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 379 records · Page 21

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗

Generative learning of densities on manifolds

A generative modeling framework is proposed that combines diffusion models and manifold learning to efficiently sample data densities on manifolds. The approach utilizes Diffusion Maps to uncover possible low-dimensional underlying (latent) spaces in the high-dimensional data (ambient) space. Two approaches for sampling from the latent data density are described. The first is a score-based diffusion model, which is trained to map a standard normal distribution to the latent data distribution using a neural network. The second one involves solving an Itô stochastic differential equation in the latent space. Additional realizations of the data are generated by lifting the samples back to the ambient space using Double Diffusion Maps , a recently introduced technique typically employed in studying dynamical system reduction; here the focus lies in sampling densities rather than system dynamics. The proposed approaches enable sampling high dimensional data densities restricted to low-dimensional, a priori unknown manifolds. The efficacy of the proposed framework is demonstrated through a benchmark problem and a material with multiscale structure.

Double diffusion maps↗

A novel physics-based and data-supported microstructure model for part-scale simulation of laser powder bed fusion of Ti-6Al-4V

The elasto-plastic material behavior, material strength and failure modes of metals fabricated by additive manufacturing technologies are significantly determined by the underlying process-specific microstructure evolution. In this work a novel physics-based and data-supported phenomenological microstructure model for Ti-6Al-4V is proposed that is suitable for the part-scale simulation of laser powder bed fusion processes. The model predicts spatially homogenized phase fractions of the most relevant microstructural species, namely the stable β-phase, the stable α s -phase as well as the metastable Martensite α m -phase, in a physically consistent manner. In particular, the modeled microstructure evolution, in form of diffusion-based and non-diffusional transformations, is a pure consequence of energy and mobility competitions among the different species, without the need for heuristic transformation criteria as often applied in existing models. The mathematically consistent formulation of the evolution equations in rate form renders the model suitable for the practically relevant scenario of temperature- or time-dependent diffusion coefficients, arbitrary temperature profiles, and multiple coexisting phases. Due to its physically motivated foundation, the proposed model requires only a minimal number of free parameters, which are determined in an inverse identification process considering a broad experimental data basis in form of time-temperature transformation diagrams. Subsequently, the predictive ability of the model is demonstrated by means of continuous cooling transformation diagrams, showing that experimentally observed characteristics such as critical cooling rates emerge naturally from the proposed microstructure model, instead of being enforced as heuristic transformation criteria. Eventually, the proposed model is exploited to predict the microstructure evolution for a realistic selective laser melting application scenario and for the cooling/quenching process of a Ti-6Al-4V cube of practically relevant size. Numerical results confirm experimental observations that Martensite is the dominating microstructure species in regimes of high cooling rates, e.g., due to highly localized heat sources or in near-surface domains, while a proper manipulation of the temperature field, e.g., by preheating the base-plate in selective laser melting, can suppress the formation of this metastable phase.

Inverse parameter identification↗

Deep learning assisted jet tomography for the study of Mach cones in QGP

Abstract Mach cones are expected to form in the expanding quark-gluon plasma (QGP) when energetic quarks and gluons traverse the hot medium at a velocity faster than the speed of sound in high-energy heavy-ion collisions. The shape of the Mach cone and the associated diffusion wake are sensitive to the initial jet production location and the propagation direction of the parton shower relative to the radial flow because of the distortion caused by the collective expansion of the QGP and the large density gradient. The shape of jet-induced Mach cones and their distortions in heavy-ion collisions provide a unique and direct probe of the dynamical evolution and the equation of state of QGP. However, it is difficult to identify the Mach cone and the diffusion wake in current experimental measurements of final hadron distributions because they are averaged over all possible initial jet production locations and parton-shower propagation directions. To overcome this difficulty, we develop a deep learning assisted jet tomography which uses the full information of the final hadrons from jets to localize the initial jet production positions. This method can help to constrain the initial regions of jet production in heavy-ion collisions and enable a differential study of Mach-cones with different path lengths and orientations relative to the radial flow of the QGP in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗

Investigations into the Ternary NaF-KF-UF4 Salt System – Phase A

A knowledge gap exists in the data and understanding of fresh fuel salt and irradiated multicomponent fuel salt systems thermophysical properties. Quantifying these properties is necessary for the design and construction of test reactors, as well as the licensing of future commercial molten-salt reactors. To facilitate thermal property determination on a proposed fuel salt composition for Seaborg Technologies, several samples containing depleted uranium tetrafluoride (UF4), sodium fluoride (NaF), and potassium fluoride (KF) were blended, and a melt temperature analysis was performed. From the melting temperature analysis, it was determined that sample Seaborg-7, a ternary salt composition of 26.4UF4-24.7KF-48.9NaF (mol%), was very near a ternary eutectic point. Therefore, thermal properties such as melting temperature, salt stability, density, heat capacity, thermal diffusivity, and viscosity were experimentally determined on the Seaborg-7 salt. These measurements document the baseline properties of fresh fuel salt as a function of temperature, where future experiments on irradiated fuel salt will provide a holistic perspective on the change of thermophysical properties during reactor operations. Several precision instruments were used to collect property data, and instrument calibrations and data collection were performed and documented in a standardized and reproducible manner with meticulous detail. This process ensured that the measurement procedures and resulting data can readily be duplicated elsewhere. The Seaborg-7 salt was shown to be stable at temperatures up to 900°C, as no mass change was observed upon repeated heating and cooling. The peak melting temperature was determined to be 547°C (557°C endset). The enthalpy of fusion (??H?_fus^o) was determined to be 167.5 ± 2.7 J/g while the enthalpy of crystallization (??H?_c^o) was determined to be -147.8 ± 13.3 J/g. In addition to the eutectic melting peak, upon heating, several pre eutectic peaks were observed, occurring at 470°C (onset) and 499°C (peak). Specific heat capacity measurements showed a slightly increasing trend with respect to temperature in the solid phase, while the liquid-specific heat capacity showed a somewhat flat trend with an average value of 106.1 ± 1.24 J/mol·K between 600 to 800°C. Three independent trials using the Seaborg-7 salt determined the density to be ?(T) = 4.908 – 0.000363·T(°C), validated between 32 to 200°C, and ?(T) = 4.808 – 0.00113·T(°C), validated between ~575 to 850°C. Thermal diffusivity was determined for the liquid state and is represented by the linear equation y = 0.1581 + 0.000207·T(°C) between 550 to 850°C. The viscosity was determined from 600 to 800°C and is represented by the exponential fit equation, ? (mPa·s) = 736.58e^(-0.006·T(°C)). This report documents the conclusion of fuel salt thermophysical property measurements for the Seaborg SPP, Phase A project.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fully Implicit Monte Carlo (FIMC) (September 2019-September 2020 Progress Report)

This research is focused on applying the NLem method to the full thermal radiative transfer (TRT) equations in the context of a Monte Carlo solution strategy. A framework has been developed for obtaining a solution to the TRT equations using the nonlinear elimination method (NLem) method. Current efforts are focused of the development and verification of the required Monte Carlo algorithm, diffusion algorithm, and material energy solver. Once all three pieces have been verified, the algorithms will then be coupled into a fully functioning method for applying NLem to the TRT equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Risk-Averse Control of Fractional Diffusion with Uncertain Exponent

In this paper, we introduce and analyze a new class of optimal control problems constrained by elliptic equations with uncertain fractional exponents. We utilize risk measures to formulate the resulting optimization problem. We develop a functional analytic framework, study the existence of solution and rigorously derive the first-order optimality conditions. Additionally, we employ a sample-based approximation for the uncertain exponent and the finite element method to discretize in space. Further, we prove the rate of convergence for the optimal risk neutral controls when using quadrature approximation for the uncertain exponent and conclude with illustrative examples.

finite element method↗

Comment on “Hermitian stochastic methodology for x-ray superfluorescence”

This comment critically examines the paper “Hermitian stochastic methodology for x-ray superfluorescence” by Chuchurka, Sukharnikov, and Rohringer [Phys. Rev. A 109, 063705 (2024)]. Here, the paper proposes a Hermitian framework based on the positive-𝑃 stochastic equations (PPSEs), claimed to be derived “from first principles,” through a redefinition of physical variables without a clear rationale. However, because the original PPSEs cannot satisfy Hermiticity under any stochastic-diffusion-gauge transformation, the framework necessarily employs representations beyond this gauge freedom, thereby contradicting the original PPSEs and their ensemble-average predictions. The resulting equations have close structural similarities with the previously established Hermitian formalism of stochastic Maxwell-Bloch equations [arXiv:2305.04653 and Phys. Rev. A 110, 023724 (2024)].

Park, Jeong-Wan [Argonne National Laboratory (ANL)↗

Finite domain solution of a KGD hydraulic fracture in the viscosity-dominated regime

This paper describes a numerical algorithm for solving the classic problem of a plane strain (KGD) fracture propagating in an impermeable elastic medium with zero toughness. The method, which takes advantage of the self-similar nature of the solution, combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation. This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects, noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method. To enhance the efficiency and accuracy of the numerical scheme, the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening, while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation. The proposed technique yields results that closely match the analytical solution, even with a coarse mesh. This approach offers potential for addressing more complex hydraulic fracturing problems in the future.

Domain-based method↗

A model for droplet heating and evaporation of water-in-oil emulsified fuel

Water-in-oil (W/O) emulsified fuel is a promising alternative fuel by inducing flash boiling of water at high temperature that can improve the atomization of fuel spray. The heating and evaporation process of emulsified fuel droplet is affected by the diffusion and coalescence of small dispersed water droplets in oil during the heating. In this study, a model is developed for the droplet heating and evaporation of W/O emulsified fuel with these key physics considered. The diffusivity of dispersed water droplets in oil is calculated using the Stokes-Einstein equation. The deactivation temperature of surfactant is selected as a criterion for the water coalescence since the water coalescence is due to the deactivation of surfactant. The process of water coalescence is simplified that the dispersed water droplets coalesce into one single water sub-droplet at the center of the oil droplet instantly as the droplet temperature reaches the deactivation temperature of surfactant. The model is validated against experimental data of single droplet under different heating temperatures, surfactant concentrations, and sizes of dispersed water droplets. Based on the proposed model, the effects of fuel properties and heating conditions on the droplet heating and evaporation are analyzed.

Droplet coalescence↗

Toward data assimilation of ship-induced aerosol–cloud interactions

Satellite imagery can detect temporary cloud trails or ship tracks formed from aerosols emitted from large ships traversing our oceans, a phenomenon that global climate models cannot directly reproduce. Ship tracks are observable examples of marine cloud brightening, a potential solar climate intervention that shows promise in helping combat climate change. In this paper, we demonstrate a simulation-based approach in learning the behavior of ship tracks based upon a novel stochastic emulation mechanism. Our method uses wind fields to determine the movement of aerosol–cloud tracks and uses a stochastic partial differential equation (SPDE) to model their persistence behavior. This SPDE incorporates both a drift and diffusion term which describes the movement of aerosol particles via wind and their diffusivity through the atmosphere, respectively. We first present our proposed approach with examples using simulated wind fields and ship paths. We then successfully demonstrate our tool by applying the approximate Bayesian computation method-sequential Monte Carlo for data assimilation.

54 ENVIRONMENTAL SCIENCES↗

Moisture-induced autonomous surface potential oscillations for energy harvesting

A variety of autonomous oscillations in nature such as heartbeats and some biochemical reactions have been widely studied and utilized for applications in the fields of bioscience and engineering. Here, we report a unique phenomenon of moisture-induced electrical potential oscillations on polymers, poly([2-(methacryloyloxy)ethyl] dimethyl-(3-sulfopropyl) ammonium hydroxide-co-acrylic acid), during the diffusion of water molecules. Chemical reactions are modeled by kinetic simulations while system dynamic equations and the stability matrix are analyzed to show the chaotic nature of the system which oscillates with hidden attractors to induce the autonomous surface potential oscillation. Using moisture in the ambient environment as the activation source, this self-excited chemoelectrical reaction could have broad influences and usages in surface-reaction based devices and systems. As a proof-of-concept demonstration, an energy harvester is constructed and achieved the continuous energy production for more than 15,000 seconds with an energy density of 16.8 mJ/cm 2 . A 2-Volts output voltage has been produced to power a liquid crystal display toward practical applications with five energy harvesters connected in series.

42 ENGINEERING↗

Generalized impurity pinch in partially magnetized multi-ion plasma

In a two-ion-species plasma with disparate ion masses, heavy ions tend to concentrate in the low-temperature region of collisionally magnetized plasma and in the high-temperature region of collisionally unmagnetized plasma, respectively. Moreover, collisional magnetization can be determined as the ratio of the light ion gyrofrequency to the collision frequency of light and heavy ion species, and the behavior of this effect in the intermediate regime of partially magnetized plasma is predominantly dependent on this Hall parameter. Multi-ion cross-field transport has been described before in the collisionally magnetized plasma regime, and generalized pinch relations, which describe densities of ion species in equilibrium in that plasma, are found in the literature. Here, the role of collisional magnetization and Larmor magnetization in multi-ion collisional transport is clarified, and generalized pinch relations are extended to the partially magnetized regime in which the ion Hall parameter may be small, as long as electrons remain collisionally magnetized. Equilibrium ion density profiles have the same dependence on external forces and on each other regardless of collisional magnetization of ions. The expansion of the range of validity of multi-ion collisional transport models makes them applicable to a wider range of laboratory plasma conditions. In particular, ion density profiles evolve sufficiently fast for radial impurity transport to be observable around stagnation on MagLIF, leading to expulsion of heavy ion impurities from the hotspot as long as plasma becomes sufficiently collisionally magnetized during the implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

jaxhps: An elliptic PDE solver built with machine learning in mind

Elliptic partial differential equations (PDEs) can model many physical phenomena, such as electrostatics, acoustics, wave propagation, and diffusion. In scientific machine learning settings, a high-throughput PDE solver may be required to generate a training dataset, run in the inner loop of an iterative algorithm, or interface directly with a deep neural network. To provide value to machine learning users, such a PDE solver must be compatible with standard automatic differentiation frameworks, scale efficiently when run on graphics processing units (GPUs), and maintain high accuracy for a large range of input parameters. We have designed the jaxhps package with these use-cases in mind by implementing a highly efficient and accurate solver for elliptic problems with native hardware acceleration and automatic differentiation support.

97 MATHEMATICS AND COMPUTING↗

Analysis of Tempered Fractional Operators

Tempered fractional operators provide an improved predictive capability for modeling anomalous effects that cannot be captured by standard partial differential equations. These effects include subdiffusion and superdiffusion (i.e. the mean square displacement in a diffusion process is proportional to a fractional power of the time), that often occur in, e.g., geoscience and hydrology. We analyze tempered fractional operators within the nonlocal vector calculus framework in order to assimilate them to the rigorous mathematical structure developed for nonlocal models. First, we show they are special instances of generalized nonlocal operators by means of a proper choice of the nonlocal kernel. Then, we present a plan for showing tempered fractional operators are equivalent to truncated fractional operators. These truncated operators are useful because they are less computationally intensive than the tempered operators.

97 MATHEMATICS AND COMPUTING↗

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗