Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic simulations”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

How eigenmode self-interaction affects zonal flows and convergence of tokamak core turbulence with toroidal system size

Self-interaction is the process by which a microinstability eigenmode that is extended along the direction parallel to the magnetic field interacts non-linearly with itself. This effect is particularly significant in gyrokinetic simulations accounting for kinetic passing electron dynamics and is known to generate stationary $E\times B$ zonal flow shear layers at radial locations near low-order mode rational surfaces (Weikl et al. Phys. Plasmas, vol. 25, 2018, 072305). Here, we find that self-interaction, in fact, plays a very significant role in also generating fluctuating zonal flows, which is critical to regulating turbulent transport throughout the radial extent. Unlike the usual picture of zonal flow drive in which microinstability eigenmodes coherently amplify the flow via modulational instabilities, the self-interaction drive of zonal flows from these eigenmodes are uncorrelated with each other. It is shown that the associated shearing rate of the fluctuating zonal flows therefore reduces as more toroidal modes are resolved in the simulation. In simulations accounting for the full toroidal domain, such an increase in the density of toroidal modes corresponds to an increase in the toroidal system size, leading to a finite system size effect that is distinct from the well-known profile shearing effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Ethanol Conversion over La 0.7 Sr 0.3 MnO 3–x (100): Autocatalysis, Adjacent O-Vacancies, Disproportionation, and Dehydrogenation

The mechanism for catalytic conversion of ethanol over La 0.7 Sr 0.3 MnO 3– x (100) surface to acetaldehyde and ethene was investigated. Pre-exposure temperature-programmed reaction (PE-TPR) experiments were performed in which ethanol was introduced to oxidized or reduced surfaces followed by heating. In particular, sequential PE-TPR experiments were conducted to incrementally and gradually reduce the surface. The products and their ratios were investigated as a function of surface reduction. The data show that acetaldehyde and ethene production is catalyzed with hydrogen abstraction and oxygen abstraction reactions occurring by intermediates in vacancies at various temperatures >400 K. Adsorption of acetaldehyde followed by a temperature-programmed reaction does not produce ethene, indicating that acetaldehyde is not an intermediate to ethene and that the hydrogen and oxygen abstraction from ethanol to ethene are decoupled steps. Further evidence for this mechanistic nuance was obtained using isotopically labeled ethanol (CD 3 CH 2 OH), which produces CD 3 CHO and CD 2 CH 2 . Additionally, the ratio of aldehyde production to alkene production increases with reduction, suggesting that aldehyde is produced from a disproportionation reaction between ethoxy species in adjacent O-vacancies, while ethene is produced from a dehydrogenation reaction with ethoxy species in vacancies without requiring adjacent O-vacancies. Counterintuitively, this finding indicates that the more oxygenated product (aldehyde vs ethene) is favored with more vacancies and that the net alcohol conversion is autocatalytic. Density functional theory calculations were able to find the previously unknown disproportionation pathway between ethoxies in adjacent O-vacancies, and kinetic Monte Carlo simulations support this interpretation by reproducing experimental selectivities. The activation energies for these pathways are estimated as 132 ± 10 kJ/mol for the disproportionation reaction (when occurring between ethoxies in adjacent vacancies) and as 148 ± 11 kJ/mol for the direct dehydrogenation reaction of an ethoxy in a vacancy. Based on these results, a mechanism with operative pathways based on elementary steps in O-vacancies is reported.

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

Kinetic Model of Acid Gas Induced Defect Propagation in Zeolitic Imidazolate Frameworks

Understanding the degradation of nanoporous materials under exposure to common acid gas contaminants (e.g., SO 2 , CO 2 , NO 2 , and H 2 S) is essential to elongate their lifetime and thus enable their practical applications in separations and catalysis. Previous theoretical investigations have focused on the formation of isolated point defects, which are insufficient to provide direct insights into the long-term evolution of the bulk properties of materials such as zeolitic imidazolate frameworks (ZIFs) under sustained acid gas exposure. Here, to bridge this divide in both length and time scales, we developed a first-principles lattice-based kinetic model to simulate the defect propagation and bulk material breakdown in ZIFs. This model closely reproduces the experimentally measured macroscopic evolution of the time-dependent bulk materials proprieties and also yields important new insights regarding the autocatalytic nature of ZIF degradation and the spatial distribution of defects. Our results suggest new experimental directions to identify nascent defect clusters in degraded ZIFs and avenues to mitigate degradation under challenging conditions of acid gas exposure.

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

Nonlinear saturation of reverse shear Alfvén eigenmodes induced by kinetic electrons

A strong nonlinear saturation mechanism of the reverse shear Alfvén eigenmode (RSAE) induced by thermal electrons is observed in gyrokinetic $\delta$$f$ particle-in-cell simulation. This nonlinear effect occurs at moderate mode amplitude, $\delta$$B$ ⟂ /$B$ < 10 -3 , and is associated with the electron parallel streaming along the perturbed magnetic field lines, the term v $\|$ $\frac{δΒ_{⟂}}{B}$ ∙ ∇$\delta$$f$, in the electron drift-kinetic equation. In the case of an n = 4 RSAE, where n is the toroidal mode number, this magnetic fluttering nonlinearity leads to strong mode coupling and broadening of the saturation spectrum. The nonlinear components, which are highly damped through the interaction between thermal electrons and the parallel electric field, modifies the RSAE mode structure and suppresses energetic particle drive significantly. The split-weight scheme (Chen and Parker 2007 J. Comput. Phys. 220 839–55) is used to simulate kinetic electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Atomic kinetics of laboratory photoionized plasmas relevant to astrophysics. Final report

The major goal of this project was to study the atomic and radiation physics of photoionized plasmas driven by a broadband x-ray flux and its impact on plasma x-ray heating and ionization, emissivity and opacity, and radiation-hydrodynamics. This problem is relevant to the physics of accretion disks surrounding black holes, x-ray binaries, and active galactic nuclei in astrophysics. In addition, it is also important for understanding non-equilibrium laboratory plasmas driven by a distribution of photons. To this end, and motivated by previous and ongoing gas-cell experiments performed by the PI at the Z facility of Sandia National Laboratories, we developed a new experimental platform at the 1MA Zebra pulsed power accelerator of the University of Nevada, Reno based on using a gas jet produced by a supersonic nozzle. The gas jet is turned into a photoionized plasma driven by the x- ray flux from a wire-array z-pinch. Establishing this new experimental platform at Zebra was a major goal of this project. In addition, we have also dedicated effort to continue complementary experiments at the Z facility using the gas-cell photoionized plasma set up with support from the Fundamental Science Program at Z as well as to perform data processing and analysis, and theory and modeling simulations, including electron kinetics and radiation-hydrodynamics. Since we seek to test and establish what physics models are needed to describe photoionized plasma guided by experimental observation, these complementary experimental efforts were central to achieving the goals of this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Shape Stability of Truncated Octahedral fcc Metal Nanocrystals

Metallic nanocrystals (NCs) can be synthesized with tailored nonequilibrium shapes to enhance desired properties, e.g., octahedral fcc metal NCs optimize catalytic activity associated with {111} facets. However, maintenance of optimized properties requires stability against thermal reshaping. Thus, we analyze the reshaping of truncated fcc metal octahedra mediated by surface diffusion using a stochastic atomistic-level model with energetic input parameters for Pd. Furthermore, the model describes NC thermodynamics by an effective nearest-neighbor interaction and includes a realistic treatment of diffusive hopping for undercoordinated surface atoms. Kinetic Monte Carlo simulation reveals that the effective barrier, E eff , for the initial stage of reshaping is strongly tied to the degree of truncation of the vertices in the synthesized initial octahedral shapes. This feature is elucidated via exact analytic determination of the energy variation along the optimal reshaping pathway at low-temperature (T), which involves transfer of atoms from truncated {100} vertex facets to form new layers on {111} side facets. Deviations from predictions of the low-T analysis due to entropic effects are more prominent for higher T and larger NC sizes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Applying the Gyrokinetic Formulation of Magnetized Plasma Dynamics to Z-pinch plasmas

This project adapts and applies the gyrokinetic Eulerian code COGENT to an emerging area of highly-magnetized Z-pinch plasma systems. Simulation efforts for these plasmas currently use either magneto-hydrodynamic (MHD) fluid codes, which suffer from fidelity being restricted to narrow plasma conditions, or fully kinetic particle-in-cell (PIC) methods, which suffer from the run-time needed to resolve the very small spatial and temporal scales required for the plasma conditions. The progress with achieving detailed theoretical understanding of Z-pinch plasma systems can be greatly expedited using some combination of gyrokinetic and drift-ordered MHD models, developed in the tokamak community and implemented in the COGENT code. These formulations are applicable to a much larger parameter space than standard MHD models, while providing a substantial speed-up in run-time over fully kinetic codes. In particular, the gyrokinetic formalism is successfully applied to a flow-stabilized Z pinch (FSZP) system – a promising fusion concept where a sheared axial flow is used to provide stability against large-scale ideal-MHD modes. Results of gyrokinetic simulations with COGENT agree well with those from fully kinetic simulations, while requiring ~500x less run-time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Exascale granular microstructure reconstruction in 3D volumes of arbitrary geometries with generative learning

Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method’s capability to efficiently handle complex geometries. In conclusion, the proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Radiation-reduced diffusion of Nd in bcc Fe

Here, this study combines ab initio calculations and kinetic Monte Carlo simulations to investigate the bulk diffusion of Nd, an exemplar lanthanide element, in bcc Fe. Contrary to the conventional belief that irradiation accelerates diffusion processes, our study reveals that irradiation can unexpectedly reduce the diffusion of Nd due to the strong binding between substitutional Nd atoms and vacancies in bcc Fe. Combined with the high concentration of nonequilibrium vacancies, this promotes the abundant formation of Nd+3Va defect complexes containing a single Nd substitutional atom and three vacancies, accompanied by the depletion of Nd-vacancy pairs that mediate the diffusion of Nd.

Density functional theory (DFT)↗
DOE OSTI2020

Turbulent kinetic energy over large offshore wind farms observed and simulated by the mesoscale model WRF (3.8.1)

Abstract. Wind farms affect local weather and microclimates;hence, parameterizations of their effects have been developed for numerical weather prediction models.While most wind farm parameterizations (WFPs) include drag effects of wind farms, models differ on whether or not an additionalturbulent kinetic energy (TKE) source should be included in theseparameterizations to simulate the impact of wind farms on the boundary layer.Therefore, we use aircraft measurements above large offshore wind farms in stable conditions to evaluateWFP choices. Of the three case studies we examine, we find the simulated ambient backgroundflow to agree with observations of temperature stratification and winds. This agreement allows us toexplore the sensitivity of simulated wind farm effects with respect to modeling choicessuch as whether or not to include a TKE source, horizontal resolution, vertical resolutionand advection of TKE. For a stably stratified marine atmospheric boundary layer (MABL), a TKE source and a horizontal resolutionon the order of 5 km or finer are necessary to represent the impact of offshore wind farms on the MABL.Additionally, TKE advection results in excessively reduced TKE over the wind farms, which in turn causesan underestimation of the wind speed deficit above the wind farm. Furthermore, using fine verticalresolution increases the agreement of the simulated wind speed with satellite observations of surface wind speed.

17 WIND ENERGY↗
DOE OSTI2024

Mineral Deposition on the Rough Walls of a Fracture

Modeling carbonate growth in fractures and pores is important for understanding carbon sequestration in the environment or when supersaturated solutions are injected into rocks. Here, we study the simple but nontrivial problem of calcite growth on fractures with rough walls of the same mineral using kinetic Monte Carlo simulations of attachment and detachment of molecules and scaling approaches. First, we consider wedge-shaped fracture walls whose upper terraces are in the same low-energy planes and show that the valleys are slowly filled by the propagation of parallel monolayer steps in the wedge sides. The growth ceases when the walls reach these low-energy configurations so that a gap between the walls may not be filled. Second, we consider fracture walls with equally separated monolayer steps (vicinal surfaces with roughness below 1 nm) and show that growth by step propagation will eventually clog the fracture gap. In both cases, scaling approaches predict the times to attain the final configurations as a function of the initial geometry and the step-propagation velocity, which is set by the saturation index. The same reasoning applied to a random wall geometry shows that step propagation leads to lateral filling of surface valleys until the wall reaches the low-energy crystalline plane that has the smallest initial density of molecules. Thus, the final configurations of the fracture walls are much more sensitive to the crystallography than to the roughness or the local curvature. The framework developed here may be used to determine those configurations, the times to reach them, and the mass of deposited mineral. Effects of transport limitations are discussed when the fracture gap is significantly narrowed.

calcite↗
DOE OSTI2022

Thermal ion kinetic effects and Landau damping in fishbone modes

We report the kinetic–magnetohydrodynamic (MHD) hybrid simulation approach for macroscopic instabilities in plasmas can be extended to include the kinetic effects of both thermal ions and energetic ions. The new coupling scheme includes synchronization of the density and parallel velocity between thermal ions and MHD, in addition to pressure coupling, to ensure the quasineutrality condition and avoid numerical errors. The new approach has been implemented in the kinetic-MHD code M3D-C1-K, and was used to study the thermal ion kinetic effects and Landau damping in fishbone modes in both DIII-D and NSTX. It is found that the thermal ion kinetic effects can cause an increase of the frequencies of the non-resonant n = 1 fishbone modes driven by energetic particles for q min >1, and Landau damping can provide additional stabilization effects. A nonlinear simulation for n = 1 fishbone mode in NSTX is also performed, and the perturbation on magnetic flux surfaces and the transport of energetic particles are calculated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Design, scaling and cost evaluations of circulating fluidized-bed systems for biomass pyrolysis

To generate updated and transparent capital cost estimates for biomass fast pyrolysis equipment, refinery fluidized catalytic cracking design and sizing principles are examined and extended to pyrolysis of woody biomass. Capital costs for the sized equipment are estimated with process-industry software. A one-dimensional flow simulation with pyrolysis kinetics is leveraged to validate the fluidization conditions and thermal energy balance. After successful sizing and a system cost estimate of $\$$2.8 M (in 2016 US$\$$) at the biorefinery scale of 1000 metric tons per day (MTD), these methods were exercised for even smaller equipment at the distributed pyrolysis scale with modifications to the process design constraints, and not directly comparable with the 1000 MTD case, arriving at capital cost estimates of $\$$1.2 M for a 500 ton/day system and $\$$0.9 M for a 200 ton/day system. Here in this work, It is noted that this work only estimates purchased equipment costs at the ±50% accuracy level; there are significant other custom factors applicable to each installation based on location, maturity, scale, complexities during installation, engineering and licensing costs, etc. that need to be added on to these estimates to derive investment costs.

09 BIOMASS FUELS↗
DOE OSTI2020

Non-compact oxide-island growth induced by surface phase transition of the intermetallic NiAl during vacuum annealing

Crystal structure and composition are inter-dependent and decoupling their effects on surface reactivity is challenging. Using low-energy electron microscopy to spatially and temporally resolve the oxide film growth during the oxidation of NiAl(100), we differentiate such coupled effects by monitoring oxide growth while simultaneously fine-tuning the surface structure and composition during oxidation. We demonstrate that the oxidation of chemically ordered surfaces results in compact oxide island growth whereas non-compact oxide growth during the surface phase transition. By incorporating the surface phase transition induced chemical disordering into kinetic Monte Carlo simulations, we illustrate that the non-compact oxide growth is induced by the composition effect on the surface diffusion of oxygen, which can be described by the concept of “ant in the labyrinth”.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Novel effects of grain size and ion implantation on grain boundary segregation in ion irradiated austenitic steel

Nuclear reactor environments are extreme in nature due to the combination of exposure to corrosive coolant, mechanical stresses, and neutron irradiation damage. This combination of environments can lead to irradiation induced/assisted precipitation, segregation, and embrittlement coupled with enhanced creep and stress corrosion cracking. While nanostructured materials are of interest in a variety of applications due to their unique properties including high mechanical strength, for nuclear applications nanostructured materials have the added benefit of irradiation resistance due to the high volume density of sinks for irradiation-induced defects. In this study, we show unique radiation-induced segregation (RIS) behavior in nanocrystalline 304 stainless steel (SS304). Typical RIS in austenitic steels includes Ni enrichment and Cr depletion within the grain boundaries. In nanostructured SS304 this RIS is reduced. Furthermore, interestingly, some grain boundaries in the nanocrystalline SS304 in this study were found to be enriched in Cr after Fe 2+ ion irradiation at 500 °C up to 50 displacements per atom (dpa). To understand this unique observation, lattice-based atomic kinetic Monte Carlo simulations were performed to investigate the effects of grain size, thermal segregation, as well as injected interstitials (introduced as Fe self-ions). The results indicate that the Cr enrichment at grain boundaries may be caused by combined effects of grain size, self-ion injection and preferential Cr diffusion via the interstitial mechanism. In conclusion, the findings also suggest a possible difference in RIS between ion and neutron irradiations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The neutron number probability distribution in coupled lumped assemblies

Here, the validity of the gamma distribution in describing the neutron number probability distribution function for both isolated and coupled multiplying assemblies when constrained to reproduce the true mean and variance is investigated in lumped geometry by numerical comparison with kinetic Monte Carlo simulations. The mean and variance are obtained from numerical solution of moment equations constructed from the relevant forward Master equation with assembly coupling coefficients obtained from a view factor method. Numerical results for a two-group, two coupled assemblies model, with static and dynamic reactivity insertion, show that except for subcritical assemblies, the gamma distribution well-approximates the number distribution. Differences in the fast and thermal neutron population shapes are explained in terms of effective source strengths due to downscatter and coupling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Role of Short-Range Order on Diffusion Coefficients in the Li–Mg Alloy

Li–Mg alloys are important because of their beneficial role in fostering uniform plating and stripping of lithium in all-solid-state batteries. The alloy Li x Mg 1–x forms a solid solution on the BCC crystal structure when the lithium content is greater than x ≈ 0.3. The activation barriers of lithium and magnesium exchanges with a vacancy, crucial for substitutional diffusion, are predicted to be exceptionally low and almost identical, with negligible dependence on the alloy composition. The equilibrium vacancy concentration at room temperature is predicted to be very low, and it also remains almost constant with no dependence on Mg content in the alloy (for x Li ≥ 0.5). Nevertheless, both experiments and kinetic Monte Carlo simulations indicate that the tracer diffusion coefficients decrease by almost an order of magnitude with the addition of Mg to the alloy. In this contribution, the crucial role that chemical short-range order plays in affecting the diffusion coefficients is studied. Chemical short-range order is found to increase the effective activation barrier for lithium and magnesium diffusion by making successive atomic hops with vacancies correlated.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantifying the Decomposition Kinetics of Linear C 2 –C 4 Perfluoroalkanes

The pyrolysis of C 2 –C 4 linear perfluoroalkanes was investigated using laser schlieren densitometry (LS) in a diaphragmless shock tube. Mixtures of 1, 2, and 4% perfluoroethane (PFE), perfluoropropane (PFP), and perfluorobutane (PFB) in excess krypton were shock-heated to 1400–2500 K at pressures of 60, 120, and 240 Torr. The initiation step for all perfluoroalkanes involved carbon–carbon bond fission. The rate constants for unimolecular decomposition were computed using microcanonical rate theory, and these rate constants were combined into a detailed chemical kinetic mechanism. For PFE and PFP, there was one dissociation channel, while PFB dissociated via two competing bond-fission reactions. The perfluoroalkyl radicals subsequently recombined or underwent further unimolecular decomposition. In conclusion, the kinetic model accurately simulated the density gradient profiles of all perfluoroalkanes.

Guzman, Eduardo H. [Brown Univ., Providence, RI (U↗