Oscillatory combustion in rockets Semiannual report, 1 Jun. - 30 Nov. 1968
Vaporization of droplets near critical point
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Vaporization of droplets near critical point
A three-dimensional analysis of combustion instabilities in liquid rocket engines is presented based on a mixed finite difference/spectral solution methodology for the gas phase and a discrete droplet tracking formulation for the liquid phase. Vaporization is treated by a simplified model based on an infinite thermal conductivitiy assumption for spherical liquid droplets of fuel in a convective environment undergoing transient heating. A simple two parameter phenomenological combustion response model is employed for validation of the results in the small amplitude regime. The computational procedure is demonstrated to capture the phenomena of wave propagation within the combustion chamber accurately. Results demonstrate excellent amplitude and phase agreement with analytical solutions for properly selected grid resolutions under both stable and unstable operating conditions. Computations utilizing the simplified droplet model demonstrate stable response to arbitrary pulsing. This is possibly due to the assumption of uniform droplet temperature which removes the thermal inertia time-lag response of the vaporization process. The mixed-character scheme is sufficiently efficient to allow solutions on workstations at a modest increase in computational time over that required for two-dimensional solutions.
The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.
Blending petroleum fuels with oxygenates is a common approach to stem the depletion of crude oil while also mitigating the impact of their combustion on the environment. It has recently been considered that butyl acetate (BA, C6H12O2 , boiling point of 399K) could be a viable oxygenate additive to diesel fuel. The platform used to study combustion of commercial and synthesized BA was an isolated droplet burning with spherical symmetry. The experiments were carried out in a drop tower to minimize buoyant convection. Restricting motion by supporting droplets on 14 micrometer fibers eliminated forced convection. Precise control of the initial diameter was achieved by first deploying large droplets and letting them evaporate to the size of interest, here being between 0.3 mm to 0.8 mm. The experimental results show that commercial and the newly synthesized BA burned almost identically for the evolution of droplet and flame diameters. Droplets burned slower as the initial diameter was increased and soot formation was noted for the larger sizes studied due to a residence size effect. Burning in a buoyant convective field showed higher burn rates due to the convective effect and that the asymmetric flames became nearly spherical as diameter decreased. The study illustrates the significant value that the microgravity environment brings to evaluating performance of oxygenated fuels synthesized by new production methods.
Hafnium droplets burning in ultrapure oxygen at various pressures, considering combustion induced natural mode periodic oscillations
Proposed test setup would measure combustion performance of new fuels by rapidly heating a droplet of coal/oil/water mixture and recording resulting explosion. Such mixtures are being considered as petroleum substitutes in oil-fired furnaces.
The combustion of electrosprays in a laminar counterflow diffusion flame has been experimentally studied by measuring droplet size and velocity distributions and gas-phase temperature. Detailed examination of the evolution of droplet size distribution as droplets approach the flame shows that, if substantial evaporation occurs before droplets 'interact' with the flame, the size distribution becomes bimodal. A secondary, sharp peak, in fact, develops in correspondence of diameters about one order of magnitude smaller than the mean. No evaporation mechanism can account for the development of such bimodality, that can be explained only in terms of a disintegration of droplets into finer fragments of size much smaller than the parent ones. This fission is of electric nature and it occurs when the repulsion of electric charges overcomes the surface tension cohesive force ultimately leading to a disintegration into finer fragments at or about the so-called Rayleigh limit. We here report on the first observation in combustion environments of such 'explosions'. If, on the other hand, droplets enter the very high temperature region before exploding, there appears to be no evidence of bimodality in their size distribution. In this case, in fact, flame chemi-ions may neutralize the charge on the droplets and thus prevent disruption.
A free jet configuration was chosen for measuring noncombusting spray fields and hydrocarbon-air spray flames in an effort to develop computational models of the dynamic interaction between droplets and the gas phase and to verify and refine numerical models of the entire spray combustion process. The development of a spray combustion facility is described including techniques for laser measurements in spray combustion environments and methods for data acquisition, processing, displaying, and interpretation.
An experimental and numerical study of combustion of a gasoline certification fuel (‘indolene’), and four (S4) and five (S5) component surrogates for it, is reported for the configurations of an isolated droplet burning with near spherical symmetry in the standard atmosphere, and a single cylinder engine designed for advanced compression ignition of pre-vaporized fuel. The intent was to compare performance of the surrogate for these different combustion configurations and to assess the broader applicability of the kinetic mechanism and property database for the simulations. A kinetic mechanism comprised of 297 species and 16797 reactions was used in the simulations that included soot chemistry and accounted for unsteady transport, liquid diffusion inside the droplet, radiative heat transfer, and variable properties. The droplet data showed a clear preference for the S5 surrogate in terms of burning rate. The simulations showed generally very good agreement with measured droplet, flame, and soot shell diameters. Measurements of combustion timing, in-cylinder pressure, and mass-averaged gas temperature were also well predicted with a slight preference for the S5 surrogate. Preferential vaporization was not evidenced from the evolution of droplet diameter but was clearly revealed in simulations of the evolution of mixture fractions inside the droplets. The influence of initial droplet diameter (Do) on droplet burning was strong, with S5 burning rates decreasing with increasing Do due to increasing radiation losses from the flame. Flame extinction was predicted for Do =3.0 mm as a radiative loss mechanism but not predicted for smaller Do for the conditions of the simulations.
The combustion of electrostatic sprays of heptane in laminar counterflow diffusion flames was experimentally studied by measuring droplet size and velocity distributions, as well as the gas-phase temperature. A detailed examination of the evolution of droplet size distribution as droplets approach the flame shows that, if substantial evaporation occurs before droplets interact with the flame, an initially monodisperse size distribution becomes bimodal. A secondary sharp peak in the size histogram develops in correspondence of diameters about one order of magnitude smaller than the mean. No evaporation mechanism can account for the development of such bimodality, that can be explained only in terms of a disintegration of droplets into finer fragments of size much smaller than that of the parent. Other evidence in support of this interpretation is offered by the measurements of droplet size-velocity correlation and velocity component distributions, showing that, as a consequence of the ejection process, the droplets responsible for the secondary peak have velocities uncorrelated with the mean flow. The fission is induced by the electric charge. When a droplet evaporates, in fact, the electric charge density on the droplet surface increases while the droplet shrinks, until the so-called Rayleigh limit is reached at which point the repulsion of electric charges overcomes the surface tension cohesive force, ultimately leading to a disintegraton into finer fragments. We report on the first observation of such fissions in combustion environments. If, on the other hand, insufficient evaporation has occurred before droplets enter the high temperature region, there appears to be no significant evidence of bimodality in their size distribution. In this case, in fact, the concentration of flame chemi-ions or, in the case of positively charged droplets, electrons may be sufficient for them to neutralize the charge on the droplets and to prevent disruption.
This presentation will describe results of an experimental and numerical study of combustion of a Tier II gasoline certification fuel and its surrogate blended with iso-butanol as a potential biofuel additive. The experiments were carried out under conditions where gas transport originated entirely by fuel evaporation with external convection (both forced and natural) were minimized. The resulting near spherically symmetric droplet flames were amenable to detailed numerical modeling. The droplet initial diameters were in the narrow range of 0.55 mm to 0.62 mm while and burning is performed under standard atmosphere. The combustion characteristics such as evolution of droplet, flame and soot shell diameters for indolene/iso-butanol and surrogate/iso-butanol mixtures are compared. The presence of iso-butanol is experimentally found to reduce the formation of soot as the iso-butanol concentration was increased. Simulations are also reported using a detailed kinetic chemistry comprised of 324 species and 17608 reactions using the CRECK POLIMI mechanism. The detailed numerical model included unsteady gas and liquid transport, soot chemistry and radiative transport. Droplet burn rates were well predicted while flame diameters showed more variability depending on how the flame was defined in the numerical simulation (peak gas temperature and OH concentration). Simulations also show the influence of iso-butanol concentration on flame temperature, greenhouse gas emissions and particulate emissions.
Fuel cells are a promising next-generation energy-conversion technology designed to replace internal combustion engines in transportation applications. However, much work remains to optimize them. Operation at high humidities causes liquid water droplet formation on Pt catalyst particles during oxygen reduction, potentially impeding reactant arrival to the reactive electrode. In this work, four different cases of water droplet growth in fuel-cell catalyst layers are considered: pinned or advancing droplets on a bare Pt surface, advancing droplets on a Nafion film, and water-layer growth in carbon nanopores. Transient drop growth is captured with a combination of mass, species mass, and momentum balances, and the subsequent limiting current is determined via oxygen diffusion and Tafel kinetics. Further, water droplets are found not to be mass-transfer limiting due to the relatively large liquid-gas area compared to the Pt nanoparticle. Mass-transfer-limited behavior is calculated in carbon nanopores.
X-ray phase contrast imaging is conducted at 90,517 Hz on a fuel spray to characterize the breakup of the bulk fuel into ligaments and individual droplets in a single-sector, swirl-stabilized combustor under reacting conditions. Phase contrast imaging is used to qualitatively assess spray breakup for two fuels: Jet-A, which represents a normal viscosity fuel, and a blend of JP-5 (64%) and farnesane (36% by volume), which represents a high-viscosity fuel, at three fuel-flow rates and two air inlet temperatures. Time-series images of the fuels show that atomization occurs much more rapidly at higher fuel pressure with thinner ligaments being rapidly shredded. The high-viscosity fuel is qualitatively observed to have markedly longer, thicker ligaments and larger droplet diameters after the primary breakup than the normal viscosity fuel. Using an image processing and droplet quantification strategy, the mean droplet velocity magnitudes were found to decrease with decreasing fuel pressure and fuel-flow rate. The increase in preheat temperature from 323 to 370 K results in an increase in the mean droplet velocity magnitude and a decrease in numerical average and Sauter mean diameter, indicating improved atomization. The high-viscosity fuel was found to have larger average droplet diameters at each condition than the normal viscosity fuel. The increased preheat is found to enhance the combustion efficiency, which is hypothesized to couple with the improved atomization due to the higher combustion temperatures. This indicates the importance of examining spray atomization and breakup under realistic combusting conditions for relevant fuels as the heat release directly couples with the atomization process.
Many thermoplastic materials in common use for a wide range of applications, including spacecraft, develop bubbles internally as they burn due to chemical reactions taking place within the bulk. These bubbles grow and migrate until they burst at the surface, forceably ejecting volatile gases and, occasionally, molten fuel. In experiments in normal gravity, Kashiwagi and Ohlemiller observed vapor jets extending a few centimeters from the surface of a radiatively heated polymethylmethacrylate (PMMA) sample, with some molten material ejected into the gas phase. These physical phenomena complicated the combustion process considerably. In addition to the non-steady release of volatiles, the depth of the surface layer affected by oxygen was increased, attributed to the roughening of the surface by bursting events. The ejection of burning droplets in random directions presents a potential fire hazard unique to microgravity. In microgravity combustion experiments on nylon Velcro fasteners and on polyethylene wire insulation, the presence of bursting fuel vapor bubbles was associated with the ejection of small particles of molten fuel as well as pulsations of the flame. For the nylon fasteners, particle velocities were higher than 30 cm/sec. The droplets burned robustly until all fuel was consumed, demonstrating the potential for the spread of fire in random directions over an extended distance. The sequence of events for a bursting bubble has been photographed by Newitt et al.. As the bubble reaches the fluid surface, the outer surface forms a dome while the internal bubble pressure maintains a depression at the inner interface. Liquid drains from the dome until it breaks into a cloud of droplets on the order of a few microns in size. The bubble gases are released rapidly, generating vortices in the quiescent surroundings and transporting the tiny droplets. The depression left by the escaping gases collapses into a central jet, which rises with a high velocity and may break up, releasing one or more relatively large drops (on the order of a millimeter in these experiments). A better understanding of bubble development and bursting processes, the effects of bursting behavior on burning rate of the bulk material, and the circumstances under which large droplets are expelled, as well as their trajectories, sizes, and burning rates, is sought through computer modeling compared with experiment.
The Variable Thrust Engine (VTE), developed by TRW, for the Orbit Maneuvering Vehicle (OMV) uses a hypergolic propellant combination of Monomethyl Hydrazine (MMH) and Nitrogen Tetroxide (NTO) as fuel and oxidizer, respectively. The propellants are pressure fed into the combustion chamber through a single pintle injection element. The performance of this engine is dependent on the pintle geometry and a number of complex physical phenomena and their mutual interactions. The most important among these are (1) atomization of the liquid jets into fine droplets; (2) the motion of these droplets in the gas field; (3) vaporization of the droplets (4) turbulent mixing of the fuel and oxidizer; and (5) hypergolic reaction between MMH and NTO. Each of the above phenomena by itself poses a considerable challenge to the technical community. In a reactive flow field of the kind occurring inside the VTE, the mutual interactions between these physical processes tend to further complicate the analysis. The objective of this work is to develop a comprehensive mathematical modeling methodology to analyze the flow field within the VTE. Using this model, the effect of flow parameters on various physical processes such as atomization, spray dynamics, combustion, and radiation is studied. This information can then be used to optimize design parameters and thus improve the performance of the engine. The REFLEQS CFD Code is used for solving the fluid dynamic equations. The spray dynamics is modeled using the Eulerian-Lagrangian approach. The discrete ordinate method with 12 ordinate directions is used to predict the radiative heat transfer in the OMV combustion chamber, nozzle, and the heat shield. The hypergolic reaction between MMH and NTO is predicted using an equilibrium chemistry model with 13 species. The results indicate that mixing and combustion is very sensitive to the droplet size. Smaller droplets evaporate faster than bigger droplets, leading to a well mixed zone in the combustion chamber. The radiative heat flux at combustion chamber and nozzle walls are an order of negligible less than the conductive heat flux. Simulations performed with the heat shield show that a negligible amount of fluid is entrained into the heat shield region. However, the heat shield is shown to be effective in protecting the OMV structure surrounding the engine from the radiated heat.
Research designed to develop detailed knowledge of the secondary atomization and ignition characteristics of aluminum slurry propellants was started. These processes are studied because they are the controlling factors limiting the combustion efficiency of aluminum slurry propellants in rocket applications. A burner and spray rig system allowing the study of individual slurry droplets having diameters from about 10 to 100 microns was designed and fabricated. The burner generates a near uniform high temperature environment from the merging of 72 small laminar diffusion flames above a honeycomb matrix. This design permits essentially adiabatic operation over a wide range of stoichiometries without danger of flashback. A single particle sizing system and velocimeter also were designed and assembled. Light scattered from a focused laser beam is related to the particle (droplet) size, while the particle velocity is determined by its transit time through the focal volume. Light from the combustion of aluminum is also sensed to determine if ignition was achieved. These size and velocity measurements will allow the determination of disruption and ignition times as functions of drop sizes and ambient conditions.
Characterization of liquid rocket engine injection elements is an important part of the development process for rocket engine combustion devices. Modern nonintrusive instrumentation for flow velocity and spray droplet size measurement and automated, computer-controlled test facilities allow rapid, low-cost evaluation of injection element performance and behavior. Application of these methods in rocket engine development, paralleling their use in gas turbine engine development, will reduce rocket engine development cost and risk. Two types of liquid rocket engine injection elements have been characterized using such methods.
Characterization of liquid rocket engine injector elements is an important part of the development process for rocket engine combustion devices. Modern nonintrusive instrumentation for flow velocity and spray droplet size measurement, and automated, computer-controlled test facilities allow rapid, low-cost evaluation of injector element performance and behavior. Application of these methods in rocket engine development, paralleling their use in gas turbine engine development, will reduce rocket engine development cost and risk. The Alternate Turbopump (ATP) Hot Gas Systems (HGS) preburner injector elements were characterized using such methods, and the methodology and some of the results obtained will be shown.