Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemically reactive 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Reaction Kernel Structure of a Slot Jet Diffusion Flame in Microgravity

Diffusion flame stabilization in normal earth gravity (1 g) has long been a fundamental research subject in combustion. Local flame-flow phenomena, including heat and species transport and chemical reactions, around the flame base in the vicinity of condensed surfaces control flame stabilization and fire spreading processes. Therefore, gravity plays an important role in the subject topic because buoyancy induces flow in the flame zone, thus increasing the convective (and diffusive) oxygen transport into the flame zone and, in turn, reaction rates. Recent computations show that a peak reactivity (heat-release or oxygen-consumption rate) spot, or reaction kernel, is formed in the flame base by back-diffusion and reactions of radical species in the incoming oxygen-abundant flow at relatively low temperatures (about 1550 K). Quasi-linear correlations were found between the peak heat-release or oxygen-consumption rate and the velocity at the reaction kernel for cases including both jet and flat-plate diffusion flames in airflow. The reaction kernel provides a stationary ignition source to incoming reactants, sustains combustion, and thus stabilizes the trailing diffusion flame. In a quiescent microgravity environment, no buoyancy-induced flow exits and thus purely diffusive transport controls the reaction rates. Flame stabilization mechanisms in such purely diffusion-controlled regime remain largely unstudied. Therefore, it will be a rigorous test for the reaction kernel correlation if it can be extended toward zero velocity conditions in the purely diffusion-controlled regime. The objectives of this study are to reveal the structure of the flame-stabilizing region of a two-dimensional (2D) laminar jet diffusion flame in microgravity and develop a unified diffusion flame stabilization mechanism. This paper reports the recent progress in the computation and experiment performed in microgravity.

Takahashi, F.↗
NASA NTRS2001

Grid Convergence of High Order Methods for Multiscale Complex Unsteady Viscous Compressible Flows

Grid convergence of several high order methods for the computation of rapidly developing complex unsteady viscous compressible flows with a wide range of physical scales is studied. The recently developed adaptive numerical dissipation control high order methods referred to as the ACM and wavelet filter schemes are compared with a fifth-order weighted ENO (WENO) scheme. The two 2-D compressible full Navier-Stokes models considered do not possess known analytical and experimental data. Fine grid solutions from a standard second-order TVD scheme and a MUSCL scheme with limiters are used as reference solutions. The first model is a 2-D viscous analogue of a shock tube problem which involves complex shock/shear/boundary-layer interactions. The second model is a supersonic reactive flow concerning fuel breakup. The fuel mixing involves circular hydrogen bubbles in air interacting with a planar moving shock wave. Both models contain fine scale structures and are stiff in the sense that even though the unsteadiness of the flows are rapidly developing, extreme grid refinement and time step restrictions are needed to resolve all the flow scales as well as the chemical reaction scales.

Sjoegreen, B.↗
DOE OSTI2022

Chemical Kinetics Underlying the Sooting Tendency and Auto-Ignition Characteristics of Linear, Branched, and Cyclic Ether Compounds

Biofuels present opportunities for improving the performance and reducing emissions from internal combustion engines by incorporating oxygenated functional groups to the fuels. Among various oxygenates, ethers have been recognized as promising candidates for an alternative to conventional diesel fuel owing to their higher reactivity and lower sooting tendency. The detailed guidelines for designing ethers, however, have not been fully discussed, even though their combustion characteristics are sensitive to the molecular structure. This study was devoted to exploring the structure-property relationships, particularly focusing on the cetane number and yield sooting index, using five linear, branched, and cyclic ethers: di-amyl-ether, 4-butoxy-heptane, 3,3-dimethyl-oxetane, 2-ethyl-4-methyl-1,3-dioxolane, and 2-isopropyl-4-methyl-1,3-dioxolane. First, we examined the chemical kinetics underlying the sooting tendency of the test fuels. The combustion product distribution was measured from flow reactor experiments at 750-1100 K, F=3, at atmospheric pressure. As a result, it was revealed that the sooting tendency is closely related to the size of hydrocarbon intermediates in the high-temperature regime (>1000 K); that is, larger hydrocarbons lead to more soot precursor formation. The underlying chemistry determining the size of the hydrocarbons from the tested fuels was analyzed using reaction pathway analysis and quantum mechanics calculations, which showed that the branched and cyclic ether structures form abundant C3-C4 compounds. Moreover, the auto-ignition characteristics of the test fuels were studied using the flow reactor at low-temperature (400-700 K) and F=1. We found a clear difference in the combustion-product distribution from high and low reactivity fuels, which was then correlated to the systematic analysis of the key reaction energy barriers with the varying molecular structure.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2020

Predicting Water Cycle Characteristics from Percolation Theory and Observational Data

The fate of water and water-soluble toxic wastes in the subsurface is of high importance for many scientific and practical applications. Although solute transport is proportional to water flow rates, theoretical and experimental studies show that heavy-tailed (power-law) solute transport distribution can cause chemical transport retardation, prolonging clean-up time-scales greatly. However, no consensus exists as to the physical basis of such transport laws. In percolation theory, the scaling behavior of such transport rarely relates to specific medium characteristics, but strongly to the dimensionality of the connectivity of the flow paths (for example, two- or three-dimensional, as in fractured-porous media or heterogeneous sediments), as well as to the saturation characteristics (i.e., wetting, drying, and entrapped air). In accordance with the proposed relevance of percolation models of solute transport to environmental clean-up, these predictions also prove relevant to transport-limited chemical weathering and soil formation, where the heavy-tailed distributions slow chemical weathering over time. The predictions of percolation theory have been tested in laboratory and field experiments on reactive solute transport, chemical weathering, and soil formation and found accurate. Recently, this theoretical framework has also been applied to the water partitioning at the Earth’s surface between evapotranspiration, ET, and run-off, Q, known as the water balance. A well-known phenomenological model by Budyko addressed the relationship between the ratio of the actual evapotranspiration (ET) and precipitation, ET/P, versus the aridity index, ET 0 /P, with P being the precipitation and ET 0 being the potential evapotranspiration. Existing work was able to predict the global fractions of P represented by Q and ET through an optimization of plant productivity, in which downward water fluxes affect soil depth, and upward fluxes plant growth. In the present work, based likewise on the concepts of percolation theory, we extend Budyko’s model, and address the partitioning of run-off Q into its surface and subsurface components, as well as the contribution of interception to ET. Using various published data sources on the magnitudes of interception and information regarding the partitioning of Q, we address the variability in ET resulting from these processes. The global success of this prediction demonstrated here provides additional support for the universal applicability of percolation theory for solute transport as well as guidance in predicting the component of subsurface run-off, important for predicting natural flow rates through contaminated aquifers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2021

Core-flood Effluent and Shale Surface Chemistries in Predicting Interaction between Shale, Brine, and Reactive Fluid

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1995

Single-drop reactive extraction/extractive reaction with forced convective diffusion and interphase mass transfer

An algorithm has been developed for the forced convective diffusion-reaction problem for convection inside and outside a droplet by a recirculating flow field hydrodynamically coupled at the droplet interface with an external flow field that at infinity becomes a uniform streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet or reactions can take place in both phases. The algorithm has been implemented and results are shown here for the case of no reaction and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Kleinman, Leonid S.↗
NASA NTRS1999

The Mars Environmental Compatibility Assessment (MECA) Wet Chemistry Experiment on the Mars 2001 Lander

The Mars Environmental Compatibility Assessment (MECA) is an instrument suite that will fly on the Mars Surveyor 2001 Lander Spacecraft. MECA is sponsored by the Human Exploration and Development of Space (HEDS) program and will evaluate potential hazards that the dust and soil of Mars might present to astronauts and their equipment on a future human mission to Mars. Four elements constitute the integrated MECA payload: a microscopy station, patch plates, an electrometer, and the wet chemistry laboratory (WCL). The WCL consists of four identical cells, each of which will evaluate a sample of Martian soil in water to determine conductivity, pH, redox potential, dissolved C02 and 02 levels, and concentrations of many soluble ions including sodium, potassium, magnesium, calcium and the halides. In addition, cyclic voltammetry will be used to evaluate reversible and irreversible oxidants present in the water/soil solution. Anodic stripping voltammetry will be used to measure concentrations of trace metals including lead, copper, and cadmium at ppb levels. Voltammetry is a general electrochemical technique that involves controlling the potential of an electrode while simultaneously measuring the current flowing at that electrode. The WCL experiments will provide information on the corrosivity and reactivity of the Martian soil, as well as on soluble components of the soil which might be toxic to human explorers. They will also guide HEDS scientists in the development of high fidelity Martian soil simulants. In the process of acquiring information relevant to HEDS, the WCL will assess the chemical composition and properties of the salts present in the Martian soil.

Grannan, S. M.↗
DOE OSTI2021

Microscale Mechanical‐Chemical Modeling of Granular Salt: Insights for Creep

Abstract Microscale numerical modeling is potentially an effective approach for understanding salt creep, but it is quite challenging because natural salt grain boundaries can have complex and evolving geometry, and because contacts between these deformable salt mineral grains are dynamic as a result of compaction, chemical reaction, fluid flow, and heat transfer. In this study, we have overcome these challenges and developed a new microscale mechanical‐chemical (MC) model to analyze creep of salt at the microscale, accounting for coupled deformation, dynamic contacts, and chemical reaction in granular systems with realistic geometric representations. The MC model was realized by linking a new microscale mechanical code based on the numerical manifold method (NMM) to a reactive transport code named Crunch. We simulated the processes of reorganization of the salt grains, microfracturing, and pressure solution that contribute to creep at larger scales. Based on this first quantitative microscale model, we found that sharp corners of mineral grains can dominate the contact dynamics, microfracturing, and pressure solution when salt is compacted, thus governing the structural changes and porosity loss of the system. We found that pressure solution, which preferentially dissolves sharp corners and edges, can lead to relatively high porosity loss in the system, thus playing an important role in the creep of salt. Our analysis shows that the dynamic changes of salt granular systems involving grain relocation and pressure solution can occur repeatedly and continuously, thus contributing to longer‐term creep of salt at larger scales.

58 GEOSCIENCES↗
DOE OSTI2023

Protection against Chemical Warfare Agents and Biological Threats Using Metal–Organic Frameworks as Active Layers

The SARS-CoV-2 pandemic outbreak and the unfortunate misuse of toxic chemical warfare agents (CWAs) highlight the importance of developing functional materials to protect against these chemical and pathogen threats. Metal–organic frameworks (MOFs), which comprise a tunable class of crystalline porous materials built from inorganic nodes and organic linkers, have emerged as a class of heterogeneous catalysts capable of rapid detoxification of multiple classes of these harmful chemical or biological hazards. In particular, zirconium-based MOFs (Zr-MOFs) feature Lewis acidic nodes that serve as active sites for a wide range of catalytic reactions, including the hydrolysis of organophosphorus nerve agents within seconds in basic aqueous solutions. In addition, postsynthetic modification of Zr-MOFs enables the release of active species capable of reacting with and deactivating harmful pathogens. Despite this impressive performance, utilizing Zr- MOFs in powder form is not practical for application in masks or protective uniforms. To address this challenge, our team sought to develop MOF/fiber composite systems that could be adapted for use under realistic operating conditions to protect civilians, military personnel, and first responders from harmful pathogens and chemical warfare agents. Over the last several years, our group has designed and fabricated reactive and biocidal MOF/fiber composites that effectively capture and deactivate these toxic species. In this Account, we describe the evolution of these porous and reactive MOF/fiber composites and focus on key design challenges and considerations. First, we devised a scalable method for the integration of Zr-MOFs onto textile substrates using aqueous precursor solutions and without using pretreated textiles, highlighting the potential scalability of this method. Moving beyond standard textiles, we also developed a microbial synthesis strategy to prepare hierarchically porous MOF/bacterial cellulose nanofiber composite sponges that can both capture and detoxify nerve agents when exposed to contaminated gas flows. The mass loading of the MOF in the nanofibrous composite sponge is up to 90%, affording higher work capacities compared to those of textile-fiber-based composites with relatively lower MOF loadings. Next, we demonstrated that heterogeneous polymeric bases are suitable replacements for volatile liquid bases typically used in solution-phase reactions, and we showed that these composite systems are capable of effectively hydrolyzing nerve agents in the solid state by using only water that is present as humidity. Moreover, incorporating a reactive dye precursor into the composite affords a dual function sensing and detoxifying material that changes color from white to orange upon reaction with the byproduct following nerve agent hydrolysis, demonstrating the versatility of this platform for use in decontamination applications. We then created chlorine-loaded MOF/fiber composites that act as biocidal and reactive textiles that are capable of not only detoxifying sulfur-mustard-based chemical warfare agents and simulants but also deactivating both bacteria and the SARS-CoV-2 virus within minutes of exposure. Lastly, we synthesized a mixed-metal Ti/Zr-MOF coating on cotton fibers to afford a photoactive biocidal cloth that shows fast and broad-spectrum biocidal performance against viruses and Gram-positive and Gram-negative bacteria under visible light irradiation. Given the tunable, multifunctional nature of these MOF/fiber composites, we believe that this Account will offer new insights for the rational design and preparation of functional MOF/fiber composites and pave the way toward the development of next-generation reactive and protective textiles.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Silicon Microleaks for Inlets of Mass Spectrometers

Microleaks for inlets of mass spectrometers used to analyze atmospheric gases can be fabricated in silicon wafers by means of photolithography, etching, and other techniques that are commonly used in the manufacture of integrated circuits and microelectromechanical systems. The microleaks serve to limit the flows of the gases into the mass-spectrometer vacuums to specified very small flow rates consistent with the capacities of the spectrometer vacuum pumps. There is a need to be able to precisely tailor the dimensions of each microleak so as to tailor its conductance to a precise low value. (As used here, "conductance" signifies the ratio between the rate of flow in the leak and the pressure drop from the upstream to the downstream end of the leak.) To date, microleaks have been made, variously, of crimped metal tubes, pulled glass tubes, or frits. Crimped-metal and pulled-glass-tube microleaks cannot readily be fabricated repeatably to precise dimensions and are susceptible to clogging with droplets or particles. Frits tend to be differentially chemically reactive with various gas constituents and, hence, to distort the gas mixtures to be analyzed. The present approach involving microfabrication in silicon largely overcomes the disadvantages of the prior approaches.

Harpold, Dan↗
NASA NTRS1995

Single-drop reactive extraction/extractive reaction with forced convective diffusion and interphase mass transfer

An algorithm has been developed for time-dependent forced convective diffusion-reaction having convection by a recirculating flow field within the drop that is hydrodynamically coupled at the interface with a convective external flow field that at infinity becomes a uniform free-streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet, or reactions can take place in both phases. The algorithm has been implemented, and for comparison results are shown here for the case of no reaction in either phase and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Kleinman, Leonid S.↗
DOE OSTI2025

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗
NASA NTRS1997

Science Support for Space-Based Droplet Combustion: Drop Tower Experiments and Detailed Numerical Modeling

This program supports the engineering design, data analysis, and data interpretation requirements for the study of initially single component, spherically symmetric, isolated droplet combustion studies. Experimental emphasis is on the study of simple alcohols (methanol, ethanol) and alkanes (n-heptane, n-decane) as fuels with time dependent measurements of drop size, flame-stand-off, liquid-phase composition, and finally, extinction. Experiments have included bench-scale studies at Princeton, studies in the 2.2 and 5.18 drop towers at NASA-LeRC, and both the Fiber Supported Droplet Combustion (FSDC-1, FSDC-2) and the free Droplet Combustion Experiment (DCE) studies aboard the shuttle. Test matrix and data interpretation are performed through spherically-symmetric, time-dependent numerical computations which embody detailed sub-models for physical and chemical processes. The computed burning rate, flame stand-off, and extinction diameter are compared with the respective measurements for each individual experiment. In particular, the data from FSDC-1 and subsequent space-based experiments provide the opportunity to compare all three types of data simultaneously with the computed parameters. Recent numerical efforts are extending the computational tools to consider time dependent, axisymmetric 2-dimensional reactive flow situations.

Marchese, Anthony J.↗
DOE OSTI2022

Non-Idealities in Lab-Scale Kinetic Testing: A Theoretical Study of a Modular Temkin Reactor

The Temkin reactor can be applied for industrial relevant catalyst testing with unmodified catalyst particles. It was assumed in the literature that this reactor behaves as a cascade of continuously stirred tank reactors (CSTR). However, this assumption was based only on outlet gas composition or inert residence time distribution measurements. The present work theoretically investigates the catalytic CO2 methanation as a test case on different catalyst geometries, a sphere, and a ring, inside a single Temkin reaction chamber under isothermal conditions. Axial gas-phase species profiles from detailed computational fluid dynamics (CFD) are compared with a CSTR and 1D plug-flow reactor (PFR) model using a sophisticated microkinetic model. In addition, a 1D chemical reactor network (CRN) model was developed, and model parameters were adjusted based on the CFD simulations. Whereas the ideal reactor models overpredict the axial product concentrations, the CRN model results agree well with the CFD simulations, especially under low to medium flow rates. This study shows that complex flow patterns greatly influence species fields inside the Temkin reactor. Although residence time measurements suggest CSTR-like behavior, the reactive flow cannot be described by either a CSTR or PFR model but with the developed CRN model.

Wehinger, Gregor D. (ORCID:0000000217743391)↗
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
NASA NTRS1981

Review of literature surface tension data for molten silicon

Measurements of the surface tension of molten silicon are reported. For marangoni flow, the important parameter is the variation of surface tension with temperature, not the absolute value of the surface tension. It is not possible to calculate temperature coefficients using surface tension measurements from different experiments because the systematic errors are usually larger than the changes in surface tension because of temperature variations. The lack of good surface tension data for liquid silicon is probably due to its extreme chemical reactivity. A material which resists attack by molten silicon is not found. It is suggested that all of the sessile drip surface tension measurements are probably for silicon which is contaminated by the substrate materials.

Hardy, S.↗
NASA NTRS1992

A subgrid mixing model for LES of nonpremixed turbulent combustion

A new sub-grid mixing model for use in large eddy simulations of turbulent combustion is presented and applied to a hydrogen-air diffusion flame. The sub-grid model is based on Kerstein's Linear Eddy Model (1988). The model is first used to predict the mixing of a conserved scalar in a turbulent shear flow. The model correctly predicts the behavior of the pdf of the scalar field and displays a nonmarching peak at the preferred mixture fraction as the shear layer is traversed. It is then illustrated how a reduced chemical mechanism can be implemented within the linear eddy subgrid model formulation. The model is used to predict NO formation in hydrogen-air diffusion flame using a reduced chemical mechanism involving nine reactive scalars.

Mcmurtry, P. A.↗