Search NASA⌕ Search

SEARCH · Search NASA

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Diffusion of cold magnetospheric ions

The paper investigates conditions under which a diffusion equation is a valid description of the transport across magnetic field lines of low energy magnetospheric plasma. The case in which the guiding center drift and its fluctuating component are comparable in size to the particle thermal speed is considered. It is concluded that a two-dimensional diffusion equation is valid in the alpha-beta space of the Euler potentials, provided that the fluctuations have a short enough correlation time so that the guiding center drifts remain negligible, and ions in Jupiter's Io plasma torus can thus be described by a diffusion equation if the correlation time is no more than a few rotation periods of the planet.

Birmingham, T. J.↗
NASA NTRS1993

Multilevel methods for transport equations in diffusive regimes

We consider the numerical solution of the single-group, steady state, isotropic transport equation. An analysis by means of the moment equations shows that a discrete ordinate S(sub N) discretization in direction (angle) with a least squares finite element discretization in space does not behave properly in the diffusion limit. A scaling of the S(sub N) equations is introduced so that the least squares discretization has the correct diffusion limit. For the resulting discrete system a full multigrid algorithm was developed.

Manteuffel, Thomas A.↗
DOE OSTI2025

Space-Time Finite Element Tensor Network Approach for the Time-Dependent Convection–Diffusion–Reaction Equation with Variable Coefficients

In this paper, we present a new space-time Galerkin-like method, where we treat the discretization of spatial and temporal domains simultaneously. This method utilizes a mixed formulation of the tensor-train (TT) and quantized tensor-train (QTT) (please see Section Tensor-Train Decomposition), designed for the finite element discretization (Q1-FEM) of the time-dependent convection–diffusion–reaction (CDR) equation. We reformulate the assembly process of the finite element discretized CDR to enhance its compatibility with tensor operations and introduce a low-rank tensor structure for the finite element operators. Recognizing the banded structure inherent in the finite element framework’s discrete operators, we further exploit the QTT format of the CDR to achieve greater speed and compression. Additionally, we present a comprehensive approach for integrating variable coefficients of CDR into the global discrete operators within the TT/QTT framework. The effectiveness of the proposed method, in terms of memory efficiency and computational complexity, is demonstrated through a series of numerical experiments, including a semi-linear example.

convection–diffusion–reaction equation↗
DOE OSTI2023

1-D coupled surface flow and transport equations revisited via the physics-informed neural network approach

The de Saint-Venant equation (SVE) and advection–diffusion equation (ADE) are commonly employed to solve solute transport problems in surface water. In this work, we propose a mesh-free method based on the physics-informed neural network (PINN) to solve the one dimensional (1-D) SVE, ADE, and the coupled SVE and ADE (SVE-ADE) under various initial and boundary conditions. The PINN model extends the architecture of deep neural networks (DNNs) with implementation of loss function, which are additionally subject to constraints imposed by the physical laws of SVE and ADE, along with their initial and boundary conditions. In such a manner, PINNs can be quickly steered to the true solution while obeying the physical laws. The results of PINN model are compared with the analytical and/or numerical solutions under various conditions to investigate its accuracy and efficiency in solving the SVE, ADE, and SVE-ADE. Our results indicate PINN can accurately simulate the shock wave morphology and avoid numerical dissipation in unsteady flow condition. The PINN method outweighed traditional numerical methods in several aspects, including its ability to function with small amounts of data, no grid discretization, and random selection of sampling points, etc. Additionally, the PINN method is also suitable for solving inverse problems with sparse and noisy data. With 1% noise and 2000 initial and boundary condition points (N u ), the errors of the estimated flow rate (v) and diffusion coefficient (D) are 0.003% and 0.105%, respectively, which indicate the accuracy and robustness of the proposed method. Finally, our results indicate the capability and robustness of the proposed PINN methodology for solving multi-physics problems, irrespective of the presence of sparse and noisy data in the training dataset.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2018

An Investigation of Chaotic Diffusion in a Family of Hamiltonian Mappings Whose Angles Diverge in the Limit of Vanishingly Action

The chaotic diffusion for a family of Hamiltonian mappings whose angles diverge in the limit of vanishingly action is investigated by using the solution of the diffusion equation. The system is described by a two-dimensional mapping for the variables action, I, and angle, θ and controlled by two control parameters: (i) ϵ, controlling the nonlinearity of the system, particularly a transition from integrable for ϵ=0 to non-integrable for ϵ≠0 and; (ii) γ denoting the power of the action in the equation defining the angle. For ϵ≠0 the phase space is mixed and chaos is present in the system leading to a finite diffusion in the action characterized by the solution of the diffusion equation. The analytical solution is then compared to the numerical simulations showing a remarkable agreement between the two procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗