Abstract
The main goal of this paper is to analyze the exergy losses during the shock and rarefaction wave of hydrogen-air mixture. First, detonation parameters (pressure, temperature, density, and species mass fraction) are calculated for three cases where the hydrogen mass fraction in air is 1.5%, 2.5%, and 5%. Then, exergy efficiency is used as objective criteria of performance evaluation. A two-dimensional computational fluid dynamic code is developed using Finite volume discretization method coupled with implicit scheme for the time discretization (Euler system equations). A seven-species and five-step global reactions mechanism is used. Implicit total variation diminishing (TVD) algorithm, based on Riemann solver, is solved. The typical diagrams of exergy balances of hydrogen detonation in air are calculated for each case. The energy balance shows a successive conversion of kinetic energy, and total enthalpy, however, does not indicate consequent losses. On the other hand, exergy losses increase with the augment of hydrogen concentration in air. It obtained an exergetic efficiency of 77.2%, 73.4% and 69.7% for the hydrogen concentrations of 1.5%, 2.5%, and 5%, respectively.
1. Introduction
The Ramjet propulsion principle was invented at the beginning of the 20th century by Rene Lorin [1, 2], who published it in the technical review “Aérophile” in 1913. A Ramjet is an air-breathing jet engine which uses the engine’s forward motion to compress incoming air, without a rotary compressor. The development of this propulsion method was slow and was practically started only from the 1950s. Some years later appeared the pulsed detonation engine (PDE), which is a type of propulsion system that uses detonation waves to combust the fuel and oxidizer mixture [3, 4].
The excitement about pulse detonation engines (PDEs) stems from the fact that, theoretically at least, they can operate at thermal efficiencies about 30% higher than conventional combustion processes. This is due to the high detonation velocities between 1500 m·s−1 and 2200 m·s−1 in a hydrocarbon fuel/air mixture; for example, its thermodynamic cycle is effectively a constant volume or Humphrey cycle [5–7]. The Humphrey cycle can be considered a modification of the Brayton cycle in which the constant-pressure heat addition process of the Brayton cycle is replaced by a constant-volume heat addition process [8].
Due to this thermodynamic efficiency and mechanical simplicity, the PDE technology is more efficient when compared with current engine types. Additionally, the PDE can also provide static thrust for a ramjet or scramjet engine or operate in combination with turbofan systems showing an enormous potential to be applied in many sectors of the aerospace, aeronautic, and military industries. However, there are still engineering challenges that must be overcome, such as to improve methods for initiating the detonation process and develop materials able to withstand with extreme conditions of heat and pressure.
In the PDE technology the combustion takes place in an open-ended tube in which fuel is mixed with air and detonated. As the detonation wave travels down the tube at supersonic speed, a refraction wave propagates into the combustor, and exhaust products leave the chamber. Pressure within the chamber jumps down to charging conditions, which draws in fresh fuel and air, and the cycle is repeated. Each pulse lasts only milliseconds [9, 10]. Detonation is initiated by a predetonator, containing a more readily detonable mixture of fuel air, or by DDT (deflagration detonation transition), a process by which a flame, initiated at the closed end of a duct by a weak spark, accelerates to speeds on the order of 1000 m·s−1 at which point a detonation wave is initiated within the flame-shock complex [11–13].
The existence of flames for problems with complex chemistry in the combustion approximation has been considered by many authors; see, for example, [14–18]. It has been pointed out in [19], that reaction mechanisms including chain branching can lead to significant qualitative differences in the initiation of detonations in comparison to single step reactions. The use of the exergy concept can be used with advantage to best understand the detonation problem [20]. The exergy gives a quantitative and illustrative description of the convertibility of different energy forms and is a function of the system and environment, whereas entropy is a function only of the system. Total exergy, thermomechanical and chemical, has been defined as the work which can be obtained by taking the system, by means of reversible processes, from its initial state to a state of total equilibrium (dead state), thermal, mechanical, and chemical (“unrestricted equilibrium”), with the environment. This definition entails that the system is in a state of internal equilibrium (uniformity of all intensive properties) and that the environment is in a state which is in internal equilibrium and does not change during time. The exergy of a fluid stream depends upon the choice of the reference environment. This choice, which seems obvious according to some authors [20–22], is, on the contrary, problematic, since every usual environment is a chemical state of nonequilibrium.
A simplified model of adiabatic detonation process in gases was considered by Petela [20] using the concept of exergy. The exergy loss for high pressure shock is larger than that for reactions of combustion. On the other hand, the total exergy loss for deflagration is also larger than that for detonation. Huntchins and Metghalchi [21] analyzing the thermal efficiency show that the Humphrey cycle reveals good advantage compared with Brayton cycle using the same fuel. They also compare the effectiveness of both cycles and conclude that the results using the Humphrey cycle have significant advantage comparing with a Brayton cycle based on the same fuel (methane). Also, the research work of Wintenberger and Shepherd [22] finds that this efficiency cannot be precisely translated into propulsive efficiency. Indeed, the results are only useful in comparing detonations with other combustion modes. They also find that the efficiency of cycles based on detonation and constant-volume combustion is very similar and superior to a constant-pressure combustion (Brayton) cycle when compared on the basis of pressure at the start of the combustion process.
This paper reports an exergy analysis on the detonation wave between the pulse shock wave to the end of the rarefaction wave for various concentrations of the hydrogen in air. Numerical simulation of a two-dimensional detonation of the hydrogen-air mixture is considered. The coupled hydrodynamical model and chemical model integrating the Euler system equations are combined with a detailed chemical reaction model.
2. Numerical Method
Supersonic combustion, called detonation, occurs coupled to a shock wave traveling at supersonic speed. If premixed gases inside a tube closed at both ends are ignited at one end, a laminar flame first develops, traveling pushed by the expanding hot products behind. The inverted small pressure jump across the flame generates local pressure pulses that wrinkle the flame, create turbulence, extend the burning area, and increase the burning rate, with a positive feedback that, if positively combined with pressure pulses reflected from the other end, might compress the fresh mixture to the autoignition temperature.
2.1. Chemical Model
Consider the reaction of hydrogen with oxygen to yield water that is usually written as
The detonation of hydrogen in air must start by dissociation of a hydrogen molecule (it requires less energy than oxygen dissociation): where M stands for an unspecified molecule (hydrogen, oxygen, a hot wire). The following basic reactions take place:
The chemical kinetic scheme involves seven species, H2, O2, , , , , and N2 and the above-described 5 elementary reversible reactions. The mechanism is the one proposed by Balakrishnan and Williams [23], and the empirical parameters from each one of the reactions [24] are depicted in Table 1. is the “rate coefficient,” which generally depends on temperature , and activation energy and can be defined as follows: where A is the preexponential factor as in Arrhenius law and is the activation energy. The three parameters , , and are determined in practice experimentally.
The chemical production rates are computed using the methodology suggested by Warnatz [25]. In general, for the reaction products not just to species , it can be defined:
2.2. Euler System Equations
The contributions of pressure-driven diffusion, bulk viscosity, and radiative heat transport can be neglected; then the balance equations of mass, momentum, and energy for the two-dimensional unsteady flow for a multicomponent chemically reacting gas mixture are written in the integral form as The state equation can be defined as follows: In these system equations the density is the total density of the whole gas in the position , and at time . is calculated as sum of mass fraction and density of each species: The above equation system is rewritten in a control volume as follows: where
In these equations and are the velocity vector components in the directions and , the mixture temperature, ρ the mixture density, the pressure, and are the internal and total energies, respectively. For each of the chemical species, is the mass fraction, is the molecular, and is the specific enthalpy. Then, the following equations of state can be also defined:
2.3. The Solution Procedure
The numerical methods that have been developed in this code for hydrogen and combustion are detailed in [26, 27]. These include lumped parameters, CFD-type codes, and also the CAST3M code [28] developed at the French Atomic Centre Energy (CEA). CAST3M is a general finite-element/finite-volume code for structural and fluid mechanics and heat transfer.
The CAST3M code is used to model both hydrogen dispersion and combustion phenomena over a very wide range of flow regimes, from nearly incompressible flow to compressible flow with shock waves. This could be done by developing a most suitable (i.e., accurate and efficient) numerical method for each type of flow or developing a single method suitable for all flow regimes. We are not aware of any single numerical method able to treat accurately and efficiently flow regimes which range from slowly evolving nearly incompressible buoyant flow to fast transient supersonic flow. However, in practice, for detonation simulation, the use of shock-capturing methods is preferred. In the hydrogen risk analysis tools CEA is developing; the choice was made to develop in the same computational platform CAST3M code for distribution calculations, an efficient pressure-based solver using a semi-implicit incremental projection algorithm which allows the use of “large” time steps, for combustion calculations, a robust and accurate density-based solver using a shock-capturing conservative method.
A brief description of the method we used here, typical for such hyperbolic system equations, follows. Consider the system of hyperbolic conservation laws, (2.10), with the initial conditions:
Given a uniform grid with time step and spatial mesh size , we define an approximation of at the point , , by the finite volume formula:
Figure 1 shows the structured mesh used in this problem. Finding the exact solution to the Riemann problem for the equation system (2.9) is complex. Even if we would find it, it would be too complex to use it at each cell boundary, as it needed in (2.13). Therefore, we have to look for some approximate Riemann solver. Here, we adopt the approach described in [29], which we have already considered in case of conservation laws. Let us rewrite the homogeneous hyperbolic system (2.9) in primitive variables for one direction: Consider the Riemann problem for (2.14) at each cell boundary , that is, with the following initial data: Following [10, 11], we calculate the Jacobian matrix in the average state: The intermediate state in the solution of the Riemann problem is where the eigenvalues and the eigenvectors of the Jacobian matrix and are the coefficients of eigenvector decomposition of : With the choice of the primitive variables, they are given by the following expressions: where are the components of and is the k th component of . Recalculating into the conservative vector , we fully determine the Godunov-type scheme (2.17) for the homogeneous system (2.14).
3. Exergy Calculation
A practical definition of exergy (E) for a closed system is the maximal work that can be obtained when a thermodynamic system is allowed to attain equilibrium with an environment defined as the “dead state,” typically defined as the system ambient. If that ambient is at temperature and pressure, , the exergy is where , and denote, respectively, the internal energy, volume, and entropy of the system, and , and are the values of the same properties in the system which were called dead state.
The exergy expression can be generalized considering as a heat transfer across a system boundary where the temperature is constant at and taking in account the mass flow across the boundary of a control. The exergy rate balance for a control can be derived using these described approaches, where the control volume forms of mass and energy rate balances are obtained by transforming the closed system forms. The exergy accompanying mass flow and heat transfer can be written as where and are the exit and inlet specific exergy and is the heat transport rate. The sum of the thermomechanical and chemical exergies is the total exergy associated with a given system at a specified state, relative to a specified exergy reference environment. The chemical exergy is defined as or in the following form, when the entropy at () is calculated: where , , are the mass fraction of reactants and products of involved reactions and , , , and are the specific enthalpy and specific entropy for the reactants and products, respectively. The specific enthalpy terms are determined using the enthalpy of formation for respective substances. The specific entropy appearing in the above equation is absolute entropy. Even if the logarithmic term normally contributes a few percent, it will not be neglected in the present paper. The Gibbs function for each one of the components could be expressed as follows:
The energy associated with a specific state of a system is the sum of two contributions: the thermomechanical contribution and the chemical contribution. On unit mass, the total exergy is: where (2.1) is the thermo mechanical contribution and (2.2) is the chemical contribution.
The molar exergy of each reactant and product, at temperature , is given by The standard chemical exergies of the reactants and products are given in Table 2.
The specific chemical exergy of a material stream is that part of the exergy that results from the sum of the chemical potential differences between the pure reference substances and the pure material stream components at and . The specific mixing exergy of a material stream is that part of the exergy resulting from the sum of two contributions. (i)The differences in potential between the exergies of the environmental components in environmental concentrations and pure environmental components at and . (ii)The difference in potential of the exergy due to mixing of the pure components into a mixed material stream at and .
The different exergy components are listed in Table 3.
Finally, it can be defined as a “performance factor” the ratio of outlet exergy to inlet exergy:
4. Results and Discussion
The structure of the pulse detonation and the shock tube is depicted in Figures 2 and 3, respectively.
The detonation field parameters in this work were computed by using the Euler multispecies equations of detonation of H2-air mixture. They were solved using the first author’s solution procedure [26, 27]. Analyses were carried out for different hydrogen concentrations in air optimum pressure with corresponding temperatures; species concentrations were determined and given first case in figures as function of time. As an example, for first case where the mass fraction of H2 is equal to 1.5%, peak detonation and end rarefaction density of the selected location were taken at 1.4 kg·m−3 and 1.1 kg·m−3, respectively. As a result, optimum shock and rarefaction temperatures were found to be 1930 K and 700 K, respectively. At the same time, the optimum pressures were found to be 18 bar and 6 bar, respectively for the same concentration of hydrogen in air.
As seen in Figure 4 representing the pressure distribution along the tube at 0.003 ms, the optimum pressure at the peak of detonation was found to be 18 bar, 22 bar, and 36 bar for 1.5%, 2.5%, and 5% of hydrogen mass fractions, respectively. Figure 5 shows the result of the temperature at 0.003 ms corresponding to the end of reactions for the case of hydrogen concentration equal to 1.5%. Optimum values vary from 1930 K for the peak of detonation to 700 K for the end of rarefaction. These results show that the structure of the detonation wave is dependent on hydrogen concentration.
Figure 6 shows the shock pressure variation with time for the three different concentrations of hydrogen in air.
The first case considers a hydrogen-air mixture at an initial pressure = 1 bar and initial temperature = 298 K. On the left side of the tube (Figure 3) a pressure pulse (10 bars) is considered to initialize the shock wave. At s, the shock wave moves from the left to the right side of the tube and initiates the reactions. Shock waves generated due to the formation of a new reaction front are called “reaction shocks” [20]. The reaction shock overtakes the shock front at around 27 μs, strengthening it and sharply increasing its propagation speed; that is, a detonation wave is formed. Subsequently, the detonation wave speed decreases gradually and a high-frequency, low-amplitude propagation mode is established. In this study no reaction delay is considered, so when the shock wave is initialized, the reactions produced energy. Generally, chemical reactions behind the shock front start generating weak shock waves after approximately 20 μs [31]. The plots show the variation in shock front speed with time. Initially the shock wave travels at a constant speed (determined by the initial conditions) slightly greater than 1500 m·s−1, 1800 m·s−1, and 2200 m·s−1 for respective hydrogen concentrations of 1.5%, 2.5%, and 5%.
Figures 6 and 7 show the variation with time of pressure and density as function of time and for mass fraction of hydrogen in air: 1.5%, 2.5%, and 5%. As can be seen from these figures, the peak of detonation is characterized by the high-pressure values followed by the rarefaction wave. During the rarefaction wave, the pressure, temperature, and density decrease; this decrease is stopped by the reflected wave from the end of the tube. The most detonation energy is absorbed by the shock when the wave is reflected. This is because by assuming that the detonation is instantaneous and all hydrogen is consumed in the first way of the shock wave, then the energy is completely kinetic, and, as a result, the density and pressure decrease.
Figure 8 shows the mass fraction of H2O, calculated as function of time for the three studied cases: 1.5%, 2.5%, and 5% of hydrogen concentration. This figure gives us a good idea about mass fraction of H2O during the detonation for each of hydrogen concentration in air. In fact, the contribution of H2O in exergy calculation is not negligible since the standard chemical exergy of H2O, O2, and N2 are, respectively, 8.63 kJ·mol−1, 3.90 kJ·mol−1, and 0.60 kJ·mol−1.
Mass fraction for the H2, OH, and H2O system is plotted in Figure 9. Here, the energy levels were chosen to be consistent with those realized in the detonation calculation. As seen in this figure the hydrogen, radical OH, and product H2O as function of time for the referent element shown on Figure 3, optimum values at the end of reactions, were found to be 0.04%, 0.7%, and 1.2% carried out for species H2, OH, and H2O. These results are presented for the case when the initial mass fraction of hydrogen in air is assumed to be equal to 2.5%.
It is possible to consider the shock and rarefaction processes as the separate successive steps. Therefore, the exergy loss due to the irreversible detonation and the end of rarefaction can be calculated as the product of the initial temperature and the difference between entropies in the start of the shock wave and the end of rarefaction wave. The typical diagrams of exergy balances of hydrogen detonation in air are shown in Figure 10 as function of hydrogen concentration (1.5%, 2.5%, and 5%). Data were calculated for the point of steady detonation. It was shown that exergy losses increase with the augment of hydrogen concentration in air; consequently low hydrogen concentrations in air achieve slight higher efficiencies. It was also verified that the exergy efficiency decreases for all the hydrogen concentrations as a function of the time. It was obtained an exergetic efficiency of 77.2%, 73.4%, and 69.7% for the hydrogen concentrations of 1.5%, 2.5% and 5%, respectively.
On the other hand, the energy balance shows a successive conversion of kinetic energy and total enthalpy, however, does not indicate consequent losses.
5. Conclusions
The CFD code was developed using finite volume discretization method coupled with implicit scheme for the time discretization. The code solves the Euler system equations including a chemistry model, with 5 reactions and 7 species, using an implicit total variation diminishing (TVD) algorithm based on Riemann solver.
There exists a large literature about the numerical discretization of the detonation of hydrogen-air mixture problems: Eulerian schemes, Lagrangian schemes, front tracking techniques, and level set methods, among others. Our main goal, in the present work, is to use detonation parameters, from the beginning to the end of detonation of hydrogen-air mixture to calculate and analyze the exergy loss and the efficiency. The main steps of H2-air mixture detonation calculation were explained. The developed program is able to calculate the pressure, temperature, density, and mass fraction, for each species from the chemical model. The exergy efficiency was used as objective criteria of performance evaluation.
We conclude in this work that the exergy losses increase with the augment of hydrogen concentration in air, and the efficiency of hydrogen detonation, with low concentration in air, is slightly elevated than the detonation using high hydrogen concentration.
Acknowledgments
The authors would like to express our gratitude to the Portuguese Foundation for Science and Technology (FCT) for the given support to the grant SFRH/BPD/71686 and to the project PTDC/AAC-AMB/103119/2008 and also to The French Atomic Centre Energy (CEA) for all the provided help.