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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Pickup protons and water ions at Comet Halley - Comparisons with Giotto observations

The cometary ion pickup process along the sun-comet line at Comet Halley is investigated using a quasi-linear diffusion model including both pitch angle and energy diffusion, adiabatic compression, and convective motion with the solar wind flow. The model results are compared with energetic ion distributions observed by instruments on board the Giotto spacecraft. The observed power spectrum index of magnetic turbulence (gamma) is 2-2.5. The present simulation shows that when gamma was 2, the calculated proton distributions were much more isotropic than the observed ones. The numerical solutions of the quasi-linear diffusion equations show that the isotropization of the pickup ion distribution, particularly at the pickup velocity, is not complete even close to the bow shock. Given the observed turbulence level, quasi-linear theory yields pickup ion energy distributions that agree with the observed ones quite well and easily produces energetic ions with energies up to hundreds of keV.

Ye, G.↗
DOE OSTI2023

Tree-based solvers for adaptive mesh refinement code FLASH – III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources

Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here, we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, TREERAY/RADPRESSURE, is an extension to the novel radiative transfer method TREERAY implemented in the grid-based Magneto-Hydrodynamics code FLASH. In TREERAY/RADPRESSURE, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust, and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 M ⊙ . We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust–gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time-scale in this work.

(ISM:) HII regions↗
NASA NTRS2009

BRDF Calibration of Sintered PTFE in the SWIR

Satellite instruments operating in the reflective solar wavelength region often require accurate and precise determination of the Bidirectional Reflectance Distribution Function (BRDF) of laboratory-based diffusers used in their pre-flight calibrations and ground-based support of on-orbit remote sensing instruments. The Diffuser Calibration Facility at NASA's Goddard Space Flight Center is a secondary diffuser calibration standard after NEST for over two decades, providing numerous NASA projects with BRDF data in the UV, Visible and the NIR spectral regions. Currently the Diffuser Calibration Facility extended the covered spectral range from 900 nm up to 1.7 microns. The measurements were made using the existing scatterometer by replacing the Si photodiode based receiver with an InGaAs-based one. The BRDF data was recorded at normal incidence and scatter zenith angles from 10 to 60 deg. Tunable coherent light source was setup. Broadband light source application is under development. Gray-scale sintered PTFE samples were used at these first trials, illuminated with P and S polarized incident light. The results are discussed and compared to empirically generated BRDF data from simple model based on 8 deg directional/hemispherical measurements.

Georgiev, Georgi T.↗
NASA NTRS2004

Quantitative Rainbow Schlieren Deflectometry as a Temperature Diagnostic for Spherical Flames

Numerical analysis and experimental results are presented to define a method for quantitatively measuring the temperature distribution of a spherical diffusion flame using Rainbow Schlieren Deflectometry in microgravity. First, a numerical analysis is completed to show the method can suitably determine temperature in the presence of spatially varying species composition. Also, a numerical forward-backward inversion calculation is presented to illustrate the types of calculations and deflections to be encountered. Lastly, a normal gravity demonstration of temperature measurement in an axisymmetric laminar, diffusion flame using Rainbow Schlieren deflectometry is presented. The method employed in this paper illustrates the necessary steps for the preliminary design of a Schlieren system. The largest deflections for the normal gravity flame considered in this paper are 7.4 x 10(-4) radians which can be accurately measured with 2 meter focal length collimating and decollimating optics. The experimental uncertainty of deflection is less than 5 x 10(-5) radians.

Feikema, Douglas A.↗
NASA NTRS1995

A Hybrid Parallel Preconditioning Algorithm For CFD

A new hybrid preconditioning algorithm will be presented which combines the favorable attributes of incomplete lower-upper (ILU) factorization with the favorable attributes of the approximate inverse method recently advocated by numerous researchers. The quality of the preconditioner is adjustable and can be increased at the cost of additional computation while at the same time the storage required is roughly constant and approximately equal to the storage required for the original matrix. In addition, the preconditioning algorithm suggests an efficient and natural parallel implementation with reduced communication. Sample calculations will be presented for the numerical solution of multi-dimensional advection-diffusion equations. The matrix solver has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The full paper will show numerous examples in CFD to demonstrate the efficiency and robustness of the method.

Barth,Timothy J.↗
NASA NTRS1991

Combined energy and pitch angle diffusion of pickup ions at Comet Halley

It is well known that cometary pickup ions, e.g., H2O(+), OH(+), O(+), CO(+), H(+), initially form a ring-beam distribution in the solar wind reference frame, which is highly unstable to the growth of MHD waves (such as ion-cyclotron waves). The low-frequency magnetic fluctuations (or waves), which were observed upstream of Comet Halley, cannot only pitch-angle-scatter the pickup ions so that the distribution becomes at least partially isotropized, but also stochastically accelerate the ions, resulting in the energetic ion populations observed in the vicinity of Comet Halley. Here, numerical solutions of the quasi-linear diffusion equation were used to investigate the cometary ion pickup process at Comet Halley. Both pitch angle and energy diffusion are taken into account. Many quasi-linear models of cometary pickup ions exist which involve one type of diffusion or the other but not both types at once. It is found that the pitch angle scattering occurs faster than the energy diffusion, as expected. Moreover, the results demonstrate that the distribution of accelerated energetic ions is more isotropic than that of ions which have just been picked up. In fact, the ion distribution function on the initial pickup shell is quite anisotropic, even close to the Comet-Halley bow shock.

Ye, Gang↗
DOE OSTI2022

Spatiotemporal Route to Understanding Metal Halide Perovskitoid Crystallization

Here, a spatiotemporal experimental route is reported for the antisolvent vapor diffusion crystal growth of metal halide perovskitoids. A computational analysis combining automated image capture and diffusion modeling enables the determination of the critical concentrations required for nucleation and crystal growth from a single experiment. Five different solvent systems and ten distinct organic ammonium iodide salts were investigated with lead iodide, from which nine previously unreported compounds were discovered. Automated image capture of the mother liquor and antisolvent vials was used to determine changes in solution meniscus positions and detect the nucleation event location. Matching the observations to a numerical solution of Fick's second law diffusion model enables the calculation of reactant, solvent, and antisolvent concentrations at both the time and position of the first stable nucleation and crystal growth. A machine learning model was trained on the resulting data, and it reveals solvent- and amine-specific crystallization tendencies. Solvent systems that interact more weakly with dissolved lead species promote crystallization, while those with stronger interactions can prevent crystallization through increased solubilities. Organic amines that interact more strongly with inorganic components and exhibit greater rigidity are more likely to be incorporated into crystalline products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Diffuse energy transport and coda-wave interferometry for resonant transmission between reverberant structures

Approximate analytic and numerical theories are developed with which to model, and compare with laboratory measurements, the diffuse ultrasonic field transmitted from one elastic body to another through a resonant coupling. Particular focus is on the sensitivity of the transmitted field to perturbations in the natural frequency of the coupling, and on the manifestation of such perturbations as time delays in the second body. Immediate application is to measurements of slow dynamic elastic nonlinearity. It is found that the developed approximate theories do a good job of predicting the time dependence of mean square transmitted signal and signal spectra. They also predict the time delays, both their erratic character and their means. Analytic dependence of these time shifts on the fractional frequency perturbation df/f is derived.

42 ENGINEERING↗
DOE OSTI2025

Extended molecular eigenmodes treatment of dipole–dipole NMR relaxation in real fluids

Traditional models of NMR relaxation fail to account for the complex, multi-exponential behavior of the autocorrelation function in realistic systems characterized by soft-interactions and molecules that are chemically and physically complex. Here, in this study, we describe the relative diffusion of the spin dipoles by means of a Fokker–Planck equation that includes an interaction potential of mean force to account for the response of the physical/chemical environment around the dipoles. By numerically solving the Fokker–Planck equation for the diffusion propagator, we estimate dipole–dipole NMR relaxation for like- and unlike-spin systems via its eigenmode solution. We test the model against molecular simulations of diffusing dipoles with harmonic potentials and also validate using experimental longitudinal relaxation data from real systems, including Gd(III)–aqua and Gd(III)–DO3A–butrol complexes, the latter being an important MRI contrast agent. Using this novel approach, we predict both the inner- and outer-shell contributions to the relaxivity rates with excellent accuracy at frequencies relevant to MRI. We also show that, under the appropriate assumptions, our framework naturally recovers the Bloembergen–Purcell–Pound, the Solomon–Bloembergen–Morgan, and the Hwang–Freed models. Our implementation is general and publicly available for application to a broad range of systems.

Pinheiro dos Santos, Thiago J. [Rice Univ., Housto↗
DOE OSTI2023

Code for the manuscript titled "Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces"

We would like to disclose two scripts, written in Jupyter notebook, in which we implement the "Blackout Diffusion Process" described in the manuscript "Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces" (LA-UR-23-20509), to be submitted to the International Conference of Machine Learning (ICML). The abstract of the manuscript is append below. == Typical generative diffusion models rely on a Gaussian diffusion process for training the backward transformations, which can then be used generate samples from Gaussian noise. However, real world data often takes place in discrete-state spaces, which occur in many scientific applications. Here, we develop a theoretical formulation for arbitrary discrete-state Markov processes in the forward diffusion process. We relate the theory to the existing continuous-state Gaussian diffusion and identify the corresponding reverse-time stochastic process and score function in the continuous-time setting, and the reverse-time mapping in the discrete-time setting. As an example of this framework, we introduce "Blackout Diffusion", which learns to produce samples from an empty image instead of from noise. Numerical experiments on the CIFAR-10 dataset confirm the feasibility of generative diffusion modeling in a discrete space. Generalizing from specific (Gaussian) forward processes to a more general framework also sheds light on how to interpret generative diffusion models and their mathematical structure, which we comment on.

Lin, Yen Ting↗
NASA NTRS1984

Evaluation of effective eddy diffusive coefficients using radar observations of turbulence in the stratosphere

Radar observations show that thin, persistent layers of turbulence occur sporadically in the troposphere and stratosphere. Two probabilistic approaches are used to show that the vertical eddy diffusivity due to such layers is of the order of 0.2-0.3 sq m/sec in the lower stratosphere. An actual realization of turbulent layers, derived from the radar observations at Arecibo, is used in a numerical approach to obtain a profile of eddy diffusivity. It is suggested that turbulence plays a significant role in the vertical transport of trace constituents in the stratosphere.

Woodman, R. F.↗
NASA NTRS1975

On upstream blocking in a viscous diffusive stratified flow

The effect of diffusion of specie upon the flow about a transverse flat plate moving horizontally in a viscous stratified medium is considered. Asymptotic expansions are used to define a parameter regime where a viscous-diffusive-buoyancy balance is dominant. The solution, expressed in terms of an inverse Fourier transform, is numerically integrated. The results show that, as in the non-diffusive problem, a region of closed streamlines exists ahead of the body. However, unlike the case where diffusion is neglected, the density field within this recirculating region is uniquely determined and found to be statically stable. It is also found that varying the relative amount of diffusion affects not only the density distribution, but the velocity profile as well, indicating a strong coupling between the vorticity and specie equation.

Koop, C. G.↗
NASA NTRS1976

Prediction of the diffuse far infrared flux from the galactic plane

A basic model and simple numerical relations useful for future far infrared studies of the galaxy are presented. Making use of recent CO and other galactic surveys, the diffuse far infrared flux distribution from the galactic plane as a function of galactic longitude alternate theta for 4 deg or = alternate theta or = 90 and the far infrared emissivity as a function of galactocentric distance is predicted. Future measurements of the galactic far infrared flux would yield valuable information on the physical properties and distribution of dust and molecular clouds in the galaxy, particularly the inner region.

Fazio, G. G.↗
NASA NTRS2001

Nonlinearity Effects of Lateral Density Diffusion Coefficient on Gain-Guided VCSEL Performance

Electron and hole diffusions in the plane of semiconductor quantum wells play an important part in the static and dynamic operations of semiconductor lasers. In this paper, we apply a hydrodynamic model developed from the semiconductor Bloch equations to numerically study the effects of nonlinearity in the diffusion coefficient on single mode operation and direct modulation of a gain-guided InGaAs/GaAs multiple quantum well laser, operating not too far from threshold. We found that a small diffusion coefficient is advantageous for lowering the threshold current and increasing the modulation bandwidth. Most importantly, the effects of nonlinearity in the coefficient can be approximately reproduced by replacing the coefficient with an effective constant diffusion coefficient, which corresponds roughly to the half height density of the density distribution.

Li, Jian-Zhong↗
DOE OSTI2023

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

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

Galactic cosmic ray modulation from 1965-1970.

Numerical solutions of the cosmic-ray equation of transport within the solar cavity and including the effects of diffusion, convection, and energy losses due to adiabatic deceleration, have been used to reproduce the modulation of galactic electrons, protons, and helium nuclei observed during the period from 1965 to 1970. Kinetic energies between 10 and 10,000 MeV/nucleon are considered. Computed and observed spectra are given for the years 1965, 1968, 1969, and 1970 together with the diffusion coefficients. These diffusion coefficients are assumed to be of separable form in rigidity and radial dependence, and are consistent with the available magnetic-field power spectra. The force-field solutions are given for these diffusion coefficients and galactic spectra and are compared with the numerical solutions. It is shown that the energy losses and convection lead to near-earth nuclei spectra at kinetic energies less than or equal to 100 MeV/nucleon in which the differential intensity is proportional to the kinetic energy with little dependence on the form of the galactic spectrum. This dependence is in agreement with the observed spectra of all species of atomic nuclei and it is argued that this provides strong observational evidence for the presence of energy losses in the propagation process, and for the exclusion of low-energy galactic nuclei from near earth.

Urch, I. H.↗
NASA NTRS1991

Turbulence and stall in plane diffusers - Computational study

The effect of free-stream turbulence level on stall conditions in plane diffusers is investigated with a finite-volume numerical method which utilizes a k-epsilon turbulence model or a differential Reynolds stress model for closure. Computed results show qualitative agreement with measured data. The ability to eliminate stall and thereby increase pressure recovery in diffusers at large included angles by utilizing high free-stream turbulence is correctly predicted.

Demuren, A. O.↗
NASA NTRS1991

A detailed numerical investigation of Burke-Schumann gaseous and spray flames

The classical Burke-Schumann gaseous diffusion flame, and the effect of introducing fuel in the form of liquid sprays are investigated numerically. The time-dependent two-dimensional, axisymmetric conservation equations are solved for the gas phase. A Lagrangian approach is used for the dispersed phase. The chemical reactions are modeled through the one-step global reaction scheme. The numerical results show that the computed flame structure is significantly different from that given by the analytical solution. The computed flame is about 15 percent longer and 5 percent narrower than the classical Burke-Schumann flame. The predicted flame shape is, however, in better agreement with experimental observations than the Burke-Schumann flame. The effects of introducing liquid sprays is found to produce thinner and longer flames than its counterpart gaseous flame. A detailed comparison of the structures of gaseous and spray flames is made, and the differences are discussed.

Mawid, M. A.↗