Rocket launch

Will an L class rocket fly straight, reach 11,000 ft and land softly in a 10 mph wind?

You ask

Will this L class rocket fly straight, reach 11,000 ft and land softly in a 10 mph wind?

Wobbl runs

Terminal
$ wobbl motion rocket.toml

Modeled6 state variables

  • gravity falling off with altitude
  • 1976 standard atmosphere
  • air drag
  • transonic drag rise
  • wind
  • launch rail
  • thrust curve
  • propellant burn
  • weathercocking
  • drogue parachute
  • main parachute

Checks

  • passstep convergence
  • passdelta v budget
  • passparachute descent
  • passenergy balance
  • passunits and sanity
  • passstability margin
  • passrail exit speed
  • infomax dynamic pressure
  • infomax acceleration
  • passapogee
  • passlanding speed
  • infodrift

Events

EventTimeAltitudeSpeed
Liftoff0.0014 s0 m0 m/s
Rail exit0.173 s2.4 m30.6 m/s
Burnout2.65 s647 m419 m/s
Apogee23.5 s3,439 m23.1 m/s
Drogue deploy23.5 s3,439 m23.1 m/s
Main deploy178 s150 m20.4 m/s
Landing203 s0 m7.26 m/s
Altitude from liftoff to landing.

OKReaches 3439 m (11284 ft) at Mach 1.27 and lands 29.3 m upwind of the pad at 5.72 m/s, with a 1.88 caliber margin and 30.6 m/s off the rail.

Apogee
3439 m
Max Mach
1.27
Stability
1.88 cal
Confidence
high
summary.txt
summary.txt
OK. Reaches 3439 m (11284 ft) at Mach 1.27 and lands 29.3 m upwind of the pad at 5.72 m/s, with a 1.88 caliber margin and 30.6 m/s off the rail. Confidence: high. passed step convergence, delta v budget, parachute descent, energy balance, units and sanityapogee: 3439 mmax Mach: 1.27stability: 1.88 caloff the rail: 30.6 m/slanding: 5.72 m/sdrift: 29.3 m Checks:  [pass] step convergence: flown again with a tolerance ten times tighter, the flight time 203 s, apogee 3439 m, burnout speed 419 m/s, max Mach 1.27, landing speed 5.72 m/s and drift 29.3 m moved by at most 0.0098% (the drift)  [pass] delta v budget: Tsiolkovsky gives 541 m/s from a 2007 m/s exhaust speed and a mass ratio of 1.31. Gravity takes 25.7 m/s, drag 96.9 m/s and steering 0.0863 m/s, which leaves 419 m/s, and the flight reaches 419 m/s at burnout (7.2e-10% apart)  [pass] parachute descent: Under the drogue sqrt(2 m g / (rho Cd A)) gives 19.9 m/s at 150 m and the flight sinks at 19.9 m/s there. Under the drogue and main it gives 5.72 m/s at the ground against 5.72 m/s from the run (0.097% apart at most)  [pass] energy balance: thrust did 867 kJ of work, drag took 798 kJ, gravity took 3.43 kJ and the burned propellant carried off 65.6 kJ, which leaves 160 J of kinetic energy at landing, and the landing speed carries 160 J (apart by 4.4e-06% of the largest term)  [pass] units and sanity: speeds stay under Mach 1.27, all 1.88 kg of propellant burned, an average thrust 18.3 times the liftoff weight, 7.95 kg in a 98 mm by 2.1 m tube, 502 kg/m3 on average  [pass] stability margin: 1.88 calibers at liftoff, the center of pressure 1.54 m from the nose tip and the balance point 1.36 m, rising to 3.38 at burnout as the propellant burns away. Barrowman's method holds below about Mach 0.8  [pass] rail exit speed: 30.6 m/s after 2.44 m of rail, above the 15 m/s the fins need to steer  [info] max dynamic pressure: 106 kPa at Mach 1.27 and 574 m, 2.48 s after ignition, the load the fins and their joints have to carry without fluttering  [info] max acceleration: 20.8 g (204 m/s2) at 0.03 s under thrust. The drogue opens with 1.3 g of drag, pulling 79.2 N on its harness. The main opens with 3.6 g of drag, pulling 214 N on its harness  [pass] apogee: 3439 m (11284 ft) above the ground at 23.5 s, 2.6% over the 3353 m target  [pass] landing speed: descends at 5.72 m/s under the drogue and main, under the 6.1 m/s limit, and with the wind it moves 7.26 m/s over the ground  [info] drift: lands 29.3 m upwind of the pad in a 4.47 m/s wind, after reaching 757 m upwind at apogee. The descent from apogee takes 179 s Assumptions:  - A point mass in a vertical plane, so it does not spin, tumble or move sideways, and its drag acts against the air rushing past it.  - Drag coefficient 0.45 below Mach 0.6, raised through the transonic range to 0.765 near Mach 1.1 by a generic curve for a slender finned rocket, the same whether the motor burns or not.  - The 1976 US standard atmosphere from 0 m above sea level, a standard day.  - Gravity falls off with the square of the distance to the Earth's center, over a flat Earth that does not rotate.  - A steady wind of 4.47 m/s at every height, blowing back toward the start.  - Thrust follows the 75 mm L class curve in straight lines between its points, 3773 N s over 2.65 s, and propellant burns in proportion to thrust, so all 1.88 kg is gone at burnout.  - The rocket travels the full 2.44 m of rail before it leaves it, tilted 5 deg from vertical.  - Off the rail it turns into the wind it feels at its weathercock frequency, from Barrowman's normal force and margin and the pitch inertia of a uniform rod of its length and mass, with no lift at angle of attack and no overshoot.  - Each parachute opens over the time it takes to travel eight of its diameters, its drag area growing with the square of time, and the drogue stays on under the main.  - Under the parachutes only their drag counts, so the descent is slightly faster than with the airframe's drag added, which errs on the safe side for landing speed. Images: flight: images/flight.png, path: images/path.pngReport: report.html

Open the full report

About this example

This is a level 2 high power rocket, a 98 mm fiberglass airframe 2.1 m long with an ogive nose and three clipped delta fins, flying on a 75 mm L class motor of 3773 N s that burns for 2.65 s. It weighs 7.95 kg on the pad and leaves an 8 ft rail tilted 5 degrees into a 10 mph wind. An altimeter fires an 18 inch drogue at apogee and a 60 inch main at 150 m on the way down. Everything is in rocket.toml, with comments on the values that need them.

Result Value
Apogee 3439 m (11284 ft) at 23.5 s, 2.6% over the 11,000 ft target
Top speed 427 m/s, Mach 1.27, just before burnout at 2.65 s
Stability 1.88 calibers off the rail by Barrowman’s method, rising to 3.38 as the propellant burns
Off the rail 30.6 m/s, twice the 15 m/s the fins need to steer
Peak loads 106 kPa of dynamic pressure at top speed and 20.8 g on the thrust peak
Descent 19.9 m/s under the drogue, then 5.72 m/s under both parachutes
Landing 29.3 m upwind of the pad, 203 s after ignition
Altitude and speed over the flight, colored by phase

The rocket passes Mach 1 at 1.76 s. Through the transonic range its drag coefficient climbs from 0.45 to 0.77, and by burnout drag has taken 96.9 m/s of the 541 m/s the motor could give it. Off the rail it turns into the wind it feels, so it climbs upwind and is 757 m from the pad at apogee, and the wind carries it back during the three minute descent until it lands near where it started.

The path through the air, upwind on the way up and drifting back under the parachutes

Wobbl checks itself on the way. Tsiolkovsky’s ideal delta v less what gravity, drag and steering took matches the speed at burnout, the work done by thrust, drag and gravity matches the change in kinetic energy, the sink rate under each parachute matches its terminal speed, and every number held when the flight was integrated again with a tolerance ten times tighter. The model is a point mass in a vertical plane that turns into the relative wind at its weathercock frequency, which is the usual depth for planning a flight but not a six degree of freedom simulation, and Barrowman’s method is subsonic, so the margin past Mach 1 needs a supersonic method. At 106 kPa it is also worth checking the fins for flutter.

The full flight is in output/trajectory.csv, at most 500 rows with one at every event, each with time, downrange distance, altitude, speed, acceleration, mass, Mach number, dynamic pressure and phase.

Run it yourself

Terminal
$ git clone https://github.com/sohan-bhat/wobbl$ cd wobbl/examples/rocket-launch && ./run.sh

The part, the script and every output are in examples/rocket-launch on GitHub.