SEARCH · Search NASA

Results for “Numerical diffusion”

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 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Electron transport in a tokamak scrape-off layer: impact of toroidal nonuniformities of divertor targets

An analysis of electron dynamics is carried out in the tokamak scrape-off layer (SOL) region. Small non-axisymmetric perturbations of the divertor target sheath potential affect electron drift orbits in the SOL and may lead to non-ambipolar electron radial transport. The resulting rates of electron convection, diffusion, and mobility are calculated analytically, and the analytic results are compared with direct numerical simulations of electron drift orbits. The proposed mechanism of non-ambipolar electron transport may be relevant to the sustainment of the SOL plasma quasi-neutrality in the ‘heuristic’ model of the SOL width (Eich et al 2011 Phys. Rev. Lett. 107 215001).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

A Least-Squares-Based Weak Galerkin Finite Element Method for Second Order Elliptic Equations

Here, in this article, we introduce a least-squares-based weak Galerkin finite element method for the second order elliptic equation. This new method is shown to provide very accurate numerical approximations for both the primal and the flux variables. In contrast to other existing least-squares finite element methods, this new method allows us to use discontinuous approximating functions on finite element partitions consisting of arbitrary polygon/polyhedron shapes. We also develop a Schur complement algorithm for the resulting discretization problem by eliminating all the unknowns that represent the solution information in the interior of each element. Optimal order error estimates for both the primal and the flux variables are established. An extensive set of numerical experiments are conducted to demonstrate the robustness, reliability, flexibility, and accuracy of the least-squares-based weak Galerkin finite element method. Finally, the numerical examples cover a wide range of applied problems, including singularly perturbed reaction-diffusion equations and the flow of fluid in porous media with strong anisotropy and heterogeneity.

nite element methods↗

DOE OSTI2020

Retention of Long-period Gas Giant Planets: Type II Migration Revisited

During their formation, emerging protoplanets tidally interact with their natal disks. Proto–gas giant planets, with Hill radii larger than the disk thickness, open gaps and quench gas flow in the vicinity of their orbits. It is usually assumed that their type II migration is coupled to the viscous evolution of the disk. Although this hypothesis provides an explanation for the origin of close-in planets, it also encounters a predicament on the retention of long-period orbits for most gas giant planets. Moreover, numerical simulations indicate that the planets' migrations are not solely determined by the viscous diffusion of their natal disk. In this work, we carry out a series of hydrodynamic simulations combined with analytic studies to examine the transition between different paradigms of type II migration. We find a range of planetary mass for which gas continues to flow through a severely depleted gap so that the surface density distribution in the disk region beyond the gap is maintained in a quasi-steady state. The associated gap profile modifies the location of corotation and Lindblad resonances. In the proximity of the planet's orbit, high-order Lindblad and corotation torque are weakened by the gas depletion in the gap, while low-order Lindblad torques near the gap walls preserve their magnitude. Consequently, the intrinsic surface density distribution of the disk delicately determines both the pace and direction of the planets' type II migration. We show that this effect might stall the inward migration of giant planets and preserve them in disk regions where the surface density is steep.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1987

Ion and electron velocity distributions within flux transfer events

The detailed nature of the thermal and suprathermal ion and electron distributions within magnetic flux transfer events (FTEs) is examined. Examples of both magnetosheath FTEs and magnetospheric FTEs are discussed. The detailed distributions confirm that FTEs contain a mixture of magnetosheath and magnetospheric plasmas. To lowest order, the distributions are consistent with a simple superposition of the two interpenetrating populations, with no strong interactions between them. To first order, some interesting differences appear, especially in the electron distributions, suggesting that considerable pitch angle scattering and some electron energy diffusion are also occurring. These observations should provide a useful test of analytical and numerical studies of interpenetrating plasmas.

Thomsen, M. F.↗

NASA NTRS1993

A study of the compressible flow through a diffusing S-duct

Benchmark aerodynamic data are presented for compressible flow through a representative S-duct configuration. A numerical prediction of the S-duct flow field, obtained from a subsonic parabolized Navier-Stokes algorithm, is also shown. The experimental and numerical results are compared. Measurements of the three-dimensional velocity field, total pressures, and static pressures were obtained at five cross-sectional planes. Aerodynamic data were gathered with calibrated pneumatic probes. Surface static pressure and surface flow visualization data were also acquired. All reported tests were conducted with an inlet centerline Mach number of 0.6. The Reynolds number, based on the inlet centerline velocity and duct inlet diameter, was 2.6 x 10(exp 6). Thin inlet turbulent boundary layers existed. The collected data should be beneficial to aircraft inlet designers and the measurements are suitable for the validation of computational codes. The results show that a region of streamwise flow separation occurred within the duct. Details about the separated flow region, including mechanisms which drive this complicated flow phenomenon, are discussed. Results also indicate that the duct curvature induces strong pressure driven secondary flows. The cross flows evolve into counter-rotating vortices. These vortices convect low momentum fluid of the boundary layer toward the center of the duct, degrading both the uniformity and magnitude of the total pressure profile.

Wellborn, Steven R.↗

DOE OSTI2020

Comparison of plug flow and multi-node stratified tank modeling approaches regarding computational efficiency and accuracy

Residential water heaters contain water stratified by temperature-driven density differences. This implies that a water tank can reach a state in which the top and bottom sections have different temperatures, unless mixing happens. A high degree of thermal stratification can improve the efficiency of some water heaters, by saving the amount of energy required for the heat-up process. Studies of stratification became popular in the 1970s and it remains an active research topic today. The research has led to the development of different models and techniques to better predict and define a stratified tanks behavior. By comparing these models and techniques used previously to describe thermal stratification, the phenomenon could be better understood, exploited, and used to increase efficiency and thermal energy capacity in modern water tanks. From the existing models, we found the one-dimensional standard plug-flow and a multi node model to be appropriate for analyzing the processes of the heat up and cool-down in a water tank. These two models are based on energy balances. This work involved comparing the accuracy and computational effort needed to implement these models. To assess accuracy, we compared both types of existing models to experimental data (also collected in this work) which included a heat up process using an external heat pump. This external process included a layering process that has an eddy diffusivity at five times the rate of thermal diffusion. For this project, we implemented the models in MATLAB, the multi-paradigm numerical computing environment. We quantified model accuracy using the root mean squared error between modeled data and experimental data for six measured tank temperatures. Comparing the accuracy and the computational time taken to run the simulation provides a method to contrast the performance of each model and a way to rate it. The multi node model was run using from 6 to 96 spatial nodes; the plug flow model was run using 1 to 0.001 º C temperature bin sizes. Additionally, timesteps were varied from 4 to 236 s. The results quantify the tradeoff between accuracy and computational time, providing guidance for simulations to intelligently select the best model type and simulation parameters. This research can be used to validate the pre-existing models and possibly improve the modern water tank.

Bulnes, Fernando Karg↗

DOE OSTI2018

Fluctuating Hydrodynamic Methods for Fluid-Structure Interactions in Confined Channel Geometries

Computational methods are developed for the study of fluid-structure interactions subject to thermal fluctuations when conned within channels having slit-like geometry. The methods take into account the hydrodynamic coupling and diusivity of microstructures when influenced by their proximity to no-slip walls. We develop stochastic numerical methods subject to no-slip boundary conditions using a staggered nite volume discretization. We show that by imposing an exact fluctuation-dissipation condition to discretize the stochastic driving elds combined with using an exact projection method to enforce incompressibility is sucient to ensure results consistent with statistical mechanics. We demonstrate our methods by investigating how the proximity of ellipsoidal colloids to the channel wall eects their active hydrodynamic responses and passive diusivity. We also study the collective dynamics of a large number of particles by considering the intermediate scattering functions for the relaxation of density fluctuations. We expect our introduced stochastic computational methods to be applicable broadly to applications in which confinement effects play an important role in the dynamics of microstructures subject to hydrodynamic coupling and thermal fluctuations. uctuations.

Wang, Yaohong↗

NASA NTRS1992

Laboratory scaled simulation of lidar cloud sounding experiments

The results of lidar measurements carried out on laboratory scale models of clouds are presented. Measurements on laboratory scale models are important since one has the knowledge of the relevant parameters of the diffusing medium, such as: scattering and absorption coefficients, phase function, homogeneity, shape, etc. Knowledge of these parameters enables one to use the results to test the reliability of theoretical and numerical investigations. To obtain a laboratory scaled model of a lidar system sounding a cloud, it is necessary to scale down all the geometrical quantities by the same factor to reduce distances of the order of kilometers to the order of meters, keeping the size and the optical depth of the diffusers unchanged. If a time resolution of the order of nanoseconds is necessary for a lidar sounding actual clouds, the corresponding time resolution for the laboratory model should be of the order of picoseconds. It is possible to obtain this resolution by using picosecond laser systems and fast electrooptical detectors like the streak camera. The results of the laboratory measurements showed that the multiple scattering effect strongly depends on the size of the diffusers, as well as on the concentration. The experimental results were compared with the numerical results of a Monte Carlo code. A generally good agreement was obtained.

Zaccanti, G.↗

DOE OSTI2025

Diffusion Quantum Monte Carlo Benchmarking of Magnetic Moments in MnBi 2 Te 4

The intrinsically antiferromagnetic topological insulator, MnBi 2 Te 4 (MBT), has garnered significant attention recently due to its potential to host numerous exotic topological quantum states. Unfortunately, their consistent realization has been hindered by intrinsic antisite defects among the Mn and Bi sublattices. In this work, we establish Mn magnetization of pristine MBT through high level diffusion Monte Carlo calculations, which can serve as a precise starting point for various models to estimate antisite defect concentrations in actual MBT samples. The benchmark quality of DMC calculations is further identified from out model estimating antisite defect concentrations, which combines the benchmarked Mn magnetization with data from magnetic susceptibility and intermediate field magnetization measurements. This reproduces well Bi Mn and Mn Bi concentrations measured in the experiments. Here, we anticipate these theoretically based magnetic purity measures may be used as minimization targets in cycles of refinement to synthesize MBT with low antisite defect concentrations and more reproducible topological properties.

Defects↗

NASA NTRS1980

Magnetic pumping of particles in the outer Jovian magnetosphere

The mechanism of magnetic pumping consists of two processes, the adiabatic motion of charged particles in a time varying magnetic field and their pitch-angle diffusion. The result is a systematic increase in the energy of charged particles trapped in mirror (and particularly, magnetospheric) magnetic fields. A numerical model of the mechanism is constructed, compared with analytic theory where possible, and, through elementary exercises, is used to predict the consequences of the process for cases that are not tractable by analytical means. For energy dependent pitch angle diffusion rates, characteristic 'two temperature' distributions are produced. Application of the model to the outer Jovian magnetosphere shows that beyond 20 Jupiter radii in the outer magnetosphere, particles may be magnetically pumped to energies of the order of 1 - 2 MeV. Two temperature distribution functions with "break points" at 1 - 4 KeV for electrons and 8 - 35 KeV for ions are predicted.

Borovsky, J. E.↗

NASA NTRS1990

A turbulent two-phase flow model for nebula flows

A new and very efficient turbulent two-phase flow numericaly model is described to analyze the environment of a protoplanetary nebula at a stage prior to the formation of planets. Focus is on settling processes of dust particles in flattened gaseous nebulae. The model employs a perturbation technique to improve the accuracy of the numerical simulations of such flows where small variations of physical quantities occur over large distance ranges. The particles are allowed to be diffused by gas turbulence in addition to settling under gravity. Their diffusion coefficients is related to the gas turbulent viscosity by the non-dimensional Schmidt number. The gas turbulent viscosity is determined by the means of the eddy viscosity hypothesis that assumes the Reynolds stress tensor proportional to the mean strain rate tensor. Zero- and two-equation turbulence models are employed. Modeling assumptions are detailed and discussed. The numerical model is shown to reproduce an existing analytical solution for the settling process of particles in an inviscid nebula. Results of nebula flows are presented taking into account turbulence effects of nebula flows. Diffusion processes are found to control the settling of particles.

Champney, Joelle M.↗

NASA NTRS2017

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗

DOE OSTI2021

Spectral Denoising for Accelerated Analysis of Correlated Ionic Transport

Computation of correlated ionic transport properties from molecular dynamics in the Green-Kubo formalism is expensive, as one cannot rely on the affordable mean square displacement approach. We use spectral decomposition of the short-time ionic displacement covariance to learn a set of diffusion eigenmodes that encode the correlation structure and form a basis for analyzing the ionic trajectories. This allows systematic reduction of the uncertainty and accelerate computations of ionic conductivity in systems with a steady-state correlation structure. We provide mathematical and numerical proofs of the method’s robustness and demonstrate it on realistic electrolyte materials.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2018

Сoncentric GaAs Nanorings Growth Modelling

The nanostructures formation process using the droplet epitaxy technique was investigated by Monte Carlo simulation. The simulation was fulfilled for two-dimensional and three-dimensional geometry substrates. The nanostructures morphology dependence on the growth temperature was presented. Crystal clusters, single and double rings were observed. The nanostructures shape was shown to be determined by the gallium diffusion length. The conditions of double rings formation during the droplet epitaxy were considered using analytical and numerical approaches. The factors that determine the rings location and shape were analyzed. The growth morphology was demonstrated to be dependent on the initial distance L between the droplets. The double ring formation was possible at a low droplet density only, when the As-stabilized region could be created between the droplets.

36 MATERIALS SCIENCE↗

DOE OSTI2018

A reaction–diffusion model for phenotypic evolution

We present a reaction–diffusion mathematical model for the evolutionary dynamics of phenotypic evolution. A detailed deduction of the equations is presented for the one-dimensional version, from which a more general model is proposed. Particular cases are studied using analytical approximations and numerical simulations. Results indicate that the approach proposed produces results that are coherent with mainstream models in evolutionary dynamics, suggesting that the reaction–diffusion model could be an alternative tool in the analysis of evolutionary dynamics.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1988

SINDA - Systems Improved Numerical Differencing Analyzer

Capability to model fluid flow added. SINDA is software system developed to solve physical problems governed by diffusion-type equations mathematically modeled by lumped-parameter representations. Used as general thermal analyzer with resistor-and-capacitor-network representations. Adapted to wide range of problems represented by Fourier, Possion, Laplace, or other differential equations. Written in FORTRAN and Assembler.

Damico, Steve J.↗

NASA NTRS1995

An Experimental and Theoretical Study of Radiative Extinction of Diffusion Flames

The objective of this research was to experimentally and theoretically investigate the radiation-induced extinction of gaseous diffusion flames in microgravity. The microgravity conditions were required because radiation-induced extinction is generally not possible in 1-g but is highly likely in microgravity. In 1-g, the flame-generated particulates (e.g. soot) and gaseous combustion products that are responsible for flame radiation, are swept away from the high temperature reaction zone by the buoyancy-induced flow and a steady state is developed. In microgravity, however, the absence of buoyancy-induced flow which transports the fuel and the oxidizer to the combustion zone and removes the hot combustion products from it enhances the flame radiation due to: (1) transient build-up of the combustion products in the flame zone which increases the gas radiation, and (2) longer residence time makes conditions appropriate for substantial amounts of soot to form which is usually responsible for most of the radiative heat loss. Numerical calculations conducted during the course of this work show that even non-radiative flames continue to become "weaker" (diminished burning rate per unit flame area) due to reduced rates of convective and diffusive transport. Thus, it was anticipated that radiative heat loss may eventually extinguish the already "weak" microgravity diffusion flame. While this hypothesis appears convincing and our numerical calculations support it, experiments for a long enough microgravity time could not be conducted during the course of this research to provide an experimental proof. Space shuttle experiments on candle flames show that in an infinite ambient atmosphere, the hemispherical candle flame in microgravity will burn indefinitely. It was hoped that radiative extinction can be experimentally shown by the aerodynamically stabilized gaseous diffusion flames where the fuel supply rate was externally controlled. While substantial progress toward this goal was made during this project, identifying the experimental conditions for which radiative extinction occurs for various fuels requires further study. Details concerning this research which are discussed in published articles are included in the appendices.

Atreya, Arvind↗

DOE OSTI2018

A Study of the von Neumann Stability Analysis of a Semi-implicit Numerical Method Applied Over the Radial Impurity Transport Equation in Tokamak Plasma

The radial impurity transport equation for tokamak plasma is a form of diffusion–convection–reaction equation. The impurity transport equation is solved to determine the distribution of impurity (non-fuel) ion species with different ionization states perpendicular to magnetic surfaces of tokamak plasma. The equation for each charge (ionization) state Z is a non-linear, second-order in space, first-order in time, parabolic partial differential equation coupled to the previous Z − 1 and the next Z + 1 charge states of the impurity species through its reaction term. The number of differential equations to be solved simultaneously is hence determined by the number of ionization states of the impurity species studied. The solution to the set of these coupled equations can be obtained using a semi-implicit numerical method applied on it. The present study describes the application of von Neumann stability analysis over the semi-implicit numerical method applied over the radial impurity transport equation and determines a generic stability criterion for the method. The stability analysis is further illustrated using the geometry of Aditya tokamak installed at the Institute for Plasma Research Gandhinagar, India as an example. The impurity species considered is oxygen (Atomic number = 8). This leads to a set of eight coupled equations for charge states Z = 1 to 8 over which von Neumann analysis is illustrated in present study.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗