SEARCH · Search NASA

Results for “anisotropic 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Transient anisotropic kernel for probabilistic learning on manifolds

PLoM (Probabilistic Learning on Manifolds) is a method introduced in 2016 for handling small training datasets by projecting an Itô equation from a stochastic dissipative Hamiltonian dynamical system, acting as the MCMC generator, for which the KDE-estimated probability measure with the training dataset is the invariant measure. PLoM performs a projection on a reduced-order vector basis related to the training dataset, using the diffusion maps (DMAPS) basis constructed with a time-independent isotropic kernel. In this paper, we propose a new ISDE projection vector basis built from a transient anisotropic kernel, providing an alternative to the DMAPS basis to improve statistical surrogates for stochastic manifolds with heterogeneous data. The construction ensures that for times near the initial time, the DMAPS basis coincides with the transient basis. For larger times, the differences between the two bases are characterized by the angle of their spanned vector subspaces. The optimal instant yielding the optimal transient basis is determined using an estimation of mutual information from Information Theory, which is normalized by the entropy estimation to account for the effects of the number of realizations used in the estimations. Consequently, this new vector basis better represents statistical dependencies in the learned probability measure for any dimension. Three applications with varying levels of statistical complexity and data heterogeneity validate the proposed theory, showing that the transient anisotropic kernel improves the learned probability measure.

Diffusion maps↗

DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1978

Two-stream theory of spectral reflectance of snow

Spectral reflectance of snow under diffuse illumination is studied using the two-stream approximation of the radiative transfer equation. The scattering and absorption within the snowcover due to the randomly distributed ice grains are characterized by the single scattering albedo and anisotropic phase function. Geometric optics calculations are used to relate the scattering and absorption parameters to grain size and density of snow. Analytical expressions for the intensity within the snowpack and the asymptotic flux extinction coefficient are also obtained. Good agreement is shown between the theory and available experimental data on visible and near-infrared reflectance and asymptotic flux extinction coefficient. The theory also may be used to explain the observed effect of aging on the snow reflectance.

Choudhury, B. J.↗

NASA NTRS2015

High-Order Hyperbolic Residual-Distribution Schemes on Arbitrary Triangular Grids

In this paper, we construct high-order hyperbolic residual-distribution schemes for general advection-diffusion problems on arbitrary triangular grids. We demonstrate that the second-order accuracy of the hyperbolic schemes can be greatly improved by requiring the scheme to preserve exact quadratic solutions. We also show that the improved second-order scheme can be easily extended to third-order by further requiring the exactness for cubic solutions. We construct these schemes based on the LDA and the SUPG methodology formulated in the framework of the residual-distribution method. For both second- and third-order-schemes, we construct a fully implicit solver by the exact residual Jacobian of the second-order scheme, and demonstrate rapid convergence of 10-15 iterations to reduce the residuals by 10 orders of magnitude. We demonstrate also that these schemes can be constructed based on a separate treatment of the advective and diffusive terms, which paves the way for the construction of hyperbolic residual-distribution schemes for the compressible Navier-Stokes equations. Numerical results show that these schemes produce exceptionally accurate and smooth solution gradients on highly skewed and anisotropic triangular grids, including curved boundary problems, using linear elements. We also present Fourier analysis performed on the constructed linear system and show that an under-relaxation parameter is needed for stabilization of Gauss-Seidel relaxation.

Mazaheri, Alireza↗

DOE OSTI2024

A fully implicit, asymptotic-preserving, semi-Lagrangian algorithm for the time dependent anisotropic heat transport equation

In this paper, we extend the operator-split asymptotic-preserving, semi-Lagrangian algorithm for time dependent anisotropic heat transport equation proposed in Chacón et al. (2014) [18] to use a fully implicit time integration with backward differentiation formulas. The proposed implicit method can deal with arbitrary heat-transport anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ $\ggg$ 1 (with $\mathcal{X}$∥, $ \mathcal{X}$⟂ the parallel and perpendicular heat diffusivities, respectively) in complicated magnetic field topologies in an accurate and efficient manner. Further, the implicit algorithm is second-order accurate temporally and demonstrates an accurate treatment at boundary layers (e.g., island separatrices), which was not ensured by the operator-split implementation. The condition number of the resulting algebraic system is independent of the anisotropy ratio, and is inverted with preconditioned GMRES. We propose a simple preconditioner that renders the finite-dimensional linear operator compact, resulting in mesh-independent convergence rates for topologically simple magnetic fields, and convergence rates scaling as ~ (NΔt) 1/4 (with N the total mesh size and Δt the timestep) in topologically complex magnetic-field configurations. We demonstrate the accuracy and performance of the approach with test problems of varying complexity, including an analytically tractable boundary-layer problem in a straight magnetic field, and a topologically complex magnetic field featuring magnetic islands with extreme anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ = 10 10 ) .

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Inertial and anisotropic pressure effects on cross-field electron transport in low-temperature magnetized plasmas

In this paper, a one-dimensional (1D) particle-in-cell Monte Carlo collision (PIC-MCC) model is developed to investigate the effects of anisotropic pressure and inertial terms due to non-Maxwellian velocity distribution functions on cross-field electron transport. The conservation of momentum is evaluated by taking the moments of the first-principles gas-kinetic equation. A steady-state discharge is obtained without any low-frequency ionization oscillations by considering an anomalous electron scattering profile. Further, the results obtained from the 1D PIC-MCC model are compared with fluid models, including the quasi-neutral drift-diffusion (DD), non-neutral DD, and full fluid moment models. The discharge current obtained from the PIC-MCC model is in good agreement with the fluid models. The cross-field electron transport due to the inertial terms, i.e. the gradient of axial and azimuthal drift, is evaluated. Moreover, PIC-MCC simulation results show non-zero, anisotropic, off-diagonal pressure tensor terms due to asymmetric non-Maxwellian electron velocity distribution function, potentially contributing to cross-field electron transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1979

Evaluation of models to predict insolation on tilted surfaces

An empirical study was performed to evaluate the validity of various insolation models which employ either an isotropic or an anisotropic distribution approximation for sky light when predicting insolation on tilted surfaces. Data sets of measured hourly insolation values were obtained over a 6-month period using pyranometers which received diffuse and total solar radiation on a horizontal plane and total radiation on surfaces tilted toward the equator at 37 degrees and 60 degrees angles above the horizon. Data on the horizontal surfaces were used in the insolation models to predict insolation on the tilted surface; comparisons of measured vs calculated insolation on the tilted surface were examined to test the validity of the sky light approximations. It was found that the Liu-Jordan isotropic distribution model provides a good fit to empirical data under overcast skies but underestimates the amount of solar radiation incident on tilted surfaces under clear and partly cloudy conditions.

Klucher, T. M.↗

DOE OSTI2022

Mesoscopic finite-size effects of unconventional electron transport in PdCoO 2

A wide range of unconventional transport phenomena has recently been observed in single-crystal delafossite metals. Here, we present a theoretical framework to elucidate electron transport using a combination of first-principles calculations and numerical modeling of the anisotropic Boltzmann transport equation. Using PdCoO 2 as a model system, we study different microscopic electron and phonon scattering mechanisms and establish the mean free path hierarchy of quasiparticles at different temperatures. We treat the anisotropic Fermi surface explicitly to numerically obtain experimentally-accessible transport observables, which bridge between the “diffusive,” “ballistic,” and “hydrodynamic” transport regime limits. Further, we illustrate that the distinction between the “quasiballistic” and “quasihydrodynamic” regimes is challenging and often needs to be quantitative in nature. From first-principles calculations, we populate the resulting transport regime plots and demonstrate how the Fermi surface orientation adds complexity to the observed transport signatures in micrometer-scale devices. Our work provides key insights into microscopic interaction mechanisms on open hexagonal Fermi surfaces and establishes their connection to the macroscopic electron transport in finite-size channels.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Aggregation of solutes in bosonic versus fermionic quantum fluids

Quantum fluid droplets made of helium-3 ( 3 He) or helium-4 ( 4 He) isotopes have long been considered as ideal cryogenic nanolabs, enabling unique ultracold chemistry and spectroscopy applications. The droplets were believed to provide a homogeneous environment in which dopant atoms and molecules could move and react almost as in free space but at temperatures close to absolute zero. Here, we report ultrafast x-ray diffraction experiments on xenon-doped 3 He and 4 He nanodroplets, demonstrating that the unavoidable rotational excitation of isolated droplets leads to highly anisotropic and inhomogeneous interactions between the host matrix and enclosed dopants. Superfluid 4 He droplets are laced with quantum vortices that trap the embedded particles, leading to the formation of filament-shaped clusters. In comparison, dopants in 3 He droplets gather in diffuse, ring-shaped structures along the equator. The shapes of droplets carrying filaments or rings are direct evidence that rotational excitation is the root cause for the inhomogeneous dopant distributions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

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

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

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Mechanochemical topological defects in an active nematic

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

42 ENGINEERING↗

NASA NTRS1976

A finite difference method of solving anisotropic scattering problems

A new method of solving radiative transfer problems is described including a comparison of its speed with that of the doubling method, and a discussion of its accuracy and suitability for computations involving variable optical properties. The method uses a discretization in angle to produce a coupled set of first-order differential equations which are integrated between discrete depth points to produce a set of recursion relations for symmetric and anti-symmetric angular sums of the radiation field at alternate depth points. The formulation given here includes depth-dependent anisotropic scattering, absorption, and internal sources, and allows arbitrary combinations of specular and non-Lambertian diffuse reflection at either or both boundaries. Numerical tests of the method show that it can return accurate emergent intensities even for large optical depths. The method is also shown to conserve flux to machine accuracy in conservative atmospheres

Barkstrom, B. R.↗

DOE OSTI2022

A model of thermal creep and annealing in finite domains based on coupled dislocation climb and vacancy diffusion

Here, we develop a framework to investigate thermal creep and annealing in finite domains, where the climb motion of discrete dislocations is coupled to the diffusion of a continuum vacancy field. The model is first formulated in a continuum finite-deformation setting. All governing equations and boundary conditions are obtained from a unified irreversible thermodynamics principle. The resulting model couples a mechanical boundary value problem (BVP), a vacancy diffusion BVP, and the climb and glide motion of the discrete dislocation network within the crystal. The framework is then linearized for implementation in three-dimensional (3D) discrete dislocation dynamics (DDD) simulations for arbitrary anisotropic crystals. A solution scheme is developed based on the superposition principle, which is imposed weakly on the dislocation network to obtain a Galerkin solution for the nodal climb velocities. The framework includes diffusional (Nabarro–Herring) creep deformation as well as dislocation creep by climb-assisted-glide. The method is applied to simulate the annealing of vacancy loops in Al, with good agreement to experimental measurements by Silcox and Hirsch. We further consider the effects of annealing under stress, and of the proximity of the vacancy loops to loaded and free boundaries Simulations in polycrystalline materials are carried out to highlight the effects of the grain size on dislocation climb and vacancy loop annealing. The method is also applied to estimate the creep rate due to climb-assisted glide of jogged-screw dislocations in γ-TiAl, and results are compared to experiments by Viswanathan et al. Finally, we discuss the effects of uniaxial and hydrostatic stresses on the two diffusive deformation pathways of the material, namely Nabarro–Herring creep and dislocation climb.

36 MATERIALS SCIENCE↗

NASA NTRS1984

Turbulent solutions of the equations of fluid motion

Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. Initial three-dimensional cosine velocity fluctuations and periodic boundary conditions are used in most of the work considered. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations, the initially nonrandom flow develops into an apparently random turbulence. Thus randomness or turbulence can arise as a consequence of the structure of the Navier-Stokes equations. The cases considered include turbulence which is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A mean shear is present in some cases. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components, and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Deissler, R. G.↗

DOE OSTI2022

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗

NASA NTRS1987

Hydrodynamical constraints on cosmic-ray acceleration in relativistic shocks

A two-fluid hydrodynamical model governing the transport of cosmic rays in a relativistically moving background plasma is developed. The equations are used to discuss the time asymptotic structure of a relativistic, plane-parallel shock wave significantly modified by first-order Fermi acceleration of cosmic rays. The model allows for an anisotropic cosmic-ray pressure tensor with pressure components parallel and perpendicular to the shock normal. It is surmised that the well-known energy threshold problem for diffusive shock acceleration that downstream particles have a large enough velocity to overtake the shock and subsequently scatter in the upstream medium before returning to the downstream medium may be expressed hydrodynamically by the condition V less than V(cr), where V is the fluid velocity relative to the shock and V(cr) is the relativistic cosmic-ray sound speed. Astrophysical implications of the results are briefly discussed.

Webb, G. M.↗

NASA NTRS2001

2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.

Svizhenko, Alexei↗