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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multigrid Algorithms with Projection and Prolongation over Elements of the Phase Space for K-Eigenvalue Transport Problems

This paper describes new multilevel acceleration methods for solving the multigroup neu- tron transport eigenvalue problems. These multilevel algorithms use different projection and prolongation operators in the phase space. The Nonlinear Diffusion Acceleration (NDA) method with multiple grids in energy is formulated with the prolongation oper- ator based on multiplication iterative correction and linear-in-energy mapping. Another multilevel NDA method uses the projection operator with coarsening in energy between the high-order transport and low-order NDA equations. The third algorithm is formu- lated with the partial-current based CMFD low-order equations and applies projection operators in space and energy. The numerical results are presented.

Cornejo, Luke↗
DOE OSTI2025

Computing Nonequilibrium Responses with Score-Shifted Stochastic Differential Equations

Using equilibrium fluctuations to understand the response of a physical system to an externally imposed perturbation is the basis for linear response theory, which is widely used to interpret experiments and shed light on microscopic dynamics. For nonequilibrium systems, perturbations cannot be interpreted simply by monitoring fluctuations in a conjugate observable and general response results rely on path ensemble averaging. Furthermore, these techniques do not apply to perturbations that affect the diffusion tensor in a stochastic system. Here, we introduce an “effective” physical process that represents the diffusion perturbed dynamics and enables accurate calculations of responses to a change in the diffusion. Interestingly, the effective dynamics contain an additional drift involving the instantaneous “score” of the system, and we leverage score matching algorithms to carry out nonequilibrium response calculations on systems for which the exact stationary distribution is unknown.

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

Learning functional priors and posteriors from data and physics

In this work, we develop a new Bayesian framework based on deep neural networks to be able to extrapolate in space-time using historical data and to quantify uncertainties arising from both noisy and gappy data in physical problems. Specifically, the proposed approach has two stages: (1) prior learning and (2) posterior estimation. At the first stage, we employ the physics-informed Generative Adversarial Networks (PI-GAN) to learn a functional prior either from a prescribed function distribution, e.g., Gaussian process, or from historical data and physics. At the second stage, we employ the Hamiltonian Monte Carlo (HMC) method to estimate the posterior in the latent space of PI-GANs. In addition, we use two different approaches to encode the physics: (1) automatic differentiation, used in the physicsinformed neural networks (PINNs) for scenarios with explicitly known partial differential equations (PDEs), and (2) operator regression using the deep operator network (DeepONet) for PDE-agnostic scenarios. We then test the proposed method for (1) meta-learning for one-dimensional regression, and forward/inverse PDE problems (combined with PINNs); (2) PDE-agnostic physical problems (combined with DeepONet), e.g., fractional diffusion as well as saturated stochastic (100-dimensional) flows in heterogeneous porous media; and (3) spatial-temporal regression problems, i.e., inference of a marine riser displacement field using experimental data from the Norwegian Deepwater Programme (NDP). The results demonstrate that the proposed approach can provide accurate predictions as well as uncertainty quantification given very limited scattered and noisy data, since historical data could be available to provide informative priors. In summary, the proposed method is capable of learning flexible functional priors, e.g., both Gaussian and non-Gaussian process, and can be readily extended to big data problems by enabling mini-batch training using stochastic HMC or normalizing flows since the latent space is generally characterized as low dimensional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{∞}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Parallel-in-Time Integration for Nonlinear Hyperbolic Problems (Final Report)

The work for the subcontract was situated in the area of parallel-in-time integration for hyperbolic partial differential equations (PDEs). Parallel-in-time integration is an active area of research due to its ability to enable faster numerical simulations for applications throughout many areas of science. Over the past two decades, much progress has been made in this area; however, this progress has largely been limited to diffusion-dominated PDEs, with some recent success in scalar linear hyperbolic PDEs. Given the ubiquity of numerical simulations of hyperbolic PDEs throughout the sciences, in particular, nonlinear hyperbolic systems, there is a strong need to develop efficient parallel-in-time techniques for hyperbolic problems beyond simple scalar and linear cases, which is the main topic of this subcontract. The main focus of the work was to further develop and perfect coarse-grid operators for the Multigrid Reduction-inTime (MGRIT) method applied to hyperbolic PDEs that were recently proposed in PhD thesis, based on a modified semi-Lagrangian approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Theory for the Balance between Warm Rain and Ice Crystal Processes of Precipitation in Mixed-Phase Clouds

Abstract Mixed-phase clouds contain both supercooled cloud liquid and ice crystals. In principle, precipitation may be initiated either by the liquid phase or by the ice phase. Ice crystals may grow by vapor diffusion to become snow (“ice crystal process”), forming “cold” precipitation. Equally, cloud droplets, when large enough, coalesce to form “warm” precipitation by the “warm rain process.” Warm rain could be supercooled and freeze as “warm” graupel. In the present paper, a new simplified theoretical analysis is provided to examine the microphysical system consisting of three species of hydrometeor, namely, cloud liquid, “cold ice” (crystals, snow), and “warm rain” (frozen or supercooled). This is obtained by nondimensionalizing and simplifying the evolution equations for the mass of each species. Analytical formulas are given for equilibria. Feedback analysis shows that the sign of the feedback is linked to the abundance of precipitation, with a neutral surface in the 3D phase space. The system’s precipitation amount explodes while in the initial unstable regime, crossing the neutral surface and approaching the equilibrium point that is a stable attractor. Positive and negative feedbacks are elucidated. In a standard case, the cold ice mass is about 1000 times larger than the warm rain mass. To illustrate the physical behavior of the theory, sensitivity tests are performed with respect to environmental conditions (e.g., aerosol, updraft speed) and microphysical parameters (e.g., riming and sedimentation rates for cold ice). Cold ice prevails, especially in fast ascent, due to its low bulk density, favoring slow sedimentation and a wide cross-sectional area for riming. Significance Statement The theory elucidates how the ice phase can prevail in the precipitation from any mixed-phase clouds with supercooled cloud liquid and crystals. The ice phase radically suppresses cloud liquid by riming when active and “wins” the competition against coalescence. This prevalence of ice is shown to arise from the low bulk density of snow. The cloud is viewed as a system of negative and positive feedbacks that prevail in realms of stability and instability in a 3D phase space.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Cost-efficient finite-volume high-order schemes for compressible magnetohydrodynamics

We present an efficient dimension-by-dimension finite-volume method which solves the adiabatic magnetohydrodynamics equations at high discretization order, using the constrained-transport approach on Cartesian grids. Results are presented up to tenth order of accuracy. The algorithmic architecture of this method is very close to that of commonly employed second-order schemes: it requires only one reconstructed value per face for each computational cell, independently of the scheme's order. This property is highly beneficial for the numerical efficiency. It results from reusing the required values already available in neighboring grid cells, in contrast to standard algorithms that require a number of reconstructions and evaluations which increases with the scheme's order of accuracy. At a given resolution, these high-order schemes present significantly less numerical dissipation than commonly employed lower-order approaches. Thus, results of comparable accuracy are achievable at a substantially coarser resolution, yielding overall performance gains. We also present a way to include physical dissipative terms: viscosity, magnetic diffusivity and cooling functions, respecting the finite-volume and constrained-transport frameworks. Benefits of this method are shown through applications in turbulent flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Analysis of ducted fuel injection at high-pressure transcritical conditions using large-eddy simulations

Ducted fuel injection (DFI) is a proposed fuel injection concept for achieving substantial reductions in emissions. In this concept, the fuel is injected through a coannular duct, resulting in increased fuel-air mixing and minimized formation of soot and other unwanted combustion products. Apart from comprehensive experimental investigations on DFI, so far computational studies have been limited to single-point Reynolds-averaged Navier Stokes simulations. Therefore, the objective of this work is to complement these studies by performing large-eddy simulations using a diffuse-interface method to examine the physical mechanisms and combustion processes of DFI, specifically focusing on the mixing process and the effect of fuel-ducting on combustion and pollutant emissions. To this end, finite-rate chemistry simulations are performed of the DFI configuration corresponding to the Engine Combustion Network Spray A injector at transcritical conditions (n-dodecane fuel, 60 bar pressure and 1000 K temperature chamber conditions). A two-equation soot model is employed for the qualitative analysis of soot emissions. Direct comparisons of averaged and instantaneous flow field results with the Spray A configuration are performed to assess the effect of DFI on the first- and second-stage ignition and soot formation. Compared to the free-spray condition, the results show that the DFI case exhibits a combination of (i) increased mass flow rate and entrained air, (ii) larger pressure drop magnitude and flow velocity, and (iii) a closer-to-stoichiometric mixture composition (both globally and locally), each of which is conjectured to contribute toward reduced soot production.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Statistical distributions for transient transport

Here, this paper introduces the use of statistical distributions based on transport differential equations for clear distinction of transport modes within transient kinetic experiments. More specifically, novel techniques are developed for the transient data obtained through the Temporal Analysis of Products (TAP) reactor and are applicable to experiments where pulse response into a gas flow is used. The methodology allows distinguishing between two domains of diffusion transport in heterogeneous catalytic systems, i.e., Knudsen and non-Knudsen diffusion, using statistical fingerprints, and finding the transition domain. Two distribution parameters were obtained that directly result in coefficients that correspond to the concentration and the rate of transport. Using a linear relationship between the rate and concentration coefficients, Knudsen diffusion is revealed when the rate of transport is constant and non-Knudsen diffusion is confirmed when the rate of transport coefficient is a function of the concentration coefficient. As a result, accurate transport information can be extracted from experimental data even in the presence of comprising instrument drift or noise particularly when analyzing higher pressure pulse responses with complex transport. This enables more direct investigation of experiments influenced by gas-phase reactions.

TAP reactor↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

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

Dimensionless parameters for cloudy Rayleigh-Bénard convection: Supersaturation, Damköhler, and Nusselt numbers

In steady-state Rayleigh-Bénard convection, heat is transported by turbulent thermal convection from the bottom, hot surface to the top, cold surface, leading to a height-independent sensible heat flux. When water vapor is present and cloud formation occurs, there is also an additional latent heat flux. Heat transport in cloudy Rayleigh-Bénard convection depends on turbulent flow as well as the microphysical state of the clouds: specifically, whether substantial supersaturations exist and whether cloud liquid water is removed through sedimentation/precipitation. In this article we bridge between the Rayleigh-Bénard convection literature and the atmospheric literature. We express the governing equations for cloudy convection in dimensionless form, thereby explicitly identifying the governing parameters relevant to the cloudy case, including Schmidt, Damköhler, supersaturation, and sedimentation numbers. We further connect to the atmospheric literature by obtaining a Nusselt number (dimensionless heat flux) for a cloud-convection system, directly from the conservation equations for temperature and water vapor. This flux has the same form as that identified by Zhang et al. [L. Zhang, K. L. Chong, and K.-Q. Xia, J. Fluid Mech. 874, 1041 (2019)] for convection with water vapor, but is extended to the cloudy case. For equal thermal and water vapor diffusivities, the flux corresponds to the widely used atmospheric quantities equivalent temperature and moist static energy. Using large eddy simulation (LES) of an idealized cloudy Rayleigh-Bénard convection system with fixed boundary conditions, we find that the equivalent heat flux (Nusselt number) is only weakly dependent on the microphysical details of the system, such as liquid water mixing ratio and cloud droplet number concentration. Finally, from the results, we show the vertical profiles of sensible and latent heat fluxes depend on the liquid water content, whereas the equivalent heat flux remains a constant throughout the height of the chamber.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Mechanochemical topological defects in an active nematic

We propose a reaction-diffusion system that converts topological information of an active nematic into chemical signals. We show that a curvature-activated reaction dipole is sufficient for creating a system that dynamically senses topology by producing a concentration field possessing local extrema coinciding with ±$\frac{1}{2}$ defects. The enabling term is analogous to polarization charge density seen in dielectric materials. We demonstrate the ability of this system to identify defects in both passive and active nematics. Our results illustrate that a relatively simple feedback scheme, expressed as a system of partial differential equations, is capable of producing chemical signals in response to inherently nonlocal structures in anisotropic media. Here, we posit that such coarse-grained systems can help generate testable hypotheses for regulated processes in biological systems, such as morphogenesis, and motivate the creation of bio-inspired materials that utilize dynamic coupling between nematic structure and biochemistry.

42 ENGINEERING↗
DOE OSTI2021

Spatial coupling of gyrokinetic simulations, a generalized scheme based on first-principles

In this paper, we present a scheme that spatially couples two gyrokinetic codes using first-principles. Coupled equations are presented and a necessary and sufficient condition for ensuring accuracy is derived. This new scheme couples both the field and the particle distribution function. The coupling of the distribution function is only performed once every few time-steps, using a five-dimensional (5D) grid to communicate the distribution function between the two codes. This 5D grid interface enables the coupling of different types of codes and models, such as particle and continuum codes, or delta-f and total-f models. Transferring information from the 5D grid to the marker particle weights is achieved using a new resampling technique. Demonstration of the coupling scheme is shown using two XGC gyrokinetic simulations for both the core and edge. We also apply the coupling scheme to two continuum simulations for a one-dimensional advection–diffusion problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Demystifying the Stern layer at a metal–electrolyte interface: Local dielectric constant, specific ion adsorption, and partial charge transfer

Electric double layer (EDL) represents one of the most basic concepts in electrochemistry and is pertinent to diverse engineering applications ranging from electrocatalysis to energy storage. Whereas phenomenological and coarse-grained models have been long established to describe ionic distributions in the diffuse layer, a faithful prediction of the physicochemical properties of the electrode–electrolyte interface from a molecular perspective remains a daunting challenge. In this work, we investigate the charging behavior of an Ag (111) electrode in NaF aqueous solutions leveraging experimental results and theoretical calculations based on the classical density functional theory for ion distributions in the diffuse layer and on the joint density functional theory (JDFT) for the electronic structure. When the Ag electrode is applied with a negative voltage, the surface charge density can be reasonably described by assuming a neutral Stern layer with the dielectric constant dependent on the local electric field as predicted by the Kirkwood equation. However, the specific adsorption of F− ions must be considered when the electrode is positively charged and the fluoride adsorption can be attributed to both physical and chemical interactions. Qualitatively, F− binding and partial charge transfer are supported by JDFT calculations, which predict an increased binding energy as the voltage increases. Our findings shed insight on the molecular characteristics of the Stern layer and the charge behavior of adsorbed species not specified by conventional EDL models.

Wang, Xuepeng (ORCID:0000000285704546)↗
DOE OSTI2023

Superior discretizations and AMG solvers for extremely anisotropic diffusion via hyperbolic operators [Slides]

Diffusion in magnetic confinement fusion is extremely anisotropic in the direction of field lines. Rewrote diffusion system based on directional gradients, apply discretization and solver techniques developed for advection. Orders of magnitude decrease in error and solve wallclock time vs. traditional methods. The next steps include: (1) incorporate into larger MHD simulations, (2) better solvers for closed field lines or mixed regimes, and (3) possibly other extremely anisotropic equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗
DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗