Search NASA⌕ Search

SEARCH · Search NASA

Results for “combustion dynamics”

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.

At least 757 records · Page 42

Analysis of Hydrodynamic (Landau) Instability in Liquid-Propellant Combustion at Normal and Reduced Gravity

The burning of liquid propellants is a fundamental combustion problem that is applicable to various types of propulsion and energetic systems. The deflagration process is often rather complex, with vaporization and pyrolysis occurring at the liquid/gas interface and distributed combustion occurring either in the gas phase or in a spray. Nonetheless, there are realistic limiting cases in which combustion may be approximated by an overall reaction at the liquid/gas interface. In one such limit, the gas flame occurs under near-breakaway conditions, exerting little thermal or hydrodynamic influence on the burning propellant. In another such limit, distributed combustion occurs in an intrusive regime, the reaction zone lying closer to the liquid/gas interface than the length scale of any disturbance of interest. Finally, the liquid propellant may simply undergo exothermic decomposition at the surface without any significant distributed combustion, such as appears to occur in some types of HydroxylAmmonium Nitrate (HAN)-based liquid propellants at low pressures. Such limiting models have recently been formulated,thereby significantly generalizing earlier classical models that were originally introduced to study the hydrodynamic stability of a reactive liquid/gas interface. In all of these investigations, gravity appears explicitly and plays a significant role, along with surface tension, viscosity, and, in the more recent models, certain reaction-rate parameters associated with the pressure and temperature sensitivities of the reaction itself. In particular, these parameters determine the stability of the deflagration with respect to not only classical hydrodynamic disturbances, but also with respect to reactive/diffusive influences as well. Indeed, the inverse Froude number, representing the ratio of buoyant to inertial forces, appears explicitly in all of these models, and consequently, in the dispersion relation that determines the neutral stability boundaries beyond which steady, planar burning is unstable to nonsteady, and/or nonplanar (cellular) modes of burning. These instabilities thus lead to a number of interesting phenomena, such as the sloshing type of waves that have been observed in mixtures of HAN and TriEthanolAmmonium Nitrate (TEAN) with water. Although the Froude number was treated as an O(1) quantity in these studies, the limit of small inverse Froude number corresponding to the microgravity regime is increasingly of interest and can be treated explicitly, leading to various limiting forms of the models, the neutral stability boundaries, and, ultimately, the evolution equations that govern the nonlinear dynamics of the propagating reaction front. In the present work, we formally exploit this limiting parameter regime to compare some of the features of hydrodynamic instability of liquid-propellant combustion at reduced gravity with the same phenomenon at normal gravity.

Margolis, Stephen B.↗

Supersonic Free-Jet Combustion in a Ramjet Burner

A new dual-mode ramjet combustor concept intended for operation over a wide flight Mach number range is described. Subsonic combustion mode is similar to that of a traditional ram combustor which allows operation at higher efficiency, and to lower flight Mach numbers than current dual-mode scramjets. High speed mode is characterized by supersonic combustion in a free-jet that traverses the subsonic combustion chamber to a variable nozzle. The maximum flight Mach number of this scheme is governed largely by the same physics as its classical counterpart. Although a variable combustor exit aperture is required, the need for fuel staging to accommodate the combustion process is eliminated. Local heating from shock-boundary-layer interactions on combustor walls is also eliminated. Given the parallel nature of the present scheme, overall flowpath length is less than that of present dual-mode configurations. Cycle analysis was done to define the flowpath geometry for computational fluid dynamics (CFD) analysis, and then to determine performance based on the CFD results. CFD results for Mach 5, 8, and 12 flight conditions indicate stable supersonic free-jet formation and nozzle reattachment, thereby establishing the basic feasibility of the concept. These results also reveal the structure of, and interactions between the free-jet and recirculating combustion chamber flows. Performance based on these CFD results is slightly less than that of the constant-pressure-combustion cycle analysis primarily due to these interactions. These differences are quantified and discussed. Additional CFD results at the Mach 8 flight condition show the effects of nozzle throat area variation on combustion chamber pressure, flow structure, and performance. Calculations with constant temperature walls were also done to evaluate heat flux and overall heat loads. Aspects of the concept that warrant further study are outlined. These include diffuser design, ramjet operation, mode transition, loss mechanisms, and the effects of secondary flow for wall cooling and combustion chamber pressurization. Also recommended is an examination of system-level aspects such as weight, thermal management and rocket integration as well as alternate geometries and variable geometry schemes.

Trefny, Charles J.↗

Combustion Stability of the Gas Generator Assembly from J-2X Engine E10001 and Powerpack Tests

Testing of a powerpack configuration (turbomachinery and gas generator assembly) and the first complete engine system of the liquid oxygen/liquid hydrogen propellant J-2X rocket engine have been completed at the NASA Stennis Space Center. The combustion stability characteristics of the gas generator assemblies on these two systems are of interest for reporting since considerable effort was expended to eliminate combustion instability during early development of the gas generator assembly with workhorse hardware. Comparing the final workhorse gas generator assembly development test data to the powerpack and engine system test data provides an opportunity to investigate how the nearly identical configurations of gas generator assemblies operate with two very different propellant supply systems one the autonomous pressure-fed test configuration on the workhorse development test stand, the other the pump-fed configurations on the powerpack and engine systems. The development of the gas generator assembly and the elimination of the combustion instability on the pressure-fed workhorse test stand have been reported extensively in the two previous Liquid Propulsion Subcommittee meetings 1-7. The powerpack and engine system testing have been conducted from mid-2011 through 2012. All tests of the powerpack and engine system gas generator systems to date have been stable. However, measureable dynamic behavior, similar to that observed on the pressure-fed test stand and reported in Ref. [6] and attributed to an injection-coupled response, has appeared in both powerpack and engine system tests. As discussed in Ref. [6], these injection-coupled responses are influenced by the interaction of the combustion chamber with a branch pipe in the hot gas duct that supplies gaseous helium to pre-spin the turbine during the start transient. This paper presents the powerpack and engine system gas generator test data, compares these data to the development test data, and provides additional combustion stability analyses of the configurations.

Hulka, J. R.↗

Vortical interactions in turbulent thermoacoustic systems

This study examines the dynamics of vortical interactions and their implications for mitigating thermoacoustic instability in a turbulent combustor. The regions of intense vortical interactions are identified as vortical communities in the network space of weighted directed vortical networks constructed from two-dimensional experimental velocity data. One can expect vortical interactions in the combustor to be strongest near the moment of vortex shedding, as the shed vortices gradually weaken due to dissipation while convecting downstream. However, we show that, during the state of thermoacoustic instability, there is a non-trivial consistent phase lag of approximately 52° between the shedding of the coherent structures from the backward-facing step and the time instant when the vortical interactions attain their local maximum value. We explain this phase lag by investigating the correlation between acoustic pressure fluctuations, spatio-temporal dynamics of coherent structures and vortical interactions in the reaction field of the combustor. We also show the aperiodic variation of vortical interactions during the states of combustion noise and aperiodic epochs of intermittency. Furthermore, the spatio-temporal evolution of pairs of vortical communities with the maximum inter-community interactions provides insight into explaining the critical regions detected in the reaction field during the states of intermittency and thermoacoustic instability, also identified in previous studies. As a result, we further show that the most efficient suppression of thermoacoustic instability via air microjet injection is achieved when steady air jets are introduced to disrupt the maximum inter-community interactions present during the state of thermoacoustic instability.

Sahay, Ankit [Indian Institute of Technology Madra↗

Mass-loading and parallel magnetized shocks

Recent observations at comets Giacobini-Zinner and Halley suggest that simple nonreacting gas dynamics or MHD is an inappropriate description for the bow shock. The thickness of the observed (sub)shock implies that mass-loading is an important dynamical process within the shock itself, thereby requiring that the Rankine-Hugoniot conditions possess source terms. This leads to shocks with properties similar to those of combustion shocks. The paper considers parallel magnetized shocks subjected to mass-loading, describes some properties which distinguish them from classical MHD parallel shocks, and establishes the existence of a new kind of MHD compound shock. These results will be of importance both to observations and numerical simulations of the comet-solar wind interaction.

Zank, G. P.↗

Computational Fluid Dynamic Modeling of Rocket Based Combined Cycle Engine Flowfields

Computational Fluid Dynamic techniques are used to study the flowfield of a fixed geometry Rocket Based Combined Cycle engine operating in rocket ejector mode. Heat addition resulting from the combustion of injected fuel causes the subsonic engine flow to choke and go supersonic in the slightly divergent combustor-mixer section. Reacting flow computations are undertaken to predict the characteristics of solutions where the heat addition is determined by the flowfield. Here, adaptive gridding is used to improve resolution in the shear layers. Results show that the sonic speed is reached in the unheated portions of the flow first, while the heated portions become supersonic later. Comparison with results from another code show reasonable agreement. The coupled solutions show that the character of the combustion-based thermal choking phenomenon can be controlled reasonably well such that there is opportunity to optimize the length and expansion ratio of the combustor-mixer.

Daines, Russell L.↗

Beagle 2: The Next Exobiology Mission to Mars

Beagle 2 is a 60 kg probe (with a 30 kg lander) developed in the United Kingdom for inclusion on the European Space Agency's 2003 Mars Express. Beagle 2 will deliver to the Martian surface a payload which consists of a high percentage of science instruments to landed spacecraft mass. Beagle 2 will be launched in June, 2003 with Mars Express on a Soyuz-Fregat rocket from the Baikonur Cosmodrome in Kazakhstan. Beagle 2 will land on Mars on December 26, 2003 in the Isidis Planitia basin (approximately 10 degrees N and 275 degrees W), a large sedimentary basin that overlies the boundary between ancient highlands and northern plains. Isidis Planitia, the third largest basin on Mars, which is possibly filled with sediment deposited at the bottom of long-standing lakes or seas, offers an ideal environment for preserving traces of life. Beagle 2 was developed to search for organic material and other volatiles on and below the surface of Mars in addition to the study of the inorganic chemistry and mineralogy. Beagle 2 will utilize a mechanical mole and grinder to obtain samples from below the surface, under rocks and inside rocks. A pair of stereo cameras will image the landing site along with a microscope for examination of surface and rock samples. Analyses will include both rock and soil samples at various wavelengths, X-ray spectrometer and Mossbauer spectrometer as well as a search for organics and other light element species (e.g. carbonates and water) and measurement of their isotopic compositions. Beagle 2 has as its focus the goal of establishing whether evidence for life existed in the past on Mars at the Isidis Planitia site or at least establishing if the conditions were ever suitable. Carbonates and organic components were first recognized as existing on Mars when they were found in the Martian meteorite Nakhla. Romanek et al showed the carbonates in ALH84001 were formed at low temperatures. McKay et al noted possible evidence of early life on Mars within the ALH84001 meteorite. Thomas-Keprta et al showed the magnetite biomarkers in ALH84001's carbonates are indistinguishable from those formed by magnetotactic bacteria found on Earth. Gibson et al showed there was significant evidence for liquid water and biogenic products present on Mars across a 3.9 billion year period. A mechanical arm (PAW) operates from the lander and is used for science operations along with sample acquisition. Instruments attached to the PAW include: stereo cameras, Mossbauer instrument, X-ray fluorescence instrument, microscope, environmental sensors, rock corer/grinder, a spoon, mirror, brushes, a mole attachment for acquisition of subsurface to depths of 1 to 2 meters and an illumination device. Each camera has 14 filters which have been optimized for mineralogy composition, dust and water vapor detection. The microscope's camera is designed for viewing the size and shape of dust particles, rock surfaces, microfossils, and characteristics of the samples prior to introduction into the gas analysis package (GAP). The camera has a resolution of 4 microns/pixel, features four color capability (red, green, blue and UV (ultraviolet) fluoresence), a depth of focus of 40 micrometers and translation stage of +3 millimeters. The heart of the Beagle 2's life detection package is the gas analysis package (GAP), which consists of a mass spectrometer with collectors at fixed masses for precise isotopic ratio measurements and voltage scanning for spectral analysis. Primary aim of the GAP is to search for the presence of bulk constituents, individual species, and isotopic fractionations for both extinct and extent life along with studying the low-temperature geochemistry of the hydrogen, carbon, nitrogen and oxygen components on Mars from both the surface and atmosphere. GAP is a magnetic sector mass spectrometer with the range of 1 to 140 amu which can be operated in both the static and dynamic modes. A triple Faraday collector array will be used for C, N and O ratios along with a double Faraday array for H/D. Pulse counting electron multiplier will be utilized for noble gases and selected organics. Anticipated detection limits are at the picomole level for operation in the static mode of operation and high precision isotopic measurements will be made in the dynamic mode. Sample processing and preparation system consists of reaction vessels along with references. Sample ovens capable of being heated are attached to the manifold for sample combustion. Surface, subsurface materials and interior rock specimens will be combusted in pure oxygen gas at various temperature intervals to release organic matter and volatiles. Combustion process will permit detection of all forms and all atoms of carbon present in the samples. A chemical processing system is capable of a variety of conversion reactions. Gases are manipulated either by cryogenic or chemical reactions and passed through the gas handling portion of the vacuum system. There are two modes of operation: quantitative analysis and precise isotopic measurements. Additional information is contained in the original extended abstract.

Gibson, Everett K., Jr.↗

Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations

The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing nu- merical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten (1978) but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch to the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992) used by Gerritsen and Olsson (1996) in an adaptive mesh refinement method, to determine regions where re nement should be done. The other is the modification of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand alone options for numerical algorithm other than the Yee et al. scheme.

Sjoegreen, B.↗

Analytical Redundancy Using Kalman Filters for Rocket Engine Sensor Validation

The use of sensor redundancy is crucial in aerospace systems to maintain safe, reliable operation. While hardware redundancy is more common in application, analytical redundancy can provide a viable alternative in systems where the installation of multiple redundant sensors is not viable. To this end, the use of Kalman filters to analytically validate sensor measurements within rocket engines was explored. First, a dynamic model of the RS 25 engine, a derivative of the Space Shuttle Main Engine (SSME), was reduced to a subset of relations, focused around the main combustion chamber pressure. These relations were used within the Kalman filter algorithm to generate an estimate of sensor measurements to be compared with true measurements for data validation purposes. By using a bank of Kalman filters, the residuals between the estimated and true measurements were used to detect and isolate sensor faults. Through fault simulations, the sensor validation performance of this Kalman filter bank design was compared to a hardware redundancy check. Sensor bias and drift faults of various magnitudes were injected into nominal RS 25 engine test data. Results for both approaches show comparable fault detection with most bias faults found nearly instantaneously by both algorithms. Drift fault detection results show certain cases where one algorithm is faster than the other. The key advantage of the Kalman filter algorithm is shown in fault isolation performance where it can isolate faults between two redundant sensors while the hardware redundancy comparisons cannot.

sensors↗

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗

Splashing Droplets

Current data on droplet breakup is scarce for the sizes and velocities typical of practical applications such as in spray combustion processes and coating processes. While much more representative of practical applications, the small spatial scales and rapid time-scales prevent detailed measurement of the internal fluid dynamics and liquid property gradients produced by impinging upon surfaces. Realized through the extended spatial and temporal scales afforded by a microgravity environment, an improved understanding of drop breakup dynamics is sought to understand and ultimately control the impingement dynamics of droplets upon surfaces in practical situations. The primary objective of this research will be to mark the onset of different 'splashing modes' and to determine their temperature, pressure and angle dependence for impinging droplets representative of practical fluids. In addition, we are modeling the evolution of droplets that do not initially splash but rather undergo a 'fingering' evolution observed on the spreading fluid front and the transformation of these fingers into splashed products. An example of our experimental data is presented below. These images are of Isopar V impacting a mirror-polished surface. They were acquired using a high-speed camera at 1000 frames per second. They show the spreading of a single droplet after impact and ensuing finger instabilities. Normal gravity experimental data such as this will guide low gravity measurements in the 2.2 second drop tower and KC-135 aircraft as available. Presently we are in the process of comparing the experimental data of droplet shape evolution to numerical models, which can also capture the internal fluid dynamics and liquid property gradients such as produced by impingement upon a heated surface. To-date isothermal numerical data has been modeled using direct numerical simulations of representative splashing droplets. The data obtained so far indicates that the present model describes well the droplet wall interactions to a point in time just before splash. Additional information is included in the original extended abstract.

VanderWal, Randall L.↗

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗

Simulation of Unsteady Hypersonic Combustion Around Projectiles in an Expansion Tube

The temporal evolution of combustion flowfields established by the interaction between wedge-shaped bodies and explosive hydrogen-oxygen-nitrogen mixtures accelerated to hypersonic speeds in an expansion tube is investigated. The analysis is carried out using a fully implicit, time-accurate, computational fluid dynamics code that we developed recently for solving the Navier-Stokes equations for a chemically reacting gas mixture. The numerical results are compared with experimental data from the Stanford University expansion tube for two different gas mixtures at Mach numbers of 4.2 and 5.2. The experimental work showed that flow unstart occurred for the Mach 4.2 cases. These results are reproduced by our numerical simulations and, more significantly, the causes for unstart are explained. For the Mach 5.2 mixtures, the experiments and numerical simulations both produced stable combustion. However, the computations indicate that in one case the experimental data were obtained during the transient phase of the flow; that is, before steady state had been attained.

Yungster, S.↗

Modeling of Transient Flow Mixing of Streams Injected into a Mixing Chamber

Ignition is recognized as one the critical drivers in the reliability of multiple-start rocket engines. Residual combustion products from previous engine operation can condense on valves and related structures thereby creating difficulties for subsequent starting procedures. Alternative ignition methods that require fewer valves can mitigate the valve reliability problem, but require improved understanding of the spatial and temporal propellant distribution in the pre-ignition chamber. Current design tools based mainly on one-dimensional analysis and empirical models cannot predict local details of the injection and ignition processes. The goal of this work is to evaluate the capability of the modern computational fluid dynamics (CFD) tools in predicting the transient flow mixing in pre-ignition environment by comparing the results with the experimental data. This study is a part of a program to improve analytical methods and methodologies to analyze reliability and durability of combustion devices. In the present paper we describe a series of detailed computational simulations of the unsteady mixing events as the cold propellants are first introduced into the chamber as a first step in providing this necessary environmental description. The present computational modeling represents a complement to parallel experimental simulations' and includes comparisons with experimental results from that effort. A large number of rocket engine ignition studies has been previously reported. Here we limit our discussion to the work discussed in Refs. 2, 3 and 4 which is both similar to and different from the present approach. The similarities arise from the fact that both efforts involve detailed experimental/computational simulations of the ignition problem. The differences arise from the underlying philosophy of the two endeavors. The approach in Refs. 2 to 4 is a classical ignition study in which the focus is on the response of a propellant mixture to an ignition source, with emphasis on the level of energy needed for ignition and the ensuing flame propagation issues. Our focus in the present paper is on identifying the unsteady mixing processes that provide the propellant mixture in which the ignition source is to be placed. In particular, we wish to characterize the spatial and temporal mixture distribution with a view toward identifying preferred spatial and temporal locations for the ignition source. As such, the present work is limited to cold flow (pre-ignition) conditions

Voytovych, Dmytro M.↗

Dynamic projection of anthropogenic emissions in China: methodology and 2015–2050 emission pathways under a range of socio-economic, climate policy, and pollution control scenarios

Abstract. Future trends in air pollution and greenhouse gas (GHG)emissions for China are of great concern to the community. A set of globalscenarios regarding future socio-economic and climate developments, combiningshared socio-economic pathways (SSPs) with climate forcing outcomes asdescribed by the Representative Concentration Pathways (RCPs), was createdby the Intergovernmental Panel on Climate Change (IPCC). Chinese researchers have also developed various emission scenarios by considering detailed local environmental and climate policies. However, a comprehensive scenario set connecting SSP–RCP scenarios with local policies and representing dynamic emission changes under local policies is still missing. In this work, to fill this gap, we developed a dynamic projection model, the Dynamic Projection model for Emissions in China (DPEC), to explore China'sfuture anthropogenic emission pathways. The DPEC is designed tointegrate the energy system model, emission inventory model, dynamicprojection model, and parameterized scheme of Chinese policies. The modelcontains two main modules, an energy-model-driven activity rate projectionmodule and a sector-based emission projection module. The activity rateprojection module provides the standardized and unified future energyscenarios after reorganizing and refining the outputs from the energy systemmodel. Here we use a new China-focused version of the Global ChangeAssessment Model (GCAM-China) to project future energy demand and supply inChina under different SSP–RCP scenarios at the provincial level. Theemission projection module links a bottom-up emission inventory model, theMulti-resolution Emission Inventory for China (MEIC), to GCAM-China andaccurately tracks the evolution of future combustion and production technologiesand control measures under different environmental policies. We developedtechnology-based turnover models for several key emitting sectors (e.g.coal-fired power plants, key industries, and on-road transportationsectors), which can simulate the dynamic changes in the unit/vehicle fleetturnover process by tracking the lifespan of each unit/vehicle on an annualbasis. With the integrated modelling framework, we connected five SSP scenarios(SSP1–5), five RCP scenarios (RCP8.5, 7.0, 6.0, 4.5, and 2.6), and threepollution control scenarios (business as usual, BAU; enhanced controlpolicy, ECP; and best health effect, BHE) to produce six combined emissionscenarios. With those scenarios, we presented a wide range of China's futureemissions to 2050 under different development and policy pathways. We foundthat, with a combination of strong low-carbon policy and air pollutioncontrol policy (i.e. SSP1-26-BHE scenario), emissions of major airpollutants (i.e. SO 2 , NO x , PM 2.5 , and non-methane volatile organic compounds – NMVOCs) in China willbe reduced by 34%–66% in 2030 and 58%–87% in 2050 compared to 2015. End-of-pipe control measures are more effective for reducing air pollutant emissions before 2030, while low-carbon policy will play a more important rolein continuous emission reduction until 2050. In contrast, China's emissionswill remain at a high level until 2050 under a reference scenario without activeactions (i.e. SSP3-70-BAU). Compared to similar scenarios set from theCMIP6 (Coupled Model Intercomparison Project Phase 6), our estimates ofemission ranges are much lower than the estimates from the harmonized CMIP6 emissions dataset in2020–2030, but their emission ranges become similar in the year 2050.

54 ENVIRONMENTAL SCIENCES↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion at Normal and Reduced Gravity

The classical Landau/Levich models of liquid-propellant combustion, despite their relative simplicity, serve as seminal examples that correctly describe the onset of hydrodynamic instability in reactive systems. Recently, these two separate models have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed either numerically or analytically in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity and other parameters on the hydrodynamic instability of the propagating liquid/gas interface. In particular, an analytical expression is derived for the neutral stability boundary A(p)(k), where A(p) is the pressure sensitivity of the burning rate and k is the wavenumber of the disturbance. The results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for negative values of A(p). In the limiting case of weak gravity, it is shown that hydrodynamic instability in liquid-propellant combustion is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumber disturbances. it is also demonstrated that, in general, surface tension and the viscosity of both the liquid and gas phases each produce comparable stabilizing effects in the long-wavenumber regime, thereby providing important modifications to previous analyses in which one or more of these effects were neglected.

Margolis, S. B.↗

Pre-Ionization Controlled Laser Plasma Formation for Ignition Applications (Final Scientific/Technical Report)

This report consists of three grand funded papers. Abstracts for each follow. 1) In this work, we present experimental and modeling studies of air pre-ionization using ultraviolet (UV) laser pulses and its effect on laser breakdown of an overlapped near-infrared (NIR) pulse. Experimental studies are conducted with a 266 nm beam (fourth harmonic of Nd:YAG) for UV pre-ionization and an overlapped 1064 nm NIR beam (fundamental of Nd:YAG), both having pulse duration of 10 ns. Results show that the UV beam produces a pre-ionized volume which assists in breakdown of the NIR beam, leading to reduction in NIR breakdown threshold by factor of >2. Numerical modeling is performed to examine the ionization and breakdown of both beams. The modeled breakdown threshold of the NIR, including assist by pre-ionization, is in reasonable agreement with the experimental results. 2) The present contribution compares the energy absorption, optical emission, temperature, and fluid dynamics of ultraviolet (UV) λ = 266 nm and near infrared (NIR) λ = 1064 nm nanosecond laser induced plasmas in ambient air. For UV pulses at the conditions studied, energy absorption by the plasmas increases relatively gradually with laser pulse energy starting at delivered energy of E~8 mJ. Corresponding measurements of plasma luminosity show that the absorption of UV radiation does not necessarily result in visible plasma emission. For the NIR induced plasmas, the energy absorption profile is far more abrupt and begins at ~55 mJ. In contrast with UV, the absorption of NIR radiation is always accompanied by intense optical emission. The temperatures of both types of plasma have been measured with Rayleigh scattering thermometry (at times after the Thomson signal sufficiently diminishes). The UV plasmas can attain a wider range of temperatures, including lower temperatures, depending on the pulse energy (e.g., T ~400–2000K for E ~7–35 mJ at Δt = 10 μs after the pulse) while the NIR plasmas show only hotter temperatures (e.g., T ~ 12 000 K for E = 75 mJ at Δt = 10 μs after the pulse) as is consistent with the literature. Differences in the fluid dynamics for UV versus NIR pulses are shown with Schlieren imaging. In this work, the contrast in the UV and NIR plasma threshold behavior is attributed to differing roles of avalanche ionization and multiphoton ionization as is also illustrated by a simple numerical model. 3) The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

Mechanical engineering↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion at Normal and Reduced Gravity

The classical Landau/Levich models of liquid-propellant combustion, despite their relative simplicity, serve as seminal examples that correctly describe the onset of hydrodynamic instability in reactive systems. Recently, these two separate models have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed either numerically or analytically in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity and other parameters on the hydrodynamic instability of the propagating liquid/gas interface. In particular, an analytical expression is derived for the neutral stability boundary A(sub p)(k), where A(sub p) is the pressure sensitivity of the burning rate and k is the wavenumber of the disturbance. The results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limiting case of weak gravity, it is shown that hydrodynamic instability in liquid-propellant combustion is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumber disturbances. It is also demonstrated that, in general, surface tension and the viscosity of both the liquid and gas phases each produce comparable stabilizing effects in the large-wavenumber regime, thereby providing important modifications to previous analyses in which one or more of these effects were neglected.

Margolis, Stephen B.↗