Case Study, Two-Phase Blowdown, 2026

Compressed-Air / Water Mixture Rocket Performance

Scope

Launch rod based flight performance divided into 3 thrust phases

Result

Wall shear costs 0.58 to 8.1 percent

Method

Transient water column momentum

The Question

What does the archived model of a toy contain, and what did its own author leave unfinished?

A 2 litre bottle part filled with water and pressurised with air is a two phase blowdown with three thrust regimes. The launch tube as a piston, the water impulse driven by the expanding gas, and the gas impulse after water exhaustion. The physics is Wheeler 2002 in the bottle fixed frame1, and the archived Simulink implementation of it ran on a quasi steady Bernoulli velocity and dry air5.

The mechanism

The water phase is momentum limited and the gas phase is energy limited, and that asymmetry decides every finding below.

The unfinished parts

Wheeler's Appendix B, the condensation of the water vapour carried in the pressurising gas, and Appendix A, the nozzle boundary layer, were both left incomplete in the source. I completed both, here1.

The Framework

What waterrocketperf computes

1 · Properties

Bottle

The axisymmetric profile with closed form volume and inertia integrals.

Atmosphere and water

Ambient state and waterProperties.

cases

Named parameter sets verbatim from their sources.

Knows no solvers.
2 · Closures

Gas processes

DryAdiabatic, the parity baseline. WheelerCondensation, Appendix B all at once1. GommesMoistAdiabat, the continuous moist adiabat2.

Nozzle

Choked and subsonic compressible outflow.

BoundaryLayer

Appendix A completed34.

Knows no rockets.
3 · Solver

WaterRocket

Through launch tube, water impulse and gas impulse.

  • out: thrust history
  • out: phase impulses
  • out: state traces
4 · Diagnostics

designMaps

Fill fraction and pressure sweeps.

payloadRatio

Flight performance.

Verification, waterrocketperf v1.1.0, 45 tests
Cross check, phase by phase

against the archived trace5, launch tube, water, gas

0.64, 0.001 and 7.84 percent, impulse over the archived window 1.14 percent
Verified on geometry

the closed form volume round trip

exact to 1e-12, and 11 times faster
Cross check, condensation

Wheeler's own worked example1

reproduced to 0.007 percent
Gas impulse gap

the archive truncated its own gas tail at 0.088 s with 6.33 N still flowing

7.84 percent, stated, not smoothed
FIG 4Framework architecture. Property sets know no solvers, closures know no rockets, and the solver takes every closure by injection. The largest deviation on the page sits in the verification row, not in a footnote.
The Cases

Four named runs

The cross-check case

The framework run on the legacy geometry against reference_thrust.csv, 1764 points at 50 microseconds, like for like5.

The reference bottle

The standard 2 litre case behind every sweep.

The condensation sweep

Dry baseline against both moist closures12.

The friction sweep

Bore diameter and neck length swept to find where Appendix A starts to matter3.

Finding 1, Figure 1

The rebuild reproduces the archive phase by phase, and says where they part

Thrust against time, archive and framework overlaid, three phases shaded. The water phase matches to 0.001 percent while the gas phase reads 7.84 percent high against an archive that its own Stop block truncated at 0.088 s with 6.33 N still flowing.

Phase zoom
FIG 1Thrust against time for the archived Simulink trace and the framework, three phases shaded. The 2 codes agree on the water phase, which carries 75 percent of the impulse.

Drawn from the page's own march of the same equations. The exported parity run against reference_thrust.csv will replace it.

Finding 2, Figure 2

Condensation is a pressure effect the impulse barely sees

Two treatments of the same vapour, all at once against the continuous moist adiabat12. Measured on the reference case, water impulse moves 0.802 percent and gas impulse 0.155 percent.

Vapour treatment
FIG 2Gas temperature as a rotatable surface over expansion and initial humidity, dry baseline drawn as the dashed edge, with the treatment toggle kept. The vapour changes the pressure trace and moves the impulse less than 1 percent. Measured on the reference case, condensation moves water impulse 0.802 percent and gas impulse 0.155 percent.
Finding 3, Figure 3

Wheeler predicted the friction and mispriced it

The completed Appendix A splits the effect. The displacement thickness area reduction on its own is impulse neutral at +0.000 percent, because narrowing the exit lowers mass flow and raises velocity in almost equal measure, while the wall shear costs 0.579 percent on the standard bottle, 1.8 percent at a 4 mm bore and 8.1 percent at a 300 mm neck34. The laminar method carries a caveat at the measured nozzle Reynolds number of 8.8e5.

Anchors from the sweep: +0.000 percent from the area reduction alone, 0.579 percent on the standard bottle, 1.8 percent at a 4 mm bore, 8.1 percent at a 300 mm neck. Laminar method, nozzle Reynolds number 8.8e5, past transition.
FIG 3Impulse penalty as a rotatable surface over bore diameter and neck length. The wall shear, not the area reduction, carries the whole friction cost. Every number is a lower bound, because the method is laminar and the measured Reynolds number sits past transition.
Discussion

One asymmetry organises everything

Restrictions and vapour physics move the water phase, which is momentum limited, and barely touch the gas phase, which is energy limited, because the same stored gas energy eventually leaves whatever the restriction. The 7.84 percent gas phase gap against the archive and the laminar caveat are the two open edges, and both are stated where they occur.

A bench thrust trace

With a pressure transducer and controlled initial humidity, which separates the two condensation treatments.

A turbulent completion of Appendix A

The measured Reynolds number sits past transition, so the laminar result is a lower bound on the friction.

A high speed video of water exhaustion

The instant where the archive and the framework disagree about the handover.

Limits

Model limitations

Laminar boundary layer only, so the friction numbers are lower bounds. No spray, foaming or two phase expulsion at water exhaustion. Rigid bottle, no compliance. One dimensional nozzle flow. No validation against a measured thrust trace.

Sources

Principal sources

  1. Wheeler, D. Modeling the thrust phase of a water rocket. 2002. The three phase model and both appendices.
  2. Gommes, C. J. A more thorough analysis of water rockets, moist adiabats, transient flows, and inertial forces in a soda bottle. American Journal of Physics 78(3), 2010.
  3. Thwaites, B. Approximate calculation of the laminar boundary layer. Aeronautical Quarterly 1, 1949.
  4. Cebeci, T., Bradshaw, P. Momentum Transfer in Boundary Layers. Hemisphere, 1977.
  5. Archived Simulink implementation of the Wheeler model. The parity artefact, reference_thrust.csv, 1764 points at 50 microseconds.
Back to start
Thermoacoustics of Hydrogen Micromix Flames
View project →
Konstantinos N. Kyriakos · © 2026All workBack to top ↑