Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equations”

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 1,477 records · Page 82

Regularization by denoising diffusion models for solving inverse PDE problems with application to full waveform inversion

Partial differential equation (PDE)-governed inverse problems are fundamental across various scientific and engineering applications; yet they face significant challenges due to nonlinearity, ill-posedness, and sensitivity to noise. Here, we introduce a computational framework, regularization by denoising using diffusion models for partial differential equations (RED-DiffEq), by integrating physics-driven inversion and data-driven learning. RED-DiffEq leverages pretrained diffusion models as a regularization mechanism for PDE-governed inverse problems. We apply RED-DiffEq to solve the full waveform inversion problem in geophysics, a challenging seismic imaging technique that seeks to reconstruct high-resolution subsurface velocity models from seismic measurement data. Our method shows enhanced accuracy and robustness compared to benchmark methods. Additionally, it exhibits strong generalization and domain decomposition capacity, enabling the inversion of more complex velocity models with larger domains than those used in training the diffusion model. Our framework can also be directly applied to diverse PDE-governed inverse problems.

Shan, Siming [Yale University, New Haven, CT (Unit↗

A Fokker–Planck approach to a moment closure for mixing in variable-density turbulence

We develop a theory for the cascade mixing terms in a moment closure approach to binary active scalar mixing in variable-density turbulence. To address the variable-density complications we apply, as a principle and constraint, the conservation of the probability density function (PDF) through a Fokker–Planck equation with bounded sample space whose attractor is the beta PDF with skewness. Mixing is related to a single-point PDF as a realisability principle to provide mathematically rigorous expressions for the small scale statistics in terms of largescale moments. The problem of the unknown small-scale mixing is replaced with the determination of the drift and diffusion terms of a Fokker–Planck equation in a beta-PDF-convergent stochastic process. We find that realisability of a beta-convergent process requires the mixing time-scale ratio, taken as a constant in passive scalar mixing, to be a function of the mean mass fraction, mean fluid density, the Atwood number, the density-volume correlation and moments of the density field. We develop and compare the new model with direct numerical simulations data of non-stationary homogeneous variable-density turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Challenging Practices of Algebraic Battery Life Models through Statistical Validation and Model Identification via Machine-Learning

Various modeling techniques are used to predict the capacity fade of Li-ion batteries. Algebraic reduced-order models, which are inherently interpretable and computationally fast, are ideal for use in battery controllers, technoeconomic models, and multi-objective optimizations. For Li-ion batteries with graphite anodes, solid-electrolyte-interphase (SEI) growth on the graphite surface dominates fade. This fade is often modeled using physically informed equations, such as square-root of time for predicting solvent-diffusion limited SEI growth, and Arrhenius and Tafel-like equations predicting the temperature and state-of-charge rate dependencies. In some cases, completely empirical relationships are proposed. However, statistical validation is rarely conducted to evaluate model optimality, and only a handful of possible models are usually investigated. This article demonstrates a novel procedure for automatically identifying reduced-order degradation models from millions of algorithmically generated equations via bi-level optimization and symbolic regression. Identified models are statistically validated using cross-validation, sensitivity analysis, and uncertainty quantification via bootstrapping. On a LiFePO 4 /Graphite cell calendar aging data set, automatically identified models utilizing square-root, power law, stretched exponential, and sigmoidal functions result in greater accuracy and lower uncertainty than models identified by human experts, and demonstrate that previously known physical relationships can be empirically "rediscovered" using machine learning.

25 ENERGY STORAGE↗

Experimental study of a three-dimensional, shear-driven, turbulent boundary layer using a three-dimensional laser Doppler velocimeter

Experimental data have been obtained in a three-dimensional, shear-driven, turbulent boundary layer subjected to sudden transverse strain. Measurements made with a newly developed, three-dimensional laser Doppler velocimeter include all components of mean flow, turbulent Reynolds stresses, and triple-product correlations. Eddy viscosities, production, convection, turbulent diffusion, and dissipation (balance of kinetic energy equation) are extracted from the data. It was shown experimentally that eddy viscosity is strongly direction-dependent (not a scalar), indicating a need for full Reynolds stress modeling. Turbulence models for eddy viscosity, stress diffusion, pressure rate-of-strain, and anisotropy of Reynolds stresses are tested using the data. Calculations using a Reynolds stress equation model are compared with the data.

Driver, D. M.↗

A near-wall four-equation turbulence model for compressible boundary layers

A near-wall four-equation turbulence model is developed for the calculation of high-speed compressible turbulent boundary layers. The four equations used are the k-epsilon equations and the theta(exp 2)-epsilon(sub theta) equations. These equations are used to define the turbulent diffusivities for momentum and heat fluxes, thus allowing the assumption of dynamic similarity between momentum and heat transport to be relaxed. The Favre-averaged equations of motion are solved in conjunction with the four transport equations. Calculations are compared with measurements and with another model's predictions where the assumption of the constant turbulent Prandtl number is invoked. Compressible flat plate turbulent boundary layers with both adiabatic and constant temperature wall boundary conditions are considered. Results for the range of low Mach numbers and temperature ratios investigated are essentially the same as those obtained using an identical near-wall k-epsilon model. In general, the numerical predictions are in very good agreement with measurements and there are significant improvements in the predictions of mean flow properties at high Mach numbers.

Sommer, T. P.↗

Near-wall variable-Prandtl-number turbulence model for compressible flows

A near-wall four-equation turbulence model is developed for the calculation of high-speed compressible turbulent boundary layers. The four equations used are the k-epsilon equations and the theta(exp 2)-epsilon (sub theta) equations. These equations are used to define the turbulent diffusivities for momentum and heat fluxes, thus allowing the assumption of dynamic similarity between momentum and heat transport to be relaxed. The Favre-averaged equations of motion are solved in conjunction with the four transport equations. Calculations are compared with measurements and with another model's predictions where the assumption of the constant turbulent Prandtl number is invoked. Compressible flat plate turbulent boundary layers with both adiabatic and constant temperature wall boundary conditions are considered. Results for the range of low Mach numbers and temperature ratios investigated are essentially the same as those obtained using an identical near-wall k-epsilon model. In general, there are significant improvements in the predictions of mean flow properties at high Mach numbers.

Sommer, T. P.↗

Steady-state thermal-solutal diffusion in a float zone

A model for a float zone in a thin vertical sheet is used to study the formation of the zone configuration for both pure and binary systems of Si and GeGa. Equations describing the steady two-dimensional diffusion of solute in the melt and diffusion of heat in the feed material and the product crystal are presented. The material properties, asymptotic solutions for the temperature and concentration profiles, and melting, solidifying, and melt/gas interfacial shapes are determined in the small aspect ratio limit.

Young, G. W.↗

Reacting viscous-shock-layer solutions with multicomponent diffusion and mass injection

Numerical solutions of the viscous-shock-layer equation where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium are presented. Also the effects of the diffusion model, surface catalysis, and mass injection on surface transport and flow parameters are considered. The flow is treated as a mixture of five inert and thermally perfect species. The viscous-shock-layer equations are solved by using an implicit-difference scheme. All calculations are for hyperboloids with included angles of 20 and 45. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented to show the effects of the chemical models; diffusion models; surface catalysis; and mass injection of air on heat transfer; skin friction; shock standoff distance; wall pressure distribution; and tangential velocity, temperature, and species profiles. The results show that an equilibrium analysis can substantially overpredict the heat-transfer rates for flow conditions experienced by earth-orbital entry vehicles. Moreover, at such conditions surface catalysis significantly influences heat-transfer and flow-field properties. If a binary rather than a multicomponent diffusion model is assumed, negligible errors in most flow properties result. Quantitative results are presented that show the effect of mass injection on flow properties within and downstream of the injection region.

Moss, J. N.↗

Analysis of Multipoint-Multitime Correlations and Diffusion in Decaying Homogeneous Turbulence

Two-point, two-time correlation equations are obtained by considering the Navier-Stokes equations for two points in a fluid at two time. By neglecting the triple correlations in the equations, a solution is obtained for the final period of decay. The analysis is extended to earlier times by considering three points at three different times. The set of equations is made determinate by neglecting the quadruple correlations in comparison with the triple correlations. The diffusion of particles from a source in a decaying turbulent field is calculated approximately by assuming that the velocity fluctuations are small.

Deissler, Robert G.↗

Nonstationary Problem of Complex Heat Transfer in a System of Semitransparent Bodies with Boundary-Value Conditions of Diffuse Reflection and Refraction of Radiation

We consider a nonstationary initial boundary-value problem describing complex (radiative-conductive) heat transfer in a system of semitransparent bodies. To describe radiation propagation, we use the radiation transfer equation with boundary-value conditions of diffuse reflection and refraction of radiation. We take into account that the radiation intensity and optical properties of bodies depend on the radiation frequency. The unique solvability of a weak solution is established. The comparison theorem is proved. A priori estimates of a weak solution as well as its regularity are obtained.

97 MATHEMATICS AND COMPUTING↗

WE-AB-204-07: Spatiotemporal Distribution of the FDG PET Tracer in Solid Tumors: Contributions of Diffusion and Convection Mechanisms

Purpose: In this study, a mathematical model is utilized to simulate FDG distribution in tumor tissue. In contrast to conventional compartmental modeling, tracer distributions across space and time are directly linked together (i.e. moving beyond ordinary differential equations (ODEs) to utilizing partial differential equations (PDEs) coupling space and time). The diffusion and convection transport mechanisms are both incorporated to model tracer distribution. We aimed to investigate the contributions of these two mechanisms on FDG distribution for various tumor geometries obtained from PET/CT images. Methods: FDG transport was simulated via a spatiotemporal distribution model (SDM). The model is based on a 5K compartmental model. We model the fact that tracer concentration in the second compartment (extracellular space) is modulated via convection and diffusion. Data from n=45 patients with pancreatic tumors as imaged using clinical FDG PET/CT imaging were analyzed, and geometrical information from the tumors including size, shape, and aspect ratios were classified. Tumors with varying shapes and sizes were assessed in order to investigate the effects of convection and diffusion mechanisms on FDG transport. Numerical methods simulating interstitial flow and solute transport in tissue were utilized. Results: We have shown the convection mechanism to depend on the shape and size of tumors whereas diffusion mechanism is seen to exhibit low dependency on shape and size. Results show that concentration distribution of FDG is relatively similar for the considered tumors; and that the diffusion mechanism of FDG transport significantly dominates the convection mechanism. The Peclet number which shows the ratio of convection to diffusion rates was shown to be of the order of 10−{sup 3} for all considered tumors. Conclusion: We have demonstrated that even though convection leads to varying tracer distribution profiles depending on tumor shape and size, the domination of the diffusion phenomenon prevents these factors from modulating FDG distribution.

60 APPLIED LIFE SCIENCES↗

Multi-resolution partial differential equations preserved learning framework for spatiotemporal dynamics

Traditional data-driven deep learning models often struggle with high training costs, error accumulation, and poor generalizability in complex physical processes. Physics-informed deep learning (PiDL) addresses these challenges by incorporating physical principles into the model. Most PiDL approaches regularize training by embedding governing equations into the loss function, yet this depends heavily on extensive hyperparameter tuning to weigh each loss term. To this end, we propose to leverage physics prior knowledge by “baking” the discretized governing equations into the neural network architecture via the connection between the partial differential equations (PDE) operators and network structures, resulting in a PDE-preserved neural network (PPNN). This method, embedding discretized PDEs through convolutional residual networks in a multi-resolution setting, largely improves the generalizability and long-term prediction accuracy, outperforming conventional black-box models. The effectiveness and merit of the proposed methods have been demonstrated across various spatiotemporal dynamical systems governed by spatiotemporal PDEs, including reaction-diffusion, Burgers’, and Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗

Effect of eddy diffusivity on wind-driven currents in a two-layer stratified lake

The steady state wind-driven circulation was numerically calculated in a rectangular stratified lake. The lake is composed of two layers having uniform but unequal densities and eddy diffusivities. The position in thermocline and the three-dimensional velocities in both layers calculated using shallow lake equations. The results show that, as the eddy diffusivity in the hypolimnion is increased, the thermocline tilt and hypolimnetic velocities increase. The effect of the other variables such as wind stress, density, basin length, and mean thermocline depth are also shown.

Gedney, R. T.↗

A stochastic diffusion process for Lochner's generalized Dirichlet distribution

The method of potential solutions of Fokker-Planck equations is used to develop a transport equation for the joint probability of N stochastic variables with Lochner’s generalized Dirichlet distribution as its asymptotic solution. Individual samples of a discrete ensemble, obtained from the system of stochastic differential equations, equivalent to the Fokker-Planck equation developed here, satisfy a unit-sum constraint at all times and ensure a bounded sample space, similarly to the process developed in for the Dirichlet distribution. Consequently, the generalized Dirichlet diffusion process may be used to represent realizations of a fluctuating ensemble of N variables subject to a conservation principle. Compared to the Dirichlet distribution and process, the additional parameters of the generalized Dirichlet distribution allow a more general class of physical processes to be modeled with a more general covariance matrix.

97 MATHEMATICS AND COMPUTING↗

Two critical issues in Langevin simulation of gas flows

A stochastic algorithm based on the Langevin equation has been recently proposed to simulate rarefied gas flows. Compared with the direct simulation Monte Carlo (DSMC) method, the Langevin method is more efficient in simulating small Knudsen number flows. While it is well-known that the cell sizes and time steps should be smaller than the mean free path and the mean collision time, respectively, in DSMC simulations, the Langevin equation uses a drift term and a diffusion term to describe molecule movements, so no direct molecular collisions have to be modeled. This enables the Langevin simulation to proceed with a much larger time step than that in the DSMC method. Two critical issues in Langevin simulation are addressed in this paper. The first issue is how to reproduce the transport properties as that described by kinetic theory. Transport coefficients predicted by Langevin equation are obtained by using Green-Kubo formulae. The second issue is numerical scheme with boundary conditions. We present two schemes corresponding to small time step and large time step, respectively. For small time step, the scheme is similar to DSMC method as the update of positions and velocities are uncoupled; for large time step, we present an analytical solution of the hitting time, which is the crucial factor for accurate simulation. Velocity-Couette flow, thermal-Couette flow, Rayleigh-Bénard flow and wall-confined problem are simulated by using these two schemes. Our study shows that Langevin simulation is a promising tool to investigate small Knudsen number flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magneto-hydrodynamics simulation study of deflagration mode in co-axial plasma accelerators

Experimental studies by Poehlmann et al. [Phys. Plasmas 17(12), 123508 (2010)] on a coaxial electrode magnetohydrodynamic (MHD) plasma accelerator have revealed two modes of operation. A deflagration or stationary mode is observed for lower power settings, while higher input power leads to a detonation or snowplow mode. A numerical modeling study of a coaxial plasma accelerator using the non-ideal MHD equations is presented. The effect of plasma conductivity on the axial distribution of radial current is studied and found to agree well with experiments. Lower conductivities lead to the formation of a high current density, stationary region close to the inlet/breech, which is a characteristic of the deflagration mode, while a propagating current sheet like feature is observed at higher conductivities, similar to the detonation mode. Results confirm that plasma resistivity, which determines magnetic field diffusion effects, is fundamentally responsible for the two modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large deviations in stochastic heat-conduction processes provide a gradient-flow structure for heat conduction

We consider three one-dimensional continuous-time Markov processes on a lattice, each of which models the conduction of heat: the family of Brownian Energy Processes with parameter m (BEP(m)), a Generalized Brownian Energy Process, and the Kipnis-Marchioro-Presutti (KMP) process. The hydrodynamic limit of each of these three processes is a parabolic equation, the linear heat equation in the case of the BEP(m) and the KMP, and a nonlinear heat equation for the Generalized Brownian Energy Process with parameter a (GBEP(a)). We prove the hydrodynamic limit rigorously for the BEP(m), and give a formal derivation for the GBEP(a). We then formally derive the pathwise large-deviation rate functional for the empirical measure of the three processes. These rate functionals imply gradient-flow structures for the limiting linear and nonlinear heat equations. We contrast these gradient-flow structures with those for processes describing the diffusion of mass, most importantly the class of Wasserstein gradient-flow systems. The linear and nonlinear heat-equation gradient-flow structures are each driven by entropy terms of the form -log ρ; they involve dissipation or mobility terms of order ρ² for the linear heat equation, and a nonlinear function of ρ for the nonlinear heat equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The solution of the laminar-boundary-layer equation for the flat plate for velocity and temperature fields for variable physical properties and for the diffusion field at high concentration

In connection with Pohlhausen's solution for the temperature field on the flat plate, a series of formulas were indicated by means of which the velocity and temperature field for variable physical characteristics can be computed by an integral equation and an iteration method based on it. With it, the following cases were solved: On the assumption that the viscosity simply varies with the temperature while the other fluid properties remain constant, the velocity and temperature field on the heated and cooled plate, respectively, was computed at the Prandtl numbers 12.5 and 100 (viscous fluids). A closer study of these two cases resulted in general relations: The calculations for a gas of Pr number 0.7 (air) were conducted on the assumption that all fluid properties vary with the temperature, and the velocities are low enough for the heat of friction to be discounted. The result was a thickening of the boundary layers, but no appreciable modification in shearing stress or heat-transfer coefficient.

BOUNDARY LAYER, LAMINAR - PLATES, FLAT↗