Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic simulation”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Implementation of a 9-point stencil in SOLPS-ITER and implications for Alcator C-Mod divertor plasma simulations

The SOLPS-ITER code suite is used worldwide for plasma edge modeling, the interpretation of experiments, as well as for the design of the ITER divertor. The numerical scheme of the plasma solver of the code, B2.5, is based on the assumption of perfectly field-aligned grids, while in practice grids are often strongly distorted to match divertor target shapes. Neglecting these grid distortion leads to qualitatively and quantitatively incorrect results for fluid neutral simulations, and may affect results in cold (detached) divertors even when using kinetic neutral simulations. In this contribution, we present the first results of a newly implemented 9-point stencil in B2.5 to properly handle misaligned grids. The new scheme is then applied to fluid neutral simulations of a well-diagnosed and previously modeled Alcator C-Mod discharge. Results are compared with the original 5-point scheme neglecting grid distortion effects, as well as with simulations including a full kinetic neutral model. We conclude that the 9-point stencil is essential to correctly model the transport of fluid neutrals on distorted grids, and to capture the effects of divertor closure on the fluid neutral behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling of a Generic Pebble Bed High-temperature Gas-cooled Reactor (PB-HTGR) with SAM

This report presents the modeling of the core of a generic pebble-bed reactor (PBR) at the system level using the System Analysis Module (SAM) code. This work is an extension of a previous work by the authors (Ooi et al. (2021)) that used the so-called 2-D ring model approach to model the PBMR-400. With the new approach, the pebble bed of the reactor is modeled with multiple PBCoreChannel components with spherical heat structures which allows the code to calculate thermal fluid parameters with built-in closure relations. The new core-channel approach is an improvement to the 2-D ring model approach as it does not introduce geometric distortions to the model and thus reduces the uncertainties of the predictions. In addition to thermal fluid simulations, point kinetics (PKE) are included to the model. Simulations are performed under a steady-state normal operation condition and a load-following transient scenario. This particular transient scenario is chosen as it tests both the thermal fluid and neutronics aspects of the model. The predicted results from both the steady-state and transient scenarios are compared with the results by Stew- art et al. (2021) who performed similar simulations with a Griffin-Pronghorn coupled tool. Despite the differences between the codes, with SAM being a system-analysis code and Pronghorn being a porous-medium code, both sets of results compare favorably. The over- all profiles and trends of the predicted temperatures and reactivities from the SAM and Griffin/Pronghorn simulations are similar, with some differences in their predicted values. The first part of the report covers the significance of a relatively fast-running approach that is capable of modeling the pebble bed reactor at the system-level while simultaneously capturing the radial thermal behavior of the core. Then, the modeling approach used in this work is discussed in details. Lastly, the results and comparisons with the Griffin/Pronghorn simulation are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2015

Next-Generation Ion Thruster Design Tool

Computational tools that accurately predict the performance of electric propulsion devices are highly desirable and beneficial to NASA and the broader electric propulsion community. The current state of the art in electric propulsion modeling relies heavily on empirical data and numerous computational "knobs." In Phase I of this project, Tech-X Corporation developed the most detailed ion engine discharge chamber model that currently exists. This kinetic model simulates all particles in the discharge chamber along with a physically correct simulation of the electric fields. In addition, kinetic erosion models are included for modeling the ion-impingement effects on thruster component erosion. In Phase II, Tech-X developed a user-friendly computer program for NASA and other governmental and industry customers. Tech-X has implemented a number of advanced numerical routines to bring the computational time down to a commercially acceptable level. NASA now has a highly sophisticated, user-friendly ion engine discharge chamber modeling tool.

Stolz, Peter↗
DOE OSTI2013

CURRENT SHEETS AND COLLISIONLESS DAMPING IN KINETIC PLASMA TURBULENCE

We present the first study of the formation and dissipation of current sheets at electron scales in a wave-driven, weakly collisional, three-dimensional kinetic turbulence simulation. We investigate the relative importance of dissipation associated with collisionless damping via resonant wave-particle interactions versus dissipation in small-scale current sheets in weakly collisional plasma turbulence. Current sheets form self-consistently from the wave-driven turbulence, and their filling fraction is well correlated to the electron heating rate. However, the weakly collisional nature of the simulation necessarily implies that the current sheets are not significantly dissipated via Ohmic dissipation. Rather, collisionless damping via the Landau resonance with the electrons is sufficient to account for the measured heating as a function of scale in the simulation, without the need for significant Ohmic dissipation. This finding suggests the possibility that the dissipation of the current sheets is governed by resonant wave-particle interactions and that the locations of current sheets correspond spatially to regions of enhanced heating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
DOE OSTI2017

In Situ TEM Multi-Beam Ion Irradiation as a Technique for Elucidating Synergistic Radiation Effects

Materials designed for nuclear reactors undergo microstructural changes resulting from a combination of several environmental factors, including neutron irradiation damage, gas accumulation and elevated temperatures. Typical ion beam irradiation experiments designed for simulating a neutron irradiation environment involve irradiating the sample with a single ion beam and subsequent characterization of the resulting microstructure, often by transmission electron microscopy (TEM). This method does not allow for examination of microstructural effects due to simultaneous gas accumulation and displacement cascade damage, which occurs in a reactor. Sandia’s in situ ion irradiation TEM (I3TEM) offers the unique ability to observe microstructural changes due to irradiation damage caused by concurrent multi-beam ion irradiation in real time. This allows for time-dependent microstructure analysis. A plethora of additional in situ stages can be coupled with these experiments, e.g., for more accurately simulating defect kinetics at elevated reactor temperatures. This work outlines experiments showing synergistic effects in Au using in situ ion irradiation with various combinations of helium, deuterium and Au ions, as well as some initial work on materials utilized in tritium-producing burnable absorber rods (TPBARs): zirconium alloys and LiAlO2.

Taylor, Caitlin↗
NASA NTRS1998

Rayleigh-Benard Simulation using Gas-Kinetic BGK Scheme in the Incompressible Limit

In this paper, a gas-kinetic BGK model is constructed for the Rayleigh-Benard thermal convection in the incompressible flow limit, where the flow field and temperature field are described by two coupled BGK models. Since the collision times and pseudo-temperature in the corresponding BGK models can be different, the Prandtl number can be changed to any value instead of a fixed Pr=1 in the original BGK model. The 2D Rayleigh-Benard thermal convection is studied and numerical results are compared with theoretical ones as well as other simulation results.

Xu, Kun↗
DOE OSTI2018

DOE Final Report: A Unified Understanding of Residual Stress in Thin Films: Kinetic Models, Experiments and Simulations

Thin films are critical for a wide range of advanced technologies. However, the deposited films often have high levels of residual stress that can limit their performance or lead to failure. The stress is known to depend on many variables, including the processing conditions, type of material, deposition technique and the film’s microstructure. The goal of this DOE program was to develop a fundamental understanding of how the different processes that control thin film growth under different conditions can be related to the development of stress. In the program, systematic experiments were performed or analyzed that related the stress to the processing conditions that were used. Measurements of stress were obtained for films that were grown at different rates, different solutions (for electrodeposition), different particle energies (for sputter deposition) and different microstructures. Based on this data, models were developed to explain the observed dependence on the different parameters. The models were based on considering the balance among different stress-inducing mechanism occurring as the film grows (for both non-energetic and energetic deposition). Comparison of the model predictions with the experiments enabled the kinetic parameters to be determined for different materials. The resulting model equations provide a comprehensive picture of how stress changes with the processing conditions that can be used to optimize the growth of thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2016

CFD Simulation of a 10 kW Chemical Looping with Oxygen Uncoupling System: Effects Of Process And Fuel Parameters

Chemical looping with oxygen uncoupling (CLOU) is a carbon-capture technology that utilizes a metal oxide as an oxygen carrier to separate oxygen from air. The metal oxide is oxidized in an air reactor and releases gaseous oxygen in a fuel reactor where a fuel, such as coal, is fed. CFD simulations will expedite the scale-up of CLOU systems to demonstration and utility size plants by resolving technical concerns and showing the CLOU technology at larger scales. In this study, Barracuda was used as the platform for simulation. The CLOU kinetics were incorporated into the simulation and the coal kinetics were explored for one-step versus two-step combustion kinetics. Two cases were run, one for a pure carbon-char and one with the inclusion of ash. Using the two-step kinetics and inclusion of ash in the simulation yielded reasonable temperatures in both the gas and particle phases. Also studied were the effects of heat input, via increasing the coal feed, and coal particle diameter. These simulations showed that the increasing of the coal heat input had little effect as compared to changing the particle diameter. Increased particle diameter increased the particle temperatures and decreased the carbon capture efficiency.

Hamilton, Matthew A.↗
DOE OSTI2018

A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons

A conservative scheme has been formulated and confirmed for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from the drift kinetic equation is used to time advance the electron density perturbation by using the perturbed mechanical flow calculated from the parallel vector potential, and the parallel vector potential is solved by using the perturbed canonical flow from the perturbed distribution function. In gyrokinetic particle simulations using this new scheme, the shear Alfvén wave dispersion relation in the shearless slab and continuum damping in the sheared cylinder have been recovered. The new scheme overcomes the stringent requirement in the conventional perturbative simulation method that perpendicular grid size needs to be as small as electron collisionless skin depth even for the long wavelength Alfvén waves. The new scheme also circumvents the problem in the typical approach that an unphysically large parallel electric field arises due to the inconsistency between electrostatic potential calculated from the perturbed density and vector potential calculated from the perturbed canonical flow. Lastly, the gyrokinetic particle simulations of the Alfvén waves in sheared cylinder have superior numerical properties compared with the fluid simulations, which suffer from numerical difficulties associated with singular mode structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Investigation of dynamics of ELM crashes and their mitigation techniques

The accurate prediction of H-mode pedestal dynamics is critical for planning experiments in existing tokamaks and in the design of future tokamaks such as ITER and DEMO. The main objective of the proposed research is to advance the understanding of the physics of H-mode pedestal. Through advances in coupled kinetic-MHD simulations, a new model for H-mode pedestal and ELM crashes as well as an improved model for the bootstrap current will be developed. ELMmitigation techniques will also be investigated. The proposed research will help design efficient confinement scenarios and reduce transient heat loads on the divertor and plasma facing components. During the last two years, the principal investigator (PI) of this proposal actively participated in physics studies related to the DOE Joint Research Targets. These studies include the modeling of divertor heat load in the DIII-D, Alcator C-Mod, and NSTX tokamaks in 2010, and the modeling of H-mode pedestal structure in the DIII-D tokamak in 2011. It is proposed that this close collaboration with experimentalists from major US tokamaks continue during the next funding period. Verification and validation will be a strong component of the proposed research. During the course of the project, advances will be made in the following areas; Dynamics of the H-mode pedestal buildup and recovery after ELM crashes – The effects of neutral fueling, particle and thermal pinches will be explored; Dynamics of ELM crashes in realistic tokamak geometries – Heat loads associated with ELM crashes will be validated against experimental measurements. An improved model for ELM crashes will be developed; ELM mitigation – The effect of resonant magnetic perturbations on ELMs stability and their evolution will be investigated; Development of a new bootstrap current model – A reduced model for will be developed through careful verification of existing models for bootstrap current against first-principle kinetic neoclassical simulations. The primary focus here will be on plasma regimes where the existing reduced models deviate from each other, such as high-beta and low aspect ratio plasmas; and Results will be verified against analytical models and simulation results from other codes. Simulation results will be compared with diagnostic experimental data such as Thomson measurements (HRTC) or flux measurements if available. We have arrangements with experimentalists from DIII-D and JET to use diagnostic edge data for this validation. In summary, we propose to study the H-mode pedestal buildup, the stability and evolution of ELMs using a model that will account for dynamically evolving neoclassical sources.

42 ENGINEERING↗
DOE OSTI2026

Scaling kinetic Monte-Carlo simulations of grain growth with combined convolutional and graph neural networks

Graph neural networks (GNN) have emerged as a promising machine learning method for microstructure simulations such as grain growth. However, accurate modeling of realistic grain boundary networks requires large simulation cells, which GNN has difficulty scaling up to. To alleviate the computational costs and memory footprint of GNN, we suggest a hybrid architecture combining a convolutional neural network (CNN) based bijective autoencoder to compress the spatial dimensions, and a GNN that evolves the microstructure in the latent space of reduced spatial sizes. Our results demonstrate that the new design significantly reduces computational costs with using fewer message passing layer (from 12 down to 3) compared with GNN alone. The reduction in computational cost becomes more pronounced as the spatial size increases, indicating strong computational scalability. For the largest mesh evaluated (160 3 ), our method reduces memory usage and runtime in inference by 117× and 115×, respectively, compared with GNN-only baseline. More importantly, it shows higher accuracy and stronger spatiotemporal capability than the GNN-only baseline, especially in long-term testing. Such combination of scalability and accuracy is essential for simulating realistic material microstructures over extended time scales. The improvements can be attributed to the bijective autoencoder’s ability to compress information losslessly from spatial domain into a high dimensional feature space, thereby producing more expressive latent features for the GNN to learn from, while also contributing its own spatiotemporal modeling capability. Training data are generated from stochastic grain growth simulations, providing realistic variability for learning robust microstructure evolution. Comprehensive system validation confirms that the model is accurate, robust, and scalable.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Characterization of laser-accelerated proton beams from a 0.5 kJ sub-picosecond laser for radiography applications

Laser-accelerated ion beams show promise for many applications, including high-resolution flash imaging of static or dynamic objects in next-generation radiography to probe materials and plasmas in extreme environments and inertial confinement fusion. To scale up ion beam production for radiography applications, we conducted experiments using sub-picosecond lasers up to 0.5 kJ at the OMEGA-EP facility to characterize proton beams from solid targets, primarily CH/CD sub-micron thin films from which ion beams were also used for static and dynamic radiography for the first time. For standalone sub-micron thin CH films, the highest detected proton energy is in the range of 72–97 MeV. Proton beams with highest energy near or above 60 MeV at full laser energy and similar beam profiles are also measured from low-density CD foams and flat CH foil target of micrometer-scale thickness. The ~ 700–800 nm CH/CD foils achieve the highest ion yield among the targets tested. For sub-micron thin films, the laser prepulse can expand the target and lead to complex interactions, which is simulated using coupled hydrodynamic and two-step kinetic models. Simulations suggest the presence of a micrometer-scale preplasma plateau with near-critical density and further indicate that target normal sheath acceleration, electron heating from Relativistic transparency in the preplasma plateau, and background proton reflection from carbon ion front at the rear side contribute to the resulting proton spectrum from these sub-micron thin targets at various stages. These proton beams show strong potential for radiography and for production of secondary sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Real-Time Optical Monitoring and Simulations of Gas Phase Kinetics in InN Vapor Phase Epitaxy at High Pressure

Understanding the kinetics of nucleation and coalescence of heteroepitaxial thin films is a crucial step in controlling a chemical vapor deposition process, since it defines the perfection of the heteroepitaxial film both in terms of extended defect formation and chemical integrity of the interface. The initial nucleation process also defines the film quality during the later stages of film growth. The growth of emerging new materials heterostructures such as InN or In-rich Ga(x)In(1-x)N require deposition methods operating at higher vapor densities due to the high thermal decomposition pressure in these materials. High nitrogen pressure has been demonstrated to suppress thermal decomposition of InN, but has not been applied yet in chemical vapor deposition or etching experiments. Because of the difficulty with maintaining stochiometry at elevated temperature, current knowledge regarding thermodynamic data for InN, e.g., its melting point, temperature-dependent heat capacity, heat and entropy of formation are known with far less accuracy than for InP, InAs and InSb. Also, no information exists regarding the partial pressures of nitrogen and phosphorus along the liquidus surfaces of mixed-anion alloys of InN, of which the InN(x)P(1-x) system is the most interesting option. A miscibility gap is expected for InN(x)P(1-x) pseudobinary solidus compositions, but its extent is not established at this point by experimental studies under near equilibrium conditions. The extension of chemical vapor deposition to elevated pressure is also necessary for retaining stoichiometric single phase surface composition for materials that are characterized by large thermal decomposition pressures at optimum processing temperatures.

Dietz, Nikolaus↗