Transverse Momentum Differential Global Analysis of Heavy-Ion Collisions
Not Available
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
SEARCH · Search NASA
Results for “plasma global models”
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.
Electromagnetic energy transport by tearing fluctuations in a self-organized reversed-field pinch plasma
Fluctuation measurements reveal the outward electromagnetic energy flux needed to drive the dynamo electromotive force supporting magnetic self-organization in a reversed-field pinch plasma. The radial Poynting flux due to tearing mode fluctuations is measured with an insertable probe during magnetic relaxation. This flux corresponds to transient power levels much larger than the input power and comparable to the global equilibrium magnetic energy transient loss rate. The probe measurements of this flux are roughly as predicted by a simple Poynting's theorem model upon substitution of equilibrium measurement data.
Physics↗
An Experimental and Numerical Study of Plasma Treatment Effects on Adhesive Joints of Metals, CFRTPs, and Metal–CFRTP Combinations
This work investigates the effects of plasma surface treatment on both the local adherend–adhesive interfacial strength and the global lap shear strength (LSS) of adhesively-bonded joints composed of metals, CFRTPs, and metal-CFRTP combinations, using AA6061 and CFRPPA as representative examples. Experimental and computational results show that plasma treatment improved the ratio of CFRPPA–adhesive interfacial normal strength to adhesive tensile strength (sint,n/scri,t,adhesive) from 0.59–1.00 to 0.91–1.73, leading to an average 90% increase in LSS for CFRPPA-related joints. In contrast, the LSS and sint,n/scri,t,adhesive values of AA6061–AA6061 joints (>=2.27) remained nearly unchanged after plasma treatment. The enhanced joint bonding was further supported by post-failure surface morphologies from three-dimensional profilometry and simulations, showing a transition from dominant interfacial failure to mixed modes involving interfacial, cohesive, and slight fiber-peeling failures. The LSS evolution with respect to sint,n/scri,t,adhesive was numerically established until the saturation of the LSS of various joints, and the behavior of AA6061–CFRPPA T-joints with varying interfacial strengths was also modeled. This study highlights the importance of local interfacial bonding strength, rather than global LSS, for understanding and optimizing the mechanical performance of adhesively-bonded composite–metal and composite-composite structures via different surface modification strategies.
Qiao, Yao (ORCID:0000000265867604)↗
Automated Fitting of a Semi-empirical Multiphase Equation of State for Carbon
The equation of state (EOS) of carbon is important in high explosive, geophysical, and inertial confinement fusion applications. Within the semi-empirical Sesame framework, the EOS of each phase is represented by a sum of cold, vibrational, and thermal electronic Helmholtz free energy contributions. Each phase has ~5-10 independent parameters that are adjusted to reproduce single-phase data (e.g., thermal expansion, isothermal compression) as well as experimentally- and computationally- derived phase boundaries. Manual calibration of the full multiphase EOS is arduous. We present our progress in development of automated EOS parameterization based on minimization of an objective function. Here, this function encodes deviation of model EOS results from experimental/computational benchmarks. Optimization is implemented as a combination of global (particle swarm) and local optimization techniques.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Toroidal modeling of runaway electron loss due to 3D fields in ITER
Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.
3D fields↗
Drift kinetic electrostatic simulations of the edge localized mode heat pulse
In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Scattering by whistler-mode waves during a quiet period perturbed by substorm activity
We study the dynamics of radiation belt electrons during a 10-day quiet period perturbed by substorm activity and preceding a high-speed stream (HSS), aiming at a global description of the radiation belts in L-shell, L in [2, 6], and energy [0.1, 10] MeV. We combine Van Allen Probes observations and Fokker-Planck numerical simulations of pitch-angle diffusion. The Fokker-Planck model uses event-driven pitch angle diffusion coefficients from whistler-mode waves, built from the wave properties and the ambient plasma density measurements from the Van Allen Probes. We first find this event has some similar characteristics to regular quiet times previously studied; a widely extended plasmasphere within which we observe strong and varying whistler-mode waves. These ambient conditions lead to strong pitch-angle scattering, which contributes to the creation of a wide slot region as well as a significant decay of the outer radiation belts, which are observed and qualitatively well simulated. In addition, we find the substorm activity causes short duration (within ± 4h) decay of the plasma density and a lowering amplitude of the whistler-mode waves within the plasmasphere, both causing opposite effects in terms of pitch angle diffusion. This leads to a diminution of pitch-angle diffusion at the time of the main substorm activity. Conversely, whistler-mode waves become enhanced in the time periods between the substorm injections. All effects cumulated, we find an enhancement of pitch angle diffusion by whistler-mode waves above L~4.7 during the 10-day period. This directly relates to the combination of quietness and substorm activity which allows pitch angle diffusing of up to 1 MeV electrons in the outer belt. Relativistic electrons of 1–2 MeV remain trapped in the outer belt, from L~4.7 to L~5.2, forming, in both the observations and the simulations, a distinct pocket of remnant electrons.
79 ASTRONOMY AND ASTROPHYSICS↗
Magnetogenesis from early structure formation due to Yukawa forces
Yukawa interactions can mediate relatively long-range attractive forces between fermions in the early universe. Such a globally attractive interaction creates an instability that can result in the growth of structure in the affected species even during the radiation dominated era. The formation and collapse of fermionic microhalos can create hot fireballs at the sites of the collapsing halos which inject energy into the cosmic plasma. In this paper we study a new phenomena which can take place in such models. We show that the injected energy can be partially converted into primordial magnetic fields and we estimate the correlation scale and the power spectrum of these fields. We show that they may be the seeds of the observed astrophysical magnetic fields.
79 ASTRONOMY AND ASTROPHYSICS↗
Coupling constant corrections in a holographic model of heavy ion collisions with nonzero baryon number density
Sufficiently energetic collisions of heavy ions result in the formation of a droplet of a strongly coupled liquid state of QCD matter known as quark-gluon plasma. By using gauge-gravity duality (holography), a model of a rapidly hydrodynamizing and thermal- izing process like this can be constructed by colliding sheets of energy density moving at the speed of light and tracking the subsequent evolution. In this work, we consider the dual gravitational description of such collisions in the most general bulk theory with a four-derivative gravitational action containing a dynamical metric and a gauge field in five dimensions. Introducing the bulk gauge field enables the analysis of collisions of sheets which carry nonzero “baryon” number density in addition to energy density. Introducing the four-derivative terms enables consideration of such collisions in a gauge theory with finite gauge coupling, working perturbatively in the inverse coupling. While the dynamics of energy and momentum in the presence of perturbative inverse-coupling corrections has been analyzed previously, here we are able to determine the effect of such finite coupling corrections on the dynamics of the density of a conserved global charge, which we take as a model for the dynamics of nonzero baryon number density. In accordance with expectations, as the coupling is reduced we observe that after the collisions less baryon density ends up stopped at mid-rapidity and more of it ends up moving near the lightcone.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Large-scale tearing-mode hazard function analysis with standard matched equilibrium reconstructions
The association between features from standard tokamak equilibrium reconstructions and the onset of n = 1 tearing modes (TMs) is analyzed at scale. The TM onset rate is directly modeled with a ‘hazard’ function which gives the expected number of onsets (per unit time spent) in a given equilibrium parameter region. In particular the different statistical modeling performance achieved for magnetics-only reconstructions and motional Stark effect (MSE) enhanced reconstructions is studied. It is observed that a better hazard model for the TM onset rate can be built with the MSE-enhanced equilibria compared to the matched magnetics-only situation. This advantage disappears if internal profile details are withheld from the matched analysis. Plausibility of the hazard function is further demonstrated with visualizations of global trends in the operational space, and time-traces from specific tokamak discharges. As a result, TMs typically degrade tokamak plasma performance and may lead to plasma termination, motivating this statistical study.
equilibrium↗
Pressure–Strain Interaction as the Energy Dissipation Estimate in Collisionless Plasma
The dissipative mechanism in weakly collisional plasma is a topic that pervades decades of studies without a consensus solution. We compare several energy dissipation estimates based on energy transfer processes in plasma turbulence and provide justification for the pressure–strain interaction as a direct estimate of the energy dissipation rate. The global and scale-by-scale energy balances are examined in 2.5D and 3D kinetic simulations. We show that the global internal energy increase and the temperature enhancement of each species are directly tracked by the pressure–strain interaction. The incompressive part of the pressure–strain interaction dominates over its compressive part in all simulations considered. The scale-by-scale energy balance is quantified by scale filtered Vlasov–Maxwell equations, a kinetic plasma approach, and the lag dependent von Kármán–Howarth equation, an approach based on fluid models. We find that the energy balance is exactly satisfied across all scales, but the lack of a well-defined inertial range influences the distribution of the energy budget among different terms in the inertial range. Therefore, the widespread use of the Yaglom relation in estimating the dissipation rate is questionable in some cases, especially when the scale separation in the system is not clearly defined. In contrast, the pressure–strain interaction balances exactly the dissipation rate at kinetic scales regardless of the scale separation.
79 ASTRONOMY AND ASTROPHYSICS↗
Global X-Ray Properties of the Vela and Puppis A Supernova Remnants
The Vela and Puppis A supernova remnants (SNRs) comprise a large emission region of ∼8° diameter in the soft X-ray sky. The HaloSat CubeSat mission provides the first soft X-ray (0.4–7 keV) observation of the entire Vela SNR and Puppis A SNR region with a single pointing and moderate spectral resolution. HaloSat observations of the Vela SNR are best fit with a two-temperature thermal plasma model consisting of a cooler component with kT{sub 1}=0.19{sub −0.01}{sup +0.01} keV in collisional ionization equilibrium and a hotter component with kT{sub 2}=1.06{sub −0.27}{sup +0.45} keV in nonequilibrium ionization. Observations of the Puppis A SNR are best fit with a single-component plane-parallel shocked plasma model with kT=0.86{sub −0.05}{sup +0.06} keV in nonequilibrium ionization. For the first time, we find the total X-ray luminosities of both components of the Vela SNR spectrum in the 0.5–7 keV energy band to be L{sub X}=4.4{sub −1.4}{sup +1.4}×10{sup 34} erg s{sup −1} for the cooler component and L{sub X}=4.1{sub −1.5}{sup +1.8}×10{sup 34} erg s{sup −1} for the hotter component. We find the total X-ray luminosities of the Vela and Puppis A SNRs to be L{sub X}=8.4×10{sup 34} erg s{sup −1} and L{sub X}=6.7{sub −0.9}{sup +1.1}×10{sup 36} erg s{sup −1}.
79 ASTRONOMY AND ASTROPHYSICS↗
Low- n global ideal MHD instabilities in the CFETR baseline scenario
This article reports an evaluation on the linear ideal magnetohydrodynamic (MHD) stability of the China Fusion Engineering Test Reactor (CFETR) baseline scenario for various first-wall locations. The initial-value code NIMROD and eigen-value code AEGIS are employed in this analysis. Despite the distinctly different approaches in modeling the scrape off layer region, the dominant growth in each of the low-n (n=1–10) modes are consistent between the two codes. We report the higher-n modes are dominated by ballooning modes and localized in the pedestal region, while the lower-n modes have more prominent external kink components and broader mode profiles. The influence of the plasma-vacuum profile and wall shape are also examined using NIMROD. In the presence of a resistive wall, the low-n ideal MHD instabilities are further studied using AEGIS. For the designed first-wall location, the n=1 resistive wall mode is found to be unstable, which could be fully stabilized by uniform toroidal rotation above the 2.9% core Alfven speed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Toroidal phase of post-disruption runaway electron loss to wall in presence of applied 3-D fields in DIII-D
The toroidal distribution of runaway electrons (REs) striking the centre post in DIII-D post-disruption RE final loss events is simulated utilizing the MARS-F code combined with the REORBIT module, by tracing the guiding-centre drift orbits of test REs in the presence of external (applied) plus internal (intrinsic instability) 3-D fields. To better recover the experimental results, three different equilibria with different safety factor profiles are adopted. Dominant resistive magnetohydrodynamic instabilities are found to be m/n=2/1, 3/2 and 1/1 modes (where m and n are the poloidal and toroidal mode numbers). Externally applied n=1 resonant magnetic perturbations (RMPs), with both even and odd parity configuration, are considered in each equilibrium. Locking of MHD instability to the applied RMPs is simulated by minimizing the total perturbed magnetic energy in the plasma. Here, the simulated toroidal locations, where most REs strike in the presence of an edge m/n= 2/1 instability and the n=1 RMPs in even parity, are found to agree well with experiments. Modeling also captures the measured n = 2 RE loss pattern on the wall for experiments where an edge m/n=3/2 instability is dominant instead of the m/n=2/1 mode. The key overall findings include (a) the measured global RE toroidal impact pattern in experiments can be explained by a locked resistive kink instability, (b) the toroidal phase of the locked mode is well predicted by minimizing the total perturbed magnetic energy, and (c) this edge locked mode determines the toroidal phase of the peak RE wall impact.
Bai, Xue [University of California San Diego, La J↗
NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra
This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Effect of neutral interactions on parallel transport and blob dynamics in gyrokinetic scrape-off layer simulations
The effect of neutral interactions on scrape-off layer (SOL) turbulence is investigated in a continuum gyrokinetic code that has been coupled to a continuum kinetic model of neutral transport. This extends the work of a previous paper, which compared two NSTX SOL simulations in simple helical geometry, one with neutrals and one without. The former included electron-impact ionization, charge exchange, and wall recycling. Here, the case with neutrals is compared to a gyrokinetic-only simulation that includes an effective ionization source to separate the effect of sourcing from charge exchange collisions. It is observed that sourcing accounts for many features of the simulated SOL with neutrals, including density and temperature magnitudes and reduced normalized density fluctuations, but differences persist. In particular, a flatter density profile results due to changes in parallel transport when neutral collisions are included, illustrating the importance of neutral drag on global plasma properties. An analysis of coherent turbulent structures, or blobs, in these simulations demonstrates the case with neutrals has slower and larger blobs. Here, a series of seeded blob simulations corroborates the blob velocity observation. In general, the blob motion does not contribute significantly to radial transport in these simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Momentum conservation in current drive and alpha-channeling-mediated rotation drive
Alpha channeling uses waves to extract hot ash from a fusion plasma, transferring energy from the ash to the wave. It has been proposed that this process could create a radial electric field, efficiently driving E × B rotation. However, existing theories ignore the nonresonant particles, which play a critical role in enforcing momentum conservation in quasilinear theory. Because cross field charge transport and momentum conservation are fundamentally linked, this non-consistency throws the rotation drive into question. This paper has two main goals. First, we provide a pedantic and cohesive introduction to the recently developed simple, general, self-consistent quasilinear theory for electrostatic waves that explains the torques which allow for current drive parallel to the magnetic field, and charge extraction across it; a theory that has largely resolved the question of rotation drive by alpha channeling. We show how the theory reveals a fundamental difference between the reaction of nonresonant particles to plane waves that grow in time vs steady-state waves that have a nonuniform spatial structure, allowing rotation drive in the latter case while precluding it in the former, and we review the local and global conservation laws that lead to this result. Second, we provide two new results in support of the theory. First, we provide a novel two-particle Hamiltonian model that rigorously establishes the relationship between charge transport and momentum conservation. Second, we compare the new quasilinear theory to the oscillation-center theories of ponderomotive forces, showing how the latter often obscure the time-dependent nonresonant recoil, but ultimately lead to similar results.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation
Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.
97 MATHEMATICS AND COMPUTING↗