The next frontier in predictive simulation: implicit kinetic simulation at the system scale [Slides]
Abstract not provided.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
SEARCH · Search NASA
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.
Fully kinetic simulations of megajoule-scale dense plasma focus
Dense plasma focus (DPF) Z-pinch devices are sources of copious high energy electrons and ions, x-rays, and neutrons. Megajoule-scale DPFs can generate 10{sup 12} neutrons per pulse in deuterium gas through a combination of thermonuclear and beam-target fusion. However, the details of the neutron production are not fully understood and past optimization efforts of these devices have been largely empirical. Previously, we reported on the first fully kinetic simulations of a kilojoule-scale DPF and demonstrated that both kinetic ions and kinetic electrons are needed to reproduce experimentally observed features, such as charged-particle beam formation and anomalous resistivity. Here, we present the first fully kinetic simulation of a MegaJoule DPF, with predicted ion and neutron spectra, neutron anisotropy, neutron spot size, and time history of neutron production. The total yield predicted by the simulation is in agreement with measured values, validating the kinetic model in a second energy regime.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Kinetic simulations and reduced modeling of longitudinal sideband instabilities in non-linear electron plasma waves
Kinetic Vlasov simulations of one-dimensional finite amplitude Electron Plasma Waves are performed in a multi-wavelength long system. A systematic study of the most unstable linear sideband mode, in particular its growth rate γ and quasi- wavenumber δk, is carried out by scanning the amplitude and wavenumber of the initial wave. Simulation results are successfully compared against numerical and analytical solutions to the reduced model by Kruer et al. [Phys. Rev. Lett. 23, 838 (1969)] for the Trapped Particle Instability (TPI). A model recently suggested by Dodin et al. [Phys. Rev. Lett. 110, 215006 (2013)], which in addition to the TPI accounts for the so-called Negative Mass Instability because of a more detailed representation of the trapped particle dynamics, is also studied and compared with simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Kinetic simulation of neutral/ionized gas and electrically charged dust in the coma of comet 67P/Churyumov-Gerasimenko
The cometary coma is a unique phenomenon in the solar system being a planetary atmosphere influenced by little or no gravity. As a comet approaches the sun, the water vapor with some fraction of other gases sublimate, generating a cloud of gas, ice and other refractory materials (rocky and organic dust) ejected from the surface of the nucleus. Sublimating gas molecules undergo frequent collisions and photochemical processes in the near-nucleus region. Owing to its negligible gravity, comets produce a large and highly variable extensive dusty coma with a size much larger than the characteristic size of the cometary nucleus.The Rosetta spacecraft is en route to comet 67P/Churyumov-Gerasimenko for a rendezvous, landing, and extensive orbital phase beginning in 2014. Both, interpretation of measurements and safety consideration of the spacecraft require modeling of the comet's dusty gas environment.In this work we present results of a numerical study of multispecies gaseous and electrically charged dust environment of comet Chyuryumov-Gerasimenko. Both, gas and dust phases of the coma are simulated kinetically. Photolytic reactions are taken into account. Parameters of the ambient plasma as well as the distribution of electric/magnetic fields are obtained from an MHD simulation of the coma connected to the solar wind. Trajectories of ions and electrically charged dust grains are simulated by accounting for the Lorentz force and the nucleus gravity.
79 ASTRONOMY AND ASTROPHYSICS↗
THERMAL ESCAPE IN THE HYDRODYNAMIC REGIME: RECONSIDERATION OF PARKER's ISENTROPIC THEORY BASED ON RESULTS OF KINETIC SIMULATIONS
The one-dimensional steady-state problem of thermal escape from a single-component atmosphere of mon- and diatomic gases is studied in the hydrodynamic (blow-off) regime using the direct simulation Monte Carlo method for an evaporative-type condition at the lower boundary. The simulations are performed for various depths into an atmosphere, indicated by a Knudsen number, Kn{sub 0}, equal to the ratio of the mean free path of molecules to the radial position of the source surface, ranging from 10 to 10{sup -5}, and for the range of the source Jeans parameter, {lambda}{sub 0}, equal to the ratio of gravitational and thermal energies, specific to blow-off. The results of kinetic simulations are compared with the isentropic model (IM) and the Navier-Stokes model. It is shown that the IM can be simplified if formulated in terms of the local Mach number and Jeans parameter. The simulations predict that at Kn{sub 0} < {approx} 10{sup -3} the flow includes a near-surface non-equilibrium Knudsen layer, a zone where the flow can be well approximated by the IM, and a rarefied far field. The corresponding IM solutions, however, only approach Parker's critical solution as {lambda}{sub 0} approaches the upper limit for blow-off. The IM alone is not capable for predicting the flow and requires boundary conditions at the top of the Knudsen layer. For small Kn{sub 0}, the scaled escape rate and energy loss rate are found to be independent of {lambda}{sub 0}. The simulation results can be scaled to any single-component atmosphere exhibiting blow-off if the external heating above the lower boundary is negligible, in particular, to sublimation-driven atmospheres of Kuiper belt objects.
79 ASTRONOMY AND ASTROPHYSICS↗
Test of Shi et al. Method to Infer the Magnetic Reconnection Geometry from Spacecraft Data: MHD Simulation with Guide Field and Antiparallel Kinetic Simulation
When analyzing data from an array of spacecraft (such as Cluster or MMS) crossing a site of magnetic reconnection, it is desirable to be able to accurately determine the orientation of the reconnection site. If the reconnection is quasi-two dimensional, there are three key directions, the direction of maximum inhomogeneity (the direction across the reconnection site), the direction of the reconnecting component of the magnetic field, and the direction of rough invariance (the "out of plane" direction). Using simulated spacecraft observations of magnetic reconnection in the geomagnetic tail, we extend our previous tests of the direction-finding method developed by Shi et al. (2005) and the method to determine the structure velocity relative to the spacecraft Vstr. These methods require data from four proximate spacecraft. We add artificial noise and calibration errors to the simulation fields, and then use the perturbed gradient of the magnetic field B and perturbed time derivative dB/dt, as described by Denton et al. (2010). Three new simulations are examined: a weakly three-dimensional, i.e., quasi-two-dimensional, MHD simulation without a guide field, a quasi-two-dimensional MHD simulation with a guide field, and a two-dimensional full dynamics kinetic simulation with inherent noise so that the apparent minimum gradient was not exactly zero, even without added artificial errors. We also examined variations of the spacecraft trajectory for the kinetic simulation. The accuracy of the directions found varied depending on the simulation and spacecraft trajectory, but all the directions could be found within about 10 for all cases. Various aspects of the method were examined, including how to choose averaging intervals and the best intervals for determining the directions and velocity. For the kinetic simulation, we also investigated in detail how the errors in the inferred gradient directions from the unmodified Shi et al. method (using the unperturbed gradient) depended on the amplitude of the calibration errors. For an accuracy of 3 for the maximum gradient direction, the calibration errors could be as large as 3% of reconnection magnetic field, while for the same accuracy for the minimum gradient direction, the calibration errors could only be as large as 0.03% of the reconnection magnetic field. These results suggest that the maximum gradient direction can normally be determined by the unmodified Shi et al. method, while the modified method or some other method must be used to accurately determine the minimum gradient direction. The structure velocity was found with magnitude accurate to 2% and direction accurate to within 5%.
Denton, R.↗
Energy dynamics and current sheet structure in fluid and kinetic simulations of decaying magnetohydrodynamic turbulence
Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfvén waves, which interact and rapidly generate strong MHD turbulence. The total energy is conserved and the rate of turbulent energy decay is very similar in both codes, although the fluid code has numerical dissipation, whereas the kinetic code has kinetic dissipation. The inertial range power spectrum index is similar in both the codes. The fluid code shows a perpendicular wavenumber spectral slope of k{sub ⊥}{sup −1.3}. The kinetic code shows a spectral slope of k{sub ⊥}{sup −1.5} for smaller simulation domain, and k{sub ⊥}{sup −1.3} for larger domain. We estimate that collisionless damping mechanisms in the kinetic code can account for the dissipation of the observed nonlinear energy cascade. Current sheets are geometrically characterized. Their lengths and widths are in good agreement between the two codes. The length scales linearly with the driving scale of the turbulence. In the fluid code, their thickness is determined by the grid resolution as there is no explicit diffusivity. In the kinetic code, their thickness is very close to the skin-depth, irrespective of the grid resolution. This work shows that kinetic codes can reproduce the MHD inertial range dynamics at large scales, while at the same time capturing important kinetic physics at small scales.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Diffusion in confinement: kinetic simulations of self- and collective diffusion behavior of adsorbed gases
Explore the source record for details and available documents.
membrane, carbon capture, materials and chemistry ↗
The linear tearing instability in three dimensional, toroidal gyro-kinetic simulations
Linear gyro-kinetic simulations of the classical tearing mode in three-dimensional toroidal geometry were performed using the global gyro-kinetic turbulence code, GKW. The results were benchmarked against a cylindrical ideal MHD and analytical theory calculations. The stability, growth rate, and frequency of the mode were investigated by varying the current profile, collisionality, and the pressure gradients. Both collisionless and semi-collisional tearing modes were found with a smooth transition between the two. A residual, finite, rotation frequency of the mode even in the absence of a pressure gradient is observed, which is attributed to toroidal finite Larmor-radius effects. When a pressure gradient is present at low collisionality, the mode rotates at the expected electron diamagnetic frequency. However, the island rotation reverses direction at high collisionality. The growth rate is found to follow a η{sup 1∕7} scaling with collisional resistivity in the semi-collisional regime, closely following the semi-collisional scaling found by Fitzpatrick. The stability of the mode closely follows the stability analysis as performed by Hastie et al. using the same current and safety factor profiles but for cylindrical geometry, however, here a modification due to toroidal coupling and pressure effects is seen.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Kinetic Simulations of Microstructural Evolution (KSOME)
"Kinetic simulations of microstructural evolution" (KSOME) is an object Kinetic Monte Carlo (OMKC) simulation code. It is developed to simulate microstructural evolution in materials during irradiation. KSOME is developed from ground up with priority given to the flexibility, ease of upgradability, and computation efficiency, in that order. With previously developed simulation codes, various processes (more precisely, the diffusion-reaction processes) that can occur in a system were hardwired into the code. Hence, to perform higher fidelity simulations, a significant coding effort was required every time to incorporate either a more accurate description of existing processes or new processes. With KSOME, the necessity for frequent code upgrades whenever there is a need to perform higher fidelity simulations is either eliminated or reduced significantly. In KSOME, the categories of various diffusion-reaction processes and the data-management system are hardwired, while the execution of the diffusion-reaction processes specific to a system is specified via input files. In the present version of KSOME, the analytical expression required to calculate defect interaction radii are inputted via a text file. However, the ability to parse mathematical expressions from a text file can be extended to input the analytical expression to calculate other defect properties; including the expressions to calculate the change in the migration barriers of mobile defects due to long-range interaction with extended defects
Nandipati, Giridhar [Pacific Northwest National La↗
Jovian Plasmas Torus Interaction with Europa. Plasma Wake Structure and Effect of Inductive Magnetic Field: 3D Hybrid Kinetic Simulation
The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa moon-magnetosphere system with respect to a variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo Orbiter mission, and for planning flyby and orbital measurements (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy et al., 2007; Shematovich et al., 2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyroradius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions). Photoionization, electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider the models with Oþ þ and Sþ þ background plasma, and various betas for background ions and electrons, and pickup electrons. The majority of O2 atmosphere is thermal with an extended non-thermal population (Cassidy et al., 2007). In this paper, we discuss two tasks: (1) the plasma wake structure dependence on the parameters of the upstream plasma and Europa's atmosphere (model I, cases (a) and (b) with a homogeneous Jovian magnetosphere field, an inductive magnetic dipole and high oceanic shell conductivity); and (2) estimation of the possible effect of an induced magnetic field arising from oceanic shell conductivity. This effect was estimated based on the difference between the observed and modeled magnetic fields (model II, case (c) with an inhomogeneous Jovian magnetosphere field, an inductive magnetic dipole and low oceanic shell conductivity).
Jovian magnetosphere↗
Modeling Solar Wind Flow with the Multi-Scale Fluid-Kinetic Simulation Suite
Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) is a package of numerical codes capable of performing adaptive mesh refinement simulations of complex plasma flows in the presence of discontinuities and charge exchange between ions and neutral atoms. The flow of the ionized component is described with the ideal MHD equations, while the transport of atoms is governed either by the Boltzmann equation or multiple Euler gas dynamics equations. We have enhanced the code with additional physical treatments for the transport of turbulence and acceleration of pickup ions in the interplanetary space and at the termination shock. In this article, we present the results of our numerical simulation of the solar wind (SW) interaction with the local interstellar medium (LISM) in different time-dependent and stationary formulations. Numerical results are compared with the Ulysses, Voyager, and OMNI observations. Finally, the SW boundary conditions are derived from in-situ spacecraft measurements and remote observations.
79 ASTRONOMY AND ASTROPHYSICS↗
Hybrid Parallelization of Adaptive MHD-Kinetic Module in Multi-Scale Fluid-Kinetic Simulation Suite
The Multi-Scale Fluid-Kinetic Simulation Suite has a computational tool set for solving partially ionized flows. In this paper we focus on recent developments of the kinetic module which solves the Boltzmann equation using the Monte-Carlo method. The module has been recently redesigned to utilize intra-node hybrid parallelization. We describe in detail the redesign process, implementation issues, and modifications made to the code. Finally, we conduct a performance analysis.
97 MATHEMATICS AND COMPUTING↗
Comparisons of dense-plasma-focus kinetic simulations with experimental measurements
Dense-plasma-focus (DPF) Z-pinch devices are sources of copious high-energy electrons and ions, x rays, and neutrons. The mechanisms through which these physically simple devices generate such high-energy beams in a relatively short distance are not fully understood and past optimization efforts of these devices have been largely empirical. Previously we reported on fully kinetic simulations of a DPF and compared them with hybrid and fluid simulations of the same device. Here we present detailed comparisons between fully kinetic simulations and experimental data on a 1.2 kJ DPF with two electrode geometries, including neutron yield and ion beam energy distributions. A more intensive third calculation is presented which examines the effects of a fully detailed pulsed power driver model. We also compare simulated electromagnetic fluctuations with direct measurement of radiofrequency electromagnetic fluctuations in a DPF plasma. These comparisons indicate that the fully kinetic model captures the essential physics of these plasmas with high fidelity, and provide further evidence that anomalous resistivity in the plasma arises due to a kinetic instability near the lower hybrid frequency.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Kinetic Simulations of Cu Doping in Chlorinated CdSeTe PV Absorbers
This paper reports on 1D kinetic modeling of p-type doping formation in Cu-doped chlorinated CdSeTe photovoltaic (PV) absorbers with graded Se profiles. Following the recent progress in kinetic simulations of the defect chemistry in stoichiometric CdTe films, this work extends simulation capabilities to a domain of semiconductor alloy films with graded stoichiometries. The defect formation energies and ionization levels, as well as the diffusion and reaction barriers used in the reaction-diffusion equations are calculated from the first principles. A new formalism is employed to account for the uncertainty of the defect formation energies in semiconductor alloys caused by unknown local surrounding. Finally, the developed methodology is used to model a practical fabrication process of p-type doping formation in graded CdSeTe absorbers and to study the effect of alloy stoichiometry on the doping activation.
14 SOLAR ENERGY↗
Kinetic simulations of anomalous resistivity in high-temperature current carrying plasmas
Anomalous resistivity caused by lower hybrid drift instabilities (LHDIs) has been extensively studied in the literature, and has been invoked to explain the rates of magnetic diffusion and plasma profile evolution in the low-density plasma periphery of theta and screw pinch plasmas. In this article, we present a fully kinetic simulation study of LHDIs with a realistic mass ratio using the particle-in-cell code Chicago over a large range of drift speeds, 0.5 < v$^{2}_{di}$/v$^{2}_{ti}$ < 14.5, plasma β, β ≥ 2×10 –5 , and temperature ratios, 0.1 < T e /T i < 10. The resistivity quantified from the simulations is compared with the analytic estimates from the literature and a generalized resistivity expression is presented, which can be applied over the entire range of plasma environments that were simulated in this study. This expression uses the amplitude of the LHDI fluctuations resulting from instability saturation by electron resonance broadening rather than plateau formation, current relaxation, or ion trapping. The generalized resistivity expression is appropriate for full-scale fluid simulations that encounter a large range of plasma conditions, and has better agreement with the resistivity quantified from kinetic simulations than previous expressions in the literature, particularly at low plasma β. Anomalous resistivity has also been considered as an explanation of the rapid magnetic field dissipation during magnetic reconnection events at large plasma β. At a large β, the magnetic field fluctuations of LHDIs can dominate the electrostatic fluctuations and this article investigates enhanced collisionality due to the nonlinear magnetic force, δJ×δB/c.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A new approach to the evaluation and solution of the relativistic kinetic dispersion relation and verification with continuum kinetic simulation
Here, the present work describes a new approach to evaluation and root finding for the kinetic dispersion relation of Langmuir waves, which is central to the analytical understanding of collisionless damping in plasmas. The plasma dispersion function is solved to machine precision using direct integration in the complex plane in combination with an analytic evaluation of the residue to account for the deformation along the Landau contour. To efficiently attain machine precision, the contour is displaced in the complex plane prior to integration, and numerical subtleties related to the placement of the contour are discussed. The approach is generic in that it applies to arbitrary distribution functions, with the present manuscript focused on relativistic cases. Detailed verification of results via direct kinetic simulation in a variety of configuration space dimensions is also presented. Finally, the technique is applied to the challenging case of highly relativistic (i.e. extremely hot) plasmas. Here we show both qualitative agreement with prior work, as well as the disappearance of the Landau root which would have significant implication for real-life observation or experiment.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗