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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Q2 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q2: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $\circ$ Run CESOL with BOUT++/Hermes-3 and immersed boundary condition to directly map to wall: • Run BOUT++/Hermes-3 through the IPS workflow to find radial particle and energy diffusivities to match either the Eich or the physics-based scaling of the SOL heat flux width, and • Expand source of first wall heat flux to include charged particles, neutrals, and radiation from the core+edge. 2. Generate medium fidelity parametrized CAD: $\circ$ Develop the TRACER tool to read an existing CAD, regenerate the geometry based on vertex location and connectivity information, define vertex translation and parameters needed for scaling the CAD, and $\circ$ Utilize the FreeGS code to determine CAT PF coil placement, including minimizing the number of coils, coil current, and electromechanical stresses. 3. Utilize plasma loading for engineering analysis: $\circ$ Couple the plasma loading to input for OpenFOAM and demonstrate initial test of thermal analysis of CAT first wall loading with typical DCLL blanket component cooling boundary conditions. 4. Demonstrate nuclear analysis: $\circ$ Apply initial analysis of tritium transport in DCLL blanket by evaluating spatially resolved tritium generation rates, tritium diffusion and convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗
DOE OSTI2019

On Thermodynamic and Kinetic Mechanisms for Stabilizing Surface Solid Solutions

Many processes for energy storage rely on transformations between phases with strong separation tendencies. In these systems, performance limitations can arise from undesirable chemical and mechanical factors associated with the phase separation behavior. Solid solutions represent a desirable alternative, provided the conditions for their formation are known. In the present work, we invoke linear stability theory and diffuse-interface mesoscopic simulations to demonstrate that solid solutions can be stabilized near surface layers of phase-separating systems. Two factors are found to drive surface solid-solution formation: surface relaxation of solution self-strain energy and anisotropy of diffusion mobility. Using a strongly phase-separating LiXFePO4 particle as a model system, we show that the relaxation of the solution self-strain energy competes against the relaxation of the coherency strain energy to stabilize surface solid solutions. Our theoretical understanding also suggests that highly anisotropic diffusion mobility can provide an alternative kinetic route to achieve the same aim, with stabilizing behavior strongly dependent on the specific alignment of the surface orientation. Our findings provide fundamental guidance for manipulating solid-solution behavior in nanoscale structures, in which surface effects become especially significant. Beyond energy storage materials, our findings have important implications for understanding solid-solution formation in other phase-separating systems from metal alloys to ceramics.

25 ENERGY STORAGE↗
DOE OSTI2018

Improving accuracy of electrochemical capacitance and solvation energetics in first-principles calculations

Here, reliable first-principles calculations of electrochemical processes require accurate prediction of the interfacial capacitance, a challenge for current computationally efficient continuum solvation methodologies. We develop a model for the double layer of a metallic electrode that reproduces the features of the experimental capacitance of Ag(100) in a non-adsorbing, aqueous electrolyte, including a broad hump in the capacitance near the potential of zero charge and a dip in the capacitance under conditions of low ionic strength. Using this model, we identify the necessary characteristics of a solvation model suitable for first-principles electrochemistry of metal surfaces in non-adsorbing, aqueous electrolytes: dielectric and ionic nonlinearity, and a dielectric-only region at the interface. The dielectric nonlinearity, caused by the saturation of dipole rotational response in water, creates the capacitance hump, while ionic nonlinearity, caused by the compactness of the diffuse layer, generates the capacitance dip seen at low ionic strength. We show that none of the previously developed solvation models simultaneously meet all these criteria. We design the nonlinear electrochemical soft-sphere solvation model which both captures the capacitance features observed experimentally and serves as a general-purpose continuum solvation model.

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

AB INITIO SIMULATIONS FOR MATERIAL PROPERTIES ALONG THE JUPITER ADIABAT

We determine basic thermodynamic and transport properties of hydrogen-helium-water mixtures for the extreme conditions along Jupiter's adiabat via ab initio simulations, which are compiled in an accurate and consistent data set. In particular, we calculate the electrical and thermal conductivity, the shear and longitudinal viscosity, and diffusion coefficients of the nuclei. We present results for associated quantities like the magnetic and thermal diffusivity and the kinematic shear viscosity along an adiabat that is taken from a state-of-the-art interior structure model. Furthermore, the heat capacities, the thermal expansion coefficient, the isothermal compressibility, the Grueneisen parameter, and the speed of sound are calculated. We find that the onset of dissociation and ionization of hydrogen at about 0.9 Jupiter radii marks a region where the material properties change drastically. In the deep interior, where the electrons are degenerate, many of the material properties remain relatively constant. Our ab initio data will serve as a robust foundation for applications that require accurate knowledge of the material properties in Jupiter's interior, e.g., models for the dynamo generation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Electrode Mesoscale as a Collection of Particles: Coupled Electrochemical and Mechanical Analysis of NMC Cathodes

Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.

42 ENGINEERING↗
DOE OSTI2023

Kernel Buffer Volume Fraction Margin of the AGR Designed Fuel Particle

Modeling results used to assess the fuel performance of the TRISO-coated fuel particles as a function of kernel/buffer volume fraction include SiC tangential stress, formation of the buffer/IPyC gap, particle temperature profile, internal particle pressure, fission gas released from the kernel, probability of fuel particle failure, and fission product diffusion. These results were evaluated at two burnup levels and irradiation temperatures to bound expected steady-state irradiation conditions. In general, increasing the kernel/buffer volume fraction increases the SiC stress and subsequently the failure probability of a fuel particle when compared to the AGR designed particle. There was little impact on the fission product diffusion through the particle as the kernel/buffer volume fraction increased.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗
DOE OSTI2024

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2014

Oxygen distribution in tumors: A qualitative analysis and modeling study providing a novel Monte Carlo approach

Purpose: To construct a Monte Carlo (MC)-based simulation model for analyzing the dependence of tumor oxygen distribution on different variables related to tumor vasculature [blood velocity, vessel-to-vessel proximity (vessel proximity), and inflowing oxygen partial pressure (pO{sub 2})]. Methods: A voxel-based tissue model containing parallel capillaries with square cross-sections (sides of 10 μm) was constructed. Green's function was used for diffusion calculations and Michaelis-Menten's kinetics to manage oxygen consumption. The model was tuned to approximately reproduce the oxygenational status of a renal carcinoma; the depth oxygenation curves (DOC) were fitted with an analytical expression to facilitate rapid MC simulations of tumor oxygen distribution. DOCs were simulated with three variables at three settings each (blood velocity, vessel proximity, and inflowing pO{sub 2}), which resulted in 27 combinations of conditions. To create a model that simulated variable oxygen distributions, the oxygen tension at a specific point was randomly sampled with trilinear interpolation in the dataset from the first simulation. Six correlations between blood velocity, vessel proximity, and inflowing pO{sub 2} were hypothesized. Variable models with correlated parameters were compared to each other and to a nonvariable, DOC-based model to evaluate the differences in simulated oxygen distributions and tumor radiosensitivities for different tumor sizes. Results: For tumors with radii ranging from 5 to 30 mm, the nonvariable DOC model tended to generate normal or log-normal oxygen distributions, with a cut-off at zero. The pO{sub 2} distributions simulated with the six-variable DOC models were quite different from the distributions generated with the nonvariable DOC model; in the former case the variable models simulated oxygen distributions that were more similar to in vivo results found in the literature. For larger tumors, the oxygen distributions became truncated in the lower end, due to anoxia, but smaller tumors showed undisturbed oxygen distributions. The six different models with correlated parameters generated three classes of oxygen distributions. The first was a hypothetical, negative covariance between vessel proximity and pO{sub 2} (VPO-C scenario); the second was a hypothetical positive covariance between vessel proximity and pO{sub 2} (VPO+C scenario); and the third was the hypothesis of no correlation between vessel proximity and pO{sub 2} (UP scenario). The VPO-C scenario produced a distinctly different oxygen distribution than the two other scenarios. The shape of the VPO-C scenario was similar to that of the nonvariable DOC model, and the larger the tumor, the greater the similarity between the two models. For all simulations, the mean oxygen tension decreased and the hypoxic fraction increased with tumor size. The absorbed dose required for definitive tumor control was highest for the VPO+C scenario, followed by the UP and VPO-C scenarios. Conclusions: A novel MC algorithm was presented which simulated oxygen distributions and radiation response for various biological parameter values. The analysis showed that the VPO-C scenario generated a clearly different oxygen distribution from the VPO+C scenario; the former exhibited a lower hypoxic fraction and higher radiosensitivity. In future studies, this modeling approach might be valuable for qualitative analyses of factors that affect oxygen distribution as well as analyses of specific experimental and clinical situations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Model-informed, Adaptive Physical Vapor Deposition to Fabricate Hierarchical Binary-alloy Thin-films

Designing next generation thin films, tailor-made for specific applications, relies on the availability of robust processing-structure-property relationships. Traditional structure zone diagrams are limited to low-dimensional mappings, with machine-learning methods only recently attempting to relate multiple processing parameters to the final microstructure. Despite this progress, structure-processing relationships are unknown for processing conditions that vary during thin-film deposition, limiting the range of microstructures and properties achievable. In this project, we employed a phase-field computational model combined with a genetic algorithm (GA) to identify and design time-dependent processing protocols that achieve tailor-made microstructures. We simulate the physical vapor deposition of a binary-alloy thin film by employing a phase-field model, where deposition rates and diffusivities are controlled via the genetic algorithm. Our GA-guided protocols achieve targeted microstructures with lateral and vertical concentration modulations, as well as more complex, hierarchical microstructures previously not described in simple structure zone diagrams. Our algorithm provides insight to experimentalists looking for additional avenues to design novel thin-film microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Water vapor oxidation of SiC layer in surrogate TRISO fuel particles

Under accidental conditions for high temperature gas-cooled reactors (HTGR), the SiC layer in tri-structural-isotropic (TRISO) fuel particles can be exposed to water vapor. In this study, oxidation behaviors of surrogate TRISO fuel particles were investigated in a He-20 vol% water vapor mixed atmosphere at temperatures up to 1600 °C. The growth of the crystalline oxide passivation layer with cracks and pores followed a parabolic law with time, where the maximum was 2.3 μm at 1600 °C. The oxide layer and the SiC surface under the oxide became flattened with increasing temperature, as a function of the silica viscosity and the diffusion path of water vapor. Volatilization of the oxide layer was analyzed using a mechanistic model that an inert gas in oxidizing atmospheres could influence the magnitude of volatilization. The fracture load and strength of the oxidized and thinned SiC layer were numerically estimated to decrease from 2.27 to 1.69 N and from 317 to 299 MPa, respectively, with the SiC thickness decrease from 35 to 32 μm. This prediction indicates that the oxidized SiC layer should retain fission products. Additionally, the mechanical integrity of each layer in the TRISO fuel particle after oxidation was evaluated. The results in this work provide important data for the safety analysis of accidental scenarios in HTGRs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hydrogen Diffusivity Measurements of YH 1.87 Moderator Material with Incoherent Quasielastic Neutron Scattering

Yttrium hydride is an optimal choice of moderator material for thermal neutron spectrum reactors requiring small core volumes. The Transformation Challenge Reactor (TCR) program at Oak Ridge National Laboratory (ORNL) aims to develop an additively manufactured microreactor core by leveraging recent advances in materials, manufacturing, data analytics, and high-fidelity modeling and simulation. Yttrium hydride was selected as the moderator material for TCR due to its high moderating power, which enables the spatially efficient reactor design and thermal stability that is superior when compared with zirconium hydride. Hydrogen desorption from yttrium hydride is expected at elevated temperatures, posing safety and performance concerns. As a lifetime component in advanced reactors, appropriate strategies to mitigate hydrogen release from yttrium hydride over long-term reactor operation are required. Understanding hydrogen mobility in yttrium hydride provides the scientific basis for developing a hydrogen barrier for hydride moderator and predicting the hydrogen redistribution within the material during various operational conditions. This study investigated the hydrogen motion in YH 1.87 as a function of temperature ranging from 200 to 1,173 K using incoherent quasielastic neutron scattering (IQNS) at the Spallation Neutron Source (SNS). The results are presented and discussed in terms of hydrogen self-diffusion coefficients, activation energy for translational hydrogen motion, and hydrogen jump distances. Based on the IQNS data, YH 1.87 is found to be stable at high temperatures with no translational H motion below 1,023 K within the IQNS detection limits. In the temperature range of 1,073–1,173 K, hydrogen diffusivity in YH 1.87 is $D = 4.57 X 10^{⁻3} [m^2/s] exp (⁻\frac{1.73[eV]}{RT}$). The hydrogen concentration’s dependence on hydrogen diffusivity is also discussed in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

A Straightforward Model for Quantifying Local pH Gradients Governing the Oxygen Evolution Reaction

The production and consumption of protons by an electrocatalyst will, under certain conditions, generate localized microenvironments with properties distinct from those of the bulk solution. These local properties are particularly impactful for reactions involving proton-coupled electron transfer, where the generation of locally basic or acidic environments may significantly influence the energy efficiency and reaction selectivity of the electrocatalyst. Whereas local pH environments have been observed and characterized in reductive half-reactions, including the CO 2 reduction and hydrogen evolution reactions, the incompatibility of conventional techniques and materials has limited studies in oxidative half-reactions, including the oxygen evolution reaction (OER), which provides the reducing equivalents for solar-to-fuels electrolysis. With the straightforward parameters bulk pH, buffer composition and pK a , and mass transport, we develop a model for describing local pH as a function of current density regardless of the microscopic details of the mechanism. Using an acid-stable PbO x OER catalyst, we observe the formation and dissipation of pH gradients during the OER and validate the model with voltammetric and potentiometric studies. Here, the model predicts how local acidic environments can develop over a narrow OER current density window, thus providing further motivation for the development of OER catalysts that are stable to acid, even when operating in basic aqueous conditions. More generally, the model is not restricted to the OER and is useful for determining the onset of local pH gradients for other electrocatalytic reactions that involve the consumption or generation of protons in energy conversion reactions.

Anions↗
DOE OSTI2017

Treatment of the polar coordinate singularity in axisymmetric wave propagation using high-order summation-by-parts operators on a staggered grid

In this paper, we develop a high-order finite difference scheme for axisymmetric wave propagation in a cylindrical conduit filled with a viscous fluid. The scheme is provably stable, and overcomes the difficulty of the polar coordinate singularity in the radial component of the diffusion operator. The finite difference approximation satisfies the principle of summation-by-parts (SBP), which is used to establish stability using the energy method. To treat the coordinate singularity without losing the SBP property of the scheme, a staggered grid is introduced and quadrature rules with weights set to zero at the endpoints are considered. Finally, the accuracy of the scheme is studied both for a model problem with periodic boundary conditions at the ends of the conduit and its practical utility is demonstrated by modeling acoustic-gravity waves in a magmatic conduit.

42 ENGINEERING↗
DOE OSTI2022

Failure analysis of solid oxide fuel cells nickel-yttria stabilized zirconia anode under siloxane contamination

In this study, the failure process of the solid oxide fuel cell (SOFC) Ni-yttria-stabilized zirconia (YSZ) anode is investigated with D4 siloxane (octamethylcyclotetrasiloxane) contamination. In order to evaluate the influence of the electrochemical reaction on the siloxane deposition process, the SOFC experiments were operated at open circuit voltage (OCV) and 50 mA cm -2 conditions at 800°C. During the failure process, electrochemical, morphology and exhaust gas component analysis testing are conducted at the critical points. An equivalent circuit model and corresponding microstructure parameter calculations for separated physicochemical processes were utilized for the quantitative analysis of the failure process. Further, the results confirm the siloxane chemical adsorption deposition mechanism proposed in previous work. As a result, the failure of the anode was attributed to the gas diffusion blockage by dense silicon dioxide layer formation. The anode failure process with siloxane contamination is faster when the anode is operated under polarization.

25 ENERGY STORAGE↗
DOE OSTI2017

Hydrogen self-diffusion in single crystal olivine and electrical conductivity of the Earth’s mantle

Nominally anhydrous minerals formed deep in the mantle and transported to the Earth’s surface contain tens to hundreds of ppm wt H 2 O, providing evidence for the presence of dissolved water in the Earth’s interior. Even at these low concentrations, H 2 O greatly affects the physico-chemical properties of mantle materials, governing planetary dynamics and evolution. The diffusion of hydrogen (H) controls the transport of H 2 O in the Earth’s upper mantle, but is not fully understood for olivine ((Mg, Fe) 2 SiO 4 ) the most abundant mineral in this region. Here we present new hydrogen self-diffusion coefficients in natural olivine single crystals that were determined at upper mantle conditions (2 GPa and 750–900 °C). Hydrogen self-diffusion is highly anisotropic, with values at 900 °C of 10 -10.9 , 10 -12.8 and 10 -11.9 m 2 /s along [100], [010] and [001] directions, respectively. Combined with the Nernst-Einstein relation, these diffusion results constrain the contribution of H to the electrical conductivity of olivine to be σH = 10 2.12 S/m·C H2O ·exp -187kJ/mol/(RT) . Comparisons between the model presented in this study and magnetotelluric measurements suggest that plausible H 2 O concentrations in the upper mantle (≤250 ppm wt) can account for high electrical conductivity values (10 -2 –10 -1 S/m) observed in the asthenosphere.

58 GEOSCIENCES↗
DOE OSTI2025

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪⁡(100) simulations rather than 𝒪⁡(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.

general physics↗