Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear stability analysis”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Pit rim decomposition into multiple quantum dots on surfaces of epitaxial thin films grown on pit-patterned substrates

Here, we report results of dynamical simulations according to an experimentally validated surface morphological evolution model on the formation of regular arrays of quantum dot molecules (QDMs) consisting of 1D arrays of smaller interacting quantum dots (QDs). These QD arrays form along the sides of each pit rim on the surface of a coherently strained thin film epitaxially deposited on a semiconductor substrate, the surface of which consists of a periodic pattern of inverted prismatic pits with rectangular pit openings. We find that this complex QDM pattern results from the decomposition of the pit rim from a “quantum fortress” with four elongated QDs into four 1D arrays of multiple smaller QDs arranged along each side of the pit rim. Systematic parametric analysis indicates that varying the pit opening dimensions and the pit wall inclination directly impacts the number of QDs in the resulting QDM pattern, while varying the pit depth only affects the dimensions of the QDs in the nanostructure pattern. Therefore, the number, arrangement, and sizes of QDs in the resulting pattern of QDMs on the film surface can be engineered precisely by proper tuning of the pit design parameters. Our simulation results are supported by predictions of morphological stability analysis, which explains the pit rim decomposition into multiple QDs as the outcome of a tip-splitting instability and provides a fundamental characterization of the post-instability nanostructure pattern. Our theoretical findings can play a vital role in designing optimal semiconductor surface patterns toward enabling future nanofabrication technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Prediction of the energetic particle redistribution by an improved critical gradient model and analysis of the transport threshold

Based on the theory of critical gradient model (CGM) and following the simulation method proposed by Waltz et al., a combination of TGLFEP and EPtran code is employed to predict the energetic particle (EP) transport induced by Alfvén eigenmodes (AEs). To be consistent with the experiment, recent improvements to the simulation method include consideration of threshold evolution and orbit loss due to finite orbit width. The revised CGM is applied to simulate two DIII-D experimental discharges (#142111 and #153071). It well reproduces the experimental profiles with multiple unstable AEs and large-scale EP transport. Discharge #142111 had previously been simulated using a nonlinear MHD-kinetic code MEGA with a transport mechanism based on stochasticity induced by overlapping AE. By comparing the simulated EP profiles, we find that the AE transport threshold is approximated by both the MEGA nonlinear stability threshold and the proposed CGM threshold (error <5% for single n and <17% for multiple n simulation). Both of them are larger than the linear stability threshold of the most unstable AE mode by a quantity of the order of the flux needed to sustain EP transport by the background turbulence. Additionally, we have also applied the improved CGM to simulate the α particle redistribution for a China Fusion Engineering Test Reactor steady state scenario. Because of the clear separation between the AE unstable region and the loss cone, only a moderate α particle loss of ~9.6% is predicted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Resistive wall tearing mode disruptions in DIII-D and ITER tokamaks

Disruptions are a serious problem in tokamaks, in which thermal and magnetic energy confinement is lost. This paper uses data from the DIII-D experiment, theory, and simulations to demonstrate that resistive wall tearing modes (RWTMs) produce the thermal quench (TQ) in a typical locked mode shot. Analysis of the linear RWTM dispersion relation shows the parameter dependence of the growth rate, particularly on the resistive wall time. Linear simulations of the locked mode equilibrium show that it is unstable with a resistive wall and stable with an ideally conducting wall. Nonlinear simulations demonstrate that the RWTM grows to sufficient amplitude to cause a complete thermal quench. The RWTM growth time is proportional to the thermal quench time. Here, the nonlinearly saturated RWTM magnetic perturbation amplitude agrees with experimental measurements. The onset condition is that the q = 2 rational surface is sufficiently close to the resistive wall. Collectively, this identifies the RWTM as the cause of the TQ. In ITER, RWTMs will produce long TQ times compared to present-day experiments. ITER disruptions may be significantly more benign than previously predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Identification of a window for quiescent H-mode operation in MHD stability diagram of DIII-D plasmas

A window for quiescent H-mode (QH-mode) operation in DIII-D plasmas in the diagram of pedestal magnetohydrodynamic (MHD) stability was identified through linear MHD stability analysis, accounting for the effects of plasma rotation and ion diamagnetic drift. The operation window lies between the stability boundaries of the kink/peeling mode (K/PM) identified with and without the inclusion of plasma rotation effects. In this region, the mode remains unstable unless rotation effects are considered alongside the ion diamagnetic drift, which is consistently included in the analysis. The stabilization of the MHD mode, facilitated by the coupled effects of plasma rotation and ion diamagnetic drift, plays a crucial role in widening the window, enabling the attainment of the QH-mode state. Furthermore, the suppression of edge localized modes (ELM) can be achieved by controlling the pedestal structure to ensure the plasma state remains within the operation window. The location of the operation window in the stability diagram depends on the K/PM stability properties. Therefore, optimizing conditions for QH-mode requires adjustments based on changes in stability characteristics. A pressure pedestal and its associated bootstrap current density near the last closed flux surface are advantageous for situating the plasma state within the window. However, excessive current density can trigger ELMs. This trend was confirmed through comparisons of MHD stability diagrams between QH-mode and ELMy H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Nonlinear error field response in the presence of plasma rotation and real frequencies due to favorable curvature

This works describes the interaction of an phenomenon known to arise in tokamak fusion experiments with the magnetic perturbations that can arise at the edge of such experiments. We present nonlinear NIMROD resistive MHD simulations of the response of a rotating plasma to an error field when the plasma has weakly damped linear tearing modes (TMs), stabilized by a pressure gradient and favorable curvature. The favorable curvature leads to the Glasser effect: the occurrence of real frequencies and stabilization with positive stability index Δ'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A comprehensive conductivity model for drift and micro-tearing modes

The parallel electrical conductivity is a crucial parameter in the study of the linear stability of drift-modes like the micro-tearing mode (MTM). The conductivity enters by closing the electromagnetic tearing layer equations. Recent progress in the understanding of the pedestal suggests that the MTM could play an important role in its structure and evolution. For this reason, we revisit and improve previous model conductivities. This parameter is generally derived from the linearized drift kinetic equation. In the past literature, it has been computed using either simplified collision operators or neglecting the spatial dependence away from the rational surface. Here, a fully consistent expression for the conductivity that would accurately model the pedestal has not been available. By applying a novel variational procedure and with the full Fokker Plank collision operator, including electron–electron collisions, we compute a closed expression for the parallel electrical conductivity in the form of a rational function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Consistent coupling algorithms for coupled core-edge simulations of plasma turbulence

Two schemes for coupling gyrokinetic simulations of microturbulence in tokamaks are proposed. The first scheme is based on an additive Schwarz domain decomposition. In this work, we show that, because the goal of turbulence is long-time averages of the dynamics rather than temporal accuracy, the iteration to self-consistency across domains, which is typically required by Schwarz schemes, can be avoided, thereby accelerating the computation. Second, we propose a coupling scheme that relies entirely on the addition of source terms, leaving the boundary conditions arbitrary. The practical motivations for such a scheme are discussed, and forms of the source terms that ensure consistency and stability are derived. The schemes are tested on a nonlinear, one-dimensional model problem, and the first scheme is further tested on the Hasegawa–Wakatani model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Microtearing instabilities and electron thermal transport in low and high collisionality NSTX discharges

Microtearing mode (MTM) real frequency, growth rate, magnetic fluctuation amplitude and resulting electron thermal transport are studied in systematic NSTX scans of relevant plasma parameters. Here, the dependency of the MTM real frequency and growth rate on plasma parameters, suitable for low and high collision NSTX discharges, is obtained by using the reduced MTM transport model. The plasma parameter dependencies are compared and found to be consistent with the results obtained from MTM using the Gyrokinetic GYRO code. The scaling trend of collision frequency and plasma beta is found to be consistent with the global energy confinement trend observed in the NSTX experiment. The strength of the magnetic fluctuation is found to be consistent with the gyrokinetic estimate. In earlier studies, it was found that the version of the Multi-Mode (MM) anomalous transport model, which did not contain the effect of MTMs, provided an appropriate description of the electron temperature profiles in standard tokamak discharges and not in spherical tokamaks. When the MM model, which involves transport associated with MTMs, is incorporated in the TRANSP code and is used in the study of electron thermal transport in NSTX discharges, it is observed that the agreement with the experimental electron temperature profile is substantially improved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Achieving record hot spot energies with large HDC implosions on NIF in HYBRID-E

HYBRID-E is an inertial confinement fusion implosion design that increases energy coupled to the hot spot by increasing the capsule scale in cylindrical hohlraums while operating within the current experimental limits of the National Ignition Facility. HYBRID-E reduces the hohlraum scale at a fixed capsule size compared to previous HYBRID designs, thereby increasing the hohlraum efficiency and energy coupled to the capsule, and uses the cross-beam energy transfer (CBET) to control the implosion symmetry by operating the inner (23° and 30°) and outer (44° and 50°) laser beams at different wavelengths (Δλ> 0). Small case to capsule ratio designs can suffer from insufficient drive at the waist of the hohlraum. We show that only a small amount of wavelength separation between the inner and outer beams (Δλ1-2 Å) is required to control the symmetry in low-gas-filled hohlraums (0.3 mg/cm3 He) with enough drive at the waist of the hohlraum to symmetrically drive capsules 1180 μm in outer radius. This campaign is the first to use the CBET to control the symmetry in 0.3 mg/cm3 He-filled hohlraums, the lowest gas fill density yet fielded with Δλ> 0. We find a stronger sensitivity of hot spot P2 in μm per Angstrom (40–50 μm/Å wavelength separation) than observed in high-gas-filled hohlraums and previous longer pulse designs that used a hohlraum gas fill density of 0.6 mg/cm3. There is currently no indication of transfer roll-off with increasing Δλ, indicating that even longer pulses or larger capsules could be driven using the CBET in cylindrical hohlraums. We show that the radiation flux symmetry is well controlled during the foot of the pulse, and that the entire implosion can be tuned symmetrically in the presence of the CBET in this system, with low levels of laser backscatter out of the hohlraum and low levels of hot electron production from intense laser–plasma interactions. Radiation hydrodynamic simulations can accurately represent the early shock symmetry and be used as a design tool, but cannot predict the late-time radiation flux symmetry during the peak of the pulse, and semi-empirical models are used to design the experiments. Deuterium–tritium (DT)-layered tests of 1100 μm inner radius implosions showed performance close to expectations from simulations at velocities up to ~360 km/s, and record yields at this velocity, when increasing the DT fuel layer thickness to mitigate hydrodynamic mixing of the ablator into the hot spot as a result of defects in the ablator. However, when the implosion velocity was increased, mixing due to these defects impacted performance. The ratio of measured to simulated yield for these experiments was directly correlated with the level of observed mixing. These simulations suggest that reducing the mixing, e.g., by improving the capsule defects, could result in higher performance. In addition, future experiments are planned to reduce the coast time at this scale, delay between the peak compression and the end of the laser, to increase the hot spot convergence and pressure. To reduce the coast time by several hundred ps compared to the 1100 μm inner radius implosions, HYBRID-E has also fielded 1050 μm inner radius capsules, which resulted in higher hot spot pressure and a fusion energy yield of ~170 kJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Resonant amplification of enzymatic chemical oscillations by oscillating flow

Using theory and simulation, we analyzed the resonant amplification of chemical oscillations that occur due to externally imposed oscillatory fluid flows. The chemical reactions are promoted by two enzyme-coated patches located sequentially on the inner surface of a pipe that transports the enclosed chemical solution. In the case of diffusion-limited systems, the period of oscillations in chemical reaction networks is determined by the rate of the chemical transport, which is diffusive in nature and, therefore, can be effectively accelerated by the imposed fluid flows. In this work, we first identify the natural frequencies of the chemical oscillations in the unperturbed reaction–diffusion system and, then, use the frequencies as a forcing input to drive the system to resonance. We demonstrate that flow-induced resonance can be used to amplify the amplitude of the chemical oscillations and to synchronize their frequency to the external forcing. In particular, we show that even 10% perturbations in the flow velocities can double the amplitude of the resulting chemical oscillations. Particularly, effective control can be achieved for the two-step chemical reactions where during the first half-period, the fluid flow accelerates the chemical flux toward the second catalytic patch, while during the second half-period, the flow amplifies the flux to the first patch. The results can provide design rules for regulating the dynamics of coupled reaction–diffusion processes and can facilitate the development of chemical reaction networks that act as chemical clocks.

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

Integrated modeling of advanced tokamaks for tearing mode avoidance on DIII-D

Understanding the physics of the onset and evolution of tearing modes (TMs) in tokamak plasmas is crucial for high-performance steady-state operations. The onset condition of (m,n) = (2,1) tearing stability, $Δ′$ > $Δ_c^′$> 0, has been studied in DIII-D steady-state hybrid discharges with accurate equilibrium reconstruction and well-measured plasma profiles. Here, m is the poloidal mode number and n is the toroidal mode number. The onset of n = 1 TMs is carefully identified using the Mirnov array to determine in-phase signals for the onset time and using electron cyclotron emission measurements to trace back from the saturated island width in the experiments for both initial island width and its growth rate. The tearing stability index $Δ′$ at TM onset is calculated in multiple ways, using the modified Rutherford equation and using the PEST-III or resistive DCON magnetohydrodynamic codes, and the results are in reasonable agreement with each other. The calculated $Δ′$ is compared with an analytical formula of the tearing stability threshold $Δ_c^′$ to determine the condition for the mode onset, $Δ′$ > $Δ_c^′$. Finally, the local gradient effect of the plasma current and the pressure profiles on the tearing stability is investigated in terms of the difference between $Δ′$ and $Δ_c^′$ through a series of modeled equilibria with the parametric variations. A database of experimental mode instability onsets is compatible with the simulated regions of instability obtained from the parametric variations.

Linear stability analysis↗
DOE OSTI2025

Prevention of resistive wall tearing mode major disruptions with feedback

Resistive wall tearing modes (RWTMs) can cause major disruptions. A signature of RWTMs is that the rational surface is sufficiently close to the wall to interact with it. For (m,n) = (2,1) modes, a RWTM requires normalized minor radius of the rational surface ρ q2 ≥ 0.75, which can also be expressed as q 75 ≤ 2. Major disruptions can occur when the criterion is satisfied. This is confirmed in simulations and theory and in a DIII-D locked mode disruption database. The q 75 < 2 criterion is valid at high β as well as at low β. A very important feature of RWTMs is that they can be feedback stabilized. If the ρ q2 criterion is not satisfied, or if the wall is ideally conducting, then the mode does not produce a major disruption, although it can produce a minor disruption. Feedback, or rotation of the mode at the wall by complex feedback, can emulate an ideal wall, preventing major disruptions. The ρ q2 criterion depends weakly on the wall radius. A simple geometric model of its dependence on wall radius is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

On resistive interchange, double tearing, and resonant and non-resonant infernal modes in spherical tokamak negative central shear discharges

We examine the linear and nonlinear stability of a sequence of reversed shear NSTX equilibria with the same toroidal current density and pressure profile but with different toroidal field strengths. All equilibria have two q = 2 surfaces and have q > 1 everywhere. For this sequence, which all have β P ∼ 0.5, as the minimum value of q decreases below about 1.4, the localized resistive interchange criterion is strongly violated, and the unstable mode with toroidal mode number n=1 changes its dominant poloidal mode number m from a (m,n)=(2,1) double tearing mode to a non-resonant (1,1) infernal mode. This (1,1) mode nonlinearly flattens both the current density and the pressure near the magnetic axis. The higher-n unstable modes are all resistive infernal modes localized around the surface where q has a minimum and the magnetic shear vanishes. They generally saturate nonlinearly at low amplitude but can cause magnetic surface breakup near the axis.

Spitzer resistivity↗
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2026

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Analysis of the linear and nonlinear stability of Alfven eigenmodes and fish-bones in JET DT discharges: mode identification and shear flows generation

The plasma in future nuclear fusion reactors will be heated by neutral beam injectors (NBIs) and high frequency electromagnetic waves as well as fusion born alpha particles. Energetic particles (EPs), with energies up to two orders of magnitude larger than the thermal plasma, can trigger EP driven modes and induce harmful EP losses, reducing the plasma heating efficiency and the economical viability of the reactor. The present study is dedicated to analyze the Alfven Eigenmode (AE) activity in JET D–T discharges, the closest experiment to reactor-like operation performed until now. There, EP driven modes are induced by the combined effect of tangential NBIs and ion cyclotron resonance heating (ICRH) driven EP. Linear and nonlinear simulations are performed with the gyro-fluid FAR3d code to analyze the AE activity observed in the discharge 99896. The linear simulations reproduce the unstable n = 3 to 5 toroidal AEs (TAE) at the inner plasma region observed in the experiment, triggered by highly energetic passing deuterium populations injected by the tangential NBIs, further accelerated by the effect of the ICRH up to 1 MeV. In addition, fish-bones triggered by energetic trapped hydrogen induced by the ICRH are also reproduced. On the other hand, the alpha particles density is too small to destabilize AEs in the experiment. Nonetheless, increasing artificially the alpha density by one order of magnitude, an n = 1 beta induced AE can be destabilized in the inner plasma region. Nonlinear simulations indicate the generation of zonal structures during the AE/fish-bone saturation phase. TAE and fish-bones causes a rather weak increase of the passing D and trapped H EP (around 2%), respectively. Shear flows and zonal currents are generated during the saturation of TAE and fish-bones. Nonlinear simulations performed for D–T and pure deuterium thermal plasma indicate AE/fish-bone activity is weaker and shear flows are less intense in the pure deuterium case, trends consistent with the experimental observations that also indicates a deterioration of the thermal plasma confinement. Therefore, both numerical studies and experimental evidence indicate the generation of shear flows by AE/fish-bones could be connected with an improvement of the thermal plasma confinement.

AE↗