Search NASA⌕ Search

SEARCH · Search NASA

Results for “global instability”

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 OSTI2023

The r -modes of slowly rotating, stratified neutron stars

ABSTRACT The only r-modes that exist in a globally barotropic, rotating, Newtonian star are the fundamental l = |m| solutions, where l and m are the indices of the spherical harmonic $Y_l^m$ that describe the mode’s angular dependence. This is in stark contrast to a stellar model that is non-barotropic throughout its interior, which hosts all the l ≥ |m| perturbations including radial overtones. In reality, neutron stars are stratified with locally barotropic regions. Therefore, we explore how stratification alters a star’s ability to support r-modes. We consider the globally stratified case and examine the behaviour of the modes as the star gets close to barotropicity. In this limit, we find that all but the fundamental l = |m| perturbations change character and become generic inertial modes. Restricting the analysis to l = |m| perturbations, we develop the r-mode equations in order to consider stellar models that exhibit local barotropicity. Our results for such models show that the r-mode overtones diverge and join the inertial modes. In order to see which r-modes persist and retain their character in realistic neutron stars, these calculations will need to be brought into full general relativity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A BOUT++ extension for full annular tokamak edge MHD and turbulence simulations

For tokamak edge plasma simulation, a plasma simulation framework BOUT++ employs a dual coordinate system to simulate moderate-n and high-n plasma instability with reasonable computational cost, where n is the toroidal mode number. This coordinate system however limits the computational domain to the toroidal wedge (full torus divided into N parts in the toroidal direction) for computational efficiency and the use of flute-ordering approximation in the field solver calculating the flow potential from the vorticity which may not be valid for low-n modes. Improving numerical treatment of low-n modes is however indispensable to address simulations of low-n current-driven edge localized mode (ELM), ELM control by resonant magnetic perturbations (RMPs), edge turbulence with RMPs and so on. In this work, BOUT++ is extended to simulate the interplay between $n=0$, low-n and high-n plasma components in a full annular tokamak edge domain through hybrid modeling of the flow potential and the vorticity. Low-n modes of flow potential are calculated in an orthogonal flux surface coordinate and high-n modes in the dual coordinate system separately in Fourier space. Finally, the proposed scheme can capture an interplay between $n=1$ global modes and high-n turbulence during pedestal collapse in a full annular torus domain with a circular cross section.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Arctic Soil Patterns Analogous to Fluid Instabilities: Supporting Data

This dataset characterizes solifluction lobe morphology and spatial patterns using pre-existing LiDAR-derived digital elevation models of 25 sites across Norway with accompanying long term climate data for each site. Data were collected as part of an effort to better understand controls on the formation of solifluction patterns and to test the idea that they are analogous to fluid instabilities. We also provide soil velocity profiles and estimates of effective viscosity from across the world, drawn from literature. They were collected to improve our understanding of the rheology of soliflucting soil. See this article for more information on the theoretical motivation behind this dataset see "Arctic soil patterns analogous to fluid instabilities" (Glade et al., 2021). Data files arranged in a hierarchy and include image files *.tif and *.png, GIS shapefiles and geopackages (*.gpkg), *.csv (with same file as *.xlsx), and *.py (Python scripts readable with a text editor). Files also bundled into *.zip files. Note (2021-10-20): unit corrections made on two files: RR.csv and snowfall.csv. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Anomalous hydrodynamics effective actions from holography

We derive an effective action for charged plasmas with an anomalous (abelian) global current charge current using holography. The holographic description is captured by the dynamics of an Einstein-Maxwell-Chern-Simons theory in an asymptotically AdS spacetime. The ’t Hooft anomaly contribution, which is encoded in the Chern-Simons term, contributes at the Gaussian order in the effective action only in the momentum diffusion sector, where it leads to chiral shear waves. However, as we demonstrate, beyond the Gaussian order, there is non-trivial imprint of the parity-violating anomaly term in sound and charge diffusion dynamics as well.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Momentum transfer and plasma rotation caused by destabilized eigenmodes in tokamaks

In this work, the influence of magnetohydrodynamic eigenmodes destabilized by energetic ions on the momentum of these ions and concomitant sheared plasma rotation are studied. Two mechanisms affecting rotation are revealed: (i) spatial channelling (SC) – radially separated emission and absorption of the momentum; (ii) mode induced redistribution (MIR) across the magnetic field of the momentum of energetic ions. Forces arising during SC and MIR produce both toroidal and poloidal rotations. In addition, the momentum emission during SC leads to a radial flux of fast ions and generation of a radial electric field. Using the developed theory, estimates were made for the ITER (International Thermonuclear Experimental Reactor) 15 MA baseline scenario. They show that a global toroidicity-induced Alfvén eigenmode destabilized by alpha particles and neutral beam injection can result in significant radial electric field and forces applied to plasma. However, available data are not sufficient for a reliable prediction of the effects of SC and MIR in ITER. In general, one can expect that sheared rotation arising after destabilization of Alfvén modes and fast magnetoacoustic modes by energetic ions will tend to suppress the turbulence and improve plasma performance. The importance of plasma rotation is supported, in particular, by the fact that during the JET DTE1 experimental campaign the best parameters were achieved in a deuterium–tritium discharge where the rotation frequency was largest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Achieving a robust grassy-ELM operation regime in CFETR

We have identified a robust grassy-edge localized mode (ELM) operation regime for future tokamak reactors. The regime exists within a pedestal top electron collisionality (ν*) window at high global poloidal beta (β p ). The existence of an upper ν* limit for grassy-ELMs is consistent with results previously reported in experiments (Oyama et al 2010 Nucl. Fusion 50 064014), while the existence of a lower ν* limit has not been reported previously. Using EPED and BOUT++, a theoretical model that quantitatively explains the physics of the grassy-ELMs within the window, which distinguishes them from the small mixed-ELMs at lower ν*, is presented for the first time. A peeling-ballooning stability boundary is obtained by scanning the operating density space. The change in density corresponds to a change in ν* that affects the pedestal bootstrap current. High β p leads to a strong Shafranov shift, which affects the flux surface averaged pressure drive. The two effects combine to create a peeling-dominated window in intermediate ν* buffered by ballooning-dominated regimes. Only the peeling-dominated regime shows a cyclic behavior in the perturbed pressure during the nonlinear simulation of an ELM crash, reminiscent of grassy-ELM dynamics. Similarly, the energy released across the separatrix is demonstrated to be significantly smaller. The quick recovery of the ELM crash is explainable by the rapid rise of a low n kink-peeling instability when the pedestal current I ped exceeds a threshold at high β p . It minimizes the excursion beyond marginal stability and is absent in the ballooning-dominated regime. Comparison with recent experiments over a range of β p and ν* strongly supports the physical picture proposed by the modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

U-Pb zircon age constraints on the earliest eruptions of the Deccan Large Igneous Province, Malwa Plateau, India

Climate instability driven by emission of volatiles during emplacement of large igneous provinces (LIPs) is frequently invoked as a potential cause of mass extinctions. However, documenting this process in the geologic record requires a holistic understanding of eruption rates, the location of eruptive centers, and potential sources of climate-changing volatiles. We present new chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb zircon geochronology from Malwa Plateau basalts on the northern margin of the Deccan LIP, India. These basalts have been previously interpreted as either an extension of the province's main volcanic stratigraphy or as an independent eruptive center active up to millions of years prior to the main eruptive phase. Our data instead demonstrate that the lower Malwa Plateau basalts are temporally correlative with the first pulse of Deccan volcanism and provide new constraints on its initiation and duration. Paleomagnetic data further indicate that upper Malwa Plateau basalts may be age-equivalent to the second, third, and fourth pulses of Deccan volcanism. The relative thicknesses of age-equivalent packages of basalt are consistent with eruption of the Deccan LIP from a southward-migrating eruptive center. The first eruptive pulse is coeval with a ~200 kyr Late Maastrichtian warming event preserved globally in contemporaneous stratigraphic sections. We propose that the first pulse of Deccan magmatism was more voluminous in the north, where it erupted through organic-rich sedimentary rocks of the Narmada-Tapti rift basin. Thermal metamorphism of these sediments could have been a source of sufficient CO 2 to drive the Late Maastrichtian warming event, which if true would reconcile the apparent dampened warming signals associated with later Deccan eruptive pulses.

58 GEOSCIENCES↗
DOE OSTI2021

Liquid-induced topological transformations of cellular microstructures

The fundamental topology of cellular structures—the location, number and connectivity of nodes and compartments—can profoundly affect their acoustic, electrical5, chemical, mechanical and optical properties, as well as heat, fluid and particle transport. Approaches that harness swelling, electromagnetic actuation and mechanical instabilities in cellular materials have enabled a variety of interesting wall deformations and compartment shape alterations, but the resulting structures generally preserve the defining connectivity features of the initial topology. Achieving topological transformation presents a distinct challenge for existing strategies: it requires complex reorganization, repacking, and coordinated bending, stretching and folding, particularly around each node, where elastic resistance is highest owing to connectivity. Here we introduce a two-tiered dynamic strategy that achieves systematic reversible transformations of the fundamental topology of cellular microstructures, which can be applied to a wide range of materials and geometries. Our approach requires only exposing the structure to a selected liquid that is able to first infiltrate and plasticize the material at the molecular scale, and then, upon evaporation, form a network of localized capillary forces at the architectural scale that ‘zip’ the edges of the softened lattice into a new topological structure, which subsequently restiffens and remains kinetically trapped. Reversibility is induced by applying a mixture of liquids that act separately at the molecular and architectural scales (thus offering modular temporal control over the softening–evaporation–stiffening sequence) to restore the original topology or provide access to intermediate modes. Guided by a generalized theoretical model that connects cellular geometries, material stiffness and capillary forces, we demonstrate programmed reversible topological transformations of various lattice geometries and responsive materials that undergo fast global or localized deformations. In conclusion, we then harness dynamic topologies to develop active surfaces with information encryption, selective particle trapping and bubble release, as well as tunable mechanical, chemical and acoustic properties.

42 ENGINEERING↗
DOE OSTI2023

Changing Characteristics of Tropical Extreme Precipitation–Cloud Regimes in Warmer Climates

In this study, we investigated the changing characteristics of climatic scale (monthly) tropical extreme precipitation in warming climates using the Energy Exascale Earth System Model (E3SM). The results are from Atmospheric Model Intercomparison Project (AMIP)-type simulations driven by (a) a control experiment with the present-day sea surface temperature (SST) and CO 2 concentration, (b) P4K, the same as in (a) but with a uniform increase of 4K in the SST globally, and (c) the same as in (a), but with an imposed SST and CO 2 concentration from the outputs of the coupled E3SM forced by a 4xCO 2 concentration. We found that as the surface warmed under P4K and 4xCO 2 , both convective and stratiform rain increased. Importantly, there was an increasing fractional contribution of stratiform rain as a function of the precipitation intensity, with the most extreme but rare events occurring preferentially over land more than the ocean, and more so under 4xCO 2 than P4K. Extreme precipitation was facilitated by increased precipitation efficiency, reflecting accelerated rates of recycling of precipitation cloud water (both liquid and ice phases) in regions with colder anvil cloud tops. Changes in the vertical profiles of clouds, condensation heating, and vertical motions indicate increasing precipitation–cloud–circulation organization from the control and P4K to 4xCO 2 . The results suggest that large-scale ocean warming, that is, P4K, was the primary cause contributing to an organization structure resembling the well-known mesoscale convective system (MCS), with increased extreme precipitation on shorter (hourly to daily) time scales. Additional 4xCO 2 atmospheric radiative heating and dynamically consistent anomalous SST further amplified the MCS organization under P4K. Analyses of the surface moist static energy distribution show that increases in the surface moisture (temperature) under P4K and 4xCO 2 was the key driver leading to enhanced convective instability over tropical ocean (land). However, a fast and large increase in the land surface temperature and lack of available local moisture resulted in a strong reduction in the land surface relative humidity, reflecting severe drying and enhanced convective inhibition (CIN). It is argued that very extreme and rare “record-breaking” precipitation events found over land under P4K, and more so under 4xCO 2 , are likely due to the delayed onset of deep convection, that is, the longer the suppression of deep convection by CIN, the more severe the extreme precipitation when it eventually occurs, due to the release of a large amount of stored surplus convective available potential energy in the lower troposphere during prolonged CIN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Co-Occurring Atmospheric Features and Their Contributions to Precipitation Extremes

Object-based identification algorithms for atmospheric features are commonly utilized to attribute global precipitation. This study employs a systematic approach to examine feature co-occurrences and their relationships to mean and extreme precipitation. Four features are identified using existing data sets for atmospheric rivers (ARs), mesoscale convective systems (MCSs), low-pressure systems (LPSs), and fronts (FTs). Often, a single atmospheric phenomenon satisfies the criteria set by multiple feature identification algorithms, yielding an association between precipitation and multiple features. Over the extra-tropics, the number of features attributed to a single event typically increases with precipitation intensity. Over two-thirds of the precipitation is from co-occurring features, with a considerable fraction related to AR-FT co-occurrences. Over the tropics, about one-quarter of precipitation is associated with co-occurring features, with LPS-MCS co-occurrences contributing substantially in monsoon regions. MCSs are the leading single-feature contributors over tropical land and oceans. In the extra-tropics, FTs, ARs, and their co-occurrences account for over half of the total precipitation over oceans. AR-FT-MCS and FT-MCS co-occurrences contribute to extremes (precipitation exceeding the 95th percentile) over both oceans (over 30%) and land (over 20%). Any combination of features involving MCSs shows a larger contribution to high percentiles of precipitation intensity. A case analysis indicates that AR-FT-MCS co-occurrences exhibit convective instability and deep vertical motion, suggesting that the feature trackers and reanalysis are capturing physics relevant to both convective and frontal systems. The results here emphasize the need for simultaneous identifications of multiple features when attributing precipitation to atmospheric phenomena.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impact of various DIII-D diagnostics on the accuracy of neural network surrogates for kinetic EFIT reconstructions

Abstract Kinetic equilibrium reconstructions make use of profile information such as particle density and temperature measurements in addition to magnetics data to compute a self-consistent equilibrium. They are used in a multitude of physics-based modeling. This work develops a multi-layer perceptron (MLP) neural network (NN) model as a surrogate for kinetic Equilibrium Fitting (EFITs) and trains on the 2019 DIII-D discharge campaign database of kinetic equilibrium reconstructions. We investigate the impact of including various diagnostic data and machine actuator controls as input into the NN. When giving various categories of data as input into NN models that have been trained using those same categories of data, the predictions on multiple equilibrium reconstruction solutions (poloidal magnetic flux, global scalars, pressure profile, current profile) are highly accurate. When comparing different models with different diagnostics as input, the magnetics-only model outputs accurate kinetic profiles and the inclusion of additional data does not significantly impact the accuracy. When the NN is tasked with inferring only a single target such as the EFIT pressure profile or EFIT current profile, we see a large increase in the accuracy of the prediction of the kinetic profiles as more data is included. These results indicate that certain MLP NN configurations can be reasonably robust to different burning-plasma-relevant diagnostics depending on the accuracy requirements for equilibrium reconstruction tasks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

First observations of edge instabilities in strongly shaped negative triangularity plasmas on DIII-D

Throughout the 2023 DIII-D negative triangularity (NT) campaign, unexpected magnetic fluctuations indicative of an edge localized instability were observed in the majority of strongly shaped (–0.55 ≤ δ ≤ –0.45) discharges. Moreover, these same fluctuations were found to be absent in weaker NT plasma shapes. Based on a database study, three separate mode behaviors are observed based primarily on the edge safety factor and normalized plasma beta β N . The main regime is the bursty state, in which short (~5 ms) magnetic bursts occur every 5–10 ms. Each burst drives a small reduction in the edge electron temperature (T e ) as measured by electron cyclotron emissions (ECE), which then rapidly recovers between the bursts. The quiescent regime appears to have governing mechanics similar to the bursty regime, but with multiple unique mode signatures observed in the magnetics. These additional modes greatly reduce the bursty nature of the plasma. Finally, a saturated state is observed in a number of discharges, particularly at early time slices. The saturated regime is often short-lived (maximum 500 ms), but can have a large impact on the dynamics of the plasma edge region, causing significant drops in electron temperature as observed in the ECE diagnostic. 2D ECE imaging (ECE-I) is utilized, showing large kink-like perturbations localized near the last closed flux surface that extend into the SOL indicative of an external kink/peeling behavior. Additionally, evidence of phase mixing and observed tearing-parity perturbations localized near the edge rational surfaces is observed in the quiescent regime plasmas. As a result, the various modes do not appear to have a large impact on global confinement but may play a role in regulating the transport of the NT edge and the observed edge gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Should We Conserve Entropy or Energy when Computing CAPE with Mixed-Phase Precipitation Physics?

Abstract The rapidly increasing resolution of global atmospheric reanalysis and climate model datasets necessitates finding methods for computing convective available potential energy (CAPE) both efficiently and accurately. To this end, this article compares two common methods for computing CAPE which conserve either energy or entropy. Inaccuracies in these computations arise from both physical and numerical errors. For instance, computing CAPE with entropy conserved results in physical errors from nonequilibrium phase transitions but minimizes numerical errors because solutions are analytic at each height. In contrast, computing CAPE with energy conserved avoids these physical errors, but accumulates numerical errors that are grid-resolution-dependent because the numerical integration of a differential equation is required. Analysis of CAPE computed with large databases of soundings from the tropical Amazon and midlatitude storm environments shows that physical errors from the entropy method are typically 1%–3% as large as CAPE, which is comparable to the numerical errors from conserving energy with grid spacing of 25 and 250 m using explicit first-order and second-order integration schemes, respectively. Errors in entropy-based CAPE calculations are also insensitive to vertical grid spacing, in contrast to energy-based calculations whose error strongly scales with the grid spacing. It is shown that entropy-based methods are advantageous when intercomparing datasets with differing vertical resolution because they produce accurate and reasonably fast results that are insensitive to grid resolution, whereas a second-order energy-based method is advantageous when analyzing data with a consistent vertical resolution because of its superior computational efficiency. Significance Statement Convective available potential energy (CAPE) is a measure of instability in the atmosphere that helps forecasters and researchers understand when and where thunderstorms will form. The purpose of this article is to identify the most efficient and accurate methods for computing CAPE. Two methods are considered here, one that relates to the entropy (a measure of thermodynamic disorder) of an air parcel and one that relates to the energy of an air parcel. Results indicate that the entropy method is most accurate and insensitive to the resolution of the data used for the calculation (which can vary considerably), whereas the energy method uses the least computation time.

Peters, John M.↗
DOE OSTI2025

Turbulent burning velocity of lean premixed hydrogen/air flames at engine conditions: Effects of turbulence intensity and length scale

For turbulent lean premixed hydrogen flames with strong thermodiffusively instabilities, most previous studies have focused on the influence of turbulence intensity, whereas the role of turbulence length scale is less well understood. Here, this study addresses this gap by conducting direct numerical simulations (DNS) of statistically planar turbulent premixed flames for a lean (ϕ=0.35) hydrogen/air mixture under independently varied turbulence intensity (u') and length scale (l T ) at engine-relevant thermodynamics conditions. Results show that as u' increases, the flame front becomes increasingly wrinkled, forming smaller cellular structures. In contrast, l T variations do not significantly alter the size of these structures. For the turbulent burning velocity (S T ), the normalized S T (i.e., S T /S L , where S L is the laminar flame speed) increases linearly with u', driven by both enhanced flame surface wrinkling (i.e., increased A T /A L ) and enhanced local burning rate (i.e., increased I 0 ). However, increasing l T reduces I 0 , despite a continued increase in A T /A L , resulting in only a marginal increase in S T /S L . To reveal the underlying mechanisms, especially the decreasing trend of I 0 with l T , local flame dynamics analyses are performed. It is found that as l T increases, the interaction between thermodiffusive effects and turbulence weakens due to the reduced tangential strain rate, while the flame curvature remains largely unchanged. This suppresses local reactivity enhancement and thus decreases I 0 , In contrast, an increase in u' enhances the interaction by amplifying both curvature fluctuation and tangential strain rate, leading to increased local reactivity (increased I 0 ). Finally, based on the DNS data, several new scaling models are proposed for the three global properties, S T /S L , A T /A L , and I 0 , and showed improvements compared to existing models. These findings provide new insights into the flame-turbulence interactions in thermodiffusively unstable hydrogen flames. The DNS dataset is also useful for the development of turbulent combustion models applicable to practical engine simulations.

Engine-relevant condition↗
DOE OSTI2020

Event Horizon Telescope imaging of the archetypal blazar 3C 279 at an extreme 20 microarcsecond resolution

3C 279 is an archetypal blazar with a prominent radio jet that show broadband flux density variability across the entire electromagnetic spectrum. We use an ultra-high angular resolution technique – global Very Long Baseline Interferometry (VLBI) at 1.3 mm (230 GHz) – to resolve the innermost jet of 3C 279 in order to study its fine-scale morphology close to the jet base where highly variable γ-ray emission is thought to originate, according to various models. The source was observed during four days in April 2017 with the Event Horizon Telescope at 230 GHz, including the phased Atacama Large Millimeter/submillimeter Array (ALMA), at an angular resolution of ∼20 μas (at a redshift of z = 0.536 this corresponds to ∼0.13 pc ∼ 1700 Schwarzschild radii with a black hole mass MBH = 8 × 108 M⊙). Imaging and model-fitting techniques were applied to the data to parameterize the fine-scale source structure and its variation. We find a multicomponent inner jet morphology with the northernmost component elongated perpendicular to the direction of the jet, as imaged at longer wavelengths. The elongated nuclear structure is consistent on all four observing days and across different imaging methods and model-fitting techniques, and therefore appears robust. Owing to its compactness and brightness, we associate the northern nuclear structure as the VLBI “core”. This morphology can be interpreted as either a broad resolved jet base or a spatially bent jet. We also find significant day-to-day variations in the closure phases, which appear most pronounced on the triangles with the longest baselines. Our analysis shows that this variation is related to a systematic change of the source structure. Two inner jet components move non-radially at apparent speeds of ∼15 c and ∼20 c (∼1.3 and ∼1.7 μas day−1, respectively), which more strongly supports the scenario of traveling shocks or instabilities in a bent, possibly rotating jet. The observed apparent speeds are also coincident with the 3C 279 large-scale jet kinematics observed at longer (cm) wavelengths, suggesting no significant jet acceleration between the 1.3 mm core and the outer jet. The intrinsic brightness temperature of the jet components are ≲1010 K, a magnitude or more lower than typical values seen at ≥7 mm wavelengths. The low brightness temperature and morphological complexity suggest that the core region of 3C 279 becomes optically thin at short (mm) wavelengths.Key words: galaxies: active / galaxies: jets / galaxies: individual: 3C 279 / techniques: interferometric⋆ The data are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/640/A69 and at https://eventhorizontelescope.org/for-astronomers/data⋆⋆ Deceased.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Comparison of edge turbulence characteristics between DIII-D and C-Mod simulations with XGC1

The physical processes taking place at the separatrix and scrape-off layer regions are crucial for the operation of tokamaks as they govern the interaction of hot plasma with the vessel walls. Numerical modeling of the edge with state-of-the-art codes attempts to elucidate the complex interactions between neoclassical drifts, turbulence, poloidal and parallel flows that control the physical set-up of the SOL region. Here, we present post-processing analysis of simulation results from the gyrokinetic code XGC1, comparing and contrasting edge turbulence characteristics from a simulation of the DIII-D tokamak against a simulation of the Alcator C-Mod tokamak. We find that the equilibrium ExB flux across the separatrix has a similar poloidal pattern in both discharges which can be explained by rB-drifts and trapped ion excursions. Furthermore, collisionality is noted to play a major role in the way that it prevents local charge accumulations from having more global effects in the C-Mod case. In both cases, turbulent electron heat flux is observed to be higher than the ion one and is possibly related to the need of electrons to maintain quasineutrality through the only channel available to them for exiting the confinement. By Fourier analysis, we identify turbulent frequencies and growth rates of the dominant mode in both simulations. In the case of C-Mod, these numbers point to the presence of a drift wave. In the DIII-D case, further linear simulations with the Gene code reveal a trapped electron mode. Furthermore, using a blob detection and tracking tool, we present the amplitude and size distributions of the blobs from both simulations. The amplitude distributions are in qualitative agreement with experimental observations while the size distributions are consistent with the fact that most of the blobs are not connecting to the divertor plates and suggest that they are generated by the shearing of the turbulent modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gyrokinetic simulation of low-frequency Alfvénic modes in DIII-D tokamak

Global gyrokinetic simulations find that a beta-induced Alfvén eigenmode (BAE) and a Low frequency mode (LFM) coexist in the DIII-D tokamak experiments. The LFM can be excited without fast ions and has a frequency inside the gap of the beta-induced Alfvén-acoustic eigenmode (BAAE). The mode structure of the LFM from simulations and many of its parametric dependencies are consistent with experimental observations. For the BAE, the simulated mode structure is consistent with experiment but the frequency is lower than experiment. An antenna scan shows that the LFM is NOT the conventioanl BAAE; rather, it is an interchange-like electromagnetic mode excited by non-resonant drive of pressure gradients. The compressible magnetic perturbations significantly increase the growth rates of the BAE and LFM. On the other hand, effects of trapped electrons and equilibrium current have modest effects on the BAE and LFM. Finally, radial transfer of wave energy is different between the BAE and LFM, resulting in their respective triangular and elliptic mode shapes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗