Calculating the Path of Intercontinental Ballistic Missiles

Initial Conditions and Projectile Motion

The trajectory of an ICBM can be initially described using the equations of projectile motion, assuming constant gravitational acceleration (g) acting downwards. The initial velocity vector (v₀), launch angle (θ), and atmospheric conditions are crucial inputs. The horizontal component of velocity remains relatively constant if air resistance is neglected, while the vertical component changes due to gravity.

For a simplified case neglecting drag, the position of the missile as a function of time can be described by: x(t) = v₀ * cos(θ) * t and y(t) = v₀ * sin(θ) * t - 0.5 * g * t².

Aerodynamic Drag – A Significant Force

At the high speeds characteristic of ICBMs, aerodynamic drag becomes a dominant force affecting trajectory. The drag force (Fd) is proportional to the square of the velocity and the drag coefficient (Cd), which depends on the missile’s shape and Reynolds number. The equation for drag is Fd = 0.5 * ρ * Cd * A * v², where ρ is air density, A is the cross-sectional area, and v is the velocity.

Impact of Drag on Trajectory

The inclusion of drag significantly alters the trajectory. The horizontal component of velocity decreases over time due to the drag force acting opposite to motion. The vertical component is also affected, though typically to a lesser extent unless the missile’s shape creates significant induced drag. Numerical integration methods are often employed to solve for the trajectory accurately when drag is present.

Consider a simplified scenario where drag is constant; the change in horizontal velocity (Δvₓ) over a small time interval (Δt) can be approximated as Δvₓ = - (Fd/m) * Δt, where m is the missile’s mass. This allows for iterative calculations of position and velocity.

Pressure Curves and Internal Ballistics

ICBMs utilize multi-stage rocket engines, each with its own pressure curve – a graph depicting the thrust generated as a function of chamber pressure. The shape of this curve is critical to performance. The change in momentum (Δp) imparted by a small mass flow rate (ṁ) over a time interval (Δt) is given by Δp = ṁ * v, where v is the exhaust velocity. Integrating this equation provides the thrust force.

The cumulative effect of these pressure curves dictates the overall trajectory, particularly during boost phases when the missile accelerates rapidly.

Coriolis Deflection – A Subtle Effect

Due to the Earth’s rotation, a force known as the Coriolis deflection acts on moving objects. This force is perpendicular to both the object's velocity and the axis of rotation. For an ICBM launched north of the equator, the deflection will be to the right; south of the equator, it will be to the left. The magnitude of this effect is proportional to the missile’s velocity and the sine of the latitude.

The Coriolis force (Fc) can be approximated as Fc = 2 * m * ω * v * sin(φ), where m is mass, ω is the angular velocity of Earth's rotation, v is the missile’s velocity, and φ is the launch latitude.

Trajectory Refinement – Iterative Solutions

Accurately predicting an ICBM trajectory requires iterative numerical methods such as Runge-Kutta or similar algorithms. These methods break down the time into small increments and repeatedly calculate the forces acting on the missile, updating its velocity and position at each step. The accuracy of the solution depends on the size of the time steps.

The process involves solving a differential equation system: d²x/dt² = Fₓ (drag, gravity, Coriolis), d²y/dt² = Fy (drag, gravity, Coriolis) where Fx and Fy are vector forces.

Frequently asked questions

Try it live

Everything above runs in your browser — open Rocket Launch Trajectory and change the parameters while it is running. Nothing is installed, nothing is uploaded, the whole model lives in one tab.