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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Optimal design of a model energy conversion device

Fuel cells, batteries, and thermochemical and other energy conversion devices involve the transport of a number of (electro-) chemical species through distinct materials so that they can meet and react at specified multi-material interfaces. Therefore, morphology or arrangement of these different materials can be critical in the performance of an energy conversion device. In this paper, we study a model problem motivated by a solar-driven thermochemical conversion device that splits water into hydrogen and oxygen. We formulate the problem as a system of coupled multi-material reaction-diffusion equations where each species diffuses selectively through a given material and where the reaction occurs at multi-material interfaces. In conclusion, we introduce a phase-field formulation of the optimal design problem and numerically study selected examples.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2019

High-order gradients with the shifted boundary method: An embedded enriched mixed formulation for elliptic PDEs

Here, we propose an extension of the embedded boundary method known as “shifted boundary method” to elliptic diffusion equations in mixed form (e.g., Darcy flow, heat diffusion problems with rough coefficients, etc.). Our aim is to obtain an improved formulation that, for linear finite elements, is at least second-order accurate for both flux and primary variable, when either Dirichlet or Neumann boundary conditions are applied. Following previous work of Nishikawa and Mazaheri in the context of residual distribution methods, we consider the mixed form of the diffusion equation (i.e., with Darcy-type operators), and introduce an enrichment of the primary variable. This enrichment is obtained exploiting the relation between the primary variable and the flux variable, which is explicitly available at nodes in the mixed formulation. The proposed enrichment mimics a formally quadratic pressure approximation, although only nodal unknowns are stored, similar to a linear finite element approximation. We consider both continuous and discontinuous finite element approximations and present two approaches: a non-symmetric enrichment, which, as in the original references, only improves the consistency of the overall method; and a symmetric enrichment, which enables a full error analysis in the classical finite element context. Combined with the shifted boundary method, these two approaches are extended to high-order embedded computations, and enable the approximation of both primary and flux (gradient) variables with second-order accuracy, independently on the type of boundary conditions applied. We also show that the primary variable is third-order accurate, when pure Dirichlet boundary conditions are embedded.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Alluvial fan morphology: A self-similar free boundary problem description

In this work, we examine approximate geometrically self-similar solutions to a parabolic free boundary value problem applied to alluvial fan surface morphology and growth. Alluvial fans are fan- or cone-shaped sedimentary deposits caused by the rapid deposition of sediment from a canyon discharging onto a flatter plain. Longitudinal, topographic profiles of fans can be readily described by a seemingly time independent dimensionless profile (DeChant et al., 1999). However, because an alluvial fan can be expected to grow over time, it is clear that this “steady” profile is certainly time dependent and can be described using a space-time self-similar solution. In an experimental and theory-based study, Guerit et al. (2014) developed a self-similar (or as they describe it a self-affine) linear solution based upon an approximate first order small parameter expansion solution for a 1-d homogeneous nonlinear diffusion equation. Direct substitution of this result into a linear diffusion equation suggests that this first order expression may not fully satisfy the associated governing equation. In contrast, we develop a more complete solution based upon a modeled approximation for the axi-symmetric formulation such that the associated temporal behavior is consistent with a 1/3 time power-law as described by Reitz and Jerolmack (2014). The resulting expression is an exact solution to a linear heat equation. Furthermore, we emphasize that a small parameter is not inherent to the resulting profile result and is not included in our model development. Though developed using rather different approaches, the formal solution developed here is in good agreement with the simple polynomial described by DeChant et al. (1999) suggesting that this self-similar solution is a suitable time dependent representation of alluvial fan longitudinal profile form and improves on earlier work.

58 GEOSCIENCES↗
DOE OSTI2019

PVRD-FASP: A Unified Solver for Modeling Carrier and Defect Transport in Photovoltaic Devices

In our paper, we present a simulator for modeling transport of charge carriers and electrically active defect centers in solar cells by treating them on an equal footing, which allows us to address metastability and reliability issues. The exact nonlinear differential equations set solved by our solver is presented. The formulation of such differential equations, namely the continuity equations, drift-diffusion equation, and Poisson equation, for studying charge and defect transport is explained. The parameters needed for forming the differential equations are taken from first principle calculations. The solver is verified with test cases built on PN heterojunctions, Cu diffusion in single crystal CdTe and comparing Silvaco simulations with our numerical results.

14 SOLAR ENERGY↗
DOE OSTI2018

Initial-Boundary Value Problem for the Nonstationary Radiative Transfer Equation with Diffuse Reflection and Refraction Conditions

We consider the initial-boundary value problem for the nonstationary radiative transfer equation in a system of semitransparent bodies with the conditions of diffuse reflection and refraction of radiation. We establish the unique solvability of the problem with boundary and initial data in the complete scale of Lebesgue spaces. We obtain estimates for solutions.

42 ENGINEERING↗
DOE OSTI2024

Tensor Network Space-Time Spectral Collocation Method for Time-Dependent Convection-Diffusion-Reaction Equations

Emerging tensor network techniques for solutions of partial differential equations (PDEs), known for their ability to break the curse of dimensionality, deliver new mathematical methods for ultra-fast numerical solutions of high-dimensional problems. Here, we introduce a Tensor Train (TT) Chebyshev spectral collocation method, in both space and time, for the solution of the time-dependent convection-diffusion-reaction (CDR) equation with inhomogeneous boundary conditions, in Cartesian geometry. Previous methods for numerical solution of time-dependent PDEs often used finite difference for time, and a spectral scheme for the spatial dimensions, which led to a slow linear convergence. Spectral collocation space-time methods show exponential convergence; however, for realistic problems they need to solve large four-dimensional systems. We overcome this difficulty by using a TT approach, as its complexity only grows linearly with the number of dimensions. We show that our TT space-time Chebyshev spectral collocation method converges exponentially, when the solution of the CDR is smooth, and demonstrate that it leads to a very high compression of linear operators from terabytes to kilobytes in TT-format, and a speedup of tens of thousands of times when compared to a full-grid space-time spectral method. These advantages allow us to obtain the solutions at much higher resolutions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

Ionospheric disturbances and gravity waves

The response of ionization to a gravity wave moving through the ionosphere is studied. Hydrodynamic equations are used, and local thermodynamic equilibrium is imposed for simplicity. The treatment involves a perturbation analysis, and the background medium is assumed to be time stationary, horizontally stratified, and known. It is shown that ionization may be locally resonant at each level for certain frequencies and directions, for which condition neutral and ionized particles are considered closely or critically coupled. The phase direction for this critical coupling is always downward in the absence of a magnetic field. A magnetic field results in two resonant directions for the same frequency, and these directions are mostly downward. Observed TID's associated with gravity waves may be indicative of such resonances. It is also noted that strong coupling may occur to neutral acoustic waves at high altitudes. Previous investigations restrict their use of momentum equations to the diffusion equation. The analysis also shows that such restrictions result in the neglect of terms arising from momentum transport due to any background ambipolar diffusion velocity and wave motion. These terms are mostly relevant at higher altitudes.

Eun, H.↗
DOE OSTI2022

Extension of the PINN diffusion model to k-eigenvalue problems

This paper extends our recent work on the Physics-Informed Neural Networks (PINN) approach for the fixed source diffusion models and applies it to the diffusion theory based k-eigenvalue problems. To make the PINN equitable for the eigenvalue problems, we introduce a novel integral regularization term to the loss function in the framework, and allow the direct inference of the principal eigenvalue and the associated eigenfunction. The regularization term enforces a pre-defined value on the integration of the model predictions, and this value can be directly related to a physical property of the system. We also introduce an additional learnable parameter to approximate the principal eigenvalue. As a proof of principle, we solve the one-group two-dimensional k-eigenvalue neutron diffusion equation in this work. We then provide two numerical examples to demonstrate the applicability of the PINN approach. In each example, we solve the k-eigenvalue diffusion equation in a multi-region configuration constrained with a set of Robin boundary conditions for generality. We use a FEM solution based on the power-iteration method to verify the results of the PINN solution. The results showed relative percentage error in the predicted eigenvalue of about 0.77% and about 1.2% for example 1 and example 2, respectively. The mean absolute error in the predicted flux for example 1 is ∼ 0.002 and for example 2 is ∼ 0.0024. These results indicate some preliminary successes of the PINN application to k-eigenvalue problems. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Optical-property coefficient estimation of bulky medium in experiments with a succinctly analytical calculation

Diffuse optical imaging enables to reconstruct distribution of optical properties, absorption and scattering coefficients, in tissue for breast cancer detection based on diffusion equation with the help of the initial guess obtained from measured data. To estimate the initial guess of the optical-property coefficients, an analytical solution of diffusion equation can be used and compared with the measured data. The analytical solution for a homogeneous infinite medium can be obtained in the frequency domain, expressing that the photon intensity and the phase lag relative to the distance between source and detector in a linear relationship. In this study, a succinct calculation using the trigonometric relation is proposed to estimate the optical-property coefficients. A tank-type and a cylinder Lipovenoes phantoms with two concentrations of 1.25% and 2.5% are employed and measured for verification. It is found that the method proposed here shows better results and results in estimation errors of 0–14.81% for μ{sub a} and 20–42% for μ{sub s}′.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1986

Multi-dimensional initial-boundary value problems with strong nonlinearities

Initial boundary value problems for the general scalar singularly perturbed parabolic equation in a cylindrical domain are studied. The one-dimensional reaction-diffusion equation is considered in the absence of convection. Then, the general one-dimensional and N-dimensional reaction-convection-diffusion equations are analyzed.

Howes, F. A.↗
DOE OSTI2023

Effective dissipation rate in a Liouvillian-graph picture of high-temperature quantum hydrodynamics

At high temperature, generic strongly interacting spin systems are expected to display hydrodynamics: local transport of conserved quantities, governed by classical partial differential equations like the diffusion equation. I argue that the emergence of this dissipative long-wavelength dynamics from the system's unitary microscopic dynamics is controlled by the structure of the Liouvillian graph of the system's Hamiltonian, that is, the graph induced on Pauli strings by commutation with that Hamiltonian. The Liouvillian graph decomposes naturally into subgraphs of Pauli strings of constant diameter, and the coherent dynamics of these subgraphs determines the rate at which operator weight spreads to long operators. In conclusion, this argument provides a quantitative theory of the emergence of a dissipative effective dynamics from unitary microscopic dynamics; it also leads to an effective model with Hilbert space dimension linear in system size and exponential in the UV cutoff for diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Theory of disordered superconductors with applications to nonlinear current response

I present a review of the theory and basic equations for charge transport in superconducting alloys starting from the Keldysh formulation of the quasiclassical transport equations developed by Eilenberger, Larkin and Ovchinnikov, and Eliashberg. This formulation is the natural extension of Landau’s theory of normal Fermi liquids to the superconducting state of strongly correlated metals. For dirty metals the transport equations reduce to equations for charge diffusion, with the current response given by the Drude conductivity at low temperatures. The extension of the diffusion equation for the charge and current response of a strongly disordered normal metal to the superconducting state yields Usadel’s equations for the nonequilibrium quasiclassical Keldysh propagator. The conditions for the applicability of the Usadel equations are discussed, the pair-breaking effect of disorder on the current response, including the nonlinear current response to an electromagnetic field in the dirty limit, τ ≪ ℏ/Δ, are reported. The same nonlinearity is shown to lead to source currents for photon generation and nonlinear Kerr rotation driven by the nonlinear response to excitation of the superconductor by a multi-mode electromagnetic field. The potential relevance of the nonlinear source currents to superconducting radio-frequency cavities as detectors of axion-like dark matter candidates is briefly discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1984

Chaotic ion motion in magnetosonic plasma waves

The motion of test ions in a magnetosonic plasma wave is considered, and the 'stochasticity threshold' of the wave's amplitude for the onset of chaotic motion is estimated. It is shown that for wave amplitudes above the stochasticity threshold, the evolution of an ion distribution can be described by a diffusion equation with a diffusion coefficient D approximately equal to 1/v. Possible applications of this process to ion acceleration in flares and ion beam thermalization are discussed.

Varvoglis, H.↗
NASA NTRS1988

Real-time optical laboratory solution of parabolic differential equations

An optical laboratory matrix-vector processor is used to solve parabolic differential equations (the transient diffusion equation with two space variables and time) by an explicit algorithm. This includes optical matrix-vector nonbase-2 encoded laboratory data, the combination of nonbase-2 and frequency-multiplexed data on such processors, a high-accuracy optical laboratory solution of a partial differential equation, new data partitioning techniques, and a discussion of a multiprocessor optical matrix-vector architecture.

Casasent, David↗
NASA NTRS1994

A Comparison of Metallographic Cooling Rate Methods Used in Meteorites

The primary objective of this study was to test the postulate that cooling rates acquired from metal grains in chondrites are consistent with those from iron meteorites. Both types of metal occur in some Group IAB meteorites, which are mixtures of massive metal with well-developed Widmanstatten structures and chondritic inclusions with dispersed metal grains. The grains have textures and compositions similar to chondritic metal, including negligible P. The meteorites studied show little or no sign of shock reheating and textural evidence indicates that silicates and metal were mixed before Widmanstatten patterns formed during cooling. Cooling rates were obtained by comparing measured to modeled taenite grain or lamellae dimensions and central Ni contents. Modeling entails solving diffusion equations using experimental diffusion coefficients, phase relations, and bulk or local Ni and P contents, taking into account geometry, undercooling, and impingement. There is one set of parameters for grains and another, quite different set for Widmanstatten lamellae, including a factor of 30 difference in diffusion coefficients. Yet cooling rates obtained from Widmanstatten structures and metal grains in chondritic inclusions of the same meteorite are consistent; uncertainties in the best data are +/- 10 K/Ma, equivalent to a factor of 1 +/- 0.25. This agreement implies that the data and models are correct or contain fortuitously offsetting errors, which is quite unlikely. Cooling rates range from 40 K/Ma to 70 K/Ma in IAB meteorites that contain both grains and Widmanstatten structures. Rates based on grains in Ni-poor and Ni-rich meteorites lacking Widmansatten patterns expand the range from 30 K/Ma to perhaps 200 K/Ma. Cooling rates correlate with Ni content; Ni-poor meteorites have slower rates than Ni-rich ones. Evidently, IAB meteorites were radially distributed over greater than 30km in a body with a radius less than 50km. A comparison of the available Ar ages with cooling times inferred from the cooling rates suggests that the parent body cooled more slowly after the metallographic cooling rates were established.

Herpfer, Marc A.↗
NASA NTRS2004

Numerical Simulation of Nanostructure Growth

Nanoscale structures, such as nanowires and carbon nanotubes (CNTs), are often grown in gaseous or plasma environments. Successful growth of these structures is defined by achieving a specified crystallinity or chirality, size or diameter, alignment, etc., which in turn depend on gas mixture ratios. pressure, flow rate, substrate temperature, and other operating conditions. To date, there has not been a rigorous growth model that addresses the specific concerns of crystalline nanowire growth, while demonstrating the correct trends of the processing conditions on growth rates. Most crystal growth models are based on the Burton, Cabrera, and Frank (BCF) method, where adatoms are incorporated into a growing crystal at surface steps or spirals. When the supersaturation of the vapor is high, islands nucleate to form steps, and these steps subsequently spread (grow). The overall bulk growth rate is determined by solving for the evolving motion of the steps. Our approach is to use a phase field model to simulate the growth of finite sized nanowire crystals, linking the free energy equation with the diffusion equation of the adatoms. The phase field method solves for an order parameter that defines the evolving steps in a concentration field. This eliminates the need for explicit front tracking/location, or complicated shadowing routines, both of which can be computationally expensive, particularly in higher dimensions. We will present results demonstrating the effect of process conditions, such as substrate temperature, vapor supersaturation, etc. on the evolving morphologies and overall growth rates of the nanostructures.

Hwang, Helen H.↗
NASA NTRSDate not supplied

Composite Bond Line Measurements Based on A Bayesian Analysis of Flash Thermography Data

For bonded composite materials, an accurate characterization of the adhesive bond line is needed to predict failure modes and fracture toughness. In this paper, bond line thickness was estimated from data obtained using through transmission flash thermography. The forward model that predicts back surface temperature is based on a three layer heat diffusion equation with varying diffusivity and flux boundary conditions. The corresponding inverse problem of estimating bond line thickness from measurement data was solved using a Bayesian approach that assumed Gaussian priors for the bond line thickness and thermal diffusivity of the adherends. Finally, the outputs of the thermography based method were compared to measurements that were collected using a micrometer and ultrasound testing.

Flash thermography↗