Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Rotational memory function of SPC/E water

Memory effects are essential for the dynamics of condensed materials and are responsible for non-exponential relaxation of correlation functions of dynamic variables through the memory function. Memory functions of dipole rotations for water have never been calculated directly from molecular dynamics simulations. We present here calculations of memory functions for single-dipole rotations and for the overall dipole moment of the sample for SPC/E water. The normalized memory functions for single-particle and collective dipole dynamics turn out to be nearly identical. This result validates theories of dielectric spectroscopy in terms of single-particle time correlation functions and the connection between the collective and single-particle relaxation times through the Kirkwood factor. The dielectric function in this formalism contains no new dynamic information that does not exist in the single-dipole correlation function. A short memory time, ≲1 fs, justifies the use of the mathematics of rotational diffusion to describe the dynamics of a single molecular dipole moment in bulk water. Here, an analytical equation for the rotational memory time is derived.

Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗
DOE OSTI2024

Relativistic diffusion model for hadron production in p -Pb collisions at energies available at the CERN Large Hadron Collider

We investigate charged-hadron production in relativistic heavy-ion collisions of asymmetric systems within a nonequilibrium-statistical framework. Calculated centrality-dependent pseudorapidity distributions for p-Pb collisions at s NN = 5.02 and 8.16 TeV are compared with data from the Large Hadron Collider. Our approach combines a relativistic diffusion model with formulations based on quantum chromodynamics while utilizing numerical solutions of a Fokker-Planck equation to account for the shift and broadening of the fragmentation sources for particle-production with respect to the stopping (net-baryon) rapidity distributions. To represent the centrality dependence of charged-hadron production in asymmetric systems over a broad region of pseudorapidities, the consideration and precise modeling of the fragmentation sources, along with the central gluon-gluon source, is found to be essential. Specifically, this results in an inversion of the particle-production amplitude from backward to forward dominance when transitioning from central to peripheral collisions, in agreement with recent ATLAS and ALICE p-Pb data at s NN = 5.02 T e V . Published by the American Physical Society 2024

Schulz, Philipp (ORCID:0009000266642148)↗
DOE OSTI2023

Final-stage densification kinetics of direct current–sinteredZrB2

Final-stage sintering was analyzed for nominally phase pure zirconium diboride synthesized by borothermal reduction of high-purity ZrO2. Analysis was conducted on ZrB2 ceramics with relative densities greater than 90% using the Nabarro–Herring stress–directed vacancy diffusion model. Temperatures of 1900°C or above and an applied uniaxial pressure of 50 MPa were required to fully densify ZrB2 ceramics by direct current sintering. Ram travel data were collected and used to determine the relative density of the specimens during sintering. Specimens sintered between 1900 and 2100°C achieved relative densities greater than 97%, whereas specimens sintered below 1900°C failed to reach the final stage of sintering. The average grain size ranged from 1.0 to 14.7 μm. The activation energy was calculated from the slope of an Arrhenius plot that used the Kalish equation. The activation energy was 162 ± 34 kJ/mol, which is consistent with the activation energy for dislocation movement in ZrB2. The diffusion coefficients for dislocation motion that controls densification were 5.1 × 10−6 cm2/s at 1900°C and 5.1 × 10−5 cm2/s at 2100°C, as calculated from activation energy and average grain sizes. This study provides evidence that the dominant mechanism for final-stage sintering of ZrB2 ceramics is dislocation motion.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hydrodynamics of higher-rank gauge theories

We extend recent work on hydrodynamics with global multipolar symmetries — known as “fracton hydrodynamics” — to systems in which the multipolar symmetries are gauged. We refer to the latter as “fracton magnetohydrodynamics”, in analogy to conventional magnetohydrodynamics (MHD), which governs systems with gauged charge conservation. We show that fracton MHD arises naturally from higher-rank Maxwell’s equations and in systems with one-form symmetries obeying certain constraints; while we focus on “minimal” higher-rank generalizations of MHD that realize diffusion, our methods may also be used to identify other, more exotic hydrodynamic theories (e.g., with magnetic subdiffusion). In contrast to semi-microscopic derivations of MHD, our approach elucidates the origin of the hydrodynamic modes by identifying the corresponding higher-form symmetries. Being rooted in symmetries, the hydrodynamic modes may persist even when the semi-microscopic equations no longer provide an accurate description of the system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Microstructural Coarsening Kinetics and Mechanical Property Changes in Long-Term Aged Sn–Pb–Sb Solder Joints

Tin-lead-antimony (50Sn–47Pb–3Sb wt.%) soldered assemblies were mechanically tested approximately 30 years after initial production and found to have solder joints of reduced strength. The microstructure of this solder alloy exhibits a ternary eutectic structure with Sn-rich, Pb-rich, and SnSb phases. Accelerated aging was performed to evaluate solder microstructural coarsening and associated strength of laboratory solder joints to correlate these properties to the “naturally aged” solder joints. Isothermal aging was conducted at room temperature, 55, 70, 100, and 135 °C and aging times that ranged from 0.1 to 365 days. The coarsening kinetics of the Pb-rich phase were determined through optical microscopy and image analysis methods established in previous studies on binary Sn–Pb solder. A kinetic equation was developed with time exponent n of 0.43 and activation energy of 24000 J/mol, suggesting grain boundary diffusion or other fast diffusion pathways controlling the microstructural evolution. Compression testing and Vickers microhardness showed significant strength loss within the first 20–30 days after soldering; then, the microstructure and mechanical properties changed more slowly over long periods of time. Further, by combining accelerated aging data and the microstructure-based kinetics, strength predictions were made that match well with the properties of the actual soldered assemblies naturally aged for 30 years. However, aging at the highest temperature of 135 °C produced anomalous behavior suggesting that extraneous aging mechanisms are active. Therefore, data obtained at this temperature or higher should not be used. Overall, the combined microstructural and mechanical property methods used in this study confirmed that the observed reduction in strength of ~ 30-year-old solder joints can be accounted for by the microstructural coarsening that takes place during long-term solid-state aging.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Learning effective SDEs from Brownian dynamic simulations of colloidal particles

We construct a reduced, data-driven, parameter dependent effective stochastic differential equation (eSDE) for electric-field mediated colloidal crystallization using data obtained from Brownian dynamics simulations. We use diffusion maps (a manifold learning algorithm) to identify a set of useful latent observables. In this latent space we identify an eSDE using a deep learning architecture inspired by numerical stochastic integrators and compare it with the traditional Kramers–Moyal expansion estimation. We show that the obtained variables and the learned dynamics accurately encode the physics of the Brownian dynamic simulations. We further illustrate that our reduced model captures the dynamics of corresponding experimental data. Further, our dimension reduction/reduced model identification approach can be easily ported to a broad class of particle systems dynamics experiments/models.

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

Nonintrusive manufactured solutions for ablation

The study of heat transfer and ablation plays an important role in many problems of scientific and engineering interest. As with the computational simulation of any physical phenomenon, the first step toward establishing credibility in ablation simulations involves code verification. Code verification is typically performed using exact and manufactured solutions. However, manufactured solutions generally require the invasive introduction of an artificial forcing term within the source code such that the code solves a modified problem for which the solution is known. In this paper, we present a nonintrusive method for manufacturing solutions for a non-decomposing ablation code, which does not require the addition of a source term.

42 ENGINEERING↗
DOE OSTI2021

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

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

Decoherence and Brownian motion of a polarizable particle near a medium

Optically levitated nanoparticles are ideal experimental testbeds for investigating macroscopic superpositions and microscopic thermodynamics. Integrating such levitated nanoparticles with photonic structures can enable strong coupling between their center-of-mass motion and guided photonic modes, facilitating enhanced control and probing of their motion. When coupling a particle to a photonic structure, such as a waveguide, the effects of fluctuations become prominent at nanoscales. In this work, we analyze the classical and quantized center-of-mass motion of a polarizable particle interacting with the fluctuations of the electromagnetic field in the presence of a medium. We derive a position localization master equation for the particle's quantized center of mass, and examine its classical center-of-mass momentum diffusion, elucidating correspondences between classical and quantum Brownian motion of polarizable particles near media. We study the decoherence rate of the particle in the presence of a planar surface as a function of temperature and distance from the surface, comparing it to common sources of decoherence. Our results are pertinent to experiments aimed at preparing levitated nanospheres in macroscopic quantum states and investigating their Brownian dynamics.

Brownian motion↗
DOE OSTI2020

Influence of Irradiation-induced Microstructural Defects on Thermal Transport in Single Crystal Thorium Dioxide

Microstructural defects formed as a result of fission fragment damage are known to drastically alter thermal conductivity, a fuel property that governs the efficiency of a nuclear reactor. A fundamental understanding of radiation-induced effects on thermal transport is critical for the development of advanced fuels. Ion-irradiation has widely been used to simulate the effects of neutron-irradiation by seeding atomic-to-nanoscale defects of controllable size and density. In this work, we investigate the influence of microstructural defects, induced by ion-irradiation, on the thermal properties of single crystal thorium dioxide samples grown using the hydrothermal technique. The samples were irradiated using 2 MeV protons at room temperature. A laser-based modulated thermoreflectance technique was used to measure the thermal diffusivity and conductivity within the damaged ThO2 region. The experimentally measured thermal properties were compared with Boltzmann Transport Equation predictions, to determine the impact of various defect types on the thermal properties of ThO2.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Comparison of Some RANS Solvers

We will take a look at solving the Reynolds-averaged Navier-Stokes (RANS) equations that are encountered in the context of wind farm performance simulations and optimizations. We will compare some of the more popular ways to solve these equations with a focus on using iterative solvers for the linear solve. We will compare their performance and reliability to a direct solve as we scale the problem both by adding more and parallel resources and by increasing the size of the domain, both in two and three dimensions. There are many strategies that can be applied to solving the RANS equations, some are very efficient, while others are very insensitive to, for example, the Reynolds number. The first contender we will consider is the Pressure-Convection-Diffusion (PCD) preconditioner. Early results suggest that PCD is indeed a very efficient solver, in particular in two dimensions, as long as the Reynold's number remains small. Next we will try to reorder our degrees of freedom such that we can use GMRES with ILU for our linear solve. Another popular choice we will consider for solving the RANS equations is SIMPLE (and its derivatives). For all of our implementations we make use of either FEniCS or Firedrake, basing our work on both existing implementations of some of these solvers while also writing new extensions for others.

CFD↗
DOE OSTI2021

Poisson_CCD: A dedicated simulator for modeling CCDs

A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson’s equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs that are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel–pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A FETI approach to domain decomposition for meshfree discretizations of nonlocal problems

We propose a domain decomposition method for the efficient simulation of nonlocal problems. Our approach is based on a multi-domain formulation of a nonlocal diffusion problem where the subdomains share “nonlocal” interfaces of the size of the nonlocal horizon. This system of nonlocal equations is first rewritten in terms of minimization of a nonlocal energy, then discretized with a meshfree approximation and finally solved via a Lagrange multiplier approach in a way that resembles the finite element tearing and interconnect method. Specifically, we propose a distributed projected gradient algorithm for the solution of the Lagrange multiplier system, whose unknowns determine the nonlocal interface conditions between subdomains. Several two-dimensional numerical tests on problems as large as 191 million unknowns illustrate the strong and the weak scalability of our algorithm, which outperforms the standard approach to the distributed numerical solution of the problem. Finally, this work is the first rigorous numerical study in a two-dimensional multi-domain setting for nonlocal operators with finite horizon and, as such, it is a fundamental step towards increasing the use of nonlocal models in large scale simulations.

42 ENGINEERING↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2020

Fundamental gap of fluorographene by many-body GW and fixed-node diffusion Monte Carlo methods

Fluorographene (FG) is a promising graphene-derived material with a large bandgap. Currently existing predictions of its fundamental gap (Δ f ) and optical gap (Δ opt ) significantly vary when compared with experiment. We provide here an ultimate benchmark of Δ f for FG by many-body GW and fixed-node diffusion Monte Carlo (FNDMC) methods. Both approaches independently arrive at Δf ≈ 7.1 ± 0.1 eV. In addition, the Bethe–Salpeter equation enabled us to determine the first exciton binding energy, Eb = 1.92 eV. We also point to the possible misinterpretation problem of the results obtained for gaps of solids by FNDMC with single-reference trial wave functions of Bloch orbitals. We argue why instead of Δ opt , in the thermodynamic limit,

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

Deuterium and tritium anomalous transport in the generalized Hasegawa–Wakatani resistive drift wave turbulence model with finite ion Larmor radius

Anomalous transport of mixed deuterium–tritium plasma in the edge of magnetic fusion reactors is investigated using numerical solutions of resistive drift wave turbulence model equations, including finite Larmor radius effects, that are derived within the generalized Hasegawa–Wakatani framework. The anomalous cross field diffusivities of deuterium and tritium are compared in turbulence regimes with different values of the electron adiabaticity parameter controlling the existence of zonal flow. The dependence of the tritium-to-deuterium diffusivity ratio on the deuterium and tritium densities and the logarithmic density gradients is analyzed, and a scaling relation is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Reproducibility of fixed-node diffusion Monte Carlo across diverse community codes: The case of water–methane dimer

Fixed-node diffusion quantum Monte Carlo (FN-DMC) is a widely trusted many-body method for solving the Schrödinger equation, known for its reliable predictions of material and molecular properties. Furthermore, its excellent scalability with system complexity and near-perfect utilization of computational power make FN-DMC ideally positioned to leverage new advances in computing to address increasingly complex scientific problems. Even though the method is widely used as a computational gold standard, reproducibility across the numerous FN-DMC code implementations has yet to be demonstrated. This difficulty stems from the diverse array of DMC algorithms and trial wave functions, compounded by the method’s inherent stochastic nature. Here, this study represents a community-wide effort to assess the reproducibility of the method, affirming that yes, FN-DMC is reproducible (when handled with care). Using the water–methane dimer as the canonical test case, we compare results from eleven different FN-DMC codes and show that the approximations to treat the non-locality of pseudopotentials are the primary source of the discrepancies between them. In particular, we demonstrate that, for the same choice of determinantal component in the trial wave function, reliable and reproducible predictions can be achieved by employing the T-move, the determinant locality approximation, or the determinant T-move schemes, while the older locality approximation leads to considerable variability in results. These findings demonstrate that, with appropriate choices of algorithmic details, fixed-node DMC is reproducible across diverse community codes—highlighting the maturity and robustness of the method as a tool for open and reliable computational science.

Della Pia, Flaviano [Univ. of Cambridge (United Ki↗