Preliminary Study, Liquid Propellant Engine, 2026

Ichor

Scope

O/F 5 to 8, Pc 1 to 4 MPa

Result

Design point O/F 6.4, 3.5 MPa

Method

Equilibrium thermochemistry, 2 codes

The Question

The choice of mixture ratio and chamber pressure

Ichor is an early stage propulsion study for a pressure fed, bipropellant liquid rocket engine targeting 3.5 kN sea level thrust, with hydrogen peroxide as oxidiser and kerosene as fuel.

Before any hardware is drawn, two numbers fix most of the engine: the oxidiser to fuel ratio and the chamber pressure. This page maps how specific impulse, chamber temperature, characteristic velocity and thrust coefficient respond to both, to find which variable earns its cost.

OxidiserH₂O₂
FuelC₁₂H₂₆
The Framework

A toolchain, cross checked, not yet a framework

Thermochemical equilibrium calculations with RPA1 and NASA CEA2, post processed in MATLAB. The sweep spans O/F 5 to 8 and chamber pressure 1 to 4 MPa, with the nozzle expansion ratio as a free parameter. Each point in the maps below is one converged equilibrium state.

The states are chemical equilibrium compositions at each O/F and Pc pair, expanded through an ideal nozzle. Shifting equilibrium is assumed through the expansion. The two tools were run on the same sweep and post processed together.

1 · Sweep

The sweep definition

  • O/F 5 to 8
  • Pc 1 to 4 MPa
  • H₂O₂ and kerosene
  • expansion ratio free
Defines every run.
2 · Solvers

RPA1

Rocket Propulsion Analysis, one independent implementation of the equilibrium problem.

NASA CEA2

Gordon and McBride, the second implementation.

No shared code between the two.
3 · Post processing

MATLAB

Common units, sweep assembly, the maps of Figures 1 to 3 and the design point selection.

4 · Diagnostics

Performance maps

  • Isp against O/F and Pc
  • Cf against O/F and Pc
  • Tc and c star
  • the design point
Verification, cross checked, 2 tools
Cross check

RPA against NASA CEA over the full sweep

the two tools agree, and nothing more is claimed
Analytic verification

no closed form exists for these mixtures

none was run
FIG 4The toolchain. This page is cross checked between two independent tools, not verified against an analytic solution, and the diagram says so.
The Cases

The sweep and the design point

The sweep

A 4 by 13 grid at minimum over O/F and chamber pressure, run through both tools and post processed to common units.

The design point

O/F 6.4 and Pc 3.5 MPa against a 3.5 kN sea level target, chosen for temperature margin against peak Isp. It feeds the stat cards under Finding 3.

Finding 1, Figure 1

Specific impulse is set by mixture ratio, not by pressure

Isp peaks near O/F 6.4 at every chamber pressure, and above about 2 MPa more pressure buys almost no Isp. The ★ marks the global maximum in each space.

FIG 1Specific impulse as a rotatable surface over mixture ratio and chamber pressure, with the design point marked. Drag to rotate. Mixture ratio sets specific impulse, and pressure above 2 MPa adds almost none.
Finding 2, Figure 2

Thrust coefficient is the pressure dependent quantity

Where Isp ignores chamber pressure, the thrust coefficient does not. Pressure decisions live in the nozzle, not in the chemistry.

FIG 2Thrust coefficient as a rotatable surface over mixture ratio and chamber pressure. Pressure acts through the thrust coefficient, not through the chemistry.
Finding 3, Figure 3

The design point trades peak performance for temperature margin

The sweep shows that chamber temperature keeps rising past the Isp optimum. The design point was set at O/F 6.4 and Pc 3.5 MPa, which maximises specific impulse while keeping chamber temperature within structural material limits.

FIG 3Chamber temperature as a rotatable surface over mixture ratio and chamber pressure, with the Isp against Tc cross plot beside it. The design point gives up peak impulse for temperature margin. C star peaks at 1627 m/s and the peak chamber temperature in the sweep is 2940 K.
6.4
Design O/F ratio
291 s
Peak Isp (shifting eq.)
1627 m/s
Peak C*
Discussion

What the maps decide

The two variables are not symmetric. Mixture ratio moves specific impulse and chamber temperature together, so the design point is a trade along one axis, peak performance against material margin, and O/F 6.4 buys that margin for a small Isp loss.

Chamber pressure is bought, not free. It costs tank mass and feed pressure in a pressure fed system, and above about 2 MPa the chemistry returns almost nothing for it. The 3.5 MPa point is set by the nozzle through the thrust coefficient, not by the chemistry.

A hot fire c star efficiency

One measured characteristic velocity from a firing closes the gap between equilibrium chemistry and a real chamber, the combustion efficiency this page assumes to be 1.

An injector pressure drop measurement

The feed system this page ignores. One measured drop across a candidate injector converts the chamber pressure choice from a chemistry statement into a system one.

Limits

Model limitations

Equilibrium only. No injector model. No combustion efficiency. No wall heat transfer. No feed system. Ideal expansion assumed. No hardware exists.

Sources

Principal sources

  1. Ponomarenko, A. RPA: Rocket Propulsion Analysis. Tool and documentation.
  2. Gordon, S., McBride, B. J. Computer Program for Calculation of Complex Chemical Equilibrium Compositions and Applications. NASA Reference Publication 1311, 1994.

The connection
The same modelling discipline pointed at a fuel the literature cannot yet size: Melt Layer Thermoacoustics →

Next Case Study
Paraffin-Based Hybrid Rocket Engines
View project →
Konstantinos N. Kyriakos · © 2026All workBack to top ↑