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 127 records · Page 7

Review of Power System Support Functions for Inverter-Based Distributed Energy Resources - Standards, Control Algorithms, and Trends

Penetration of renewable energy in power systems has been increasing in the past decades in response to increased global electricity demand and concerns for the environment. Distributed energy resources (DERs) based on renewables have experienced rapid growth thanks to the incentive programs and broad-based participation. With the growing prevalence of DERs, the risk of grid instability and vulnerability increases due to the intermittent nature of renewable energy. At the same time, the voltage and frequency deviation problems emerge more often when the reverse power flow occurs under supply-demand imbalance in distributed power systems. Standards and grid codes have been issued for DER inverters to interconnect with the distribution grid. The updated standard and grid codes expect DERs to provide a variety of power system support functions in order to incorporate higher DER penetration and to maximize DER value to the grid. This paper provides an overview of the power system support functions from renewable DER inverters, which are categorized as: voltage regulation by active/reactive power control, frequency regulation by active power control, voltage ride-through, and frequency ride-through. The benefits and drawbacks of each algorithm are presented and compared with its predecessor, manifesting the logic in the evolution of the algorithms.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Tornadoes in Southeast South America: Mesoscale to Planetary-Scale Environments

Abstract A multiscale analysis of the environment supporting tornadoes in southeast South America (SESA) was conducted based on a self-constructed database of 74 reports. Composites of environmental and convective parameters from ERA5 were generated relative to tornado events. The distribution of the reported tornadoes maximizes over the Argentine plains, while events are rare close to the Andes and south of Sierras de Córdoba. Events are relatively common in all seasons except in winter. Proximity environment evolution shows enhanced instability, deep-layer vertical wind shear, storm-relative helicity, reduced convective inhibition, and a lowered lifting condensation level before or during the development of tornadic storms in SESA. No consistent signal in low-level wind shear is seen during tornado occurrence. However, a curved hodograph with counterclockwise rotation is present. The Significant Tornado Parameter (STP) is also maximized prior to tornadogenesis, most strongly associated with enhanced CAPE. Differences in the convective environment between tornadoes in SESA and the U.S. Great Plains are discussed. On the synoptic scale, tornado events are associated with a strong anomalous trough crossing the southern Andes that triggers lee cyclogenesis, subsequently enhancing the South American low-level jet (SALLJ) that increases moisture advection to support deep convection. This synoptic trough also enhances vertical shear that, along with enhanced instability, sustains organized convection capable of producing tornadic storms. At planetary scales, the tornadic environment is modulated by Rossby wave trains that appear to be forced by convection near northern Australia. Madden–Julian oscillation phase 3 preferentially occurs 1–2 weeks ahead of tornado occurrence. Significance Statement The main goal of this study is to describe what atmospheric conditions (from local to global scales) are present prior to and during tornadic storms impacting southeast South America (SESA). Increasing potential for deep convection, wind shear, and potential for rotating updrafts, as well as reducing convective inhibition and cloud-base height, are predominant a few hours before and during the events in connection to low-level northerly winds enhancing moisture transport to the region. Remote convective activity near northern Australia appears to influence large-scale atmospheric circulation that subsequently triggers convective storms supporting tornadogenesis 1–2 weeks later in SESA. Our findings highlight the importance of accounting for atmospheric processes occurring at different scales to understand and predict tornado occurrences.

Meteorology & Atmospheric Sciences↗

Linear simulation of kinetic electromagnetic instabilities in a tokamak plasma with weak magnetic shear

Gyrokinetic simulation and eigenvalue calculation of electromagnetic instabilities are carried out for an experimentally observed low-n mode in weak magnetic shear discharge. With different magnetic shear values, the Gyrokinetic Toroidal Code simulation of the ion temperature gradient mode is consistent with the eigenvalue code (HD7) calculation. Due to the sensitivity of the kinetic ballooning mode (KBM) to global equilibrium, the simulation of the KBM deviates from the eigenvalue results, for the ballooning representation used in HD7 satisfies its spatial scale separations. Under a flat safety factor profile, the KBM is more unstable and its mode structure tends to move with the peak of the ion temperature drive. Further simulation of the KBM in an HL-2A-like equilibrium shows that the β excitation threshold of the mode is lower than 0.2% and the dominant toroidal mode number is n = 4, which is consistent with the measured experimental spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Achieving American Leadership in Cybersecurity and Digital Components Factsheet

The Biden Administration’s efforts to meet 100% clean electricity by 2035 and net-zero greenhouse gas emissions by 2050 has created opportunities to rebuild American supply chains create new jobs, strengthen community engagement, and spur U.S. economic growth. DOE’s strategy for success in the transition to a clean energy economy hinges on building and maintaining technology supply chains that are advanced, secure, and resilient to cyber threats. As the energy sector grows increasingly globalized, complex, and digitized, the supply chain for digital components of energy systems – including software, virtual platforms and services, and data – is facing greater threats. Nearly all digital components of U.S. energy sector systems are vulnerable to cyber supply chain instability, stemming from a variety of causes and shared among a broad set of interdependent stakeholders. Overall, supply chain risks for digital components in energy sector systems will continue to evolve and likely increase as these systems are increasingly interconnected, digitized, and remotely operated.

CCE↗

Implementation and Exploration of Parameterizations of Large-Scale Dynamics in NCAR's Single Column Atmosphere Model SCAM6

A single column model with parameterized large-scale (LS) dynamics is used to better understand the response of steady-state tropical precipitation to relative sea surface temperature under various representations of radiation, convection, and circulation. The large-scale dynamics are parametrized via the weak temperature gradient (WTG), damped gravity wave (DGW), and spectral weak temperature gradient (Spectral WTG) method in NCAR's Single Column Atmosphere Model (SCAM6). Radiative cooling is either specified or interactive, and the convective parameterization is run using two different values of a parameter that controls the degree of convective inhibition. Results are interpreted in the context of the Global Atmospheric System Studies -Weak Temperature Gradient (GASS-WTG) Intercomparison project. Using the same parameter settings and simulation configuration as in the GASS-WTG Intercomparison project, SCAM6 under the WTG and DGW methods produces erratic results, suggestive of numerical instability. However, when key parameters are changed to weaken the large-scale circulation's damping of tropospheric temperature variations, SCAM6 performs comparably to single column models in the GASS-WTG Intercomparison project. The Spectral WTG method is less sensitive to changes in convection and radiation than are the other two methods, performing qualitatively similarly across all configurations considered. Under all three methods, circulation strength, represented in 1D by grid-scale vertical velocity, is decreased when barriers to convection are reduced. This effect is most extreme under specified radiative cooling, and is shown to come from increased static stability in the column's reference radiative-convective equilibrium profile. This argument can be extended to interactive radiation cases as well, though perhaps less conclusively.

54 ENVIRONMENTAL SCIENCES↗

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global Gyrokinetic Simulations of Isotope Effects under Ambipolar Electric Fields and Advances Toward Whole-Volume Modeling

We review global gyrokinetic simulation studies on plasma transport in the Large Helical Device using XGC-S. XGC-S is an extended version of X-point Gyrokinetic Code for stellarators and has been progressively verified throughout the code development process. Verification tests of neoclassical transport successfully demonstrate the generation of an ambipolar electric field due to ripple-trapped particles. We perform quasi-linear analyses of the ion temperature gradient mode under the influence of the ambipolar electric field. The results reveal that the ambipolar electric field and the heavy hydrogen component in mixed isotope plasmas can lead to the favorable isotope effect observed in recent deuterium experiments. We also present recent efforts in code development toward whole-volume simulations, including the helical divertor region. A mesh generation scheme based on field-line tracing and the construction of curved surfaces perpendicular to the magnetic field would be promising for global field calculations in the whole-volume simulations.

Basic Plasma Phenomena and Gas Discharges↗

Magnetopause Reconnection and Indents Induced by Foreshock Turbulence

Abstract Based on global hybrid simulation results, we predict that foreshock turbulence can reach the magnetopause and lead to reconnection as well as Earth‐sized indents. Both the interplanetary magnetic field (IMF) and solar wind are constant in our simulation, and hence, all dynamics are generated by foreshock instabilities. The IMF in the simulation is mostly Sun‐Earth aligned with a weak northward and zero dawn‐dusk component, such that subsolar magnetopause reconnection is not expected without foreshock turbulence modifying the magnetosheath fields. We show a reconnection example to illustrate that the turbulence can create large magnetic shear angles across the magnetopause to induce local bursty reconnection. Magnetopause reconnection and indents developed from the impact of foreshock turbulence can potentially contribute to dayside loss of planetary plasmas.

58 GEOSCIENCES↗

Fractional Laplacian spectral approach to turbulence in a dusty plasma monolayer

This work presents an analytical investigation of anomalous diffusion and turbulence in a dusty plasma monolayer, where energy transport across scales leads to the spontaneous formation of spatially disordered patterns. Many-body simulations of 10,000-particle dusty plasma monolayers are used to demonstrate how the global dynamics depend on the statistical properties of the dust assembly for realistic laboratory conditions. We find that disorder due to variations in the dust size distribution and charge-driven nonlocal interactions resulting in anomalous dust diffusion are key factors for the onset of instabilities. The resulting dynamics exhibit features of inertial turbulence over slightly more than half a decade of scales proportional or smaller than the Debye shielding length. These processes are examined analytically using a recently developed Fractional Laplacian Spectral (FLS) technique, which identifies the active energy channels as a function of scale, disorder concentration, and features of the nonlocal-interactions. The predictions from the theoretical (spectral) analysis demonstrate agreement with the results from the many-body (kinetic) simulations, thus providing a powerful tool for the study of active turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetopause Reconnection as Influenced by the Dipole Tilt Under Southward IMF Conditions: Hybrid Simulation and MMS Observation

Abstract Using a three‐dimensional (3‐D) global‐scale hybrid code, the Magnetospheric Multiscale (MMS) reconnection event around 02:13 UT on 18 November 2015, highlighted in the Geospace Environment Modeling (GEM) Dayside Kinetic Challenge, is simulated, in which the interplanetary magnetic field (IMF) points southward and the geomagnetic field has a −27° dipole tilt angle. Strong southward plasma jets are found near the magnetopause as a result of the dayside reconnection. Our results indicate that the subsolar magnetopause reconnection X line shifts from the subsolar point toward the Northern Hemisphere due to the effect of the tilted geomagnetic dipole angle, consistent with the MMS observation. Subsequently, the reconnection X lines or sites and reconnection flux ropes above the equator propagate northward along the magnetopause. The formation and global distribution of the X lines and the structure of the magnetopause reconnection are investigated in detail with the simulation. Mirror mode waves are also found in the middle of the magnetosheath downstream of the quasi‐perpendicular shock where the plasma properties are consistent with the mirror instability condition. As a special outcome of the GEM challenge event, the spatial and temporal variations in reconnection, the electromagnetic power spectra, and the associated D‐shaped ion velocity distributions in the simulated reconnection event are compared with the MMS observation.

Guo, Zhifang↗

Data Assimilation with Machine Learning Surrogate Models: A Case Study with FourCastNet

Modern data-driven surrogate models for weather forecasting provide accurate short-term predictions but inaccurate and nonphysical long-term forecasts. This paper investigates online weather prediction using machine learning surrogates supplemented with partial and noisy observations. We empirically demonstrate and theoretically justify that, despite the long-time instability of the surrogates and the sparsity of the observations, filtering estimates can remain accurate in the long-time horizon. As a case study, we integrate the Fourier Forecasting Neural Network (FourCastNet), a weather surrogate model, within a variational data assimilation framework using partial, noisy ERA5 global reanalysis data from the European Centre for Medium-Range Weather Forecasts (ECMWF). Here, our results show that filtering estimates remain accurate over a year-long assimilation window and provide effective initial conditions for forecasting tasks, including extreme event prediction.

Data assimilation↗

Self-organized confinement in whole-device modeling of laboratory magnetospheres

Turbulent self-organization driven by global ideal interchange modes in a dipole-confined plasma is explored with self-consistent, whole-plasma simulations using a flux-tube averaged magnetohydrodynamic model in dipole magnetic geometry. We show the existence of robust particle pinch driven by ideal interchange-mode fluctuations, in which the particles are transported up the density gradient. It is found that the plasma profiles in a dipole field spontaneously relax to a marginally stable state as centrally peaked pressure and density are created by global interchange-mode transport.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2)

Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

54 ENVIRONMENTAL SCIENCES↗

Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry

For the research of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Empirical Evaluation of GPS Clock Accuracy for Isochronous Droop-based Inverters

Grid-forming inverters rely upon an internal reference to regulate the grid 's voltage and frequency. Without a reliable reference, issues such as loss of synchronism, frequency instability, and unequal power sharing can occur. One proposed internal phase reference for grid-forming inverters is generated by the controller's clock, but this is affected by nonidealities in the crystal oscillator leading to time drift in the phase reference. To improve the controller's clock accuracy, the Global Positioning System (GPS) can be used to calibrate the internal clock. Previously, GPS-based clock calibration was explored theoretically and various methods were developed for mitigating frequency/phase drift at the power system level. This paper experimentally validates the direct performance of GPS time-calibration at the controller level and justifies its use in maintaining an accurate phase reference for droop-based grid-forming inverters.

clock drift↗

Empirical Evaluation of GPS Clock Accuracy for Isochronous Droop-Based Inverters: Preprint

Grid-forming inverters rely upon an internal reference to regulate the grid’s voltage and frequency. Without a reliable reference, issues such as loss of synchronism, frequency instability, and unequal power sharing can occur. One proposed internal phase reference for grid-forming inverters is generated by the controller’s clock, but this is affected by nonidealities in the crystal oscillator leading to time drift in the phase reference. To improve the controller’s clock accuracy, the Global Positioning System (GPS) can be used to calibrate the internal clock. Previously, GPS-based clock calibration was explored theoretically and various methods were developed for mitigating frequency/phase drift at the power system level. This paper experimentally validates the direct performance of GPS time-calibration at the controller level and justifies its use in maintaining an accurate phase reference for droop-based grid-forming inverters.

clock drift↗

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗