Search NASA⌕ Search

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.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Global stochastic optimization of stellarator coil configurations

In the construction of a stellarator, the manufacturing and assembling of the coil system is a dominant cost. These coils need to satisfy strict engineering tolerances, and if those are not met the project could be cancelled as in the case of the National Compact Stellarator Experiment (NCSX) project. Therefore, our goal is to find coil configurations that increase construction tolerances without compromising the performance of the magnetic field. In this paper, we develop a gradient-based stochastic optimization model which seeks robust stellarator coil configurations in high dimensions. In particular, we design a two-step method: first, we perform an approximate global search by a sample efficient trust-region Bayesian optimization; second, we refine the minima found in step one with a stochastic local optimizer. To this end, we introduce two stochastic local optimizers: BFGS applied to the sample average approximation; and Adam, equipped with a control variate for variance reduction. Numerical simulations performed on a W7-X-like coil configuration demonstrate that our global optimization approach finds a variety of promising local solutions at less than 0.1% of the cost of previous work, which considered solely local stochastic optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Linear simulation of magnetohydrodynamic plasma response to three-dimensional magnetic perturbations in high-β P plasmas

In this work, we report the numerical analyses of linear magnetohydrodynamics (MHD) plasma response to applied three-dimensional magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β P (poloidal beta) plasmas, utilizing the extended-MHD code M3D-C1, with the purpose of realizing a better understanding of the existing experiment in which the n=3 MPs were applied to such high-β P plasmas attempting to control large amplitude type-I ELMs. Such high-β P plasmas obtained at the DIII-D tokamak feature an upper-biased double null configuration, a high edge safety factor q 95 ~6.4, and a stable internal transport barrier (ITB) leading to relatively high core pressures. Single-fluid simulations show that the plasma response to n=3 MPs, including both non-resonant/kinking and resonant components, is significantly weaker than that to n=1 or 2 MPs. To survey the impact of q 95 on plasma response to applied MPs, the SEGWAY (Self-consistent Equilibrium Generating Workflow for AnalYsis) module, developed in the OMFIT integrated modelling framework, is employed to generate a series of equilibria with a wide range of q 95 while other key parameters including the normalized beta, electron density at pedestal top, and plasma shape are kept fixed. Compared to the vacuum response, single-fluid M3D-C1 simulations predict a much more significant decrease of resonant plasma response to the applied n=3 MPs at the maximum penetration radii as q 95 increases. In contrast to single-fluid simulation results showing resonant penetration occurs only near the pedestal top where the E×B toroidal rotation frequency is zero, two-fluid simulations show two comparable resonant penetrations locating near the pedestal top and the ITB foot, where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot may be responsible for the observed formation of a staircase structure in both electron density and temperature profiles and thereby a considerable deterioration of global plasma performance when MPs are applied in high-β P plasmas. Motivated by this numerical work, we provide some ideas for the future research, with the purpose of realizing effective ELM control in such high-β P plasmas on the DIII-D and EAST devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Predicting nonresonant pressure-driven MHD modes in equilibria with low magnetic shear

Nonresonant internal modes can be difficult to anticipate as there is no resonant surface in the plasma. However, equilibria that are unstable to multiple nonresonant magnetohydrodynamic (MHD) modes may be more prone to global loss of confinement since these instabilities generate spatially extended linear displacements, potentially enhancing magnetic field line chaos via nonlinear interactions. Here, we successfully predict the unstable nonresonant pressure-driven modes for equilibria with zero shear in the plasma core, irrational q on axis, and a central pressure gradient, which is consistent with pre-crash profiles in sawtoothing tokamak plasmas in the large-aspect-ratio limit. A criterion for identifying nonresonant modes most likely to be unstable is developed from the convergents of the continued fraction representation of q 0 . A higher-order analysis of the standard Energy Principle reveals the conditions under which these modes are expected to dominate. Linear growth rate spectra, as a function of toroidal mode number (up to n = 30), calculated using the initial-value extended-MHD code, M3D-C1, recover the characteristic dependence observed for ideal infernal modes. Nonresonant modes have also been invoked in some ideal sawtooth crash models. Finally, this work provides a mechanism to predict the mode numbers of infernal modes and, potentially, the width of some post-sawtooth-crash profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Efficient facemask decontamination via forced ozone convection

The COVID-19 crisis has taken a significant toll on human life and the global economy since its start in early 2020. Healthcare professionals have been particularly vulnerable because of the unprecedented shortage of Facepiece Respirators (FPRs), which act as fundamental tools to protect the medical staff treating the coronavirus patients. In addition, many FPRs are designed to be disposable single-use devices, creating an issue related to the generation of large quantities of non-biodegradable waste. In this contribution, we describe a plasma-based decontamination technique designed to circumvent the shortages of FPRs and alleviate the environmental problems posed by waste generation. The system utilizes a Dielectric Barrier Discharge (DBD) to generate ozone and feed it through the fibers of the FPRs. The flow-through configuration is different than canonical ozone-based sterilization methods, in which the equipment is placed in a sealed ozone-containing enclosure without any flow through the mask polymer fibers. We demonstrate the rapid decontamination of surgical masks using Escherichia coli ( E. coli ) and Vesicular Stomatitis Virus (VSV) as model pathogens, with the flow-through configuration providing a drastic reduction in sterilization time compared to the canonical approach. We also demonstrate that there is no deterioration in mask structure or filtration efficiency resulting from sterilization. Finally, we show that this decontamination approach can be implemented using readily available tools, such as a plastic box, a glass tube, few 3D printed components, and the high-voltage power supply from a plasma globe toy. The prototype assembled for this study is portable and affordable, with effectiveness comparable to that of larger and more expensive equipment.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Global Characteristics of the IBEX Ribbon and Its Temporal Variation: Implications for Relationships with Solar Cycle Variability and the Interstellar Magnetic Field

The Interstellar Boundary Explorer (IBEX) has revealed the existence of the energetic neutral atom (ENA) ribbon, a narrow and nearly circular band of enhanced flux that enables probing plasma and magnetic field conditions outside the heliopause (HP). In this study, we investigate the global structure and temporal variation of the IBEX ribbon using 14 yr of IBEX-Hi ENA observations. We derive geometric parameters of the ribbon: center location, angular radius, eccentricity, and major axis orientation. We find that the ribbon exhibits systematic variations tied to the solar cycle. The ribbon center shows a solar cycle–modulated latitudinal motion at low energies. The angular radius increases with ENA energy and displays a pronounced enlargement between 2013 and 2018, coinciding with changes in high-latitude fast solar wind from polar coronal holes. The eccentricity remains consistently around 0.3 with some temporal evolution, while the major axis rotates gradually over time, becoming increasingly perpendicular to the ecliptic meridional plane. A notable hemispheric asymmetry has emerged over time: the ribbon in the northern hemisphere remains remarkably stable in position, whereas the southern ribbon exhibits temporal migration. By combining these results with solar wind observations, global heliosphere, and ribbon modeling, we suggest that this asymmetry is explained by a coupled effect of magnetic field draping and the latitude-dependent timing of fast solar wind from polar coronal holes. These results establish the most detailed long-term geometric characterization of the IBEX ribbon to date and provide observational support for models of secondary ENA production and magnetic field draping outside the HP.

Noh, Sung Jun [New Mexico Consortium, Los Alamos, ↗
DOE OSTI2025

Imaging nuclear shape through anisotropic and radial flow in high-energy heavy-ion collisions

Most atomic nuclei exhibit ellipsoidal shapes characterized by quadrupole deformation β 2 and triaxiality γ, and sometimes even a pear-like octupole deformation β 3 . The STAR experiment introduced a new ‘imaging-by-smashing’ technique to image the nuclear global shape by colliding nuclei at ultra-relativistic speeds and analyzing outgoing debris. Features of nuclear shape manifest in collective observables like anisotropic flow vn and radial flow via mean transverse momentum [p T ]. We present new measurements of the variances of v n (n = 2, 3, and 4) and [p T ], and the covariance of $v$$^{2}_{n}$ with [p T ], in collisions of highly deformed 238 U and nearly spherical 197 Au. Ratios of these observables between the two systems effectively suppress common final-state effects, isolating the strong impact of uranium’s deformation. By comparing results with state-of-the-art hydrodynamic model calculations, we extract β 2U and γ U values consistent with those deduced from low-energy nuclear structure measurements. Measurements of v 3 and its correlation with [p T ] also provide the first experimental suggestion of a possible octupole deformation for 238 U. These findings provide significant support for using high-energy collisions to explore nuclear shapes on femtosecond timescales, with implications for both nuclear structure and quark-gluon plasma studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving the linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the resistive DCON code [A. H. Glasser, Z. R. Wang, and J.-K. Park, Phys. Plasmas 23, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Additionally, the solution of the ideal outer region is critical to the asymptotic matching and is challenging in toroidal geometry due to the singular power series solution at the resonant surfaces. Thus, a systematic verification of the outer region Δ' matrix is made by reproducing the well-known analytical Δ' result in Furth et al. [Phys. Fluids 16, 1054–1063 (1073)] and by making a quantitative benchmark with the PEST3 code [A. Pletzer and R. L. Dewar, J. Plasma Phys. 45, 427–451 (1991)]. Finally, the reconstructed numerical solution of the resistive plasma response from the toroidal matching matrix is presented. Compared with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Metal-induced oxidative stress and human plasma protein oxidation after SARS-CoV-2 infection

Pathogenesis of COVID-19 by SARS-CoV-2 resulted in a global pandemic and public health emergency in 2020. Viral infection can induce oxidative stress through reactive oxygen species (ROS). Inflammation and environmental stress are major sources of oxidative stress after infection. Micronutrients such as iron, copper, zinc, and manganese play various roles in human tissues and their imbalance in blood can impact immune responses against pathogens including SARS CoV-2. We hypothesized that alteration of free metal ions during infection and metal-catalyzed oxidation plays a critical role towards pathogenesis after infection. We analyzed convalescent and hospitalized COVID-19 patient plasma using orthogonal analytical techniques to determine redox active metal concentrations, overall protein oxidation, oxidative modifications, and protein levels via proteomics to understand the consequences of metal-induced oxidative stress in COVID-19 plasma proteins. Metal analysis using ICP-MS showed significantly greater concentrations of copper in COVID-19 plasma compared to healthy controls. We demonstrate significantly greater total protein carbonylation, other oxidative modifications, and deamidation of plasma proteins in COVID-19 plasma compared to healthy controls. Proteomics analysis showed that levels of redox active proteins including hemoglobulin were elevated in COVID-19 plasma. Molecular modeling concurred with potential interactions between iron binding proteins and SARS CoV-2 surface proteins. Overall, increased levels of redox active metals and protein oxidation indicate that oxidative stress-induced protein oxidation in COVID-19 may be a consequence of the interactions of SARS-CoV-2 proteins with host cell metal binding proteins resulting in altered cellular homeostasis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the Resistive DCON code [A.H. Glasser, Z.R. Wang and J.-K. Park, Physics of Plasmas, \textbf{23}, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Since the solution of ideal outer region is critical to the asymptotic matching and is challenging to be solved in the toroidal geometry due to the singular power series solution at the resonant surfaces, systematic verification of the outer region $\Delta^\prime$ matrix is made by reproducing the well known analytical $\Delta^{\prime}$ result in [H.P. Furth, P.H. Rutherford and H. Selberg, The Physics of Fluids, \textbf{16}, 1054-1063 (1073)] as well as by making a quantitative benchmark with the PEST3 code [A. Pletzer and R.L. Dewar, J. Plasma Physics, \textbf{45}, 427-451 (1991)]. Finally, the reconstructed numerical solution of resistive plasma response from the toroidal matching matrix is presented. Comparing with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

MHD simulation↗
DOE OSTI2022

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: B. A reduced model based on reversed-direction two point modeling

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter–erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding (Pitts et al 2019). Part A of the present study reports new results for largely the same ITER cases, confirming the strong correlation reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not as yet been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne (Ne s –1 ), and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In this part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction two point modeling, Rev2PM. The input to the reduced model is a value of the variable pair (T e,t @q ⊥,pk , Inj Ne )for a chosen case and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , the electron density at the X-point, n e,Xpt , etc. The reduced model was quantitatively characterized using one half of the code cases; it was then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Description of global EGAM in the maximum of local frequency during current ramp-up discharges in DIII-D

Energetic-particle-induced geodesic acoustic modes, EGAMs (Fu, Phys. Rev. Let., vol. 101, 2008, pp. 185002), driven by neutral beam injection (NBI), have been observed in many DIII-D tokamak experiments (Nazikian et al., Phys. Rev. Lett., vol. 101, 2008, pp. 185001). This mechanism has been theoretically investigated in (Qiu et al., Plasma Phys. Control. Fusion, vol. 52, 2010, pp. 095003), using a sharp energetic particle distribution function, and in (Qu et al., Plasma Phys. Control. Fusion, vol. 59, 2017, pp. 055018), where the dispersion relation and eigenmode behaviour were obtained for the situation of early beam scenario, that is, for times smaller than the beam slowing down time. In this work, we extend these studies determining the eigenmode for beyond the slowing down time, in a scenario with reverse safety factor q profile, where a small concentration of energetic ions can produce an off-axis maximum in the GAM dispersion relation. Here, the characteristics of EGAM are analytically studied with the drift kinetic equation together with the MHD code NOVA. The toroidal energetic ion transit frequency, coupled with the GAM frequency, produces the maximum in the dispersion relation where the eigenmode can be found. The quantitative correspondence of experimental results with the predictions of the proposed model is analysed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Generation of Subauroral Longitudinally Extended Emissions Following Intensifications of the Poleward Boundary of the Substorm Bulge and Streamer Production

Observations showing the development of Subauroral Longitudinally Extended Emigerssions on a global scale are presented. It is demonstrated that they occur as a separation of an east-west arc-like band of luminosity that detaches away from the equatorward edge of the auroral distribution following episodes of auroral streamer production. They persist for time intervals on the order of ~30 min and devolve into patchy disjointed segments before they fade. Emissions are seen in both 557.7 nm OI and 391.4 nm N + 2 1NG lines, and 630.0 nm emissions are often observed equatorward. The 391.4 nm emissions are typically weaker and fade away more quickly than the 557.7 nm emissions, suggesting that an auroral precipitation source is initially present but is rapidly depleted as the forms age. All cases are associated with enhanced Subauroral Polarization Stream (SAPS) flows, and one event shows clear association with large-scale Giant Undulations (GUs) and the formation of spur-like forms that fold equatorward under the main oval. A model is proposed for the production of subauroral longitudinally extended forms in which nonlinear growth of SAPS-induced surface waves on the plasmapause results in a disruption of the boundary. It is suggested that hot plasma of plasma sheet origin becomes entrained in the plasmasphere to produce transient precipitation-associated auroral emissions that may decay into STEVE emissions and that cold dense plasma from the plasmasphere becomes entrained onto open drift paths to feed long-lived drainage plumes. This process may occur quasi periodically during intense substorms and storm-time conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Interplay between beam-driven chirping modes and plasma confinement transitions in spherical tokamak ST40

Experiments on the high field spherical tokamak ST40 have led to the recent observation of interplay between beam-driven modes of sweeping frequency (chirping modes) and transitions to the enhanced global confinement regime (H-mode) and back to the low confinement regime (L-mode). The H-modes of plasma confinement are identified from decreased intensity of D α signal and from clear distinctions in the edge gradients of the visible plasma boundary (observed as a sharp plasma edge in camera images). The beam-driven chirping modes, identified as ideal magnetohydrodynamics beta-induced Alfvén acoustic eigenmodes modes, are observed in Mirnov coil signals, interferometry, and soft x-ray diagnostics. A moderate amplitude ‘primer’ chirping mode usually precedes an H–L transition. This is followed by a ‘dominant’ chirping mode with higher amplitude during the L-mode. The L–H transition back to the improved confinement occurs on a longer time scale of tens of ms, consistent with the slowing down time scale of fast beam ions. A dramatic decrease in toroidal plasma rotation is systematically observed associated with chirping modes sweeping down to zero frequency. Resonance maps built for the beam-driven chirping modes with the ASCOT (accelerated simulation of charged particle orbits in torodoial devices) code show that the resonant beam ions have orbits near the trapped-passing boundary. The ASCOT modelling assesses how losses of the resonant fast ions caused by the chirping modes with high enough amplitude modify the torque, potentially affecting the plasma rotation.

magnetic confinement fusion↗
DOE OSTI2021

Hybrid simulations of sub-cyclotron compressional and global Alfvén eigenmode stability in spherical tokamaks

A comprehensive numerical study has been conducted in order to investigate the stability of beam-driven, sub-cyclotron-frequency compressional Alfvén eigenmodes (CAEs) and global Alfvén eigenmodes (GAEs) in low-aspect-ratio plasmas for a wide range of beam parameters. The presence of CAEs and GAEs has previously been linked to anomalous electron temperature profile flattening at high beam powers in NSTX experiments, prompting a further examination of the conditions necessary for their excitation. Linear simulations have been performed with the hybrid MHD–kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Additionally, three distinct types of modes were found in the simulations—co-CAEs, cntr-GAEs, and co-GAEs—with differing spectral and stability properties. The simulations revealed that unstable GAEs are more ubiquitous than unstable CAEs, which is consistent with experimental observations, as they are excited at lower beam energies and generally have larger growth rates. Local analytic theory is used to explain key features of the simulation results, including the preferential excitation of different modes based on beam injection geometry and the growth rate dependence on the beam injection velocity, critical velocity, and degree of velocity space anisotropy. The background damping rate is inferred from the simulations and analytically estimated for relevant sources absent from the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Global gyrokinetic simulation of neoclassical ambipolar electric field and its effects on microturbulence in W7-X stellarator

Global neoclassical simulations of a model equilibrium of the W7-X stellarator find an ambipolar electric field with either an ion root or an electron root solution depending on the temperature ratio between electrons and ions. The ambipolar electric field is then used as an equilibrium field in the turbulence simulations of ion temperature gradient (ITG) instability. The shear of the ambipolar electric field has modest effects on the ITG linear instability, nonlinear saturation amplitude, and turbulent transport in the ion root case. However, in the electron root case, the ambipolar electric field significantly reduces the linear ITG growth rate, turbulence intensity, and radial correlation length, resulting in an ion heat conductivity comparable to the neoclassical transport level in the strong shear region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evolution of Λ polarization in the hadronic phase of heavy-ion collisions

Using the AMPT+MUSIC+URQMD hybrid model, we study the global and local spin polarizations of Λ hyperons as functions of the freeze-out temperature of the spin degree of freedom in the hadronic phase of Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV. Including contributions from both the thermal vorticity and thermal shear of the hadronic matter, we find that with the spin freeze-out temperature dropping from the hadronization temperature of 160 MeV to 110 MeV at the kinetic freeze-out, both the global and local spin polarizations of Λ hyperons due to the thermal vorticity decrease by a factor of 2, while those due to the thermal shear decrease quickly and become negligibly small at 140 MeV. Here our results thus suggest that it is important to take into account the evolution of the spin degree of freedom in the hadronic stage of relativistic heavy-ion collisions when theoretically predicted global and local Λ spin polarizations are compared with the experimental measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗