Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion models”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Modeling of ion acceleration through drift and diffusion at interplanetary shocks

A test particle simulation designed to model ion acceleration through drift and diffusion at interplanetary shocks is described. The technique consists of integrating along exact particle orbits in a system where the angle between the shock normal and mean upstream magnetic field, the level of magnetic fluctuations, and the energy of injected particles can assume a range of values. The technique makes it possible to study time-dependent shock acceleration under conditions not amenable to analytical techniques. To illustrate the capability of the numerical model, proton acceleration was considered under conditions appropriate for interplanetary shocks at 1 AU, including large-amplitude transverse magnetic fluctuations derived from power spectra of both ambient and shock-associated MHD waves.

Decker, R. B.↗
NASA NTRS1991

Active unsteady aerodynamic suppression of rotating stall in an incompressible flow centrifugal compressor with vaned diffuser

A mathematical model is developed to analyze the suppression of rotating stall in an incompressible flow centrifugal compressor with a vaned diffuser, thereby addressing the important need for centrifugal compressor rotating stall and surge control. In this model, the precursor to to instability is a weak rotating potential velocity perturbation in the inlet flow field that eventually develops into a finite disturbance. To suppress the growth of this potential disturbance, a rotating control vortical velocity disturbance is introduced into the impeller inlet flow. The effectiveness of this control is analyzed by matching the perturbation pressure in the compressor inlet and exit flow fields with a model for the unsteady behavior of the compressor. To demonstrate instability control, this model is then used to predict the control effectiveness for centrifugal compressor geometries based on a low speed research centrifugal compressor. These results indicate that reductions of 10 to 15 percent in the mean inlet flow coefficient at instability are possible with control waveforms of half the magnitude of the total disturbance at the inlet.

Lawless, Patrick B.↗
DOE OSTI2020

Reduced Artifact Approach for Determining Diffusion Coefficients in Time-Resolved Microscopy

Ultrafast microscopy methods traditionally assume a Gaussian profile to extract excited state diffusivities from transport measurements. Although this fitting method recovers accurate diffusion coefficients when the point spread function is well-represented by a Gaussian, even minor spatial aberrations introduced by the imaging system cause significant errors in the determined value. To provide a more accurate measure of excited state transport in nano- and microscale materials systems, in this work an alternative analysis protocol is proposed that numerically convolves the Green’s function solution to the diffusion equation with the experimentally measured point spread function. In contrast to the Gaussian fitting approach, the numerical convolution is shown to be robust against artifacts caused by nonideal point spread functions. Furthermore, the numerical convolution approach is highly effective at resolving anisotropic diffusion in modeled data.

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

An efficient approach for treating composition-dependent diffusion within organic particles

Mounting evidence demonstrates that under certain conditions the rate of component partitioning between the gas and particle phase in atmospheric organic aerosol is limited by particle-phase diffusion. To date, however, particle-phase diffusion has not been incorporated into regional atmospheric models. An analytical rather than numerical solution to diffusion through organic particulate matter is desirable because of its comparatively small computational expense in regional models. Current analytical models assume diffusion to be independent of composition and therefore use a constant diffusion coefficient. To realistically model diffusion, however, it should be composition-dependent (e.g. due to the partitioning of components that plasticise, vitrify or solidify). This study assesses the modelling capability of an analytical solution to diffusion corrected to account for composition dependence against a numerical solution. Results show reasonable agreement when the gas-phase saturation ratio of a partitioning component is constant and particle-phase diffusion limits partitioning rate (<10% discrepancy in estimated radius change). However, when the saturation ratio of the partitioning component varies, a generally applicable correction cannot be found, indicating that existing methodologies are incapable of deriving a general solution. Until such time as a general solution is found, caution should be given to sensitivity studies that assume constant diffusivity. Furthermore, the correction was implemented in the polydisperse, multi-process Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) and is used to illustrate how the evolution of number size distribution may be accelerated by condensation of a plasticising component onto viscous organic particles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Analysis based on global model of nitrogen plasma produced by pulsed microwave at low pressure

This paper analyzes certain evolution processes in nitrogen plasmas discharged using pulsed microwaves at low pressure. Comparing the results obtained from the global model incorporating diffusion and the microwave transmission method, the temporal variation of the electron density is analyzed. With a discharge pressure of 300 Pa, the results obtained from experiments and the global model calculation show that when the discharge begins the electron density in the plasma rises quickly, to a level above the critical density corresponding to the discharge microwave frequency, but falls slowly when the discharge microwave pulse is turned off. The results from the global model also show that the electron temperature increases rapidly to a peak, then decays after the electron density reaches the critical density, and finally decreases quickly to room temperature when the discharge microwave pulse is turned off. In the global model, the electron density increases because the high electron temperature induces a high ionization rate. The decay of the electron density mainly comes from diffusion effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Kinetic Monte Carlo Simulation of Hydrogen Diffusion in Tungsten

Nuclear fusion is considered as the most promising sustainable energy source. If achieved, it could provide electricity to the mankind for millions of years. During the reaction, the fusion ash is continuously removed from the magnetic gap, the divertor region, and the divertor facing component is subjected to intense plasma bombardment. It is proposed to use tungsten as the divertor facing material, because it has low physical sputtering yield and superior thermomechanical properties. However, the high mobility of hydrogen in tungsten and the high bombardment flux bring the hydrogen isotopes deep into the tungsten tile. The presence of hydrogen in tungsten tends to weaken the gram boundaries and produce surface blistering. The ejection of tungsten particulates into the plasma core will cool the reaction and cease the fusion reactions. Therefore, it is important to have a detailed understanding of the hydrogen diffusion in tungsten bulk. The KMC method is a stochastic computer simulation to predict the time evolution of a certain system. Typically, at a given time, the system will evolve through a few possible paths with known probabilities. In each step, the KMC method analyzes the entire system and randomly determines the system evolving events based on all transition rates of all possible events and advance the physical time. This process is repeated till the simulation of the entire physical time is completed. In previous work, the authors developed a kinetic Monte Carlo (KMC) algorithm and program to simulate hydrogen diffusion process on tungsten reconstructed surface. In this article, this KMC method is extended to model the hydrogen diffusion in tungsten bulk. The diffusion coefficients are calculated and fitted to the Arrhenius equation. The result matches the experiment measurement very well. KMC diffusion simulation needs the locations of the hydrogen trapping sites, diffusion paths and the corresponding diffusion energy barriers. The hydrogen interstitial sites in tungsten crystal have been thoroughly studied by density function theory (DFT) and molecular dynamics (MD). Johnson and Carter used DFT to determine the preferred location of the hydrogen interstitial sites m tungsten crystal, which are the tetrahedral sites (T-sites) in body-centered-cubic (bcc) lattice surfaces. All T-sites are located on the surfaces of a lattice cell. There are four T-sites per surface and twelve T-sites per lattice cell. The hydrogen atoms diffuse by traveling from one T-site to one of its nearest four T-site neighbors. Thus, each T-site is surrounded by two neighbors on the same surface, and one neighbor on each of the two adjacent surfaces perpendicular to the surface holding the center T-site. The diffusion barrier between two neighboring T-sites was calculated to be 0.38 eV. Knowing the diffusion path, energy barrier, and the crystal structure, the KMC simulation is readily to be performed. Based on other researchers' DFT calculations, the hydrogen stable interstitial site in tungsten bulk is determined as the tetrahedral site on the surface of the bee lattice. Each T-site has 4 neighbors and the energy barrier between two adjacent T-sites are estimated to be 0.38 eV. Based on those facts and the KMC theory, a MATLAB program is developed to model the hydrogen diffusion in tungsten bulk. In the MATLAB program, substrate temperature, hydrogen concentration, physical time to be simulated, and the number of cycles are all customizable. To compare the KMC determined diffusion coefficients with experiment results in the literature, the simulations are performed for temperature 300 K - 2500 K and low hydrogen concentration (0.1%), mimicking the experiment condition. The calculated diffusion coefficients match the reference values measured by the experiment under low hydrogen concentration condition. The results are also fitted to the Arrhenius equation. The pre-exponential factor and the average activation energy in Arrhenius equation are obtained. The good agreement between simulation and experiment proves that the developed KMC algorithm truly reflect the well-established theory of interstitial hydrogens in tungsten. The developed program could easily be used to model the diffusion process and predict the diffusion coefficients under other concentrations and temperatures. This work is part of the ongoing development of hydrogen diffusion in tungsten integrating surface diffusion (previous work), bulk diffusion (this work), and adsorption/desorption on surface (future work)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Social Sensor Analytics: Making Sense of Network Models in Social Media

Social networks can be thought of as noisy sensor networks mapping real world information to the web. Owing to the extensive body of literature in sensor network analysis, this work sought to apply several novel and traditional methods in sensor network analysis for the purposes of efficiently interrogating social media data streams from raw data. We carefully revisit our definition of a social media signal from previous work both in terms of time-varying features within the data and the networked nature of the medium. Further, we detail our analysis of global patterns in Twitter over the months of November 2013 and June 2014, detect and categorize events, and illustrate how these analyses can be used to inform graph-based models of Twitter, namely using a recent network influence model called PhySense: similar to PageRank but tuned to behavioral analysis by leveraging a sociologically inspired probabilistic model. We ultimately identify forms of information dissemination via analysis of time series and dynamic graph spectra and corroborate these findings through manual investigation of the data as a requisite step in modeling the diffusion process with PhySense. We hope to sufficiently characterize global behavior in a medium such as Twitter as a means of learning global model parameters one may use to predict or simulate behavior on a large scale. We have made our time series and dynamic graph analytical code available via a GitHub repository https://github.com/cpatdowling/salsa and our data are available upon request.

event detection↗
NASA NTRS1974

Modulation of cosmic-ray protons and helium nuclei near solar maximum

Balloon observations of proton and helium spectra in 1970, 1971 and 1972, reveal a factor of approximately 1.85 deviation from a single valued regression at low rigidities. This deviation decreases with increasing rigidity for both species. The period 1969-1970 is unique because time variations at low and high energies were anti-correlated. When satellite observations are used to extend the balloon observations to energies below 100 MeV/Nucleon, the proton spectrum showed a steeper slope in 1970 and 1972 than the characteristic J = AT spectrum observed during 1965-1969. The slope of the helium spectrum became continuously flatter during the same period (1970-1972). Computer generated spectra based on simple two parameter modulation models describe the basic features of the observations if one of the variable parameters is used to characterize the rigidity dependence of the diffusion coefficient. Models which do not allow such a variation are not consistent with the observations.

Rygg, T. A.↗
NASA NTRS1992

Flow effects in a vertical CVD reactor

A model is presented to simulate the non-Boussinesq flow in a vertical, two-dimensional, chemical vapor deposition reactor under atmospheric pressure. Temperature-dependent conductivity, mass diffusivity, viscosity models, and reactive species mass transfer to the substrate are incorporated. In the limits of small Mach number and small aspect ratio, asymptotic expressions for the flow, temperature, and species fields are developed. Soret diffusion effects are also investigated. Analytical solutions predict an inverse relationship between temperature field and concentration field due to Soret effects. This finding is consistent with numerical simulations, assisting in the understanding of the complex interactions amongst the flow, thermal, and species fields in a chemically reacting system.

Young, G. W.↗
DOE OSTI2016

Detailed high-resolution three-dimensional simulations of OMEGA separated reactants inertial confinement fusion experiments

Here, we present results from the comparison of high-resolution three-dimensional (3D) simulations with data from the implosions of inertial confinement fusion capsules with separated reactants performed on the OMEGA laser facility. Each capsule, referred to as a “CD Mixcap,” is filled with tritium and has a polystyrene (CH) shell with a deuterated polystyrene (CD) layer whose burial depth is varied. In these implosions, fusion reactions between deuterium and tritium ions can occur only in the presence of atomic mix between the gas fill and shell material. The simulations feature accurate models for all known experimental asymmetries and do not employ any adjustable parameters to improve agreement with experimental data. Simulations are performed with the RAGE radiation-hydrodynamics code using an Implicit Large Eddy Simulation (ILES) strategy for the hydrodynamics. We obtain good agreement with the experimental data, including the DT/TT neutron yield ratios used to diagnose mix, for all burial depths of the deuterated shell layer. Additionally, simulations demonstrate good agreement with converged simulations employing explicit models for plasma diffusion and viscosity, suggesting that the implicit sub-grid model used in ILES is sufficient to model these processes in these experiments. In our simulations, mixing is driven by short-wavelength asymmetries and longer-wavelength features are responsible for developing flows that transport mixed material towards the center of the hot spot. Mix material transported by this process is responsible for most of the mix (DT) yield even for the capsule with a CD layer adjacent to the tritium fuel. Consistent with our previous results, mix does not play a significant role in TT neutron yield degradation; instead, this is dominated by the displacement of fuel from the center of the implosion due to the development of turbulent instabilities seeded by long-wavelength asymmetries. Through these processes, the long-wavelength asymmetries degrade TT yield more than the DT yield and thus bring DT/TT neutron yield ratios into agreement with experiment. Finally, we present a detailed comparison of the flows in 2D and 3D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Current-induced spin torque resonance of magnetic insulators affected by field-like spin-orbit torques and out-of-plane magnetizations

The spin-torque ferromagnetic resonance (ST-FMR) in a bilayer system consisting of a magnetic insulator such as Y{sub 3}Fe{sub 5}O{sub 12} and a normal metal with spin-orbit interaction such as Pt is addressed theoretically. We model the ST-FMR for all magnetization directions and in the presence of field-like spin-orbit torques based on the drift-diffusion spin model and quantum mechanical boundary conditions. ST-FMR experiments may expose crucial information about the spin-orbit coupling between currents and magnetization in the bilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Modeling the Formation of Hydrate-Filled Veins in Fine-Grained Sediments from in Situ Microbial Methane

Continental margin sediments are dominantly fine-grained silt and clay, and methane hydrates in these sediments are often found in semi-vertical veins and fractures. In several instances, these hydrate veins occupy discrete depth intervals that are a few tens of meters thick and are surrounded by hydrate-free sediments. As they are not connected with gas sources beneath the base of the gas hydrate stability zone (GHSZ), these isolated hydrate-bearing intervals have been interpreted as formed by in situ microbial methane. To investigate the formation of these hydrate deposits, we applied a time-dependent advection-diffusion-reaction model that includes the effects of sedimentation, compaction, solute diffusion, and microbial methane generation. Microbial methane generation depends on the amount of metabolizable organic carbon deposited at the seafloor, whose progressive degradation produces methane beneath the sulfate reduction zone. If the amount of organic carbon entering the methanogenic zone is kept constant in time, we found that the computed amounts of hydrate formed in discrete intervals within the GHSZ are well below those estimated from observations. On the other hand, if the deposition of organic carbon is higher in a given time interval, methane generation during burial is more intense in the corresponding sediment interval, resulting in enhanced hydrate formation. With variations in organic carbon deposition comparable to those generally observed in continental margins, our model was able to reproduce the methane hydrate contents that were estimated from drilling. These results support the suggestion that in situ microbial generation associated with transient organic carbon deposition is the source of methane that forms isolated intervals of hydrate-filled veins in fine-grained sediments.

58 GEOSCIENCES↗
DOE OSTI2026

Influence of Rigidity–Hydration Coupling on Size-Dependent Diffusion in Hydrated Polymer Membranes

Selective ion transport in polymer membranes depends critically on how penetrant motion couples to polymer dynamics and hydration. Yet, the mechanistic interplay between polymer rigidity, water content, and penetrant size remains poorly understood, especially in the regime where the penetrant diameter, polymer Kuhn length, and correlation length are comparable. Here, we employ coarse-grained molecular dynamics simulations to systematically investigate penetrant diffusion in hydrated polymer networks across a broad range of water volume fractions, chain rigidities, and penetrant sizes. The results reveal a transition from a decoupled regime, where small penetrants diffuse nearly independently of polymer relaxation, to a coupled regime in which large penetrants require cooperative polymer motion for transport. Increasing polymer rigidity amplifies the sensitivity of diffusivity to hydration, particularly at low water content, leading to pronounced deviations from Stokes−Einstein scaling. Comparison with scaling theories and free-volume models shows that classical nanoparticle-based frameworks fail to capture this intermediate regime. To address this gap, we extend the Yasuda model to incorporate polymer rigidity through a single parameter that quantifies the dynamic contribution of chain stiffness to free-volume fluctuations. The resulting model collapses diffusivity data across all sizes, water contents, and rigidities, providing a unified description of penetrant transport in hydrated polymer matrices. Furthermore, these findings establish polymer rigidity as a key, tunable determinant of diffusion and offer a framework for interpreting size-dependent transport in ion-selective membranes.

diffusion↗
NASA NTRS1985

Onion-shell model for cosmic ray electrons and radio synchrotron emission in supernova remnants

The spectrum of cosmic ray electrons, accelerated in the shock front of a supernova remnant (SNR), is calculated in the test-particle approximation using an onion-shell model. Particle diffusion within the evolving remnant is explicity taken into account. The particle spectrum becomes steeper with increasing radius as well as SNR age. Simple models of the magnetic field distribution allow a prediction of the intensity and spectrum of radio synchrotron emission and their radial variation. The agreement with existing observations is satisfactory in several SNR's but fails in other cases. Radiative cooling may be an important effect, especially in SNR's exploding in a dense interstellar medium.

Beck, R.↗
NASA NTRS1990

Model of Jovian F region ionosphere (Jovian ionosphere model in offset dipole approximation)

The geomagnetic control of the Earth's atmosphere is well understood. In the F-region and the topside ionosphere, non-electrical forces transport plasma along the magnetic field lines only. In consequence, the worldwide distribution of ionization is strongly dependent on the dip angle. For example, the equatorial anomaly is roughly symmetrical about the dipole equator rather than the geographic. The same appears to be the case in the Jovian ionosphere (Mahajan, 1981). The influence of the magnetic field of Jupiter on its ionization pattern is one of several outstanding topics which need to be studied. Tan (1986) investigated the formation of the equatorial anomaly in the Jovian ionosphere under a centered dipole model. Tan (1988) further studied the effect of the tilt of the Jovian dipole. The results were in broad agreement with those of a diffusive equilibrium model (Tan and Wu, 1981). An off-centered dipole model is constructed and its effects on the ionization pattern are investigated.

Tan, A.↗
NASA NTRS1996

Migrational Instabilities in Particle Suspensions

This work deals with an instability arising from the shear-induced migration of particles in dense suspensions coupled with a dependence of viscosity on particle concentration. The analysis summarized here treats the inertialess (Re = O) linear stability of homogeneous simple shear flows for a Stokesian suspension model of the type proposed by Leighton and Acrivos (1987). Depending on the importance of shear-induced migration relative to concentration-driven diffusion, this model admits short-wave instability arising from wave-vector stretching by the base flow and evolving into particle-depleted shear bands. Moreover, this instability in the time-dependent problem corresponds to loss of ellipticity in the associated static problem (Re = O, Pe = O). While the isotropic version of the Leighton-Acrivos model is found to be stable with their experimentally determined parameters for simple shear, it is known that the stable model does not give a good quantitative description of particle clustering in the core of pipe flow (Nott and Brady 1994). This leads to the conjecture that an appropriate variant on the above model could explain such clustering as a two-phase bifurcation in the base flow.

Goddard, Joe D.↗
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↗