SEARCH · Engineering Papers
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.
All of the elements used in the Reacting and Multi-Phase (RAMP2) computer code are described in detail. The code can be used to model the dominant phenomena which affect the prediction of liquid and solid rocket nozzle and orbital plume flow fields.
We report field and laboratory observations to date indicate that the efficiency of hydraulic fracturing, as it relates to hydrocarbon recovery, depends significantly on geochemical alterations to rock surfaces that diminish accessibility by partial or total plugging of the pore and fracture networks. This is caused by mineral scale deposition such as coating of fracture surfaces with precipitates, particle migration, and/or crack closure due to dissolution under stress. In reactive flow-through experiments, mineral reactions in response to acidic fluid injection significantly reduced system porosity and core permeability. The present study focuses on changes to fluid chemistry and shale surfaces (inlet and fracture walls) resulting from shale-fluid interactions and integrating these findings for an improved estimate of transport-related consequences. The reacted shale surfaces were examined by spatially-resolved scanning electron microscopy - energy dispersive spectroscopy (SEM-EDS) analysis. Importantly, inductively coupled plasma - mass spectrometry/optical emission spectroscopy (ICP-MS/OES) was utilized to probe the chemical evolution of the core-flood effluents. The three study cores selected from the Marcellus formation represent different mineralogies and structural features. In flow-through experiments, lab-generated brine and HCl-based fracture fluid (pH=2) were injected sequentially under effective stress (up to 500 psi) at reservoir temperature (80°C). SEM-EDS results confirmed by the ICP concentration trends showed significant Fe hydroxide precipitates in clay- and pyrite-rich outcrop samples due to partial oxidation of Fe-bearing phases in the case of intrusion of low salinity water-based fluids. Porosity reduction in the MSEEL (Marcellus Shale Energy and Environmental Laboratory) carbonate-rich sample is related to compaction of cores under stress due to matrix softening with substantial dissolution, and pore-filling by hydroxides, as well as barite and salts. Despite the same fluid compositions and experimental conditions used for both MSEEL samples, barite precipitation was much more intense in the MSEEL clay-rich sample due to its greater sorption capacity and additional sulfate source as well as fissile nature with multiple lengthwise cracks. ICP tests revealed time-resolved concentration trends in produced brine and reactive fluids that in turn complemented the pre-/post-reaction SEM-EDS observations.
A review of chemically reactive flow analysis is presented. Key studies are reviewed. The approaches to modeling chemically reactive flows that are discussed are PDF, conserved PDF, conserved scalar, stochastic, and mixing length. Turbulent mixing mechanisms are also discussed and a model is formulated that could be used for optimizing non-premixed gas reacting systems. Consequently, a re-examination of the NO(x) reduction concepts for nonpremixed reactants are suggested together with two new designs for low NO(x) subsonic combustor research.
A numerical study was performed to assess the effects of vitiated air on the chemical kinetics of hydrogen, ethane, and methane combustion with air. A series of calculations in static reacting systems was performed, where the initial temperature was specified and reactions occurred at constant pressure. Three different types of test flow contaminants were considered: NP, H2O, and a combustion of H2O and CO2. These contaminants are present in the test flows of facilities used for hypersonic propulsion testing. The results were computed using a detailed reaction mechanism and are presented in terms of ignition and reaction times. Calculations were made for a wide range of contaminant concentrations, temperatures and pressures. The results indicate a pronounced kinetic effect over a range of temperatures, especially with NO contamination and, to a lesser degree, with H2O contamination. In all cases studied, CO2 remained kinetically inert, but had a thermodynamic effect on results by acting as a third body. The largest effect is observed with combustion using hydrogen fuel, less effect is seen with combustion of ethane, and little effect of contaminants is shown with methane combustion.
A numerical study was performed to assess the effects of vitiated air on the chemical kinetics of hydrogen, ethane, and methane combustion with air. A series of calculations in static reacting systems was performed, where the initial temperature was specified and reactions occurred at constant pressure. Three different types of test flow contaminants were considered: NO, H2O, and a combination of H2O and CO2. These contaminants are present in the test flows of facilities used for hypersonic propulsion testing. The results were computed using a detailed reaction mechanism and are presented in terms of ignition and reaction times. Calculations were made for a wide range of contaminant concentrations, temperatures and pressures. The results indicate a pronounced kinetic effect over a range of temperatures, especially with NO contamination and, to a lesser degree, with H2O contamination. In all cases studied, CO2 remained kinetically inert, but had a thermodynamically effect on results by acting as a third body. The largest effect is observed with combustion using hydrogen fuel, less effect is seen with combustion of ethane, and little effect of contaminants is shown with methane combustion.
Work on understanding the mechanisms of mixing and reaction in high speed turbulent reacting flows was continued. Efforts, in particular, were concentrated on taking advantage of modern computational methods to simulate high speed turbulent flows. In doing so, two methodologies were used: large eddy simulations (LES) and direct numerical simulations (DNS). In the work related with LES the objective is to study the behavior of the probability density functions (pdfs) of scalar properties within the subgrid in reacting turbulent flows. The data base obtained by DNS for a detailed study of the pdf characteristics within the subgrid was used. Simulations are performed for flows under various initializations to include the effects of compressibility on mixing and chemical reactions. In the work related with DNS, a two-dimensional temporally developing high speed mixing layer under the influence of a second-order non-equilibrium chemical reaction of the type A + B yields products + heat was considered. Simulations were performed with different magnitudes of the convective Mach numbers and with different chemical kinetic parameters for the purpose of examining the isolated effects of the compressibility and the heat released by the chemical reactions on the structure of the layer. A full compressible code was developed and utilized, so that the coupling between mixing and chemical reactions is captured in a realistic manner.
Application of a partially calibrated National Combustion Code (NCC) for providing guidance in the design of the 3rd generation of the Lean-Direct Injection (LDI) multi-element combustion configuration (LDI-3) is summarized. NCC was used to perform non-reacting and two-phase reacting flow computations on several LDI-3 injector configurations in a single-element and a five-element injector array. All computations were performed with a consistent approach for mesh-generation, turbulence, spray simulations, ignition and chemical kinetics-modeling. Both qualitative and quantitative assessment of the computed flowfield characteristics of the several design options led to selection of an optimal injector LDI- 3 design that met all the requirements including effective area, aerodynamics and fuel-air mixing criteria. Computed LDI-3 emissions (namely, NOx, CO and UHC) will be compared with the prior generation LDI- 2 combustor experimental data at relevant engine cycle conditions.
The present technique measures atmospheric OH concentrations by first titrating OH into isotopically labeled H2(exp 34)SO4 by a series of gas phase reactions in an effectively wall-less flow tube reactor on a time-scale which is short in comparison to the OH lifetime. The gas phase H2SO4 concentration is then determined by reacting it with NO3(-)(center dot)HNO3 ions in a flow tube, and subsequently measuring the resulting (HSO4(-))/(NO3(-)) ion ratio. Since the determination of H2SO4 concentration does not depend on the measurement of an absolute photon flux or chemical concentration but rather only on the relative concentration measurement of two fairly similar mass ions in conjunction with a fixed reaction time and reaction rate coefficient, no in-field calibration technique should be required. The ion-assisted OH measurement technique provides a concentration measurement about once each 10 to 20 seconds and had a detection sensitivity less than or equal to 1 x 10(exp 5) molecules/cu cm (2 standard deviations) with an integration time of 5 minutes during the past field study. Typically, the total background signal level is in the low 10(exp 6) molecule/cu cm range with the scatter in this background (or background noise in 5 minutes) in the mid 10(exp 4) molecules/cu cm range.
High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.
This article presents molecular-level analysis of a reactive, near-continuum, Mach 21 nitrogen flow over a blunt wedge using the direct molecular simulation (DMS) method. The flow conditions lead to internal energy excitation and dissociation in the flow field, resulting in thermal and chemical nonequilibrium in the flow. Thermal nonequilibrium in the vibrational mode is observed to extend to the molecular level, where the vibrational energy distributions at various points in the flow field are observed to be non-Boltzmann. Furthermore, this is the first reactive DMS calculation where the wall is assumed to be isothermal and full momentum accommodation of the particles is enforced, hence incorporating viscous wall effects. Since the DMS method uses a quantum mechanically generated interaction potential as its only modeling input, all thermochemical and transport properties of the flow field can directly be attributed to the ab initio potential energy surface. Using the DMS solution as a benchmark, this article assesses the performance of Navier–Stokes computational fluid dynamics solutions using lower fidelity two-temperature models. Two models are chosen as points of comparison: the well-known Park two-temperature model and the recently developed modified Marrone and Treanor model.
This paper presents a study of the temporal evolution of the combustion flowfield established by the interaction of ram-accelerator-type projectiles with an explosive gas mixture accelerated to hypersonic speeds in an expansion tube. The Navier-Stokes equations for a chemically reacting gas are solved in a fully coupled manner using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order, total variation diminishing (TVD) scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of expansion tube flowfields are then presented for two different configurations. In particular, the development of the shock-induced combustion process is followed. In one case, designed to ensure ignition only in the boundary layer, the lateral extent of the combustion front during the initial transient phase was surprisingly large. The time histories of the calculated thrust and drag forces on the ram accelerator projectile are also presented.
A detailed numerical study was conducted on the dynamics and thermal response of inert, spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport, For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.
A detailed numerical study was conducted on the dynamics and thermal response of inert spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport. For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.
The multi-dimensional numerical model has been developed to analyze the nonlinear combustion instabilities in liquid-fueled engines. The present pressure-based approach can handle the implicit pressure-velocity coupling in a non-iterative way. The additional scalar conservation equations for the chemical species, the energy, and the turbulent transport quantities can be handled by the same predictor-corrector sequences. This method is time-accurate and it can be applicable to the all-speed, transient, multi-phase, and reacting flows. Special emphasis is given to the acoustic/vaporization interaction which may act as the crucial rate-controlling mechanism in the liquid-fueled rocket engines. The subcritical vaporization is modeled to account for the effects of variable thermophysical properties, non-unitary Lewis number in the gas-film, the Stefan flow effect, and the effect of transient liquid heating. The test cases include the one-dimenisonal fast transient non-reacting and reacting flows, and the multi-dimensional combustion instabilities encountered in the liquid-fueled rocket thrust chamber. The present numerical model successfully demonstrated the capability to simulate the fast transient spray-combusting flows in terms of the limiting-cycle amplitude phenomena, correspondence between combustion and acoustics, and the steep-fronted wave and flame propagation. The investigated parameters include the spray initial conditions, air-fuel mixture ratios, and the engine geometry. Stable and unstable operating conditions are found for the liquid-fueled combustors. Under certain conditions, the limiting cycle behavior of the combusting flowfields is obtained. The numerical results indicate that the spray vaporization processes play an important role in releasing thermal energy and driving the combustion instability.
The objective of this effort is to develop an efficient and accurate computational methodology to predict both detailed thermo-fluid environments and global characteristics of the internal ballistics for a hypothetical solid-core nuclear thermal thrust chamber assembly (NTTCA). Several numerical and multi-physics thermo-fluid models, such as real fluid, chemically reacting, turbulence, conjugate heat transfer, porosity, and power generation, were incorporated into an unstructured-grid, pressure-based computational fluid dynamics solver as the underlying computational methodology. The numerical simulations of detailed thermo-fluid environment of a single flow element provide a mechanism to estimate the thermal stress and possible occurrence of the mid-section corrosion of the solid core. In addition, the numerical results of the detailed simulation were employed to fine tune the porosity model mimic the pressure drop and thermal load of the coolant flow through a single flow element. The use of the tuned porosity model enables an efficient simulation of the entire NTTCA system, and evaluating its performance during the design cycle.
Water ices in surface crusts of Europa, Enceladus, Saturn's main rings, and Kuiper Belt Objects can become heavily oxidized from radiolytic chemical alteration of near-surface water ice by space environment irradiation. Oxidant accumulations and gas production are manifested in part through observed H2O2 on Europa. tentatively also on Enceladus, and found elsewhere in gaseous or condensed phases at moons and rings of Jupiter and Saturn. On subsequent chemical contact in sub-surface environments with significant concentrations of primordially abundant reductants such as NH3 and CH4, oxidants of radiolytic origin can react exothermically to power gas-driven cryovolcanism. The gas-piston effect enormously amplifies the mass flow output in the case of gas formation at basal thermal margins of incompressible fluid reservoirs. Surface irradiation, H2O2 production, NH3 oxidation, and resultant heat, gas, and gas-driven mass flow rates are computed in the fluid reservoir case for selected bodies. At Enceladus the oxidant power inputs are comparable to limits on nonthermal kinetic power for the south polar plumes. Total heat output and plume gas abundance may be accounted for at Enceladus if plume activity is cyclic in high and low "Old Faithful" phases, so that oxidants can accumulate during low activity phases. Interior upwelling of primordially abundant NH3 and CH4 hydrates is assumed to resupply the reductant fuels. Much lower irradiation fluxes on Kuiper Belt Objects require correspondingly larger times for accumulation of oxidants to produce comparable resurfacing, but brightness and surface composition of some objects suggest that such activity may be ongoing.
The RPLUS2D code, capable of calculating high speed reacting flows, has been adopted to design a compressible shear layer facility. In order to create reacting shear layers at high convective Mach numbers, hot air streams at supersonic speeds, render by converging-diverging nozzles, must be provided. A finite rate chemistry model is used to simulate the nozzle flows. Results are compared with one-dimensional, chemically equilibrium solutions. Additionally, a two equation turbulence model with compressibility effects has been successfully incorporated with the RPLUS code. The model has been applied to simulate a supersonic shear layer. Preliminary results show favorable comparisons with the experimental data.
Injector geometry, physical mixing, chemical processes, and engine cycle conditions together govern performance, operability and emission characteristics of aviation gas-turbine combustion systems. The present investigation explores swirl-venturi lean direct injection combustor fundamentals, characterizing the influence of key geometric injector parameters on reacting flow physics and emission production trends. In this computational study, a design space exploration was performed using a parameterized swirl-venturi lean direct injector model. From the parametric geometry, 20 three-element lean direct injection combustor sectors were produced and simulated using steady-state, Reynolds-averaged Navier-Stokes reacting computations. Species concentrations were solved directly using a reduced 18-step reaction mechanism for Jet-A. Turbulence closure was obtained using a nonlinear κ-ε model. Results demonstrate sensitivities of the geometric perturbations on axially averaged flow field responses. Output variables include axial velocity, turbulent kinetic energy, static temperature, fuel patternation and minor species mass fractions. Significant trends have been reduced to surrogate model approximations, intended to guide future injector design trade studies and advance aviation gas-turbine combustion research.