1. Introduction
In the very early development of jet-propulsion engines, it was known from the thermodynamic analysis cycle that an engine based on a constant-volume combustion process achieves higher thermodynamic efficiency than a constant pressure engine. The earliest non-piston-engine-type prime mover employing constant volume combustion with a deflagrative and not a detonative reaction was the Holzwarth gas turbine manufactured by Brown-Boveri (now ABB) in Switzerland during the early part of the last century, but its success was limited [
1]. Eidelman, Grossmann, and Lottati [
2] and Ma, Choi, and Yang [
3] have summarized that the first reported work on intermittent detonation is attributed to Hoffman in 1940, using acetylene and benzene as fuels with oxygen. After the work was terminated during World War II, Nicholls and co-workers reinitiated the effort in the 1950s by experimenting with a series of single- and multiple-cycle detonation experiments with different mixtures of hydrogen, oxygen, acetylene, and air in a six-foot tube. The Naval Postgraduate School (NPS) reexamined the pulse detonation engines (PDE) concept in the late 1980s and successfully demonstrated the self-aspirating feature of air breathing PDE using ethylene/oxygen and ethylene/air mixtures. Since then, there has been a growing interest in PDEs as a propulsion technology for both air breathing and rocket systems.
In the context of repetitive mode detonative burning to develop thrust [
4,
5,
6], PDEs represent one of the pressure-rise unsteady propulsion systems, which differ from conventional propulsion frameworks [
3,
7]. This potential alternative combustor technology introduces a rapid detonation wave as a more thermodynamically efficient means for converting chemical to mechanical energy and thus generating far higher kinetic energy [
8,
9,
10,
11] compared to the normal deflagration combustion process [
12]. PDE offers numerous potential advantages; however, it also suffers from several drawbacks too. A list of advantages and disadvantages are summarized in
Table 1. The theory, operational considerations, and research done for PDE will be described in the next section.
Awareness of environmental and energy crises has prompted tremendous efforts such as the Clean Sky JTI Projects by the European countries, The Environmentally Responsible Aviation Project (ERA) by National Aeronautics and Space Administration (NASA), and several more. The aviation industries have shifted their strategy to use alternative fuels based on biofuels. The use of drop-in fuels and blended fuels in aircraft engines has significantly attracted the attention and interest of engineers and researchers throughout the world. Drop-in fuels need minor or no modifications at all in the aircraft engine in service. It offers a future ‘greener’ aircraft and less dependency on crude oil. Following the successful flights of many commercial aircraft running with different biofuels, these have become a viable choice to sustain the environment as well as conserve energy. However, there are shortcomings associated with the use of biofuels alone in aircraft engines, such as in terms of thermodynamic efficiency and performance. Running an engine using biofuels with pressure-rise combustors would certainly be a good choice and strategy to satisfy greener technology with better performance. Indeed, it is believed that such alternative combustor technologies fueled by alternative fuels could meet the 2050 emissions targets plan for aviation.
Four biofuels, namely Jatropha Bio-synthetic Paraffinic Kerosene (JSPK), Camelina Bio-synthetic Paraffinic Kerosene (CSPK), Microalgae Biofuel, and Algal Biofuel have been evaluated as pure fuels and are compared here with conventionally used kerosene and acetylene fuels. These particular biofuels were chosen because of previously reported successful use in conventional engine test flight programs and because their fuel properties are available in the published literature, as listed in the
Appendix A. The 3rd and 4th Generation’ biofuels derived from “algae-to-biofuels” and microalgae biofuel technologies are based on algae biomass processing for biofuel production and metabolic engineering from oxygenic photosynthetic microorganisms [
27]. These are considered to provide a technically viable alternative energy resource, overcoming the major drawbacks associated with first and second generation biofuels [
28]. The other main feature of such algae and microalgae-based biofuels is they have the highest oil yield compared to other types of feedstocks because of their unique fast growing capabilities. Microalgae, in particular, have been suggested as potential candidates for fuel production, capable of meeting the global and sustainable demand for transport fuels [
29] because of a number of key advantages, including higher energy yields per hectare, higher photosynthetic efficiency, higher biomass production, higher growth rates, and a non-requirement of agricultural land compared to other energy crops [
30,
31]. (Some microalgae have also been reported as good producers of hydrogen which offers higher energy potential and almost no pollution [
32]). Microalgae biofuel has properties similar to those of petro-diesel in terms of density, viscosity, flash point, cold flow, and heating value. None of the other potential sources of biodiesel present as realistic an option for replacing petro diesel sustainablyas microalgae do [
33].
However, to date, almost no efforts have been made to study the use of such alternative fuels under detonation combustion conditions. Although studies have been made of heavy-hydrocarbon fuel such as Jet Propellant (JP10), none have been made of other commercialized alternative fuels. Since alternative fuels are targeted to be used in the near future, it is certainly worth exploring the wider capability of these fuels as well. The uniqueness of this paper is to assess the behavior of such alternative fuels in terms of physical and chemical properties for changes in the different initial conditions. This study only uses one-step chemistry reactions for a start in order to make a straight comparison between different fuels and to assess whether these might be sufficiently accurate to be useful. Thus, remaining differences within the experiment are most probably due to not using full multi-step chemistry and leave open this extension for further investigations and improvements.
2. Theory, Process and Previous Works on PDE
PDE detonation is modeled as a normal shock wave or Zel’dovich–von Neumann–Doering (ZND) detonation wave, advancing into the undisturbed fuel-air mixture of a uniform cross-sectional area tube, which is almost at rest for combustor entry conditions [
34]. This is then followed by Rayleigh type combustion [
35]. The whole process satisfies the Chapman-Jouguet (CJ) condition, which requires that the local Mach number at the termination of the heat expansion region be choked [
34]. CJ theory requires chemical reactions to be represented by heat discharge in an infinitesimally thin shock front that brings the material from a starting state on the inert Hugoniot line to a subsequent CJ point state [
20]. The CJ point also forms a tangent from the initial to final state on a Pressure-Volume diagram (p–v diagram) equivalent to the Rayleigh heating process. It is difficult to evaluate the relative performance of air-breathing PDEs with respect to conventional steady-flow propulsion systems without performing a full unsteady computational analysis because of the intrinsically unsteady nature of the flow field due to the detonation process [
16,
36,
37].
In a conventional Brayton cycle, the heat injection process has the maximum exergy, which is fixed by the compressor’s delivered pressure and the maximum temperature allowed by the cycle. Therefore, the exergy can be increased if the heat injection process follows a different thermodynamic cycle path [
38]. The thermodynamic cycle of ideal PDE is similar to the ideal Brayton cycle, while the Humphrey cycle is considered a modification to the Brayton cycle in which the constant-pressure heat addition process is replaced by a constant-volume heat addition process [
34]. The Humphrey cycle is much more efficient than the Brayton cycle, in which a very rapid burning takes place. Due to the rapidity of the process, there is not enough time for pressure equilibration, and the overall process is thermodynamically closer to a constant volume process than to the constant pressure process typical of conventional propulsion systems [
39]. The thermodynamic efficiency of Chapmen–Jouget detonation has minimum entropy generation along the Hugoniot curve as compared to other combustion modes, which appear to have a potential thermodynamic advantage [
13,
40]. In general, the Humphrey cycle consists of four processes. The first is an isentropic compression. This compression occurs ahead of the detonation wave in PDEs. Compression is followed by constant volume combustion. Another isentropic process expands the combustion products back to atmospheric pressure. In the PDE, the rarefaction waves cause this expansion process. Finally, an isobaric process brings the cycle back the start of the cycle.
Figure 1 illustrates the detonation process in order.
Researchers throughout the world have clearly already done quite extensive work on PDE. In addition Ma, Choi, and Yang [
41] have summarized the findings of both numerical and experimental work on air breathing PDE using hydrogen fuel in a review article.
Figure 2 attempts to summarise by key words all the topics covered across their identified themes and illustrates the wide breadth of PDE research. Roy et al. [
42] and Kailasanath [
39,
43,
44] have also presented detailed review discussions of PDE work, so we will not repeat these here. The present work can be categorized as a first examination of the applicability and feasibility of the selected alternative fuels for detonation combustion. It should be noted, however, that two previous numerical studies have been conducted investigating the detonation characteristics of biofuel and the feasibility of biogas; by Shimada et al. [
45], using bio-ethanol, and Dairobi et al. [
46], using biogas. Shimada at al. utilized STANJAN for 2D bio-ethanol chemical reactions to study the two-phase detonation of bio-ethanol/air, which showed that the biofuel resulted in a smaller cell size. Meanwhile, the biogas studies suggested that this requires supplementary additives for higher detonation pressure. Our own numerical work utilises ZND Theory and CJ Theory in a zero dimensional analysis under a few basic assumptions, which will be discussed in the theoretical model framework. As a first attempt, the approach was to employ single tube, single phase, and single cycle processes. The theoretical formulation and numerical framework will be discussed in the following section.
3. Methodology: Theoretical Formulation and Numerical Framework
The model adopted here uses an open-ended constant-area tube geometry in a single cycle operation. It incorporates appropriate expressions, including the Rankine-Hugoniot Equation, Rayleigh Line Equation, species mole and mass fraction of the reactants, enthalpies-of-formation, and ideal-gas normal shock equations. The computational results from our analyses have been verified using the available limited published data from the literature to ensure the consistency of our model across an acceptable range of cases. Five key simplifying assumptions have been made in our work: upstream and downstream boundaries are included in the control volume, with no temperature or species concentration gradients; there is uniform one-dimensional flow under adiabatic conditions; body forces, dissociation of products, and atomization of fuel are neglected; and only the normal shock relation is considered. In addition, although there are many variations in the molecular structures of these alternative fuels, consideration of such variability in the characteristics of the fuels is also neglected, and the analysis is based solely on the information given in the
Appendix A.
Figure 3 illustrates the different stages for the calculations below.
One-dimensional analysis with the variation of mass flux, initial temperature, and pressure is calculated from the conservation of mass and momentum. Thus the Rayleigh line yields the following relationship:
Combining conservation of mass, momentum, and energy with heat addition yields:
From the Rayleigh Line,
is:
Substituting into the Rankine-Hugonoit Curve yields:
Expanding and converting to the quadratic equation:
where
for every
are calculated accordingly. Next, detonation velocity in stoichiometric conditions and the gas-mixture properties at the shock front (state 2’) are estimated by applying a stoichiometric relation.
Every species mole and mass fraction are calculated, and thermochemical properties such as specific heat, gas constant, and specific heat ratio are obtained using these relations:
Heat formation,
q, is calculated using enthalpies-of-formation in the tabulated table, which is converted to a mass basis.
Detonation velocity and temperature at State 2 are determined using:
Using ideal-gas normal-shock and knowing the mixture specific-heat ratio and Mach number at the initial state, these relations are used to find State 2’:
is calculated using the conservation of mass. The State 2 Mach number should be equal to one (upper CJ point).