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,351 records · Page 75

On the modeling of the scalar correlations necessary to construct a second-order closure description of turbulent reacting flows

A tentative scheme is proposed for constructing second-order closure models of turbulent reacting flows, with particular reference to mixing processes in turbulent diffusion flows. The model is based on the concepts of typical structure and typical eddy, and treats various hypothetical turbulent flame situations, including the time dependence of species mass fractions passing a point in the flame, the most typical occurrence of species mass fractions, and the case of a rather bizarre typical eddy which shows the effects of nonmixedness.

Donaldson, C. DUP.↗

Vibration characteristics of the HPOTP (High-Pressure Oxygen Turbopump) of the SSME (Space Shuttle Main Engine)

Attention is given to rotor dynamic problems that have been encountered and eliminated in the course of Space Shuttle Main Engine (SSME) development, as well as continuing, subsynchronous problems which are being encountered in the development of a 109-percent power level engine. The basic model for the SSME's High Pressure Oxygen Turbopump (HPOTP) encompasses a structural dynamic model for the rotor and housing, and component models for the liquid and gas seals, turbine clearance excitation forces, and impeller diffuser forces. Linear model results are used to examine the synchronous response and stability characteristics of the HPOTP, with attention to bearing load and stability problems associated with the second critical speed. Differences between linear and nonlinear model results are discussed and explained in terms of simple models. Simulation results indicate that while synchronous bearing loads can be reduced, subsynchronous motion is not eliminated by seal modifications.

Childs, D. W.↗

CaloTrilogy: Toward a Breakthrough in One-Step, End-to-End, Physics-Guided Shower Generation for Modern Calorimeters

High-precision calorimeter simulation at current and future colliders imposes rapidly growing computational demands, motivating the development of machine-learning surrogates for traditional Monte Carlo tools such as Geant4. Flow matching and diffusion-based generative models have become leading approaches for high-dimensional fast simulation because of their sample quality, but typically require ${\cal O}(100)$ function evaluations at inference and often rely on auxiliary networks to constrain global observables, compromising streamlined end-to-end generation. We introduce a unified framework that improves the balance between speed, shower quality, and physics fidelity. The method combines: (i) an average velocity field integrator that enables sampling in one or a few evaluations; (ii) a learned generative prior in shower space, constructed from data rather than random noise; and (iii) physics-guided loss terms that impose inductive biases on key observables during training. These elements are training time regularizers, preserving end-to-end inference with no additional cost. With only one or a few evaluation steps, the model achieves shower quality competitive with state-of-the-art flow and diffusion approaches, tested on several public high granularity calorimeter datasets. The results demonstrate inter-layer shower structure consistent with the underlying physics, providing a strong candidate for future fast simulation workflows.

Jiang, Cheng [Edinburgh U.]↗

Radial diffusion in Saturn's radiation belts - A modeling analysis assuming satellite and ring E absorption

A modeling analysis is carried out of six experimental phase space density profiles for nearly equatorially mirroring protons using methods based on the approach of Thomsen et al. (1977). The form of the time-averaged radial diffusion coefficient D(L) that gives an optimal fit to the experimental profiles is determined under the assumption that simple satellite plus Ring E absorption of inwardly diffusing particles and steady-state radial diffusion are the dominant physical processes affecting the proton data in the L range that is modeled. An extension of the single-satellite model employed by Thomsen et al. to a model that includes multisatellite and ring absorption is described, and the procedures adopted for estimating characteristic satellite and ring absorption times are defined. The results obtained in applying three representative solid-body absorption models to evaluate D(L) in the range where L is between 4 and 16 are reported, and a study is made of the sensitivity of the preferred amplitude and L dependence for D(L) to the assumed model parameters. The inferred form of D(L) is then compared with that which would be predicted if various proposed physical mechanisms for driving magnetospheric radial diffusion are operative at Saturn.

Hood, L. L.↗

Li-7 abundances in halo stars: Testing stellar evolution models and the primordial Li-7 abundance

A large number of stellar evolution models with (Fe/H) = -2.3 and -3.3 have been calculated in order to determine the primordial Li-7 abundance and to test current stellar evolution models by a comparison to the extensive database of accurate Li abundances in extremely metal-poor halo stars observed by Thorburn (1994). Standard models with gray atmospheres do a very good job of fitting the observed Li abundances in stars hotter than approximately 5600 K. They predict a primordial. Li-7 abundance of log N(Li) = 2.24 +/- 0.03. Models which include microscopic diffusion predict a downward curvature in the Li-7 destruction isochrones at hot temperatures which is not present in the observations. Thus, the observations clearly rule out models which include uninhibited microscopic diffusion of Li-7 from the surface of the star. Rotational mixing inhibits the microscopic diffusion and the (Fe/H) = -2.28 stellar models which include both diffusion and rotational mixing provide an excellent match to the mean trend in T(sub eff) which is present in the observations. Both the plateau stars and the heavily depleted cool stars are well fit by these models. The rotational mixing leads to considerable Li-7 depletion in these models and the primordial Li-7 abundance inferred from these models is log N(Li) = 3.08 +/- 0.1. However, the (Fe/H) = -3.28 isochrones reveal problems with the combined models. These isochrones predict a trend of decreasing log N(Li) with increasing T(sub eff) which is not present in the observations. Possible causes for this discrepancy are discussed.

Chaboyer, Brian↗

Modeling the impact of radiation-enhanced diffusion on implanted ion profiles

We report that ion irradiations are a valuable research tool for exploring radiation effects in materials. However, it is well recognized that the implanted ions can artificially modify the radiation damage evolution, e.g., enhancing amorphization processes at low irradiation temperatures and suppressing void swelling at elevated temperatures. Therefore, the implanted ion region should be avoided for most studies of ion irradiated materials. Due to increased interest in high damage, high temperature ion irradiations studying radiation effects in materials for proposed high dose Generation IV fission and fusion energy applications, it is crucial to quantify the extent of diffusional broadening of the implanted ion profile for a variety of temperatures, irradiation fluxes, and sink strengths. The present study summarizes computational analyses of thermal and depth-dependent radiation enhanced difusion (RED) on diffusion broadening of the implanted ion profiles in Fe and Ni for a variety of irradiation conditions. For a low assumed RED coefficient (10 -20 -10 -19 m 2 /s and 10 -4 and 10 -3 peak dpa/s, respectively) or high flux, broadening of the as-implanted ion profile is very small and a suitable artifact-free midrange region for analysis exists for ion energies above 5-6 MeV at 100 peak dpa. At high RED coefficients (10 -19 -10 -18 m 2 /s and 10 -4 and 10 -3 peak dpa/s, respectively) broadening is much more significant, and no valid region for investigation exists below 6-8 MeV ion energies at any damage level >50 dpa. While increasing flux decreases irradiation time, it also increases the RED coefficient; these effects mostly offset except under recombination-dominant conditions. For typical high dose irradiation studies of void swelling that exceed ~100 displacement per atom (dpa), ion energies > 6-8 MeV must be employed in order to achieve suitable artifact-free midrange irradiation analysis regions, depending on the material system. Lastly, analysis of amorphization/precipitation require even higher energies (>10 MeV).

07 ISOTOPE AND RADIATION SOURCES↗

Diffusion in random networks

The ensemble averaging technique is applied to model mass transport by diffusion in random networks. The system consists of an ensemble of random networks, where each network is made of pockets connected by tortuous channels. Inside a channel, fluid transport is assumed to be governed by the one-dimensional diffusion equation. Mass balance leads to an integro-differential equation for the pocket mass density. The so-called dual-porosity model is found to be equivalent to the leading order approximation of the integration kernel when the diffusion time scale inside the channels is small compared to the macroscopic time scale. As a test problem, we consider the one-dimensional mass diffusion in a semi-infinite domain. Because of the required time to establish the linear concentration profile inside a channel, for early times the similarity variable is xt $-$1/4 rather than xt $-$1/2 as in the traditional theory. We found this early time similarity can be explained by random walk theory through the network.

97 MATHEMATICS AND COMPUTING↗

Matrix Diffusion in Fractured Media: New Insights Into Power Law Scaling of Breakthrough Curves

Here, we develop a theoretical model for power law tailing behavior of transport in fractured rock based on the relative dominance of the decay rate of the advective travel time distribution, modeled using a Pareto distribution (with tail decaying as ~ time -(1+α) ), versus matrix diffusion, modeled using a Lévy distribution. The theory predicts that when the advective travel time distribution decays sufficiently slowly (α<1), the late-time decay rate of the breakthrough curve is -(1+α/2) rather than the classical -3/2. Yet, if α>1, the -3/2 decay rate is recovered. For weak matrix diffusion or short advective first breakthrough times, we identify an early-time regime where the breakthrough curve follows the Pareto distribution, before transitioning to the late-time decay rate. The theoretical predictions are validated against particle tracking simulations in the three-dimensional discrete fracture network simulator dfnWorks, where matrix diffusion is incorporated using a time domain random walk.

58 GEOSCIENCES↗

Construction and Modeling of a Coculture Microplate for Real-Time Measurement of Microbial Interactions

The dynamic structures of microbial communities emerge from the complex network of interactions between their constituent microorganisms. Quantitative measurements of these interactions are important for understanding and engineering ecosystem structure. Here, we present the development and application of the BioMe plate, a redesigned microplate device in which pairs of wells are separated by porous membranes. BioMe facilitates the measurement of dynamic microbial interactions and integrates easily with standard laboratory equipment. We first applied BioMe to recapitulate recently characterized, natural symbiotic interactions between bacteria isolated from the Drosophila melanogaster gut microbiome. Specifically, the BioMe plate allowed us to observe the benefit provided by two Lactobacillus strains to an Acetobacter strain. We next explored the use of BioMe to gain quantitative insight into the engineered obligate syntrophic interaction between a pair of Escherichia coli amino acid auxotrophs. We integrated experimental observations with a mechanistic computational model to quantify key parameters associated with this syntrophic interaction, including metabolite secretion and diffusion rates. This model also allowed us to explain the slow growth observed for auxotrophs growing in adjacent wells by demonstrating that, under the relevant range of parameters, local exchange between auxotrophs is essential for efficient growth. The BioMe plate provides a scalable and flexible approach for the study of dynamic microbial interactions.

3D printed device↗

Morphological evolution of growing crystals - A Monte Carlo simulation

The combined effects of nutrient diffusion and surface kinetics on the crystal morphology were investigated using a Monte Carlo model to simulate the evolving morphology of a crystal growing from a two-component gaseous nutrient phase. The model combines nutrient diffusion, based on a modified diffusion-limited aggregation process, with anisotropic surface-attachment kinetics and surface diffusion. A variety of conditions, ranging from kinetic-controlled to diffusion-controlled growth, were examined. Successive transitions from compact faceted (dominant surface kinetics) to open dendritic morphologies (dominant volume diffusion) were obtained.

Xiao, Rong-Fu↗

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point-source and extended-source studies rely on the modeling of this diffuse emission for accurate characterization. Here, we describe the development of the Galactic Interstellar Emission Model (GIEM),which is the standard adopted by the LAT Collaboration and is publicly available. This model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse-Compton emission produced in the Galaxy. In the GIEM, we also include large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20deg and we observe an enhanced emission toward their base extending in the north and south Galactic directions and located within approximately 4deg of the Galactic Center.

gamma rays: diffuse background – gamma rays: gen↗

Superresolution Diffuse Optical Imaging by Localization of Fluorescence

The multiple scattering of light presents major challenges in realizing useful in vivo imaging at tissue depths of more than about one millimeter, where many answers to health questions lie. Visible through near-infrared photons can be readily and safely detected through centimeters of tissue; however, limited information is available for image formation. One strategy for obtaining images is to model the photon transport and a simple incoherent model is the diffusion equation approximation to the Boltzmann transport equation. Such an approach provides a prediction of the mean intensity of heavily scattered light and hence provides a forward model for optimization-based computational imaging. While diffuse optical imaging methods have received substantial attention, they remain restricted in terms of resolution because of the loss of high-spatial-frequency information that is associated with the multiple scattering of photons. Consequently, only relatively large inhomogeneities, such as tumors or organs in small animals, can be effectively resolved. We introduce a superresolution imaging approach based on point localization in a diffusion framework that enables over two orders of magnitude improvement in the spatial resolution of diffuse optical imaging. The method is demonstrated experimentally by localizing a fluorescent inhomogeneity in a highly scattering slab and characterizing the localization uncertainty. The approach allows imaging through centimeters of tissue with a resolution of tens of microns, thereby enabling cells or cell clusters to be resolved. More generally, this high-resolution imaging approach could be applied with any physical transport or wave model and hence to a broad class of physical problems. Finally, paired with a suitable optical contrast mechanism, as can be realized with targeted fluorescent molecules or genetically modified animals, superresolution diffuse imaging should open alternative dimensions for in vivo applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experiment and analysis of directional effects on radiant heat transfer.

Critical examination of the validity of the commonly used simplified models for predicting local and overall radiant heat interchange among real surfaces. The spectral local irradiation was measured at wavelengths of 3.08 and 4.51 microns in a configuration consisting of three plane parallel surfaces of finite extent. The test surfaces were gold with mechanical roughness and temperature varying from 0.02 to 7.1 microns and 77 to 760 K, respectively. Comparisons between experimental data and predictions for six different models ranging from the simple diffuse to the most detailed diffuse-plus-specular-directional-property model with specular component calculated according to Beckmann are presented. It is shown that the calculation of the overall irradiation using appropriate constant-property models agrees well with the experimental data to within combined experimental and analytical uncertainty. In general, it is concluded that for accurate prediction of local irradiation the directional characteristics of the surfaces must be considered.

Toor, J. S.↗

SkyFACT: high-dimensional modeling of gamma-ray emission with adaptive templates and penalized likelihoods

We present SkyFACT (Sky Factorization with Adaptive Constrained Templates), a new approach for studying, modeling and decomposing diffuse gamma-ray emission. Like most previous analyses, the approach relies on predictions from cosmic-ray propagation codes like GALPROP and DRAGON. However, in contrast to previous approaches, we account for the fact that models are not perfect and allow for a very large number (∼> 10{sup 5}) of nuisance parameters to parameterize these imperfections. We combine methods of image reconstruction and adaptive spatio-spectral template regression in one coherent hybrid approach. To this end, we use penalized Poisson likelihood regression, with regularization functions that are motivated by the maximum entropy method. We introduce methods to efficiently handle the high dimensionality of the convex optimization problem as well as the associated semi-sparse covariance matrix, using the L-BFGS-B algorithm and Cholesky factorization. We test the method both on synthetic data as well as on gamma-ray emission from the inner Galaxy, |ℓ|<90{sup o} and | b |<20{sup o}, as observed by the Fermi Large Area Telescope. We finally define a simple reference model that removes most of the residual emission from the inner Galaxy, based on conventional diffuse emission components as well as components for the Fermi bubbles, the Fermi Galactic center excess, and extended sources along the Galactic disk. Variants of this reference model can serve as basis for future studies of diffuse emission in and outside the Galactic disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Hydrocarbon photochemistry in the upper atmosphere of Jupiter

The hydrocarbon photochemistry in the upper atmosphere of Jupiter is investigated using a one-dimensional, photochemical-diffusive, and diurnally averaged model. The important chemical cycles and pathways among the major species are outlined and a standard model for the North Equatorial Belt region is examined in detail. It is found that several traditionally dominant chemical pathways among the C and C2 species are replaced in importance by cycles involving C-C4 species. The pressure and altitude profiles of mixing ratios for several observable hydrocarbon species are compared with available ultraviolet- and infrared-derived abundances. The results of sensitivity studies on the standard model with respect to variations in eddy diffusion profile, solar flux, atomic hydrogen influx, latitude, temperature, and important chemical reaction rates are presented. Measured and calculated airglow emissions of He at 584 angstroms and H at 1216 angstroms are also used to provide some constraints on the range of model parameters. The relevance of the model results to the upcoming Galileo mission is briefly discussed. The model is subject to considerable improvement; there is a great need for laboratory measurements of basic reaction rates and photodissociation quantum yields, even for such simple species as methylacetylene and allene. Until such laboratory measurements exist there will be considerable uncertainty in the understanding of the C3 and higher hydrocarbons in the atmospheres of the jovian planets.

NASA Discipline Number 52-20↗

Nanoparticle dynamics in the spatial afterglows of nonthermal plasma synthesis reactors

Nonthermal plasma flow tube reactors are industrially scalable systems for the production of nanocrystal (NC) based materials and coatings. One key advantage of nonthermal plasma synthesis is the ability to both synthesize NCs and deposit films in a single reactor, as at the reactor outlet, NCs can be inertially deposited onto a target substrate. The size and morphology of deposited particles can substantially influence the film structure and function. Though NCs are typically near-spherical and monodispersed as-produced in plasma synthesis reactors, NC charge and growth dynamics can be altered substantially when NCs are sampled out of the plasma and through the spatial afterglow region, affecting deposition. Experiments have demonstrated changes of NC size and charge in the spatial afterglow; however, these dynamics remain unexplored and unexplained via theory and simulation. To address this, we developed a constant number Monte Carlo (CNMC) simulation model to examine the mechanisms of NC decharging and growth in the spatial afterglow of plasma flow tube reactors. Collisions between NC and plasma species, diffusive deposition, and electron desorption from NCs are incorporated in the CNMC simulation. The simulation results are specifically compared with previous experiments on Si NCs synthesized from a low pressure Ar-SiH 4 nonthermal plasma reactor. Furthermore, the experiment-model comparison shows that CNMC models can be implemented which accurately model NC size distribution evolution in a spatial afterglow. Simultaneously, results show that improved collision models, energetic species diffusion models, and electron desorption models will be necessary to accurately depict NC dynamics in spatial afterglows.

42 ENGINEERING↗

Understanding the structure and mechanism of Na + diffusion in NASICON solid-state electrolytes and the effect of Sc- and Al/Y-substitution

NASICON (sodium superionic conductor) based ceramics are one of the most promising classes of solid-state electrolytes for all-solid-state batteries. However, the mechanism of sodium ion diffusion is not understood in great detail since there is still a discrepancy between reported average structure models, local structures, and the number and position of sodium sites. To close this gap, we investigate the underlying diffusion mechanism and structural changes governing the Na + transport in Na 3.4 Zr 2 Si 2.4 P 0.6 O 12 using quasielastic neutron scattering (QENS) and powder X-ray diffraction (XRD). In the temperature range from 298 K to 640 K, the correlations between structural changes of a monoclinic C2/c to rhombohedral R $\bar{3}$c phase transition and the result of ion diffusion are investigated. The analysis of the quasielastic neutron scattering data reveals two quasielastic components corresponding to the Chudley-Elliott jump-diffusion model. It clearly shows two different Na + diffusion processes, local and long-range, on two different time and length scales and allows calculations of their corresponding activation energies. Additionally, the effects of Sc 3+ and Al 3+ /Y 3+ aliovalent substitution of Zr 4+ ions on the crystal structure and Na + diffusion are also studied. We can distinguish a local, chain, and cross-chain diffusion mechanism based on correlated QENS and XRD comparison of relevant nearest crystallographic Na–Na distances. The results reveal that the Na + diffusion in these NASICONs is three-dimensional and can provide guidelines on how dopants and changes in the crystal structure can affect the Na + conductivity.

25 ENERGY STORAGE↗