Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Directional Effect”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Advanced Instrumentation for Metal Additive Manufacturing

Laser powder bed fusion (LPBF) is the most widely used process for metal additive manufacturing (AM), particularly where complex geometries provide performance advantages unattainable with traditional manufacturing techniques. However, LPBF is highly sensitive to innate variability in both the powder spreading and fusion steps, often leading to defects such as pores that are difficult to detect yet significantly impair component mechanical properties and fatigue life. This thesis presents a range of novel instruments enabling both precise assessment of powder layer characteristics and in-situ thermal metrology of metal AM to advance the quality control of LPBF. First, leveraging a custom X-ray microscope and a radiation-transport model developed through this work, transmission X-ray imaging is used to study spreading of thin metal powder layers. Effective layer depth is directly mapped at a process-relevant size scale, surpassing optical techniques that can only estimate local deposition from layer surface topography. Layer packing density and quality are shown to be influenced by powder flowability and particle size relative to nominal powder layer thickness. Layer quality is additionally connected to the geometry of the spreading implement and its velocity. This technique and its presented findings enable pairing feedstocks with spreading strategies that create layers with consistent packing density and uniformity. Second, a twofold approach is employed to optically interrogate the laser fusion step of LPBF for observing signatures of defect formation. Aperture division multiplexing is conceptualized, providing for simultaneous laser delivery and high-fidelity infrared (IR) process monitoring through a common optic. In-situ microscopy at 50 μm spatial resolution and at mid-wave IR wavelengths is proven readily achievable with the first purpose-built optic of this type. Next, a bespoke imaging spectrometer, along with a temperature-emissivity separation technique, is used to retrieve accurate process temperatures over a 1000 K range. Data from these instruments are correlated to porosity as fine as 4.3 μm in two LPBF test artifacts, as verified using computed tomography (CT), establishing the viability of robust optically-based component qualification.

Penny, Ryan↗
DOE OSTI2024

Economic and Jobs Impacts of Point-Source Carbon Capture in Cement Industry – Case Study

The cement industry accounts for an estimated 8% of global CO2 emissions, which surpasses that of the entire aviation sector. In contrast with other industries, where CO2 emissions can be drastically reduced via electrification or fuels substitution, cement production releases CO2 as part of its process, during the calcination of carbonates to yield oxides. Thus, point-source carbon capture has become a key technology in the cement industry’s decarbonization. Apart from the expected environmental benefits, point-source carbon capture in the cement industry can yield important economic benefits and create jobs. The objective of this study was to perform a preliminary assessment of the economic and workforce impacts associated with the construction and operation of a point-source carbon capture retrofit of an existing cement production facility, using as a basis the data from a front-end engineering design (FEED) study to install a 3.9 million metric tons per year (Mtpy) CO2 capture facility at Holcim Ste Genevieve cement plant in Missouri, United States of America. The advanced carbon capture technology used in this FEED study was Air Liquide’s Cryocap™ FG carbon capture technology. The study evaluated the direct, indirect, and induced economic impacts of the construction, operation, and maintenance activities of the project over its lifespan. It also covered how the project will generate new jobs, their nature, and quantity, along with strategies to prepare the workforce. To perform this study, construction, operation, and maintenance cost estimates, as well as construction and operation staffing plans from the FEED study were input into IMPLAN version 7.5 software, licensed by IMPLAN Group LLC (Huntersville, VC), to predict the direct, indirect and induced economic impacts of the project using industry multipliers from the software. Additionally, recruitment strategies were developed for hiring individuals who belong to groups that are historically underserved or underrepresented, as well as anticipated recruitment of workers from the local community (whether training will be required or if the skills are associated with an existing labor force). The analysis estimated that the construction and operation of the carbon capture at Holcim Ste. Genevive will result in over 24 thousand work-years of job opportunities, close to USD 10 billion of economic impacts, including over USD 460 million of tax revenue. These results encompass the direct, indirect, and induced effects. A strategy to maximize hiring from the project and neighboring counties was developed, leveraging training agreements with local trade groups and universities. The result of this study can be used for a strategic preliminary assessment of the potential regional economic and job impacts of retrofitting existing cement plants with point source carbon systems, and its methodology can be replicated to individual projects to aid in planning and workforce development.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Conservation laws and effective hadronization models

Hadronization models based on local string-breaking dynamics are typically Markovian by construction, yet the physical ensemble of final states is shaped by global constraints that couple the entire fragmentation trajectory. Recasting hadronization as a conditioned stochastic diffusion process provides a precise mathematical resolution to this tension. In particular, this language reveals explicitly that constraints stemming from conservation laws induce non-Markovian correlations between otherwise independent fragmentation steps, and that these correlations can be absorbed exactly into a renormalization of the local dynamics through a Doob $h$-transform. We develop this formalism for a $q\bar{q}$ string in the chiral limit, where the longitudinal-transverse factorization of the Lund kernel becomes exact, enabling systematic power counting and clean ultraviolet (UV)/infrared (IR) separation. The dynamics organize naturally into a tower of effective theories distinguished by the remaining string mass, spanning a UV fixed point with scale-invariant transport coefficients, an intermediate regime where transverse phase space induces controlled running, and an IR boundary layer where non-local effects enter at leading order. The tower exhibits genuine Wilsonian structure, including $β$-functions, anomalous dimensions, and systematic matching conditions. The resulting framework achieves a clean factorization of universal microscopic fragmentation dynamics from infrared constraint effects, and opens new directions for both the theoretical analysis and practical simulation of hadronization.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗
DOE OSTI2021

Exposure, health effects, sensing, and remediation of the emerging PFAS contaminants – Scientific challenges and potential research directions

Per- and polyfluoroalkyl substances (PFAS) make up a large group of persistent anthropogenic chemicals which are difficult to degrade and/or destroy. PFAS are an emerging class of contaminants, but little is known about the long-term health effects related to exposure. In addition, technologies to identify levels of contamination in the environment and to remediate contaminated sites are currently inadequate. In this opinion-type discussion paper, a team of researchers from the University of Connecticut and the University at Albany discuss the scientific challenges in their specific but intertwined PFAS research areas, including rapid and low-cost detection, energy-saving remediation, the role of T helper cells in immunotoxicity, and the biochemical and molecular effects of PFAS among community residents with measurable PFAS concentrations. Potential research directions that may be employed to address those challenges and improve the understanding of sensing, remediation, exposure to, and health effects of PFAS are then presented. Here, we hope our account of emerging problems related to PFAS contamination will encourage a broad range of scientific experts to bring these research initiatives addressing PFAS into play on a national scale.

Challenges and research directions↗
DOE OSTI2026

Infrared triggered dwell and active cooling thermal control effects on microstructural uniformity in DED

Directed Energy Deposition (DED) offers rapid large scale fabrication, but difficulty in delivering consistent microstructures and properties hinders the use of DED fabricated components in safety or performance critical applications. Variability stems from the complex thermal cycles generated by the toolpath used to print the required geometry. Several practical methods have become established in DED to regulate overheating, such as active cooling of the baseplate structure or the use of an infrared camera to inject interlayer pauses to ensure the top layer of the component cools to a set temperature, which have been shown to affect microstructure. However, no critical assessment has been performed as to how effective these controls are in promoting microstructural uniformity in the context of complex layer timing commonly generated by non-prismatic geometries. Here we show how controls influence the thermal field, phase transformations, and dynamic annealing of a low-temperature transformation steel using infrared imaging and operando neutron diffraction. Counterintuitively, common thermal homogenization process controls can reduce microstructural uniformity because these approaches stabilize peak temperature while overlooking temperatures near the solid-state phase transformation fronts. Instead, the cyclic reheating induces spatially-variant dynamically annealed regions which can be modulated via control parameters. We show that these controls have spatially linked effects centimeters away from the active weld, which implies that microstructure control must co-optimize thermal input across many subsequent layers. In conclusion, our results demonstrate the pressing need for higher order controls that integrate predictive elements of simulation data to stabilize printed properties for future qualification of DED components.

Directed energy deposition↗
DOE OSTI2024

Nuclear Dependence of Charged Current Inclusive $\bar{\nu}_\mu -$A Cross Section Measurements at 6 GeV at MINERvA

The main advantage of the MINERvA detector is that it uses differently sized nuclear targets, in the same neutrino beam, allowing for a direct study of nuclear medium effects, with reduced flux and detector uncertainties. This thesis presents the first high statistics charged current antineutrino-nucleus inclusive two-dimensional cross section measurements using an antineutrino beam with $\langle E_{\bar\nu}\rangle\sim$ 6 GeV, on iron, lead, carbon and hydrocarbon targets, as a function of bjorken $x$ and four-momentum transfer squared $Q^2$, invariant mass $W$ and $Q^2$ and muon longitudinal and transverse momenta, $p_t$ and $p_z$. These nuclei of different sizes were utilized to measure the nuclear dependence of the cross section. The antineutrino interactions that were utilized for the analysis presented in this thesis, were having an antimuon in the final state with energy 2-20GeV and with an antimuon scattering angle less than $17^o$. This analysis represents one of the largest antineutrino datasets analyzed to date in the energy region of $\langle E_{\bar\nu} \rangle \sim$ 6 GeV. This direct evidence for the nuclear medium effects, can help benchmarking models, which can be used by future neutrino experiments like DUNE, by extrapolating these results to Argon target.

Gaur, Prameet Kumar↗
DOE OSTI2022

Fabrication of Polymer Membrane-Suspended Microstructures on Printed Circuit Boards

When microelectromechanical systems are manufactured in high densities, the necessary routing layer structures can become complex and thus expensive if executed in silicon technology. If, additionally, the microstructures are also to be fabricated suspended on membranes, complex thinning processes become necessary to form these often rather fragile membranes. Here, a novel fabrication concept is demonstrated that allows for the fabrication of membrane-suspended microstructures using standard microfabrication techniques directly integrated with printed circuit boards (PCB) by laminating polyimide (PI) films over cavities in the PCB. This results in cost-effective routing layer fabrication and direct compatibility with common electronic component standards such as SMT (surface-mount technology) and THT (through-hole technology). To illustrate the capabilities of the fabrication process, we fabricated an array of membrane suspended microheaters. Individual platinum microheater elements are suspended on PI membranes that are 25 μm thick and 500 μm in diameter. The fabricated PCB carries a 10x10 array of suspended microheaters on an area of 26,mm x 26,mm. Here, the characterization of the microheater array shows that a yield of 86% functional microheaters was achieved. An individual microheater shows an almost linear heating characteristic at least up to 200°C demonstrating the viability of the technique for creating functional membrane-suspended microstructures on a PCB. [2021-0226]

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Synergizing a Large Ordinary Nernst Effect and Axis‐Dependent Conduction Polarity in Flat Band KMgBi Crystals

The exploration of quantum materials in which an applied thermo/electrical/magnetic field along one crystallographic direction produces an anisotropic response has led to unique functionalities. Along these lines, KMgBi is a layered, narrow gap semiconductor near a critical state between multiple Dirac phases due to the presence of a flat band near the Fermi level. The valence band is highly anisotropic with minimal cross-plane dispersion, which, in combination with an isotropic conduction band, enables axis-dependent conduction polarity. Thermopower and Hall measurements indicate dominant p-type conduction along the cross-plane direction, and n-type conduction along the in-plane direction, leading to a significant zero-field transverse thermoelectric response when the heat flux is at an angle to the principal crystallographic directions. Additionally, a large Ordinary Nernst effect (ONE) is observed with an applied field. It arises from the ambipolar term in the Nernst effect, whereby the Lorentz force on electrons and holes makes them drift in opposite directions so that the resulting Nernst voltage becomes a function of the difference between their partial thermopowers, greatly enhancing the ONE. It is proven that axis-dependent polarity can synergistically enhance the ONE, in addition to leading to a zero-field transverse thermoelectric performance.

36 MATERIALS SCIENCE↗
DOE OSTI2022

CO 2 fertilization of terrestrial photosynthesis inferred from site to global scales

Significance The magnitude of the CO 2 fertilization effect on terrestrial photosynthesis is uncertain because it is not directly observed and is subject to confounding effects of climatic variability. We apply three well-established eco-evolutionary optimality theories of gas exchange and photosynthesis, constraining the main processes of CO 2 fertilization using measurable variables. Using this framework, we provide robust observationally inferred evidence that a strong CO 2 fertilization effect is detectable in globally distributed eddy covariance networks. Applying our method to upscale photosynthesis globally, we find that the magnitude of the CO 2 fertilization effect is comparable to its in situ counterpart but highlight the potential for substantial underestimation of this effect in tropical forests for many reflectance-based satellite photosynthesis products.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Reducing Sample Loss in Measurement of Heat of Vaporization of Ethanol/Gasoline Blends by Differential Scanning Calorimetry/Thermogravimetric Analysis

Higher gasoline Heat of Vaporization (HOV) can enable higher compression-ratio direct injection spark ignition engines by providing evaporative cooling that effectively increases fuel knock resistance. Methods to directly measure this fuel property in complex gasoline samples are not well developed. The objective of the present study was to further improve a Differential Scanning Calorimetry/Thermogravimetric Analyzer (DSC/TGA) method to measure the total and partial HOV of gasoline. Ten market gasoline samples were chosen to have a wide range of properties to assess the method's capability across the entire volatility range with an emphasis on understanding how well the method captures the initial 10 percent (%) of sample evaporation and how much sample is left unevaporated at the end of the experiment. Modifications to the sample preparation/introduction method as well as to the instrument itself were made to reduce initial sample losses which included measurements of the HOV at 10 degrees C and 5 degrees C (in a cold chamber) as well as under (uncontrolled) ambient conditions. Experimental results from the DSC/TGA were compared to calculated total HOV results based on Detailed Hydrocarbon Analysis (DHA). Results from the two methods agreed very well with the difference being 5% or less in almost every case. In addition, the repeatability of the method was investigated by analyzing samples in triplicate at the three temperatures investigated. One valuable conclusion from the study was that the lower temperatures of 10 degrees C and 5 degrees C enabled more reproducible measurements for both total and partial HOV. This improved precision may be caused by the fact that temperature control in the cold chamber was more reliable than the ambient laboratory temperature control. Additionally, cold chamber experiments, due to ergonomic limitations, did not allow for the use of a lid on the sample pan. The reproducibility of the evaporation rate was found to be highly dependent on the pan/lid fit which can vary significantly such that elimination of the lid improved measurement precision while operation at sub-ambient temperature slowed the evaporation rate. Results detailing sample preparation and instrument modifications as well as a detailed comparison of total and partial HOV results are presented.

ADVANCED PROPULSION SYSTEMS,BIOMASS FUELS↗
DOE OSTI2021

Magnetoelastic transition and magnetocaloric effect in induction melted Fe 100-x Rh x bulk alloys with x = 50, 51

Magnetoelastic transitions (METs) in bulk in nearly equiatomic Fe-Rh alloys produced by arc melting may show poor reproducibility related to insufficient chemical homogeneity and presence of impurity phases in variable concentrations. To better understand the synthesis conditions that reliably yield bulk FeRh materials with reproducible MET characteristics, Fe 100-x Rh x alloys with x = 50, 50.5 and 51 at. % were prepared by induction melting and thermal annealing under identical conditions. The fabricated samples were cut into several slices, followed by characterization of METs in each of the slices using isothermal and isofield magnetization measurements, differential scanning calorimetry, and direct measurements of the magnetocaloric effect. All of the slices exhibit METs between the AFM and FM states, but the transitions are abrupt with nearly the same change of magnetization, ΔM, when x = 50.5 and 51, whereas for the x = 50 alloy the transition spreads over a wide temperature interval and ΔM may fluctuate by as much as 10 % from one specimen to another. A comparison of the magnetocaloric responses of x = 50 and 51 materials is presented. The clearly different effect of the magnetic field on the transition in both directions leads to significant differences in the reversibility and maximum values of the magnetic field-induced entropy and adiabatic temperature changes, as well as average hysteresis losses. In terms of reproducibility, our results suggest that induction melting is a more appropriate technique to prepare these binary alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A combined experimental and crystal plasticity study of grain size effects in magnesium alloys

This work presents a method to incorporate the micro Hall-Petch equation into the crystal plasticity finite element (CPFE) framework accounting for the microstructural features to understand the coupling between grain size, texture, and loading direction in magnesium alloys. The effect of grain size and texture is accounted for by modifying the slip resistances of individual basal and prismatic systems based on the micro Hall-Petch equation. The modification based on the micro Hall-Petch equation endows every slip system at each microstructural point with a slip system-level grain size and maximum compatibility factor, which are in turn used to modify the slip resistance. While the slip-system level grain size is a measure of the grain size, the maximum compatibility factor encodes the effect of the grain boundary on the slip system resistance modification and is computed based on the Luster-Morris factor. The model is calibrated using experimental stress-strain curves of Mg-4Al samples with three different grain sizes from which the Hall-Petch coefficients are extracted and compared with Hall-Petch coefficients predicted using original parameters from previous work. The predictability of the model is then evaluated for a Mg-4Al sample with different texture and three grain sizes subjected to loading in different directions. The calibrated parameters are then used for some parametric studies to investigate the variation of Hall-Petch slope for different degrees of simulated spread in basal texture, variation of Hall-Petch slope with loading direction relative to basal poles for a microstructure with strong basal texture, and variation of yield strength with change in grain morphology. The proposed approach to incorporate the micro Hall-Petch equation into the CPFE framework provides a foundation to quantitatively model more complicated scenarios of coupling between grain size, texture and loading direction in the plasticity of Mg alloys.

Crystal plasticity↗
DOE OSTI2025

Supernova pointing capabilities of DUNE

The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electron-neutrino charged-current absorption on 40 Ar and elastic scattering of neutrinos on electrons. Procedures to reconstruct individual interactions, including a newly developed technique called “brems flipping,” as well as the burst direction from an ensemble of interactions are described. Performance of the burst direction reconstruction is evaluated for supernovae happening at a distance of 10 kpc for a specific supernova burst flux model. The pointing resolution is found to be 3.4 degrees at 68% coverage for a perfect interaction-channel classification and a fiducial mass of 40 kton, and 6.6 degrees for a 10 kton fiducial mass respectively. Assuming a 4% rate of charged-current interactions being misidentified as elastic scattering, DUNE’s burst pointing resolution is found to be 4.3 degrees (8.7 degrees) at 68% coverage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Supernova pointing capabilities of DUNE

The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electron-neutrino charged-current absorption on Ar 40 and elastic scattering of neutrinos on electrons. Procedures to reconstruct individual interactions, including a newly developed technique called “brems flipping,” as well as the burst direction from an ensemble of interactions are described. Performance of the burst direction reconstruction is evaluated for supernovae happening at a distance of 10 kpc for a specific supernova burst flux model. The pointing resolution is found to be 3.4 degrees at 68% coverage for a perfect interaction-channel classification and a fiducial mass of 40 kton, and 6.6 degrees for a 10 kton fiducial mass respectively. Assuming a 4% rate of charged-current interactions being misidentified as elastic scattering, DUNE’s burst pointing resolution is found to be 4.3 degrees (8.7 degrees) at 68% coverage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗
DOE OSTI2026

Relativistic collapse of the classical triple bond in the CBi − molecular ion

The conventional framework for chemical bonding between main-group elements involves separate σ and π orbitals to describe multiple bonds. However, relativistic effects mix these orbitals in molecules containing heavy elements through spin–orbit coupling, leaving the total angular-momentum projection (ω) as the only good quantum number. Direct experimental evidence that relativistic effects change the σ-π bonding framework has remained elusive. Here, we probe the carbon-bismuth triple bond in the CBi − anion using high-resolution cryogenic photoelectron spectroscopy, coupled with relativistic four-component Dirac-Coulomb coupled-cluster calculations. Furthermore, even though the CBi − anion is isovalent to the well-known CN − species, we demonstrate that the traditional σ + 2π triple-bond picture collapses into a pure π-like |ω| = 3/2 and two |ω| = 1/2 Kramers pairs containing substantial σ/π mixing.

Kahraman, Deniz [Brown University, Providence, RI ↗
DOE OSTI2020

Enabling fast charging of lithium-ion batteries through secondary-/dual- pore network: Part II - numerical model

To increase the market share of electric vehicles, it is desirable to reduce the battery charge times, which are significantly limited by poor electrolyte transport. A high rate charging is achievable by using expensive and low energy density cells with thin electrodes. For higher energy density cells, new electrolytes with improved conductivity and diffusivity and/or electrodes with advanced architecture are required to boost the electrolyte transport, leading to a more uniform utilization of active materials. In our previous work, an analytical model was developed to investigate the effect of secondary pore network (SPN) on electrolyte transport and the configuration of SPN was optimized by enforcing equal characteristic diffusion times in through-plane and in-plane directions. Here, to evaluate the effect of SPN on the fast-charging capability of lithium-ion batteries, a 2D physics-based electrochemical model is developed with SPN in either one or both electrodes. Additionally, the effect of SPN on cell energy density and lithium plating is investigated for cells with different loadings and electrode porosities. Combining SPN with elevated charging temperatures, the model predicts that the volumetric discharge energy density of a 3 mA h/cm2 cell can reach 270 Wh/L after a 6C constant-current charging.

25 ENERGY STORAGE↗
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↗