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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,747 records · Page 97

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Selective Uptake of Ethane/Ethylene Mixtures by UTSA-280 is Driven by Reversibly Coordinated Water Defects

The metal–organic framework (MOF) UTSA-280 has been shown previously to separate ethane/ethylene mixtures. We show that this separation is kinetic; if full equilibrium is reached, the MOF is not selective. The time scales associated with this kinetic separation are sufficient to make the material potentially interesting for practical separations. The observation of kinetic separation is surprising because the MOF’s one-dimensional (1D) channels are small enough that molecular motion along the channels would be expected to take place by single-file diffusion, which is not observed. We show using a combination of experimental and modeling techniques that diffusion of molecules between adjacent 1D channels through local defects associated with coordinated water molecules allows net mixing of molecules and, moreover, that this process is potentially responsible for the selective nature of the observed kinetic separation. We show how activation and regeneration conditions play an essential role in controlling the density of these local defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

CTF Void Drift Validation Study

This milestone report is a summary of work performed in support of expansion of the validation and verification (V&V) matrix for the thermal-hydraulic subchannel code, CTF. The focus of this study is on validating the void drift modeling capabilities of CTF and verifying the supporting models that impact the void drift phenomenon. CTF uses a simple turbulent-diffusion approximation to model lateral cross-flow due to turbulent mixing and void drift. The void drift component of the model is based on the Lahey and Moody model. The models are a function of two-phase mass, momentum, and energy distribution in the system; therefore, it is necessary to correctly model the ow distribution in rod bundle geometry as a first step to correctly calculating the void distribution due to void drift.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

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

NASA NTRS1984

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

DOE OSTI2026

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

NASA NTRS1983

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

DOE OSTI2016

GRIZZLY MODEL OF FULLY COUPLED HEAT TRANSFER, MOISTURE DIFFUSION, ALKALI-SILICA REACTION AND FRACTURING PROCESSES IN CONCRETE

Idaho National Laboratory (INL) has been developing a fully coupled thermo-hydro-mechanical-chemical (THMC) concrete model, referred as Grizzly, with the ultimate goal of reliably simulating and predicting the long-term performance of concrete structures subjected to various aging and degradation mechanisms. In this paper, we present a fully coupled concrete model implemented in our Grizzly code for simulating coupled moisture/heat transfer, stress-dependent ASR swelling and elastic/inelastic mechanics. Example simulation results and preliminary model validation results against experimental data reported in the literature will be provided also in this paper. The close matches between the simulation results and experiments clearly reveal the potential of the fully coupled solution strategy for solving tightly coupled THMC concrete models. We also present simulation results of using a rate-dependent damaging mechanics model for simulating the ASR swelling induced degradation of concrete structures and its feedback to moisture and heat transport processes.

99 GENERAL AND MISCELLANEOUS↗

NASA NTRS1994

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↗

DOE OSTI2014

Transport of an interacting Bose gas in 1D disordered lattices

We use ultracold atoms in a quasiperiodic lattice to study two outstanding problems in the physics of disordered systems: a) the anomalous diffusion of a wavepacket in the presence of disorder, interactions and noise; b) the transport of a disordered superfluid. a) Our results show that the subdiffusion, observed when interaction alone is present, can be modelled with a nonlinear diffusion equation and the peculiar shape of the expanding density profiles can be connected to the microscopic nonlinear diffusion coefficients. Also when noise alone is present we can describe the observed normal diffusion dynamics by existing microscopic models. In the unexplored regime in which noise and interaction are combined, instead, we observe an anomalous diffusion, that we model with a generalized diffusion equation, where noise- and interaction-induced contributions add each other. b) We find that an instability appearing at relatively large momenta can be employed to locate the fluid-insulator crossover driven by disorder. By investigating the momentum-dependent transport, we observe a sharp crossover from a weakly dissipative regime to a strongly unstable one at a disorder-dependent critical momentum. The set of critical disorder and interaction strengths for which such critical momentum vanishes, can be identified with the separation between a fluid regime and an insulating one and can be related to the predicted zero-temperature superfluid-Bose glass transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

SU-E-T-333: Towards Customizable Radiotherapy Enhancement (CuRE) for Prostate Cancer Using Cisplatin Nanoparticles

Purpose: Replacing routinely used brachytherapy spacers with multifunctional ones loaded with cisplatin nanoparticles (CNP), which can be released into the tumor after implantation, could enable customizable radiation boosting to the prostate tumor in addition to chemotherapy effect. This study investigates the feasibility of customizing the intra-tumor biodistribution and corresponding dose enhancement (DEF) over time for the released CNP as a function of nanoparticle size. Methods: Dose enhancement factors (DEF) due to photon-induced emission of photo-/Auger electrons from CNPs were calculated as a function of concentration using previously published analytical calculation method. An experimentally determined diffusion coefficient (D) for 10 nm nanoparticles in mouse tumor model was employed to estimate D for other sizes using the Stoke- Einstein equation. The error function diffusion model in the experimental study was applied to generate the intra-tumor concentration profile for a burst release of CNPs from the spacer over time. The corresponding DEF profiles were then determined for brachytherapy using Pd-103 and I-125 sources. Results: As expected, the generated profiles showed greater DEF over time for smaller CNP sizes at sample distances from the spacer. For example, for a centrally located spacer, clinically significant DEF (> 20%) could be achieved near the tumor periphery (ca. 0.85 cm distance from the spacer for average PCa tumor size) after 20, and 100 days, respectively for CNPs sizes of 2 nm, and 10 nm, using I-125. Meanwhile for Pd-103, clinically significant DEF could be achieved at the same position after 22 and 108 days, respectively, for same size particles. Conclusion: Our preliminary results demonstrate the feasibility of customizing dose enhancement to prostate tumors as a function of spacer location, brachytherapy source type or size of CNPs released from multifunctional spacers. Such an approach could enable customizable radiation boosting to tumor sub-volumes, while minimizing dose to healthy tissues.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

Interstellar simulations using a unified microscopic-macroscopic Monte Carlo model with a full gas-grain network including bulk diffusion in ice mantles

We have designed an improved algorithm that enables us to simulate the chemistry of cold dense interstellar clouds with a full gas-grain reaction network. The chemistry is treated by a unified microscopic-macroscopic Monte Carlo approach that includes photon penetration and bulk diffusion. To determine the significance of these two processes, we simulate the chemistry with three different models. In Model 1, we use an exponential treatment to follow how photons penetrate and photodissociate ice species throughout the grain mantle. Moreover, the products of photodissociation are allowed to diffuse via bulk diffusion and react within the ice mantle. Model 2 is similar to Model 1 but with a slower bulk diffusion rate. A reference Model 0, which only allows photodissociation reactions to occur on the top two layers, is also simulated. Photodesorption is assumed to occur from the top two layers in all three models. We found that the abundances of major stable species in grain mantles do not differ much among these three models, and the results of our simulation for the abundances of these species agree well with observations. Likewise, the abundances of gas-phase species in the three models do not vary. However, the abundances of radicals in grain mantles can differ by up to two orders of magnitude depending upon the degree of photon penetration and the bulk diffusion of photodissociation products. We also found that complex molecules can be formed at temperatures as low as 10 K in all three models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

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↗

DOE OSTI2017

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↗

DOE OSTI2019

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↗

DOE OSTI2016

A GENERALIZED TWO-COMPONENT MODEL OF SOLAR WIND TURBULENCE AND AB INITIO DIFFUSION MEAN-FREE PATHS AND DRIFT LENGTHSCALES OF COSMIC RAYS

We extend a two-component model for the evolution of fluctuations in the solar wind plasma so that it is fully three-dimensional (3D) and also coupled self-consistently to the large-scale magnetohydrodynamic equations describing the background solar wind. The two classes of fluctuations considered are a high-frequency parallel-propagating wave-like piece and a low-frequency quasi-two-dimensional component. For both components, the nonlinear dynamics is dominanted by quasi-perpendicular spectral cascades of energy. Driving of the fluctuations by, for example, velocity shear and pickup ions is included. Numerical solutions to the new model are obtained using the Cronos framework, and validated against previous simpler models. Comparing results from the new model with spacecraft measurements, we find improved agreement relative to earlier models that employ prescribed background solar wind fields. Finally, the new results for the wave-like and quasi-two-dimensional fluctuations are used to calculate ab initio diffusion mean-free paths and drift lengthscales for the transport of cosmic rays in the turbulent solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

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↗

DOE OSTI2016

HARD X-RAY EMISSIONS FROM CASSIOPEIA A OBSERVED BY INTEGRAL

Cassiopeia A (Cas A), as the nearby young remnant of a core-collapse supernova, is the best candidate for astrophysical studies in supernova explosion and its environment. We studied the hard X-ray emission from Cas A using the 10 year data of INTEGRAL observations, and first detected non-thermal continuum emission from the source up to 220 keV. The {sup 44}Ti line emission at 68 and 78 keV is confirmed by our observations with a mean flux of ∼(2.2 ± 0.4) × 10{sup 5} ph cm{sup 2} s{sup 1}, corresponding to a {sup 44}Ti yield in Cas A of (1.3 ± 0.4) × 10{sup 4} M {sub ⊙}. The continuum emission from 3 to 500 keV can be fit with a thermal bremsstrahlung of kT ∼ 0.79 ± 0.08 keV plus a power-law model of Γ ∼ 3.13 ± 0.03. The non-thermal emission from Cas A is well fit by a power-law model without a cutoff up to 220 keV. This radiation characteristic is inconsistent with diffusive shock acceleration models with a remnant shock velocity of only 5000 km s{sup 1}. The central compact object in Cas A cannot significantly contribute to the emission above 80 keV. Some possible physical origins of the non-thermal emission above 80 keV from the remnant shock are discussed. We deduce that the asymmetrical supernova explosion scenario of Cas A is a promising scenario for the production of high-energy synchrotron radiation photons, where a portion of the ejecta with a velocity of ∼0.1 c and opening angle of ∼10° can account for the 100 keV emission, as is consistent with the “jet” observed in Cas A.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1988

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↗