Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Bounce- and MLT-averaged diffusion coefficients in a physics-based magnetic field geometry obtained from RAM-SCB for the March 17 2013 storm

Local acceleration via whistler wave and particle interaction plays a significant role in particle dynamics in the radiation belt. In this work we explore gyro-resonant wave-particle interaction and quasi-linear diffusion in different magnetic field configurations related to the March 17 2013 storm. We consider the Earth's magnetic dipole field as a reference and compare the results against non-dipole field configurations corresponding to quiet and stormy conditions. The latter are obtained with the ring current-atmosphere interactions model with a self-consistent magnetic field RAM-SCB, a code that models the Earth's ring current and provides a realistic modeling of the Earth's magnetic field. By applying quasi-linear theory, the bounce- and MLT-averaged electron pitch angle, mixed term, and energy diffusion coefficients are calculated for each magnetic field configuration. For radiation belt (~1 MeV) and ring current (~100 keV) electrons, it is shown that at some MLTs the bounce-averaged diffusion coefficients become rather insensitive to the details of the magnetic field configuration, while at other MLTs storm conditions can expand the range of equatorial pitch angles where gyro-resonant diffusion occurs and significantly enhance the diffusion rates. When MLT average is performed at drift shell L = 4.25 (a good approximation to drift average), the diffusion coefficients become quite independent of the magnetic field configuration for relativistic electrons, while the opposite is true for lower energy electrons. These results suggest that, at least for the March 17 2013 storm and for L ≲ 4.25, the commonly adopted dipole approximation of the Earth's magnetic field can be safely used for radiation belt electrons, while a realistic modeling of the magnetic field configuration is necessary to describe adequately the diffusion rates of ring current electrons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Computational modeling study of the radial line slot antenna microwave plasma source with comparisons to experiments

The radial line slot antenna plasma source is a high-density microwave plasma source comprising a high electron temperature source region within the plasma skin depth from a coupling window and low electron temperature diffusion region far from the window. The plasma is typically comprised of inert gases like argon and mixtures of halogen or fluorocarbon gases for etching. Following the experimental study of Tian et al.[J. Vac. Sci. Technol. A 24, 1421 (2006)], a two-dimensional computational model is used to describe the essential features of the source. A high density argon plasma is described using the quasi-neutral approximation and coupled to a frequency-domain electromagnetic wave solver to describe the plasma-microwave interactions in the source. The plasma is described using a multispecies plasma chemistry mechanism developed specifically for microwave excitation conditions. The plasma is nonlocal by nature with locations of peak power deposition and peak plasma density being very different. The spatial distribution of microwave power coupling depends on whether the plasma is under- or over-dense and is described well by the model. The model predicts the experimentally observed low-order diffusion mode radial plasma profiles. The trends of spatial profiles of electron density and electron temperature over a wide range of power and pressure conditions compare well with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

THE ROLE OF THE MAGNETOROTATIONAL INSTABILITY IN MASSIVE STARS

The magnetorotational instability (MRI) is key to physics in accretion disks and is widely considered to play some role in massive star core collapse. Models of rotating massive stars naturally develop very strong shear at composition boundaries, a necessary condition for MRI instability, and the MRI is subject to triply diffusive destabilizing effects in radiative regions. We have used the MESA stellar evolution code to compute magnetic effects due to the Spruit-Tayler (ST) mechanism and the MRI, separately and together, in a sample of massive star models. We find that the MRI can be active in the later stages of massive star evolution, leading to mixing effects that are not captured in models that neglect the MRI. The MRI and related magnetorotational effects can move models of given zero-age main sequence mass across ''boundaries'' from degenerate CO cores to degenerate O/Ne/Mg cores and from degenerate O/Ne/Mg cores to iron cores, thus affecting the final evolution and the physics of core collapse. The MRI acting alone can slow the rotation of the inner core in general agreement with the observed ''initial'' rotation rates of pulsars. The MRI analysis suggests that localized fields ∼10{sup 12} G may exist at the boundary of the iron core. With both the ST and MRI mechanisms active in the 20 M {sub ☉} model, we find that the helium shell mixes entirely out into the envelope. Enhanced mixing could yield a population of yellow or even blue supergiant supernova progenitors that would not be standard SN IIP.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Effects of various parameters of different porous transport layers in proton exchange membrane water electrolysis

Porous transport layers (PTLs) play an important role in proton exchange membrane water electrolysis (PEMWE) cells. The PTL facilitates water and gas transport, as well as thermal and electrical conduction, and is required to sustain good contact with adjacent components. It is expected that using PTLs with variations in material properties such as structure, composition, thickness and wettability results in performance changes of the PEMWE. Here, a general mathematical PEMWE model is developed that separates and analyzes the contributions of ohmic, activation, diffusion and Nernst potentials. For model validation, three inherently different anode PTL structures (carbon paper, sintered titanium particles, and titanium felt) are operated over a range of conditions. Additionally, the effects of PTL wettability were used to verify the model using Polytetrafluoroethylene (PTFE) treated Toray papers with PTFE loading ranging from 0% to 20%. The modeling results of both PTFE treated and untreated materials show good agreement with the experimental data. Mass transport or diffusion loss is the primary reason for performance differences between PTFE treated and untreated PTLs. Sintered titanium PTLs with thicknesses above 1 mm suffer from up to 33% increased ohmic losses without indicating any obvious changes in activation and diffusion losses when compared to untreated PTLs. The losses of the cell increase when using PTFE treated Toray paper. Individual contributions are quantified and assigned to increased ohmic, activation, and diffusion losses. In conclusion, the proposed model offers insights into the overpotential contributions of a PEMWE. It is a useful tool for predicting performance of various PTL materials and can be applied for PTL development and optimization efforts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗
DOE OSTI2013

A coupled mechanical and chemical damage model for concrete affected by alkali–silica reaction

To model the complex degradation phenomena occurring in concrete affected by alkali–silica reaction (ASR), we formulate a poro-mechanical model with two isotropic internal variables: the chemical and the mechanical damage. The chemical damage, related to the evolution of the reaction, is caused by the pressure generated by the expanding ASR gel on the solid concrete skeleton. The mechanical damage describes the strength and stiffness degradation induced by the external loads. As suggested by experimental results, degradation due to ASR is considered to be localized around reactive sites. The effect of the degree of saturation and of the temperature on the reaction development is also modeled. The chemical damage evolution is calibrated using the value of the gel pressure estimated by applying the electrical diffuse double-layer theory to experimental values of the surface charge density in ASR gel specimens reported in the literature. The chemo-damage model is first validated by simulating expansion tests on reactive specimens and beams; the coupled chemo-mechanical damage model is then employed to simulate compression and flexure tests results also taken from the literature. -- Highlights: •Concrete degradation due to ASR in variable environmental conditions is modeled. •Two isotropic internal variables – chemical and mechanical damage – are introduced. •The value of the swelling pressure is estimated by the diffuse double layer theory. •A simplified scheme is proposed to relate macro- and microscopic properties. •The chemo-mechanical damage model is validated by simulating tests in literature.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Modeling the Ionic Strength Effect on Diffusion in Clay. The DR-A Experiment at Mont Terri

Solute diffusion in compacted clays depends on ionic strength through its control on the thickness of the electrical double layer (EDL) on the charged clay surfaces. In the DR-A field experiment (Mont Terri, Switzerland), synthetic porewater was circulated through a borehole for 189 days, leading to the out-diffusion of a variety of tracers into the Opalinus Clay. The borehole solution was then replaced with a higher-salinity solution for an additional 540 days, leading to the diffusion of Cs + , Ca 2+ , Mg 2+ , and Sr 2+ back into the borehole and to an increase in the out-diffusion of anions (I - , Br - ) and 3 H. The experimental results were interpreted using the CrunchClay code, which includes a mean electrostatic potential model for the EDL. The EDL corresponds to a second continuum in addition to bulk electrically neutral porewater. Species-specific diffusion (Nernst-Planck equation) occurs through both domains. A 1D radial model considered a single pore diffusion coefficient (Dp = 10-9 m2/s) for cations and 3H in the bulk porosity, and a smaller D p (3 × 10 -10 m 2 /s) for anions. Dp values in the EDL were smaller (10 -11 m 2 /s), except for Cs + and K + (5 × 10 -10 and 2 × 10 -10 m 2 /s, respectively). The model reproduced well the experimental results and showed the capability to consider temporal changes in geochemical conditions affecting the transport and retention of potentially important radionuclide contaminants (e.g., 137 Cs + , 90 Sr 2+ , 129 I - ) in underground geological nuclear waste repositories. Finally, coupled multicomponent diffusion together with the electrostatic properties of the charged surfaces are essential in the development of predictive models for ion transport in clays.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A Comparison of Adverse Meteorology Inputs for Emergency Management and Safety Basis Dispersion Modeling Applications at SRS

The MELCOR Accident Consequence Calculation System (MACCS), used for safety basis accident analyses at SRS, implements adjustment factors to consider the increased vertical turbulence associated with local surface roughness and the increased lateral turbulence associated with time-based plume meander. The two models used for SRS emergency planning radiological consequence assessment, i.e., Hotspot and Puff/Plume, currently do not have an option to input or apply diffusion adjustments for meander or surface roughness. Adverse meteorology conditions at SRS have been determined as E stability class and 1.3 meters per second wind speed for ground level releases. This is the input used for emergency planning models. Contrastingly safety basis analyses use a different input, representing adverse meteorology values as F stability class and wind speeds of 1.3 meters per second (for a ground level release). This results in an apparent discrepancy between the two approaches. In this report we document the differing modeling assumptions, discuss the logic behind each modeling approach used to represent local roughness and plume meander conditions, and numerically compare dispersion coefficients and concentration estimates obtained from each set of model inputs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Physics-Informed Deep Neural Networks for Learning Parameters and Constitutive Relationships in Subsurface Flow Problems

In this paper, we present a physics informed deep neural network (DNN) method for estimating parameters and unknown physics (constitutive relationships) in partial differential equation (PDE) models. We use PDEs in addition to measurements to train DNNs to approximate unknown parameters and constitutive relationships as well as states. The proposed approach increases the accuracy of DNN approximations of partially known functions when a limited number of measurements is available and allows for training DNNs when no direct measurements of the functions of interest are available. We employ physics informed DNNs to estimate the unknown space-dependent diffusion coefficient in a linear diffusion equation and an unknown constitutive relationship in a non-linear diffusion equation. For the parameter estimation problem, we assume that partial measurements of the coefficient and states are available and demonstrate that under these conditions, the proposed method is more accurate than state-of-the-art methods. For the non-linear diffusion PDE model with a fully unknown constitutive relationship (i.e., no measurements of constitutive relationship are available), the physics informed DNN method can accurately estimate the non-linear constitutive relationship based on state measurements only. Finally, we demonstrate that the proposed method remains accurate in the presence of measurement noise.

42 ENGINEERING↗
DOE OSTI2026

Intrinsic kinetics of water-inhibited ultra-lean methane oxidation over PtPd-Mg/ θ -Al 2 O 3 catalyst

Here, this study develops intrinsic methane oxidation kinetics for ultra-lean methane conditions in the presence of water over a highly active and stable PtPd–Mg/θ-Al 2 O 3 catalyst. Comprehensive laboratory experiments were conducted over a wide range of methane concentrations (150–1200 ppm CH 4 ), water contents (1–5% H 2 O), and industrially relevant space velocities (80,000 ≤ GHSV ≤ 160,000 h -1 ). These systematic experiments informed a two-dimensional, axisymmetric, multiscale reactor model that was used to develop and validate methane oxidation kinetics under practically relevant conditions, including non-isothermal operation and high conversion regimes. Combined experimental and modeling results revealed significant intraparticle diffusion resistance and transport-induced reaction exotherm at elevated temperatures, which limited catalyst utilization despite high intrinsic activity. These transport effects were explicitly incorporated into the reactor model, enabling accurate estimation of intrinsic kinetic parameters without reliance on conventional effectiveness-factor corrections. The resulting kinetic model successfully captured both kinetically controlled and mass-transfer-limited regimes and reliably predicted CH 4 conversion across broad ranges of temperature, methane concentration, and water content

Heat and mass transfer limitations↗
DOE OSTI2017

Surface transport processes in charged porous media

Surface transport processes are important in chemistry, colloidal sciences, engineering, biology, and geophysics. Natural or externally produced charges on surfaces create electrical double layers (EDLs) at the solid-liquid interface. The existence of the EDLs produces several complex processes including bulk and surface transport of ions. In this work, a model is presented to simulate bulk and transport processes in homogeneous porous media comprising big pores. It is based on a theory for capacitive charging by ideally polarizable porous electrodes without Faradaic reactions or specific adsorption of ions. A volume averaging technique is used to derive the averaged transport equations in the limit of thin electrical double layers. Description of the EDL between the electrolyte solution and the charged wall is accomplished using the Gouy-Chapman-Stern (GCS) model. The surface transport terms enter into the average equations due to the use of boundary conditions for diffuse interfaces. Two extra surface transports terms appear in the closed average equations. One is a surface diffusion term equivalent to the transport process in non-charged porous media. The second surface transport term is a migration term unique to charged porous media. The effective bulk and transport parameters for isotropic porous media are calculated solving the corresponding closure problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Atomistic modeling of dropwise condensation

The basic aim of the atomistic modeling of condensation of water is to determine the size of the stable cluster and connect phenomena occurring at atomic scale to the macroscale. In this paper, a population balance model is described in terms of the rate equations to obtain the number density distribution of the resulting clusters. The residence time is taken to be large enough so that sufficient time is available for all the adatoms existing in vapor-phase to loose their latent heat and get condensed. The simulation assumes clusters of a given size to be formed from clusters of smaller sizes, but not by the disintegration of the larger clusters. The largest stable cluster size in the number density distribution is taken to be representative of the minimum drop radius formed in a dropwise condensation process. A numerical confirmation of this result against predictions based on a thermodynamic model has been obtained. Results show that the number density distribution is sensitive to the surface diffusion coefficient and the rate of vapor flux impinging on the substrate. The minimum drop radius increases with the diffusion coefficient and the impinging vapor flux; however, the dependence is weak. The minimum drop radius predicted from thermodynamic considerations matches the prediction of the cluster model, though the former does not take into account the effect of the surface properties on the nucleation phenomena. For a chemically passive surface, the diffusion coefficient and the residence time are dependent on the surface texture via the coefficient of friction. Thus, physical texturing provides a means of changing, within limits, the minimum drop radius. The study reveals that surface texturing at the scale of the minimum drop radius does not provide controllability of the macro-scale dropwise condensation at large timescales when a dynamic steady-state is reached.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Simulating Uranium(VI) Diffusion in Sodium-Montmorillonite as a Function of Chemical Solution Conditions

Most nuclear waste disposal options include compacted bentonite, predominately consisting of montmorillonite clay, as part of a barrier system to minimize contaminant mobility. Uranium (U) is the primary element in spent nuclear fuel, and a potential contaminant of water resources. Uranium sorption onto clay and its slow diffusive transport away from waste canisters are expected to limit U(VI) mobility in these systems. However, their prediction is complicated by the complex mineralogical structure of montmorillonite clay, leading to two types of surface sites and clay porosities, and by the dependence of U(VI) solution speciation on chemical conditions. For instance, depending on pH and ionic strength, U(VI) may sorb onto montmorillonite due to weak ion exchange reactions or the formation of stable surface complexes. Furthermore, at the high degrees of clay compaction expected in future barrier systems, the ability of U(VI) species to access small, highly negatively charged interlayer spaces may ultimately control diffusive fluxes. In this modeling study, we investigate how the charge of U(VI) solution species affects U(VI) sorption and diffusion behavior. For this purpose, we simplified complex chemical solution speciation diagrams to charge speciation diagrams. Furthermore, we coupled an existing U(VI) sorption model (in PHREEQC) with U(VI) diffusion calculations (in GLE) to simulate U(VI) mobility in lab-scale through-diffusion experiments over various chemical conditions. Our results suggest that there are two important parameters that govern U(VI) diffusion in these systems. Uranium(VI) surface complexation on clay edge surfaces leads to contaminant retardation, while the access of U(VI) solution species to diffuse layers in clay interlayer spaces controls diffusive fluxes. Hence on the field scale, these parameters will determine how early a plume may cross site boundaries and pose risks to drinking water resources, and the levels of U(VI) concentrations to be expected in contaminated waters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Helium diffusivity in oxide nuclear fuel: Critical data analysis and new correlations

Helium is relevant in determining nuclear fuel behaviour. It affects the performance of nuclear fuel both in reactor and in storage conditions. Helium becomes important in reactor conditions when high burnups are targeted or MOX fuel is used, whereas for storage conditions it can represent a threat to the fuel rods integrity. The accurate knowledge of helium behaviour combined with predictive model capabilities is fundamental for the safe management of nuclear fuel, with helium diffusivity being a critical property. For this reason, a considerable number of separate effect experiments in the last fifty years investigated helium diffusivity in nuclear fuel. The aim of this work is to critically review and assess the experimental results concerning the helium diffusivity. Experimental results are critically analysed in terms of the helium introduction technique used (either infusion, implantation or doping) and of sample characteristics (single crystal, poly-crystal or powder). Accordingly, we derived two different correlations for the diffusivity. Clearly, each of the new correlations corresponds to a limited range of application conditions, depending on the experimental data used to derive it. We provide recommendations regarding the proper application conditions for each correlation (e.g., in reactor or storage conditions).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

SUNSPOT AND STARSPOT LIFETIMES IN A TURBULENT EROSION MODEL

Quantitative models of sunspot and starspot decay predict the timescale of magnetic diffusion and may yield important constraints in stellar dynamo models. Motivated by recent measurements of starspot lifetimes, we investigate the disintegration of a magnetic flux tube by nonlinear diffusion. Previous theoretical studies are extended by considering two physically motivated functional forms for the nonlinear diffusion coefficient D : an inverse power-law dependence D ∝ B {sup −ν} and a step-function dependence of D on the magnetic field magnitude B . Analytical self-similar solutions are presented for the power-law case, including solutions exhibiting “super fast” diffusion. For the step-function case, the heat-balance integral method yields approximate solutions, valid for moderately suppressed diffusion in the spot. The accuracy of the resulting solutions is confirmed numerically, using a method which provides an accurate description of long-time evolution by imposing boundary conditions at infinite distance from the spot. The new models may allow insight into the differences and similarities between sunspots and starspots.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Derivation of an Analytical Solution to a Reaction-Diffusion Model for Autocatalytic Degradation and Erosion in Polymer Microspheres

A mathematical reaction-diffusion model is defined to describe the gradual decomposition of polymer microspheres composed of poly(D,L-lactic-co-glycolic acid) (PLGA) that are used for pharmaceutical drug delivery over extended periods of time. The partial differential equation (PDE) model treats simultaneous first-order generation due to chemical reaction and diffusion of reaction products in spherical geometry to capture the microsphere-size dependent effects of autocatalysis on PLGA erosion that occurs when the microspheres are exposed to aqueous media such as biological fluids. The model is solved analytically for the concentration of the autocatalytic carboxylic acid end groups of the polymer chains that comprise the microspheres as a function of radial position and time. The analytical solution for the reaction and transport of the autocatalytic chemical species is useful for predicting the conditions under which drug release from PLGA microspheres transitions from diffusion-controlled to erosion-controlled release, for understanding the dynamic coupling between the PLGA degradation and erosion mechanisms, and for designing drug release particles. The model is the first to provide an analytical prediction for the dynamics and spatial heterogeneities of PLGA degradation and erosion within a spherical particle. The analytical solution is applicable to other spherical systems with simultaneous diffusive transport and first-order generation by reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Enhanced mobility of cations and anions in the redox state: The polaronium mechanism

Hugely enhanced slow-ion diffusivity has been widely observed under extreme redox conditions and for unclear reasons. Aided by first-principles calculations on model systems of ZrO 2 , CeO 2 , BaTiO 3 and Li 4/3 Mn 2/3 O 2 , here we successfully explained the intriguing phenomenon by a polaronium mechanism. We found a polaronium, defined as a transitory complex of a polaronic electron or hole and a migrating counterion, becomes highly mobile when the counterion comes from a d 0 or f 0 cation (e.g., Zr 4+ , Ti 4+ and Ce 4 ) or a p 6 anion (e.g., O 2- ) in the host compound. Upon a redox reaction, the complex attains a d 1 /f 1 or p 5 configuration, which spontaneously forms because it is favored by an electron-phonon interaction (manifest as the Jahn-Teller effect in high symmetry systems) that enables local relaxation and lowers the system energy. Here our calculations found such interaction reaching its peak at the saddle point where the local environment is softest, so soft that it allows a reorientation of the anisotropic d/f/p orbital to minimize the electron repulsion locally. Since the complex may dissolve after a successful ion-migration event, the redox electron/hole can be recycled to form another free-radical-like polaronium elsewhere, thereby enhancing ion migration repeatedly. The proposed polaronium mechanism, which also operates in ceramics doped with mixed-valence cations, is most relevant under dynamic and extreme thermal/field/irradiation conditions where extra electrons/holes are abundantly generated by non-equilibrium redox reactions. For such operations, some with emerging applications, our diffusion-enhancing mechanism may provide new theoretical insight to help understand their material/microstructure stability and performance.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A DNS evaluation of mixing models for transported PDF modelling of turbulent nonpremixed flames

Transported probability density function (TPDF) approaches are well suited to modelling turbulent, reacting, variable density flows. One of the main challenges to the successful deployment of TPDF methods is accurately modelling the unclosed molecular mixing term. This study examines three of the most widely used mixing models: the Interaction by Exchange with the Mean (IEM), Modified Curl (MC) and Euclidean Minimum Spanning Tree (EMST) models. Direct numerical simulation (DNS) data-sets were used to provide both initial conditions and inputs needed over the course of the runs, including the mean flow velocities, mixing frequency, and the turbulent diffusion coefficient. The same chemical mechanism and thermodynamic properties were used, allowing the study to focus on the mixing model. The simulation scenario was a one-dimensional, nonpremixed, turbulent jet flame burning either a syngas or ethylene fuel stream that featured extinction and reignition. This test scenario was selected because extinction and reignition phenomena are sensitive to the mixing model. Three DNS cases were considered for both the syngas and ethylene cases with a parametric variation of Reynolds and Damköhler numbers, respectively. Extinction events became more popular with increasing Reynolds number in the syngas cases and with decreasing Damköhler number in the ethylene cases. The model was first tested with the mixing frequency defined from the dissipation rate and variance of mixture fraction. With this definition, for the syngas cases this study finds that the TPDF method is successful at predicting flame extinction and reignition using all three mixing models for the relatively lower and intermediate Reynolds number cases, but that all models under-predict reignition in the relatively higher Reynolds number case. In the ethylene fuelled cases, only the EMST mixing model correctly predicts the reignition event for the two higher Damköhler number cases, yet, in the lowest Damköhler number case the EMST model over-predicts reignition and the IEM and MC models under-predict it. Mixing frequency was then modelled based on the turbulence frequency and a model constant $C_{\phi}$, the ratio of scalar to mechanical mixing rates. The DNS cases were reexamined with this definition and the findings suggested that the optimal value for $C_{\phi}$ is mixing model and case dependent. In particular, it was found in the ethylene case considered that reignition could be achieved with the IEM and MC models by adjusting the value of $C_{\phi}$.

42 ENGINEERING↗