Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion model”

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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Kinetic and Transport Modeling for Entry Flows in Hydrogen-Helium Atmospheres

Several different thermochemical models for H2/He mixtures are used to simulate a 1-D shock in a viscous CFD code at conditions relevant for hypersonic entry into the ice and gas giants. These models range from a 6-species, 1-temperature (1-T) model, to more complex 17 and 25-species state-to-state (StS) models. Corresponding radiance profiles for each of these models are then computed using a radiation code, and results are compared to experimental data from the NASA Ames Electric Arc Shock Tube (EAST) facility. Overall, the 6-species model with a quasi-steady state (QSS) solver for atomic H and modified ionization rates, along with the 17 and 25-species StS models, are found to reproduce the magnitude of the radiance seen in the EAST experiments reasonably accurately. Two additional changes to the modeling of species diffusion fluxes are also investigated: the exclusion of a ∇W term, as well as the inclusion of Soret/ Dufour effects. It is determined that these changes to the transport modeling have a negligible impact on the flowfield and therefore the radiance predictions for the H2/He shock conditions considered in this study.

Outer planets↗
DOE OSTI2026

Efficient data-driven regression for reduced-order modeling of spatial pattern formation

We present an efficient data-driven regression approach for constructing reduced-order models (ROMs) of reaction-diffusion systems exhibiting pattern formation. The ROMs are learned non-intrusively from available training data of physically accurate numerical simulations. The method can be applied to general nonlinear systems through the use of polynomial model form, while not requiring knowledge of the underlying physical model, governing equations, or numerical solvers. The process of learning ROMs is posed as a low-cost least-squares problem in a reduced-order subspace identified via Proper Orthogonal Decomposition (POD). Numerical experiments on classical pattern-forming systems–including the Schnakenberg and Mimura–Tsujikawa models–demonstrate that higher-order surrogate models significantly improve prediction accuracy while maintaining low computational cost. The proposed method provides a flexible, non-intrusive model reduction framework, well suited for the analysis of complex spatio-temporal pattern formation phenomena.

Data-driven modeling↗
DOE OSTI2019

A Method for Including Diffusive Effects in Texture Evolution

This work describes a model for representing grain misorientation in continuum models for texture evolution. The model relies on a connection, derived here, between the fluctuation of the reorientation velocity in discrete models of texture evolution and an isotropic diffusion operator that can be used to represent misorientation effects in continuum models for texture evolution. Previous work has developed a method for computing the fluctuation of the reorientation velocity in discrete models. These values, along with the average reorientation velocity, are used to set up inhomogeneous continuum advection and diffusion operators. Here, the complete homogenized continuum model uses the discrete harmonics to represent and evolve the orientation distribution function, now including a diffusion operator representing the effects of intragranular misorientation. Results from the new model are compared to fully-resolved crystal plasticity finite element simulations and experimental data to demonstrate the accuracy of the proposed diffusion operator and the overall method for simulating texture evolution.

(B) crystal plasticity↗
NASA NTRS1979

Acceleration and transport of galactic cosmic rays

A recently proposed cosmic-ray acceleration mechanism is combined with a model of transport and loss from the Galaxy to construct a model for galactic cosmic rays. Various aspects of the cosmic-ray spectrum, including the secondary-to-primary ratio as a function of energy, are found to fit naturally into the overall picture. A solution to the dynamical-halo model with energy-dependent diffusion coefficient (kappa) is presented. This model is consistent with the observed spectra if kappa is proportional to the square root of the kinetic energy. The observed flattening of the secondary-to-primary ratio below about 1-2 GeV per nucleon is shown to be a natural consequence of the model. Similarly, the primary electron spectrum is predicted to flatten below about 1-2 GeV, which is consistent with deconvolution of the local nonthermal radio spectrum, although the flattening is not quite as much as indicated by the observations.

Jokipii, J. R.↗
NASA NTRS1978

Turbulent convection within rapidly rotating superadiabatic fluids with horizontal temperature gradients

A model is developed for small scale thermally driven turbulent convection in rapidly rotating systems. The Boussinesq equations of adiabatic and inviscid flow are perturbed about a basic state which is superadiabatic and which contains a horizontal temperature gradient balanced by a geostrophic wind. The growth rates and length scales of the most unstable modes are determined. It is found that these modes are either zonally or latitudinally symmetric, depending on the latitude and the horizontal temperature gradient. The transports of heat and momentum by the modes are computed. The linear amplitudes are assumed to be limited by shear instability. The model predicts a latitude dependent eddy diffusivity. The applicability of the model to astrophysical and planetary bodies is briefly discussed.

Flasar, F. M.↗
NASA NTRS1984

Analysis techniques for tracer studies of oxidation

Analysis techniques to obtain quantitative diffusion data from tracer concentration profiles were developed. Mass balance ideas were applied to determine the mechanism of oxide growth and to separate the fraction of inward and outward growth of oxide scales. The process of inward oxygen diffusion with exchange was theoretically modelled and the effect of lattice diffusivity, grain boundary diffusivity and grain size on the tracer concentration profile was studied. The development of the tracer concentration profile in a growing oxide scale was simulated. The double oxidation technique was applied to a FeCrAl-Zr alloy using 0-18 as a tracer. SIMS was used to obtain the tracer concentration profile. The formation of lacey oxide on the alloy was discussed. Careful consideration was given to the quality of data required to obtain quantitative information.

Basu, S. N.↗
NASA NTRS2010

Combined Heat Transfer in High-Porosity High-Temperature Fibrous Insulations: Theory and Experimental Validation

Combined radiation and conduction heat transfer through various high-temperature, high-porosity, unbonded (loose) fibrous insulations was modeled based on first principles. The diffusion approximation was used for modeling the radiation component of heat transfer in the optically thick insulations. The relevant parameters needed for the heat transfer model were derived from experimental data. Semi-empirical formulations were used to model the solid conduction contribution of heat transfer in fibrous insulations with the relevant parameters inferred from thermal conductivity measurements at cryogenic temperatures in a vacuum. The specific extinction coefficient for radiation heat transfer was obtained from high-temperature steady-state thermal measurements with large temperature gradients maintained across the sample thickness in a vacuum. Standard gas conduction modeling was used in the heat transfer formulation. This heat transfer modeling methodology was applied to silica, two types of alumina, and a zirconia-based fibrous insulation, and to a variation of opacified fibrous insulation (OFI). OFI is a class of insulations manufactured by embedding efficient ceramic opacifiers in various unbonded fibrous insulations to significantly attenuate the radiation component of heat transfer. The heat transfer modeling methodology was validated by comparison with more rigorous analytical solutions and with standard thermal conductivity measurements. The validated heat transfer model is applicable to various densities of these high-porosity insulations as long as the fiber properties are the same (index of refraction, size distribution, orientation, and length). Furthermore, the heat transfer data for these insulations can be obtained at any static pressure in any working gas environment without the need to perform tests in various gases at various pressures.

Daryabeigi, Kamran↗
DOE OSTI2020

Volatilization Behavior of Supported Au Nanoparticle Arrays under H 2 at High Temperature

In this work, the fabrication and experimental study of model Au nanoparticle arrays (nanodots) was used to allow a direct comparison of reactive evaporation at elevated temperatures and nonambient H 2 -containing gas atmospheres with theoretical models. Strong temperature and H 2 -concentration dependencies were confirmed as expected, and the observed dependence in particle diameter with time was effectively described using a reactive evaporation and gas phase diffusion based quantitative model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Discontinuous Galerkin Methods for Turbulence Simulation

A discontinuous Galerkin (DG) method is formulated, implemented, and tested for simulation of compressible turbulent flows. The method is applied to turbulent channel flow at low Reynolds number, where it is found to successfully predict low-order statistics with fewer degrees of freedom than traditional numerical methods. This reduction is achieved by utilizing local hp-refinement such that the computational grid is refined simultaneously in all three spatial coordinates with decreasing distance from the wall. Another advantage of DG is that Dirichlet boundary conditions can be enforced weakly through integrals of the numerical fluxes. Both for a model advection-diffusion problem and for turbulent channel flow, weak enforcement of wall boundaries is found to improve results at low resolution. Such weak boundary conditions may play a pivotal role in wall modeling for large-eddy simulation.

Collis, S. Scott↗
NASA NTRS1997

Shapes of Buoyant and Nonbuoyant Methane Laminar Jet Diffusion Flames

Laminar gas jet diffusion flames represent a fundamental combustion configuration. Their study has contributed to numerous advances in combustion, including the development of analytical and computational combustion tools. Laminar jet flames are pertinent also to turbulent flames by use of the laminar flamelet concept. Investigations into the shapes of noncoflowing microgravity laminar jet diffusion flames have primarily been pursued in the NASA Lewis 2.2-second drop tower, by Cochran and coworkers and by Bahadori and coworkers. These studies were generally conducted at atmospheric pressure; they involved soot-containing flames and reported luminosity lengths and widths instead of the flame-sheet dimensions which are of Greater value to theory evaluation and development. The seminal model of laminar diffusion flames is that of Burke and Schumann, who solved the conservation of momentum equation for a jet flame in a coflowing ambient by assuming the velocity of fuel, oxidizer and products to be constant throughout. Roper and coworkers improved upon this model by allowing for axial variations of velocity and found flame shape to be independent of coflow velocity. Roper's suggestion that flame height should be independent of gravity level is not supported by past or present observations. Other models have been presented by Klajn and Oppenheim, Markstein and De Ris, Villermaux and Durox, and Li et al. The common result of all these models (except in the buoyant regime) is that flame height is proportional to fuel mass flowrate, with flame width proving much more difficult to predict. Most existing flame models have been compared with shapes of flames containing soot, which is known to obscure the weak blue emission of flame sheets. The present work involves measurements of laminar gas jet diffusion flame shapes. Flame images have been obtained for buoyant and nonbuoyant methane flames burning in quiescent air at various fuel flow-rates, burner diameters and ambient pressures. Soot concentrations were minimized by selecting conditions at low flowrates and low ambient pressures; this allows identification of actual flame sheets associated with blue emissions of CH and CO2. The present modeling effort follows that of Roper and is useful in explaining many of the trends observed.

Sunderland, Peter B.↗
NASA NTRS1976

On the kinetics of the pack - Aluminization process

An investigation has been made of the aluminization of unalloyed Ni in fluoride-activated packs of varying Al activity. In packs of low Al activity, in which the ratio of Al to Ni was less than 50 at. pct, the specimen surface quickly came to equilibrium with the pack and remained close to equilibrium for the duration of normal coating runs. In these packs the kinetics of aluminization was controlled by diffusion in the solid. In packs of higher Al activity the surface of the specimen did not come to equilibrium with the pack and the kinetics of the process was governed by a combination of solid and gas diffusion rates. Under most conditions however, the surface composition was time-invariant and a steady-state appeared to exist at the pack-coating interface. By combining Levine and Caves' model for gaseous diffusion in pure-Al packs with calculations of solid diffusion rates some success has been achieved in explaining the results.

Sivakumar, R.↗
DOE OSTI2018

UEDGE modelling of detached divertor operation for long-leg divertor geometries in ARC

In this work, simulations of the ARC scrape-off layer and divertor configuration are carried out using UEDGE to perform the first power-handling assessment for the divertor design. A convective–diffusive anomalous transport model is implemented in the code to reproduce the radial plasma profiles anticipated for the device at the outer mid-plane. Levels of exhaust power from the core are varied, and a 0.5% fixed fraction neon impurity seeding is included in investigations to find steady-state detached solutions. Initial studies employing a super-X divertor configuration have shown that a stable detached divertor-operational window exists for reactor exhaust powers in the range of 72–92 MW, which is close to projected levels of exhaust power of 105 MW across the separatrix anticipated for ARC. Under these conditions, the divertor heat flux is fully dissipated by radiation and radial losses to the sidewalls of the divertor channel, leading to acceptable levels of power flux density on plasma-facing components. Simulations are being extended to study the performance of the X-point target geometry and to explore the sensitivity to the physics model and input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Monte-Carlo model of pitch-angle scattering in solar cosmic ray events

Monte Carlo simulations of the propagation of solar cosmic-rays in interplanetary space are reported, including the effects of pitch-angle scattering and adiabatic focusing. Time intensity profiles agree well with the corresponding spatial diffusion approximation for models close to those used by Gombosi and Owens (1980). Monte Carlo simulations of the same model problem yield results that disagree with those of Gombosi and Owens, and the most likely possible source of error in the Monte Carlo simulation is considered to be in the region near the 90 deg pitch angle where the finite step size introduces certain inconsistencies.

Palmer, I. D.↗
DOE OSTI2021

Single molecule characterization of anomalous transport in a thin, anisotropic film

Here, the diffusion of small, charged molecules incorporated in an anisotropic polyelectrolyte multilayers (PEMs) was tracked in three dimensions by combining single-molecule fluorescence localization (to characterize lateral diffusion) with Förster resonance energy transfer (FRET) between diffusing molecules and the supporting surface (to measure diffusion in the surface-normal direction). Analysis of the surface-normal diffusion required model-based statistical analysis to account for the inherently noisy FRET signal. Combining these distinct single-molecule methods, which are inherently sensitive to different length scales, permitted simultaneous characterization of severely anisotropic diffusion, which was more than three orders of magnitude slower in the surface-normal direction. We hypothesize that the anomalously slow surface-normal diffusion was related to the periodic distribution of charge in the PEM, which created electrostatic barriers. The motion was strongly subdiffusive, with anomalous temporal scaling exponents in lateral and normal directions, suggesting a connection to the transient, random fractal conformation of polymer chains in the film’s matrix.

36 MATERIALS SCIENCE↗
NASA NTRS1997

A Model Study of Zonal Forcing in the Equatorial Stratosphere by Convectively Induced Gravity Waves

A two-dimensional cloud-resolving model is used to examine the possible role of gravity waves generated by a simulated tropical squall line in forcing the quasi-biennial oscillation (QBO) of the zonal winds in the equatorial stratosphere. A simulation with constant background stratospheric winds is compared to simulations with background winds characteristic of the westerly and easterly QBO phases, respectively. In all three cases a broad spectrum of both eastward and westward propagating gravity waves is excited. In the constant background wind case the vertical momentum flux is nearly constant with height in the stratosphere, after correction for waves leaving the model domain. In the easterly and westerly shear cases, however, westward and eastward propagating waves, respectively, are strongly damped as they approach their critical levels, owing to the strongly scale-dependent vertical diffusion in the model. The profiles of zonal forcing induced by this wave damping are similar to profiles given by critical level absorption, but displaced slightly downward. The magnitude of the zonal forcing is of order 5 m/s/day. It is estimated that if 2% of the area of the Tropics were occupied by storms of similar magnitude, mesoscale gravity waves could provide nearly 1/4 of the zonal forcing required for the QBO.

Alexander, M. J.↗
DOE OSTI2018

Conditions for translation and scaling invariance of the neutron diffusion equation

Lie group methods are applied to the time-dependent, monoenergetic neutron diffusion equation in materials with spatial and time dependence. To accomplish this objective, the underlying 2nd order partial differential equation (PDE) is recast as an exterior differential system so as to leverage the isovector symmetry analysis approach. Some of the advantages of this method as compared to traditional symmetry analysis approaches are revealed through its use in the context of a 2nd order PDE. In this context, various material properties appearing in the mathematical model (e.g., a diffusion coefficient and macroscopic cross section data) are left as arbitrary functions of space and time. The symmetry analysis that follows is restricted to a search for translation and scaling symmetries; consequently the Lie derivative yields specific material conditions that must be satisfied in order to maintain the presence of these important similarity transformations. The principal outcome of this work is thus the determination of analytic material property functions that enable the presence of various translation and scaling symmetries within the time-dependent, monoenergetic neutron diffusion equation. The results of this exercise encapsulate and generalize many existing results already appearing in the literature. While the results contained in this work are primarily useful as phenomenological guides pertaining to the symmetry behavior of the neutron diffusion equation under certain assumptions, they may eventually be useful in the construction of exact solutions to the underlying mathematical model. Furthermore, the results of this work are also useful as a starting point or framework for future symmetry analysis studies pertaining to the neutron transport equation and its many surrogates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Along the Channel Gradients Impact on the Spatioactivity of Gas Diffusion Electrodes at High Conversions during CO 2 Electroreduction

Results of a 2-D transport model for a gas diffusion electrode performing CO 2 reduction to CO with a flowing catholyte are presented, including the concentration gradients along the flow cell, spatial distribution of the current density and local pH in the catalyst layer. The model predicts that both the concentration of CO 2 and the buffer electrolyte gradually diminish along the channels for a parallel flow of gas and electrolyte as a result of electrochemical conversion and nonelectrochemical consumption. At high single-pass conversions, significant concentration gradients exist along the flow channels leading to large local variations in the current density (>150 mA/cm 2 ), which becomes prominent when compared to ohmic losses. In addition, concentration overpotentials change dramatically with CO 2 flow rate, which results in significant differences in outlet concentrations at high conversions. The outlet concentration of CO attains a maximum of 80% along with 5% CO 2 and 15% H 2 , although the maximum single-pass conversion is limited to below 60% due to homogeneous consumption by the electrolyte. Fundamental and practical implications of our findings on electrochemical CO 2 reduction are discussed with a focus on the trade-off between high current density operation and high single-pass conversion efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗