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 865 records · Page 48

Directing Anisotropic Assembly of Metallic Nanoclusters by Exploiting Linear Trio Interactions and Quantum Size Effects: Au Chains on Ag(100) Thin Films

Discovery and understanding of mechanisms for kinetically controlled growth of metal nanoclusters can be enabled by realistic atomistic-level modeling with ab initio kinetics. KMC simulation of such a model for Au deposition on Ag(100) films reveals the formation of single-atom-wide Au chains below 275 K, even though 2D islands are thermodynamically preferred. Chain formation is shown to reflect a combination of strong linear trio attractions guiding assembly and a weak driving force and slow rate of transformation of 1D chains to 2D islands (or sometimes irreversible rounding of adatoms from chain sides to ends). Behavior can also be tuned by quantum size effects: chain formation predominates on 2-monolayer Ag(100) films supported on NiAl(100) at 250 K for low coverages but not on 1- or 3-monolayer films, and longer chains form than on bulk Ag(100). Finally, our predictive kinetic modeling shows the potential for simulation-guided discovery and analysis of novel self-assembly processes.

ab initio periphery diffusion kinetics↗

Ephemeral collision complexes mediate chemically termolecular transformations that affect system chemistry [Ephemeral collision complexes induce chemically termolecular transformations that affect global chemistry]

Termolecular association reactions involve ephemeral collision complexes—formed from the collision of two molecules—that collide with a third and chemically inert ‘bath gas’ molecule that simply transfers energy to/from the complex. These collision complexes are generally not thought to react chemically on collision with a third molecule in the gas-phase systems of combustion and planetary atmospheres. Such ‘chemically termolecular’ reactions, in which all three molecules are involved in bond making and/or breaking, were hypothesized long ago in studies establishing radical chain branching mechanisms, but were later concluded to be unimportant. Here, with data from ab initio master equation and kinetic-transport simulations, we reveal that reactions of H+O 2 collision complexes with other radicals constitute major kinetic pathways under common combustion situations. These reactions are also found to influence flame propagation speeds, a common measure of global reactivity. As a result, analogous chemically termolecular reactions mediated by ephemeral collision complexes are probably of significance in various combustion and planetary environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electron Injection via Modified Diffusive Shock Acceleration in High-Mach-number Collisionless Shocks

The ability of collisionless shocks to efficiently accelerate nonthermal electrons via diffusive shock acceleration (DSA) is thought to require an injection mechanism capable of preaccelerating electrons to high enough energy where they can start crossing the shock front potential. We propose, and show via fully kinetic plasma simulations, that in high-Mach-number shocks electrons can be effectively injected by scattering in kinetic-scale magnetic turbulence produced near the shock transition by the ion Weibel, or current filamentation, instability. We describe this process as a modified DSA mechanism where initially thermal electrons experience the flow velocity gradient in the shock transition and are accelerated via a first-order Fermi process as they scatter back and forth. The electron energization rate, diffusion coefficient, and acceleration time obtained in the model are consistent with particle-in-cell simulations and with the results of recent laboratory experiments where nonthermal electron acceleration was observed. This injection model represents a natural extension of DSA and could account for electron injection in high-Mach-number astrophysical shocks, such as those associated with young supernova remnants and accretion shocks in galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of point kinetics, improved quasistatic and theta method as space-time kinetics solvers in DONJON-3 simulations

To ensure the safety of nuclear reactors, we have to simulate accurately their normal operation and also accident cases. To perform transient calculations, coupled neutronic and fhermo-hydraulic codes are used. This article compares three neutronic solvers. The first one is the point kinetic approach where the flux shape is constant during all the transient. For the second method (the improved quasistatic method), the flux shape is constant but only during small time steps. Finally, we used the theta approach where both flux and precursors distributions vary with time and space. Transients of Lost Of Coolant Accident in CANDU-6 reactors have been simulated with DONJON and the outputs of a thermalhydraulic system code. Results show that the point kinetics is inappropriate for transient with large distortion of the flux shape. Improved quasistatic and theta methods give relatively similar results. However, the improved quasistatic approach is less stable and a little bit more sensitive on time-step and spatial discretization than the theta method is. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Review of Molten Salt Reactor Kinetics Models

Interest in circulating fuel reactors (CFRs), particularly molten salt reactors (MSRs) of the fluid fuel type, has been growing in the last two decades. Starting with a resurgence of interest in Europe, there have been a growing number of methods proposed and codes developed to model the kinetics of CFRs, which is a capability essential to the design and evaluation of such reactors. This work first reviews the physical phenomena unique to CFRs in light of current research and how CFR kinetics are impacted by these considerations. In general, it is found that the movement of delayed neutron precursors (DNPs) through the primary loop has significant impacts on transients at low reactor powers or those with significant spatial components such as a change in the primary loop mass flow rate. Effects on the neutron flux are exceedingly minimal and entirely negligible. An extensive review of published models and methods for simulating CFR kinetics is presented, along with transient simulations in fast and thermal neutron flux systems using representative codes from each of the main modeling categories. Comparisons among methods are presented as are recommendations for their use or nonuse in various transient and work-flow scenarios. In general, it is recommended that time-resolved, multigroup neutron diffusion approaches be used to establish ranges of applicability for point reactor kinetics (PRK)–based approaches that themselves may not be applicable for all modeling situations. In such cases, it is suggested that quasi-static approaches be used where PRK-based approaches cannot be used. Finally, a review of common assumptions used in these models is presented, along with an evaluation of their impact on model performance. It is found that neglecting turbulent diffusion in open core–type CFRs is a poor assumption that leads to an underestimation of the reduction of the delayed neutron fraction. Furthermore, it is seen that exclusion of secondary heat transfer loops in models leads to underestimation of transient peaks and troughs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A low-frequency variational model for energetic particle effects in the pressure-coupling scheme

Energetic particle effects in magnetic confinement fusion devices are commonly studied by hybrid kinetic-fluid simulation codes whose underlying continuum evolution equations often lack the correct energy balance. While two different kinetic-fluid coupling options are available (current coupling and pressure coupling), this paper applies the Euler–Poincaré variational approach to formulate a new conservative hybrid model in the pressure-coupling scheme. In our case the kinetics of the energetic particles are described by guiding center theory. The interplay between the Lagrangian fluid paths and phase space particle trajectories reflects an intricate variational structure which can be approached by letting the four-dimensional guiding center trajectories evolve in the full six-dimensional phase space. Then, the redundant perpendicular velocity is integrated out to recover a four-dimensional description. A second equivalent variational approach is also reported, which involves the use of phase space Lagrangians. Not only do these variational structures confer on the new model a correct energy balance, but also they produce a cross-helicity invariant which is lost in the other pressure-coupling schemes reported in the literature.

Physics↗

Sheath-induced distortions in particle distributions near enhanced polar outflow probe particle sensors

We discuss sheath and kinetic effects on ion and electron distribution functions at the aperture of enhanced Polar Outflow Probe particle sensors. For this purpose, the interaction between the CASSIOPE spacecraft and space environment is simulated fully kinetically using the electrostatic Particle In Cell code PTetra. The simulations account for the geometry of the main features of the spacecraft body, the booms, and the sensors. In addition to the background plasma, the model also accounts for Earth magnetic field. The plasma parameters assumed in the simulations are obtained from the latest version of the International Reference Ionosphere (IRI) model and the value of magnetic field is obtained from the International Geophysical Reference Field model. Our analysis shows significant distortions in the ion distribution function in the plane of the sensor aperture, as well as in the direction along the boom holding the sensor. We argue that significant distortions and asymmetries should also occur at the aperture of the suprathermal electron imager when suprathermal electrons are detected, with energies of 5 eV or more.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of Kinetic Mechanisms for Next-Generation Fuels and CFD Simulation of Advanced Combustion Engines

Predictive chemical kinetic models are needed to represent next-generation fuel components and their mixtures with conventional gasoline and diesel fuels. These kinetic models will allow the prediction of the effect of alternative fuel blends in CFD simulations of advanced spark-ignition and compression-ignition engines. Enabled by kinetic models, CFD simulations can be used to optimize fuel formulations for advanced combustion engines so that maximum engine efficiency, fossil fuel displacement goals, and low pollutant emission goals can be achieved.

30 DIRECT ENERGY CONVERSION↗

Pedestal Fueling Simulations with a Coupled Kinetic-kinetic Plasma-neutral Transport Code

A Monte Carlo neutral transport routine, based on DEGAS2, has been coupled to the guiding center ion-electron-neutral neoclassical PIC code XGC0 to provide a realistic treatment of neutral atoms and molecules in the tokamak edge plasma. The DEGAS2 routine allows detailed atomic physics and plasma-material interaction processes to be incorporated into these simulations. The spatial pro le of the neutral particle source used in the DEGAS2 routine is determined from the uxes of XGC0 ions to the material surfaces. The kinetic-kinetic plasma-neutral transport capability is demonstrated with example pedestal fueling simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic (particle-in-cell) simulation of nonlinear laser absorption in a finite-size plasma with a background inhomogeneous magnetic field

In this paper, the effect of an external inhomogeneous magnetic field on the high intensity laser absorption rate in a sub-critical plasma has been investigated by employing a relativistic electromagnetic 1.5 dimensional particle-in-cell code. Relying on the effective nonlinear phenomena such as phase-mixing and scattering, this study shows that in a finite-size plasma the laser absorption increases with inhomogeneity of the magnetic field (i.e., reduction of characteristic length of inhomogeneous magnetic field, λ{sub p}) before exiting a considerable amount of laser energy from the plasma due to scattering process. On the other hand, the presence of the external inhomogeneous magnetic field causes the maximum absorption of laser to occur at a shorter time. Moreover, study of the kinetic results associated with the distribution function of plasma particles shows that, in a special range of the plasma density and the characteristic length of inhomogeneous magnetic field, a considerable amount of laser energy is transferred to the particles producing a population of electrons with kinetic energy along the laser direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrated Multiscale Modeling of Molecular Computing Devices. Final Report

The general theme of this research has been to expand the capabilities of a simulation technique, Kinetic Monte Carlo (KMC) and apply it to study self-assembled nano-structures on epitaxial thin films. KMC simulates thin film growth and evolution by replacing the detailed dynamics of the system's evolution, which might otherwise be studied using molecular dynamics, with an appropriate stochastic process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry

For the research of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Iodide oxidation by ozone at the surface of aqueous microdroplets

The kinetics of iodide oxidation by ozone are studied in single microdroplets to explore chemical reactivity at the air–water interface. Molecular simulations and kinetic modeling are employed to understand the underlying multiphase mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigation of thermolytic hydrogen generation rate of tank farm simulated and actual waste

To support resolution of Potential Inadequacies in the Safety Analysis for the Savannah River Site (SRS) Tank Farm, Savannah River National Laboratory conducted research to determine the thermolytic hydrogen generation rate (HGR) with simulated and actual waste. Gas chromatography methods were developed and used with air-purged flow systems to quantify hydrogen generation from heated simulated and actual waste at rates applicable to the Tank Farm Documented Safety Analysis (DSA). Initial simulant tests with a simple salt solution plus sodium glycolate demonstrated the behavior of the test apparatus by replicating known HGR kinetics. Additional simulant tests with the simple salt solution excluding organics apart from contaminants provided measurement of the detection and quantification limits for the apparatus with respect to hydrogen generation. Testing included a measurement of HGR on actual SRS tank waste from Tank 38. A final series of measurements examined HGR for a simulant with the most common SRS Tank Farm organics at temperatures up to 140 °C. The following conclusions result from this testing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Investigation of thermolytic hydrogen generation rate of tank farm simulated and actual waste

To support resolution of Potential Inadequacies in the Safety Analysis for the Savannah River Site (SRS) Tank Farm, Savannah River National Laboratory conducted research to determine the thermolytic hydrogen generation rate (HGR) with simulated and actual waste. Gas chromatography methods were developed and used with air-purged flow systems to quantify hydrogen generation from heated simulated and actual waste at rates applicable to the Tank Farm Documented Safety Analysis (DSA). Initial simulant tests with a simple salt solution plus sodium glycolate demonstrated the behavior of the test apparatus by replicating known HGR kinetics. Additional simulant tests with the simple salt solution excluding organics apart from contaminants provided measurement of the detection and quantification limits for the apparatus with respect to hydrogen generation. Testing included a measurement of HGR on actual SRS tank waste from Tank 38. A final series of measurements examined HGR for a simulant with the most common SRS Tank Farm organics at temperatures up to 140 °C. The following conclusions result from this testing. HGR was measured for an actual-waste Tank 38 sample. The HGR measurements at 80, 90, 100, and 110 °C were 1.7×10 -7 , 4.7×10 -7 , 5.9×10 -7 , and 2.6×10 -6 ft 3 hr -1 gal -1 , respectively. The Limit of Detection and Limit of Quantification for the HGR measurements were 2.6×10-8 ft3 hr-1 gal-1 and 5.1×10 -8 ft 3 hr -1 gal -1 , respectively. HGR was measured in a high boiling point (HBP) simulant with major known SRS Tank Farm organic compounds. Testing used fresh organics and did not simulate the long-term exposure of these organics already experienced in the waste. Hence, testing is expected to provide high bias to the hydrogen generation rates. The final HGR measurements at 75, 100, and 140 °C were 6.2×10 -7 , 3.0×10 -6 , and 8.1×10 -5 ft 3 hr -1 gal -1 , respectively. Decreases in the HGR were noted at isothermal conditions with steady gas tracer concentration, which clouds development of a temperature dependent correlation for generation rate. While the simulant contains high concentrations of major SRS organics tested at high temperatures found in the 3H Evaporator, this limited data set may not fully bound values under all possible operating conditions. The HGR measurement apparatus was successfully employed to replicate hydrogen production rate experiments. Simulant without deliberately added organic carbon was found to contain approximately 200 mg/L of organic carbon from tramp organics in the reagents used to prepare the solution. The tests using such simulant generated hydrogen at rates consistent with, but somewhat lower than, those predicted at equivalent carbon and aluminum loadings by methods used at the Hanford Site. Prior to identifying the organic contaminants, potential sources of hydrogen in simulants without deliberately added organic carbon were examined. Testing determined that this hydrogen is not due to caustic reaction of organic materials of construction within the reaction vessel, electrolysis due to current leak, anode/cathode reactions, corrosion, or the method of heating employed in the experiment. During the HBP and Tank 38 testing, a peak was identified at a retention time greater than the retention time of the inert tracer of krypton. This peak has been identified as methane, and the data have been reexamined to obtain methane concentrations. For the HBP testing, methane concentration was constant at the highest temperature studied (140 °C) and higher than hydrogen. At the lower temperatures (100 and 75 °C), methane and hydrogen tracked closely, with methane concentration generally lower than hydrogen. For Tank 38, methane concentration was approximately 1/3 of the hydrogen concentration during testing (100 and 110 °C). At lower temperatures for Tank 38 testing (90, 80, and 27 °C), methane fell below the quantification range.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Dissipation measures in weakly collisional plasmas

ABSTRACT The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell vpic, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov–Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic and Potential Sputtering of Lunar Regolith: Contribution of Solar-Wind Heavy Ions

Sputtering of lunar regolith by protons as well as solar-wind heavy ions is considered. From preliminary measurements of H+, Ar+1, Ar+6 and Ar+9 ion sputtering of JSC-1A AGGL lunar regolith simulant at solar wind velocities, and TRIM simulations of kinetic sputtering yields, the relative contributions of kinetic and potential sputtering contributions are estimated. An 80-fold enhancement of oxygen sputtering by Ar+ over same-velocity H+, and an additional x2 increase for Ar+9 over same-velocity Ar+ was measured. This enhancement persisted to the maximum fluences investigated is approximately 1016/cm (exp2). Modeling studies including the enhanced oxygen ejection by potential sputtering due to the minority heavy ion multicharged ion solar wind component, and the kinetic sputtering contribution of all solar wind constituents, as determined from TRIM sputtering simulations, indicate an overall 35% reduction of near-surface oxygen abundance. XPS analyses of simulant samples exposed to singly and multicharged Ar ions show the characteristic signature of reduced (metallic) Fe, consistent with the preferential ejection of oxygen atoms that can occur in potential sputtering of some metal oxides.

Meyer, F. W.↗

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗