Case Study · Hybrid Rocket Combustion Dynamics · 2026

Melt Layer Thermoacoustics

Corpus

57 papers · 18 read in full

Measured FTFs

Zero, for any liquefying fuel

Method

Fixed-column extraction table

The Question

Why does changing nothing but the fuel destabilise a rocket engine?

Paraffin makes a hybrid rocket burn three times faster. It may also be the reason the engine shakes itself apart. Two published models say why, and both are built on three numbers nobody has ever measured.

Every published stability prediction for a paraffin hybrid rests on three numbers that have never been measured for a liquefying fuel, at any condition.

18 primary sources read57 papers in the corpus0 measured response functionsAugust 2026

The Problem

One engine, two fuels, one result

In 2023 Giuseppe Gallo and Carmine Carmicino ran a 100 N hybrid rocket motor twice. Same hardware. Same oxygen mass flow. They changed only the fuel.

With high density polyethylene, the motor burned stably from ignition to shutdown.

With a paraffin grain, chamber pressure built to six atmospheres and held. Then, 1.7 seconds into the burn, it entered a limit cycle. Pressure began swinging between six and nine and a half atmospheres peak to peak, and the mean pressure climbed from six to eight point four.

A fifty eight per cent pressure oscillation and a forty per cent shift in mean chamber pressure, caused by changing nothing but the fuel.

Their explanation points at the thing that makes paraffin attractive in the first place.

The Mechanism

Why paraffin burns fast

Regression rate is how fast the fuel surface recedes. It is the single number that decides whether a hybrid motor is worth building.

Conventional hybrid fuels pyrolyse straight from the solid. Paraffin does something else. It melts into a thin, low viscosity, low surface tension film on the grain surface, and the oxidiser gas tearing across that film destabilises it. Waves grow. Droplets are torn off and carried into the flow, where they vaporise and burn.

That droplet entrainment is an extra mass transfer path stacked on top of ordinary evaporation, and it is worth a factor of three to four in regression rate. It is the reason anyone puts wax in a rocket.

Karabeyoglu, Altman and Cantwell wrote the governing relation in 2002. Entrainment scales with the gas dynamic pressure raised to a power, times film thickness raised to another power, divided by melt viscosity and surface tension.

Two features of that relation matter more than they look. The mechanism is mechanical, driven by shear, with no heat needing to soak into a solid first. And it has a hard threshold: below a critical value the entrainment mass transfer is not small, it is exactly zero. The mechanism switches on.

Finding 1

Three numbers, all assumed

Two independent groups have now modelled paraffin hybrid instability. Gallo and Carmicino built an acoustic model of the coupled feed line and chamber and ran a linear stability analysis. A Turin group, publishing with Avio and the Italian Space Agency, built a multiphysics model with a quasi one-dimensional chamber solver and a neural network standing in for the chemistry.

Both rest on the same three quantities. A dynamic pressure exponent. A film thickness exponent. And the delay between a droplet leaving the surface and its heat arriving in the flow.

QuantityWhat it controlsStatus
αHow hard fuel production responds to a flow perturbationFitted to 1970s wind tunnel data on kerosene films
βSensitivity to melt film thicknessSame origin, range 1 to 2
τPhase between fuel production and chamber pressureA free parameter in one model. “Arbitrarily chosen” in the other

The Turin and Avio paper is unusually candid about this. It sets the exponents to 1.5 and 2, notes in the same paragraph that these values have been shown to produce unrealistic regression curves, and states that finding better approximations is future work. Of the delay term it says plainly that the value is arbitrarily chosen.

Written out, their entrainment term scales with oxidiser mass flux cubed. A cubed sensitivity, published, with the authors flagging that the exponent may be wrong.

Every published stability prediction for a paraffin hybrid rests on three numbers that have never been measured for a liquefying fuel, at any condition.

Finding 2 · Figure 1

The pressure nobody plots

Paraffin wax of carbon number around thirty has a critical pressure near 6.5 bar and a critical temperature near 840 K.

Above roughly 6.5 bar the wax goes supercritical. Surface tension and heat of vaporisation approach zero, and there is no clean distinction between the liquid and the gas. The droplet entrainment picture, the one all of this rests on, is a sub-critical picture.

Every flight engine runs an order of magnitude above that line.

Anna Petrarolo filmed paraffin burning at 1, 13 and 40 bar in an optical chamber rated to fifty. Below the critical point, periodic Kelvin-Helmholtz waves at 300 to 340 hertz, wavelengths of a few millimetres, droplets shedding as the theory says. Above it, no periodic wave structure could be recognised in the flame at all. What she saw instead were intermittent, violent blowing events and flame bursts, getting more frequent as pressure rose.

Her thesis closes that thread with one sentence: some of the chamber pressure oscillations could also derive from these unsteady phenomena. And then it stops.

Sorting the studies by chamber pressure has not, as far as I can find, been done. Here is the first attempt. Hover any marker for the campaign behind it.

UNSTABLE STABLE SLAB RIG no chamber 6.5 bar — wax critical pressure 16.513203040 Chamber pressure, bar Gallo & Carmicino, paraffin, 100 N onset near 6, mean 8.5, entrainment–acoustic Bertoldi et al., 2 × 1 kN, N₂O feed-system coupled HDPE Gallo & Carmicino 1 kN · Kobald 16 & 19 bar paraffin, stable Petrarolo, optical slab, 1 / 13 / 40 bar
FIG 1Every test campaign in the source set for which chamber pressure is reported, against stability outcome. Filled circles are entrainment-acoustic mechanisms, squares are feed-system coupled, the open circle is the non-liquefying control. Open ember circles are two-dimensional slab rigs, which have no chamber acoustics and therefore cannot be scored. Six of thirteen sources do not state chamber pressure at all, which for a fuel whose critical point sits in the middle of the operating range is itself a finding.

The picture is not the one I expected. Feed-coupled instability appears well above the critical pressure. Entrainment-driven acoustic instability has so far only been observed with onset in the neighbourhood of it. What nobody has done is look for the second kind, at scale, above the line.

The Lineage · Figure 2

Where the numbers came from

Two branches, one root. The place they would meet is the thing nobody has measured.

Every model in current use traces to one correlation published in 2002, which was itself fitted to wind tunnel measurements of non-reacting liquid films made in 1970. From there the field splits. One branch went to the fuel surface and photographed it. The other went to the chamber and modelled its acoustics. Neither branch has come back for the other. Hover any node.

FUEL SURFACE ENTRAINMENT THEORY CHAMBER MODELS Gater andL'Ecuyer1970 KarabeyogluParts 1 and 22002 Carmicino2009 Kobaldviscosity2017 Petrarolooptical KHI2018 Petraroloto 40 bar2020 HyEnD STERNand HyImpulse2019 Gallo andCarmicino2023 Casalinoand Avio2024 MEASUREDRESPONSEFUNCTIONdoes not exist fits entrainment law low-frequency formula film stability mode formulas acoustic model exponents 1.5 and 2 inherited adds pressure no acoustics 1 kHz sampling ceiling delay assumed delay arbitrary
FIG 2The lineage of the entrainment correlation and its two descendant branches. Solid arrows are inheritance, labelled with what is passed down. Dashed ember arrows mark where each branch stops short. The fuel-surface branch reached engine pressure but carries no chamber acoustics. The engine branch measured chamber pressure but at a sampling rate that could not resolve its own computed acoustic modes. Both modelling branches close their loop with a time delay that is assumed rather than measured. The 1970 wind tunnel data underneath all of it was taken on non-reacting films of kerosene and methanol.
Finding 3

The two models disagree

Read the papers separately and each is coherent. Read them together and they contradict each other.

Gallo and Carmicino say instability is caused by the combustion delay of entrained droplets. No entrainment, no instability. That is why the polyethylene test was quiet.

The Turin and Avio paper says, in its conclusions, that an increase in the entrainment phenomenon raises regression rate but has a beneficial effect on instability. Their design recommendation is to cut additive content so that entrainment goes up.

One says entrainment destabilises. The other says more of it stabilises. There is a related tension underneath: the optical work shows higher viscosity gives a more stable film, while the model says lower viscosity gives a more stable engine. Film stability and engine stability are being treated as the same quantity by a literature that has never said they are different.

The Hypothesis

Droplet size reconciles them

This section is inference, not a published finding. It is the part of this piece most likely to be wrong, which is why it is labelled.

The two claims fit together if droplet size is the mediating variable.

Lower melt viscosity gives more entrainment. By the same shear mechanism it also gives finer droplets. Finer droplets vaporise and burn sooner, so the delay between entrainment and heat release shortens. A shorter delay moves the phase between fuel production and chamber pressure away from the driving condition.

So more entrainment destabilises through the sheer quantity of delayed heat release, and stabilises through the shortening of the delay. Which effect wins is a question about how that delay depends on droplet diameter, and droplet diameter is set by viscosity, surface tension and dynamic pressure, which is the same group that sets entrainment in the first place.

Supporting evidence, indirect: the optical work reports that droplet size is expected to fall as chamber pressure rises, and that above the critical point droplet surfaces become indistinct altogether.

Nobody measures droplet size in a firing motor. The one dedicated droplet sizing rig in this literature runs cold, with no combustion at all.

Finding 4 · Figure 3

At flight scale, the modes converge

Every number in this section is published by the engine builders themselves. Only the arithmetic putting them on one axis is new.

The largest paraffin engine programme in Europe grew out of the University of Stuttgart's HyEnD project, whose test campaigns are the best-documented stability record in this literature. In their 500 N class engines an intrinsic low-frequency instability appeared in most tests, at 70 to 140 hertz. A Helmholtz bulk mode sat at 350 to 400 hertz, vortex shedding near 500. The first longitudinal acoustic mode of the chamber was computed at 2000 to 2500 hertz.

It was never measured. Their own campaign report says why: the acoustic modes could not be resolved because the pressure transducer sampled at 1000 hertz. Everything above 500 hertz was invisible to the instrumentation by construction.

In a small chamber that hardly matters, because the acoustic mode sits a full decade above every driving mechanism in the engine. But the first longitudinal frequency falls as the reciprocal of chamber length, and the same group has published the chamber length of their 75 kN flight engine: about 3.5 metres. Run their own formula at their own length and the mode arrives at roughly 145 to 170 hertz.

The acoustic mode that was ten times above the instability band in the test engines lands on top of it at flight scale, in the band where the instability was measured in most tests.

Move across the chart to slide a chamber length along the mode line.

5k2k50030014070 0.10.5123.55 CHAMBER LENGTH, m — log scale FREQUENCY, Hz — log MEASURED INSTABILITY 70–140 Hz NYQUIST CEILING 1 kHz probe FIRST LONGITUDINAL MODE, f = c / 2L Test chambers, 0.3–0.4 m 2000–2500 Hz computed, never measured 75 kN engine, published 3.5 m the mode arrives at 145–170 Hz
FIG 3The first longitudinal mode against chamber length, using the engine builders' own relation and their own published numbers. In the test engines the mode sat a decade above the measured instability band, and above the recording ceiling of the instrumentation. At the published 3.5 metre chamber of the 75 kN engine the mode descends into the measured band. Whether the two mechanisms couple there is exactly the kind of question a fuel response function answers, and it is unanswerable from the existing test data because the acoustic band was never recorded.

This is not a claim that the flight engine is unstable. It is a claim about evidence: the separation that made the acoustic band ignorable in the test campaigns does not survive the scale-up, and the instrumentation that produced the existing stability record could not have seen the collision if it had happened. For a sounding rocket, stability can be observed flight by flight. For a licensed launcher it has to be evidenced, at sampling rates chosen for the modes the chamber actually has.

What Would Settle It

Three measurements, each closing a named gap

A measured fuel response function

Gain and phase against forcing frequency, for a liquefying fuel. This is the closure both models assume. In gas turbine and liquid rocket combustion it is routine practice: force the system, sweep the frequency, phase-lock average against a heat release marker, extract the transfer function. It has never been done for a fuel that entrains.

The exponents, at engine pressure

Vary dynamic pressure and film thickness independently and measure the entrainment response. This is the item the Turin and Avio paper explicitly names as future work.

Droplet size under combustion

Optical sizing in a firing motor, against viscosity and pressure. The step past the cold-flow rigs, and the test of the hypothesis above.

None of these is exotic. All three are standard practice in adjacent fields. They have simply never been pointed at this fuel, at the pressure that matters.

Method

How this was built

This is a case study, so the method is part of the output.

Thirteen primary sources were read in full or by targeted section: two theory papers, four experimental campaigns, two doctoral and masters theses, two modelling papers, and three on additives and injector geometry. Each was entered into an extraction table with a fixed column set, including chamber pressure, regime against the wax critical point, instrumentation sampling rate, attributed acoustic cavity, claimed mechanism, and whether the response parameters were measured or assumed.

A second tranche of five programme papers from the HyEnD project and its successor company was read after the table was built: the STERN project report, the HEROS launch campaign paper, the LOX transition paper and two systems papers. Finding 4 is built entirely from numbers published in that tranche.

That last column is where Figure 1 and the central claim come from. Sorted, it reads: defined in one paper, assumed in two, modelled from motor data in one, not addressed in nine, measured in none.

Two database sweeps established the corpus size. The main query returns 57 records since 2015. The intersection with the supercritical question returns 15. A third query, probing for any measured response function or describing function applied to a hybrid, returns effectively nothing. A documented empty search is a finding rather than a failure, and it is why this piece can say zero rather than few.

What this cannot support

Thirteen papers is not a systematic review, and no completeness is claimed. Four sources cited inside these papers remain unread, and one of them, on film instability under supercritical conditions, sits directly under the supercritical section. The supercritical argument rests on three tests compared qualitatively in one thesis. It is the best evidence available and it is thin. The lineage reading and the hypothesis are inference, labelled as such. None of it has been checked against proprietary test data, because none was available.

Sources

Principal sources

  1. Karabeyoglu, M. A., Altman, D., Cantwell, B. J. Combustion of Liquefying Hybrid Propellants, Part 1, General Theory. Journal of Propulsion and Power 18(3), 2002.
  2. Karabeyoglu, M. A., Cantwell, B. J. Combustion of Liquefying Hybrid Propellants, Part 2, Stability of Liquid Films. Journal of Propulsion and Power 18(3), 2002.
  3. Gallo, G., Carmicino, C. Model for Combustion Instability in Hybrid Rocket Engines Burning Liquefying Fuels. AIAA Journal 61(8), 2023.
  4. Casalino, L., Ferrero, A., Folcarelli, L., Masseni, F., Muscarà, L., Pastrone, D., Frezzotti, M. L., Cretella, A., Pellegrini, R. C., Cavallini, E. Multiphysics Modeling for Combustion Instability in Paraffin-Fueled Hybrid Rocket Engines. Journal of Spacecraft and Rockets 61(3), 2024.
  5. Petrarolo, A. Liquid Layer Combustion Instabilities in Paraffin-Based Hybrid Rocket Fuels. Doctoral thesis, University of Stuttgart, 2020.
  6. Petrarolo, A., Kobald, M., Schlechtriem, S. Understanding Kelvin–Helmholtz Instability in Paraffin-Based Hybrid Rocket Fuels. Experiments in Fluids 59:62, 2018.
  7. Kobald, M., Schmierer, C., Ciezki, H. K., Schlechtriem, S., Toson, E., De Luca, L. T. Viscosity and Regression Rate of Liquefying Hybrid Rocket Fuels. Journal of Propulsion and Power 33(5), 2017.
  8. Bertoldi, A. E. M., Bouziane, M., Lee, J., Gurgel Veras, C. A., Hendrick, P., Simone, D. Theoretical and Experimental Study of Combustion Instability in Hybrid Rocket Motors. EUCASS, 2019.
  9. Mengu, D., Kumar, R. Evaluation of combustion stability of single and multi-protrusion inserted hybrid rocket motor. Acta Astronautica 219, 2024.
  10. Breitinger, J., Schmierer, C., Kobald, M., Schlechtriem, S. Launch Campaign of the Hybrid Sounding Rocket HEROS. 2017.
  11. Kobald, M., Schmierer, C., Fischer, U., Tomilin, K., Petrarolo, A., Rehberger, M. The HyEnD Stuttgart Hybrid Sounding Rocket Project. Progress in Propulsion Physics, EUCASS, 2019.
  12. Schmierer, C., Kobald, M., Fischer, U., Tomilin, K., Petrarolo, A., Hertel, F. Advancing Europe's Hybrid Rocket Engine Technology with Paraffin and LOX. EUCASS, 2019.

The measurement this calls for
I have built and run the forced-response chain this case study says is missing, on a different fuel: Thermoacoustic Response of Hydrogen Micromix Flames →

Back to start
Thermoacoustics of Hydrogen Micromix Flames
View project →
Konstantinos N. Kyriakos · © 2026All workBack to top ↑