Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability theory”

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 613 records · Page 34

Anomalous entropy-driven kinetics of dislocation nucleation

The kinetics of dislocation reactions, such as dislocation multiplication, controls the plastic deformation in crystals beyond their elastic limit, therefore critical mechanisms in a number of applications in materials science. We present a series of large-scale molecular dynamics simulations that shows that one such type of reactions, the nucleation of dislocation at free surfaces, exhibit unconventional kinetics, including unexpectedly large nucleation rates under compression, very strong entropic stabilization under tension, as well as strong non-Arrhenius behavior. These unusual kinetics are quantitatively rationalized using a variational transition state theory approach coupled with an efficient numerical scheme for the estimation of vibrational entropy changes. These results highlight the need for a variational treatment of the kinetics to quantitatively capture dislocation reaction kinetics, especially at low-to-moderate strains where large deformations are required to activate reactions. These observations suggest possible explanations to previously observed unconventional deformation kinetics in both molecular dynamics simulations and experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Oxidative dehydrogenation of propane on the oxygen adsorbed edges of boron nitride nanoribbons

Metal-free boron nitride has recently been reported to exhibit high catalytic activity and selectivity towards oxidative dehydrogenation of propane (ODHP). Despite several experimental studies, the exact nature and function of the active sites in this emerging catalyst are unknown. In the present work, density functional theory calculations are combined with microkinetic modeling to systematically explore the ODHP on oxygen passivated boron nitride nanoribbons (BNNRs). Here, the relative stabilities of different edge structures of BNNR were examined using ab initio thermodynamic analysis, and the most stable oxygen passivated edge structure was used for a mechanistic study. Microkinetic analysis revealed that a N 2 O or NO x -type active site is active and selective for the conversion of propane to propene. In addition to a heterogeneous catalytic cycle, the proposed N 2 O/NO x -type active sites are expected to be able to generate gas phase C 3 H 7 ˙, which in turn can trigger gas phase reactions of ODHP as experimentally speculated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Novel intrinsic helical cores and MHD dynamo flux pumping evidence in DIII-D

Evidence of a MHD dynamo brought about by intrinsic helical core structures in negative triangularity (NT) and ITER baseline scenario (IBS) plasmas demonstrates paths to operating tokamaks without seeds for deleterious global MHD modes. A stationary sawtooth- and edge localized mode-free state maintains high performance for long durations in a strong NT shape in the DIII-D tokamak. Disparity between minimum safety factor values predicted by magnetic flux diffusion theory and derived from experiment suggests anomalous flux diffusion by an MHD dynamo, driven by flows associated with quasi-interchange modes. Stability calculations with the GATO code find NT and IBS experimental equilibria unstable to quasi-interchange modes, in agreement with observed saturated activity. IBS discharges with similar magnetic winding structure occupy robust helical core states following a transition from sawtooth instabilities. These states agree with prior modeling of helical core bifurcation thresholds with the VMEC equilibrium code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantified limits of the nuclear landscape

The chart of the nuclides is limited by particle drip lines beyond which nuclear stability to proton or neutron emission is lost. Predicting the range of particle-bound isotopes poses an appreciable challenge for nuclear theory as it involves extreme extrapolations of nuclear masses well beyond the regions where experimental information is available. Still, quantified extrapolations are crucial for a wide variety of applications, including the modeling of stellar nucleosynthesis. We use microscopic nuclear global mass models, current mass data, and Bayesian methodology to provide quantified predictions of proton and neutron separation energies as well as Bayesian probabilities of existence throughout the nuclear landscape all the way to the particle drip lines. Here, we apply nuclear density functional theory with several energy density functionals. We also consider two global mass models often used in astrophysical nucleosynthesis simulations. To account for uncertainties, Bayesian Gaussian processes are trained on the separation-energy residuals for each individual model, and the resulting predictions are combined via Bayesian model averaging. This framework allows to account for systematic and statistical uncertainties and propagate them to extrapolative predictions. We establish and characterize the drip-line regions where the probability that the nucleus is particle- bound decreases from 1 to 0. In these regions, we provide quantified predictions for one- and two-nucleon separation energies. According to our Bayesian model averaging analysis, 7759 nuclei with Z ≤ 119 have a probability of existence ≥ 0.5. The extrapolation results obtained in this study will be put through stringent tests when new experimental information on existence and masses of exotic nuclei becomes available. In this respect, the quantified landscape of nuclear existence obtained in this study should be viewed as a dynamical prediction that will be fine-tuned when new experimental information and improved global mass models become available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Upward band gap bowing and negative mixing enthalpy in multi-component cubic halide perovskite alloys

Physical properties intermediate between constituents of alloys can be achieved as downward convex positive bowing, upward concave negative bowing, or zero bowing. Such bowing effects are essential for band gap engineering in semiconductor alloys. Upward band gap bowing effects are rather rare, hindering the exploration on half of the available physical property space of alloys. Part of this being a rare event is related to the need to stabilize an alloy with low mixing enthalpy, so it does not phase separate. Here, in this paper, we find via density functional theory that one can satisfy the simultaneous conditions of negative mixing enthalpy and upward band gap bowing in four-component AB⁢X 3 halide perovskite alloys in the cubic perovskite structure. Such perovskite alloys have the B-site occupied by a mixture of group IVB and IIB elements that have the IVB-𝑠 and IIB-s states in the valence bands and conduction bands, respectively, leaving the delocalized 𝑠 states to strongly repel each other. This 𝑠−𝑠 repulsion leads to the upward band gap bowing and negative mixing enthalpies simultaneously. Remarkably, we identify a perovskite alloy that has a band gap much larger than all its components. Analogous trends of upward band gap bowing and negative mixing enthalpy also appear in the corresponding three-component and two-component AB⁢X 3 halide perovskite alloys. These observations of upward band gap bowing and negative mixing enthalpy will significantly accelerate the design of stable upward band gap bowing alloys in a broad range of material families.

14 SOLAR ENERGY↗

Development of an Ultrahigh-bandwidth Phase Contrast Imaging System for detection of electron scale turbulence and Gigahertz Radio-Frequency Waves

The study of waves and turbulence is vital to the development of future reactor-grade plasma devices developed in the quest for fusion energy. These fluctuations are responsible for moving heat and particles across the magnetic field, and a predictive understanding of them is needed to achieve the density and temperature required to sustain a plasma fusion reaction. While many techniques have been developed for measuring waves and fluctuations, every measurement method has limitations. There are relatively few techniques for measuring very high frequency fluctuations, such as radio frequency waves injected to heat the plasma, unstable waves driven by suprathermal particles, or short wavelength electrostatic waves driven by electron temperature or density gradients. The present project builds upon the proven phase contrast imaging (PCI) technique to extend the response of the diagnostic by orders of magnitude in frequency and almost a factor of ten in spatial resolution. PCI provides a measurement of electron density based on small angle scattering of a CO-2 laser beam by using optical techniques to render a phase shift as an intensity change on a detector. Due to the telecommunications revolution, technological development by manufacturers has focused on components in the near -infrared, so that high- quality lasers and detectors at 1.55 µm are readily available. Shifting PCI design to a new, shorter wavelength has numerous advantages and challenges. Similar detector performance is available with room-temperature arrays with GHz bandwidth, while the detector arrays for 10.6 µm required liquid nitrogen cooling and were therefore limited to a bandwidth of about 1 MHz. Shifting to a shorter wavelength reduces the angle at which the laser beam scatters off of plasma waves, allowing more such scattered components to pass through the aperture of the vacuum vessel port and be collected by the PCI, which increases the spatial resolution. Concomitant with these benefits, various questions of performance arise. At shorter laser wave- length, the sensitivity to mirror and lens quality is increased, the contribution of the laser to the overall system noise is increased, and the sensitivity to vibrations is increased. The custom optical components at the heart of the PCI technique were required to be properly scaled for the shorter wavelength, so fabrication technologies needed to be explored. This project was designed to show that a low noise, high response PCI system at 1.55 µm could be constructed and operated, and then to quantify potential issues to allow extrapolation to a full-size production system providing physics measurements on a large plasma device. The first stage, producing the custom optical component called a Phase Plate, was successfully achieved using two methods. First, an easily reproduced masking and coating technique was able to produce good phase plates with the required parameters. Second, a nanofabrication technique was found to produce extremely high quality phase plates at a competitive price. The PCI constructed with the new phase plates and 1.55 µm laser was found to provide excellent wavelength measurements with the theoretically expected response. The sensitivity to optical surface quality was found to be in line with previous measurements at 10.6 µm. The observed signal-to-noise ratio was similar to the theoretically expected value. The effect of vibrations on PCI was studied with the first measurement of the effect of beam motion on PCI response and comparison to theory, allowing for a quantitative prediction of the effect of vibrations on a production PCI system and the requirements for improved beam stabilization. PCI is an extremely cost-effective method to provide a low noise, absolutely calibrated measurement of plasma fluctuations across a wide spatial scale. This project has shown that a 1.55 µm PCI using modern techniques and components is less expensive than the 10.6 µm PCI of a few years ago, with the largest savings in phase plate fabrication and the infrared detector array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine learning prediction of the mechanical properties of refractory multicomponent alloys based on a dataset of phase and first principles simulation

In this work, a dataset including structural and mechanical properties of refractory multicomponent alloys was developed by fusing computations of phase diagram (CALPHAD) and density functional theory (DFT). The refractory multicomponent alloys, also named refractory complex concentrated alloys (CCAs) which contain 2–5 types of refractory elements were constructed based on Special Quasi-random Structure (SQS). The phase of alloys was predicted using CALPHAD and the mechanical property of alloys with stable and single body-centered cubic (BCC) at high temperature (over 1,500°C) was investigated using DFT-based simulation. As a result, a dataset with 393 refractory alloys and 12 features, including volume, melting temperature, density, energy, elastic constants, mechanical moduli, and hardness, were produced. To test the capability of the dataset on supporting machine learning (ML) study to investigate the property of CCAs, CALPHAD, and DFT calculations were compared with principal components analysis (PCA) technique and rule of mixture (ROM), respectively. It is demonstrated that the CALPHAD and DFT results are more in line with experimental observations for the alloy phase, structural and mechanical properties. Furthermore, the data were utilized to train a verity of ML models to predict the performance of certain CCAs with advanced mechanical properties, highlighting the usefulness of the dataset for ML technique on CCA property prediction.

36 MATERIALS SCIENCE↗

Ab Initio Molecular Dynamics Insights into Stress Corrosion Cracking and Dissolution of Metal Oxides

Oxide phases such as α-Fe2O3 (hematite) and α-Al2O3 (corundum) are highly insoluble in water; however, subcritical crack growth has been observed in humidity nonetheless. Chemically induced bond breaking at the crack tip appears unlikely due to sterically hindered molecular transport. The molecular mechanics of a crack in corundum with a reactive force field reveal minimal lattice trapping, leading to bond breaking before sufficient space opens for water transport. To address this, we model a pre-built blunt crack with space for H2O molecule adsorption at the tip and show that it reduces fracture toughness by lowering the critical J-integral. Then, we explore stress-enhanced dissolution to understand the mechanism of crack tip blunting in the oxide/water system. Density functional theory combined with metadynamics was employed to describe atomic dissolution from flat hematite and corundum surfaces in pure water. Strain accelerates dissolution, stabilizing intermediate states with broken bonds before full atom detachment, while the free energy profile of unstrained surfaces is almost monotonic. The atomistic calculations provided input for a kinetic model, predicting the shape evolution of a blunt crack tip, which displays three distinct regimes: (i) dissolution primarily away from the tip, (ii) enhanced blunting near but not at the apex, and (iii) sharpening near the apex. The transition between regimes occurs at a low strain, highlighting the critical role of water in the subcritical crack growth of oxide scales, with dissolution as the fundamental microscopic mechanism behind this process.

Chemistry↗

Phase and Morphology of Water-ice Grains Formed in a Cryogenic Laboratory Plasma

Abstract Grains of ice are formed spontaneously when water vapor is injected into a weakly ionized laboratory plasma in which the background gas has been cooled to cryogenic temperatures comparable to those of deep space. These ice grains are levitated indefinitely within the plasma so that their time evolution can be observed under free-floating conditions. Using microscope imaging, ice grains are shown to have a spindle-like fractal structure and grow over time. Both crystalline and amorphous phases of ice are observed using Fourier transform infrared spectroscopy. A mix of crystalline and amorphous grains coexists under certain thermal conditions, and a linear mixing model is used on the ice absorption band surrounding 3.2 μ m to examine the ice phase composition and its temporal stability. The extinction spectrum is also affected by inelastic scattering as grains grow, and characteristic grain radii are obtained from Mie scattering theory and compared to size measurements from direct imaging. Observations are used to compare possible ice nucleation mechanisms, and it is concluded that nucleation is likely catalyzed by ions, as ice does not nucleate in the absence of plasma and impurities are not detected. Ice grain properties and infrared extinction spectra show similarity to observations of some astrophysical ices observed in protoplanetary disks, implying that the fractal morphology of the ice and observed processes of homogeneous ice nucleation could occur as well in such astrophysical environments with weakly ionized conditions.

Astronomy & Astrophysics↗

Empty perovskites as Coulomb floppy networks: Entropic elasticity and negative thermal expansion

Floppy networks (FNs) provide valuable insight into the origin of anomalous mechanical and thermal properties in soft matter systems, from polymers, rubber, and biomolecules to glasses and granular materials. Here, we use the same FN concept to construct a quantitative microscopic theory of empty perovskites, a family of crystals with ReO 3 structure, which exhibit a number of unusual properties. One remarkable example is ScF 3 , which shows a near-zero-temperature structural instability and large negative thermal expansion (NTE). Furthermore, we trace these effects to an FN-like crystalline architecture formed by strong nearest-neighbor bonds, which is stabilized by net electrostatic repulsion that plays a role similar to osmotic pressure in polymeric gels. NTE in these crystalline solids, which we conceptualize as Coulomb floppy networks, emerges from the tension effect of Coulomb repulsion combined with the FN's entropic elasticity and has the same physical origin as in gels and rubber. Our theory provides an accurate, quantitative description of phonons, thermal expansion, compressibility, and structural phase diagram, all in excellent agreement with experiments. The entropic stabilization of critical soft modes, which play only a secondary role in NTE, explains the observed phase diagram. Significant entropic elasticity resolves the puzzle of a marked, ≈50% discrepancy between the experimentally observed bulk modulus and ab initio calculations. The Coulomb FN approach is potentially applicable to other important materials with markedly covalent bonds, from perovskite oxides to iron chalcogenides, whose anomalous vibrational and structural properties are still poorly understood.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effects of carbon concentration on the local atomic structure of amorphous GST

Ge-Sb-Te (GST) alloys are leading phase-change materials for data storage due to the fast phase transition between amorphous and crystalline states. Ongoing research aims at improving the stability of the amorphous phase to improve retention. This can be accomplished by the introduction of carbon as a dopant to Ge 2 Sb 2 Te 5 , which is known to alter the short- and mid-range structure of the amorphous phase and form covalently bonded C clusters, both of which hinder crystallization. The relative importance of these processes as a function of C concentration is not known. Here we used molecular dynamics simulation based on density functional theory to study how carbon doping affects the atomic structure of GST-C. Carbon doping results in an increase in tetrahedral coordination, especially of Ge atoms, and this is known to stabilize the amorphous phase. We observe an unexpected, non-monotonous trend in the number of tetrahedral bonded Ge with the amount of carbon doping. Our simulations show an increase in the number of tetrahedral bonded Ge up to 5 at.% C, after which the number saturates and begins to decrease above 14 at.% C. The carbon atoms aggregate into clusters, mostly in the form of chains and graphene flakes, leaving less carbon to disrupt the GST matrix at higher carbon concentrations. Different degrees of carbon clustering can explain divergent experimental results for recrystallization temperature for carbon doped GST.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Disruption avoidance via radio frequency current condensation in magnetic islands produced by off-normal events

This article discusses the use of radio frequency (RF) current drive to stabilize large islands, focusing on nonlinear effects that appear when relatively high powers are used to stabilize large islands. We are interested in developing a capability to stabilize large islands via RF driven currents to avoid the need for mitigation to the extent possible. As tokamaks are designed and built with increasing levels of stored energy in the plasma, disruptions become increasingly dangerous. It has been reported that 95% of the disruptions in the Joint European Torus tokamak with the ITER-like wall are preceded by the growth of large locked islands. These large islands are mostly produced by off-normal events other than neoclassical tearing modes. This paper presents theory and modeling for a nonlinear “RF current condensation” effect that can concentrate the RF driven current near the center of a large island, thereby increasing the efficiency of the stabilization. A nonlinear shadowing effect can hinder the stabilization of islands if the aiming of the ray trajectories does not properly consider the nonlinear effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Establishing defect-property relationships for 2D-nanomaterials (Final technical report)

Studies of new families of two-dimensional nanomaterials (2DNMs) have established that they possess unique properties that diverge from those of their bulk counterparts. This project aimed to determine how tolerant 2DNMs are to extreme photon and particle fluxes and to identify the mechanisms governing their radiation response. Early work indicated that graphene is not as representative of other 2DNMs as previously assumed; however, the origin of this difference remained unclear. To address these knowledge gaps, this project investigated the physical processes occurring at multiple length scales in transition metal dichalcogenides (TMDs) and quantified their structural stability and property evolution under far-from-equilibrium conditions. In-situ and ex-situ ion and electron irradiations were performed to directly control and monitor defect formation in TMDs. Complementary density functional theory (DFT) and molecular dynamics (MD) simulations were employed to elucidate the mechanisms of defect generation and evolution. The outcomes of this work established mechanistic understanding of irradiation-induced defect formation in 2DNMs, quantified the radiation tolerance of TMDs, and built a fundamental knowledge base for correlating defect structures with material properties. Collectively, these results provide new insights into the stability of low-dimensional materials in extreme environments and enable the predictive design of radiation-tolerant 2DNMs.

2D materials↗

Inducing Fe moment in LaFeSi with p-block element substitution

We studied the LaFeSi phase using density functional theory based full-potential linearized augmented plane wave (FP-LAPW) method. Specifically we examined the effect of p-element substitution on the stability and magnetic properties focusing on pathways to induce Fe magnetic moments. We demonstrate that either partial or complete substitutions at the Si 2 c-site by several p-block elements lead to non-zero Fe moments. Our theoretical study shows that partially substituted LaFeSi 1-y Ga y and LaFeSi 1-y Al y are more thermodynamically stable than the fully substituted LaFeGa and LaFeAl. We also found that the P substitutions, either partial or complete, have the most negative formation energies, however the compounds containing phosphorus are non-magnetic. Our work highlights a pathway toward manipulations of the otherwise quenched Fe-moments in LaFeSi-related compounds and help in advancing the potential magnetic functionalities of LaFeX compounds.

36 MATERIALS SCIENCE↗

High Performance SOFCs with a Superior Stability for Reliable and Durable Power Systems

Next generation of fuel cells, electrolyzers, and batteries requires higher power, faster kinetics, and larger energy density, which necessitate the use of compositionally complex oxides to achieve multifunctionalities and activity. These compositionally complex oxides may change their phases and structures during an electrochemical process – a so-called “electrochemically driven phase transformation”. The origin for such a phase change has remained obscure. More importantly, there is a need to develop high performance solid oxide fuel cells with an enhanced stability. In this work, the La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3 (LSCF) cathode surface is modified by infiltration of Pr 6 O 11 and the power density at 0.8V and 750 °C is improved by 21%. In addition, by replacing the traditional barrier layer Gd 0.2 Ce 0.8 O 1.9 with mixed conducting Pr 0.1 Gd 0.1 Ce 0.8 O 1.9 , the power density increases by 38%. The different mechanism of promotions was investigated by electrochemical impedance spectroscopy. The ohmic resistance is dramatically reduced by applying the PGCO interlayer, and the distribution of relaxation time was used to analyze the mechanism for which the polarization resistance was decreased attributing to the mixed conduction nature in PrO x . An increase of power density at 0.8 V of 0.358 W/cm 2 (71%) is achieved with the implementation of both surface modification and buffer layer engineering. An experimental study and a theoretical analysis were then carried out on phase evolution in praseodymium nickelates. Nickelate-based electrodes show up to 60× greater phase transformation during operation when compared to thermally annealed ones. Theoretical analysis suggests that the presence of a reduced oxygen partial pressure at the interface between the oxygen electrode and the electrolyte is the origin for the phase change in an oxygen electrode. Guided by the theory, an addition of the electronic conduction in the interface layer leads to the significant suppression of phase change, while improving cell performance and performance stability. When an oxygen electrode is under polarization, the oxygen partial pressure at the interface between the oxygen electrode and the electrolyte is lower than that of incoming oxidant. Under a high polarization, the environment at the aforementioned interface may lead to phase transformation of the oxygen electrode. The local oxygen partial pressure is determined by the transport properties at the interfaces. An addition of the electronic conduction in the interface layer, for instance using (Pr,Gd)-doped ceria to replace Gd-doped ceria, results in improved cell performance and performance stability, while the phase transformation is significantly suppressed. This work provides a fundamental understanding of the origin for phase transformation in oxygen electrodes during operation and use this knowledge to develop a high-performance electrode that exhibits improved performance stability.

30 DIRECT ENERGY CONVERSION↗

Stabilization of Dinuclear Rhodium and Iridium Clusters on Layered Titanate and Niobate Supports

Atomically dispersed organometallic clusters can provide well-defined nuclearity of active sites for both fundamental studies as well as new regimes of activity and selectivity in chemical transformations. More recently, dinuclear clusters adsorbed onto solid surfaces have shown novel catalytic properties resulting from the synergistic effect of two metal centers to anchor different reactant species. Difficulty in synthesizing, stabilizing, and characterizing isolated atoms and clusters without agglomeration challenges allocating catalytic performance to atomic structure. Here, we explore the stability of dinuclear rhodium and iridium clusters adsorbed onto layered titanate and niobate supports using molecular precursors. Both systems maintain their nuclearity when Statistical analysis of HAADF-STEM images revealed that rhodium and iridium dimers had mean cluster-to-cluster distances very similar to what is expected from a random distribution of atoms over a large area, indicating that they are dispersed without aggregation. The stability of dinuclear rhodium clusters supported on titanate nanosheets was also investigated by X-ray absorption fine structure (EXAFS), DRIFTS, and first-principles calculations. Both X-ray absorption spectroscopy and HAADF-STEM simulations, guided by density functional theory (DFT)-optimized structure models, suggested that rhodium dimers adsorb onto the nanosheets in an end-on binding mode that is stable up to 100 degrees C under reducing conditions. Finally, this study highlights that crystalline nanosheets derived from layered metal oxides can be used as model supports to selectively stabilize dinuclear clusters, which could have implications for heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical Insights on the Fragmentation of Cellulosic Oligomers to Form Hydroxyacetone and Hydroxyacetaldehyde

The presence of heavy unknown oligomeric sugar products in bio-oil is evidenced in experimental results reported in the literature. In this paper, we study the fragmentation reactions yielding acetol and glycolaldehyde from oligomeric sugars following previous work on dehydration reactions to propose structures of these oligomers. Acetol and glycolaldehyde are primary products of cellulose fast pyrolysis but the fragmentation reaction mechanism of these compounds from oligomers merits further study. The density functional theory (DFT) approach was employed to study this reaction. Results revealed that acetol and glycolaldehyde fragments are favorably removed from the non-reducing end based on their thermodynamic stabilities. Theoretical FTIR and NMR spectra were calculated to aid in understanding the structures of the oligomeric sugars. Also, the thermodynamics and physical properties of these compounds were estimated using the Group Contribution Method (GCM). In conclusion, these properties are essential in the design of processing technologies and the upgrading of products.

10 SYNTHETIC FUELS↗

Hybrid simulations of sub-cyclotron compressional and global Alfven 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 (CAE) and global (GAE) Alfven Eigenmodes 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 power in NSTX experiments, prompting further examination of the conditions for their excitation. Linear simulations are performed with the hybrid MHD-kinetic initial value code HYM in order to capture the general Doppler-shifted cyclotron resonance that drives the modes. Three distinct types of modes are found in simulations -- co-CAEs, cntr-GAEs, and co-GAEs -- with differing spectral and stability properties. The simulations reveal that unstable GAEs are more ubiquitous than unstable CAEs, 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 simulations and estimated analytically for relevant sources not present in the simulation model, indicating that co-CAEs are closer to marginal stability than modes driven by the cyclotron resonances.

Compressional Alfven Eigenmodes↗