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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Database-wide hazard modelling of the onset of DIII-D tearing modes with field features

The rate of onset (hazard) of tearing modes is modelled probabilistically using statistical learning algorithms. Axisymmetric energy-density equilibrium fields are taken as raw high-dimensional input features which are reduced with principal component analysis. Signal processing of non-axisymmetric magnetics fluctuation array data provides the target information from which to learn. Model selection, visualization and calibration assessment procedures are detailed. Here, the analysis is deployed at large scale across the DIII-D tokamak database. Standard model selection criteria suggest that the energy-density post-processed feature is a better choice for modelling the onset rate compared to the non-processed equilibrium reconstruction solution. Two example applications of the learned rate function are demonstrated: (i) proximity-to-onset discharge monitoring and (ii) database analysis showing an (expected) observational global trend that the general hazard increases as a plasma performance metric increases. An important connection between the hazard function and its use as a conditional probability generator is reviewed in the Appendix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Improved CO 2 Capture Performance of Electrochemically Mediated Amine Regeneration Processes with Ionic Surfactant Additives

For the effective reduction of global CO 2 emissions, it is essential to develop and deploy efficient and cost-effective technologies for CO 2 capture, especially from large point sources. We recently developed an electrochemically mediated amine regeneration (EMAR) system to replace traditional thermal desorption for the capture of CO 2 from post-combustion flue gases. Despite EMAR effectiveness on a laboratory scale, concerns regarding the high gas-to-liquid ratio in the electrochemical cell and long-term instability of the electrodes need to be addressed before further scale-up of the process to a pilot plant and beyond can be entertained. Accordingly, we investigated the effect of using sodium dodecyl sulfate (SDS) as an anionic surfactant and dodecyltrimethylammonium bromide (DTAB) as a cationic surfactant on the process operation. It was found that it is advantageous to use an anionic surfactant for a system such as EMAR that contains hydrophilic electrodes and a positively charged electrochemically active species. The overall cell resistance was notably reduced when SDS anionic surfactant was used. The precipitation of copper particles observed in the anode outlet when no surfactant was used was effectively avoided when SDS was added to the electrolyte, resulting in electrode stability. In addition, smaller gas bubbles were produced in the presence of the SDS surfactant, which resulted in less blockage of the electrode by the gas with a resultant lower cell potential under constant current conditions, driving more efficient CO 2 desorption. This led to an ~25% reduction in the electrochemical energy requirement, the lowest ever achieved experimentally for the EMAR process. Altogether, the addition of a very low concentration of SDS resulted in the successful circumvention of the important problems faced by the EMAR system regarding further scale-up.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

High-order dimensionally-split Cartesian embedded boundary method for non-dissipative schemes

Centered finite-difference schemes are commonly used for high-fidelity turbulent flow simulations in canonical configurations because of their non-dissipative property and computational efficiency. However, their use in flow simulations over complex geometries is limited by the requirements of a structured grid and a stable boundary treatment in the absence of artificial (numerical) dissipation. Cartesian embedded boundary (EB) approaches provide an efficient structured-grid framework to apply difference schemes over complex domains. However, they are often restricted to low orders of accuracy because of numerical instabilities at the embedded boundaries and the issues of small-cell problem that are difficult to address with high-order accuracy. The present work discusses a systematic approach to obtain high-order EB methods with non-dissipative centered schemes in the interior. This approach, based on satisfying the primary and secondary conservation conditions, is employed to derive EB schemes that are up to sixth-order accurate in the interior and fourth-order accurate globally for hyperbolic, parabolic as well as incompletely parabolic problems. The proposed finite-difference discretization is, by construction, dimensionally split and addresses the small-cell problem without any cell/geometry transformations, thus, highly simplifying implementation in a flow solver. Various linear and non-linear numerical tests are performed to evaluate the stability and the accuracy of the proposed EB schemes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Nonlocal Coulomb interaction and spin-freezing crossover as a route to valence-skipping charge order

Abstract Multiorbital systems away from global half-filling host intriguing physical properties promoted by Hund’s coupling. Despite increasing awareness of this regime dubbed Hund’s metal, effect of nonlocal interaction is still elusive. Here we study a three-orbital model with 1/3 filling (two electrons per site) including the intersite Coulomb interaction ( V ). Using the G W plus extended dynamical mean-field theory, the valence-skipping charge order transition is shown to be driven by V . Most interestingly, the instability to this transition is significantly enhanced in the spin-freezing crossover regime, thereby lowering the critical V to the formation of charge order. This behavior is found to be closely related to the population profile of the atomic multiplet states in the spin-freezing regime. In this regime, maximum spin states are dominant in each total charge subspace with substantial amount of one- and three-electron occupations, which leads to almost equal population of one- and the maximum spin three-electron state. Our finding unveils another feature of the Hund’s metal and has potential implications for the broad range of multiorbital systems as well as the recently discovered charge order in iron pnictides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

A study of adhesive bonding in metal–metal, metal–CFRP, and CFRP–CFRP material combinations under shear deformation: Fracture morphologies and damage mechanisms

Safe design of adhesive joining in multi-materials in engineered structures requires the accumulation of numerous experimental data on the failure behavior of various adhesively-bonded material combinations under different loading conditions. The deep understanding of mechanical performance, fracturing morphologies, and main damage mechanisms is also quintessential for accelerating the development of proper physics-based and multi-scale models for assisting the design. Towards this goal, this work presents a comprehensive characterization of the failure behavior of adhesively-bonded metal–metal, metal–CFRP, and CFRP–CFRP material combinations under global shear deformation via single lap shear testing. Thanks to a synergistic combination of measurement methods by using Digital Imaging Correlation (DIC) and 3D optical profilometry, adhesive features on the adherend after failure were quantified and the main progressive damage mechanisms were identified. The characterization performed in this work provides quantitative data that contributes to a better understanding of shear failure in adhesive bonding across different bi-material combinations. Finally, the obtained results have practical implications, including the potential to enhance adhesive bonding design, identify failure causes in adhesive joints, and develop or validate computational models capable of capturing the observed behavior in various adhesively-bonded materials under global shear deformation.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Overview and technology opportunities for thermochemically-produced bio-blendstocks

Global demand for transportation fuels is projected to increase 40% by 2040, and biomass-derived fuels (biofuels) play a crucial role in substituting fossil fuels and mitigating greenhouse gas emissions. Currently, biofuels are mainly consumed as blendstocks combined with petroleum-based fuels, and effective conversion technologies can address the quality challenges for offering standalone biofuels. Thermochemical conversion process is one of the most promising pathways among existing technologies for biofuel production. However, the major barriers are unwanted characteristics (e.g., thermal instability) of intermediate products, such as bio-oil, and required upgrading treatments for producing compatible fuels. Here, this study highlights the merits and critical challenges of thermochemical conversion and physicochemical upgrading technologies for bio-blendstock production from lignocellulosic biomass. The novelty of this study lies in potential directions for future research through both critical and systematic literature reviews, and the proposed intensified process for lignocellulosic-based fuel blendstocks production. It is concluded that recovery and fractionation strategies (e.g., quenching and stripping) can maximize process yields and add values in the efficient conversion pathways. Effective quenching can stop secondary free radical reactions and improve liquid yields over gas and solid yields. Stripping process can improve process yield, catalyst lifespan, and thermal stability. It is further concluded that physicochemical treatments are not as effective as thermochemical treatments, but have advantages of mild operating conditions and potential for integrated solutions in conjunction with other treatments.

09 BIOMASS FUELS↗

DOE OSTI2022

A High-End Estimate of Sea Level Rise for Practitioners

Sea level rise (SLR) is a long-lasting consequence of climate change because global anthropogenic warming takes centuries to millennia to equilibrate for the deep ocean and ice sheets. SLR projections based on climate models support policy analysis, risk assessment and adaptation planning today, despite their large uncertainties. The central range of the SLR distribution is estimated by process-based models. However, risk-averse practitioners often require information about plausible future conditions that lie in the tails of the SLR distribution, which are poorly defined by existing models. Here, a community effort combining scientists and practitioners builds on a framework of discussing physical evidence to quantify high-end global SLR for practitioners. The approach is complementary to the IPCC AR6 report and provides further physically plausible high-end scenarios. High-end estimates for the different SLR components are developed for two climate scenarios at two timescales. For global warming of +2°C in 2100 (RCP2.6/SSP1-2.6) relative to pre-industrial values our high-end global SLR estimates are up to 0.9 m in 2100 and 2.5 m in 2300. Similarly, for a (RCP8.5/SSP5-8.5), we estimate up to 1.6 m in 2100 and up to 10.4 m in 2300. The large and growing differences between the scenarios beyond 2100 emphasize the long-term benefits of mitigation. However, even a modest 2°C warming may cause multi-meter SLR on centennial time scales with profound consequences for coastal areas. Earlier high-end assessments focused on instability mechanisms in Antarctica, while here we emphasize the importance of the timing of ice shelf collapse around Antarctica. This is highly uncertain due to low understanding of the driving processes. Hence both process understanding and emission scenario control high-end SLR.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

A Decadal Hybrid GCM Simulation Using Deep‐Learning‐Based Cloud and Convection Parameterization Generalized to a Warm Climate

A critical challenge for machine‐learning (ML) parameterization in global climate models (GCMs) is to achieve stable, accurate simulations under climates not seen during training. Previous studies have demonstrated promising offline performance and year‐long online stability in aquaplanet simulations but have encountered difficulties in real geography and under climate warming. Here we report that a GCM with real geography configuration using neural‐network‐based cloud and convection parameterization, trained exclusively with present‐day climate data, successfully performs a stable, decade‐long simulation of a warm climate with +4 K sea surface temperature (SST). The neural network (NN) is based on Han et al. (2023, https://doi.org/10.1029/2022ms003508 ) with additional inputs. The simulation captures the global precipitation distribution, surface temperatures, vertical atmospheric structures, and extreme precipitation very well, closely matching simulations from both the superparameterized CAM (SPCAM) and the conventional CAM5 in the warm climate without accuracy degradation compared to those in the baseline climate. Moreover, it produces a climate response to +4 K SST in atmospheric thermodynamic states and circulations similar to those from SPCAM and CAM5. Prognostic ablation tests on NN input variables show that the NN without convective memory as input suffers from numerical instability, and the NN without considering radiative variables and land fraction as input, or with reduced training samples produce less accurate results. To our knowledge, this is the first time an ML parameterization successfully achieves online extrapolation to a warm climate without using additional warm‐climate data for training. It demonstrates the potential of ML‐driven parameterizations for credible long‐term climate projections.

Atmosphere model↗

DOE OSTI2023

Semi-Implicit Computation of Fast Modes in a Scheme Integrating Slow Modes by a Leapfrog Method Based on a Selective Implicit Time Filter

Abstract A scheme for integration of atmospheric equations containing terms with differing time scales is developed. The method employs a filtered leapfrog scheme utilizing a fourth-order implicit time filter with one function evaluation per time step to compute slow-propagating phenomena such as advection and rotation. The terms involving fast-propagating modes are handled implicitly with an unconditionally stable method that permits application of larger time steps and faster computations compared to fully explicit treatment. Implementation using explicit and recurrent formulation is provided. Stability analysis demonstrates that the method is conditionally stable for any combination of frequencies involved in the slow and fast terms as they approach the origin. The implicit filter used in the method damps the computational modes without noticeably sacrificing the accuracy of the physical mode. TheO[(Δt 4 )] accuracy for amplitude errors achieved by the implicitly filtered leapfrog is preserved in applications where terms responsible for fast propagation are integrated with a semi-implicit method. Detailed formulation of the method for soundproof nonhydrostatic anelastic equations is provided. Procedures for implementation in global spectral shallow-water models are also given. Examples comparing numerical and analytical solutions for linear gravity waves demonstrate the accuracy of the scheme. The performance is also shown in more practical nonlinear applications, where numerical solutions accomplished by the method are evaluated against those computed from a scheme where the slow terms are handled by the third-order Runge–Kutta scheme. It demonstrates that the method is able to accurately resolve fine-scale dynamics of Kelvin–Helmholtz shear instabilities, the evolution of density current, and nonlinear drifts of twin tropical cyclones.

Meteorology & Atmospheric Sciences↗

DOE OSTI2025

Sheath transitions in a cylindrical filament discharge: Axisymmetric 1D3V PIC-MCC simulations

We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.

Electrical hysteresis↗

DOE OSTI2022

Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients

We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Two Mechanisms Limiting the Emitted Electron Current from a Cathode to an Anode

It is known that the current of emitted electrons flowing through a plasma can saturate upon formation of a potential well adjacent to the cathode (the “space charge effect”). Here, we demonstrate another saturation mechanism that will often set a more restrictive limit on the global current. When “backflow saturation” occurs, the cathode sheath weakens to allow emitted electrons that already entered the plasma to backflow to the cathode. This effect could not be captured by studies modelling the cathode sheath by itself because its origin is coupled to processes in the interior plasma and anode sheath. By modeling a full plasma diode, we show that depending on conditions the global current can be limited in four ways; by backflow alone, by space charge alone, by both mechanisms in a stable cooperative form, or by both in a competing oscillatory form. Published by the American Physical Society 2025

Laboratory plasma↗

DOE OSTI2026

Crowdsourcing the Frontier: Advancing Hybrid Physics‐ML Climate Simulation via a $\$$50,000 Kaggle Competition

Subgrid machine-learning (machine learning [ML]) parameterizations have the potential to introduce a new generation of climate models that incorporate the effects of higher-resolution physics without incurring the prohibitive computational cost associated with more explicit physics-based simulations. However, important issues, ranging from online instability to inconsistent online performance, have limited their operational use for long-term climate projections. To more rapidly drive progress in solving these issues, domain scientists and ML researchers opened up the offline aspect of this problem to the broader ML and data science community with the release of ClimSim, a NeurIPS Data sets and Benchmarks publication, and an associated Kaggle competition. This paper reports on the downstream results of the Kaggle competition by coupling emulators inspired by the winning teams' architectures to an interactive climate model (including full cloud microphysics, a regime historically prone to online instability) and systematically evaluating their online performance. Our results demonstrate that online stability in the low-resolution real-geography setting is reproducible across multiple diverse architectures, which we consider a key milestone. All tested architectures exhibit strikingly similar offline and online biases, though their responses to architecture-agnostic design choices (e.g., expanding the list of input variables) can differ significantly. Multiple Kaggle-inspired architectures achieve state-of-the-art results on certain metrics such as zonal mean bias patterns and global Root Mean Squared Error, indicating that crowdsourcing the essence of the offline problem is one path to improving online performance in hybrid physics-AI climate simulation.

Environmental sciences↗

DOE OSTI2022

The Role of Strong Magnetic Fields in Stabilizing Highly Luminous Thin Disks

Abstract We present and analyze a set of three-dimensional, global, general relativistic radiation magnetohydrodynamic simulations of thin, radiation-pressure-dominated accretion disks surrounding a nonrotating, stellar-mass black hole. The simulations are initialized using the Shakura–Sunyaev model with a mass accretion rate of M ̇ = 3 L Edd / c 2 (corresponding to L = 0.17 L Edd ). Our previous work demonstrated that such disks are thermally unstable when accretion is driven by an α -viscosity. In the present work, we test the hypothesis that strong magnetic fields can both drive accretion through magnetorotational instability and restore stability to such disks. We test four initial magnetic field configurations: (1) a zero-net-flux case with a single, radially extended set of magnetic field loops (dipole), (2) a zero-net-flux case with two radially extended sets of magnetic field loops of opposite polarity stacked vertically (quadrupole), (3) a zero-net-flux case with multiple radially concentric rings of alternating polarity (multiloop), and (4) a net-flux, vertical magnetic field configuration (vertical). In all cases, the fields are initially weak, with a gas-to-magnetic pressure ratio ≳100. Based on the results of these simulations, we find that the dipole and multiloop configurations remain thermally unstable like their α -viscosity counterpart, in our case collapsing vertically on the local thermal timescale and never fully recovering. The vertical case, on the other hand, stabilizes and remains so for the duration of our tests (many thermal timescales). The quadrupole case is intermediate, showing signs of both stability and instability. The key stabilizing factor is the ability of specific field configurations to build up and sustain strong, P mag ≳ 0.5 P tot , toroidal fields near the midplane of the disk. We discuss the reasons why certain configurations are able to do this effectively and others are not. We then compare our stable simulations to the standard Shakura–Sunyaev disk.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak

In this paper, we report a detailed experimental study of the role of ion-temperature-gradient driven turbulence in driving ion thermal transport in neutral beam injection (NBI)-heated L-mode plasmas in Experimental Advanced Superconducting Tokamak (EAST) (Wan et al 2000 Nucl. Fusion 40 1057) coupled with linear and nonlinear gyrokinetic simulations. Significant ion-scale turbulence k ⟂ < cm −1 (k ⟂ < 1.5, where k ⟂ where is the perpendicular wavenumber and ρ s is the ion gyroradius calculated using local electron temperature T e ), measured with a microwave reflectometer, is observed in the plasma core at for r/a ≈ 0.265 and 0.5 (where r is half the diameter of the closed flux surface at a given radial position, and a is half the diameter of the last closed flux surface). Local linear stability analysis with experimental equilibrium quantities at these two radial locations using the GS2 gyrokinetic code shows that the most unstable ion-scale micro-instability is the ion temperature gradient (ITG) mode. Since the computed maximum ion-scale linear growth rates are larger than the local Waltz-Miller E × B shearing rate (Waltz and Miller 1999 Phys. Plasmas 6 4265), the E × B shear is unable to suppress ion-scale turbulence, consistent with the experimental observation of ion-scale turbulence. Experimental ion and electron thermal transports, calculated with power balance analysis, are both anomalous with ion thermal transport being the dominant channel. Nonlinear gyrokinetic simulations using the electrostatic global particle-in-cell gyrokinetic tokamak simulation (GTS) code (Weixing Wang et al 2010 Phys. Plasmas 17 072511) show good quantitative agreement between predicted ion thermal transport and that from the power balance analysis. These results demonstrate for the first time the important role of ITG turbulence in driving thermal transport in NBI-heated L-mode plasmas in EAST, showing the effectiveness of GTS code in simulating electrostatic turbulence in EAST and the need of suppressing ITG turbulence to improve energy confinement in EAST.

Sun, P. J. [Chinese Academy of Sciences (CAS), Hef↗

DOE OSTI2020

The effect of initial conditions on mixing transition of the Richtmyer–Meshkov instability

We investigate the late-time Richtmyer–Meshkov instability (RMI) growth of sinuous perturbations on an air/sulphur hexafluoride interface (Atwood number, A ~ 0.67 ) subjected to a Mach 1.2 planar shock wave at Los Alamos National Laboratory's vertical shock tube facility. Interface perturbations are established using a novel membraneless technique where cross-flowing air and SF 6 separated by an oscillating splitter plate create a perturbed density interface. The interface formed has multi-modal features and residual small perturbations, however, a dominant mode is still noticeable. The late-time perturbation growths scale with $ka_0$ initial conditions (where $k$ is the wavenumber and $a_0$ is the initial amplitude of the dominant mode) as measured at the pre-shock interface. Past nonlinear models based on potential-flow theory, heuristic/interpolation approaches, Padé approximants and numerical simulations are evaluated against present experimental results. Accounting for an explicit $ka_0$ dependence in Sadot et al.'s model, we propose an empirical rational function that captures the asymptotic behaviour of perturbation growth for a broad range of initial conditions ( 0.30 $≤ ka_0 ≤$ 0.86 ). The onset of mixing transition and its initial condition dependence are investigated with respect to the minimum state criterion ($Re =$ 1.6 $\times$ 10 5 ) for unsteady flows by Zhou. Earlier mixing transitions for higher ka0 initial conditions are noted from local and global Reynolds number estimates which are corroborated by the existence of an inertial sub-range and formation of mixing regions indicating the physical significance of the minimum state criterion in RMI flows. The transition is accompanied by the increasing teapot-like appearance of joint probability density functions of $p–q$ (invariants of the reduced velocity gradient tensor), establishing the technique as a useful tool for turbulence detection in two-dimensional diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗