Search NASA⌕ Search

SEARCH · Search NASA

Results for “edge 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Study of the electrothermal and MHD instabilities in exploding cylindrical foil liner

Abstract An experimental and numerical study of the plasma instabilities in an electrical exploding cylindrical Al liner is reported. The Al liner 3 mm in diameter and 10 μ m in thickness is exploded at the Qin-1 facility (450 ns, 400 kA). Various diagnostics, such as multi-frame laser shadowgraphy, an x-ray framing camera, and an x-ray backlighting system are developed. The different imaging systems are sensitive to plasma of different areal densities based on the comparison between the experiments and simulation, which reveal the dynamics of the exploding liner in more detail. The laser shadow images show the low-density plasma (∼1–2 × 10 −4 g cm −2 ) at the edge of the liner, and both the amplitude and wavelength of the plasma instabilities increase over time, which are considered to be magnetohydrodynamic (MHD) instabilities rather than electrothermal instabilities. During the ablation of the liner, quasi-periodic azimuthally correlated striations are directly observed in extreme ultraviolet (EUV) self-emission images. Meanwhile, the vertical filaments, which are electrothermal instabilities for plasma under the condition of ∂ η /∂ T < 0, are also observed in EUV self-emission images. The x-ray backlighting images of the exploding liner are obtained by placing an X -pinch load on the current-return path to serve as an x-ray point source (∼1 ns, ∼10 μ m). The x-ray backlighting results show the behavior of the high-density plasma (∼1.89 × 10 −3 g cm −2 ), which includes the transition from electrothermal to MHD instabilities. Finally, we realized a 2D MHD simulation of the exploding liner under experimental conditions, which shows good agreement with the results of the experimental perturbation.

Physics↗
DOE OSTI2024

Active control of Alfvén eigenmodes by external magnetic perturbations with different spatial spectra

Alfvén eigenmodes have been suppressed and excited in tokamak plasmas by (just) modifying the poloidal spectra of externally applied static magnetic perturbations. This effect is observed experimentally when toroidal spectra of n = 2, n = 4 as well as a mixed spectrum of n = 2 and n = 4 is applied. Under the n = 2 magnetic perturbations, the modes are excited or suppressed by modifying the coil phasing between the upper and the lower set of coils. Regardless of the absolute rotation, an even parity for the n = 4 perturbation is observed to reduce the amplitude of the Alfvénic instabilities, while an odd parity amplifies it. To combine the stabilizing (and destabilizing) effect of n = 2 and n = 4, a mixed spectrum is applied, finding similar reduction (and amplification) trends. However, the impact on the mode amplitude is more subtle, due to the reduced coil current required for a mixed spectrum. The signal level on the fast-ion loss detector is sensitive to the applied poloidal spectrum, which is consistent with Hamiltonian full-orbit modelling of an edge resonant transport layer activated by the 3D perturbative fields. An internal redistribution of the fast-ion population is induced, modifying the phase-space gradients driving the Alfvénic instabilities, and ultimately determining their existence. The calculated edge resonant layers for both n = 2 and n = 4 toroidal spectra are consistent with the observed suppressed and excited phases. Moreover, hybrid kinetic-magnetohydrodynamic (MHD) simulations reveal that this edge resonant transport layer overlaps in phase-space with the population responsible for the fast-ion drive. The results presented here may help to control fast-ion driven Alfvénic instabilities in future burning plasmas with a significant fusion born alpha particle population.

Alfvén waves↗
DOE OSTI2024

Krash-ESL-Final report

During the period of the current Grant, the effort of the UCSD researchers was focused on: i) The application of the gyrokinetic code COGENT code to study the kinetic drift wave instability and its impact on blob dynamics, and ii) The modeling of edge plasma transport in collisional and semi-collisional regimes, and iii) Impurity transport in edge plasma and the incorporation of plasma-impurity collision operator (hereafter unlike particle collision operator) into COGENT. Over the years four UCSD graduate students (J. Angus, W. Lee, Y. Zhang, and A. Knyazev), funded under this Grant, contributed to the COGENT development/utilization and defended their PhD theses, which are based, partly, on their research related to the COGENT. Three of them (J. Angus, W. Lee, and Y. Zhang) have received “The Best MAE Graduate Student of the Year” awards. Whereas W. Lee after defending his PhD thesis returned back to his home country (South Korea), J. Angus, and Y. Zhang are continuing their research in plasma physics in the LLNL and LANL, whereas A. Knyazev currently is a postdoctoral fellow at Columbia University working on the theoretical aspects of advanced stellarators. The results obtained in our studies were disseminated in many papers published in first-class refereed scientific journals and presented on prestigious scientific meetings including the Annual APS DPP Meeting, Sherwood Fusion Theory Conference, IAEA Fusion Energy Conferences, International Conferences on Plasma-Surface Interactions in Controlled Fusion Devices, International Workshops on Edge Plasma Theory in Fusion Devices, etc. In what follows we highlight the major accomplishments of the UCSD Team in the main areas of our research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Tearing Stable Stationary ITER Baseline Operation in DIII-D

We present a tearing stable, stationary and reproducible operating solution for DIII-D plasmas characterized by the ITER Baseline Scenario normalized parameter set and shape. In these plasmas low differential rotation (Δf1,2) between the core and edge is identified as the single root cause of the onset of the primary limiting instability, the 2,1 Neoclassical Tearing Mode (NTM), while a group of other examined parameters have a much weaker impact on the stability. Explanation is offered by drift kinetic simulations, which show a six-fold reduction in the NTM onset threshold due to diminished stabilizing ion polarization currents when the seed island drift frequency ceases in the local plasma frame. Based on the experimental observation and the explanation provided by the theory, a new method is developed to sustain Δf1,2 by modifying the edge neoclassical potential through neutral gas fueling. This method enables stationary ITER baseline operation free of disruptive 2,1 tearing modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Observation of impurity accumulation and its compatibility with high plasma performance in W7-X

At the W7-X stellarator, the bolometer system has measured an intensive radiation zone in the inner plasma region (at a normalized radius $ρ$ ~ 0.3–0.4) in the hydrogen plasma generated by electron cyclotron resonance heating; it differs from the normal plasma radiation distribution with an edge-localized emission zone. Spectroscopic diagnostics have recorded high-Z elements such as iron. This phenomenon happens in the plasma phases after gas supply turn-off, which results in all impurity relevant diagnostic signals increasing for several seconds. Despite the enhancement of the core radiation, the plasma energy confinement is improved. A transport analysis shows that this impurity radiation behavior is associated with a low diffusion coefficient ($D$~ 0.02 m 2 s -1 ) and a reversal of the convection around the radial position of the emission peak, which, under normal conditions, separates the zones of outward convection in the central (|$V$| ~ 0.1 m s -1 ) and inward convection in the outer region (|$V$| ~ 0.3 m s -1 ). An impurity accumulation around this radial position has been identified. The transport coefficients obtained are comparable with the theoretical predictions of collisional impurity transport. In the plasma phases studied, both impurity and energy confinement are enhanced. The mechanism responsible for the improvement is believed to be a reduction of micro-instabilities associated with the observed steepening of the density profile, initiated by a low edge plasma density (<1.0 × 10 19 m -3 ) after switching off the gas fueling. The normalized temperature and density gradients fulfil the condition for the suppression of ITG turbulence.

suppression of ITG turbulence↗
DOE OSTI2026

Effects of density turbulence on helicon wave propagation in the core plasmas

Radio frequency wave propagation can be significantly affected by density irregularities, such as filaments in the scrape-off layer or instabilities in the core plasma. In this study, we examine the impact of edge turbulence on helicon wave propagation using the Petra-M simulation tool. To analyze the effect of edge turbulence, we utilize a realistic background plasma derived from XGC simulations, which includes spatial density fluctuations at the edge. This focus is particularly relevant for the DIII-D configuration, characterized by edge density turbulence of the core plasmas in a wide pedestal QH-mode. We focus on helicon wave propagation in the core plasma since the slow mode cannot propagate into the core plasma when the density is higher than the lower hybrid resonance, where a wide pedestal QH-mode occurs. The simulation results indicate that edge density fluctuations have a substantial impact on helicon wave coupling due to scattering in the core. We specifically demonstrate that the toroidal mode number and the level of density fluctuations are significant factors influencing scattering. While we have minimized the excitation of slow modes from the antenna, we still observe mode-converted slow modes resulting from the incoming helicon waves. The insights derived from these simulations will inform upcoming tokamak experiments regarding helicon antenna coupling in long pulse scenarios.

Kim, Eun-Hwa [Princeton Plasma Physics Laboratory ↗
DOE OSTI2021

Investigations of LHW-plasma coupling and current drive at high density related to H-mode experiments in EAST

Two important issues in achieving lower hybrid current drive (LHCD) high confinement plasma in EAST are to improve lower hybrid wave (LHW)-plasma coupling and to drive the plasma current at a high density. Studies in different configurations with different directions of toroidal magnetic field (Bt) show that the density near the antenna is affected by both the radial electric field induced by plasma without a LHW (Er-plasma) in the scrape off layer (SOL), and the radial electric field induced by LHW power (Er-LH) near the grill. Investigations indicate that Er× Btin the SOL leads to a different effect of configuration on the LHW-plasma coupling and Er-LH× Btaccounts for the asymmetric density behaviour in the SOL observed in the experiments, where Eris the total radial electric field in the SOL. Modelling of parametric instability (PI), collisional absorption (CA) and scattering from density fluctuations (SDF) in the edge region, performed considering the parameters of high density LHCD experiments in EAST, has shown that these mechanisms could be responsible for the low current drive (CD) efficiency at high density. Radiofrequency probe spectra, useful for documenting PI occurrence, show sidebands whose amplitude in the case of the lithiated vacuum chamber is smaller than in the case of poor lithiation, consistently with growth rates from PI modeling of the respective reference discharges. Since strong lithiation is also expected to diminish the parasitic effect on the LHCD of the remaining possible mechanisms, this appears to be a useful method for improving LHCD efficiency at a high density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Energizing charged particles by an orbit instability in a slowly rotating magnetic field

In this work, the stability of charged particle motion in a uniform magnetic field with an added spatially uniform transverse rotating magnetic field (RMF) is studied analytically. It is found that the stability diagram of a single particle's orbit depends critically on the chosen boundary conditions. We show that for many boundary conditions and wide regions in the parameter space, RMFs oscillating far below the cyclotron frequency can cause linear instabilities in the motion which break μ invariance and energize particles. Such energization may appear at odds with the adiabatic invariance of μ; however, adiabatic invariance is an asymptotic result and does not preclude such heating by magnetic fields oscillating at slow frequencies. This mechanism may contribute to heating in the edge plasma of field-reversed configurations (FRCs) in rotamak-FRC experiments. Furthermore, these RMF-driven instabilities may significantly enhance azimuthal current drive during the formation of FRCs in such devices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Broadening of the Divertor Heat Flux Profile in High Confinement Tokamak Fusion Plasmas with Edge Pedestals Limited by Turbulence in DIII-D

Multimachine empirical scaling predicts an extremely narrow heat exhaust layer in future high magnetic field tokamaks, producing high power densities that require mitigation. Here, in the experiments presented, the width of this exhaust layer is nearly doubled using actuators to increase turbulent transport in the plasma edge. This is achieved in low collisionality, high confinement edge pedestals with their gradients limited by turbulent transport instead of large-scale, coherent instabilities. The exhaust heat flux profile width and divertor leg diffusive spreading both double as a high frequency band of turbulent fluctuations propagating in the electron diamagnetic direction doubles in amplitude. The results are quantitatively reproduced in electromagnetic XGC particle-in-cell simulations which show the heat flux carried by electrons emerges to broaden the heat flux profile, directly supported by Langmuir probe measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Impurity radiation seeding of neoclassical tearing mode growth

Abstract The physics of neoclassical tearing modes (NTMs) is of great concern to tokamak plasma stability and performance, especially in the burning plasma regime. Whereas in many situations the seed events can be clearly identified, such as sawteeth and edge localized modes, the potential seeding mechanism of NTMs due to the resistive tearing instability driven by the impurity radiation cooling still needs more study. Recent NIMROD simulations have demonstrated that local impurity radiation cooling can drive the seed island growth and trigger the subsequent onset of NTM instability. The seed island is mainly driven by the local helical perturbation of the diamagnetic current induced by the perturbed pressure gradient as a result of the impurity radiative cooling on the rational surface. A heuristic closure for the neoclassical viscosity is adopted, and the seed island is further driven by the perturbed bootstrap current induced from the neoclassical electron viscous stress in the extended Ohm’s law. The growth rate of the NTM in simulations is found proportional to the electron neoclassical viscosity, and a theoretical neoclassical driving term is adopted to account for the nonlinear neoclassical island growth in the simulations.

Zeng, Shiyong (ORCID:000000017098818X)↗
DOE OSTI2021

Regular and in-plane skyrmions and antiskyrmions from boundary instabilities

We formulate a theory of skyrmion and antiskyrmion generation using magnetic field and charge current pulses. We show that the topological defect can be created at an edge of a system with Dzyaloshinskii-Moriya interaction (DMI) as well as at a boundary between regions with different DMI. We consider both perpendicular and in-plane (also known as magnetic bimerons) versions of skyrmions and antiskyrmions. We show that the magnetization twist in the vicinity of an edge or a boundary is described by a kink solution, the presence of which can instigate the generation of topological defects. We study the collective excitations of magnetization analytically and numerically, and demonstrate that under application of magnetic field and charge current pulses the magnon modes localized near boundaries can develop instabilities leading to the formation of skyrmions or antiskyrmions. As a result, due to the skyrmion and antiskyrmion Hall effects, a properly chosen current direction can push the topological defects away from the boundary, thus facilitating their generation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Turbulence simulations of drift-Alfvén wave effects on edge-localized modes and divertor target heat flux width

A series of BOUT++ turbulence simulations are performed for two Pre-Fusion Power Operation (PFPO-1,2) phases of the ITER Research plan proclaimed in 2019. Linear simulation results show that PFPO-1 is unstable to Peeling-ballooning modes (PB), while the PFPO-2 is unstable to the coupling of PB and Drift-Alfvén wave (DAW) instabilities. The linear results are qualitatively consistent with the dispersion relation of DAW. To study the influence of the DAW on the edge-localized mode (ELM) crash, the PFPO-2 is adopted in this section. Different from the grassy ELM in [Li et al 2022 Nucl. Fusion 62, 096030], nonlinear simulations show that the ELM size of PFPO-2 is almost one-third of the grassy ELM, representing a distinct small ELM. However, simulations then show that if the PB instability is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crash disappears, which is in accord with the theory in [Xu et al 2010 Physical Review Letters, Vol. 105, 175005] and the results in [Li et al 2022 Nucl. Fusion 62, 096030], confirming that the PB instability is a necessary condition for an ELM crash. Furthermore, removing the DAW drive also suppresses ELM crashes, implying that PB instability is necessary but insufficient for PFPO-2 ELM and that DAW could amplify PB-driven turbulence. In addition, by integrating results of both PFPO phases by BOUT++ turbulence analysis with heat flux width (λ q ) and electron thermal diffusivity (χ e ), these simulations establish a predictive framework for heat flux width and its dominant mechanisms, thereby providing a comprehensive scheme for the prediction of the heat flux width. Moreover, simulations indicate that DAW driving can increase the transport coefficient by enhancing the turbulent transport, leading to a broadened heat flux width once the transport coefficient exceeds

Drift-Alfvén wave↗
DOE OSTI2026

Density limit in peeling-limited pedestals at and above the Greenwald value in DIII-D high poloidal beta plasmas

High pressure, peeling limited pedestals with pedestal normalized beta β N,ped >2 and pedestal top density n ped at or above the Greenwald density nG have been achieved in DIII-D high poloidal beta plasmas, with high global normalized beta β N >3 and energy confinement H 98 ~1.2-1.7. Higher β N allows higher pedestal density above the Greenwald value and higher pedestal pressure, even with a low injected torque. MHD modeling confirms that the experimental profiles lie near the peeling-mode unstable boundary with high normalized pressure gradient βMHD and high edge current density. Experimental analysis and stability calculations indicate that the high poloidal beta with strong Shafranov shift, high βMHD and weak/negative magnetic shear improves the pedestal stability by decoupling the peeling and ballooning modes and stabilizing the ballooning modes, thus facilitating access to the second stable region of peeling-ballooning mode. The access to the second stable peeling-ballooning stability region opens Super-H-like channels without extremely strong shaping or strong torque injection. The high-pressure peeling pedestal allows the pedestal density to go beyond the Greenwald limit with strong ExB shear maintained: pedestal pressure increases with pedestal density even when n ped >n G , until reaching the ideal MHD instability boundary, where giant ELMs occur. The giant ELMs are dominated by a strong n=1 component and cause a large reduction of the edge pressure, but a negligible change of the core pressure, consistent with kink/peeling-mode induced instability. The pedestal recovers from the collapse and typically sustains a high baseline density, around the Greenwald limit, during the whole discharge duration. Experiments also found that internal transport barriers and n ped ~0.9nG, peeling limited pedestals could be simultaneously achieved in high β N plasmas, while an internal feedback between ITB strength and pedestal performance is found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Lithium splashing from flowing liquid lithium limiter and its effect on high confinement plasma performance in EAST tokamak

To study the potential applications of liquid lithium (Li) as plasma facing material, the third generation of flowing liquid Li (FLiLi) limiter has been designed and successfully tested in EAST tokamak H-mode discharges in the 2018 run campaign. Here, the stability of liquid Li on FLiLi surface, and Li splashing and its effect on the plasma performance are investigated in detail. Li splashing was observed in L-mode and H-mode discharges due to Rayleigh-Taylor (R-T) instability, resulting from electromagnetic forces J → x B → . Li splashing in low-parameter plasma decreased both core and edge electron temperatures, resulting in the degradation of plasma energy confinement and even H-L back transition, however, it almost did not influence plasma confinement performance in high-parameter plasma. These findings indicate that high-parameter plasma has higher tolerance to Li impurities than low-parameter plasma. By contrast, tungsten splashing from the melted layer on the upper divertor was also induced by the R-T instability but generally resulted in plasma disruption, which suggests lower tolerance to tungsten impurities in EAST tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Alcove formation in dissolving cliffs driven by density inversion instability

We demonstrate conditions that give rise to cave-like features commonly found in dissolving cliffsides with a minimal two-phase physical model. Alcoves that are wider at the top and tapered at the bottom, with sharp-edged ceilings and sloping floors, are shown to develop on vertical solid surfaces dissolving in aqueous solvents. As evident from descending plumes, sufficiently large indentations evolve into alcoves as a result of the faster dissolution of the ceiling due to a solutal Rayleigh–Bénard density inversion instability. In contrast, defects of size below the boundary layer thickness set by the critical Rayleigh number smooth out, leading to stable planar interfaces. Furthermore, the ceiling recession rate and the alcove opening area evolution are shown to be given to first-order by the critical Rayleigh number. By tracking passive tracers in the fluid phase, we show that the alcoves are shaped by the detachment of the boundary layer flow and the appearance of a pinned vortex at the leading edge of the indentations. The attached boundary layer past the developing alcove is then found to lead to rounding of the other sides and the gradual sloping of the floor.

58 GEOSCIENCES↗
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 OSTI2020

Double X/Peanut structures in barred galaxies – insights from an N -body simulation

ABSTRACT Boxy, peanut-, or X-shaped ‘bulges’ are observed in a large fraction of barred galaxies viewed in, or close to, edge-on projection, as well as in the Milky Way. They are the product of dynamical instabilities occurring in stellar bars, which cause the latter to buckle and thicken vertically. Recent studies have found nearby galaxies that harbour two such features arising at different radial scales, in a nested configuration. In this paper, we explore the formation of such double peanuts, using a collisionless N-body simulation of a pure disc evolving in isolation within a live dark matter halo, which we analyse in a completely analogous way to observations of real galaxies. In the simulation, we find a stable double configuration consisting of two X/peanut structures associated with the same galactic bar – rotating with the same pattern speed – but with different morphology, formation time, and evolution. The inner, conventional peanut-shaped structure forms early via the buckling of the bar, and experiences little evolution once it stabilizes. This feature is consistent in terms of size, strength, and morphology, with peanut structures observed in nearby galaxies. The outer structure, however, displays a strong X, or ‘bow-tie’, morphology. It forms just after the inner peanut, and gradually extends in time (within 1–1.5 Gyr) to almost the end of the bar, a radial scale where ansae occur. We conclude that, although both structures form, and are dynamically coupled to, the same bar, they are supported by inherently different mechanisms.

Ciambur, Bogdan C.↗