Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLASMA 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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

The evolution of the ion diffusion region during collisionless magnetic reconnection in a force-free current sheet

Two-dimensional particle-in-cell simulation is performed to investigate magnetic reconnection in a force-free current sheet. The results show that the evolution of the ion diffusion region has two different phases. In the first phase, the electrons flow toward the X line along one pair of separatrices and away from the X line along the other pair of separatrices. Therefore, in the ion diffusion region, a distorted quadrupole structure of the out-of-plane magnetic field is formed, which is similar to that of a typical guide field reconnection in the Harris current sheet. In the second phase, the electrons move toward the X line along the separatrices and then flow away from the X line at the inner side of the separatrices. In the ion diffusion region, the out-of-plane magnetic field exhibits a characteristic quadrupole pattern with a good symmetry, which is similar to that of antiparallel reconnection in the Harris current sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

On the electron dynamics during island coalescence in asymmetric magnetic reconnection

We present an analysis of the electron dynamics during rapid island merging in asymmetric magnetic reconnection. We consider a doubly periodic system with two asymmetric transitions. The upper layer is an asymmetric Harris sheet of finite width perturbed initially to promote a single reconnection site. The lower layer is a tangential discontinuity that promotes the formation of many X-points, separated by rapidly merging islands. Across both layers, the magnetic field and the density have a strong jump, but the pressure is held constant. Our analysis focuses on the consequences of electron energization during island coalescence. We focus first on the parallel and perpendicular components of the electron temperature to establish the presence of possible anisotropies and non-gyrotropies. Thanks to the direct comparison between the two different layers simulated, we can distinguish three main types of behavior characteristic of three different regions of interest. The first type represents the regions where traditional asymmetric reconnections take place without involving island merging. The second type of regions instead shows reconnection events between two merging islands. Finally, the third regions identify the regions between two diverging island and where typical signature of reconnection is not observed. Electrons in these latter regions additionally show a flat-top distribution resulting from the saturation of a two-stream instability generated by the two interacting electron beams from the two nearest reconnection points. Finally, the analysis of agyrotropy shows the presence of a distinct double structure laying all over the lower side facing the higher magnetic field region. This structure becomes quadrupolar in the proximity of the regions of the third type. The distinguishing features found for the three types of regions investigated provide clear indicators to the recently launched Magnetospheric Multiscale NASA mission for investigating magnetopause reconnection involving multiple islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Magnetic reconnection mediated by hyper-resistive plasmoid instability

Magnetic reconnection mediated by the hyper-resistive plasmoid instability is studied with both linear analysis and nonlinear simulations. The linear growth rate is found to scale as S{sub H}{sup 1/6} with respect to the hyper-resistive Lundquist number S{sub H}≡L{sup 3}V{sub A}/η{sub H}, where L is the system size, V{sub A} is the Alfvén velocity, and η{sub H} is the hyper-resistivity. In the nonlinear regime, reconnection rate becomes nearly independent of S{sub H}, the number of plasmoids scales as S{sub H}{sup 1/2}, and the secondary current sheet length and width both scale as S{sub H}{sup −1/2}. These scalings are consistent with a heuristic argument assuming secondary current sheets are close to marginal stability. The distribution of plasmoids as a function of the enclosed flux ψ is found to obey a ψ{sup −1} power law over an extended range, followed by a rapid fall off for large plasmoids. These results are compared with those from resistive magnetohydrodynamic studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Dynamics of laser mass-limited foil interaction at ultra-high laser intensities

By using three-dimensional particle-in-cell simulations with synchrotron radiation damping incorporated, dynamics of ultra-intense laser driven mass-limited foils is presented. When a circularly polarized laser pulse with a peak intensity of ∼10{sup 22} W/cm{sup 2} irradiates a mass-limited nanofoil, electrons are pushed forward collectively and a strong charge separation field forms which acts as a “light sail” and accelerates the protons. When the laser wing parts overtake the foil from the foil boundaries, electrons do a betatron-like oscillation around the center proton bunch. Under some conditions, betatron-like resonance takes place, resulting in energetic circulating electrons. Finally, bright femto-second x rays are emitted in a small cone. It is also shown that the radiation damping does not alter the foil dynamics radically at considered laser intensities. The effects of the transverse foil size and laser polarization on x-ray emission and foil dynamics are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Kinetic simulations of plasmoid chain dynamics

The dynamics of a plasmoid chain is studied with three dimensional Particle-in-Cell simulations. The evolution of the system with and without a uniform guide field, whose strength is 1/3 the asymptotic magnetic field, is investigated. The plasmoid chain forms by spontaneous magnetic reconnection: the tearing instability rapidly disrupts the initial current sheet generating several small-scale plasmoids that rapidly grow in size coalescing and kinking. The plasmoid kink is mainly driven by the coalescence process. It is found that the presence of guide field strongly influences the evolution of the plasmoid chain. Without a guide field, a main reconnection site dominates and smaller reconnection regions are included in larger ones, leading to an hierarchical structure of the plasmoid-dominated current sheet. On the contrary in presence of a guide field, plasmoids have approximately the same size and the hierarchical structure does not emerge, a strong core magnetic field develops in the center of the plasmoid in the direction of the existing guide field, and bump-on-tail instability, leading to the formation of electron holes, is detected in proximity of the plasmoids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Global fluid turbulence simulations in the scrape-off layer of a stellarator island divertor

Isothermal fluid turbulence simulations have been performed in the edge and scrape-off layer (SOL) of an analytic stellarator configuration with an island divertor, thereby providing numerical insight into edge turbulence in regions around islands in a stellarator. The steady-state transport follows the a curvature drive that is inverse to the major radius ( $1/R$ ) toward the outboard side, but large fluctuations are present throughout the island divertor region, with the average wavelength of similar size to the island width. The system exhibits a prominent $m=2$ , $n=5$ mode, where m is the poloidal mode number and n is the toroidal mode number, although other modes are present. The amplitude and radial extent of the density fluctuations are similar throughout the edge and SOL, but can decrease near island O-points. The fluctuations exhibit a predominantly positive skewness on the outboard midplane, indicating blob-like perturbations for the transport into the outer SOL. It is determined that a point on the separatrix is generally more correlated with regions outside of the SOL than a nearby reference point which does not lie on the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Kinetic neoclassical transport in the H-mode pedestal

Multi-species kinetic neoclassical transport through the QH-mode pedestal and scrape-off layer on DIII-D is calculated using XGC0, a 5D full-f particle-in-cell drift-kinetic solver with self-consistent neutral recycling and sheath potentials. Quantitative agreement between the flux-driven simulation and the experimental electron density, impurity density, and orthogonal measurements of impurity temperature and flow profiles is achieved by adding random-walk particle diffusion to the guiding-center drift motion. The radial electric field (E{sub r}) that maintains ambipolar transport across flux surfaces and to the wall is computed self-consistently on closed and open magnetic field lines and is in excellent agreement with experiment. The E{sub r} inside the separatrix is the unique solution that balances the outward flux of thermal tail deuterium ions against the outward neoclassical electron flux and inward pinch of impurity and colder deuterium ions. Particle transport in the pedestal is primarily due to anomalous transport, while the ion heat and momentum transport are primarily due to the neoclassical transport. The full-f treatment quantifies the non-Maxwellian energy distributions that describe a number of experimental observations in low-collisionallity pedestals on DIII-D, including intrinsic co-I{sub p} parallel flows in the pedestal, ion temperature anisotropy, and large impurity temperatures in the scrape-off layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Particle-in-cell simulation of multipactor discharge on a dielectric in a parallel-plate waveguide

An original 2D3V (two-dimensional in coordinate space and three-dimensional in velocity space) particle-in-cell code has been developed for simulation of multipactor discharge on a dielectric in a parallelplate metal waveguide with allowance for secondary electron emission (SEE) from the dielectric surface and waveguide walls, finite temperature of secondary electrons, electron space charge, and elastic and inelastic scattering of electrons from the dielectric and metal surfaces. The code allows one to simulate all stages of the multipactor discharge, from the onset of the electron avalanche to saturation. It is shown that the threshold for the excitation of a single-surface multipactor on a dielectric placed in a low-profile waveguide with absorbing walls increases as compared to that in the case of an unbounded dielectric surface due to escape of electrons onto the waveguide walls. It is found that, depending on the microwave field amplitude and the SEE characteristics of the waveguide walls, the multipactor may operate in two modes. In the first mode, which takes place at relatively low microwave amplitudes, a single-surface multipactor develops only on the dielectric, the surface of which acquires a positively potential with respect to the waveguide walls. In the second mode, which occurs at sufficiently high microwave intensities, a single-surface multipactor on the dielectric and a two-surface multipactor between the waveguide walls operate simultaneously. In this case, both the dielectric surface and the interwall space acquire a negative potential. It is shown that electron scattering from the dielectric surface and waveguide walls results in the appearance of high-energy tails in the electron distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Numerical verification of Orbit and Nemato codes for magnetic topology diagnosis

We present the results of a benchmark study involving the Orbit and Nemato codes. The two codes have been used to compare magnetic structure in a reversed-field pinch (RFP), where conserved magnetic structures/islands appear both in the core (dominated by m = 1 modes) and in the edge (dominated by m = 0 modes). As input, a snapshot of a 3D nonlinear MHD visco-resistive simulation (produced by the SpeCyl code) has been used. The first test is given by the reconstruction via Poincaré surface of section plot of an island generated by a single mode. In this case, the magnetic field topology corresponds to a time-independent Hamiltonian and shows conserved flux-surfaces used as a reference. Both codes successfully yield field lines which follow the same flux surfaces, in both the m = 1 and m = 0 cases. The benchmark between the codes has then been extended to a more complex configuration with chaotic magnetic field, using as input a fully 3D multiple mode RFP condition, characterized by the typical chain of edge magnetic islands providing a transport barrier. Finally, a quantitative benchmark has been performed, using the same 3D input, by estimating the correlation length of the magnetic field line in a bounded stochastic domain. The conclusion is that both codes yield consistent Poincaré plot structure on one hand, and very good quantitative agreement in correlation length estimate. This gives confidence to the application of both codes to magnetic topology in the RFX-mod device, for which they are used routinely, as well as for a generic fusion device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Effects of a Guide Field on the Larmor Electric Field and Upstream Electron Temperature Anisotropy in Collisionless Asymmetric Magnetic Reconnection

We perform the first study of the properties of the Larmor electric field (LEF) in collisionless asymmetric magnetic reconnection in the presence of an out-of-plane (guide) magnetic field for different sets of representative upstream parameters at Earth’s dayside magnetopause with an ion temperature greater than the electron temperature (the ion-to-electron temperature ratio fixed at 2) using two-dimensional particle-in-cell simulations. We show that the LEF does persist in the presence of a guide field. We study how the LEF thickness and strength change as a function of guide field and the magnetospheric temperature and reconnecting magnetic field strength. We find that the thickness of the LEF structure decreases, while its magnitude increases when a guide field is added to the reconnecting magnetic field. The added guide field makes the Larmor radius smaller, so the scaling with the magnetospheric ion Larmor radius is similar to that reported for the case without a guide field. Note, however, that the physics causing the LEF is not well understood, so future work in other parameter regimes is needed to fully predict the LEF for arbitrary conditions. We also find that a previously reported upstream electron temperature anisotropy arises in the vicinity of the LEF region both with and without a guide field. We argue that the generation of the anisotropy is linked to the existence of the LEF. The LEF can be used in combination with the electron temperature anisotropy as a signature to effectively identify dayside reconnection sites in observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Microsecond ramp compression of a metallic liner driven by a 5 MA current on the SPHINX machine using a dynamic load current multiplier pulse shaping

SPHINX is a 6 MA, 1-μs Linear Transformer Driver (LTD) operated by the CEA Gramat (France) and primarily used for imploding Z-pinch loads for radiation effects studies. Among the options that are currently being evaluated to improve the generator performances are an upgrade to a 20 MA, 1-μs LTD machine and various power amplification schemes, including a compact Dynamic Load Current Multiplier (DLCM). A method for performing magnetic ramp compression experiments, without modifying the generator operation scheme, was developed using the DLCM to shape the initial current pulse in order to obtain the desired load current profile. In this paper, we discuss the overall configuration that was selected for these experiments, including the choice of a coaxial cylindrical geometry for the load and its return current electrode. We present both 3-D Magneto-hydrodynamic and 1D Lagrangian hydrodynamic simulations which helped guide the design of the experimental configuration. Initial results obtained over a set of experiments on an aluminium cylindrical liner, ramp-compressed to a peak pressure of 23 GPa, are presented and analyzed. Details of the electrical and laser Doppler interferometer setups used to monitor and diagnose the ramp compression experiments are provided. In particular, the configuration used to field both homodyne and heterodyne velocimetry diagnostics in the reduced access available within the liner's interior is described. Current profiles measured at various critical locations across the system, particularly the load current, enabled a comprehensive tracking of the current circulation and demonstrate adequate pulse shaping by the DLCM. The liner inner free surface velocity measurements obtained from the heterodyne velocimeter agree with the hydrocode results obtained using the measured load current as the input. An extensive hydrodynamic analysis is carried out to examine information such as pressure and particle velocity history profiles or magnetic diffusion across the liner. The potential of the technique in terms of applications and achievable ramp pressure levels lies in the prospects for improving the DLCM efficiency through the use of a closing switch (currently under development), reducing the load dimensions and optimizing the diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Diffusive shock re-acceleration

Here, we have performed two-dimensional hybrid simulations of non-relativistic collisionless shocks in the presence of pre-existing energetic particles (‘seeds’); such a study applies, for instance, to the re-acceleration of galactic cosmic rays (CRs) in supernova remnant (SNR) shocks and solar wind energetic particles in heliospheric shocks. Energetic particles can be effectively reflected and accelerated regardless of shock inclination via a process that we call diffusive shock re-acceleration. We determine that re-accelerated seeds can drive the streaming instability in the shock upstream and produce effective magnetic field amplification. This can later trigger the injection of thermal protons even at oblique shocks that ordinarily cannot inject thermal particles. We characterize the current in reflected seeds, finding that it tends to a universal value $J\simeq en_{\text{CR}}v_{\text{sh}}$ , where $en_{\text{CR}}$ is the seed charge density and $v_{\text{sh}}$ is the shock velocity. When applying our results to SNRs, we find that the re-acceleration of galactic CRs can excite the Bell instability to nonlinear levels in less than ${\sim}10~\text{yr}$ , thereby providing a minimum level of magnetic field amplification for any SNR shock. Finally, we discuss the relevance of diffusive shock re-acceleration also for other environments, such as heliospheric shocks, galactic superbubbles and clusters of galaxies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

He{sup 2+} HEATING VIA PARAMETRIC INSTABILITIES OF PARALLEL PROPAGATING ALFVÉN WAVES WITH AN INCOHERENT SPECTRUM

The preferential heating of heavy ions in the solar corona and solar wind has been a long-standing hot topic. In this paper we use a one-dimensional hybrid simulation model to investigate the heating of He{sup 2+} particles during the parametric instabilities of parallel propagating Alfvén waves with an incoherent spectrum. The evolution of the parametric instabilities has two stages and involves the heavy ion heating during the entire evolution. In the first stage, the density fluctuations are generated by the modulation of the pump Alfvén waves with a spectrum, which then results in rapid coupling with the pump Alfvén waves and the cascade of the magnetic fluctuations. In the second stage, each pump Alfvén wave decays into a forward density mode and a backward daughter Alfvén mode, which is similar to that of a monochromatic pump Alfvén wave. In both stages the perpendicular heating of He{sup 2+} particles occurs. This is caused by the cyclotron resonance between He{sup 2+} particles and the high-frequency magnetic fluctuations, whereas the Landau resonance between He{sup 2+} particles and the density fluctuations leads to the parallel heating of He{sup 2+} particles. The influence of the drift velocity between the protons and the He{sup 2+} particles on the heating of He{sup 2+} particles is also discussed in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

The relation between reconnected flux, the parallel electric field, and the reconnection rate in a three-dimensional kinetic simulation of magnetic reconnection

We investigate the distribution of parallel electric fields and their relationship to the location and rate of magnetic reconnection in a large particle-in-cell simulation of 3D turbulent magnetic reconnection with open boundary conditions. The simulation's guide field geometry inhibits the formation of simple topological features such as null points. Therefore, we derive the location of potential changes in magnetic connectivity by finding the field lines that experience a large relative change between their endpoints, i.e., the quasi-separatrix layer. We find a good correspondence between the locus of changes in magnetic connectivity or the quasi-separatrix layer and the map of large gradients in the integrated parallel electric field (or quasi-potential). Furthermore, we investigate the distribution of the parallel electric field along the reconnecting field lines. We find the reconnection rate is controlled by only the low-amplitude, zeroth and first–order trends in the parallel electric field while the contribution from fluctuations of the parallel electric field, such as electron holes, is negligible. The results impact the determination of reconnection sites and reconnection rates in models and in situ spacecraft observations of 3D turbulent reconnection. It is difficult through direct observation to isolate the loci of the reconnection parallel electric field amidst the large amplitude fluctuations. However, we demonstrate that a positive slope of the running sum of the parallel electric field along the field line as a function of field line length indicates where reconnection is occurring along the field line.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Experimental room temperature hohlraum performance study on the National Ignition Facility [Experimental evidence for improved performance in room temperature hohlraums on the National Ignition Facility]

Room temperature or “warm” (273 K) indirect drive hohlraum experiments have been conducted on the National Ignition Facility with laser energies up to 1.26 MJ and compared to similar cryogenic or “cryo” (~20 K) experiments. Warm experiments use neopentane (C 5 H 12 ) as the low pressure hohlraum fill gas instead of helium, and propane (C 3 H 8 ) to replace the cryogenic DT or DHe3 capsule fill. The increased average Z of the hohlraum fill leads to increased inverse bremsstrahlung absorption and an overall hotter hohlraum plasma in simulations. The cross beam energy transfer (CBET) from outer laser beams (pointed toward the laser entrance hole) to inner beams (pointed at the equator) was inferred indirectly from measurements of Stimulated Raman Scattering (SRS). These experiments show that a similar hot spot self-emission shape can be produced with less CBET in warm hohlraums. The measured inner cone SRS reflectivity (as a fraction of incident power neglecting CBET) is ~2.5× less in warm than cryo shots with similar hot spot shapes, due to a less need for CBET. The measured outer-beam stimulated the Brillouin scattering power that was higher in the warm shots, leading to a ceiling on power to avoid the optics damage. These measurements also show that the CBET induced by the flow where the beams cross can be effectively mitigated by a 1.5 Å wavelength shift between the inner and outer beams. A smaller scale direct comparison indicates that warm shots give a more prolate implosion than cryo shots with the same wavelength shift and pulse shape. Lastly, the peak radiation temperature was found to be between 5 and 7 eV higher in the warm than the corresponding cryo experiments after accounting for differences in backscatter.

42 ENGINEERING↗
DOE OSTI2018

Simulations of the Stable Expansion of a Monoelement Three-Component Explosive-Emission Plasma

The results of numerical simulations of the stable expansion stage of a micro explosive plasma plume are presented. A three-liquid hydrodynamic model was used with corrections allowing its application to the consistent description of the transfer of particles in a monoelement plasma with ions with two different charges and of the electron transfer in a changing plasma–collector vacuum gap. The evolution of concentration profiles, the drift velocity, and electron and ion temperatures in plasma are analyzed depending on conditions of a collector. It is shown in the model that the velocity separation of ions with different charges occurs efficiently in the absence of an external electric field, whereas in a current-carrying plasma, the ions with different charges move with almost identical velocities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗