Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive Flow”

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 973 records · Page 54

Dynamic Surface Incorporation of Pb 2+ Ions at the Actively Dissolving Calcite (104) Surface

The reaction of dissolved Pb 2+ with calcite surfaces at near-equilibrium conditions involves adsorption of Pb 2+ and precipitation of secondary heteroepitaxial Pb-carbonate minerals. A more complex behavior is observed under far-from-equilibrium conditions, including strong inhibition of calcite dissolution, development of microtopography, and near-surface incorporation of multiple monolayers (ML) of Pb 2+ without precipitation of secondary phases [where 1 ML ≡ 1 Ca/20.2 Å 2 , the crystallographic site density of the calcite (104) lattice plane]. However, the mechanistic controls governing far-from-equilibrium reactivity are not well understood. Here, in this study, we observe the interfacial incorporation of dissolved Pb 2+ during the dissolution of calcite (104) surfaces at pH ~3.7 in a flow-through reaction cell, revealing the formation of a ~1 nm thick Pb-rich calcite layer with a total Pb coverage of ~1.4 ML. These observations of the sorbed Pb distribution used resonant anomalous X-ray reflectivity, X-ray fluorescence, and nanoinfrared atomic force microscopy. We propose that this altered surface layer represents a novel sorption mode that is stabilized by conditions of sustained disequilibrium. This behavior may significantly impact the transport of dissolved metals during disequilibrium processes occurring in acid mine drainage and subsurface CO 2 injection and, if appropriately accounted for, could improve the predictive capability of geochemical reactive-transport models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New Instrument for Time-Resolved OH and HO 2 Quantification in High-Pressure Laboratory Kinetics Studies

Here, we have constructed a new time-resolved high-pressure fluorescence assay by gas expansion (HP-FAGE) apparatus, optimized for the detection of OH and HO 2 radicals in complex gas-phase reactions. The new instrument fills a gap in the existing experimental toolkit for chemical kinetics by enabling the quantification of two key reactive species with microsecond time resolution from high-pressure sources, which was previously not attainable. The HP-FAGE is interfaced with a flow reactor, designed for pressures up to 100 bar and temperatures up to 1000 K, in which reactions are initiated by laser photolysis of radical precursors at repetition rates of 1–10 Hz. The HP-FAGE samples gas out of the reactor into a miniature FAGE chamber, where OH is detected by resonant laser-induced fluorescence using a time-delayed probe laser pulse. HO 2 is converted to OH via reaction with NO and then detected by OH fluorescence. The novel FAGE design places the probe region very close to the gas expansion, minimizing the transport time of sampled molecules and resulting in time resolution better than 20 μs for both OH and HO 2 . We calibrate the sensitivity of HP-FAGE, validate its performance with measurements of well-known reaction kinetics (OH + CH 4 , OH + OH, OH + HO 2 , and HO 2 + HO 2 ), and discuss prospects for its future use.

calibration↗

[18F]Tosyl fluoride as a versatile [18F]fluoride source for the preparation of 18F-labeled radiopharmaceuticals

Abstract Positron emission tomography (PET) is an in vivo imaging technology that utilizes positron-emitting radioisotope-labeled compounds as PET radiotracers that are commonly used in clinic and in various research areas, including oncology, cardiology, and neurology. Fluorine-18 is the most widely used PET-radionuclide and commonly produced by proton bombardment of 18 O-enriched water in a cyclotron. The [ 18 F]fluoride thus obtained generally requires processing by azeotropic drying in order to completely remove H 2 O before it can be used for nucleophilic radiofluorination. In general, the drying step is important in facilitating the radiofluorination reactions and the preparation of 18 F-labeled PET radiotracers. In this communication, we have demonstrated the feasibility of using [ 18 F]tosyl fluoride ([ 18 F]TsF) as a versatile [ 18 F]fluoride source for radiofluorination to bypass the azeotropic drying step, and we have developed a continuous flow solid-phase radiosynthesis strategy to generate [ 18 F]TsF in a form that is excellent for radiofluorination. [ 18 F]TsF shows high reactivity in radiofluorination and provides the features suitable for preparing PET radiotracers on a small scale and exploring novel radiolabeling technologies. Thus, using [ 18 F]TsF as a [ 18 F]fluoride source is a promising strategy that facilitates radiofluorination and provides a convenient and efficient solution for the preparation of 18 F-labeled radiopharmaceuticals that is well matched to the emerging trends in PET imaging technologies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Primary-Side Power Flow Control of Wireless Power Transfer for Electric Vehicle Charging

Various noncontacting methods of plug-in electric vehicle charging are either under development or now deployed as aftermarket options in the light-duty automotive market. Wireless power transfer (WPT) is now the accepted term for wireless charging and is used synonymously for inductive power transfer and magnetic resonance coupling. WPT technology is in its infancy; standardization is lacking, especially on interoperability, center frequency selection, magnetic fringe field suppression, and the methods employed for power flow regulation. This paper proposes a new analysis concept for power flow in WPT in which the primary provides frequency selection and the tuned secondary, with its resemblance to a power transmission network having a reactive power voltage control, is analyzed as a transmission network. Analysis is supported with experimental data taken from Oak Ridge National Laboratory s WPT apparatus. Lastly, this paper also provides an experimental evidence for frequency selection, fringe field assessment, and the need for low-latency communications in the feedback path.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Assessment of CO 2 Storage Resources in Depleted Oil and Gas Fields in the Ship Shoal Area, Gulf of Mexico

The Gulf of Mexico is one of the most important regions in the United States for energy resources and infrastructure. Gulf of Mexico federal offshore oil and gas production accounts for 17% of total U.S. crude oil production and 5% of total U.S. gas production (EIA Gulf of Mexico Fact Sheet). This region presents an excellent combination of high need and significant opportunity for large scale geologic storage of CO 2 . The Ship Shoal Area is located about 20 miles offshore Louisiana within the Gulf Coast federal waters. Miocene and Pliocene sediments in the Ship Shoal Area are proven to provide excellent and secure traps for oil and gas. The Ship Shoal Area contains a large number of depleted oil and gas fields either currently abandoned, or planned for abandonment by 2025, which may provide very significant potential CO 2 storage capacity. GeoMechanics Technologies conducted a comprehensive research project to better characterize Neogene sediments in the Ship Shoal Area for high volume CO 2 storage. The data generated from literature, well data, well log, and formation evaluation were input into Rockwork 16 geologic software to create a detailed geologic model for the Ship Shoal study area. The geologic grids were then fed into TOUGH2, the gas migration model software and FLAC 3D, the geomechanical model software. The fluid flow modeling for both the Miocene and Pliocene simulation from both Block 107 and Block 84 indicated that Pliocene and Miocene are a good reservoir for the CO 2 sequestration. Thirty years of CO 2 injection (at a rate of 1 million metric tons CO 2 per year) and 30 years of observation simulations were run with the fluid flow modeling. The 3D geomechanical modeling for Ship Shoal Block 107 and Block 84 were developed to evaluate induced surface deformation and potential fault reactivation after 30 years CO 2 injection at base of Pliocene and at upper Miocene injection locations. Evaluation results indicate low to no risks for fault slips or fault reactivation. Twelve well bores and 76 well bores for Ship Shoal Blocks 84 and 107 respectively were reviewed for their cement history. Most wells (57 out of 76 and 7 out of 12 wells in Ship Shoal Block 107 and 84 respectively) have good integrity. Nineteen and 5 wells (in Ship Shoal Block 107 and 84 respectively) with no top plug, incomplete cement or Plug and Abandonment information are given yellow cautionary indicators. These cautionary wells may provide leakage paths of CO 2 through the well bores to the USDWs. Using our Quantitative Risk & Decision Analysis Tool (QRDAT) for caprock integrity evaluation, we compared Ship Shoal’s risk to that of In Salah, Sleipner, Kevin Dome, Loudon, Illinois Industrial CCS and Wilmington Graben. We found the risk at the Ship Shoal Blocks 84 and 107 fields are similar to the known CO 2 active sequestration sites, but lower than the Wilmington Graben turbidities offshore California studied site. The risk of natural seismicity in the Gulf of Mexico is relatively low. Geomechanics Technologies has documented the top 25 CO 2 emission sources within the close proximity of the Ship Shoal Block 84 and 107 fields. All the offshore and onshore pipelines have been digitized and can be viewed in an interactive website (http://www.geomechanicstech.com/shipshoal.html). We also performed a feasibility study on the potential for converting existing oil or gas pipelines for CO 2 transport. There are abandoned, idled and retired onshore and offshore pipelines. However, there is no standard specification for maximum pipeline pressures needed for CO 2 transport. It is a function of design, materials, and testing provided; thus it will be the responsibility of the pipeline operator to correctly determine, maintain and operate within the limits of the pipelines. There are a few transit corridors extending from onshore Louisiana (Cailou Bay) to offshore trunk-lines. The cost for constructing a pipeline has increased about 46% from 2015 to an average cost of $5,064,046 per mile in 2016. We use the NETL CO 2 storage resource mass estimate to calculate the potential resources for Ship Shoal Block 84 and 107 fields using the depleted oil and gas reservoir volume and sand volume generated from our geologic modeling. The estimated storage resource results are greater when using the sand volumes derived through geologic modeling than the oil and gas reservoir data. The difference is due to the depleted oil and gas reservoir information not accounting for the water-flooded sand located below the oil/gas-water contact or unproductive sand units. Resource calculation using the sand volume obtained through geologic modeling overestimate the storage capacity as the model accounts for all sand within the formation, not just the interconnected sand.

02 PETROLEUM↗

State Estimation in Low-Observable Distribution Systems Using Matrix Completion: Preprint

The need for distribution system state estimation is on the rise because of the increased penetration of distributed energy resources and flexible load. To manage the distribution systems in real time, operators need to first overcome the challenge of low observability in distribution systems. Also, because of the amount of data present from smart meters, distributed generation measurements, switches, etc., the ideal distribution state estimation methods need to be able to process heterogeneous data. In this paper, an algorithm is developed for voltage phasor estimation in low-observability distribution systems. The algorithm is based on the matrix completion approach from signal processing. The traditional matrix completion formulation is augmented with power-flow constraints to improve results while requiring less data. This method can also use all types of measurements (voltage magnitude, voltage angle, real power, reactive power) to complete the state matrix.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Three-Dimensional Ignition and Flame Propagation Above Liquid Fuel Pools: Computational Analysis

A three-dimensional unsteady reactive Navier-Stokes code is developed to study the ignition and flame spread above liquid fuels initially below the flashpoint temperature. Opposed air flow to the flame spread due to forced and/or natural convection is considered. Pools of finite width and length are studied in air channels of prescribed height and width. Three-dimensional effects of the flame front near the edge of the pool are captured in the computation. The formation of a recirculation zone in the gas phase similar to that found in two-dimensional calculations is also present in the three-dimensional calculations. Both uniform spread and pulsating spread modes are found in the calculated results.

Cai, Jinsheng↗

Initial steady-state core simulation capability or thermal and pool-type molten salt reactors, coupling reactor physics, thermal-hydraulics, and evolving chemistry

This report presents the development and validation of an initial steady-state multiphysics capability for molten salt reactors (MSRs) under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in Fiscal Year 2025. The framework integrates neutronics, thermal-hydraulics, species transport, and thermochemistry to capture the coupled dynamics of liquid-fueled systems. Implementation and testing were performed on two representative designs: the Molten Salt Reactor Experiment (MSRE), a thermal-spectrum, channeled-flow reactor, and the Lotus Molten Salt Reactor (L-MSR), a fast-spectrum, pool-type reactor. The modeling suite employs Griffin for reactor physics and depletion, Pronghorn and SAM for thermal-hydraulics, Thermochimica for chemistry, and Saline for thermophysical properties, with benchmarking and validation carried out against historical MSRE data, experimental flow-loop measurements, and reference depletion calculations from Monte Carlo codes. The framework demonstrated the ability to reproduce key reactor behaviors including temperature feedback, reactivity losses, delayed neutron precursor transport, xenon poisoning, and redox potential evolution. The results confirm the feasibility and accuracy of the coupled models in predicting steady-state and selected transient MSR behaviors. This latter ones are used in this report as a proxy indicating that the steady-state models from which the transient starts are accurate. For MSRE, validation showed good agreement with pump start-up and natural circulation tests, while for the L-MSR, benchmarking confirmed hydraulic calibration and consistency of neutronics–thermal coupling. The tools also provided new insights into species transport, noble metal deposition, and salt solidification dynamics. On the Xenon transport front, the code is validated against the steady state Xenon poisoining measurement and showed good agreement with the experimental value. Identified areas for future work include advanced void transport modeling, three-dimensional simulations, improved alloy corrosion models, and tighter integration with high-fidelity Monte Carlo codes. These developments provide a foundation for high-fidelity MSR simulations that can support reactor design optimization, safety assessments, and long-term operational strategies.

42 - ENGINEERING↗

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗

Conservation Voltage Reduction (CVR) via Two-Timescale Control in Unbalanced Power Distribution Systems

Voltage control devices are employed in power distribution systems to reduce the power consumption by operating the system closer to the lower acceptable voltage limits; this technique is called conservation voltage reduction (CVR). The different modes of operation for system’s legacy devices (with binary control) and new devices (e.g. smart inverters with continuous control) coupled with variable photovoltaic (PV) generation results in voltage fluctuations which makes it challenging to achieve CVR objective. This paper presents a two-timescale control of feeder’s voltage control devices to achieve CVR that includes (1)a centralized controller operating in a slower time-scale to coordinate voltage control devices across the feeder and (2) local controllers operating in a faster timescale to mitigate voltage fluctuations due to PV variability. The centralized controller utilizes a three-phase optimal power flow model to obtain the decision variables for both legacy devices and smart inverters. The local controllers operate smart inverters to minimize voltage fluctuations and restore nodal voltages to their reference values by adjusting the reactive power support. The proposed approach is validated using the IEEE-123 bus (medium-size) and R3-12.47-2 (large-size) feeders. It is demonstrated that the proposed approach is effective in achieving the CVR objective for unbalanced distribution systems.

optimal control, distributed generation, fluctuati↗

NO 2 oxidation reactivity and burning mode of diesel particulates

The NO 2 oxidation kinetics and burning mode for diesel particulate from light-duty and medium-duty engines fueled with either ultra low sulfur diesel or soy methyl ester biodiesel blends have been investigated and are shown to be significantly different from oxidation by O 2 . Oxidation kinetics were measured using a flow-through packed bed microreactor for temperature programmed reactions and isothermal differential pulsed oxidation reactions. The burning mode was evaluated using the same reactor system for flowing BET specific surface area measurements and HR-TEM with fringe analysis to evaluate the nanostructure of the nascent and partially oxidized particulates. The low activation energy measured, specific surface area progression with extent of oxidation, HR-TEM images and difference plots of fringe length and tortuosity paint a consistent picture of higher reactivity for NO 2 , which reacts indiscriminately immediately upon contact with the surface, leading to the Zone I or shrinking core type oxidation. In comparison, O 2 oxidation is shown to have relatively lower reactivity, preferentially attacking highly curved lamella, which are more reactive due to bond strain, and short lamella, which have a higher proportion of more reactive edge sites. Furthermore, this preferential oxidation leads to Zone II type oxidation, where solid phase diffusion of oxygen via pores contributes significantly to slowing the overall oxidation rate, by comparison.

33 ADVANCED PROPULSION SYSTEMS↗

P-wave velocity anisotropy related to sealed fractures reactivation tracing the structural diagenesis in carbonates

Fracture properties are important in carbonate reservoir characterization, as they are responsible for a large part of the fluid transfer properties at all scales. It is especially true in tight rocks where the matrix transfer properties only slightly contribute to the fluid flow. Open fractures are known to strongly affect seismic velocities, amplitudes and anisotropy. Here, we explore the impact of fracture evolution on the geophysical signature and directional V p anisotropy of fractured carbonates through diagenesis. For that purpose, we studied a meter-scale, parallelepiped quarry block of limestone using a detailed structural and diagenetic characterization, and numerous V p measurements. The block is affected by two en-échelon fracture clusters, both being formed in opening mode (mode 1) and cemented, but only one being reactivated in shear. We compared the diagenetic evolution of the fractures, which are almost all 100% filled with successive calcite cements, with the P-wave velocities measured across this meter-scale block of carbonate, which recorded the tectonic and diagenetic changes of a South Provence sedimentary basin. We found that a directional V p anisotropy magnitude as high as 8–16% correlates with the reactivated fractures’ cluster dip angle, which is explained by the complex filling sequence and softer material present inside the fractures that have been reactivated during the basin's tectonic inversion. We show that although a late karstification phase preferentially affected these reactivated fractures, it only amplified the pre-existing anisotropy due to tectonic shear. We conclude that V p anisotropy measurements may help to identify the fracture sealing/opening processes associated with polyphased tectonic history, the anisotropy being independent of the current stress-state. Finally, this case shows that velocity anisotropies induced by fractures resulted here from a cause that is different from how these features have often been interpreted: selective reactivation of sealed fractures clusters rather than direction of currently open ones.

58 GEOSCIENCES↗

Determining the biogeochemical transformations of organic matter composition in rivers using molecular signatures

Inland waters are hotspots for biogeochemical activity, but the environmental and biological factors that govern the transformation of organic matter (OM) flowing through them are still poorly constrained. Here we evaluate data from a crowdsourced sampling campaign led by the Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) consortium to investigate broad continental-scale trends in OM composition compared to localized events that influence biogeochemical transformations. Samples from two different OM compartments, sediments and surface water, were collected from 97 streams throughout the Northern Hemisphere and analyzed to identify differences in biogeochemical processes involved in OM transformations. By using dimensional reduction techniques, we identified that putative biogeochemical transformations and microbial respiration rates vary across sediment and surface water along river continua independent of latitude (18°N–68°N). In contrast, we reveal small- and large-scale patterns in OM composition related to local (sediment vs. water column) and reach (stream order, latitude) characteristics. These patterns lay the foundation to modeling the linkage between ecological processes and biogeochemical signals. We further showed how spatial, physical, and biogeochemical factors influence the reactivity of the two OM pools in local reaches yet find emergent broad-scale patterns between OM concentrations and stream order. OM processing will likely change as hydrologic flow regimes shift and vertical mixing occurs on different spatial and temporal scales. As our planet continues to warm and the timing and magnitude of surface and subsurface flows shift, understanding changes in OM cycling across hydrologic systems is critical, given the unknown broad-scale responses and consequences for riverine OM.

rivers↗

Modelling uranium dioxide corrosion under repository conditions: A pore-scale study of the chemical and thermal processes

This paper investigates the corrosion behaviors of spent uranium dioxide (UO 2 ) fuel when placed in geological repositories. We performed pore-scale reactive transport simulations of the UO 2 corrosion process in a defective fuel rod with different orientations of breach or fracture on its clad. It is found that the corrosion rate has strong dependency on the breach orientation. The highest corrosion rate of UO 2 is calculated when the angle between the groundwater flow direction and the clad fracture reaches 180°. UO 2 corrosion simulations with different flow rates are compared. The results show the higher flow rate accelerates the corrosion of UO 2 fuel. The effect of pH on the corrosion process is also determined. It is indicated that higher pH can limit the corrosion of UO 2 by reducing the reaction rate. The dependence of reaction rate and time on reactive surface area is explored. Spent fuel with lower surface area demonstrates longer lifetime under corrosion conditions. A numerical model coupling thermal conduction and chemical reactions is developed to assess the impact of temperature on the process of UO 2 corrosion. The results show that higher temperature leads to larger corrosion rates for UO 2 . Here, the predicted reaction rates are higher in comparison with the isothermal results. Hence, our results can help improve the fundamental understanding of UO 2 corrosion in geological repositories for long-term storage of spent nuclear fuels and provide guidance for the safe operations and selection of appropriate repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Real-time process monitoring for direct ink write additive manufacturing

Direct ink write (DIW) printing of reactive resins presents a unique challenge due to the time-dependent nature of the rheological and chemical properties of the ink. As a result, careful print optimization or process control is important to obtain consistent, high quality prints. The present invention uses a flow-through characterization cell for in situ chemical monitoring of a resin ink during DIW printing. Additionally, in-line extrusion force monitoring can be combined with off-line post inspection using machine vision. By combining in-line spectroscopy and force monitoring, it is possible to follow reaction kinetics (for example, curing of a reactive resin) and viscosity changes during printing, which can be used for a closed-loop process control. Additionally, the capability of machine vision to automatically identify and quantify print artifacts can be incorporated on the printing line to enable real-time, AI-assisted quality control of the printed products. Together, these techniques can form the building blocks of an optimized process control strategy when complex reactive ink must be used to produce printed hardware.

Cook, Adam W.↗

LX-17 Thermal Decomposition-Characterization of Solid Residues from Cook-Off in a Small-Scale Vessel Under Confinement

Concerns surround whether insensitive (or any) energetic materials are more dangerous to handle when exposed to abnormal thermal environments. This study characterizes the residual material remaining after LX-17 (92.5 % 1,3,5-triamino 2,4,6-trinitro benzene (TATB) and 7.5 % Kel-F) is exposed to various thermal environments in a sealed small-scale vessel cook-off test reactor (heated at 0.1 to 100 °C/min until the reactor opened at 3000 psi (20.7 MPa)). Previous work has shown no additional sensitivity of these residues as evaluated by small-scale safety analysis, but characterization on the molecular scale indicates the TATB is transformed to more reactive compounds as well as the residue could be precursors to toxic gases. The solids and chars were characterized by various analytical methods. Heat-flow measurements indicated exothermic release is due to a mixture of residual TATB and related decomposition products (which may be more energetic). The N/C and O/N ratios indicated a material much more degraded than TATB. Primarily, the solids were a network of amorphous C inter-dispersed with N and O. Types of bonding include C–C, C–N, N–H, N–C, N=C, N≡C, C–O=, and –OH. Solvent extracts of the solids showed TATB decomposition intermediates benzo-furazans and benzo-furoxans, substituted TATB (mono-nitroso, hydroxyl, and chlorinated) along with several unidentified smaller molecules. These results indicate thermal treatment produces an amorphous carbon residue with heteroatoms incorporated through differing functionality, varying depending upon the thermal severity of exposure. These structures also could further decompose producing toxic light gases (such as cyanide).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fluid age-based analysis of a lifted turbulent DME jet flame DNS [Residence Time-Based Analysis of a Lifted Turbulent DME Jet Flame]

The link between the distribution of fluid residence time and the distribution of reactive scalars is analysed using Direct Numerical Simulation data. Information about the reactive scalar distribution is needed in order to model the reaction terms that appear in Large Eddy and Reynolds-Averaged simulations of turbulent reacting flows. The lifted flame is simulated taking account of multi-step chemistry for dimethyl-ether fuel and differential diffusion. Due to autoignition and flame propagation, the reaction progress increases with residence time. The variation of fluid residence time is evaluated by solving an Eulerian transport equation for the fluid age. The fluid age is a passive scalar with a spatially-uniform source term, meaning that its moments and dissipation rates in turbulent flows can be modelled using closures already established for conserved scalars such as mixture fraction. In combination with the mixture fraction, the fluid age serves as a useful mapping variable to distinguish younger less-reacted fluid near the inlet from older more-reacted fluid downstream. The local fluctuations of mixture fraction and fluid age have strong negative correlation and, building upon established presumed-pdf models for mixture fraction, this feature can be used to construct an accurate presumed-pdf model for the joint mixture fraction/fluid age pdf. It is demonstrated that the double-conditional first-order moment closure combined with the proposed presumed model for the joint pdf of mixture fraction and fluid age gives accurate predictions for unconditional reaction rates – both for pre-ignition radical species produced by low-temperature processes upstream of the flame base, and for major species that are produced at the flame front.

42 ENGINEERING↗