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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Electron acceleration in laboratory-produced turbulent collisionless shocks

Astrophysical collisionless shocks are among the most powerful particle accelerators in the Universe. Generated by violent interactions of supersonic plasma flows with the interstellar medium, supernova remnant shocks are observed to amplify magnetic fields and accelerate electrons and protons to highly relativistic speeds. In the well-established model of diffusive shock acceleration, relativistic particles are accelerated by repeated shock crossings. However, this requires a separate mechanism that pre-accelerates particles to enable shock crossing. This is known as the ‘injection problem’, which is particularly relevant for electrons, and remains one of the most important puzzles in shock acceleration. In most astrophysical shocks, the details of the shock structure cannot be directly resolved, making it challenging to identify the injection mechanism. In this letter we report results from laser-driven plasma flow experiments, and related simulations, that probe the formation of turbulent collisionless shocks in conditions relevant to young supernova remnants. We show that electrons can be effectively accelerated in a first-order Fermi process by small-scale turbulence produced within the shock transition to relativistic non-thermal energies, helping overcome the injection problem. Our observations provide new insight into electron injection at shocks and open the way for controlled laboratory studies of the physics underlying cosmic accelerators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Structural differences between human and mouse neurons and their implementation in generative AIs

Mouse and human brains have different functions that depend on their neuronal networks. We analyzed nanometer-scale three-dimensional structures of brain tissues of the mouse medial prefrontal cortex and compared them with structures of the human anterior cingulate cortex. The obtained results indicated that mouse neuronal somata are smaller and neurites are thinner than those of human neurons. We implemented these characteristics of mouse neurons in convolutional layers of a generative adversarial network (GAN) and a denoising diffusion implicit model (DDIM), which were then subjected to image generation tasks using photo datasets of cat faces, cheese, human faces, birds, and automobiles. The mouse-mimetic GAN outperformed a standard GAN in the image generation task using the cat faces and cheese photo datasets, but underperformed for human faces and birds. The mouse-mimetic DDIM gave similar results, suggesting that the nature of the datasets affected the results. Analyses of the five datasets indicated differences in their image entropy, which should influence the number of parameters required for image generation. The preferences of the mouse-mimetic AIs coincided with the impressions commonly associated with mice. The relationship between the neuronal network and brain function should be investigated by implementing other biological findings in artificial neural networks.

generative AI↗
DOE OSTI2019

Linear-scaling ab initio simulations of spin diffusion in rotating solids

We investigated the utility of locally restricting the basis sets involved in low-order correlations in Liouville space (LCL) calculations of spin diffusion. Using well-known classical models of spin diffusion, we describe a rationale for selecting the optimal basis set for such calculations. We then show that the use of these locally restricted basis sets provides the same computational accuracy as the full LCL set while reducing the computational time by several orders of magnitude. Speeding up the calculations also enables us to use higher maximum spin orders and increase the computational accuracy. Here, unlike exact and full LCL calculations, locally restricted LCL calculations scale linearly with the system size and should thus enable the ab initio study of spin diffusion in spin systems containing several thousand spins.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Parallel algorithms for hyperdynamics and local hyperdynamics

Hyperdynamics (HD) is a method for accelerating the timescale of standard molecular dynamics (MD). It can be used for simulations of systems with an energy potential landscape that is a collection of basins, separated by barriers, where transitions between basins are infrequent. HD enables the system to escape from a basin more quickly while enabling a statistically accurate renormalization of the simulation time, thus effectively boosting the timescale of the simulation. In [Kim, Perez, Voter, J Chem Phys, 139:144110, 2013)1, a local version of HD was formulated, which exploits the intrinsic locality characteristic typical of most systems to mitigate the poor scaling properties of standard HD as the system size is increased. In this paper, we discuss how both HD and local HD can be formulated to run efficiently in parallel. We have implemented these ideas in the LAMMPS MD code, which means HD can be used with any interatomic potential LAMMPS supports. Together, these parallel methods allow simulations of any size to achieve the time acceleration offered by HD (which can be orders of magnitude), at a cost 3-5x that of standard MD. As examples, we performed two simulations of a million-atom system to model the diffusion and clustering of Pt adatoms on a large patch of Pt(100) surface for 80 and 160 μs.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

The “good,” the “bad,” and the “hidden” in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl 2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca 2+ –Cl - interactions were primarily “hidden” in the Ca-NDIS signal due to overlapping with Ca 2+ –O w and Ca 2+ –H w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models “on average,” the question “which model is the best” has not been answered but replaced by the question “which model is better for a given research.” An overall “good” model can still be inappropriate or, in some instances, “bad” and, unfortunately, produce erroneous results. Additionally, the accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

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

Electron heating in 2D: Combining Fermi–Ulam acceleration and magnetic-moment non-adiabaticity in a mirror-configuration plasma

We analyze a new mechanism for the creation and confinement of energetic electrons in a mirror-configuration plasma. A Fermi–Ulam-type process, driven by end-localized coherent electrostatic oscillations, provides axial acceleration, while a natural non-adiabaticity of μ provides phase decorrelation and energy isotropization. This novel 2D combination causes the electron energy distribution function, calculated with a diffusive-loss model, to assume a Maxwellian shape with the μ non-adiabaticity, reducing loss-cone escape and annulling the absolute-barrier energy-limiting Chirikov criterion of lower dimensional models. Here, the theoretical predictions are compared with data from an experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electron kinetics in a positive column of AC discharges in a dynamic regime

Abstract We have performed hybrid kinetic-fluid simulations of a positive column in alternating current (AC) argon discharges over a range of driving frequencies f and gas pressure p for the conditions when the spatial nonlocality of the electron energy distribution function (EEDF) is substantial. Our simulations confirmed that the most efficient conditions of plasma maintenance are observed in the dynamic regime when time modulations of mean electron energy (temperature) are substantial. The minimal values of the root mean square electric field and the electron temperature have been observed at f/p values of about 3 kHz Torr −1 in a tube of radius R = 1 cm. The ionization rate and plasma density reached maximal values under these conditions. The numerical solution of a kinetic equation allowed accounting for the kinetic effects associated with spatial and temporal nonlocality of the EEDF. Using the kinetic energy of electrons as an independent variable, we solved an anisotropic tensor diffusion equation in phase space. We clarified the role of different flux components during electron diffusion in phase space over surfaces of constant total energy. We have shown that the kinetic theory uncovers a more exciting and rich physics than the classical ambipolar diffusion (Schottky) model. Non-monotonic radial distributions of excitation rates, metastable densities, and plasma density have been observed in our simulations at pR > 6 Torr cm. The predicted off-axis plasma density peak in the dynamic regime has never been observed in experiments so far. We hope our results stimulate further experimental studies of the AC positive column. The kinetic analysis could help uncover new physics even for such a well-known plasma object as a positive column in noble gases.

Physics↗
DOE OSTI2025

Using convolutional neural networks to accelerate three-dimensional coherent synchrotron radiation computations

Calculating the effects of coherent synchrotron radiation (CSR) is one of the most computationally expensive tasks in accelerator physics. Here, we use convolutional neural networks (CNNs), along with a latent conditional diffusion (LCD) model, trained on physics-based simulations to speed up calculations. Specifically, we produce the 3D CSR wakefields generated by electron bunches in circular orbit in the steady-state condition. Two datasets are used for training and testing the models: wakefields generated by three-dimensional Gaussian electron distributions and wakefields from a sum of up to 25 three-dimensional Gaussian distributions. The CNNs are able to accurately produce the 3D wakefields ∼250–1000 times faster than the numerical calculations, while the LCD achieves a gain of a factor of ∼34. We also test the extrapolation and out-of-distribution generalization ability of the models. They generalize well on distributions with larger spreads than what they were trained on but struggle with smaller spreads.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Localized basis set for plutonium

The implementation of optimal strictly localized atomic orbitals basis for plutonium (Pu) using norm-conserving pseudopotential density-functional theory (DFT) is presented. The basis set was applied to the α, β, γ, δ, δ', and ε phases of Pu, δ-Pu surface, and δ-PuGa alloys. The computed properties of the Pu phases and δ-Pu surface were in good agreement with both available experimental data and prior DFT calculations based on plane-wave methodologies. Results for the δ-PuGa alloys were also in good agreement with experimental data. Finally, the reliability of the basis set was further demonstrated by using ab initio molecular dynamics to model the diffusion coefficient and activation barrier for atomic diffusion in a δ-PuGa alloy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Millisecond timescale reactions observed via X-ray spectroscopy in a 3D microfabricated fused silica mixer

Determination of electronic structures during chemical reactions remains challenging in studies which involve reactions in the millisecond timescale, toxic chemicals, and/or anaerobic conditions. Here, a three-dimensionally (3D) microfabricated microfluidic mixer platform that is compatible with time-resolved X-ray absorption and emission spectroscopy (XAS and XES, respectively) is presented. This platform, to initiate reactions and study their progression, mixes a high flow rate (0.50–1.5 ml min –1 ) sheath stream with a low-flow-rate (5–90 µl min –1 ) sample stream within a monolithic fused silica chip. The chip geometry enables hydro­dynamic focusing of the sample stream in 3D and sample widths as small as 5 µm. The chip is also connected to a polyimide capillary downstream to enable sample stream deceleration, expansion, and X-ray detection. In this capillary, sample widths of 50 µm are demonstrated. Further, convection–diffusion-reaction models of the mixer are presented. The models are experimentally validated using confocal epifluorescence microscopy and XAS/XES measurements of a ferricyanide and ascorbic acid reaction. The models additionally enable prediction of the residence time and residence time uncertainty of reactive species as well as mixing times. Residence times (from initiation of mixing to the point of X-ray detection) during sample stream expansion as small as 2.1 ± 0.3 ms are also demonstrated. Importantly, an exploration of the mixer operational space reveals a theoretical minimum mixing time of 0.91 ms. The proposed platform is applicable to the determination of the electronic structure of conventionally inaccessible reaction intermediates.

3D microfabrication↗
DOE OSTI2023

Lorentz factor for time-of-flight neutron Bragg and total scattering

We report the three fundamental origins of the Lorentz factor for neutron time-of-flight powder diffraction are revisited. A detailed derivation of the Lorentz factor is presented in the context of diffuse scattering modelling in reciprocal space when perfect periodicity is assumed, and the total scattering pattern is constructed in its discrete form – the factor in this case becomes 1/Q 2 (or d 2 ). Discussion is also presented with respect to practical data reduction where a vanadium measurement is usually taken as the normalization factor (to account for various factors such as detector efficiency), and it is shown that the existence of the Lorentz factor is independent of such a normalization process.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Nonisothermal Free-Molecular Flow of Gas in an Elliptic Channel with a Circular Cylindrical Element Inside

The linearized problem of free-molecular gas flow in a long elliptic channel with a circular cylindrical element inside has been solved in the kinetic approximation. The flow in the channel is caused by temperature and pressure drops between its ends. The Boltzmann kinetic equation for collisionless gas has been used as a basic equation, and the boundary condition has been set in terms of the diffuse reflection model. The distribution of the mass velocity of the gas over the cross section of the channel has been obtained. The mass flow rate of the gas in the channel versus the temperature and pressure drops between its ends has been calculated. It has been found that the mass flow of the gas substantially depends on the radius of the inner cylinder.

42 ENGINEERING↗
DOE OSTI2018

Avoiding Pitfalls in the Determination of Reliable Electrochemical Kinetics Parameters for the Cu 2+ → Cu 1+ Reduction Reaction in Deep Eutectic Solvents

Deep eutectic solvents (DES), on account of their low-cost, non-flammability and electrochemical stability, are attracting attention for their potential use in applications such as electrodeposition and energy storage. In these applications, knowledge of the transport and electrochemical kinetics properties of DES is critically important. To date, attempts to measure the kinetics parameters of the Cu 2+ + e ⇌ Cu 1+ reaction in ethaline DES have yielded a cathodic charge transfer coefficient (α) in the range of 0.2–0.3 suggesting an unexpected asymmetric polarization behavior. In the present work, we pursued a comprehensive study of the kinetics and transport properties of the aforementioned reaction. Using steady-state and transient polarization measurements on RDE and microelectrodes combined with diffusion-reaction modeling, we demonstrate that the Cu 2+ /Cu 1+ transition exhibits a charge transfer coefficient in the range of 0.49–0.54 and an exchange current density in the range 1.72–1.88 mA/cm 2 . Recommendations are provided for avoiding pitfalls in kinetics analysis of highly resistive DES electrolytes. These recommendations include the application of accurate $IR_Ω$ correction and the use of electrode configurations where the secondary current distribution is relatively uniform.

25 ENERGY STORAGE↗
DOE OSTI2023

Tools for Assessing Performance (TAP) Project: FY2023 Quarter 2 Report

LANL personnel presented TAP-sponsored work at the 103 rd Annual Meeting of the American Meteorological Society, which was held in Denver, CO in January 2023. There were four different presentations related to TAP-sponsored work that were given by LANL staff and studentsat this meeting: (1) Gloeckler, T., M. Nelson, H. Shah, M. Deshler, L. Wedell, R. Linn, N. Duboc, S. Pol, P. Bieringer, A. Annunzio, B. Martin, and H. Tinnesand, 2023. Testing Validity and Limitations of a Low-Order Diffusive Wake Model. 15 th Symposium on the Urban Environment, January 8-12, 2023, Denver, CO, U.S.A. (2) Nelson, M., H. Shah, T. Gloeckler, M. Deshler, L. Wedell, A. Ortega, R. Linn, N. Duboc, P. Conry, S. Pol, and H. Tinnesand, 2023. The Isolated Building Wake Experiment. Paper 1A.1. 24 th Symposium on Boundary Layers and Turbulence, January 8-12, 2023. Denver, CO, U.S.A. (3) Shah, H., M. Nelson, S. Pol, R. Linn, N. Duboc, and H. Tinnesand, 2023. Detailed Measurement of a Building Wake Under Real-World Atmospheric Conditions. Poster 102. 24 th Symposium on Boundary Layers and Turbulence, January 8-12, 2023, Denver, CO, U.S.A. (4) Wedell, L., M. Nelson, H. Shah, T. Gloeckler, M. Deshler, R. Linn, S. Pol, H. Tinnesand, 2023. The Isolated Building Wake Experiment: Using a Rotating Building to Mulitply the Number of Wake Measurement Locations without Moving the Instruments. Paper 1A.4A. 24 th Symposium on Boundary Layers and Turbulence, January 8-12, 2023, Denver, CO, U.S.A.

54 ENVIRONMENTAL SCIENCES↗