Projectile flight with quadratic drag

Abstract

A shell is launched from a prescribed initial state. The dashed curve is the closed-form no-air result; the red curve integrates the drag model. The difference between the two curves is the useful part of the experiment.

Experimental setup

In this run, the launch speed is m/s at degrees from height m. The shell mass is kg, diameter m, drag coefficient , wind m/s, and the RK4 step is s.

Kinematic baseline

Ignoring air, horizontal motion is uniform and vertical motion is accelerated by gravity. This baseline makes the effect of drag visually measurable.

y(t) = y0 + v0 sin(theta)t - 12gt2

Drag model

The simulated shell is treated as a circular projectile moving through still or moving air. Air resistance grows with the square of relative velocity, so wind changes the drag force even when launch conditions are unchanged.

Fd = 12rho Cd A|v-w|2, A = pi(d/2)2

Numerical method

The red trajectory is advanced with fourth-order Runge-Kutta. Smaller time steps cost more computation but reduce integration error.

The plot uses locked x/y scale, so equal distances have equal visual length.

Figure 1. Trajectory comparison

0.00m 307.75m 615.50m 923.25m 1231.00m 1538.76m 66.79m 133.58m 200.37m 267.16m
The red trajectory is the RK4 drag simulation: range 1090.74 m, flight 13.97 s, impact speed 91.52 m/s. The dashed trajectory is the no-air solution: range 1424.77 m, peak 278.29 m.
RK4 dragno-air baseline