Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemically reacting flows”

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 289 records · Page 16

Direct simulations of chemically reacting turbulent mixing layers, part 2

The results of direct numerical simulations of chemically reacting turbulent mixing layers are presented. This is an extension of earlier work to a more detailed study of previous three dimensional simulations of cold reacting flows plus the development, validation, and use of codes to simulate chemically reacting shear layers with heat release. Additional analysis of earlier simulations showed good agreement with self similarity theory and laboratory data. Simulations with a two dimensional code including the effects of heat release showed that the rate of chemical product formation, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release. Subsequent three dimensional simulations showed similar behavior, in agreement with laboratory observations. Baroclinic torques and thermal expansion in the mixing layer were found to produce changes in the flame vortex structure that act to diffuse the pairing vortices, resulting in a net reduction in vorticity. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layers were shown to result from vorticity generation by baroclinic torques.

Metcalfe, Ralph W.↗

Highly-scalable GPU-accelerated compressible reacting flow solver for modeling high-speed flows

Emerging supercomputing systems utilize a combination of central processing units (CPUs) and graphics processing units (GPUs) in an effort to reach exascale capabilities while minimizing the energy footprint of operating such systems. Such heterogeneous machines introduce new challenges for fluids solvers because the hardware architecture and operation of a GPU are fundamentally different from conventional CPUs. In this work, a general approach for efficient implementation of finite-volume based reacting flow solvers on such heterogeneous systems is presented. Three main challenges, namely, data access pattern, thread divergence, and thread safety, are addressed. Since compressible reacting flows require special methods to deal with chemical reactions, hyperbolic and nonlinear convection terms, and the presence of turbulence, specific algorithms that ensure GPU-based efficiency are developed. The approach is demonstrated on the widely available OpenFOAM open source software by modifying core algorithms for GPU accessibility. The scalability of the resulting solver, is demonstrated using practical test cases, including flow through a scramjet engine and the dynamics of a rotating detonation engine. Here, the solver provides near-ideal scaleup on a large number of GPUs (>3000), and extremely efficient use of the GPUs, with throughput nearly a constant even when processing a large number of control volumes.

42 ENGINEERING↗

AdapChem

AdapChem software enables high efficiency, low computational cost, and enhanced accuracy on computational fluid dynamics (CFD) numerical simulations used for combustion studies. The software dynamically allocates smaller, reduced chemical models instead of the larger, full chemistry models to evolve the calculation while ensuring the same accuracy to be obtained for steady-state CFD reacting flow simulations. The software enables detailed chemical kinetic modeling in combustion CFD simulations. AdapChem adapts the reaction mechanism used in the CFD to the local reaction conditions. Instead of a single, comprehensive reaction mechanism throughout the computation, a dynamic distribution of smaller, reduced models is used to capture accurately the chemical kinetics at a fraction of the cost of the traditional single-mechanism approach.

Oluwole, Oluwayemisi O.↗

Earth's Early Biosphere and the Biogeochemical Carbon Cycle

Our biosphere has altered the global environment principally by influencing the chemistry of those elements most important for life, e g., C, N, S, O, P and transition metals (e.g., Fe and Mn). The coupling of oxygenic photosynthesis with the burial in sediments of photosynthetic organic matter, and with the escape of H2 to space, has increased the state of oxidation of the Oceans and atmosphere. It has also created highly reduced conditions within sedimentary rocks that have also extensively affected the geochemistry of several elements. The decline of volcanism during Earth's history reduced the flow of reduced chemical species that reacted with photosynthetically produced O2. The long-term net accumulation of photosynthetic O2 via biogeochemical processes has profoundly influenced our atmosphere and biosphere, as evidenced by the O2 levels required for algae, multicellular life and certain modem aerobic bacteria to exist. When our biosphere developed photosynthesis, it tapped into an energy resource that was much larger than the energy available from oxidation-reduction reactions associated with weathering and hydrothermal activity. Today, hydrothermal sources deliver globally (0.13-1.1)x10(exp l2) mol yr(sup -1) of reduced S, Fe(2+), Mn(2+), H2 and CH4; this is estimated to sustain at most about (0.2-2)xl0(exp 12)mol C yr(sup -1) of organic carbon production by chemautotrophic microorganisms. In contrast, global photosynthetic productivity is estimated to be 9000x10(exp 12) mol C yr(sup -1). Thus, even though global thermal fluxes were greater in the distant geologic past than today, the onset of oxygenic photosynthesis probably increased global organic productivity by some two or more orders of magnitude. This enormous productivity materialized principally because oxygenic photosynthesizers unleashed a virtually unlimited supply of reduced H that forever freed life from its sole dependence upon abiotic sources of reducing power such as hydrothermal emanations and weathering. Communities sustained by oxygenic photosynthesis apparently thrived wherever supplies of sunlight, moisture and nutrients were sufficient. Prior to the development of oxygenic photosynthesis, the net global effect of the ancient global biosphere was to facilitate chemical equilibrium between reduced species from thermal activity and weathering and more oxidized constituents in the surface environment. But even this ancient biosphere might have been globally pervasive. The global geothermal heat flow was substantially higher during Earth's first billion years, and thus reduced chemical species might have persisted in sunlit aquatic environments. Perhaps the substantial decline in thermal activity between 4000 and 3000 Ma created a driver for oxygenic photosynthesis to develop.

DesMarais, David↗

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗

Atomic Oxygen Treatment Technique for Removal of Smoke Damage from Paintings

Soot deposits that can accumulate on surfaces of a painting during a fire can be difficult to clean from some types of paintings without damaging the underlying paint layers. A non-contact technique has been developed which can remove the soot by allowing a gas containing atomic oxygen to flow over the surface and chemically react with the soot to form carbon monoxide and carbon dioxide. The reaction is limited to the surface, so the underlying paint is not touched. The process can be controlled so that the cleaning can be stopped once the paint surface is reached. This paper describes the smoke exposure and cleaning of untreated canvas, acrylic gesso, and sections of an oil painting using this technique. The samples were characterized by optical microscopy and reflectance spectroscopy.

Rutledge, S. K.↗

Simultaneous Laser Raman-rayleigh-lif Measurements and Numerical Modeling Results of a Lifted Turbulent H2/N2 Jet Flame in a Vitiated Coflow

An experiment and numerical investigation is presented of a lifted turbulent H2/N2 jet flame in a coflow of hot, vitiated gases. The vitiated coflow burner emulates the coupling of turbulent mixing and chemical kinetics exemplary of the reacting flow in the recirculation region of advanced combustors. It also simplifies numerical investigation of this coupled problem by removing the complexity of recirculating flow. Scalar measurements are reported for a lifted turbulent jet flame of H2/N2 (Re = 23,600, H/d = 10) in a coflow of hot combustion products from a lean H2/Air flame ((empty set) = 0.25, T = 1,045 K). The combination of Rayleigh scattering, Raman scattering, and laser-induced fluorescence is used to obtain simultaneous measurements of temperature and concentrations of the major species, OH, and NO. The data attest to the success of the experimental design in providing a uniform vitiated coflow throughout the entire test region. Two combustion models (PDF: joint scalar Probability Density Function and EDC: Eddy Dissipation Concept) are used in conjunction with various turbulence models to predict the lift-off height (H(sub PDF)/d = 7,H(sub EDC)/d = 8.5). Kalghatgi's classic phenomenological theory, which is based on scaling arguments, yields a reasonably accurate prediction (H(sub K)/d = 11.4) of the lift-off height for the present flame. The vitiated coflow admits the possibility of auto-ignition of mixed fluid, and the success of the present parabolic implementation of the PDF model in predicting a stable lifted flame is attributable to such ignition. The measurements indicate a thickened turbulent reaction zone at the flame base. Experimental results and numerical investigations support the plausibility of turbulent premixed flame propagation by small scale (on the order of the flame thickness) recirculation and mixing of hot products into reactants and subsequent rapid ignition of the mixture.

Cabra, R.↗

A Computational Study of Plume Modeling for Space Launch System Abort Scenarios

A series of viscous CFD simulations depicting Space Launch System (SLS) post Mode-1abort scenarios are conducted and analyzed. The primary purpose of these simulations is to study the effect of launch abort vehicle attitude control motor (ACM) and abort motor(AM) plume modeling fidelity on the drag of the aborted core stage. The simulations are first conducted with a fully-coupled, multi-species and chemically reacting model for the freestream gases and the solid rocket abort motor combustion products. Following this, equivalent species are created for the ACM and AM exhaust gases. These equivalent species aim to provide near the same results for core stage drag as the chemically reacting model at a significantly lower computational cost. Finally, the exhaust gases are modeled as calorically perfect air. It was concluded that under most scenarios, both the equivalent species and the perfect air models predict the drag of the aborted core stage sufficiently close to that predicted by the chemically reacting simulations. In these cases, the flow and geometry conditions did not induce significant changes in the composition of the ‘true’ exhaust gas via afterburning of the combustion products. For cases where this was not the case, the properties of the essentially frozen equivalent species and perfect air models deviated from the true mixture, yielding unsatisfactory results. Finally, the accuracy improvements by using the chemically reacting model are weighed against the significant increase in associated computational costs. It was found that the simplified exhaust gas models are at least 4 times less expensive than the chemically reacting model. This may justify the judicious use of these simple models in future work.

CFD↗

A Computational Study of Plume Modeling for Space Launch System Abort Scenarios

A series of viscous CFD simulations depicting Space Launch System (SLS) post Mode-1abort scenarios are conducted and analyzed. The primary purpose of these simulations is to study the effect of launch abort vehicle attitude control motor (ACM) and abort motor(AM) plume modeling fidelity on the drag of the aborted core stage. The simulations are first conducted with a fully-coupled, multi-species and chemically reacting model for the freestream gases and the solid rocket abort motor combustion products. Following this, equivalent species are created for the ACM and AM exhaust gases. These equivalent species aim to provide near the same results for core stage drag as the chemically reacting model at a significantly lower computational cost. Finally, the exhaust gases are modeled as calorically perfect air. It was concluded that under most scenarios, both the equivalent species and the perfect air models predict the drag of the aborted core stage sufficiently close to that predicted by the chemically reacting simulations. In these cases, the flow and geometry conditions did not induce significant changes in the composition of the ‘true’ exhaust gas via afterburning of the combustion products. For cases where this was not the case, the properties of the essentially frozen equivalent species and perfect air models deviated from the true mixture, yielding unsatisfactory results. Finally, the accuracy improvements by using the chemically reacting model are weighed against the significant increase in associated computational costs. It was found that the simplified exhaust gas models are at least 4 times less expensive than the chemically reacting model. This may justify the judicious use of these simple models in future work.

CFD↗

Symbolic construction of the chemical Jacobian of quasi-steady state (QSS) chemistries for Exascale computing platforms

The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of chemical mechanisms for implementation in computational reacting flow solvers. However, for many applications, the resulting model still requires implicit methods for efficient time integration. Here, in this paper, we outline an approach to formulating the QSSA reduction that is coupled with a strategy to generate C++ source code to evaluate the net species production rates, and the chemical Jacobian. The code-generation component employs a symbolic approach enabling a simple and effective strategy to analytically compute the chemical Jacobian. For computational tractability, the symbolic approach needs to be paired with common subexpression elimination which can negatively affect memory usage. Several solutions are outlined and successfully tested on a 3D multipulse ignition problem, thus allowing portable application across chemical model sizes and GPU capabilities. The implementation of the proposed method is available at https://github.com/AMReX-Combustion/PelePhysics under an open-source license.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Mass Transfer on Aeroheating in Hypersonic Chemically Reacting Boundary Layers

Heat flux characterization of high-enthalpy boundary layer flows is key to optimize the performance and design of Thermal Protection System of next generation aerospace vehicles [1]. At atmospheric entry hypersonic speeds, ablation as well as surface catalycity impact boundary layer aeroheating. Out-gassing occurring from an ablative surface in planetary entry environment introduces a rich set of problems in thermodynamic, fluid dynamic, and material pyrolysis. Ablation leads to out-gassing and surface roughness, both of which are known to affect surface heating in hypersonic chemically reacting boundary layers via three main routes: gas blowing into the boundary layer from the wall, changing the surface heat transfer due to wall-flow chemical reactions, and modifying surface roughness via ablative processes.

Karimi, Mona↗

LES, DNS, and RANS for the Analysis of High-Speed Turbulent Reacting Flows

A filtered density function (FDF) method suitable for chemically reactive flows is developed in the context of large eddy simulation. The advantage of the FDF methodology is its inherent ability to resolve subgrid scales (SGS) scalar correlations that otherwise have to be modeled. Because of the lack of robust models to accurately predict these correlations in turbulent reactive flows, simulations involving turbulent combustion are often met with a degree of skepticism. The FDF methodology avoids the closure problem associated with these terms and treats the reaction in an exact manner. The scalar FDF approach is particularly attractive since it can be coupled with existing hydrodynamic computational fluid dynamics (CFD) codes.

Colucci, P. J.↗

Chemical timescale effects on detonation convergence

Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next- generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver-which enables resolutions up to ($\mathcal{O}$ (1000) cells-per-induction length-the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen- air mechanism, where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. Furthermore, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.

Adaptive mesh refinement↗

Stochastic modeling of turbulent reacting flows

Direct numerical simulations of a single-step irreversible chemical reaction with non-premixed reactants in forced isotropic turbulence at R(sub lambda) = 63, Da = 4.0, and Sc = 0.7 were made using 128 Fourier modes to obtain joint probability density functions (pdfs) and other statistical information to parameterize and test a Fokker-Planck turbulent mixing model. Preliminary results indicate that the modeled gradient stretching term for an inert scalar is independent of the initial conditions of the scalar field. The conditional pdf of scalar gradient magnitudes is found to be a function of the scalar until the reaction is largely completed. Alignment of concentration gradients with local strain rate and other features of the flow were also investigated.

Fox, R. O.↗

Oxidation Behavior of Carbon Fiber-Reinforced Composites

OXIMAP is a numerical (FEA-based) solution tool capable of calculating the carbon fiber and fiber coating oxidation patterns within any arbitrarily shaped carbon silicon carbide composite structure as a function of time, temperature, and the environmental oxygen partial pressure. The mathematical formulation is derived from the mechanics of the flow of ideal gases through a chemically reacting, porous solid. The result of the formulation is a set of two coupled, non-linear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations are solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of the spatial location and time. The local rate of carbon oxidation is determined at each time step using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The non-linear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual finite element method, allowing for the solution of the differential equations numerically.

Sullivan, Roy M.↗

Calculation of Propulsive Nozzle Flowfields in Multidiffusing Chemically Reacting Environments

An advanced engineering model has been developed to aid in the analysis and design of hydrogen/oxygen chemical rocket engines. The complete multispecies, chemically reacting and multidiffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and the Dufour energy transfer terms. In addition to the spectrum of multispecies aspects developed, the model developed in this study is also conservative in axisymmetric flow for both inviscid and viscous flow environments and the boundary conditions employ a viscous, chemically reacting, reference plane characteristics method. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film cooled nozzle, and a transpiration cooled plug and spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 and the 25 lbf film cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements when all of the chemical reaction and diffusion terms are considered. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. The Soret thermal diffusion term is demonstrated to have a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle in both the laminar flow 1030:1 nozzle and the turbulent flow plug and spool nozzle analysis cases performed. Further, the Soret term was shown to represent an important fraction of the diffusion fluxes occurring in a transpiration cooled rocket engine.

Kacynski, Kenneth John↗