Merge pull request #12 from Team-Deep-Impact/simulation

Simulation
This commit is contained in:
Dan
2024-05-19 03:51:32 -04:00
committed by GitHub
3 changed files with 140 additions and 60 deletions

View File

@@ -3,4 +3,4 @@ SENTRY_API_KEY=<"Your NASA Sentry API Key">
OPENAI_API_KEY=<"Your OpenAI API Key"> OPENAI_API_KEY=<"Your OpenAI API Key">
HOST_CLIENT=<"Your FrontEnd URL"> HOST_CLIENT=<"Your FrontEnd URL">
ALLOWED_HOSTS=<"Your Server URL"> ALLOWED_HOSTS=<"Your Server URL">
DEBUG=<"True for Deployment | False for Development"> DEBUG=<"True for Deployment | False for Development">

3
backend/.gitignore vendored
View File

@@ -1,2 +1,3 @@
.env .env
.venv .venv
/__pycache__/

View File

@@ -1,9 +1,8 @@
import React, { useState, useEffect, useRef } from 'react'; import React, { useState, useEffect, useRef } from 'react';
import { Stage, Layer, Circle, Text, Line } from 'react-konva'; import { Stage, Layer, Circle, Text, Line } from 'react-konva';
import { Slider, Button } from '@mui/material'; import { Slider, Button, MenuItem, Select, Typography, Box } from '@mui/material';
import axios from 'axios'; import axios from 'axios';
const EARTH_RADIUS_KM = 6371; // Earth's radius in kilometers const EARTH_RADIUS_KM = 6371; // Earth's radius in kilometers
const EARTH_DISPLAY_SCALE = 0.01; // Display scale for Earth const EARTH_DISPLAY_SCALE = 0.01; // Display scale for Earth
const ASTEROID_DISPLAY_SCALE = 0.3; // Display scale for asteroids const ASTEROID_DISPLAY_SCALE = 0.3; // Display scale for asteroids
@@ -16,6 +15,7 @@ const CANVAS_HEIGHT = window.innerHeight;
const EARTH_X = CANVAS_WIDTH - 100; const EARTH_X = CANVAS_WIDTH - 100;
const EARTH_Y = CANVAS_HEIGHT / 2; const EARTH_Y = CANVAS_HEIGHT / 2;
const SIMULATION_AREA_MULTIPLIER = 2; // Extend the simulation area const SIMULATION_AREA_MULTIPLIER = 2; // Extend the simulation area
const ROCKET_SPEED = 0.025; // Reduced speed of the rocket
function Scenario() { function Scenario() {
const [asteroid, setAsteroid] = useState({ const [asteroid, setAsteroid] = useState({
@@ -35,11 +35,21 @@ function Scenario() {
const [simulationResult, setSimulationResult] = useState(''); const [simulationResult, setSimulationResult] = useState('');
const [strategy, setStrategy] = useState(''); // State for selected strategy const [strategy, setStrategy] = useState(''); // State for selected strategy
const [timeStep, setTimeStep] = useState(100); // State for time step const [timeStep, setTimeStep] = useState(100); // State for time step
const [rocket, setRocket] = useState({
active: false,
x: EARTH_X,
y: EARTH_Y,
trajectory: [],
vx: 0,
vy: 0,
});
const [rocketDistance, setRocketDistance] = useState(10000);
const [collisionOccurred, setCollisionOccurred] = useState(false);
useEffect(() => { useEffect(() => {
const fetchAsteroids = async (self, request) => { const fetchAsteroids = async () => {
try { try {
const response = await axios.get(`http://localhost:8000/api/sentry/`); const response = await axios.get('http://localhost:8080/api/sentry/');
asteroidsData.current = response.data.near_earth_objects; asteroidsData.current = response.data.near_earth_objects;
setAsteroidData(0); // Set initial asteroid setAsteroidData(0); // Set initial asteroid
} catch (error) { } catch (error) {
@@ -72,23 +82,15 @@ function Scenario() {
}); });
setTrajectoryPoints([]); setTrajectoryPoints([]);
calculateForecastPath(size, speed, Math.PI / 9, mass); // Calculate the forecast path once calculateForecastPath(size, speed, Math.PI / 9, mass); // Calculate the forecast path once
setRocket({ active: false, x: EARTH_X, y: EARTH_Y, trajectory: [], vx: 0, vy: 0 });
setRocketDistance(10000);
setCollisionOccurred(false);
}; };
const handleAsteroidSelect = (event) => { const handleAsteroidSelect = (event) => {
setAsteroidData(event.target.value); setAsteroidData(event.target.value);
}; };
const applyDeflectionStrategy = (vx, vy) => {
if (strategy === 'Nuclear Detonation') {
vx += 0.5; // Example effect
} else if (strategy === 'Kinetic Impact') {
vx += 0.25; // Example effect
} else if (strategy === 'Gravity Tractor') {
vy -= 0.1; // Example effect
}
return { vx, vy };
};
const calculateForecastPath = (size, speed, angle, mass) => { const calculateForecastPath = (size, speed, angle, mass) => {
let x = 100; let x = 100;
let y = CANVAS_HEIGHT / 2 - 100; let y = CANVAS_HEIGHT / 2 - 100;
@@ -135,6 +137,20 @@ function Scenario() {
setForecastPoints(points); setForecastPoints(points);
}; };
const applyVelocityImpulse = (asteroid, strategy, rocketDx, rocketDy) => {
let impulseX = 0;
let impulseY = 0;
const angle = Math.atan2(rocketDy, rocketDx);
if (strategy === 'Nuclear Detonation') {
impulseX = 0.01 * Math.cos(angle);
impulseY = 0.01 * Math.sin(angle);
} else if (strategy === 'Kinetic Impact') {
impulseX = 0.005 * Math.cos(angle);
impulseY = 0.005 * Math.sin(angle);
}
return { vx: asteroid.vx + impulseX, vy: asteroid.vy + impulseY };
};
const simulate = () => { const simulate = () => {
let { x, y, speed, angle, mass } = asteroid; let { x, y, speed, angle, mass } = asteroid;
@@ -142,8 +158,6 @@ function Scenario() {
let vy = speed * Math.sin(angle); let vy = speed * Math.sin(angle);
const updateFrame = () => { const updateFrame = () => {
console.log(`Asteroid position: (${x}, ${y}), velocity: (${vx}, ${vy})`);
const dx = EARTH_X - x; const dx = EARTH_X - x;
const dy = EARTH_Y - y; const dy = EARTH_Y - y;
const distance = Math.sqrt(dx * dx + dy * dy); const distance = Math.sqrt(dx * dx + dy * dy);
@@ -162,19 +176,57 @@ function Scenario() {
vx += ax * timeStep; vx += ax * timeStep;
vy += ay * timeStep; vy += ay * timeStep;
({ vx, vy } = applyDeflectionStrategy(vx, vy));
x += vx * timeStep; x += vx * timeStep;
y += vy * timeStep; y += vy * timeStep;
// Update rocket position
setRocket((prev) => {
const rocketDx = x - prev.x;
const rocketDy = y - prev.y;
const rocketDistance = Math.sqrt(rocketDx * rocketDx + rocketDy * rocketDy);
const rocketAngle = Math.atan2(rocketDy, rocketDx);
const rocketVX = ROCKET_SPEED * Math.cos(rocketAngle);
const rocketVY = ROCKET_SPEED * Math.sin(rocketAngle);
const newX = prev.x + rocketVX * timeStep;
const newY = prev.y + rocketVY * timeStep;
if (rocketDistance <= asteroid.size * ASTEROID_DISPLAY_SCALE && !collisionOccurred) {
const newVelocity = applyVelocityImpulse({ vx, vy }, strategy, rocketDx, rocketDy);
vx = newVelocity.vx;
vy = newVelocity.vy;
setSimulationResult('Rocket intercepted the asteroid!');
setCollisionOccurred(true); // Mark collision occurred
return { ...prev, active: false, x: newX, y: newY, vx: 0, vy: 0 }; // Stop the rocket
}
return {
...prev,
x: newX,
y: newY,
vx: prev.active ? rocketVX : 0, // Stop velocity if inactive
vy: prev.active ? rocketVY : 0, // Stop velocity if inactive
trajectory: [...prev.trajectory, newX, newY],
};
});
setRocketDistance(rocketDistance);
if (rocketDistance <= asteroid.size * ASTEROID_DISPLAY_SCALE && !collisionOccurred) {
setSimulationResult('Rocket intercepted the asteroid!');
setIsSimulating(false);
return;
}
if (y > CANVAS_HEIGHT * SIMULATION_AREA_MULTIPLIER || x > CANVAS_WIDTH * SIMULATION_AREA_MULTIPLIER || y < -CANVAS_HEIGHT * SIMULATION_AREA_MULTIPLIER || x < -CANVAS_WIDTH * SIMULATION_AREA_MULTIPLIER) { if (y > CANVAS_HEIGHT * SIMULATION_AREA_MULTIPLIER || x > CANVAS_WIDTH * SIMULATION_AREA_MULTIPLIER || y < -CANVAS_HEIGHT * SIMULATION_AREA_MULTIPLIER || x < -CANVAS_WIDTH * SIMULATION_AREA_MULTIPLIER) {
setSimulationResult('Asteroid missed Earth!'); setSimulationResult('Asteroid missed Earth!');
setIsSimulating(false); setIsSimulating(false);
return; return;
} }
setAsteroid(prev => ({ ...prev, x, y })); setAsteroid((prev) => ({ ...prev, x, y }));
setTrajectoryPoints(prev => [...prev, x, y]); setTrajectoryPoints((prev) => [...prev, x, y]);
animationRef.current = window.requestAnimationFrame(updateFrame); animationRef.current = window.requestAnimationFrame(updateFrame);
}; };
@@ -184,8 +236,13 @@ function Scenario() {
}; };
const handleSimulate = () => { const handleSimulate = () => {
if (isSimulating) return; if (isSimulating) {
simulate(); setIsSimulating(false);
window.cancelAnimationFrame(animationRef.current);
setSimulationResult('Simulation stopped.');
} else {
simulate();
}
}; };
const handleTimeStepChange = (event, newValue) => { const handleTimeStepChange = (event, newValue) => {
@@ -193,51 +250,67 @@ function Scenario() {
}; };
const handleAngleChange = (event, newValue) => { const handleAngleChange = (event, newValue) => {
setAsteroid(prev => ({ ...prev, angle: newValue })); setAsteroid((prev) => ({ ...prev, angle: newValue }));
calculateForecastPath(asteroid.size, asteroid.speed, newValue, asteroid.mass); calculateForecastPath(asteroid.size, asteroid.speed, newValue, asteroid.mass);
}; };
return ( return (
<> <>
<div> <Box sx={{ padding: '20px', maxWidth: '800px', margin: '0 auto', textAlign: 'center' }}>
<h1>Asteroid Collision Scenario</h1> <Typography variant="h4" gutterBottom>
<select onChange={handleAsteroidSelect}> Asteroid Collision Scenario
{asteroidsData.current.map((ast, index) => ( </Typography>
<option key={index} value={index}>{ast.name}</option> <Box sx={{ marginBottom: '20px' }}>
))} <Select
</select> value={asteroid.name}
<select onChange={(e) => setStrategy(e.target.value)}> onChange={handleAsteroidSelect}
<option value="">Select Deflection Strategy</option> displayEmpty
<option value="Nuclear Detonation">Nuclear Detonation</option> fullWidth
<option value="Kinetic Impact">Kinetic Impact</option> >
<option value="Gravity Tractor">Gravity Tractor</option> {asteroidsData.current.map((ast, index) => (
</select> <MenuItem key={index} value={index}>
{ast.name}
</MenuItem>
))}
</Select>
</Box>
<Box sx={{ marginBottom: '20px' }}>
<Select
value={strategy}
onChange={(e) => setStrategy(e.target.value)}
displayEmpty
fullWidth
>
<MenuItem value="">Select Deflection Strategy</MenuItem>
<MenuItem value="Nuclear Detonation">Nuclear Detonation</MenuItem>
<MenuItem value="Kinetic Impact">Kinetic Impact</MenuItem>
</Select>
</Box>
<Button <Button
variant="contained" variant="contained"
color="primary" color="primary"
onClick={handleSimulate} onClick={handleSimulate}
disabled={isSimulating} style={{ marginBottom: '20px' }}
style={{ marginTop: '10px' }}
> >
Simulate {isSimulating ? 'Stop Simulation' : 'Simulate'}
</Button> </Button>
{simulationResult && ( {simulationResult && (
<div> <Box sx={{ marginBottom: '20px' }}>
<h2>Simulation Result</h2> <Typography variant="h6">Simulation Result</Typography>
<p>{simulationResult}</p> <Typography>{simulationResult}</Typography>
</div> </Box>
)} )}
{asteroid.name && ( {asteroid.name && (
<div> <Box sx={{ marginBottom: '20px' }}>
<h3>Asteroid Information</h3> <Typography variant="h6">Asteroid Information</Typography>
<p><strong>Name:</strong> {asteroid.name}</p> <Typography><strong>Name:</strong> {asteroid.name}</Typography>
<p><strong>Size:</strong> {asteroid.size.toFixed(3)} kilometers</p> <Typography><strong>Size:</strong> {asteroid.size.toFixed(3)} kilometers</Typography>
<p><strong>Speed:</strong> {(asteroid.speed * 1000).toFixed(3)} km/s</p> <Typography><strong>Speed:</strong> {(asteroid.speed * 1000).toFixed(3)} km/s</Typography>
<p><strong>Mass:</strong> {asteroid.mass.toExponential(3)} kg</p> <Typography><strong>Mass:</strong> {asteroid.mass.toExponential(3)} kg</Typography>
</div> </Box>
)} )}
<div> <Box sx={{ marginBottom: '20px' }}>
<h3>Time Step: {timeStep}</h3> <Typography>Time Step: {timeStep}</Typography>
<Slider <Slider
value={timeStep} value={timeStep}
min={10} min={10}
@@ -246,9 +319,9 @@ function Scenario() {
onChange={handleTimeStepChange} onChange={handleTimeStepChange}
aria-labelledby="time-step-slider" aria-labelledby="time-step-slider"
/> />
</div> </Box>
<div> <Box sx={{ marginBottom: '20px' }}>
<h3>Initial Angle: {(asteroid.angle * 180 / Math.PI).toFixed(2)}°</h3> <Typography>Initial Angle: {(asteroid.angle * 180 / Math.PI).toFixed(2)}°</Typography>
<Slider <Slider
value={asteroid.angle} value={asteroid.angle}
min={0} min={0}
@@ -257,8 +330,8 @@ function Scenario() {
onChange={handleAngleChange} onChange={handleAngleChange}
aria-labelledby="angle-slider" aria-labelledby="angle-slider"
/> />
</div> </Box>
</div> </Box>
<Stage width={CANVAS_WIDTH} height={CANVAS_HEIGHT}> <Stage width={CANVAS_WIDTH} height={CANVAS_HEIGHT}>
<Layer> <Layer>
<Text text="Earth" fontSize={20} x={EARTH_X - EARTH_RADIUS_KM * EARTH_DISPLAY_SCALE} y={EARTH_Y - 10} /> <Text text="Earth" fontSize={20} x={EARTH_X - EARTH_RADIUS_KM * EARTH_DISPLAY_SCALE} y={EARTH_Y - 10} />
@@ -270,6 +343,12 @@ function Scenario() {
{forecastPoints.length > 0 && ( {forecastPoints.length > 0 && (
<Line points={forecastPoints} stroke="gray" strokeWidth={2} dash={[10, 10]} /> <Line points={forecastPoints} stroke="gray" strokeWidth={2} dash={[10, 10]} />
)} )}
{rocket.active && (
<Circle x={rocket.x} y={rocket.y} radius={5} fill="red" />
)}
{rocket.trajectory.length > 0 && (
<Line points={rocket.trajectory} stroke="green" strokeWidth={2} />
)}
</Layer> </Layer>
</Stage> </Stage>
</> </>