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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s), and passive scalars K(sub c). First, we took K(sub c) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub c) = K(sub h). Indeed, the turbulence model yelds a salt diffusivity that depends on R(sub i) and R(sub rho) = ((beta)(delta)S/(alpha)(delta)T/(delta)z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗
NASA NTRS1999

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in General Circulation Models (GCM'S). The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s). and passive scalars K(sub c). First, we took K(sub s) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derived using the latest model, the KPP model by the NCAR group. Presently, we are running the same Geophysical Fluid Dynamics Laboratory (GFDL) code relaxing the assumptive K(sub s) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R rho (= Beta delta S/delta z/alpha delta T/delta z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗
NASA NTRS1999

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of FLuids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusivities K for moment K(sub m), temperature K(sub h), salt K(sub g) and passive scalars K(sub c). First, we took K(sub g) = K(sub h) and tested the model using the Geophysical Fluid Dynamics Laboratory (GFDL) ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub g) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R Rho (= Beta Sigma S/Sigma z/Alpha Sigma T/Sigma z) in such a way that K(sub g) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, Vittorio↗
DOE OSTI2022

Diffusion of cesium in oxidized and unoxidized IG-110 nuclear graphite

Time-release diffusion measurements of cesium have been conducted over the temperature range 1073 K – 1973 K on oxidized and unoxidized IG-110 graphite. Four cesium concentrations were tested to investigate the concentration dependence of the diffusion coefficient. Two levels of oxidation were tested and compared to unoxidized concentration-matched sets to explore the effects of graphite oxidation. Here, the results demonstrate that cesium diffusion coefficient in unoxidized IG-110 graphite is independent of concentration within the range 34 – 163 µg Cs/g graphite . Above this, the effective cesium diffusion coefficient changes with concentration. The diffusion coefficient was increased by a factor of 2–12 in the oxidized set with 7.8% mass loss. These results can be used to aid predictive modeling of cesium diffusion in HTGR cores.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Delayed Gamma Signature Accounting for Field Line Transport of Debris Ions after a Nuclear Explosion

Delayed gammas (DGs) are produced by the decay of radioactive fission fragments after a nuclear explosion. Currently, the time and energy profile of the delayed gamma signal is propagated in DIORAMA to the GEO satellite sensors by assuming a point source at the burst location for the production of the delayed gammas. However, in reality the production of delayed gammas is spread over some spatial extent. First, the debris ions are confined within an expanding region by the magnetic field and/or surrounding plasma. This expanding confinement bubble eventually collapses, at which point the debris ions can be trapped on magnetic field lines and transported away from the original burst location. In this work the effects of debris ion transport on the associated time profile of delayed gammas at a GEO satellite is presented. The first step was the development of the Delayed Gamma Source Model (DGSM), which is based on the analytical Electron Source Model (ESM) code developed at LANL that models the spatial distribution of energetic electrons produced after a high-altitude nuclear explosion (HANE) for use as a source term in LANL artificial belt diffusion code. This model creates a spatial grid in latitude, longitude, and altitude of the distribution of radioactive debris ions. The output of this code is combined with a separate analysis code also described, which produces the DG time profile signature at a GEO sensor faceplate when accounting for the spatial spread of delayed gamma production during field line transport. Some simple case studies are presented to show how the DG signature at GEO calculated with this higher-fidelity model and a simple point-source model, as is currently assumed in DIORAMA, compare.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1986

Solvability condition for needle crystals at large undercooling in a nonlocal model of solidification

It is explicitly shown that, in a realistic model of diffusion-controlled dendritic solidification, Ivantsov's continuous family of steady-state needle crystals is destroyed by the addition of surface tension. The starting point is in the exact integro-differential equation for the one-sided model, in two dimensions, in a moving frame of reference. In the limit of large undercooling, where the range of the diffusion field is much smaller than the radius of curvature of the tip of the needle, this problem is reduced to a linear, inhomogeneous differential equation of infinite order. A solvability condition for this equation is derived and it is shown that solutions cease to exist for arbitrarily small but finite isotropic surface tension.

Caroli, B.↗
NASA NTRS1986

Third-moment closure of turbulence for predictions of separating and reattaching shear flows: A study of Reynolds-stress closure model

A numerical study of computations in backward-facing steps with flow separation and reattachment, using the Reynolds stress closure is presented. The highlight of this study is the improvement of the Reynold-stress model (RSM) by modifying the diffusive transport of the Reynolds stresses through the formulation, solution and subsequent incorporation of the transport equations of the third moments, bar-u(i)u(j)u(k), into the turbulence model. The diffusive transport of the Reynolds stresses, represented by the gradients of the third moments, attains greater significance in recirculating flows. The third moments evaluated by the development and solution of the complete transport equations are superior to those obtained by existing algebraic correlations. A low-Reynolds number model for the transport equations of the third moments is developed and considerable improvement in the near-wall profiles of the third moments is observed. The values of the empirical constants utilized in the development of the model are recommended. The Reynolds-stress closure is consolidated by incorporating the equations of k and e, containing the modified diffusion coefficients, and the transport equations of the third moments into the Reynolds stress equations. Computational results obtained by the original k-e model, the original RSM and the consolidated and modified RSM are compared with experimental data. Overall improvement in the predictions is seen by consolidation of the RMS and a marked improvement in the profiles of bar-u(i)u(j)u(k) is obtained around the reattachment region.

Amano, R. S.↗
DOE OSTI2021

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Comprehensive modeling of turbulent particulate flows using Eulerian and Lagrangian schemes

This paper addresses turbulent particle dispersion and modulation effects in dilute gas-particle turbulent flows using the Eulerian and Lagrangian modeling approaches. Gradient diffusion approximations are employed in the Eulerian formulation, while a stochastic procedure is utilized to simulate turbulent dispersion in the Lagrangian formulation. The k-epsilon turbulence model is used to characterize the time and length scales of the continuous phase turbulence. For the particle size and loading considered, the turbulence transport equations must be modified to account for the modulation effects. Models are proposed for both Eulerian and Lagrangian schemes. Comparisons and predictions are made in fully developed gas-solid pipe flow and confined coaxial jets laden with particles. For the monodispersed system investigated, Eulerian approach is less expensive and gives more consistent results than the Lagrangian approach. The Lagrangian technique should be further developed to eliminate current inherent inconsistencies especially with regard to the symmetry boundary condition.

Adeniji-Fashola, A.↗
NASA NTRS1995

Structure and dynamics of diffusion flames in microgravity

The objectives of this project are to gain insight into diffusion flames by modeling various configurations related to ongoing experimental investigations in the microgravity combustion science program. The emphasis of the work is to understand the structure and dynamics of diffusion flames. Improving our fundamental understanding of diffusion flames is most relevant to issues related to fire safety and fire prevention because most fires consist of diffusion flames.

Matalon, Moshe↗
DOE OSTI2024

Preserving Tracer Correlations in Moment-Based Atmospheric Transport Models

A linear non-diffusive algorithm for advective transport is developed that greatly improves the detail at which aerosols and clouds can be represented in atmospheric models. Linear advection schemes preserve tracer correlations but the most basic linear scheme is rarely used by atmospheric modelers on account of its excessive numerical diffusion. Higher-order schemes are in widespread use, but these present new problems as nonlinear adjustments are required to avoid occurrences of negative concentrations, spurious oscillations, and other non-physical effects. Generally successful at reducing numerical diffusion during the advection of individual tracers, for example, particle number or mass, the higher-order schemes fail to preserve even the simplest of correlations between interrelated tracers. As a result, important attributes of aerosol and cloud populations including radial moments of particle size distributions, molecular precursors related through chemical equilibria, aerosol mixing state, and distribution of cloud phase are poorly represented. We introduce a new transport scheme, minVAR, that is both non-diffusive and preservative of tracer correlations, thereby combining the best features of the basic and higher-order schemes while enabling new features such as the tracking of sub-grid information at arbitrarily fine scales with high computational efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Simulations of divertor heat flux width using transport code with cross-field drifts under the BOUT++ framework

The fluid transport code [trans-electric field (Er) module] under the BOUT++ framework has been used to simulatedivertor heat flux width and boundary Er with all drifts and the sheath potential in the scrape-off layer. The calculated steady state radial Er in the pedestal region has been compared with that of experimental measurements from the Alcator C-Mod tokamak. The magnitude and shape of Er are similar to those of the experimental data. In order to understand the relative role of cross-field drifts vs turbulent transport in setting the heat flux width, four C-Mod enhanced Dα H-mode discharges with a lower single null divertor configuration should be simulated. BOUT++ transport simulations with cross-field drifts included yield similar heat flux width λq to that of experimental measurements (within a factor of 2) from both the probe and the surface thermocouple diagnostics and show a similar trend with plasma current to that of the Eich experimental scaling. The simulations show that both drifts and turbulent transport compete to determine the heat flux width. The magnetic drifts play a dominant role in setting the divertor heat-flux width, while the E × B drift decreases the heat flux width by 10%–25%, leading to improved agreement with the experiment relative to Goldston’s model. A turbulence diffusivity scan (χ) identifies two distinct regimes: a drift dominant regime when χ is small and a turbulence dominant regime when χ is large. The Goldston heuristic drift model yields a lower limit of the width λq.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Rocket exhaust effluent modeling for tropospheric air quality and environmental assessments

The various techniques for diffusion predictions to support air quality predictions and environmental assessments for aerospace applications are discussed in terms of limitations imposed by atmospheric data. This affords an introduction to the rationale behind the selection of the National Aeronautics and Space Administration (NASA)/Marshall Space Flight Center (MSFC) Rocket Exhaust Effluent Diffusion (REED) program. The models utilized in the NASA/MSFC REED program are explained. This program is then evaluated in terms of some results from a joint MSFC/Langley Research Center/Kennedy Space Center Titan Exhaust Effluent Prediction and Monitoring Program.

Stephens, J. B.↗
NASA NTRS1978

Compressor and fan wake characteristics

A triaxial probe and a rotating conventional probe, mounted on a traverse gear operated by two step motors were used to measure the mean velocities and turbulence quantities across a rotor wake at various radial locations and downstream stations. The data obtained was used in an analytical model developed to study how rotor flow and blade parameters and turbulence properties such as energy, velocity correlations, and length scale affect the rotor wake characteristics and its diffusion properties. The model, includes three dimensional attributes, can be used in predicting the discrete as well as broadband noise generated in a fan rotor, as well as in evaluating the aerodynamic losses, efficiency and optimum spacing between a rotor and stator in turbomachinery.

Reynolds, B.↗
NASA NTRS1996

Hard Gamma Ray Emission from the Starburst Galaxy NGC 253

We have completed the study to search for hard gamma ray emission from the starburst galaxy NGC 253. Since supernovae are thought to provide the hard gamma ray emission from the Milky Way, starburst galaxies, with their extraordinarily high supernova rates, are prime targets to search for hard gamma ray emission. We conducted a careful search for hard gamma ray emission from NGC 253 using the archival data from the EGRET experiment aboard the CGRO. Because this starburst galaxy happens to lie near the South Galactic Pole, the Galactic gamma ray background is minimal. We found no significant hard gamma ray signal toward NGC 253, although a marginal signal of about 1.5 sigma was found. Because of the low Galactic background, we obtained a very sensitive upper limit to the emission of greater than 100 MeV gamma-rays of 8 x 10(exp -8) photons/sq cm s. Since we expected to detect hard gamma ray emission, we investigated the theory of gamma ray production in a dense molecular medium. We used a leaky-box model to simulate diffusive transport in a starburst region. Since starburst galaxies have high infrared radiation fields, we included the effects of self-Compton scattering, which are usually ignored. By modelling the expected gamma-ray and synchrotron spectra from NGC 253, we find that roughly 5 - 15% of the energy from supernovae is transferred to cosmic rays in the starburst. This result is consistent with supernova acceleration models, and is somewhat larger than the value derived for the Galaxy (3 - 10%). Our calculations match the EGRET and radio data very well with a supernova rate of 0.08/ yr, a magnetic field B approx. greater than 5 x 10(exp -5) G, a density n approx. less than 100/sq cm, a photon density U(sub ph) approx. 200 eV/sq cm, and an escape time scale tau(sub 0) approx. less than 10 Myr. The models also suggest that NGC 253 should be detectable with only a factor of 2 - 3 improvement in sensitivity. Our results are consistent with the standard picture of gamma-ray acceleration by supernovae.

Jackson, James M.↗
DOE OSTI2017

Sorption Behavior of Cationic Radionuclides on Biotite Flakes at pH 5 - 17326

Biotite contained in crystalline rock (granite) is a notable mineral because of its high sorption effects for cationic radionuclides. To date, studies on the sorption effects of rocks have focused on the powder form of minerals. However, the form of a mineral will seriously affect its sorption of radionuclides. This study examined the sorption of Cs, Sr, and Eu by biotite in the flake form. Biotite flakes (5 mm x 6 mm) were prepared for the sorption experiments. The initial concentrations of the tracer ions were set to 1 mM for Sr and Eu and 2 mM for Cs. The sorption experiments for Cs, Sr, and Eu were conducted by mixing 0.6 g of biotite flakes and 6 mL of trace solution at pH 5 for 7 days. As a result, the sorption of tracer ions gradually increased with continuous diffusion into the biotite flakes. The diffusion of the tracer ions into the layer structure of biotite was also observed by secondary ion-microprobe mass spectrometry (SIMS) analysis. The apparent diffusion coefficient of Eu into the biotite flakes was estimated to be 3.0 x 10-12 m{sup 2}/s by applying the experimental data to a theoretical model. This estimated diffusion coefficient is smaller than that for matrix diffusion in granite. Therefore, these diffusion and sorption behaviors may be important factors to consider to more reliably estimate the migration rates of radionuclides in host rock. (authors)

58 GEOSCIENCES↗
DOE OSTI2022

On the parameter combinations that matter and on those that do not: data-driven studies of parameter (non)identifiability

Abstract We present a data-driven approach to characterizing nonidentifiability of a model’s parameters and illustrate it through dynamic as well as steady kinetic models. By employing Diffusion Maps and their extensions, we discover the minimal combinations of parameters required to characterize the output behavior of a chemical system: a set of effective parameters for the model. Furthermore, we introduce and use a Conformal Autoencoder Neural Network technique, as well as a kernel-based Jointly Smooth Function technique, to disentangle the redundant parameter combinations that do not affect the output behavior from the ones that do. We discuss the interpretability of our data-driven effective parameters, and demonstrate the utility of the approach both for behavior prediction and parameter estimation. In the latter task, it becomes important to describe level sets in parameter space that are consistent with a particular output behavior. We validate our approach on a model of multisite phosphorylation, where a reduced set of effective parameters (nonlinear combinations of the physical ones) has previously been established analytically.

Evangelou, Nikolaos↗
NASA NTRS1996

Studies using IMP, Voyager and Pioneer cosmic ray data to determine the size of the heliosphere

The purpose of this project was to use the cosmic ray data from the IMP, Voyager and Pioneer spacecraft in the heliosphere out to approximately 65 AU to estimate the size of the heliosphere. We used several techniques to develop a consistent picture of the size of the heliosphere. The first method used a response function approach which determined the intensity as a function of time by scaling the modulation effect as they move outward and eventually reach the boundary of the heliosphere. In this model the effects of transient cosmic ray disturbances is included. A second approach using the perturbation method in which drifts are considered as a perturbation to the standard diffusion-convection modulation models was not fully developed. In a third approach the location of the modulation boundary beyond the termination shock was estimated using observations of the intensity and radial gradients between Voyager 2 and Pioneer 10 along with new estimates of the interstellar intensity of more than 70 MeV galactic cosmic rays. Using this method we found that for 7 years, from 1983 to 1990, the modulation boundary remained constant at 83 +/- 5 AU. We infer from these studies that a modulation boundary can be estimated only by extrapolating the observed radial gradients when the solar magnetic field polarity is such that cosmic-ray particles are drifting in the heliosphere inward toward the Earth along the neutral sheet. The boundary distance is larger than the estimates of the location of the termination shock at 67 +/- 5 AU using the same method. Two other studies partially supported by this grant are attached. The first deals with the recovery period of the greater than 70 MeV cosmic rays in the outer heliosphere from 1992-1995. In the second paper we compare the rigidity dependence of the 11-year cosmic ray variation at the Earth in two cycles of opposite solar magnetic field polarity.

Lockwood, John A.↗