2.1. The Schmidt-Appleman Criterion for Fuel Cells
The physical process behind the formation of contrails is the mixing of two airmasses, one warm and moist (the exhaust gases), the other cold and drier (ambient air). This mixing may lead to a supersaturated state and then droplets can form, which may freeze if they get sufficiently cold. An aircraft exhaust plume mixes isobarically with the ambient air, such that its state-point in a temperature vs. water vapour partial pressure diagram follows a straight line, see, for instance, Figure 3 of [
3]. The condition for condensation of droplets, supersaturation with respect to liquid supercooled water, is met if the mixing line crosses the curve of saturation pressure (saturation with respect to liquid supercooled water). In such a case, a contrail is formed.
The mixing trajectory has the slope
where the indices “
p” and “
a” mean “plume” and “ambient”. The denumerator is related to mass conservation of water, where
is the partial pressure of water vapour. The denominator expresses energy conservation, and
are the static temperatures.
Engineers may skip the rest of this section and jump directly to
Section 2.5, because the calculations that follow are in a certain sense academic; they consider the unaltered exhaust from a fuel cell. This leads to a form of the G-factor that is similar to its traditional form and it yields physical insight. However, technical manipulations of the exhaust, like heat exchanger and condenser for water vapour differ from case to case; they cannot be covered by a common formula of the traditional form. In this case, an engineer can still use the basic definition of
G above if only the parameters at exhaust exit (
) and their atmospheric counterparts are known. For the unmodified exhaust gases these parameters can be computed with basic physical considerations, which will be done in the following.
Starting with mass conservation, the flows of matter into and out of the fuel cell (FC) are considered in terms of molar fluxes (units mole per second, mole s
). The fed-in air and hydrogen fuel can be considered ideal gases, such that pressure ratios equal the number or molar ratios. Let the input rate of the FC consist of
moles s
of H
and
moles s
of air. The air consists mainly of nitrogen, N
, oxygen, O
, with a molar fraction
, and ambient water vapour with a molar fraction
. Hydrogen and oxygen react to water vapour in the following reaction:
which yields one mole of water per one mole of hydrogen and half a mole of oxygen. Thus the output rate at the FC exit in moles s
is lower than that at the input, since
moles are converted into one mole in the reaction.
If all hydrogen molecules are consumed in the reaction, the reaction rate equals the rate at which H
is fed in,
. The rate at which O
from the ambient air is consumed is half as large,
, and the rate at which water molecules are added to the exhaust gases equals the rate at which H
is fed in,
. The FC output thus has a total gas flow of
moles s
, since the input flux of hydrogen equals the output flux of water vapour produced in the reaction (
3). The water vapour fraction of the output consists of the produced water plus the water from the input air, which flows in and out with the rate
. The molar ratio of water vapour at the FC exit is thus
where the assumption is made that the exhaust gas obtains immediately the ambient pressure (which is a reasonable assumption since any pressure contrasts are relaxed at the speed of sound; see the chapter on elastic waves in any physics textbook). After a short calculation, one finds
Since
one may neglect this contribution to find the desired expression for
as:
For the consideration of energy conservation, it will be useful to begin with a few thermodynamic considerations of the involved quantities. Assuming that the FC operates such that the water is produced as steam, not as liquid water, implies that one needs to choose the so-called lower heating values of the thermodynamic variables. At standard conditions (i.e., pressure 1 bar or
Pa and temperature 25
C or
K), the formation enthalpy of water vapour from its constituent elements is
kJ mol
and the corresponding free enthalpy (or Gibbs free enthalpy) is
kJ mol
(Table A 1.1 of [
18]). The free enthalpy (instead of the free energy) is relevant here as the reaction occurs at constant pressure and temperature in the FC. It is the maximum non-expansive work, that is, the maximum electric work, the FC can perform. In the reaction (
3) the entropy of the reactants is higher than the entropy of the product by
J
K
at standard conditions. Thus it is necessary that a part of the reaction enthalpy is used to heat the exhaust gases in order to guarantee that the overall entropy change is not negative. This “entropy tax” is expressed by the difference between
and
. It is heat that increases the temperature of the exhaust gases.
Generally, an FC will work under conditions that are different from the standard conditions and the thermodynamic quantities
and
will differ from
and
. The differences are however small (see for instance [
18], Chapter 2.7, [
19], Chapter 2.4, or [
20], Chapter 4.4.2). Thus only small errors incur in numerical calculations if one uses the enthalpy and free enthalpy in standard conditions. In order to present the results for general conditions,
and
(without the superscript 0) are used in the following derivations.
The fuel cell is characterized by its electromotive force
, which depends on the free enthalpy of the reaction (
3):
where Faraday’s constant
F is the number of negative elementary charges per mole of electrons, expressed in C mol
:
C mol
.
F is multiplied by 2 since 2 electrons flow from the anode to the cathode [
18,
19,
20]. The electromotive force of the hydrogen fuel cell is
V (if the water is released in gaseous state). The electromotive force is an ideal voltage of the cell which requires an infinite resistance between the electrodes, thus zero current. Once the cell operates and produces work, its current is given by the rate of the reaction
its actual voltage
U is less than
, typically in the range
to 1 V (see Figure 4.3 of [
20]), and thus the electric power that can be used to drive an engine is
. The remaining power
contributes to heat the exhaust gases, which adds to the “entropy tax”,
mentioned above.
The energy balance of the working FC is thus
That is, the reaction enthalpy of the oxyhydrogen reaction is split into electric work and heating of the exhaust due to a lower than ideal voltage of the FC and due to the entropy tax.
It will turn out useful to introduce two efficiencies, a basic efficiency
and an electric efficiency
With these efficiencies the energy balance reads
The power that heats the exhaust is thus . The product can be conceived to be the efficiency of the fuel cell.
The heat capacity of the input air is
J mol
K
and that of gaseous water is
J mol
K
[
18]. The waste heat and entropy tax warms the air as follows:
Solving the energy balance for
yields
In principle, the ingredients to calculate
are now at hand, but the presence of different heat capacities would lead to an unwieldy expression. One can achieve a more concise expression defining a mean heat capacity of the exhaust gases as follows:
It is clear that this heat capacity depends on the air/fuel ratio, but it allows now to write the expression for
in a familiar form:
a quite simple result with a form that resembles the one of the G-factor for traditional jet engines. While the traditional form explicitly shows a factor for the amount of water vapour emission, namely its emission index, such a factor does not appear here. The simple reason for this is that this factor is unity: one mole water vapour is emitted for each mole of hydrogen fuel. The only essential difference to the traditional formula is that the heat capacity is not a constant. Since the heat capacities of air and water vapour are moderately different, the mixture of the gases (i.e., the air/fuel ratio) has a moderate effect on the resulting G-factor (see below).
2.2. Which Efficiency Is Relevant?
The basic efficiency is almost at standard conditions. It is larger than the value of which refers to the upper heating values of the enthalpy and free enthalpy, while the refers to the lower heating values that we use here.
Not all of the hydrogen fed into the FC is actually burnt; there may be a fraction of unburnt H. In principle, this loss fraction can be minimised if the unburnt H is fed back into the FC input pipe. One may assume therefore that this efficiency loss can by minimised to a negligible degree.
Thus, for the fuel cell itself the relevant efficiency is the product of the basic efficiency and the electric efficiency, . The electric efficiency varies with the current. If there is no current, U is close to its theoretical maximum , but once the current flows, U decreases with increasing I. The U vs. I dependence is called the operational characteristic of the fuel cell. It depends on the type and construction of the fuel cell. Ideally there is a range of currents, where the voltage and thus varies little.
For a stand-alone fuel cell running idle this consideration would suffice. However, for FCs driving airplanes there are external boundary conditions. The airplane needs a certain power
to overcome drag and friction.
F is its thrust and
V its speed. Only a fraction of the original power
produced by the fuel cell (or stack of fuel cells) will be converted eventually into the mechanical power
, thus
, where
stands for an “installation” efficiency, which obviously depends on how the power produced by the FC is transferred to the propeller or fan and from there to the momentum change of the air around the propeller or fan. It can itself be a product of a number of partial efficiencies. It’s exact value is certainly important to know for the estimation of how much H
needs to be carried on a given flight, but it is not relevant to know for the contrail formation theory, since the heat produced by less-than-ideal installation components will not heat the exhaust gases of the FC. It will eventually heat the air inside the wake vortex of the plane, but long after condensation of the exhaust water vapour has been occurred. This argument is borrowed from the first paper that introduced the factor
into the Schmidt-Appleman criterion [
21]. In that paper, it was argued that only the heat that directly warms the exhaust gases is relevant for contrail formation. Eventually all the work that is performed by the engines is converted to heat as well, but that warms the air in the aircraft’s wake only long after contrail formation has happened. Thus, only
has to appear in the equation for
. This differs from the situation of jet engines, where the overall propulsion efficiency is used in the criterion. The reason for this is obviously that in a jet engine power generation and the major fraction of power consumption occur at the same place, which is an important difference to the case of FCs.
2.4. Simple Applications
In order to get an impression of the magnitude of
in various situations, a couple of simple applications of Equation (
16) are presented assuming stoichiometric conditions for which
J mol
K
. The effect of a higher air/fuel ratio does not exceed a factor
, as demonstrated above, so that the application can be restricted to stoichiometric conditions. Further one needs assumptions for ambient temperature and pressure. For this an average atmospheric profile
is used, given by the US-Standard Atmosphere [
22]. Two versions of fuel cell operation are generally assumed, one at
V and one at
V.
In 10 km altitude, K and hPa. This yields a G-factor of Pa K for a FC operated at V and a lower value of Pa K for a FC at V. The same fuel cells working at ground level: K and hPa yield Pa K for 1 and V, respectively. These G-factors are considerably larger than those of traditional jet engines, for which a value of 2 Pa K is typical.
From the simple examples it is clear that the G-factors of FCs are much larger than those of jet engines, even when the latter are driven with liquid hydrogen. This implies a large potential for condensation of water vapour in the expanding exhaust plume. A large air/fuel ratio lowers the G-factor only little. One consequence of this is that the threshold temperature for contrail formation is much higher than for jet engines, such that contrails may form even at the ground.
2.6. Atmospheric Implications
Figure 2 shows
for
V (solid red) and
V (solid blue) together with the temperature profile of the US standard atmosphere (black) [
22]. The stoichiometric air/fuel ratio of
has been assumed. Lines for higher air/fuel ratios would appear very similar. In all cases, contrail formation is possible at much lower altitudes than for jet engines. The minimum altitude necessary for contrail formation decreases with increasing FC voltage, because higher voltage implies higher efficiency to drive the engine and hence a colder exhaust follows, which is more susceptible to contrail formation. With the higher voltage, contrail formation is possible even at ground if the ambient temperature falls to values below about 5
C.
There is another threshold,
(dashed lines), at which contrail formation occurs even in absolutely dry air. This is given as
That is, at contrail formation occurs always, at never, and between these boundaries it depends on the relative humidity (see below). The lower thresholds, that is, below which contrails are always formed, are slightly above C for the 1 V FC on ground. Thus, contrail formation at the ground in winter can occur frequently for such a FC. A lower voltage leads to a substantial reduction of these thresholds at the expense of efficiency.
Figure 2 also shows that under conditions of the US-Standard Atmosphere contrails are always formed above 380 hPa (about
km). The mean temperature at that level is still a few degrees above the temperature where minute water droplets freeze spontaneously. Droplets formed under these conditions will likely evaporate. However, the limit for spontaneous freezing (ca.
C, see [
23], sect. 4a) is not far away. At the usual cruise levels, fuel cells will generally produce contrails. Below
C the droplets will freeze and the contrail may be persistent then, if the ambient relative humidity exceeds ice saturation.
For temperatures between the upper and lower thresholds contrail formation requires a minimum relative humidity of the ambient air:
Figure 3 shows for each pressure level the relative humidity required to produce a contrail (if the temperature actually is that of the assumed US standard atmosphere from
Figure 2). Furthermore it displays for the same assumed temperature profile the relative humidity at which ice saturation is reached, which is the minimum required for contrail persistence. Again the curves are displayed for voltages of 1 V (red) and
V (blue).
2.7. Comparison to Hydrogen Burning in a Gas Turbine
The burning of gases in a gas turbine delivers expansive work that finally drives the turbine. The reaction enthalpy of the combustion reaction equals the expansive work done by the reacting gases under conditions of constant pressure and temperature. (It may sound surprising to invoke the condition of constant temperature, but once the engine is switched on the temperature rises until it attains a constant high value when the engine runs continuously.) Hence the relevant energy for a gas turbine is
. The temperature in a combustion chamber is much higher than 100
C, such that the product of burning hydrogen is water vapour, and the relevant enthalpy is the lower heating value that was already applied above,
kJ mol
. The contrail factor for a hydrogen gas turbine is thus
where all quantities are given per mole. A more familiar form appears when the quantities are given per kg with
J (kg K)
and
MJ kg
,
and with the emission index of water vapour
kg(water vapour) per kg(liquid hydrogen):
The G-factor for burning liquid hydrogen in a gas turbine is times larger than that for kerosene, if the overall propulsion efficiencies are equal. The thermal efficiency of the relevant Brayton cycle depends on the heat capacity ratio, , of the exhaust gas which in turn depends on the composition of the gas and thus on the amount of water vapour in the exhaust. However, the values of differ little for water vapour (≈1.33) and air (≈1.4) and the water vapour fraction in the exhaust is small both for kerosene and for LH2. Thus the difference between the thermal efficiencies are small for the two types of fuel. Other components of the overall efficiency do not depend on the fuel.
Figure 4 shows a comparison of the contrail formation conditions for LH2 gas turbines with the hydrogen fuel cell. An overall propulsion efficiency of
has been assumed for the gas turbine. It turns out that
for the LH2 gas turbine is lower then even
for the fuel cell operated at
V, but the contrail formation conditions for the LH2 gas turbine are very similar to those of the FC at
V, as a comparison with
Figure 2 and
Figure 3 shows.