.. jannaf
.. _ref_to_jannaf_standard:
JANNAF Standard
===============
In 1975, the `John Hopkins University, Applied Physics Laboratory `_
published `CPIA 246 `_ in conjunction with
`JANNAF (Joint Army-Navy-NASA-Air Force) `_ titled
`Rocket Engine Performance Test Data Acquisition and Interpretation Manual `_
:ref:`CPIA 246 JANNAF Rocket Engine Performance Prediction and Evaluation Manual ` established a US
national standard for the analytical and experimental evaluation of the performance of
liquid propellant rocket engines and is a result of the work by the
`Interagency Chemical Rocket Propulsion Group (ICRPG) `_
formed in 1962. In 1973,
`NASA-CR-131519 `_
was published to discuss the history of the `ICRPG `_ meetings that led to
:ref:`CPIA 246 ` and to
summarize the selected model approaches.
In 1977, `Aerojet Liquid Rocket Company (ALRC) `_, now known as
`Aerojet Rocketdyne `_,
published an update to :ref:`CPIA 246 `
entitled `Analysis of rocket engine injection combustion processes `_
that reflects on-going improvements to the models and software used to implement the
"standard" JANNAF model. The `Aerojet `_ report also indicates that,
in 1977, the `JANNAF `_ model, in general,
could not yet achieve a 1% accuracy on predicting Isp without the help of test data.
There are two model approaches offered by `JANNAF `_,
a rigorous model and a simplified model.
Rigorous Model
--------------
The rigorous model is intended for the detailed analysis of an engine that is preparing for a
hardware production design; either test or flight hardware.
The rigorous model requires many more resources and much more technical expertise than
can be presented here.
Since 1977, much has changed in how aerospace companies conduct rigorous performance analysis.
The rigorous approach has evolved.
The advent of `Computational fluid dynamics (CFD) `_
has given rise to analysis approaches like the NASA Space Launch System (SLS) at
`NASA Ames Research Center `_ and
`NASA Marshall Space Flight Center (MSFC) `_.
Such models require super computers or large computer clusters in order to complete a single run
in a reasonable time (often several days of run time).
Simplified Model
----------------
The simplified model is intended more for preliminary analysis where general trends and influences are more important
than low-level design details.
The simplified model can be performed with some
first-order efficiency models and the `python `_ project
`RocketCEA `_
(`RocketCEA `_ wraps
the `NASA FORTRAN CEA code `_ that calculates
One Dimensional Equilibrium Isp)
**Rocketisp** will focus on a simplified model of performance prediction.
`NASA-CR-131519 `_
describes the simplified model for calculating delivered Isp (:math:`Isp_{del}`) as having 3 basic assumptions::
1 One-dimensional isentropic expansion (ODE)
2. Simple nozzle geometry
3. A system of correction losses to the ODE value
a. Mixture ratio distribution
b. Chamber friction
c. Chamber heat transfer
d. Energy release efficiency
e. Nozzle friction
f. Nozzle heat transfer
g. Nozzle geometry
h. Kinetics
i. Mass distribution
j. Two-phase flow
k. Nonstandard inlet conditions
In its most simplified form, the delivered Isp (:math:`Isp_{del}`) of a liquid propellant rocket engine
with a conventional nozzle
is just the One Dimensional Equilibrium Isp
(:math:`Isp_{ODE}`) generated by `RocketCEA `_ multiplied
times each of the efficiencies (:math:`\eta_i`) described by the correction losses, a through k, above.
.. math::
\Large{Isp_{del} = Isp_{ODE} * \prod_{i=a}^k \eta_i}
The above efficiencies are usually separated into two groups; those that happen in the chamber and
those that happen in the nozzle.
Nozzle Losses
-------------
The Isp losses in the nozzle are more simple to model than the chamber losses and are usually
more consequential to the delivered Isp.
**Rocketisp** calculates overall nozzle efficiency from::
effDiv = Divergence Efficiency of Nozzle
effKin = Kinetic Efficiency of Nozzle
effBL = Boundary Layer Efficiency of Nozzle
effTP = Two Phase Efficiency of Nozzle
Chamber Losses
--------------
Chamber losses are usually more difficult to model than losses in the nozzle and
they are usually more subject to detailed design features.
That said, the combined chamber losses are usually small in production engines.
For example, 98% efficiency is a very typical value for all the combined chamber losses.
**Rocketisp** calculates overall chamber efficiency from::
effMix = Inter-Element Mixing Efficiency of Injector
effEm = Intra-Element Mixing Efficiency of Injector
effVap = Vaporization Efficiency of Injector
effHL = Heat Loss Efficiency of Chamber
effFFC = Fuel Film Cooling Efficiency of Chamber
effPulse = Pulsing Efficiency of Thruster
Mixing loss results from less than ideal mixing between the fuel and the oxidizer.
Vaporization loss results when the fuel and/or oxidizer do not completely vaporize to enable full combustion.
Those losses combined with any inability of combustion to reach
equilibrium chemistry are lumped to yield the injector losses,
sometimes referred to as combustion efficiency or energy release losses.
The most universally used parameter to characterize an injector is the characteristic velocity (:math:`c^*`).
`Rocket Propulsion Elements `_,
defines :math:`c^*`, pronounced "cee-star", as:
.. math::
\Large{c^* = P_c * A_t / \dot m}
The characteristic velocity (:math:`c^*`) is easily
determined from measured data of mass flow rate (:math:`\dot m`), chamber pressure (:math:`P_c`),
and throat area (:math:`A_t`). It relates to the efficiency of
the combustion and is essentially independent of nozzle characteristics.
:math:`c^*` is often
used as a measure of the injector excellence since fully mixed and
reacted propellants generally have higher :math:`c^*` than poorly-mixed and
partially-reacted propellants.
The most direct measurement of energy release efficiency (:math:`\eta_{ERE}`)
is the measured :math:`c^*` efficiency (:math:`\eta_{c^*}`) of the thrust chamber. If the
:math:`c^*` efficiency is used to back out energy release efficiency , then the effect of the nozzle discharge
coefficient :math:`C_D` must be taken into account.
In other words, for a perfect injector: :math:`\Large{\eta_{c^*_{PI}} = 1 / C_D}`
For a real injector: :math:`\Large{\eta_{c^*} = \eta_{ERE} / C_D}`
such that the energy release efficiency can be calculated as:
.. math::
\Large{\eta_{ERE} = \eta_{c^*} * C_D}
Note that, *In general, the use of empirical data is required to predict the energy release efficiency.*
-- `D.E. Coats `_