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 883 records · Page 49

GPU-accelerated Software Library for Unsteady Flamelet Modeling of Turbulent Combustion with Complex Chemical Kinetics

Modeling and simulation of turbulent combustion are computationally demanding and challenging problems due to the complexity of turbulent flows and the complex chemical reaction kinetics. Recent advances in high performance computing (HPC) has allowed simulations to utilize increasing computer power to obtain newer insights that were not feasible before. The arrival of GPUs and many-core processors is set to revolutionize the HPC landscape. This paper presents an unsteady flamelet software library developed for utilizing GPUs as an onboard accelerator while the flow solver continues to run on conventional host processors using legacy software. The software was architected to expose multiple levels of parallelism and concurrency available in the unsteady flamelet model which is also necessary to realize the full potential of these accelerators. While the model and its coupling with the flow solver have been known for some time, the potential to target these models and libraries for newer HPC architectures provides a greater incentive for pursuing combustion model development.

Sankaran, Ramanan↗

Local reduced-order modeling for electrostatic plasmas by physics-informed solution manifold decomposition

Despite advancements in high-performance computing and modern numerical algorithms, computational cost remains prohibitive for multi-query kinetic plasma simulations. Here, in this work, we develop data-driven reduced-order models (ROMs) for collisionless electrostatic plasma dynamics, based on the kinetic Vlasov-Poisson equation. Our ROM approach projects the equation onto a linear subspace defined by the proper orthogonal decomposition (POD) modes. We introduce an efficient tensorial method to update the nonlinear term using a precomputed third-order tensor. We capture multiscale behavior with a minimal number of POD modes by decomposing the solution manifold into multiple time windows and creating temporally local ROMs. We consider two strategies for decomposition: one based on the physical time and the other based on the electric field energy. Applied to the 1D1V Vlasov–Poisson simulations, that is, prescribed E-field, Landau damping, and two-stream instability, we demonstrate that our ROMs accurately capture the total energy of the system both for parametric and time extrapolation cases. The temporally local ROMs are more efficient and accurate than the single ROM. In addition, in the two-stream instability case, we show that the energy-windowing reduced-order model (EW-ROM) is more efficient and accurate than the time-windowing reduced-order model (TW-ROM). With the tensorial approach, EW-ROM solves the equation approximately 90 times faster than Eulerian simulations while maintaining a maximum relative error of 7.5% for the training data and 11% for the testing data.

Electrostatic plasmas↗

Verification of 5D continuum gyrokinetic code COGENT: Studies of kinetic drift wave instability

COGENT (Continuum Gyrokinetic Edge New Technology) is a kinetic plasma simulation code that is being developed by the edge simulation laboratory (ESL) collaboration. The original version of the code has been developed in a 4D phase space (2D configuration space and 2D velocity space) to address kinetic plasma phenomena in a complex magnetic field geometry including the core, magnetic separatrix, and scrape-off layer regions. This paper is focused on extending the original 4D phase space to a 5D phase space (3D2V) to address full kinetic turbulences in a tokamak edge region. Here, we report on the current status of 5D COGENT, which presently operates in a shear-less, simple slab geometry. As a verification study, we use the problems of collision-less drift wave instability (universal instability) including its modification in the presence of collisional effects. The simulation model includes the gyrokinetic equations for ion and electron species coupled to the long-wavelength limit of the gyro-Poisson equation. Collisional effects are represented by the Krook collision model. Linear analytical results for the drift mode growth rate and real frequency are recovered, and the non-linear stage is modelled and analysed as well. In addition, extensive 5D runs have been performed to address the effects of the drift wave instability on blob/filamentary structures characteristic of a tokamak edge. Finally, a helical-shaped potential perturbation is observed to grow exponentially in time while spinning around the filament axis with electron drift frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensitivity analysis of large system of chemical kinetic parameters for engine combustion simulation

In this study, the authors applied the state-of-the art sensitivity methods to downselect system parameters from 4000+ to 8, (23000+ -> 4000+ -> 84 -> 8). This analysis procedure paves the way for future works: (1) calibrate the system response using existed experimental observations, and (2) predict future experiment results, using the calibrated system.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Atomic-scale understanding of oxide growth and dissolution kinetics of Ni-Cr alloys

Aqueous corrosion of metals is governed by formation and dissolution of a passivating, multi-component surface oxide. Unfortunately, a detailed atomistic description is challenging due to the compositional complexity and the need to consider multiple kinetic factors simultaneously. To this end, we combine experiments with a first-principles-derived, multiscale computational framework that transcends thermodynamic descriptions to explicitly simulate the kinetic evolution of surface oxides of Ni-Cr alloys as a function of composition, temperature, pH, and applied voltage. In the absence of pitting, we identify three distinct voltage regimes, which are kinetically dominated by oxide growth, dissolution, and competitive dissolution and reprecipitation. Evolving compositional gradients and oxide thickness are revealed, including a transition between a metastable Ni-Cr mixed oxide and a thick, porous Ni-dominated oxide. Beyond elucidating the underlying physics, we highlight the need for competing kinetics in models to properly predict the transition from passivation to corrosion. Our results provide a key step towards co-design of alloy composition alongside environmental conditions for sustainable use across a variety of critical energy and infrastructure applications.

36 MATERIALS SCIENCE↗

VO2 kinetics during submaximal exercise following simulated weightlessness

A study is presented of the effects of deconditioning following 7 days of continuous head-down (-6 degrees) bedrest on changes in steady-state VO2, O2, and recovery VO2 during the performance of constant-load exercises. The deconditioning effects of bedrest on the physical working capacity were manifested in the subjects by significant changes in VO2 kinetics following the 7 days of head-down bedrest. While the subjects demonstrated the ability to attain similar steady-state VO2, simulated weightlessness using head-down bedrest resulted in a reduction of total VO2 capacity and an increase in the O2 deficit and VO2 halftime during submaximal constant-load exercise. It is concluded that this change in VO2 kinetics was induced by reexposure of the cardiovascular system to the +1 Gz (upright) environment.

Convertino, V. A.↗

Accelerated Steam Methane Reforming by Dynamically Applied Charges

Catalyst design has traditionally focused on tuning active site properties to optimally bind reaction intermediates and balance the kinetic requirements of multiple competing chemical processes, as necessitated by the Sabatier principle. It has recently been proposed that for reactions following certain potential energy landscapes, the activity limit imposed by the Sabatier principle may be overcome by using programmed oscillations of surface electron density at the timescales of surface reactions (i.e., “catalytic resonance”). Here, we use a combination of density functional theory (DFT) simulations and transient kinetic models (TKMs) to simulate the kinetics of steam methane reforming (SMR) on Ru(211) surfaces under statically and dynamically applied charges. DFT-calculated binding energies of SMR intermediates and transition states exhibit strong sensitivity to positively applied charges and follow unique scaling relationships that deviate from linear periodic trends across transition metals. Our simulations demonstrate that applying a small positive charge to Ru dramatically enhances the steady-state turnover frequency (TOF) of SMR by up to 5 orders of magnitude above the TOF observed over neutral Ru. Thus, statically charging Ru catalysts may be an effective strategy to lower the temperature requirements for SMR. Dynamic square-wave oscillations in charge resulted in SMR catalytic resonance with an onset frequency f ∼ 106 Hz and the corresponding average TOFs exceeding the statically charged Ru surface by an additional 15%. Here, based on sensitivity analyses performed for the two end points of oscillation, we propose that dynamic TOF improvement beyond the Sabatier maximum can be expected when the system oscillates between two kinetic regimes that are uniquely controlled by distinct elementary steps.

Catalysts↗

Kinetics of chemical ordering in a Ag-Pt nanoalloy particle via first-principles simulations

The energetics and kinetic energy barriers of vacancy/atom exchange in a 37-atom truncated octahedron Ag-Pt binary cluster in the Ag-rich range of compositions are investigated via a first-principles atomistic approach. The energy of the local minima obtained considering various distributions of a single vacancy and a few Pt atoms within the cluster and the energy barriers connecting them are evaluated using accurate density-functional calculations. The effects of the simultaneous presence of a vacancy and Pt atoms are found to be simply additive when their distances are larger than first-neighbors, whereas when they can be stabilizing at low Pt content due to the release of strain by the Pt/vacancy interaction or destabilizing close to a perfect Pt(core)/Ag(shell) arrangement. It is found that alloying with Pt appreciably increases the barriers for homotops transformations, thus rationalizing the issues encountered at the experimental level in producing Ag-Pt equilibrated nanoparticles and bulk phase diagram.

36 MATERIALS SCIENCE↗

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E×B drift to cause a non-diffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should in principle relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

intrinsic rotation↗

Modeling and verification of dynamic field ionization for laser-target interactions

Integrating field ionization models into kinetic plasma simulations is required for a variety of applications, especially when field strengths vary from low to high regimes, such as in laser-target interactions. The introduction of new physics models into kinetic codes requires a rigorous verification of their accuracy through well-defined verification problems. In this work, the field ionization model that has been included in the kinetic plasma code VPIC is presented, along with the detailed approach adopted for its integration. This model includes a comprehensive range of field ionization mechanisms: multiphoton ionization, tunneling ionization, and barrier suppression ionization. New verification problems employed to evaluate the ionization model's fidelity are outlined, and the simulation parameters that affect the accuracy of simulation results are explored. Additionally, this work addresses the impact of field ionization on computational performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mesoscale model for fission-induced recrystallization in U-7Mo alloy

A mesoscale model is developed by integrating the rate theory and phase-field models and is used to study the fission-induced recrystallization in U-7Mo alloy. The rate theory model is used to predict the dislocation density and the recrystallization nuclei density due to irradiation. The predicted fission rate and temperature dependences of the dislocation density are in good agreement with experimental measurements. This information is used as input for the multiphase phase-field model to investigate the fission-induced recrystallization kinetics. The simulated recrystallization volume fraction and bubble induced swelling agree well with experimental data. The effects of the fission rate, initial grain size, and grain morphology on the recrystallization kinetics are discussed based on an analysis of recrystallization growth rate using the modified Avrami equation. Here, we conclude that the initial microstructure of the U-Mo fuels, especially the grain size, can be used to effectively control the rate of fission-induced recrystallization and therefore swelling.

36 MATERIALS SCIENCE↗

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E × B drift to cause a nondiffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine f L , the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, f L scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should, in principle, relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear $\hat{s}$ was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of f L with the electron collisionality. Still, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theory Aided Design of Active and Durable Nanoscale Cathode Catalysts (Final Report)

One of the greatest technological hurdles to deployment of fuel cells relates to the sluggish activity, low durability and the high cost of the catalysts that are currently employed. For automotive PEM fuel cells to become commercially viable, the Pt-specific power density would need to be reduced to less than 0.2gPt/kW (at cell voltages >0.65 V). This would require the Pt loadings to be less than 0.15 mg Pt /cm 2 MEA within the membrane electrode assembly. This could be achieved by enhancing the catalytic activity at the cathode, thus lowering its overpotential. Various different Pt-alloys have shown 2-4 times enhanced activities over Pt alone but still suffer some of the same durability issues as that of the pure Pt. There is a general loss of active Pt due to dissolution and sintering. While there have been a number of elegant fundamental experimental and theoretical studies on ideal single crystal Pt and Pt alloy surfaces which have helped to elucidate the factors that control the activity, there have been very few fundamental studies focused on the stability, reactivity and durability of well-defined Pt nanoparticles. We carried out ab initio density functional theory together with a novel double reference method that we developed to simulate constant potential electrochemical systems in order to model the electrocatalytic reduction of oxygen over model Pt alloy surfaces and nanoparticles. These simulations were used to probe the factors that control the electrocatalytic activity, guide the potential selection of new materials and test their stability under reaction conditions. Ab initio calculations were used to determine the reaction energies and activation barriers for a comprehensive array of different elementary adsorption, desorption, surface reaction and diffusion steps over Pt and Pt alloys involved in the electroreduction of oxygen as a function of surface coverage, the electrochemical potential and temperature. The calculations were used to simulate the potential dependent adsorption and surface reaction energies along with activation barriers in order to determine the kinetics for different surface structures and structural features (step edge and corner sites) to provide the necessary input for kinetic Monte Carlo simulations to follow the rates of reaction. More coarse-grained simulated annealing methods were used to help establish the lowest energy structures and morphologies for different Pt and Pt-alloy nanoparticles that form. The results from the DFT calculations were used to establish an ab initio-derived kinetic database that was used in both 2D and 3D kinetic Monte Carlo simulations that were used to follow the electrocatalytic performance over different particle sizes, shapes and compositions. The ab initio calculations together with the kinetic Monte Carlo simulations were used to complete the following objectives: 1) Determine the controlling elementary reaction pathways and intrinsic kinetics involved in ORR and their potential dependent behavior; 2) Establish the influence of the extrinsic reaction environment including surface structure, alloy composition and spatial arrangement and the humidity on ORR kinetics; 3) Elucidate the effects of particle size and morphology as well as the atomic structure and composition of nanoparticles of Pt, Pt 3 Co, Pt 3 Ni, Pt 3 Fe and other Pt-alloys on the electrocatalytic activity and reaction selectivity; 5) Construct surface Pourbaix phase diagrams for different Pt and Pt-alloys in order to map out the most stable surface phases for different material compositions and reaction conditions; 6) Elucidate the mechanisms that control dissolution and re-deposition of Pt.

08 HYDROGEN↗

On the evolution of the steady state in gold-silica nanocomposites under sustained irradiation

We investigate the kinetic evolution of a nanocomposite submitted to sustained irradiation. The study focuses on gold nanoparticles embedded within a silica host matrix. Irradiations with 4 MeV Au {sup ++} ions are performed in the temperature range from 300 K to 1100 K and for fluences up to 8 × 10{sup 16 }cm{sup −2}. The experimental results are combined with kinetic Monte Carlo simulations to fully characterize the kinetic evolution of the nanoparticles under irradiation. In particular, the change of the steady state with the temperature is investigated in details. First, we demonstrate the existence of a direct-coarsening regime at high temperature, e.g., above 900 K, and of a coarsening-resistant regime at low temperature, e.g., below 700 K. Second, we show that the transition between these two regimes takes place in a temperature region, where the nucleation of new nanoparticles is hindered, and the mass is redistributed among the existing nanoparticles. We also use our experimental results to critically analyze the existing theoretical models. As it has recently been pointed out that the physical principles underlying precipitation under irradiation are similar to those that govern the behavior of driven alloys, this work can be considered as an attempt to bridge the gap between two fields of experimental and theoretical studies, viz., those on driven alloys and those on metal-glass nanocomposites submitted to a sustained irradiation.

36 MATERIALS SCIENCE↗

Lower Hybrid ←→Whistler Wave Conversion [Slides]

In this presentation we discuss recent work concerning the conversion of whistler waves to lower hybrid waves (as well as the inverse process). These efforts have been motivated by the issue of attenuation of upward propagating whistler waves in the ionosphere generated by VLF transmitters on the ground, i.e., the 'Starks 20 db' problem, which affects the lifetimes of energetic electrons trapped in the geomagnetic field at low magnetic altitude (L). We discuss recent fluid and kinetic plasma simulations as well as ongoing experiments at UCLA to quantify linear and nonlinear mode conversion of lower hybrid to whistler waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of detonation structures with hydrocarbon fuels for application towards rotating detonation engines

Rotating Detonation Engines (RDEs) provide a promising approach to increasing efficiency of gas turbine combustors by utilizing detonation-driven combustion process. While RDEs have been studied extensively in the past, much of this work has focused on the use of hydrogen as fuel. In order to develop RDEs for power generation applications, it is necessary to understand the physics of hydrocarbon detonation. This study utilizes detailed chemical kinetics to simulate a sequence of hydrocarbon-based canonical RDE configurations. The cases emulate ethylene and methane detonation in air with varying degrees of hydrogen dilution as well as a range of operating conditions. The results indicate the while ethylene-based detonations are not significantly affected by the addition of hydrogen, methane mixtures exhibit large changes to the detonation structure. In particular, the critical pressure at which heat release reaches a peak changes to lower values as the hydrogen concentration in the reactant mixture is increased. Detailed comparisons with one-dimensional profiles and the impact of back pressure on the detonation and post-detonation flow are analyzed. Further, profiles of species extracted from the simulations are compared with one-dimensional detonation profiles. Comparisons with theoretical models for thrust and specific impulse are also provided. It is established that the use of detailed chemical kinetics provides a reliable approach to assessing the performance characteristics of RDEs.

Sato, Takuma↗

The reduction of magnetic reconnection outflow jets to sub-Alfvénic speeds

The outflow velocity of jets produced by collisionless magnetic reconnection is shown to be reduced by the ion exhaust temperature in fully kinetic particle in cell simulations and in situ satellite observations. We derive a scaling relationship for the outflow velocity based on the upstream Alfvén speed and the parallel ion exhaust temperature, which is verified in kinetic simulations and observations. The outflow speed reduction is shown to be due to the firehose instability criterion, and so, for large enough guide fields, this effect is suppressed and the outflow speed reaches the upstream Alfvén speed based on the reconnecting component of the magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗