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 271 records · Page 15

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.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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.↗

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↗

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.↗

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↗

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗

Ensemble‐Based, Large‐Eddy Reconstruction of Wind Turbine Inflow in a Near‐Stationary Atmospheric Boundary Layer Through Generative Artificial Intelligence

ABSTRACT To validate the second‐by‐second dynamics of turbines in field experiments, it is necessary to accurately reconstruct the winds going into the turbine. Current time‐resolved inflow reconstruction techniques estimate wind behavior in unobserved regions using relatively simple spectral‐based models of the atmosphere. Here, we develop a technique for time‐resolved inflow reconstruction that is rooted in a large‐eddy simulation model of the atmosphere. Our “large‐eddy reconstruction” technique blends observations and atmospheric model information through a diffusion model machine learning algorithm, allowing us to generate probabilistic ensembles of reconstructions for a single 10‐min observational period. Our generated inflows can be used directly by aeroelastic codes or as inflow boundary conditions in a large‐eddy simulation. We verify the second‐by‐second reconstruction capability of our technique in three synthetic field campaigns, finding positive Pearson correlation coefficient values () between ground‐truth and reconstructed streamwise velocity, as well as smaller positive correlation coefficient values for unobserved fields (spanwise velocity, vertical velocity, and temperature). We validate our technique in three real‐world case studies by driving large‐eddy simulations with reconstructed inflows and comparing to independent inflow measurements. The reconstructions are visually similar to measurements, follow desired power spectra properties, and track second‐by‐second behavior ().

17 WIND ENERGY↗

Influence of carrier trapping on radiation detection properties in $\mathrm{CVD}$ grown 4H-SiC epitaxial layers with varying thickness up to 250 µm

We report low doped n-type 4H-SiC epitaxial layers of thickness 50, 150, and 250 µm grown by hot wall chemical vapor deposition were used to fabricate Ni/4H-SiC Schottky radiation detectors. The epitaxial layers were grown on the (0001) face of highly conductive bulk 4H-SiC substrates 8° offcut towards the $\langle 11\bar{2}1\rangle$ direction. The 50, 150, and 250 µm thick epilayer detectors, under optimized settings, showed energy resolutions of 2.0%, 0.78%, and 0.63%, respectively for 5.48 MeV alpha particles. Deep level transient spectroscopy studies showed that the observed variation in the detector resolution is linked to the defect parameters in the devices. Least-squares fitting of the bias dependence of the charge collection efficiency according to a drift-diffusion model, revealed minority carrier diffusion lengths of 16, 10, and 9.2 µm, respectively, implying that the detectors are not limited by minority carrier trapping. The detector performance was observed to be primarily dependent on the concentration and capture cross-sections of the lifetime killing electron traps Z 1/2 and EH 6/7 .

36 MATERIALS SCIENCE↗

Verification methods for drift–diffusion reaction models for plasma simulations

Abstract Compared to other computational physics areas such as codes for general computational fluid dynamics, the documentation of verification methods for plasma fluid codes remains under developed. Current analytical solutions for plasma are often highly limited in terms of testing highly coupled physics, due to the harsh assumptions needed to derive even simple plasma equations. This work highlights these limitations, suggesting the method of manufactured solutions (MMSs) as a potential option for future verification efforts. To demonstrate the flexibility of MMS in verifying these highly coupled systems, the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework was utilized. Thanks to the MOOSE framework’s robustness and modularity, as well as to its physics module capabilities and ecosystem applications (i.e. Zapdos and the chemical reaction network) developed for plasma physics modeling and simulation, this report lays the groundwork for a structured method of conducting plasma fluid code verification.

DeChant, Corey (ORCID:0000000293779732)↗