Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Numerical solution of the problem of flame propagation by the use of the random element method

A numerical, grid-free algorithm is presented for one-dimensional reaction-diffusion model of laminar flame propagation in premixed gases. It is based on the random element method we developed for the analysis of diffusional processes. The effect of combustion is taken into account by applying the principle of fractional steps to separate the process of diffusion, modeled by the random walk of computational elements, from the exothermic effects of chemical reaction, monitoring their strength. The validity of the algorithm is demonstrated by application to flame propagation problems for which exact solutions exist. The flame speed evaluated by its use oscillates around the exact value at a relatively small amplitude, while the temperature and species concentration profiles are self-correcting in their convergence to the exact solution. A satisfactory resolution is obtained by the use of quite a small number of computational elements which automatically adjust their distribution of fit sharp gradients.

Ghoniem, A. F.↗
DOE OSTI2016

Phase space effects on fast ion distribution function modeling in tokamaks

Integrated simulations of tokamak discharges typically rely on classical physics to model energetic particle (EP) dynamics. However, there are numerous cases in which energetic particles can suffer additional transport that is not classical in nature. Examples include transport by applied 3D magnetic perturbations and, more notably, by plasma instabilities. Focusing on the effects of instabilities, ad-hoc models can empirically reproduce increased transport, but the choice of transport coefficients is usually somehow arbitrary. New approaches based on physics-based reduced models are being developed to address those issues in a simplified way, while retaining a more correct treatment of resonant wave-particle interactions. The kick model implemented in the tokamak transport code TRANSP is an example of such reduced models. It includes modifications of the EP distribution by instabilities in real and velocity space, retaining correlations between transport in energy and space typical of resonant EP transport. The relevance of EP phase space modifications by instabilities is first discussed in terms of predicted fast ion distribution. Results are compared with those from a simple, ad-hoc diffusive model. It is then shown that the phase-space resolved model can also provide additional insight into important issues such as internal consistency of the simulations and mode stability through the analysis of the power exchanged between energetic particles and the instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

Estimates of ground-level mercury vapor concentrations from postulated launch-pad accidents of NASA space transportation systems

NASA has considered the use of mercury ion bombardment engines as extra-terrestrial propulsion systems where low thrust must be provided for long periods of time. However, the mercury fuel in these engines may pose a hazard if it is accidentally released to the troposphere. A description is presented of diffusion-model estimates of maximum 10-minute and 24-hour ground-level concentrations that could result from the accidental release of mercury from Titan/Centaur and Space Shuttle/Tug transportation systems. Three types of accidental release modes have been considered in the concentration calculations. For all three modes, the diffusion-model calculations were made for a variety of meteorological conditions in an effort to obtain estimates of the maximum credible concentrations under 'worst-case' conditions.

Dumbauld, R. K.↗
DOE OSTI2025

Observation of kinetic mix enhancement in thin-shell OMEGA implosions

Recent separated reactant experiments for thin-shell (6 µ⁢m) shock-driven implosions on OMEGA have demonstrated significant mix from a buried deuterated layer of the shell into the hot spot. Time resolved D 3 He-p reaction history data demonstrate a (50 ± 20)⁢ ps shift earlier in peak nuclear emission for separated reactant experiments relative to control, in contrast to past experimental data for thicker, 20 µ⁢m shells with no laser burn through that show a 75 ps delay due to the time required for hydrodynamic instabilities to develop. This contrast suggests that the mix mechanism was not hydrodynamic. Ion kinetic simulations utilizing fall line analyses show much closer agreement with mix yield and temperature than diffusion models, predicting a D 3 He-p mix yield of 1.7 × 10 9 as compared to the experimental value of 9.3⁢ (±2.1) × 10 8 . This is three orders of magnitude closer than the fall line analysis from a hydrodynamic simulation with an inline diffusive mix model, which suggests minimal mix and D 3 He-p yields of 5×10 5 . This makes kinetic mechanisms the only feasible explanation for the mix seen, demonstrating impact of a non-standard mix mechanism. An analytical model of this kinetic mix mechanism suggests that it can remain significant in situations when the shell expands significantly to low densities, and diffusive models predict negligible mix. Finally, kinetic mix will impact multiple types of high energy density, laser-driven fusion experiments including high-adiabat direct drive cryoexperiments, nuclear cross section experiments, and thin-shell polar direct drive experiments used to tune heat conduction models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An Ensemble Score Filter for Tracking High-Dimensional Nonlinear Dynamical Systems

We propose an ensemble score filter (EnSF) for solving high-dimensional nonlinear filtering problems with superior accuracy. A major drawback of existing filtering methods, e.g., particle filters or ensemble Kalman filters, is the low accuracy in handling high-dimensional and highly nonlinear problems. EnSF attacks this challenge by exploiting the score-based diffusion model, defined in a pseudo-temporal domain, to characterizing the evolution of the filtering density. EnSF stores the information of the recursively updated filtering density function in the score function, in stead of storing the information in a set of finite Monte Carlo samples (used in particle filters and ensemble Kalman filters). Unlike existing diffusion models that train neural networks to approximate the score function, we develop a training-free score estimation that uses mini-batch-based Monte Carlo estimator to directly approximate the score function at any pseudo-spatial-temporal location, which provides sufficient accuracy in solving high-dimensional nonlinear problems as well as saves tremendous amount of time spent on training neural networks. Another essential aspect of EnSF is its analytical update step, gradually incorporating data information into the score function, which is crucial in mitigating the degeneracy issue faced when dealing with very high-dimensional nonlinear filtering problems. High-dimensional Lorenz systems are used to demonstrate the performance of our method. EnSF provides surprisingly impressive performance in reliably tracking extremely high-dimensional Lorenz systems (up to 1,000,000 dimension) with highly nonlinear observation processes, which is a well-known challenging problem for existing filtering methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Water transport in reverse osmosis membranes is governed by pore flow, not a solution-diffusion mechanism

We performed nonequilibrium molecular dynamics (NEMD) simulations and solvent permeation experiments to unravel the mechanism of water transport in reverse osmosis (RO) membranes. The NEMD simulations reveal that water transport is driven by a pressure gradient within the membranes, not by a water concentration gradient, in marked contrast to the classic solution-diffusion model. We further show that water molecules travel as clusters through a network of pores that are transiently connected. Permeation experiments with water and organic solvents using polyamide and cellulose triacetate RO membranes showed that solvent permeance depends on the membrane pore size, kinetic diameter of solvent molecules, and solvent viscosity. This observation is not consistent with the solution-diffusion model, where permeance depends on the solvent solubility. Motivated by these observations, we demonstrate that the solution-friction model, in which transport is driven by a pressure gradient, can describe water and solvent transport in RO membranes.

42 ENGINEERING↗
DOE OSTI2021

Spatial profiles of photon chemical potential in near-field thermophotovoltaic cells

Emitted photons stemming from the radiative recombination of electron–hole pairs carry chemical potential in radiative energy converters. This luminescent effect can substantially alter the local net photogeneration in near-field thermophotovoltaic cells. Several assumptions involving the luminescent effect are commonly made in modeling photovoltaic devices; in particular, the photon chemical potential is assumed to be zero or a constant prescribed by the bias voltage. The significance of photon chemical potential depends upon the emitter temperature, the semiconductor properties, and the injection level. Hence, these assumptions are questionable in near-field thermophotovoltaic devices with nanoscale separation distances between the emitter and the cell. In the present work, an iterative solver that combines fluctuational electrodynamics with the drift-diffusion model is developed to tackle the coupled photon and charge transport problem, enabling the determination of the spatial profile of photon chemical potential beyond the detailed balance approach. Here, the difference between the results obtained by allowing the photon chemical potential to vary spatially and by assuming a constant value demonstrates the limitations of the conventional approaches. This study is critically important for performance evaluation of near-field thermophotovoltaic systems with nanoscale vacuum gaps.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS2001

Assessment of Remote Sensing Technologies for Location of Hydrogen and Helium Leaks

In Phase 1 of this project, a hierarchy of techniques for H2 and He leak location was developed. A total of twelve specific remote sensing techniques were evaluated; the results are summarized. A basic diffusion model was also developed to predict the concentration and distribution of H2 or He resulting from a leak. The objectives of Phase 2 of the project consisted of the following four tasks: Advance Rayleigh Doppler technique from TRL 1 to TRL 2; Plan to advance Rayleigh Doppler technique from TRL 2 to TRL 3; Advance researchers and resources for further advancement; Extend diffusion model.

Sellar, R. Glenn↗
NASA NTRS2010

Development and Application of Agglomerated Multigrid Methods for Complex Geometries

We report progress in the development of agglomerated multigrid techniques for fully un- structured grids in three dimensions, building upon two previous studies focused on efficiently solving a model diffusion equation. We demonstrate a robust fully-coarsened agglomerated multigrid technique for 3D complex geometries, incorporating the following key developments: consistent and stable coarse-grid discretizations, a hierarchical agglomeration scheme, and line-agglomeration/relaxation using prismatic-cell discretizations in the highly-stretched grid regions. A signi cant speed-up in computer time is demonstrated for a model diffusion problem, the Euler equations, and the Reynolds-averaged Navier-Stokes equations for 3D realistic complex geometries.

Nishikawa, Hiroaki↗
DOE OSTI2015

SOLAR INTERACTING PROTONS VERSUS INTERPLANETARY PROTONS IN THE CORE PLUS HALO MODEL OF DIFFUSIVE SHOCK ACCELERATION AND STOCHASTIC RE-ACCELERATION

With the first observations of solar γ-rays from the decay of pions, the relationship of protons producing ground level enhancements (GLEs) on the Earth to those of similar energies producing the γ-rays on the Sun has been debated. These two populations may be either independent and simply coincident in large flares, or they may be, in fact, the same population stemming from a single accelerating agent and jointly distributed at the Sun and also in space. Assuming the latter, we model a scenario in which particles are accelerated near the Sun in a shock wave with a fraction transported back to the solar surface to radiate, while the remainder is detected at Earth in the form of a GLE. Interplanetary ions versus ions interacting at the Sun are studied for a spherical shock wave propagating in a radial magnetic field through a highly turbulent radial ray (the acceleration core) and surrounding weakly turbulent sector in which the accelerated particles can propagate toward or away from the Sun. The model presented here accounts for both the first-order Fermi acceleration at the shock front and the second-order, stochastic re-acceleration by the turbulence enhanced behind the shock. We find that the re-acceleration is important in generating the γ-radiation and we also find that up to 10% of the particle population can find its way to the Sun as compared to particles escaping to the interplanetary space.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
NASA NTRS2008

Modeling the Diffuse Cloud-Top Optical Emissions from Ground and Cloud Flashes

A number of studies have indicated that the diffuse cloud-top optical emissions from intra-cloud (IC) lightning are brighter than that from normal negative cloud-to-ground (CG) lightning, and hence would be easier to detect from a space-based sensor. The primary reason provided to substantiate this claim has been that the IC is at a higher altitude within the cloud and therefore is less obscured by the cloud multiple scattering medium. CGs at lower altitudes embedded deep within the cloud are more obscured, so CG detection is thought to be more difficult. However, other authors claim that because the CG source current (and hence luminosity) is typically substantially larger than IC currents, the greater CG source luminosity is large enough to overcome the effects of multiple scattering. These investigators suggest that the diffuse cloud top emissions from CGs are brighter than from ICs, and hence are easier to detect from space. Still other investigators claim that the detection efficiency of CGs and ICs is about the same because modern detector sensitivity is good enough to "see" either flash type no matter which produces a brighter cloud top emission. To better assess which of these opinions should be accepted, we introduce an extension of a Boltzmann lightning radiative transfer model previously developed. It considers characteristics of the cloud (geometry, dimensions, scattering properties) and specific lightning channel properties (length, geometry, location, current, optical wave front propagation speed/direction). As such, it represents the most detailed modeling effort to date. At least in the few cases studied thus far, it was found that IC flashes appear brighter at cloud top than the lower altitude negative ground flashes, but additional model runs are to be examined before finalizing our general conclusions.

Solakiewicz, Richard↗
NASA NTRS1974

Mariner 6, 7, and 9 ultraviolet spectrometer experiment - Analysis of hydrogen Lyman alpha data

Four Lyman alpha airglow measurements of the limb and disk of Mars, made by ultraviolet spectrometers on Mariner 6 and 7 in 1969 and Mariner 9 in 1971, are analyzed to determine the amount and distribution of atomic hydrogen above 80 km. The variation of atomic hydrogen with altitude is calculated by using time-independent chemical diffusion models from 80 to 250 km, and an exospheric model is used above 250 km. By employing radiative transfer theory that includes effects of pure absorption and accounts for temperature variations in the atmosphere, a spherical model of the airglow Lyman alpha emission is used to produce theoretical intensities for comparison with the data. It is found that (1) the exospheric temperature and distribution in 1971 are consistent with those determined in 1969, (2) the vertical optical depth above 80 km was 2.2 in 1969 and 5 in 1971, and (3) the derived atomic hydrogen distribution from 80 to 250 km requires a source of atomic hydrogen above 80 km. Comparison of observed profiles with chemical diffusion models implies a large downward flow of atomic hydrogen at 80 km coupled with a large upward flow of molecular hydrogen.-

Anderson, D. E., Jr.↗
DOE OSTI2023

Multiscale modeling of solute diffusion in triblock copolymer membranes

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Furthermore, our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

36 MATERIALS SCIENCE↗
NASA NTRS1977

Radio-emitting electrons and cosmic ray confinement

Evidence from radio data obtained at frequencies of 10 to 8000 MHz for steepening of the observed background-radio-continuum spectra in the directions of the galactic Anticenter (A) and Halo Minimum (H) is used to deduce quantitative information on the variation of the magnetic field with distance from the galactic plane. The radio data are analyzed in the framework of the cosmic-ray diffusion model developed by Bulanov and Dogel (1975), and some inferences are drawn regarding the injection spectrum of cosmic-ray electrons as well as their propagation in the Galaxy. The results indicate that: (1) radio-spectrum steepening is centered around 200 MHz in the direction of H and around 330 MHz in the direction of A, implying a lower mean magnetic field toward H and supporting the existence of a radio halo; (2) an injection electron spectrum with a single power law up to the lowest energies cannot explain the radio observations satisfactorily in terms of the diffusion model; and (3) the observed spectral steepening can be satisfactorily understood as being due partly to the deviation from a power-law injection spectrum below a few GeV and partly to the first break arising from electron energy losses in approximately the same energy region.

Badhwar, G. D.↗
DOE OSTI2015

Tritium release from SS316 under vacuum condition

The plasma facing surface of the ITER vacuum vessel, partly made of low carbon austenitic stainless steel type 316L, will incorporate tritium during machine operation. In this paper the kinetics of tritium release from stainless steel type 316 into vacuum and into a noble gas stream are compared and modelled. Type 316 stainless steel specimens loaded with tritium either by exposure to 1.2 kPa HT at 573 K or submersion into liquid HTO at 298 K showed characteristic thin surface layers trapping tritium in concentrations far higher than those determined in the bulk. The evolution of the tritium depth profile in the bulk during heating under vacuum was non-discernible from that of tritium liberated into a stream of argon. Only the relative amount of the two released tritium-species, i.e. HT or HTO, was different. Temperature-dependent depth profiles could be predicted with a one-dimensional diffusion model. Diffusion coefficients derived from fitting of the tritium release into an evacuated vessel or a stream of argon were found to be (1.4 ± 1.0)*10{sup -7} and (1.3 ± 0.9)*10{sup -9} cm{sup 2}/s at 573 and 423 K, respectively. Polished surfaces on type SS316 stainless steel inhibit considerably the thermal release rate of tritium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Reduced model prediction of electron temperature profiles in microtearing-dominated National Spherical Torus eXperiment plasmas

A representative H-mode discharge from the National Spherical Torus eXperiment is studied in detail to utilize it as a basis for a time-evolving prediction of the electron temperature profile using an appropriate reduced transport model. The time evolution of characteristic plasma variables such as β{sub e}, ν{sub e}{sup ∗}, the MHD α parameter, and the gradient scale lengths of T{sub e}, T{sub i}, and n{sub e} were examined as a prelude to performing linear gyrokinetic calculations to determine the fastest growing micro instability at various times and locations throughout the discharge. The inferences from the parameter evolutions and the linear stability calculations were consistent. Early in the discharge, when β{sub e} and ν{sub e}{sup ∗} were relatively low, ballooning parity modes were dominant. As time progressed and both β{sub e} and ν{sub e}{sup ∗} increased, microtearing became the dominant low-k{sub θ} mode, especially in the outer half of the plasma. There are instances in time and radius, however, where other modes, at higher-k{sub θ}, may, in addition to microtearing, be important for driving electron transport. Given these results, the Rebut-Lallia-Watkins (RLW) electron thermal diffusivity model, which is based on microtearing-induced transport, was used to predict the time-evolving electron temperature across most of the profile. The results indicate that RLW does a good job of predicting T{sub e} for times and locations where microtearing was determined to be important, but not as well when microtearing was predicted to be stable or subdominant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

LANL Meteorology Program Self-Assessment 2025 Update

The LANL Meteorological (Met) Program has been subject to several external reviews over the past 19 years. The DOE Meteorological Coordinating Council (DMCC) conducted an initial Met Program Site Assist Visit (SAV) in August 2006 (DMCC 2006). A follow-up SAV in August 2015 assessed progress (DMCC 2015), and in June 2023, the DOE Meteorological Subcommittee (DMSC), successor to the DMCC, conducted a second follow-up SAV (DMSC 2023), in which the Met Program was evaluated relative to the following 8 high-level questions: • What is the state of the meteorological services provided to its customers? • What is the quality of meteorological data provided to its customers and is it adequate and available to meet all customer needs? • What is the quality of atmospheric transport and diffusion modeling provided to its customers and is it applicable to complex wind flow patterns at LANL? • Are the current and future meteorological service customers being serviced appropriately? • Are there adequate human resources to meet present and future program customer needs and are they being appropriately leveraged? • Are existing instrumentation, facilities, and systems adequate to meet present and future customer needs? • Are LANL meteorological services conducted in an efficient, cost-effective manner? • Is meteorological data used to ensure safety & health of LANL personnel? More specific evaluations were performed relative to 23 performance objectives extracted from the ANSI/ANS-3.11- 2024 national standard and 14 separate performance objectives associated with consequence assessment and atmospheric transport and diffusion modeling in the consequence assessment element of DOE G 151.1-1B. In 2023, the DMSC SAV Team also reviewed the status of each of the 18 remaining recommendations from its 2015 SAV. Based on this review, DMSC stated in its Exit Briefing that the LANL meteorological program has gotten much stronger and more robust since 2015 and now represents one of the better managed programs within the DOE complex.

54 ENVIRONMENTAL SCIENCES↗