Search NASA⌕ Search

SEARCH · Search NASA

Results for “perturbation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An improved theory of the response of DIII-D H-mode discharges to static resonant magnetic perturbations and its implications for the suppression of edge localized modes

The plasma response to an externally generated, static, n = 2, resonant magnetic perturbation (RMP) in the pedestal region of DIII-D H-mode discharge #158115 is investigated using a toroidal generalization of the asymptotic matching model presented by Fitzpatrick [Phys. Plasmas 27, 042506 (2020)]. Just as in a recent paper [Q. M. Hu et al., Phys. Plasmas 26, 120702 (2019)], it is hypothesized that the density pump-out phenomenon is due to locked magnetic island chains induced at the bottom of the pedestal, whereas the suppression of edge localized modes is associated with a locked magnetic island chain induced at the top of the pedestal. Neutral penetration inside the last closed magnetic flux-surface is found to have a significant influence on locked magnetic island chain formation at the bottom of the pedestal. It is found that locked island formation at the top of the pedestal is only possible when q95 lies in certain narrow windows. Finally, it is inferred that, in order to successfully induce a locked island chain at the top of the pedestal, an external RMP field-coil system must generate a magnetic field that is simultaneously strongly amplified by the plasma (via the ideal kink response) and has a large resonant component at a rational surface that lies close to the top of the pedestal.

Fitzpatrick, R. (ORCID:0000000162379309)↗
DOE OSTI2020

Analytical model for quasi-linear flow response to resonant magnetic perturbation in resistive-inertial and viscous-resistive regimes

In this work, a quasi-linear model for plasma flow response to the resonant magnetic perturbation (RMP) in a tokamak has been rigorously developed in the resistive-inertial (RI) and viscousresistive (VR) regimes purely from the two-field reduced MHD model. Models for plasma response to RMP are commonly composed of equations for the resonant magnetic field response (i.e. the magnetic island) and the torque balance of plasma flow. However, in previous plasma response models, the magnetic island and the torque balance equations are often derived separately from reduced MHD and full MHD equations, respectively. By contrast, in this work we derive both the magnetic island response and the torque balance equations in a quasi-linear model for plasma flow response entirely from a set of two-field reduced MHD equations. Such a quasi-linear model can recover previous plasma flow response models within certain limits and approximations. Furthermore, the physical origins of quasi-linear forces and moments in the flow response equation are also accurately calculated and clarified self-consistently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modeling q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

A toroidal asymptotic matching model of the response of a tokamak plasma to a static resonant magnetic perturbation (RMP) is used to simulate the n = 3 RMP-induced edge-localized-mode-suppression windows in q 95 that are evident when the plasma current is slowly ramped in DIII-D discharge #145380. All quantities employed in the simulation are derived from experimental measurements, apart from the neutral particle data. Three cases are considered. In the first case, the natural frequencies of tearing modes resonant in the plasma are determined by the ion flows at the corresponding resonant surfaces, which is the prediction of nonlinear tearing mode theory. In the second case, the natural frequencies are determined by the local E x B velocities at the resonant surfaces. In the third case, the natural frequencies are determined by the electron flows at the resonant surfaces, which is the prediction of linear tearing mode theory. The second case gives the best agreement between the simulations and the experimental observations. The first and third cases only lead to partial agreement between the simulations and the observations. In the first case, the lack of complete agreement may be a consequence of using an inaccurate assumption for the neutral particle distribution in the pedestal. In the third case, the lack of complete agreement is probably due to the fact that the response of a tokamak plasma to an RMP is not accurately described by linear tearing mode theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Predicting operational windows of ELMs suppression by resonant magnetic perturbations in the DIII-D and KSTAR tokamaks

A newly developed plasma response model, combining the nonlinear two-fluid MHD code TM1 and toroidal MHD code GPEC run in ideal mode, quantitatively predicts the narrow isolated q 95 windows (Δq 95 ~ 0.1) of ELM suppression by n = 1, 2 and 3 resonant magnetic perturbations (RMPs) in both DIII-D and KSTAR tokamaks across a wide range of plasma parameters. The key physics that unites both experimental observations and our simulations is the close alignment of essential resonant q-surfaces and the location of the top of the pedestal prior to an ELM. This alignment permits an applied RMP to produce field penetration due to the lower E × B rotation at the pedestal top rather than being screened. The model successfully predicts that narrow magnetic islands form when resonant field penetration occurs at the top of pedestal, and these islands are easily screened when q 95 moves off resonance, leading to very narrow windows of ELM suppression (typically Δq 95 ~ 0.1). Furthermore, the observed reduction in the pedestal height is also well captured by the calculated classical collisional transport across the island. We recover observed q 95 , β N and plasma shape dependence of ELM suppression due to the effect of magnetic islands on pedestal transport and Peeling-Ballooning- Mode (PBM) stability. Importantly, experiments do occasionally observe wide windows of ELM suppression (Δq 95 > 0.5). Our model reveals that at low pedestal-top density multiple islands open, leading to wide operational windows of ELM suppression consistent with experiment. The model indicates that wide q95 windows of ELM suppression can be achieved at substantially higher pedestal pressure with less confinement degradation in DIII-D by operating at higher toroidal mode number (n = 4) RMPs. As a result, this can have significant implications for the operation of the ITER ELM control coils for maintaining high confinement together with ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the KSTAR tokamak

The EPEC code is employed to model the q95 windows for n = 2 and n = 1 resonant magnetic perturbation (RMP)-induced edge localized mode (ELM) suppression in typical KSTAR H-mode discharges. The plasma equilibria used in the study are derived by rescaling the experimental plasma equilibrium in KSTAR discharge #18594 measured at time t = 6450 ms. The n = 2 and n = 1 ELM suppression windows predicted by EPEC are comparatively narrow in q95 (i.e., Δq 95 ≃ 0.1) and are such that (assuming a relative pedestal pressure reduction of 15% is required to trigger ELM suppression), even in the middle of the windows, there is a threshold RMP coil current of about 1–3 kA/turn that must be exceeded before ELM suppression occurs. Furthermore, the n = 2 and n = 1 ELM suppression windows calculated by the EPEC code are consistent with the experimental observations and are also broadly similar to those recently calculated by the TM1 code [Hu et al., Phys. Plasmas 28, 052505 (2021)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Numerical study of $\delta$-function current sheets arising from resonant magnetic perturbations

We report general three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially regions of stochastic field lines when islands overlap. Numerically resolving singular current sheets in the ideal magnetohydrodynamics (MHD) limit has been a significant challenge. This work presents numerical solutions of the Hahm–Kulsrud–Taylor (HKT) problem, which is a prototype for resonant singular current sheet formation. The HKT problem is solved by two codes: a Grad–Shafranov (GS) solver and the Stepped Pressure Equilibrium Code (SPEC) code. The GS solver has built-in nested flux surfaces with prescribed magnetic fluxes. The SPEC code implements multi-region relaxed magnetohydrodynamics (MRxMHD), whereby the solution relaxes to a Taylor state in each region while maintaining force balance across the interfaces between regions. As the number of regions increases, the MRxMHD solution appears to approach the ideal MHD solution assuming a continuum of nested flux surfaces. We demonstrate agreement between the numerical solutions obtained from the two codes through a convergence study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A resistive MHD model and simulation on plasma flow evolution in the presence of resonant magnetic perturbation in a tokamak

Nonaxisymmetric magnetic fields such as the intrinsic error field and the externally applied resonant magnetic perturbation (RMP) in a tokamak are known to influence the plasma momentum transport and flow evolution through plasma response, which itself strongly depends on the plasma flow as well. The nonlinear interaction between plasma response and flow has been previously modeled in the conventional error field theory with the “no-slip” condition, which has been recently extended to allow the “free-slip” condition. In this work, we further target this specific process and numerically simulate the nonlinear plasma response and flow evolution in the presence of a single-helicity RMP in a circular-shaped model tokamak configuration, based on the full resistive MHD model in the initial-value code NIMROD. Time evolution of the parallel (to k) flow or “slip frequency” profile and its asymptotic steady state obtained from the NIMROD simulations are compared with both conventional and extended nonlinear response theories. Here, k is the wave vector of the propagating island. Good agreement with the extended theory with free-slip condition has been achieved for the parallel flow profile evolution in response to RMP in all resistive regimes, whereas the difference from the conventional theory with the no-slip condition tends to diminish as the plasma resistivity approaches zero.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Computational study of runaway electrons in MST tokamak discharges with applied resonant magnetic perturbation

A numerical study of magnetohydrodynamics (MHD) and tracer-particle evolution investigates the effects of resonant magnetic perturbations (RMPs) on the confinement of runaway electrons (REs) in tokamak discharges conducted in the Madison Symmetric Torus. In computational results of applying RMPs having a broad toroidal spectrum but a single poloidal harmonic, m = 1 RMP does not suppress REs, whereas m = 3 RMP achieves significant deconfinement but not the complete suppression obtained in the experiment. MHD simulations with the NIMROD code produce sawtooth oscillations, and the associated magnetic reconnection can affect the trajectory of REs starting in the core region. Simulations with m = 3 RMP produce chaotic magnetic topology over the outer region, but the m = 1 RMP produces negligible changes in field topology, relative to applying no RMP. Using snapshots of the MHD simulation fields, full-orbit relativistic electron test particle computations with KORC show ≈50% loss from the m = 3 RMP compared to the 10%–15% loss from the m = 1 RMP. Here, test particle computations of the m = 3 RMP in the time-evolving MHD simulation fields show correlation between MHD activity and late-time particle losses, but total electron confinement is similar to computations using magnetic-field snapshots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Microgap breakdown with floating metal rod perturbations

Here, we report the characterization of microgap breakdown with perturbations from a metal rod floating between anode and cathode electrodes. The effects of the metal rod on the electric field distribution and the field enhancement factor are evaluated by numerical simulation and the conformal mapping method, and they indicate that the field emission regime is not reached. The breakdown voltages in the Townsend discharge regime are determined based on the voltage–current characteristics, which are obtained from two-dimensional fluid simulations. It is found that the breakdown characteristics can be significantly modulated by the floating metal rod, and the breakdown curve (breakdown voltage vs the net gap distance) is no longer U-shaped, which deviates from the conventional Paschen's law. The underlying physical mechanisms are related to the electric field enhancement, curved breakdown path, and nonuniform ion flux caused by the electric shielding effect. The results provide insights into breakdown characteristics in microscale discharges, which may promote conventional investigation of simplified clean gaps toward more complex conditions (e.g., with floating microparticles) in miniaturized plasma devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Stochastic fluctuation and transport of tokamak edge plasmas with the resonant magnetic perturbation field

Here we present that a statistical method known as the complexity–entropy analysis is useful to characterize a state of plasma turbulence and flux in the resonant magnetic perturbation (RMP) edge localized mode (ELM) control experiment. The stochastic pedestal top temperature fluctuation in the RMP ELM suppression phase is distinguished from the chaotic fluctuation in the natural ELM-free phase. It is discussed that the stochastic temperature fluctuation can be originated from the narrow layer of the field penetration on the pedestal top. The forced magnetic island can emit the resonant drift wave of comparable sizes (relatively low-k) in the RMP ELM suppression phase, and it can result in the generation of stochastic higher wavenumber fluctuations coupled to tangled fields around the island. The analysis of the ion saturation current measurement around the major outer striking point on the divertor shows that it also becomes more stochastic as the stronger plasma response to the RMP field is expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of aspect ratio on plasma response to resonant magnetic perturbations in tokamak devices

A systematic numerical study is carried out, computing and comparing the plasma response to the resonant magnetic perturbation (RMP) field, applied for controlling edge localized modes (ELMs), in a series of tokamak plasmas with varying aspect ratio and utilizing the MARS-F code. The aspect ratio is scanned either by varying the plasma major radius at a fixed minor radius or by varying the latter while fixing the former. Both approaches yield similar results when compared in terms of quantities with proper normalizations. In general, a non-monotonic dependence of the resonant response field (normalized by the vacuum counterpart) near the plasma edge is found with varying aspect ratio, indicating that a given ELM control coil current configuration strongly favors plasmas with a certain aspect ratio. This optimal aspect ratio, on the other hand, depends on the toroidal as well as poloidal (i.e., coil phasing) spectra of the applied RMP field. The equilibrium (edge) safety factor, the plasma shape, and the plasma toroidal flow are all fixed to ensure that the effects identified here are predominantly due to the plasma aspect ratio.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Porting fragmentation methods to GPUs using an OpenMP API: Offloading the resolution-of-the-identity second-order Møller–Plesset perturbation method

Here, using an OpenMP Application Programming Interface, the resolution-of-the-identity second-order Møller–Plesset perturbation (RI-MP2) method has been off-loaded onto graphical processing units (GPUs), both as a standalone method in the GAMESS electronic structure program and as an electron correlation energy component in the effective fragment molecular orbital (EFMO) framework. First, a new scheme has been proposed to maximize data digestion on GPUs that subsequently linearizes data transfer from central processing units (CPUs) to GPUs. Second, the GAMESS Fortran code has been interfaced with GPU numerical libraries (e.g., NVIDIA cuBLAS and cuSOLVER) for efficient matrix operations (e.g., matrix multiplication, matrix decomposition, and matrix inversion). The standalone GPU RI-MP2 code shows an increasing speedup of up to 7.5× using one NVIDIA V100 GPU with one IBM 42-core P9 CPU for calculations on fullerenes of increasing size from 40 to 260 carbon atoms using the 6-31G(d)/cc-pVDZ-RI basis sets. A single Summit node with six V100s can compute the RI-MP2 correlation energy of a cluster of 175 water molecules using the correlation consistent basis sets cc-pVDZ/cc-pVDZ-RI containing 4375 atomic orbitals and 14 700 auxiliary basis functions in ~0.85 h. In the EFMO framework, the GPU RI-MP2 component shows near linear scaling for a large number of V100s when computing the energy of an 1800-atom mesoporous silica nanoparticle in a bath of 4000 water molecules. The parallel efficiencies of the GPU RI-MP2 component with 2304 and 4608 V100s are 98.0% and 96.1%, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Non-perturbative many-body treatment of molecular magnets

Molecular magnets have received significant attention because of their potential applications in quantum information and quantum computing. A delicate balance of electron correlation, spin–orbit coupling (SOC), ligand field splitting, and other effects produces a persistent magnetic moment within each molecular magnet unit. The discovery and design of molecular magnets with improved functionalities would be greatly aided by accurate computations. However, the competition among the different effects poses a challenge for theoretical treatments. Electron correlation plays a central role since d- or f-element ions, which provide the magnetic states in molecular magnets, often require explicit many-body treatments. SOC, which expands the dimensionality of the Hilbert space, can also lead to non-perturbative effects in the presence of strong interaction. Furthermore, molecular magnets are large, with tens of atoms in even the smallest systems. We show how an ab initio treatment of molecular magnets can be achieved with auxiliary-field quantum Monte Carlo, in which electron correlation, SOC, and material specificity are included accurately and on an equal footing. Furthermore, the approach is demonstrated by an application to compute the zero-field splitting of a locally linear Co 2+ complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Integrated control of edge localized modes and divertor flux using mixed toroidal harmonic resonant magnetic perturbations in EAST

Mixed harmonic resonant magnetic perturbations (RMPs) for integrated edge localized modes (ELMs) and divertor flux control are demonstrated in EAST target plasmas of low input torque and normalized beta β N ~1.7–1.9, which are close to the equivalent value in ITER high Q operation. The applied RMPs are designed to combine a static harmonic of the toroidal mode number n = 3 with a static or rotating harmonic of n = 2. ELM suppression is achieved without a drop of plasma energy confinement, and tungsten concentration is effectively reduced during the application of RMPs. With mixed harmonics, the toroidal varying steady state heat and particle fluxes on the divertor target can be modified with the rotating n = 2 harmonic, which agrees with the numerical modeling of three-dimensional magnetic topology, with plasma responses being taken into account. ELM suppression correlates with the times of larger n = 3 response with mixed n = 2 and n = 3 RMPs. The mixture of harmonics and the rotating n = 2 harmonic does not require additional coil current because the variation is only in the upper-lower coil current phase space. Furthermore, these results further affirm the effectiveness of integrated ELM and divertor flux control using RMPs with mixed harmonics and improve the understanding of the role of plasma responses in ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Neoclassical transport of impurities in tokamaks with non-axisymmetric perturbations

The effect of resonant magnetic perturbations (RMPs) on neoclassical transport of impurities is calculated in DIII-D and NSTX. Here, the neoclassical fluxes are evaluated using the NEO code with nonlinear one-fluid nonaxisymmetric equilibrium calculated using M3D-C1. Neoclassical fluxes of impurities show significant changes with RMPs if the impurity resides in low-collisionality regime, but are weakly affected by RMPs in the Pfirsch–Schlüter (P–S) regime. Charge number (Z) of the impurity affects the collisionality of impurity species, which then determines the effect of 3D fields on neoclassical impurity transport. This suggests that RMPs can be possibly used for exhaust of low-Z impurities in these facilities, but have little effect for impurities with charge number greater than 10 or so. Additionally, it is shown that it is the change in the convective fluxes which is a main contributor in modifying the neoclassical impurity fluxes in the presence of RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Ultrafast photochemistry and electron-diffraction spectra in n → (3s) Rydberg excited cyclobutanone resolved at the multireference perturbative level

We study the ultrafast time evolution of cyclobutanone excited to the singlet n → Rydberg state through non-adiabatic surface-hopping simulationsperformed at extended multi-state complete active space second-order perturbation (XMS-CASPT2) level of theory. These dynamics predict relaxation to the ground-state with a timescale of 822 ± 45 fs with minimal involvement of the triplets. The major relaxation path to the ground-state involves a three-state degeneracy region and leads to a variety of fragmented photoproducts. We simulate the resulting time-resolved electron-diffraction spectra, which track the relaxation of the excited state and the formation of various photoproducts in the ground state.

Chemistry↗
DOE OSTI2025

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Identifying sequence perturbations to an intrinsically disordered protein that determine its phase-separation behavior

Phase separation of intrinsically disordered proteins (IDPs) commonly underlies the formation of membraneless organelles, which compartmentalize molecules intracellularly in the absence of a lipid membrane. Identifying the protein sequence features responsible for IDP phase separation is critical for understanding physiological roles and pathological consequences of biomolecular condensation, as well as for harnessing phase separation for applications in bioinspired materials design. To expand our knowledge of sequence determinants of IDP phase separation, we characterized variants of the intrinsically disordered RGG domain from LAF-1, a model protein involved in phase separation and a key component of P granules. Based on a predictive coarse-grained IDP model, we identified a region of the RGG domain that has high contact probability and is highly conserved between species; deletion of this region significantly disrupts phase separation in vitro and in vivo. We determined the effects of charge patterning on phase behavior through sequence shuffling. We designed sequences with significantly increased phase separation propensity by shuffling the wild-type sequence, which contains well-mixed charged residues, to increase charge segregation. This result indicates the natural sequence is under negative selection to moderate this mode of interaction. We measured the contributions of tyrosine and arginine residues to phase separation experimentally through mutagenesis studies and computationally through direct interrogation of different modes of interaction using all-atom simulations. Finally, we show that despite these sequence perturbations, the RGG-derived condensates remain liquid-like. Together, these studies advance our fundamental understanding of key biophysical principles and sequence features important to phase separation.

36 MATERIALS SCIENCE↗