Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effectiveness of denoising diffusion probabilistic models for fast and high-fidelity whole-event simulation in high-energy heavy-ion experiments

Artificial intelligence (AI) generative models, such as generative adversarial networks (GANs), variational autoencoders, and normalizing flows, have been widely used and studied as efficient alternatives for traditional scientific simulations. However, they have several drawbacks, including training instability and inability to cover the entire data distribution, especially for regions where data are rare. This is particularly challenging for whole-event, full-detector simulations in high-energy heavy-ion experiments, such as sPHENIX at the Relativistic Heavy Ion Collider and Large Hadron Collider experiments, where thousands of particles are produced per event and interact with the detector. This work investigates the effectiveness of denoising diffusion probabilistic models (DDPMs) as an AI-based generative surrogate model for the sPHENIX experiment that includes the heavy-ion event generation and response of the entire calorimeter stack. DDPM performance in sPHENIX simulation data is compared with a popular rival, GANs. Results show that both DDPMs and GANs can reproduce the data distribution where the examples are abundant (low-to-medium calorimeter energies). Nonetheless, DDPMs significantly outperform GANs, especially in high-energy regions where data are rare. Additionally, DDPMs exhibit superior stability compared to GANs. The results are consistent between both central and peripheral centrality heavy-ion collision events. Moreover, DDPMs offer a substantial speedup of approximately a factor of 100 compared to the traditional Geant4 simulation method.

42 ENGINEERING↗
DOE OSTI2020

Nonadiabatic Effects on Defect Diffusion in Silicon-Doped Nanographenes

Single atom impurities in graphene, substitutional silicon defects in particular, have been observed to diffuse under electron beam irradiation. However, the relative importance of elastic and inelastic scattering in facilitating their mobility remains unclear. Here, we employ excited-state electronic structure calculations to explore potential inelastic effects, and find an electronically nonadiabatic excited-state silicon diffusion pathway involving “softened” Si–C bonding that presents an ~2 eV lower diffusion barrier than the ground-state pathway. Beam-induced transition rates to this state indicate that the excited-state pathway is accessible through irradiation of the defect site. However, even in the limit of fully elastic scattering, upward nonadiabatic transitions are also possible along the diffusion coordinate, increasing the diffusion barrier and further demonstrating the potential for electronic nonadiabaticity to influence beam-induced atomic transformations in materials. Additionally, we also propose some experimentally testable signatures of such excited-state pathways.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Effects of Moisture on Diffusion in Unmodified Wood Cell Walls: A Phenomenological Polymer Science Approach

Despite the importance of cell wall diffusion to nearly all aspects of wood utilization, diffusion mechanisms and the detailed effects of moisture remain poorly understood. In this perspective, we introduce and employ approaches established in polymer science to develop a phenomenological framework for understanding the effects of moisture on diffusion in unmodified wood cell walls. The premise for applying this polymer-science-based approach to wood is that wood polymers (cellulose, hemicelluloses, and lignin) behave like typical solid polymers. Therefore, the movement of chemicals through wood cell walls is a diffusion process through a solid polymer, which is in contrast to previous assertions that transport of some chemicals occurs via aqueous pathways in the cell wall layers. Diffusion in polymers depends on the interrelations between free volume in the polymer matrix, molecular motions of the polymer, diffusant dimensions, and solubility of the diffusant in the polymer matrix. Because diffusion strongly depends on whether a polymer is in a rigid glassy state or soft rubbery state, it is important to understand glass transitions in the amorphous wood polymers. Through a review and analysis of available literature, we conclude that in wood both lignin and the amorphous polysaccharides very likely have glass transitions. After developing and presenting this polymer-science-based perspective of diffusion through unmodified wood cell walls, suggested directions for future research are discussed. A key consideration is that a large difference between diffusion through wood polymers and typical polymers is the high swelling pressures that can develop in unmodified wood cell walls. This pressure likely arises from the hierarchical structure of wood and should be taken into consideration in the development of predictive models for diffusion in unmodified wood cell walls.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

The kinetic Sunyaev-Zel'dovich effect from the diffuse gas in the Local Group

Since the Local Group (LG) of galaxies moves with a bulk velocity with respect to the cosmic microwave background radiation (CMB), free electrons in its gaseous halo should imprint large-scale non-primordial temperature shifts in the CMB via the kinetic Sunyaev-Zel'dovich (kSZ) effect. By modeling the distribution of gas in the LG halo and using its inferred velocity with respect to the CMB, we calculate the resulting kSZ signal from the diffuse LG medium. We find that it is dominated by a hot spot ∼ 10° in size in the direction of M31, where the optical depth of free electrons is the greatest. By performing a correlation analysis, we find no statistical evidence that the kSZ signal from model of the LG halo is embedded in the CMB temperature map measured by the Planck satellite. We constrain the amount of mass in the LG halo by limiting the kSZ temperature shift around the hot spot to be smaller than the observed temperature shift in the Planck map. We find the tightest constraints for models where the halo mass is highly concentrated, with the mass limited to roughly 2.5–5 × 10{sup 12}M{sub ⊙}, but note that halos with such high concentrations are rare.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Comparison of experimental methods for estimating matrix diffusion coefficients for contaminant transport modeling

Here, diffusion cell and diffusion wafer experiments were conducted to compare methods for estimating effective matrix diffusion coefficients in rock core samples from Pahute Mesa at the Nevada Nuclear Security Site (NNSS). A diffusion wafer method, in which a solute diffuses out of a rock matrix that is pre-saturated with water containing the solute, is presented as a simpler alternative to the traditional through-diffusion (diffusion cell) method. Both methods yielded estimates of effective matrix diffusion coefficients that were within the range of values previously reported for NNSS volcanic rocks. The difference between the estimates of the two methods ranged from 14 to 30%, and there was no systematic high or low bias of one method relative to the other. From a transport modeling perspective, these differences are relatively minor when one considers that other variables (e.g., fracture apertures, fracture spacings) influence matrix diffusion to a greater degree and tend to have greater uncertainty than effective matrix diffusion coefficients. For the same relative random errors in concentration measurements, the diffusion cell method yields effective matrix diffusion coefficient estimates that have less uncertainty than the wafer method. However, the wafer method is easier and less costly to implement and yields estimates more quickly, thus allowing a greater number of samples to be analyzed for the same cost and time. Given the relatively good agreement between the methods, and the lack of any apparent bias between the methods, the diffusion wafer method appears to offer advantages over the diffusion cell method if better statistical representation of a given set of rock samples is desired.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Steady state effects in a two-pulse diffusion-weighted sequence

In conventional nuclear magnetic resonance (NMR) diffusion measurements a significant amount of experimental time is used up by magnetization recovery, serving to prevent the formation of the steady state, as in the latter case the manifestation of diffusion is modulated by multiple applications of the pulse sequence and conventional diffusion coefficient inference procedures are generally not applicable. Here, an analytical expression for diffusion-related effects in a two-pulse NMR experiment (e.g., pulsed-gradient spin echo) in the steady state mode (with repetition times less than the longitudinal relaxation time of the sample) is derived by employing a Fourier series expansion within the solution of the Bloch-Torrey equations. Considerations are given for the transition conditions between the full relaxation and the steady state experiment description. The diffusion coefficient of a polymer solution (polyethylene glycol) is measured by a two-pulse sequence in the full relaxation mode and for a range of repetition times, approaching the rapid steady state experiment. The precision of the fitting employing the presented steady state solution by far exceeds that of the conventional fitting. Additionally, numerical simulations are performed yielding results strongly supporting the proposed description of the NMR diffusion measurements in the steady state.

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

On the Degradation of the Effectiveness of Nonlinear Diffusion Acceleration with Parallel Block Jacobi Splitting

MOOSE's radiation transport module Rattlesnake comprises a first order S N solver that can be accelerated with a flexible Nonlinear Diffusion Acceleration. At each NDA iteration, the streaming and collision operator of the S{sub N} equations are (approximately) inverted using a transport sweep. When executed on multiple processors parallel block Jacobi (PBJ) splitting is utilized lagging the processor sub-domain's incoming angular fluxes making the implemented method asynchronous, i.e. the convergence properties of NDA become dependent on the number of concurrent tasks and details of the partitioning. Within this work, we study the NDA convergence properties for the two-dimensional C5G7 benchmark problem PBJ with p = 1, 2, 4, 8, 12, 16, 24, 32, 48, 60, 72, 96 and 144 concurrent tasks. Rosa studied the convergence properties of PBJ splitting methods with and without transport synthetic acceleration (TSA) in slab geometry and two-dimensional Cartesian geometry, respectively. Rosa's iteration scheme lags both the incoming angular fluxes and the scattering source making it very similar to the method described here. Rosa found by means of Fourier analysis that convergence of unaccelerated PBJ can degrade for optically thin subdomains, i.e. for a large number of processors given a fixed domain size. Depending on the choice of parameters application of TSA can lead to an effective algorithm. In contrast to Rosa's work, this work focuses on PBJ splitting for NDA accelerated solution of the multigroup eigenvalue problem. Later, Zerr and Anistratov used their integral transport matrix method in conjunction with PBJ. This iteration scheme is different as it does not lag the scattering source. Nevertheless, Fourier analysis with Diamond Difference spatial discretization showed that the method degrades with increasing the overall size of the spatial domain while keeping the sub-domain size constant. In the case of an infinite medium, the observed spectral radius becomes unity. These references establish the iterative performance of PBJ for one-group fixed source problems discretized on rectangular cells, but for the Rattlesnake first order NDA method, an emphasis lies on multigroup eigenvalue problems. In this setting, the additional complication arises that the low order equations are not exactly inverted, but only to some iterative criterion. Further, as the low order problem is solved using Jacobian-Free Newton Krylov methods, a suitable preconditioner for the GMRES iterations is required. As the effectiveness of the preconditioners may decrease when more concurrent processes are used, we expect the low order solve to impact the overall effectiveness of the NDA method in parallel execution as well. Therefore, this work focuses on the question of how detrimental the asynchronicity of the high and low order iterative solutions is for the effectiveness of NDA solutions of realistic, unstructured grid, multigroup eigenvalue problems. The iterative convergence properties of Rattlesnake's NDA method for multigroup eigenvalue problems are studied for varying number of concurrent processes. Rattlesnake's NDA method is asynchronous because the high order transport sweep uses PBJ splitting and because the low order solve's preconditioners, algebraic multigrid and additive Schwarz, are asynchronous. We found that up to 144 processes the number of Picard iterations increases by a factor of 2-3 affecting but not destroying NDA's parallel performance for the C5G7 benchmark problem, a representative LWR (light water reactor) configuration. The solution of the NDA low order eigenvalue problem posed greater problems than expected for the solvers currently available in Rattlesnake mainly because of problems to effectively precondition the modified low-order problem. As a matter of fact, the low order problem becomes an advection-diffusion problem. Non-robustness was observed for the AMG preconditioned coarse mesh acceleration for more than 60 concurrent processes, which was fixed by increasing the number of permissible V-cycles per AMG call. For fine mesh acceleration AMG proved to be ineffective and we resorted to ASM. For both AMG and ASM, increasing number of linear iterations per nonlinear are observed when increasing the number of processes, but for ASM this problem appears to be worse. The results from this work also indicate that initial development efforts should be diverted from implementation of a high-order parallel sweeping algorithm to the effective and robust iterative solution of the low order NDA multigroup eigenvalue problem. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Strengthening model development and effects of low diffusing solutes to coarsening resistance in aluminum alloys

A modified Orowan strengthening model is proposed to account for finite rod-shaped precipitates with hemispherical caps in aluminum alloy systems. A combined computational and experimental approach is used to study the influences of anisotropic Orowan looping and solute-dislocation interaction on temperature-dependent yield strength. Here, the strengthening model is validated with a dataset containing 297 experimental precipitate geometries, chemistries, temperatures, and strength measurements, and achieves a strong predictive correlation of 0.8713 with experimentally measured yield strengths. Under conditions that are applicable to coarsening, constant particle volume fraction and aspect ratio, the model predicts that short rod precipitates provide far superior strengthening effects compared to plate precipitates. A cast Al-Si-Mg-Cu alloy with rod-shaped Q-phase (Al 3 Cu 2 Mg 9 Si 7 ) precipitates was developed with a novel chemistry exploring the use of low-diffusivity elements (Mn, Ni, V, Zr) to limit precipitate coarsening. The thermodynamic behavior of Mn, Ni, V, and Zr across the Q-phase interface is examined using transmission electron microscopy (TEM), first-principles density-functional theory (DFT) calculations. and atom-probe tomography (APT). DFT calculations utilizing TEM identified Q-phase/Al-matrix interfaces show that Mn, Ni, V, and Zr preferentially segregate to the Q-phase precipitate boundaries which suggests inhibition of precipitate coarsening and higher strengths after temperature exposure. Atom-probe tomography confirms solute atom partitioning/segregation at the Q-phase/Al-matrix interface, found in the modified commercial AS7GU alloy (A356 +0.5%Cu), which supports these observations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Osmotic Pressure Gradient Effects on Water Diffusion in Porous Rock: Can They Pervert Permeability Tests?

Abstract Generation of a large network of hydraulic cracks is of key importance not only for the success of fracking of shale but also for the recent scheme of sequestration of CO2 in deep formations of basalt and peridotite, which are mafic and ultramafic rocks that combine chemically with CO2. In numerical simulation of the creation of a fracture network in porous rock, an important goal is to enhance the rock permeability. The objective of this article is to calculate the effect of osmotic pressure gradients caused by gradients of concentration of the ions of Ca, Mg, Na, etc. on the effective permeability of the rock. The basic differential equations are formulated, and their explicit solutions for appropriate initial and boundary conditions are obtained under certain plausible simplifications. The main result is explicit approximate formulas for the critical time before which no water permeation through a test specimen can be observed. Depending on various parameters, this time can be unacceptably long, which is manifested as a zero water outflow. The solution may also explain the unreasonably small permeability values reported for some shales.

Mechanics↗
DOE OSTI2014

Protection of rubbers against aging with the use of structural, diffusion and kinetic effects

Processes of rubber aging in the different kinetic modes and mathematical models of the processes including aging and destruction processes under extreme conditions have been studied. The aging process for rubbers as thermodynamically open nonlinear systems has also been considered in the paper. It has been revealed that the aging process can be controlled by the internal physical-chemical processes organization and by creation of external influences (by organization of thermodynamic forces and flows). According to the Onsager principle, in certain conditions of aging the conjugation of thermodynamic forces and flows is possible. The diffusion and structural aspects of aging in elastomeric rubbers and articles thereof have been considered. The composite antiaging systems of prolonged action based on the use of microparticles with a saturated solution of the antiagers migrating into an elastomeric matrix at the exploitation have been proposed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Effective Stochastic Model For Reactive Transport

We propose an effective stochastic advection-diffusion-reaction (SADR) model. Unlike traditional advection-dispersion-reaction models, the SADR model describes mechanical and diffusive mixing as two separate processes. In the SADR model, the mechanical mixing is driven by random advective velocity with the variance given by the coefficient of mechanical dispersion. The diffusive mixing is modeled as a fickian diffusion with the effective diffusion coefficient. Both coefficients are given in terms of Peclet number (Pe) and the coefficient of molecular diffusion. We use the experimental results of \cite{Gramling} to demonstrate that for transport and bimolecular reactions in porous media the SADR model is significantly more accurate than the traditional dispersion model, which overestimates the mass of the reaction product by as much as 25\%.

stochastic equations, mixing controlled reactions,↗
DOE OSTI2016

Investigating the age distribution of fracture discharge using multiple environmental tracers, Bedrichov Tunnel, Czech Republic

We investigated the transit time distribution (TTD) of discharge collected from fractures in the Bedrichov Tunnel, Czech Republic, using lumped parameter models and multiple environmental tracers. We then utilize time series of δ 18 O, δ 2 H and 3 H along with CFC measurements from individual fractures in the Bedrichov Tunnel of the Czech Republic to investigate the TTD, and the uncertainty in estimated mean travel time in several fracture networks of varying length and discharge. We also compare several TTDs, including the dispersion distribution, the exponential distribution, and a developed TTD which includes the effects of matrix diffusion. The effect of seasonal recharge is explored by comparing several seasonal weighting functions to derive the historical recharge concentration. We identify best fit mean ages for each TTD by minimizing the error-weighted, multi-tracer χ2 residual for each seasonal weighting function. We use this methodology to test the ability of each TTD and seasonal input function to fit the observed tracer concentrations, and the effect of choosing different TTD and seasonal recharge functions on the mean age estimation. We find that the estimated mean transit time is a function of both the assumed TTD and seasonal weighting function. Best fits as measured by the χ2 value were achieved for the dispersion model using the seasonal input function developed here for two of the three modeled sites, while at the third site, equally good fits were achieved with the exponential model and the dispersion model and our seasonal input function. The average mean transit time for all TTDs and seasonal input functions converged to similar values at each location. The sensitivity of the estimated mean transit time to the seasonal weighting function was equal to that of the TTD. These results indicated that understanding seasonality of recharge is at least as important as the uncertainty in the flow path distribution in fracture networks and that unique identification of the TTD and mean transit time is difficult given the uncertainty in the recharge function. But, the mean transit time appears to be relatively robust to the structural model uncertainty. The results presented here should be applicable to other studies using environmental tracers to constrain flow and transport properties in fractured rock systems.

54 ENVIRONMENTAL SCIENCES↗