Search NASA⌕ Search

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.

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

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Space-energy separated representations for multigroup neutron diffusion using proper generalized decompositions

Proper generalized decomposition (PGD) can be described as a numerical extension of separation of variables and thus be employed as a solution technique for multi-dimensional problems, where dimensions should be understood in the phase-space context of the governing law at hand. This paper presents a PGD approach on efficiently solving problems involving multigroup neutron diffusion. Two PGD approaches are described: space-only and space-energy decompositions. Numerical results include 2-D and 3-D examples with two-, seven-, and 145-group structures. The PGD solutions are compared against traditional multi-dimensional finite element discretization. For few-group problems, both PGD approaches prove effective for mildly heterogeneous geometries, but showed reduced performance with increasing heterogeneity. The space-energy representation was found to be slower than the space-only approach for two-group problems, but proved more effective for seven-group problems. For even larger numbers of groups, the space-energy PGD decomposition was very effective at reducing the computational time.

42 ENGINEERING↗
NASA NTRS2001

PIV Measurements in Weakly Buoyant Gas Jet Flames

Despite numerous experimental investigations, the characterization of microgravity laminar jet diffusion flames remains incomplete. Measurements to date have included shapes, temperatures, soot properties, radiative emissions and compositions, but full-field quantitative measurements of velocity are lacking. Since the differences between normal-gravity and microgravity diffusion flames are fundamentally influenced by changes in velocities, it is imperative that the associated velocity fields be measured in microgravity flames. Velocity measurements in nonbuoyant flames will be helpful both in validating numerical models and in interpreting past microgravity combustion experiments. Pointwise velocity techniques are inadequate for full-field velocity measurements in microgravity facilities. In contrast, Particle Image Velocimetry (PIV) can capture the entire flow field in less than 1% of the time required with Laser Doppler Velocimetry (LDV). Although PIV is a mature diagnostic for normal-gravity flames , restrictions on size, power and data storage complicate these measurements in microgravity. Results from the application of PIV to gas jet flames in normal gravity are presented here. Ethane flames burning at 13, 25 and 50 kPa are considered. These results are presented in more detail in Wernet et al. (2000). The PIV system developed for these measurements recently has been adapted for on-rig use in the NASA Glenn 2.2-second drop tower.

Sunderland, Peter B.↗
NASA NTRS2007

Effects of Lewis Number on Temperatures of Spherical Diffusion Flames

Spherical diffusion flames supported on a porous sphere were studied numerically and experimentally. Experiments were performed in 2.2 s and 5.2 s microgravity facilities. Numerical results were obtained from a Chemkin-based program. The program simulates flow from a porous sphere into a quiescent environment, yields both steady-state and transient results, and accounts for optically thick gas-phase radiation. The low flow velocities and long residence times in these diffusion flames lead to enhanced radiative and diffusive effects. Despite similar adiabatic flame temperatures, the measured and predicted temperatures varied by as much as 700 K. The temperature reduction correlates with flame size but characteristic flow times and, importantly, Lewis number also influence temperature. The numerical results show that the ambient gas Lewis number would have a strong effect on flame temperature if the flames were steady and nonradiating. For example, a 10% decrease in Lewis number would increase the steady-state flame temperature by 200 K. However, for these transient, radiating flames the effect of Lewis number is small. Transient predictions of flame sizes are larger than those observed in microgravity experiments. Close agreement could not be obtained without either increasing the model s thermal and mass diffusion properties by 30% or reducing mass flow rate by 25%.

Santa, K. J.↗
DOE OSTI2024

A Finite Difference informed Random Walk solver for simulating radiation defect evolution in polycrystalline structures with strongly inhomogeneous diffusivity

Diffusivity of species and defects on grain boundaries is usually several orders of magnitude larger than that inside grains. Such strongly inhomogeneous diffusivity requires prohibitively high computational demands for modeling microstructural evolution. Here, this paper presents a highly-efficient numerical solver, combining the Finite Difference method and Random Walk model, designed for accurately modeling strongly inhomogeneous diffusion within polycrystalline structures. The proposed solver, termed Finite Difference informed Random Walk (FDiRW), integrates a customized Finite Difference (cFD) scheme tailored for fast diffusion along thin grain boundaries represented by a single-layer of nodes. Numerical experiments demonstrate that the FDiRW solver achieves an impressive efficiency gain of 1560x compared to traditional Finite Difference methods while maintaining accuracy, making it feasible for personal computer machines to handle diffusional systems with strongly inhomogeneous diffusivity across static polycrystalline microstructures. The model has been successfully applied to simulate radiation defect evolution, showcasing its scalability to engineering scales in both length and time dimensions.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Advection-Diffusion Problems

Test problems are used to examine the performance of several one-dimensional numerical schemes based on the space-time conservation and solution element (CE/SE) method. Investigated in this paper are the CE/SE schemes constructed previously for solving the linear unsteady advection-diffusion equation and the schemes derived here for solving the nonlinear viscous and inviscid Burgers equations. In comparison with the numerical solutions obtained using several traditional finite-difference schemes with similar accuracy, the CE/SE solutions display much lower numerical dissipation and dispersion errors.

Wang, Xiao-Yen↗
DOE OSTI2019

Turbulent shear-layer mixing: initial conditions, and direct-numerical and large-eddy simulations

Aspects of turbulent shear-layer mixing are investigated over a range of shear-layer Reynolds numbers, $Re_{\unicode[STIX]{x1D6FF}}=\unicode[STIX]{x0394}U\unicode[STIX]{x1D6FF}/\unicode[STIX]{x1D708}$ , based on the shear-layer free-stream velocity difference, $\unicode[STIX]{x0394}U$ , and mixing-zone thickness, $\unicode[STIX]{x1D6FF}$ , to probe the role of initial conditions in mixing stages and the evolution of the scalar-field probability density function (p.d.f.) and variance. Scalar transport is calculated for unity Schmidt numbers, approximating gas-phase diffusion. The study is based on direct-numerical simulation (DNS) and large-eddy simulation (LES), comparing different subgrid-scale (SGS) models for incompressible, uniform-density, temporally evolving forced shear-layer flows. Moderate-Reynolds-number DNS results help assess and validate LES SGS models in terms of scalar-spectrum and mixing estimates, as well as other metrics, to $Re_{\unicode[STIX]{x1D6FF}}\lesssim 3.3\times 10^{4}$ . High-Reynolds-number LES investigations to $Re_{\unicode[STIX]{x1D6FF}}\lesssim 5\times 10^{5}$ help identify flow parameters and conditions that influence the evolution of scalar variance and p.d.f., e.g. marching versus non-marching. Initial conditions that generate shear flows with different mixing behaviour elucidate flow characteristics in each flow regime and identify elements that induce p.d.f. transition and scalar-variance behaviour. P.d.f. transition is found to be largely insensitive to local flow parameters, such as $Re_{\unicode[STIX]{x1D6FF}}$ , or a previously proposed vortex-pairing parameter based on downstream distance, or other equivalent criteria. The present study also allows a quantitative comparison of LES SGS models in moderate- and high- $Re_{\unicode[STIX]{x1D6FF}}$ forced shear-layer flows.

Mechanics↗
NASA NTRS2009

Thermal Transport Model for Heat Sink Design

A document discusses the development of a finite element model for describing thermal transport through microcalorimeter arrays in order to assist in heat-sinking design. A fabricated multi-absorber transition edge sensor (PoST) was designed in order to reduce device wiring density by a factor of four. The finite element model consists of breaking the microcalorimeter array into separate elements, including the transition edge sensor (TES) and the silicon substrate on which the sensor is deposited. Each element is then broken up into subelements, whose surface area subtends 10 10 microns. The heat capacity per unit temperature, thermal conductance, and thermal diffusivity of each subelement are the model inputs, as are the temperatures of each subelement. Numerical integration using the Finite in Time Centered in Space algorithm of the thermal diffusion equation is then performed in order to obtain a temporal evolution of the subelement temperature. Thermal transport across interfaces is modeled using a thermal boundary resistance obtained using the acoustic mismatch model. The document concludes with a discussion of the PoST fabrication. PoSTs are novel because they enable incident x-ray position sensitivity with good energy resolution and low wiring density.

Chervenak, James A.↗
NASA NTRS1980

Eddy diffusion coefficients and the variance of the atmosphere 30-60 km

The results of numerical models or of new observational programs are checked by comparing them with past observations. In view of the differing analysis techniques or differing data samples, the eddy diffusivities presented here agree remarkably well with past estimates. However, in the application of K-values to two-dimensional models, the actual magnitude of the diffusivities is no more important than their spatial patterns, i.e., their gradients with height and latitude. It should thus be noted that the present patterns are often much different from those of past results.

Nastrom, G. D.↗
NASA NTRS1988

On the origin of steep edges and filaments in vorticity and potential vorticity fields

High-resolution numerical calculations are shown which capture the fundamental process responsible for the intensification of vorticity gradients in an isolated vortex subject to externally imposed distrubances. Imposition of almost any weak strain or shear field in stripping away the relatively weak vorticity at the edge of the vortex and leaves it with gradients four to six orders of magnitude greater than in the initial state. Calculations displaying such enormous gradients have never been reported previously, because of the artificial eddy diffusivities that always limit such gradients in standard numerical models. The present calculations, which have no such limitations, have been made possible by the development of a novel and robust new numerical technique for vortex dynamics called contour surgery.

Dritschel, D. G.↗
NASA NTRS1982

Stagnation point heat transfer augmentation due to free stream turbulence

A model has been proposed for the momentum eddy diffusivity induced by free stream turbulence intensity and integral length scale. The eddy diffusivity model is applied to the stagnation point of a cylinder situated in a uniform crossflow in the presence of free stream turbulence. A numerical solution of the governing momentum and energy equations with the proposed eddy diffusivity model yielded results for the skin friction coefficient and the Nusselt number. The numerical predictions of the present work are compared with experimental data and the agreement between the two is seen to be very good.

Gorla, R. S. R.↗
NASA NTRS1981

A computer program for the simulation of heat and moisture flow in soils

A computer program that simulates the flow of heat and moisture in soils is described. The space-time dependence of temperature and moisture content is described by a set of diffusion-type partial differential equations. The simulator uses a predictor/corrector to numerically integrate them, giving wetness and temperature profiles as a function of time. The simulator was used to generate solutions to diffusion-type partial differential equations for which analytical solutions are known. These equations include both constant and variable diffusivities, and both flux and constant concentration boundary conditions. In all cases, the simulated and analytic solutions agreed to within the error bounds which were imposed on the integrator. Simulations of heat and moisture flow under actual field conditions were also performed. Ground truth data were used for the boundary conditions and soil transport properties. The qualitative agreement between simulated and measured profiles is an indication that the model equations are reasonably accurate representations of the physical processes involved.

Camillo, P.↗
NASA NTRS1984

Ambipolar diffusion in the middle atmosphere

In the middle atmosphere above 60 km, the electron concentration increases with altitude, reaching values of 10 to the 10th per cu m in the daytime ionospheric E region near 100 km. The electrons are more mobile than the ions and diffuse more rapidly through the neutral atmosphere. The electron diffusion polarizes the medium, causing an electric field to develop that acts to retard the electron diffusion and enhance the conduction current of ions. A global zonally averaged numerical model of atmosheric electricity from the ground to 100 km is used to examine the effect of ambipolar diffusion and the earth's geomagnetic field on the currents and fields in the middle atmosphere. The results show that above about 65 km, ambipolar diffusion generates local electric fields and conduction currents that balance electron diffusion currents. The electric fields and conduction currents are a few orders of magnitude larger than the vertical fields and currents calculated from the downward mapping of the ionospheric potential without taking electron diffusion into account. Ambipolar diffusion does not alter the total current flowing in the global circuit. It is a local effect where enhanced conduction currens flow to balance the electron diffusion current.

Tzur, I.↗
DOE OSTI2026

FROMP-from-Monomer

This repository provides the 1D reaction-diffusion simulation engine for the multi-scale framework. It takes monomer-level kinetic and thermodynamic descriptors as inputs, approximated from quantum chemical calculations, to simulate and predict the macroscopic propagation behavior of Frontal Ring-Opening Metathesis Polymerization (FROMP). By solving a coupled 1D reaction-diffusion PDE system, this tool predicts front propagation speed, front temperature, and degree of conversion. The solver explicitly captures the three fundamental ROMP steps, initiator activation/inhibition, chain initiation, and propagation, while accounting for high-temperature cycloreversion as a competing pathway. The numerical implementation uses FEniCS to handle the coupled heat diffusion and reaction kinetics.

Chua, Lauren [Massachusetts Inst. of Technology (M↗
DOE OSTI2023

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

An investigation of voids formation mechanisms and their effects on freeze and thaw processes of lithium and lithium fluoride

The mechanisms of void formation during the cooldown and freezing of lithium coolant within the primary loop of SP-100 type systems are investigated. These mechanisms are: (1) homogeneous nucleation; (2) heterogeneous nucleation; (3) normal segregation of helium gas dissolved in liquid lithium; and (4) shrinkage of lithium during freezing. To evaluate the void formation potential due to segregation, a numerical scheme that couples the freezing and mass diffusion processes in both the solid and liquid regions is developed. The results indicated that the formation of He bubbles is unlikely by either homogeneous or heterogeneous nucleation during the cooldown process. However, homogeneous nucleation of He bubbles following the segregation of dissolved He in liquid lithium ahead of the solid-liquid interface is likely to occur. Results also show that total volume of He void is insignificant when compared to that of shrinkage voids. In viewing this, the subsequent research focuses on the effects of shrinkage void forming during freezing of lithium on subsequent thaw processes are investigated using a numerical scheme that is based on a single (solid/liquid) cell approach. The cases of lithium-fluoride are also investigated to show the effect of larger volume shrinkage upon freezing on the freeze and thaw processes. Results show that a void forming at the wall appreciably reduces the solid-liquid interface velocity, during both freeze and thaw, and causes a substantial rise in the wall temperature during thaw. However, in the case of Li, the maximum wall temperature was much lower than the melting temperature of PWC-11, which is used as the structure material in the SP-100 system. Hence, it is included that a formation of hot spots is unlikely during the startup or restart of the SP-100 system.

El-Genk, Mohamed S.↗
DOE OSTI2020

Learning how structures form in drift-wave turbulence

Drift-wave turbulence produces anomalous transport via cross-correlations between fluctuations. This transport has profound implications for confinement, structure formation, and virtually all aspects of the non-linear turbulent dynamics. Herein, we use a data-driven method based on deep learning in order to study turbulent transport in the 2D Hasegawa–Wakatani system and infer a reduced mean-field model from numerical solution. In addition to the usual turbulent diffusion, we find an effect which couples the particle flux to the local gradient of vorticity, which tends to modulate the density profile. The direct coupling to the shear is relatively weak. In addition, the deep learning method finds a model for spontaneous zonal flow generation by negative viscosity, stabilized by non-linear and hyperviscous terms. We compare these results to analytic calculations using quasilinear theory and wave kinetics, finding qualitative agreement, though the calculations miss certain higher-order effects. A simplified, 1-D model for the evolution of the profile, flow, and intensity based on the deep learning results is solved numerically and compared to previous models for staircasing based on bistability. We see that the physics uncovered by the deep learning method provided simple explanations for the formation of zonal structures in the density, flow, and turbulence fields. We highlight the important role of symmetry in the deep learning method and speculate on the portability of the method to other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

A direct numerical method for predicting concentration profiles in a turbulent boundary layer over a flat plate

A numerical solution of the turbulent mass transport equation utilizing the concept of eddy diffusivity is presented as an efficient method of investigating turbulent mass transport in boundary layer type flows. A FORTRAN computer program is used to study the two-dimensional diffusion of ammonia, from a line source on the surface, into a turbulent boundary layer over a flat plate. The results of the numerical solution are compared with experimental data to verify the results of the solution. Several other solutions to diffusion problems are presented to illustrate the versatility of the computer program and to provide some insight into the problem of mass diffusion as a whole.

Dow, J. W.↗