Power Density of Hydrogen Magnetohydrodynamic (H2MHD) Generators at Different Pressures, Seed Types, Seed Levels, and Oxidizers

1. Introduction

1.1. Background

Hydrogen is the lightest element in the periodic table but it is the most plentiful element in the universe [1,2] (while oxygen is the most abundant element in the earth’s crust by mass [3,4]). Hydrogen is an important building block in several chemical processes and power systems [5,6], as listed in Table 1.

Hydrogen facilitates the transition to a low-emission economy [55,56,57], taking advantage of the absence of released carbon dioxide (CO₂) due to the consumption of hydrogen [58,59,60,61]. The production of green hydrogen is based on water electrolysis powered by electricity generated from renewable energy sources (such as solar energy [62,63,64] and wind energy [65,66,67,68]), and thus it does not involve the release of greenhouse gases (GHG) [69,70,71], making it an attractive way to mitigate GHG emissions without resorting to carbon capture technologies as needed in the case of blue hydrogen [72,73,74,75] that is originating from fossil fuels.

In an earlier study, we proposed a novel use for hydrogen (either green hydrogen or blue hydrogen) in open-cycle magnetohydrodynamic (OCMHD) power generators [76], in which direct power extraction (DPE) is realized with the aid of the Lorentz force. In this concept, hot gases formed as combustion products of hydrogen are seeded with a small amount of an alkali metal (either cesium “Cs” or potassium “K”) [77,78,79], which allows for an appreciable level of thermal ionization leading to a suitable electric conductivity of the weakly-ionized plasma gas, which is accelerated through a linear channel [80,81,82]. The high-speed plasma loses part of its energy in reaction to an applied magnetic field, where an induced electric field allows collecting an induced electric current to power an external electric load [83,84,85]. This concept of hydrogen-powered magnetohydrodynamic power generation, which we refer to here as (H2MHD), can be combined with a classical power cycle [86,87,88] for better utilization of the input heat released from the zero-emission hydrogen combustion through a dual-cycle power plant, comprised of a higher-temperature open-cycle MHD (OCMHD) generator and either a lower-temperature closed-cycle Rankine system or a lower-temperature open-cycle Brayton system [89,90,91]. A triple combined cycle that incorporates the open-cycle MHD (OCMHD) generator as a topping layer (high temperature) with an intermediate layer having an open-cycle Brayton gas turbine (intermediate temperature) and then a layer of a closed-cycle steam turbine (low temperature) is also possible, to recover part of the rejected heat containing in the still-hot plasma gas leaving the MHD generator [92,93,94]. Introducing an even lower-temperature refrigerant-based organic Rankine cycle [95,96,97] as an additional terminal downstream level (very low temperature) was proposed, yielding a quadruple combined cycle [98].

Figure 1 depicts how the proposed hydrogen magnetohydrodynamic (H2MHD) topping cycle operates. The plasma velocity vector is designated by the symbol ($\stackrel{\rightharpoonup }{u}$), the applied magnetic field vector (applied magnetic-field flux density vector) is designated by the symbol ($\stackrel{\rightharpoonup }{B}$), and the electric current-density vector is designated by the symbol ($\stackrel{\rightharpoonup }{J}$). Hydrogen is combusted with an oxidizer, which is either air (in the case of the conventional air-combustion mode) or oxygen (in the case of the intensified oxy-combustion mode). Oxy-combustion (the combustion of HHO or Brown’s gas) offers the advantage of elevating the combustion temperature due to the absence of diluent nitrogen [99,100,101]. A small amount of solid alkali salt, such as potassium carbonate (K₂CO₃) [102] is added in the form of an aqueous solution that is introduced as a fine spray by atomization, to serve as the source of the vapor alkali metal after combustion [103,104]. Either potassium element or cesium element is readily vaporized well below the combustion temperature. Potassium (K) vaporizes at 760 °C only [105], while cesium (Cs) vaporizes at a lower boiling point of 671 °C [106]. The hot-seeded hydrogen-combustion products are accelerated to supersonic speeds using a de Laval nozzle (convergent-divergent nozzle or CD nozzle) [107,108,109]. The term “channel” or “linear channel” here refers to the portion of the divergent section of the convergent-divergent nozzle where an external magnetic field is applied perpendicular to the plasma’s direction of bulk motion, and an induced electric emf (electromotive force) is produced; thus, an electric current-density can be collected through electrodes (positive cathode and negative anode) to be delivered to an external load. Like solar photovoltaic (PV) modules [110,111] and thermoelectric generators (TEG) [112,113], MHD generators operate on the principle of direct power conversion (from input energy to electricity in a single step), without an intermediate turbogenerator [114,115], and without reciprocating or rotary parts [116,117]. This is a remarkable feature in MHD generators.

The power density (or volumetric power density) is a performance metric for power generation units, describing the power output per unit volume occupied. This metric quantifies the effective utilization of space. For MHD generators, the intensive high-temperature and high rate of power extraction from the plasma gives this technology a unique favorable characteristic when compared to other power technologies. For example, the Sakhalin MHD pulsed power generator experimentally demonstrated a temporary power density of about 90.7 MW/m³ (computed as an electric power output of 510 MW divided by an MHD linear channel volume of 5.625 m³) [118,119,120]. Another shock-tunnel MHD generator experimentally demonstrated a high power density of 140 MW/m³ [121]. These power densities are much higher than typical values of automotive engines (around 15 MW/m³, the volume here refers to the displacement volume [122,123]), permanent magnet (PM) synchronous motors (SM) (around 7.5 MW/m³ [124]), and triboelectric nanogenerator (TENG) energy harvesting units (up to about 5 MW/m³ [125]), and microbial fuel cells (MFCs) (between 1 × 10^–6 MW/m³ and 5 × 10^–5 MW/m³ [126,127,128]). Numerical two-dimensional modeling of a hypothetical MHD channel with methane (CH₄) combustion suggested a lower power density of 17.9 MW/m³ (computed through dividing an electric power of 31.4 MW by a volume of 1.75 m³) [129]. This is not remarkably high, but it is still encouraging.

There is another metric for expressing the effective utilization of a power system, which is the specific power of gravimetric power density (W/kg). However, this is more challenging to estimate compared to the volumetric power density; because estimating the specific power requires more detailed and precise knowledge about the particular power system and its components, including the minor auxiliary ones if they largely affect the weight [130,131]. The specific power (or the power-to-weight ratio [132,133]) is of particular importance in aeronautical and space propulsion systems [134,135,136]; which are outside the scope of this study that is concerned with terrestrial power systems. Thus, the volumetric power density is adopted here. The specific energy or the gravimetric energy density metric (Wh/kg) [137,138], and the energy density or volumetric energy density metric (Wh/L) [139,140] pertain more to energy storage units (such as batteries), and thus are not used in our study.

1.2. Goal of the Study

Because the proposed H2MHD concept is a new technology for utilizing hydrogen, there is a need for extensive investigations regarding the technological feasibility of this technology. This study aims to partly fill this huge knowledge gap through estimating the volumetric power density of an ideal H2MHD channel, at a representative value for the applied magnetic field (the applied magnetic-field flux density), a reasonable supersonic Mach number (the ratio of the plasma speed to the speed of sound), and a suitable temperature representative of hydrogen combustion.

The performance of the volumetric power density is estimated for a wide range of operational conditions, with different total pressures, different levels of seeding, two types of the seed metal, and two types of the oxidizer.

The results of this study provide several insights into the use of OCMHD in general, rather than H2MHD solely because some of the generic findings reported here can actually be applicable to other OCMHD systems with fuels other than hydrogen. Thus, this study is considered useful for MHD power generators in general, and for hydrogen use in these power generators in particular.

2. Research Method

2.1. Assumptions

The current study is centered around zero-dimensional computational modeling (point modeling or system-level modeling) of a hypothetical supersonic H2MHD (hydrogen magnetohydrodynamic) channel. Therefore, no exact geometry is defined; instead, the analysis assumes uniform local properties everywhere. Despite losing boundary effects and temporal transience [141,142,143,144,145] that can be captured using computational fluid dynamics (CFD) methods [146,147,148], our simplifying assumption is important to allow having a closed-form expression for the power density, which is conveniently applied to estimate theoretical upper bounds. In addition, our assumption is useful in generalizing the outcomes of this study, making them broadly valid, rather than being limited to a specific H2MHD channel.

In addition, due to the large number of variables involved in the H2MHD problem, we need to fix some variables at reasonable values while changing others to explore their impact. Therefore, three quantities are fixed in the current study, and these fixed quantities are (1) the temperature, (2) the applied magnetic-field flux density, and (3) the Mach number. These fixed parameters are set to the values listed in Table 2. The significance of these fixed parameters lies also in the fact they are the same values assumed in our previous work mentioned earlier in which we proposed the H2MHD concept (but the term “M2HHD” is coined here for the first time).

2.2. Mathematical Equations

For a plasma speed ($u$), an applied magnetic-field flux density magnitude ($B$), and an electric conductivity (σ); the theoretical upper limit for the volumetric electric power density ($P_V$) that can be extracted from the moving plasma gas is [159]

$$P_V={\displaystyle \frac{1}{4}}\sigma u^2{B}^2$$

(1)

This power density is restricted to ideal conditions, where the external load is matched (optimized) such that its electric resistance is equal to the equivalent internal electric resistance of the MHD generator [160,161,162].

The plasma speed ($u$) can be expressed in terms of the speed of sound ($a$) and the Mach number ($\mathrm{M}$) as [163]

$$u=\mathrm{M} a$$

(2)

Using Equation (2) in Equation (1) allows the expressing of the power density ($P_V$) in a different form as

$$P_V=0.25 \sigma {\mathrm{M}}^2 a^2 B^2$$

(3)

The plasma electric conductivity ($\sigma$) in our study is computed according to an algorithm we described and verified in detail in a previous study [164], which is based on the robust Frost’s two-conductivity approximation that takes into account the scattering of electrons by neutrals as well as the scattering of electrons by charged particles (ions and other electrons) [165], and thus interpolates the electric conductivity for the partially-ionized plasma in OCMHD channels between the two extreme values for which an exact solution is known (namely, the extreme case of very weakly-ionized plasma, and the extreme case of fully-ionized plasma) [166,167,168]. In our algorithm, the plasma’s electric conductivity depends on the total static pressure ($p$) (the term “total” here means the sum of partial pressures of involved gaseous species, not the stagnation zero-speed pressure, and not the gauge pressure [169]), the absolute temperature ($T$), and the pre-ionization mole fractions ($X_i$, where $i$ is an index referring to each involved gaseous species in the gaseous mixture formed after adding the seed alkali metal vapor to the hydrogen combustion products). Since the temperature is fixed in the current study at 2,300 K; practically, the electric conductivity depends only on the total static pressure ($p$) and the mole fractions ($X_i$). Therefore, we can symbolically represent the electric conductivity ($\sigma$) as ($\sigma \left(p,X_i\right)$) to emphasize its functional dependence on these plasma properties.

The speed of sound in the pre-ionization plasma ($a$) is modeled here using the ideal gas approximation, making it a function of the absolute temperature ($T$) and the mole fractions ($X_i$) only [170,171,172]. This speed of sound can be expressed as

$$a=\sqrt{{\displaystyle \frac{C_{p,mix}}{C_{p,mix}-R_{mix}}} R_{mix} T}$$

(4)

where ($C_{p,mix}$) is the specific heat capacity at constant pressure of the gas mixture after seeding but before ionization, and ($R_{mix}$) is the specific gas constant for this gas mixture. The pre-ionization mixture’s specific heat capacity ($C_{p,mix}$) is computed from the specific heat capacities of the individual gaseous species ($C_{p,i}$ for the species with the index $i$) that make up the pre-ionization gas mixture, and these individual specific heat capacities are described according to the updated nine-coefficient edition of NASA (United States National Aeronautics and Space Administration) nonlinear fitting functions for computing thermodynamic properties [173,174,175,176]. Therefore, we can analytically express the speed of sound ($a$) as ($a\left(T,X_i\right)$) to emphasize its functional dependence on the plasma properties ($T$) and ($X_i$).

The nondimensional quantity ($C_{p,mix}/\left(C_{p,mix}-R_{mix}\right)$) in Equation (4) is called the specific heat ratio or the adiabatic index, and it is assigned here the symbol (${\gamma }_{mix}$). Therefore, the previous equation for the speed of sound ($a$) can be rewritten as [177]

$$a=\sqrt{{\gamma }_{mix} R_{mix} T}$$

(5)

The gas constant ($R_i$) for each gaseous species having the index $i$ is obtained using

$$R_i={\displaystyle \frac{\overline{R}}{M_i}}$$

(6)

where ($\overline{R}$) is the universal (molar) gas constant, and ($M_i$) is the molecular weight of the species with an index ($i$). In the current study, we have [178,179,180]

$$\overline{R}=8,314.462618{\displaystyle \frac{\mathrm{J}}{\mathrm{kmol}.\mathrm{K}}}$$

(7)

The molecular weights of the gaseous species that appear in the current study are listed in Table 3. These were taken values from the NIST (United States National Institute of Standards and Technology) individual gaseous species Chemistry WebBook database [181,182,183].

The molecular weight of the gas mixture ($M_{mix}$) is computed as [188]

$$M_{mix}={\displaystyle \sum }_{i=1}^{n_s}{M}_i X_i$$

(8)

where ($n_s$) is the number of gaseous species in the pre-ionization gas mixture. We have ($n_s=2$) when the oxidizer is pure oxygen (oxy-hydrogen combustion), while ($n_s=3$) when the oxidizer is air (air-hydrogen combustion). In the case of oxy-combustion, the pre-ionization mixture is treated here as a mixture of water vapor (H₂O) and the vaporized alkali metal (either Cs or K). In the case of air-hydrogen combustion, the pre-ionization mixture is treated here as a mixture of water vapor (H₂O), excess nitrogen (N₂), and the vaporized alkali metal (either Cs or K).

After obtaining ($M_{mix}$), the specific gas constant of the pre-ionization gas mixture ($R_{mix}$) can be computed as

$$R_{mix}={\displaystyle \frac{\overline{R}}{M_{mix}}}$$

(9)

Then, the mass fraction of each individual species ($Y_i$) can be computed as [189]

$$Y_i={\displaystyle \frac{X_i M_i}{M_{mix}}}$$

(10)

Then, the specific heat capacity of the mixture is computed from the specific heat capacities of the individual species through mass-weighted averaging as

$$C_{p,mix}={\displaystyle \sum }_{i=1}^{n_s}{C}_{p,i} Y_i$$

(11)

The NASA nonlinear fitting function for estimating ($C_{p,i}$) has the following general structure, with seven terms:

$$C_{p,i}=R_i\left(b_{1,i} T^{-2}+b_{2,i} T^{-1}+b_{3,i}+b_{4,i} T+b_{5,i} T^2+b_{6,i} T^3+b_{7,i} T^4\right)$$

(12)

It should be noted that the fitting coefficients ($b_{1,i}$ to b_7,i) in the previous equation correspond to the higher-temperature range (from 1,000 K to 6,000 K); where there is another set of fitting coefficients for a lower-temperature range (from 200 K to 1,000 K), but these are irrelevant in our study (too cold to be of interest).

The combustion of hydrogen with air is approximated as a complete stoichiometric reaction of molecular hydrogen (H₂) with a mixture having mole fractions of 79% for molecular nitrogen (N₂) and 21% for molecular oxygen (O₂) [190,191,192]. Therefore, the combustion products are idealized as having only water vapor (H₂O) and excess non-reacting molecular nitrogen (N₂) according to the following reaction equation [193,194]:

$$2{\mathrm{H}}_2+{\mathrm{O}}_2+3.762{\mathrm{N}}_2\to 2{\mathrm{H}}_2\mathrm{O}+3.762{\mathrm{N}}_2$$

(13)

Vapor seed metal (either cesium or potassium) is then added to the hydrogen combustion products with a predetermined target mole fraction. This seeded mixture is then analyzed under thermal equilibrium to obtain the electrochemical properties of the post-ionization plasma (such as its electric conductivity and mole fractions of ions).

The combustion of hydrogen with oxygen is approximated as a complete stoichiometric reaction of molecular hydrogen (H₂) with molecular oxygen (O₂), resulting in pure water vapor (H₂O) as the single combustion product. This is described by the following reaction equation:

$$2{\mathrm{H}}_2+{\mathrm{O}}_2\to 2{\mathrm{H}}_2\mathrm{O}$$

(14)

Vapor seed metal (either cesium or potassium) is then added to this combustion product water vapor with a predetermined target mole fraction. This seeded mixture is then analyzed under thermal equilibrium to obtain the electrochemical properties of the post-ionization plasma.

We conclude the current mathematical description of the modeling methodology by presenting a customized version of the expression governing the ideal volumetric power density in Equation (3), by setting the Mach number to (M = 2), and the magnitude of the applied magnetic-field flux density to ($B=5 \mathrm{T}$). This gives the following special expression pertaining to our study:

$$P_V\left[{\displaystyle \frac{\mathrm{W}}{{\mathrm{m}}^3}}\right]=25 \sigma \left[{\displaystyle \frac{\mathrm{S}}{\mathrm{m}}}\right] {\left(a\left[{\displaystyle \frac{\mathrm{m}}{\mathrm{s}}}\right]\right)}^2$$

(15)

To conveniently express the values of the power density ($P_V$), the larger unit of (MW/m³) is adopted rather than the basic SI unit of (W/m³), thus the previous equation is modified to

$$P_V\left[{\displaystyle \frac{\mathrm{MW}}{{\mathrm{m}}^3}}\right]={\displaystyle \frac{25}{{10}^6}} \sigma \left[{\displaystyle \frac{\mathrm{S}}{\mathrm{m}}}\right] {\left(a\left[{\displaystyle \frac{\mathrm{m}}{\mathrm{s}}}\right]\right)}^2$$

(16)

While the speed of sound ($a$) in the above equation should correspond to the post-ionization plasma condition (taking into account the small changes in the gaseous mixture because of the partial ionization of some alkali seed atoms), we use the pre-ionization value that corresponds to an all-neutral (no ions, no electrons) seeded gas mixture. This establishes consistency with the reported total pressure values as a controllable variable, where these pressures correspond to the pre-ionization condition. This matter is practically negligible because only a small portion of the seeded alkali metal can be thermally ionized at the fixed temperature of 2,300 K; causing the mole fraction of either the alkali ions or the free electrons to be below 0.1% in all cases we cover, and even much less than this in many of the analyzed cases.

2.3. Varied Conditions

In Table 5, we list the quantities to be varied in the current study while estimating the ideal volumetric power density ($P_V$) for a magnetohydrodynamic (MHD) generator channel, as well as the values assigned to these varied quantities. These quantities are

(1)

the total absolute pressure ($p$) after the vaporous seed alkali metal is mixed with the hydrogen combustion gases. The unit is the standard atmosphere (1 atm = 101,325 Pa = 1.01325 bar). Nine values are assigned to this continuous variable, sampling it between a very low value of 1/16 atm (6332.812 Pa or 6.332812 kPa) to a high value of 16 atm (1.6212 MPa or 1,621.2 kPa). These values allow exploring the impact of this variable over a wide range. The nine selected pressures form a geometric sequence with a multiplicative progression constant of two.

(2)

the mole fraction ($X_{seed}$) of the vaporous seed alkali metal after it is mixed with the hydrogen combustion gases. This fraction is expressed as a percentage. Five values are assigned to this continuous variable, sampling it between a very low value of 0.0625% or (1/16)% to a high value of 16%. These values allow exploring the impact of this variable over a wide range. The five selected seed mole fractions form a geometric sequence with a multiplicative progression constant of four.

(3)

The type of the vaporous seed alkali metal to be added to the hydrogen combustion gases. This is a binary variable, and it can be either cesium (Cs) or potassium (K). These two particular chemical elements are chosen here as the possible seed metal that serves as the source of electrons due to their easy thermal ionization compared to the combustion products gases (H₂O or N₂). The ionization energy of cesium (Cs, atomic number 55; located in the alkali metal group IA or group 1 and period 6 in the periodic table of chemical elements) is 3.893 eV/particle (375.6 kJ/mol), and this is the lowest value among all chemical elements in the periodic table [195,196,197]. We point out here that 1 eV/particle (electronvolt per particle) is equivalent to 96.49 kJ/mol (kilojoule per mole) [198,199,200]. The ionization energy of potassium (K, atomic number 19; located in the alkali metal group IA or group 1 and period 4 in the periodic table of chemical elements) is 4.34 eV/particle (419 kJ/mol), and this is the fourth lowest value among all chemical elements in the periodic table [201,202,203]. Although francium (Fr, atomic number 87; located in the alkali metal group IA or group 1 and period 7 in the periodic table) has a lower ionization energy of (3.9 eV/particle or 376 kJ/mol) than potassium [204,205]; francium is radioactive and very stable (has no practical uses) due to its small half-life of only 22 minutes [206,207,208]. Although rubidium (Rb, atomic number 37; located in the alkali metal group IA or group 1 and period 5 in the periodic table) has a lower ionization energy of (4.177 eV/particle or 403.0 kJ/mol) than potassium [209,210]; rubidium is a rare element, and its use is typically limited to scientific research [211,212,213]. The ionization energy of molecular nitrogen (N₂) is very large, being 15.6 eV/particle (1,505 kJ/mol) approximately [214,215]. The ionization energy of atomic nitrogen (N) is also large, being 14.5 eV/particle (1,399 kJ/mol) approximately [216,217]. The dissociation energy of molecular nitrogen (N₂) into two atoms is also much lower (9.76 eV/mole or 942 kJ/mol) [218,219], but it is still higher than the ionization energy of either potassium or cesium. The ionization energy of the water molecules (H₂O) is very large, being 12.6 eV/particle (1,216 kJ/mol) approximately [220]. The dissociation energy of the water molecules (H₂O) into a hydrogen atm (H) atm and a hydroxyl radical (OH) radical is approximately 5.16 eV/particle (498 kJ/mol) [221], making this dissociation less probable than the ionization of either potassium or cesium.

(4)

The type of oxidizer to be used in the combustion of hydrogen. This is a binary variable, and it can be either air or pure oxygen (O₂).

Therefore, there is a total of 180 conditions (data points) to be computationally represented in the current study (because there are nine pressures, five seed mole fractions, two seed types, and two oxidizer types).

Table 5. Varied quantities in the current study.

Counter	Varied Quantity
	$\mathit{p}$ [atm]	${\mathit{X}}_{\mathit{s}\mathit{e}\mathit{e}\mathit{d}}$ [%]	Seed metal	Oxidizer
1.	0.0625	0.0625%, (1/16)%	Cesium (Cs)	Air (79% N2, 21% O2)
2.	0.125	0.25%	Potassium (K)	Oxygen (100% O2)
3.	0.25	1%
4.	0.5	4%
5.	1	16%
6.	2
7.	4
8.	8
9.	16

Table 5. Varied quantities in the current study.

Counter	Varied Quantity
	$\mathit{p}$ [atm]	${\mathit{X}}_{\mathit{s}\mathit{e}\mathit{e}\mathit{d}}$ [%]	Seed metal	Oxidizer
1.	0.0625	0.0625%, (1/16)%	Cesium (Cs)	Air (79% N2, 21% O2)
2.	0.125	0.25%	Potassium (K)	Oxygen (100% O2)
3.	0.25	1%
4.	0.5	4%
5.	1	16%
6.	2
7.	4
8.	8
9.	16

2.4. Power Density Criterion

In our study, it can be sufficient to estimate the volumetric electric power density output under the various operational conditions for the H2MHD (hydrogen magnetohydrodynamic) concept, as a primary goal for our study and as it major contribution. However, the value of our study is boosted when a criterion for this power density is also recommended, to help in judging the feasibility of the H2MHD concept and the preferred set of conditions suitable for it.

To justify such a criterion (a power density threshold), we consider that a mid-size power plant has a power capacity on the order of 600 MW [222,223,224,225]. If an H2MHD power density of 50 MW/m³ is possible, then this 600 MW power capacity can be accomplished with a compact space of only 12 m³ in the case of a properly operated H2MHD generator. This small volume can be formed for example as a divergent channel with a width of 1 m, a length of 6 m, and a mean height of 2 m [226,227,228,229]. To account for the assumptions made in the current study, and the losses ignored in our system-level approach, we set our threshold to twice this value, making it 100 MW/m³. This means that 100 MW/m³ in our analysis may actually be reduced to half of this value in reality. To show the gigantic amount of power delineated by this criterion, it is helpful to compare this threshold (100 MW/m³) with the typical area power density of photovoltaic (PV) solar panels under good sunshine, which is only about 200 W/m² (0.0002 MW/m²) [230,231,232].

We emphasize that this suggested criterion of 100 MW/m³ here is not a precise value, but involves some arbitrariness. However, it is considered a reasonable rough guide for assessing the suitability of H2MHD power generation.

3. Results

This results section is divided into four subsections, which are mapped with the four possibilities arising from the two binary variables (namely, two alkali seed types and two oxidizer types). For each case of these four conditions, we present first the variation of the plasma electric conductivity ($\sigma$) after thermal equilibrium, as a function of the total pre-ionization absolute pressure ($p$) with the mole fraction of the neutral alkali seed metal ($X_{seed}$) being a parameter. Thus, the variation of ($\sigma$) is visualized versus ($p$) as five curves (one curve per $X_{seed}$), and this visualization is presented four times (for two seed types and two oxidizer types). Similarly, the estimated ideal volumetric power density ($P_V$) is visualized as a function of ($p$) four times in the four coming subsections, at five values of ($X_{seed}$) in each ($P_V$) plot.

To aid in interpreting the electric conductivity, it is voluntarily expressed as a multiple of the typical seawater electric conductivity, which is assigned here a value of 5 S/m [233,234,235]. The secondary (right) vertical axis in the electric conductivity ($\sigma$) plots is utilized for this purpose, while the primary (left) vertical axis is utilized for expressing the absolute (unscaled) value of ($\sigma$), expressed in the SI unit of S/m.

Regarding the pre-ionization speed of sound ($a$), and thus the plasma speed ($u$); because ($a$) is independent of the total pressure, it does not need to be visualized as a dependent variable of ($p$). Instead, we only tabulate its values for the five values of ($X_{seed}$) at the beginning of each of the four coming subsections, along with other thermochemical properties that are also independent of the total absolute pressure (such as the pre-ionization mixture’s molecular weight).

In addition to the visualized curves, we also list the numerical values for the variations in ($\sigma$) and ($P_V$). This is useful for the sake of precision, as well as for subsequent utilization of these values by interested readers in external quantitative analysis or comparisons conducted by themselves.

For effective contrast and use of the plots, the ranges of the axes are retained the same regardless of the actual range of the data being visualized. This allows for quick identification of the relative size of the data visualized in different figures, which correspond to different operational cases or the H2MHD channel.

For effective identification of the individual influence of the seed type and the oxidizer type, we intentionally order the four coming subsections such that in the transition from one subsection to the next, only one change is made. This change is either altering the seed type or altering the oxidizer type (but not altering both of them at once; which makes it difficult to separate the influence of each alteration from the other).

3.1. Cesium Seed and Air Oxidizer

The first set of results corresponds to the condition of using cesium vapor as the ionizable gas, and using air as the oxidizer for the hydrogen combustion.

In Table 6, we provide our computed thermochemical characteristics for the seeded gas mixture before ionization, according to the mathematical model presented in the previous section.

It is noticed that as the mole fraction of cesium increases, the mixture’s gas constant ($R_{mix}$) decreases (because the mixture’s molecular weight increases). This is justified by the heavier molecular weight of cesium (Cs) compared to either water vapor (H₂O) or molecular nitrogen (N₂). One mole of cesium (Cs) is 4.744 times heavier than one mole of molecular nitrogen (N₂), and 7.377 times heavier than one mole of water (H₂O). One of the implications of the decline in the mixture’s gas constant ($R_{mix}$) is a decline in the speed of sound ($a$), and thus the plasma’s bulk speed (given that it is always twice the speed of sound in the current study). However, these influences of the additive cesium remain marginal up to a mole fraction of 1%. The remarkable changes at the large mole fraction of 16% are quite virtual, because practically it is not expected to have such a big portion of the seed.

Despite the seed content, the mole fraction ratio ($X_{N2}/X_{H2O}$) is always 1.881, as implied by the stoichiometry reaction in Equation (13), which stipulates that this ratio should be $3.762/2=1.881$ in the cases of air-hydrogen combustion.

Figure 2 shows how the electric conductivity of the hydrogen plasma drops monotonically and nonlinearly as the total pre-ionization absolute pressure increases. This trend is identical for the five seed levels of cesium vapor. This drop in the electric conductivity can be explained by the larger number density of neutral heavy particles (atoms and molecules) in the plasma at higher pressures, which retards the directional movement of the liberated electrons from the ionized cesium atoms due to more frequent collisions with these neutral particles [236,237,238].

Unlike the effect of the pressure, the effect of the cesium additive is non-monotonic. Starting from a small level of (Cs), adding more (Cs) initially improves the electric conductivity. This boost slows down and reaches a peak near a cesium mole fraction of $X_{Cs}=6\%$, and then the electric conductivity starts to drop as more cesium vapor is added in place of water vapor and nitrogen. This peaking behavior of the electric conductivity ($\sigma$) can be attributed to the presence of two factors that impact the electric conductivity due to the free electrons produced by the ionization of a fraction of the cesium atoms. The first factor is the large number density of electrons, due to more cesium atoms; this improves the electric conductivity as more cesium is added. The other factor is the electron-ion interaction (attraction) and electron-electron interaction (repulsion) due to their electric charges. As the number of electrons (and thus the number of cesium ions) increases beyond an appreciable level, the attenuating effect on the electric conductivity caused by these interactions (Coulomb scattering) [239,240,241] overpowers the accentuating effect on the electric conductivity caused by the electrons’ number density; and thus the electric conductivity declines after reaching a peak value. This explains why the electric conductivity with a higher cesium mole fraction of 16% is less than the electric conductivity with a lower cesium mole fraction of 4%. From the shown figure, it may be inferred qualitatively that a pre-ionization cesium mole fraction of 1% or 2% is preferred, because the weak increase in the electric conductivity beyond these levels does not seem to justify such additions of cesium content.

Table 7 lists the numerical values of the plasma electric conductivity as visualized in the previous figure. We also add in the table the electric conductivity values at a pre-ionization cesium mole fraction of $X_{Cs}=6\%$, which approximately corresponds to the maximized electric conductivity for this case (air-hydrogen combustion plasma with cesium seed).

Considering the five selected seed levels (0.0625%, 0.25%, 1%, 4%, and 16%), the highest obtained electric conductivity is $\sigma =55.0610\mathrm{S}/\mathrm{m}$ (at $X_{Cs}=4\%$ and $p=0.0625 \mathrm{atm}$), and this is about 11 times the electric conductivity of seawater. At the higher $X_{Cs}=16\%$ (and the same total pressure of 0.0625 atm), the electric conductivity is slightly lower, with a value of 51.7766 S/m; while at the lower $X_{Cs}=1\%$ (and the same total pressure of 0.0625 atm), the electric conductivity is mildly lower, with a value of 42.2099 S/m. If we take a total pressure of 1 atm (atmospheric MHD channel) as a reference, one can see that the electric conductivity can exceed 10 S/m (twice the conductivity of seawater) with only 1% cesium level.

Using the presented values of the speed of sound (see Table 6) and the values of the electric conductivity (see Table 7) in Equation (16) gives the corresponding ideal volumetric power densities ($P_V$) for the current H2MHD case of cesium seeding and air oxidizer. The ($P_V$) results are visualized in Figure 3, and listed in Table 8. We also add in the table the power density values at a pre-ionization cesium mole fraction of $X_{Cs}=3\%$, which approximately corresponds to the maximized power density for this case (air-hydrogen combustion plasma with cesium seed).

Like the electric conductivity, the power density declines monotonically and nonlinearly as the total pressure increases. Also, like the electric conductivity, the power density exhibits a peak when the cesium level increases, and then it starts to decrease. However, unlike the case of the electric conductivity ($\sigma$), the cesium level for maximizing the power density ($P_V$) has shifted down from being close to 6% to being close to 3%. This shift is explained by the simultaneous decline in the speed of sound as the cesium level increases. This added effect of the speed of sound (which quadratically affects the power density) favors lower cesium values for maximizing the power density.

The tabulated data is helpful in quantifying potential ranges of the power density in the current case of cesium seed and air oxidizer. Considering the five selected seed levels (0.0625%, 0.25%, 1%, 4%, and 16%); the best obtained power density is 1,139.332 MW/m³ (at 4% cesium and 0.0625 atm); and the worst obtained power density is 23.854 MW/m³ (at 0.0625% cesium and 16 atm). At an atmospheric pressure, a promising power density of 362.806 MW/m³ is attainable with 4% cesium, and an attractive power density of 298.752 MW/m³ is attainable with only 1% cesium. Even with 0.25% cesium, a power density of 177.106 MW/m³ can be achieved at atmospheric combustion. These are favorably intense power densities, and all of them exceed well our voluntary threshold criterion of 100 MW/m³; which we found to require a minute amount of cesium seed (at 1 atm) that is not more than 0.08% Cs (at which $P_V=104.932\mathrm{MW}/{\mathrm{m}}^3$).

3.2. Potassium Seed and Air Oxidizer

The second set of results corresponds to the condition of using potassium vapor as the ionizable gas, and using air as the oxidizer for the hydrogen combustion. Thus, compared to the previous set of results in the preceding subsection, the only change made here is changing the seed type from cesium to potassium.

Before presenting the results of this change, it is useful to mention here that due to the higher ionization energy of potassium compared to cesium; potassium is harder to lose its valence electron in its outermost orbital (4s¹) [242,243,244] upon heating compared to cesium, whose outermost valence electron is in (6s¹) [245,246,247]. Therefore, it is expected that this change is accompanied by a drop in the H2MHD performance, primarily because of a decline in the electric conductivity. However, due to the nonlinearity of the problem and the present influence from other factors (particularly the speed of sound), the exact effect of this change in the seed type is not easy to predict analytically. Thus, our computational modeling becomes valuable in understanding quantitatively and qualitatively the implications of this single change.

In Table 9, we provide our computed thermochemical characteristics for the seeded gas mixture before ionization. Similar to the case of cesium seeding, the mixture’s gas constant ($R_{mix}$) decreases (and the mixture’s molecular weight increases) as the level of potassium increases. However, this decrease in the case of potassium additive is much weaker than in the case of cesium additive due to the weaker difference in the molecular weight between potassium (K) and either molecular nitrogen (N₂) or water vapor (H₂O). One mole of potassium (K) is only 1.396 times heavier than one mole of molecular nitrogen (N₂), and 2.170 times heavier than one mole of water (H₂O).

As a reference; the molecular weight of the combustion products without any seeding (34.710170% H₂O and 65.289830% N₂ by mole) is 24.54304 kg/kmol; the unseeded-mixture’s gas constant is 338.7707 J/(kg.K); the unseeded-mixture’s specific heat ratio is 1.243900; and the unseeded-mixture’s speed of sound (at 2,300 K) is 984.486 m/s.

Figure 4 shows how the electric conductivity of the hydrogen plasma monotonically and nonlinearly drops as the total pre-ionization absolute pressure increases. This is qualitatively similar to what was observed when cesium was the additive seed. Also, the existence of a peak electric conductivity at a certain level of potassium seed is similar to what was observed in the case of cesium. This potassium level that maximizes the electric conductivity is again (as was the case in the cesium case) approximately near 6%. It is noticeable that the electric conductivity (or any given level of seeding) dropped strongly when the seed is changed from cesium to potassium. For example, the maximum obtained electric conductivity in the case of cesium was 55.0610 S/m (at $X_{Cs}=4\%$ and 0.0625 atm) is 2.55 times the maximum obtained electric conductivity in the case of potassium, which is 21.6201 S/m (at $X_K=4\%$ and 0.0625 atm).

Table 10 lists the numerical values of the plasma electric conductivity as visualized in the previous figure. We also add in the table the electric conductivity values at a pre-ionization cesium mole fraction of $X_K=6\%$, which approximately corresponds to the maximized electric conductivity for this case (air-hydrogen combustion plasma with potassium seed).

Considering the five selected seed levels (0.0625%, 0.25%, 1%, 4%, and 16%), the highest obtained electric conductivity is $\sigma =21.6201\mathrm{S}/\mathrm{m}$ (at $X_K=4\%$ and $p=0.0625 \mathrm{atm}$). At the higher $X_K=16\%$ (and the same total pressure of 0.0625 atm), the electric conductivity is slightly lower, with a value of 19.7671 S/m; while at the lower $X_{Cs}=1\%$ (and the same total pressure of 0.0625 atm), the electric conductivity is mildly lower, with a value of 16.0343 S/m. If we take a total pressure of 1 atm (atmospheric MHD channel) as a reference, one can see that the electric conductivity is restricted to 6 S/m.

The ideal volumetric power densities ($P_V$) for the current H2MHD case of potassium seeding and air oxidizer are visualized in Figure 5, and listed in Table 11. We also add in the table the power density values at a pre-ionization cesium mole fraction of $X_K=5\%$, which approximately corresponds to the maximized power density for this case (air-hydrogen combustion plasma with potassium seed).

Unlike the case of cesium seeding, the peak power density with potassium seeding occurs near a seed level of 5% (not 3%), which is close to the potassium seeding level that maximizes the electric conductivity. This similarity of the potassium seed level ($X_K$) that maximizes the electric conductivity and the one that maximizes the power density is related to the near-neutral influence of the seed level on the speed of sound for the case of potassium seed (much weaker influence than the case of cesium seed). A nonlinear monotonic decline in the power density ($P_V$) as the total pressure ($p$) increases is noted for any potassium level, and this behavior resembles the one found earlier for the case of cesium seeding.

The power density with potassium seeding is clearly smaller than the one with cesium seeding (when other settings are kept the same). For example; considering the five selected seed levels (0.0625%, 0.25%, 1%, 4%, and 16%), the largest obtained power density at 1 atm in the case of potassium seeding is 141.808 MW/m³ (with $X_K=4\%$). This is 2.56 times the power density at 1 atm and $X_{Cs}=4\%$, which was 362.806 MW/m³. Although potassium seeding is associated with higher speeds of sound compared to cesium seeding, this small advantage for potassium seeding is largely counteracted by the big penalty in the electric conductivity due to changing the seed vapor from cesium to potassium.

Despite the drop in the power density in the case of potassium seeding, the power density can still exceed 100 MW/m³ at 1 atm, with as small potassium levels as 1%. Thus, the H2MHD concept is still viable.

3.3. Potassium Seed and Oxygen Oxidizer

The third set of results corresponds to the condition of using potassium vapor as the ionizable gas, and using pure oxygen as the oxidizer for the hydrogen combustion. Thus, compared to the previous set of results in the preceding subsection, the only change made here is changing the oxidizer type from air (oxygen-nitrogen mixture) to pure oxygen.

Before presenting the results of this change, it is useful to mention here that due to the absence of nitrogen from the combustion products, hydrogen combustion results in pure water vapor (H₂O). Water vapor has a relatively low molecular weight and thus has a relatively high specific gas constant. Therefore, the speed of sound in the current case of oxy-combustion of hydrogen is higher than the speed of sound in the previous case of air-combustion of hydrogen (for the same level of potassium seed). Therefore, it is expected that the change of the oxidizer (from air to oxygen) improves the speed of sound, and thus the plasma speed (given its fixed relation to the speed of sound in our study). It is thus the change in the electric conductivity that determines how eventually the power density is affected by the change of the oxidizer (from air to oxygen). This is difficult to predict without numerical simulations as done in our study. This emphasizes the value and contribution made by our study to the field of MHD power generation, through quantitatively and qualitatively investigating the influence of the oxidizer on the MHD power generation process.

In Table 12, we provide our computed thermochemical characteristics for the seeded gas mixture before ionization. As mentioned earlier, the specific gas constant and the speed of sound in the current case of oxy-hydrogen combustion are higher than their values in the previous case of air-hydrogen combustion. For example, with 1% potassium seed; the speed of sound with oxy-hydrogen combustion (1,114.801 m/s) is 1.135 times the speed of sound with air-hydrogen combustion (982.178 m/s) with the same level of potassium seeding. Due to the larger molecular weight of potassium compared to water waver; the more the added potassium, the higher the mixture’s molecular weight and thus the lower the mixture’s gas constant and the lower the speed of sound in that mixture of water vapor and potassium vapor. The mixture’s specific gas constant (${\gamma }_{mix}$) increases as the potassium level increases, but the decrease in the mixture’s gas constant ($R_{mix}$) is stronger; and thus, ultimately the speed of sound decreases when the potassium level increases.

Figure 6 shows the profiles of the electric conductivity of the hydrogen plasma. The monotonic and nonlinear drop as the total pre-ionization absolute pressure increases is retained, similar to the previous case of air-hydrogen combustion. However, there are two remarkable changes that can be observed from the figure when compared to its counterpart in the previous subsection of air-hydrogen combustion (with the same potassium seeding).

The first change worthy of attention is the overall suppression of the electric conductivity. This can be attributed to the larger effective electron-neutral momentum transfer collision cross-section for water vapor compared to molecular nitrogen [248,249,250,251]. Thus, when the molecular nitrogen in the combustion products is replaced with water vapor, the electrons mobility is retarded, and this suppresses the electric conductivity of the plasma.

The second change worthy of attention is that now with oxy-hydrogen combustion, the electric conductivity with 16% potassium is larger than the electric conductivity with 4% potassium. In the previous case of air-hydrogen combustion, the electric conductivity with 16% potassium was smaller than the electric conductivity with 4% potassium, due to the detrimental effect of over-seeding. In fact, while the electric conductivity with 16% potassium is larger than the electric conductivity with 4% potassium in the current case of oxy-hydrogen combustion, the peak electric conductivity has been already passed at $X_K=16\%$. This means that this peaking phenomenon of the electric conductivity is actually still present in the current case of oxy-hydrogen plasma. The main difference is that this peak was easy to notice in the previous electric conductivity figure corresponding to air-hydrogen plasma with potassium seeding (in the previous subsection) because the potassium mole fraction ($X_K$) that maximized the electric conductivity was closer to the intermediate value of 4% than the terminal value of 16%. Changing the oxidizer from air to oxygen causes this particular ($X_K$) of maximized electric conductivity ($\sigma$) to shift to a larger position, becoming closer to the terminal value of 16% (but below it; located nearly at $X_K=13\%$). Despite the gain in the electric conductivity when the potassium seed mole fraction ($X_K$) is increased from 4% to 16% (or 13%), the gain is too small to justify the need for exceeding the seed level beyond 4%.

Table 13 lists the numerical values of the plasma electric conductivity as visualized in the previous figure. We also add in the table the electric conductivity values at a pre-ionization potassium mole fraction of $X_K=13\%$, which approximately corresponds to the maximized electric conductivity for this case (oxy-hydrogen combustion plasma with potassium seed).

Considering the five selected seed levels (0.0625%, 0.25%, 1%, 4%, and 16%), the highest obtained electric conductivity is $\sigma =16.1854\mathrm{S}/\mathrm{m}$ (at $X_K=16\%$ and $p=0.0625 \mathrm{atm}$). At the lower $X_K=4\%$ (and the same total pressure of 0.0625 atm), the electric conductivity is mildly lower, with a value of 13.9154 S/m. If we take a total pressure of 1 atm as a reference, one can see that the electric conductivity is below 5 S/m (the approximate electric conductivity of seawater) regardless of the seeding level.

The ideal volumetric power densities ($P_V$) for the current H2MHD case of potassium seeding and oxygen oxidizer are visualized in Figure 7, and listed in Table 14. We also add in the table the power density values at a pre-ionization cesium mole fraction of $X_K=9\%$, which approximately corresponds to the maximized power density for this case (oxy-hydrogen combustion plasma with potassium seed).

The electric conductivity ($\sigma$) at a potassium seed level ($X_K=16\%$) was found to be larger than the electric conductivity at ($X_K=4\%$). Similarly, the volumetric power density ($P_V$) at a potassium seed level ($X_K=16\%$) is larger than the power density at ($X_K=4\%$). However; because the speed of sound (and thus the plasma speed) at ($X_K=16\%$) is larger than its value at ($X_K=4\%$), the gap in ($P_V$) is diminished when the seed levels of 4% and 16% are compared.

When comparing the range of ($P_V$) here with those in the previous subsection (air-combustion instead of oxy-combustion); we notice that although the power density dropped, this drop is not large. For example, at 1 atm and 4% potassium, the power density here with oxy-combustion is 111.573 MW/m³, while it was 141.808 MW/m³ in the previous case of air-combustion. The difference (30.2350 MW/m³) is 23.9% of the average of both values (126.6905 MW/m³).

For exceeding the threshold power density of 100 MW/m³, the potassium seed level should not be below 3% (at which we found that $a=1,103.376 \mathrm{m}/\mathrm{s}$ and $\sigma =3.4118\mathrm{S}/\mathrm{m}$; thus, $P_V=103.841\mathrm{MW}/{\mathrm{m}}^3$)

3.4. Cesium Seed and Oxygen Oxidizer

The fourth (and final) set of results corresponds to the condition of using cesium vapor as the ionizable gas, and using pure oxygen as the oxidizer for the hydrogen combustion. Thus, compared to the previous set of results in the preceding subsection, the only change made here is changing the seed type from potassium to cesium.

In Table 15, we provide our computed thermochemical characteristics for the seeded gas mixture before ionization. The speed of sound drops faster with the cesium fraction compared to the previous case of potassium seeding, due to the large influence of cesium on the mixture’s gas constant compared to potassium; this is because the cesium atom is 3.399 times heavier than the potassium atom. For example; at the same mole fraction of 4%, the speed of sound in the case of cesium seeding is 1,002.552 m/s, which was larger (1,097.826 m/s) in the case of potassium seeding. For comparison purposes, the speed of sound in pure water vapor (at 2,300 K) is 1,120.682 m/s. Thus, introducing a fraction of 4% cesium vapor causes a relative drop in this original speed of sound (1,120.682 m/s) by 118.130 m/s or 10.54%. On the other hand, introducing a fraction of 4% potassium vapor causes a relative drop in the original speed of sound (1,120.682 m/s) by only 22.856 m/s or 2.04%.

In addition, it may be useful to add that for unseeded water vapor (H₂O), the specific gas constant is 461.5223 J/(kg.K), and the specific heat ratio is 1.183162.

Figure 8 shows how the variation of the electric conductivity of the hydrogen plasma with the total absolute pressure at the five selected cesium levels. The variations show a monotonical nonlinear drop as observed in all previous cases (the three preceding subsections). However, compared to the previous subsection particularly, the electric conductivity is largely boosted here after replacing the seed from potassium (the previous subsection) to cesium (the current subsection). For example, at 1 atm and $X_{Cs}=4\%$, the electric conductivity ($\sigma$) here is 11.0877 S/m. But in the previous case, at 1 atm and $X_K=4\%$, this electric conductivity was 3.7030 S/m. The ratio of the two values of ($\sigma$) is 2.99.

Also, as in the previous set of results (oxy-hydrogen combustion with potassium seeding), the seeding level at which the electric conductivity is maximized is larger than that level in the case of air-hydrogen combustion. That seeding level of maximum ($\sigma$) is approximately 13% in the current case of oxy-hydrogen combustion with cesium seeding (while it was near 6% in the case of air-hydrogen combustion with cesium seeding).

Table 16 lists the numerical values of the plasma electric conductivity as visualized in the previous figure. We also add in the table the electric conductivity values at a pre-ionization cesium mole fraction of $X_{Cs}=13\%$, which approximately corresponds to the maximized electric conductivity for this case (oxy-hydrogen combustion plasma with cesium seed).

Considering the five selected seed levels (0.0625%, 0.25%, 1%, 4%, and 16%), the highest obtained electric conductivity is $\sigma =43.5831\mathrm{S}/\mathrm{m}$ (at $X_{Cs}=16\%$ and $p=0.0625 \mathrm{atm}$). At the lower $X_{Cs}=4\%$ (and the same total pressure of 0.0625 atm), the electric conductivity is mildly lower, with a value of 37.7620 S/m (thus, 13.36% of the $\sigma$ value of 43.5831 S/m at 16% cesium is lost). If we take a total pressure of 1 atm as a reference, an electric conductivity of 10 S/m is achievable with a cesium mole fraction of 3% (at which the electric conductivity is 10.2333 S/m).

The ideal volumetric power densities ($P_V$) for the current H2MHD case of cesium seeding and oxygen oxidizer are visualized in Figure 9, and listed in Table 17. We also add in the table the power density values at a pre-ionization cesium mole fraction of $X_{Cs}=5\%$, which approximately corresponds to the maximized power density for this case (oxy-hydrogen combustion plasma with cesium seed).

Like the case of cesium seeding with air-hydrogen combustion, but unlike the previous case of potassium seeding and oxy-hydrogen combustion; the power density here with 4% Cs is larger than the power density with 16% seed. This is a combined effect of how the electric conductivity ($\sigma$) and the speed of sound ($a$) respond to changes in the cesium mole fraction ($X_{Cs}$); which is different from their responses to changes in the potassium mole fraction ($X_K$).

The power density with cesium seeding (the current subsection) is clearly larger than the one with potassium seeding (the previous subsection).

A power density of 100 MW/m³ or more is easily attainable at 1 atm, requiring a small mole fraction of cesium seed such as 0.25% only. Even with $X_{Cs}=0.20\%$, we obtain $P_V=101.838\mathrm{MW}/{\mathrm{m}}^3$.

4. Conclusions

In the current study, we reported novel results about our proposed concept of hydrogen magnetohydrodynamic (H2MHD) power generation, particularly results about potential volumetric electric power densities under different operational conditions. We forecast the performance of a Mach-two supersonic channel, operating with hydrogen-combustion plasma at a uniform temperature of 2,300 K. The influence of the pressure, seed level, seed type, and oxidizer type was presented graphically and numerically.

The following findings can be stated:

• Hydrogen magnetohydrodynamic (H2MHD) direct power extraction is an innovative technology for large-scale commercial electricity generation. It has a major environmental advantage of zero greenhouse gas (GHG) emissions, unlike any fossil hydrocarbons.
• H2MHD generators permit concentrated power generation, with theoretical power densities as high as 300 MW/m³ under a controlled condition of 1 atm and 1% cesium seeding (mole fraction), while using air as a conventional oxidizer.
• Using cesium vapor as a seeding alkali metal for producing free electrons (as the main charge carriers in the plasma) is highly advantageous compared to potassium. The gain can be more than doubling the output power.
• Using air as an oxidizer is preferred to using oxygen in terms of the plasma electric conductivity and the output power. However, it should be noted that this preference for air as an oxidizer does not mean preferring an air-combustion process over an oxy-combustion process. The reason is that we here assumed that the temperature remains the same when the oxidizer is changed from air to oxygen. However, practically oxy-combustion leads to elevated temperatures [252,253,254], which strongly boosts the electric plasma conductivity, and thus the electric power output. Therefore, favoring air as an oxidizer is an outcome we found based on its influence on the chemical composition of the combustion products only, while factoring out its influence on the temperature of these combustion products. When the influence of the temperature due to oxy-fuel combustion is factored in, then the use of pure oxygen becomes preferred, because the large gain in the electric conductivity due to the increase in temperature is more important than the small decline in the electric conductivity due to the altered chemical composition of the plasma.
• Increasing the total static pressure monotonically decreases the electric conductivity and thus decreases the power density and performance of the H2MHD generator.
• Increasing the alkali metal vapor seed amount monotonically increases the mixture’s molecular weight, and thus monotonically decreases the mixture’s specific gas density and in turn monotonically decreases the speed of sound in the mixture.
• Increasing the alkali metal vapor seed amount increases the electric conductivity up to a certain level, where a peak electric conductivity is reached and then it declines with further seeding.
• With air-hydrogen combustion and cesium seeding; the electric conductivity is maximized at a pre-ionization seeding mole fraction near 6%, while the power density is maximized (with a fixed Mach number) at a lower mole fraction near 3%.
• With air-hydrogen combustion and potassium seeding; the electric conductivity is maximized at a pre-ionization seeding mole fraction near 6%, while the power density is maximized (with a fixed Mach number) at a lower mole fraction near 5%.
• With oxy-hydrogen combustion and potassium seeding; the electric conductivity is maximized at a pre-ionization seeding mole fraction near 13%, while the power density is maximized (with a fixed Mach number) at a lower mole fraction near 9%.
• With oxy-hydrogen combustion and cesium seeding; the electric conductivity is maximized at a pre-ionization seeding mole fraction near 13%, while the power density is maximized (with a fixed Mach number) at a lower mole fraction near 5%.