SEARCH · Engineering Papers

Results for “curvature confinement”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2002

First Results of the Gasdynamic Mirror Fusion Propulsion Experiment

An experimental Gasdynamic Mirror or GDM device has been constructed at the NASA Marshall Space Flight Center to provide an initial assessment of the applicability of this technology for propulsion systems. This paper presents the first experimental results obtained from the machine and an analysis of the types of plasma instabilities likely to be encountered. It is intended that this device operate at higher plasma densities and with much larger L/D ratios than previous mirror machines. The high L/D ratio minimizes to a large extent certain magnetic curvature effects which lead to plasma instabilities causing a loss of plasma confinement. The high plasma density results in the plasma behaving much more like a conventional fluid with a mean free path shorter than the length of the device. This characteristic helps reduce problems associated with "loss cone" microinstabilities. The device has been constructed to allow a considerable degree of flexibility in its configuration thus permitting the experiment to grow over time without necessitating a great deal of additional fabrication.

Emrich, William J., Jr.↗

NASA NTRS2003

Enhanced gravitropism of roots with a disrupted cap actin cytoskeleton

The actin cytoskeleton has been proposed to be a major player in plant gravitropism. However, understanding the role of actin in this process is far from complete. To address this problem, we conducted an analysis of the effect of Latrunculin B (Lat B), a potent actin-disrupting drug, on root gravitropism using various parameters that included detailed curvature kinetics, estimation of gravitropic sensitivity, and monitoring of curvature development after extended clinorotation. Lat B treatment resulted in a promotion of root curvature after a 90 degrees reorientation in three plant species tested. More significantly, the sensitivity of maize (Zea mays) roots to gravity was enhanced after actin disruption, as determined from a comparison of presentation time of Lat B-treated versus untreated roots. A short 10-min gravistimulus followed by extended rotation on a 1-rpm clinostat resulted in extensive gravitropic responses, manifested as curvature that often exceeded 90 degrees. Application of Lat B to the cap or elongation zone of maize roots resulted in the disruption of the actin cytoskeleton, which was confined to the area of localized Lat B application. Only roots with Lat B applied to the cap displayed the strong curvature responses after extended clinorotation. Our study demonstrates that disrupting the actin cytoskeleton in the cap leads to the persistence of a signal established by a previous gravistimulus. Therefore, actin could function in root gravitropism by providing a mechanism to regulate the proliferation of a gravitropic signal originating from the cap to allow the root to attain its correct orientation or set point angle.

NASA Discipline Plant Biology↗

DOE OSTI2024

Nonlinear gyrokinetic modelling of high confinement negative triangularity plasmas

Abstract Nonlinear gyrokinetic simulations correctly predict particle as well as ion and electron energy fluxes of high confinement plasmas with a negative triangularity cross sectional shape, showing that core transport in these plasmas is well described by standard gyrokinetic models. Experimentally inferred power balance fluxes are mostly reproduced within one standard deviation across a wide portion of the minor radius. Experimental conditions are reproduced by ion scale simulations, without the need to include density and temperature profile curvature effects. The experimental case is used as baseline to predict that the non-dimensional confinement scaling in negative triangularity plasmas increases strongly with plasma current while slightly degrading at increasing normalized pressure and decreasing collisionality. Recent experiments showed that low toroidal rotation negatively impacts confinement; consistent with the experiment, simulations predict that low rotational shear significantly affects confinement unless the plasma effective charge is maintained above a minimum level. Core confinement is predicted to significantly degrade in low aspect ratio devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1983

Aerothermal modeling program, phase 1

The combustor performance submodels for complex flows are evaluated. The benchmark test cases for complex nonswirling flows are identified and analyzed. The introduction of swirl into the flow creates much faster mixing, caused by radial pressure gradients and increase in turbulence generation. These phenomena are more difficult to predict than the effects due to geometrical streamline curvatures, like the curved duct, and sudden expansion. Flow fields with swirl, both confined and unconfined are studied. The role of the dilution zone to achieve the turbine inlet radial profile plays an important part, therefore temperature field measurements were made in several idealized dilution zone configurations.

Srinivasan, R.↗

DOE OSTI2026

Omnigenous umbilic stellarators

To better understand the dependence of the magnetic field structure in the plasma edge on the plasma boundary shape, in the context of X-point and island divertor designs, we define and develop a class of stellarators called umbilic stellarators. These equilibria are characterised by a single continuous high-curvature edge on the plasma boundary that goes around multiple times toroidally before meeting itself. We develop a technique that allows us to simultaneously optimise the plasma boundary along with a curve lying on the boundary on which we impose a high curvature while imposing omnigenity – a property of the magnetic field that ensures trapped particle confinement throughout the plasma volume. We find that umbilic stellarators naturally tend to favour piecewise omnigenity instead of omnigenity with a specific helicity. After generating omnigenous umbilic stellarators, we design coil sets for some of them and explore the fieldline structure in the edge and its sensitivity to small fluctuations in the plasma. Finally, using single-stage optimisation, we simultaneously modify the plasma and coil shape and propose an experiment to modify an existing tokamak to a finite-𝛽 stellarator using this technique and explore a potentially simpler way to convert a limited tokamak into a diverted stellarator.

fusion plasma↗

NASA NTRS2012

Equilibrium Limit of Boundary Scattering in Carbon Nanostructures: Molecular Dynamics Calculations of Thermal Transport

It is widely known that graphene and many of its derivative nanostructures have exceedingly high reported thermal conductivities (up to 4000 W/mK at 300 K). Such attractive thermal properties beg the use of these structures in practical devices; however, to implement these materials while preserving transport quality, the influence of structure on thermal conductivity should be thoroughly understood. For graphene nanostructures, having average phonon mean free paths on the order of one micron, a primary concern is how size influences the potential for heat conduction. To investigate this, we employ a novel technique to evaluate the lattice thermal conductivity from the Green-Kubo relations and equilibrium molecular dynamics in systems where phonon-boundary scattering dominates heat flow. Specifically, the thermal conductivities of graphene nanoribbons and carbon nanotubes are calculated in sizes up to 3 microns, and the relative influence of boundary scattering on thermal transport is determined to be dominant at sizes less than 1 micron, after which the thermal transport largely depends on the quality of the nanostructure interface. The method is also extended to carbon nanostructures (fullerenes) where phonon confinement, as opposed to boundary scattering, dominates, and general trends related to the influence of curvature on thermal transport in these materials are discussed.

Haskins, Justin↗

NASA NTRS1997

Magnetophoretic induction of curvature in coleoptiles and hypocotyls

Coleoptiles of barley (Hordeum vulgare) were positioned in a high gradient magnetic field (HGMF, dynamic factor gradient of H(2)/2 of 10(9)-10(10) Oe2 cm-1), generated by a ferromagnetic wedge in a uniform magnetic field and rotated on a 1 rpm clinostat. After 4 h 90% of coleoptiles had curved toward the HGMF. The cells affected by HGMF showed clear intracellular displacement of amyloplasts. Coleoptiles in a magnetic field next to a non-ferromagnetic wedge showed no preferential curvature. The small size of the area of nonuniformity of the HGMF allowed mapping of the sensitivity of the coleoptiles by varying the initial position of the wedge relative to the coleoptile apex. When the ferromagnetic wedge was placed 1 mm below the coleoptile tip only 58% of the coleoptiles curved toward the wedge indicating that the cells most sensitive to intracellular displacement of amyloplasts and thus gravity sensing are confined to the top 1 mm portion of barley coleoptiles. Similar experiments with tomato hypocotyls (Lycopersicum esculentum) also resulted in curvature toward the HGMF. The data strongly support the amyloplast-based gravity-sensing system in higher plants and the usefulness of HGMF to substitute gravity in shoots.

NASA Discipline Plant Biology↗

NASA NTRS1986

A model of impulsive acceleration and transport of energetic particles in Mercury's magnetosphere

A qualitative model of substorm processes in the Mercury magnetosphere is presented based on Mariner 10 observations obtained in 1974-1975. The model is predicated on close analogies observed with the terrestrial case. Particular emphasis is given to energetic particle phenomena as observed by Mariner on March 29, 1974. The suggestion is supported that energetic particles up to about 500 keV are produced by strong induced electric fields at 3 to about 6 Mercury radii in the Hermean tail in association with substorm neutral line formation. The bursts of energetic particles produced are, in this model, subsequently confined on closed field lines near Mercury and drift adiabatically on quasi-trapped orbits for many tens of seconds. Such gradient and curvature drift of the particles can explain prominent periodicities of 5-10 s seen in the Mariner for greater than 170-keV electron flux profiles.

Baker, D. N.↗

DOE OSTI2025

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗

DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1999

CUSP Energetic Particles: Confinement, Acceleration and Implications

The cusp energetic particle (CEP) event is a new magnetospheric phenomenon. The events were detected in the dayside cusp for hours, in which the measured helium ions had energies up to 8 MeV. All of these events were associated with a dramatic decrease and large fluctuations in the local magnetic field strength. During January 1999 - December 1999 covered by this report, I have studied the CEP events by analyzing the POLAR, GEOTAIL, and WIND particle and magnetic field data measured during the geomagnetic quiet periods in 1996 and one geomagnetic storm period in 1998. The simultaneous observations indicated that the ion fluxes in the CEP events were higher than that in both the upstream and the downstream from the bow shock. The pitch angle distribution of the helium ions in the CEP events was found to peak around 90 deg. It was found that the mirror parameter, defined as the ratio of the square root of the integration of the parallel turbulent power spectral component over the ultra-low frequency (ULF) ranges to the mean field in the cusp, is correlated with the intensity of the cusp MeV helium flux, which is a measure of the influence of mirroring interactions and an indication of local effect. It was also found that the turbulent power of the local magnetic field in the ultra-low frequency (ULF) ranges is correlated with the intensity of the cusp energetic helium ions. Such ULF ranges correspond to periods of about 0.33-500 seconds that cover the gyroperiods, the bounce periods, and the drift periods of the tens keV to MeV charged particles when they are temporarily confined in the high-altitude dayside cusp. These observations represent a discovery that the high-altitude dayside cusp is a new acceleration and dynamic trapping region of the magnetosphere. The cusp geometry is connected via gradient and curvature drift of these energized ions to the equatorial plasma sheet as close as the geostationary orbit at local midnight. It implies that the dayside cusp is potentially an important source of magnetospheric particles. The discovery of the CEP events has been recognized as one of the most significant results from POLAR. I was invited to give a talk at 1999 IUGG meeting to interpret the CEP events. This discovery has also been written into the web-based Space Physics Text Book (http://www.oulu.fi/- spaceweb/textbook/cusp.html).

Chen, Jiasheng↗

NASA NTRSDate not supplied

Flow Characterization of the NASA Langley Unitary Plan Wind Tunnel, Test Section 2: Computational Results

This is an abstract for an invited paper at the AIAA Aviation Conference, June 2021. The work described here is part of an effort of coordinated experiments in the Unitary Plan Wind Tunnel (UPWT) facility at the NASA Langley Research Center (LaRC) and matching CFD simulations. The primary goal of the work is to assess the productivity and true predictive accuracy of CFD, absent any guidance from experiment, in the high supersonic speed range as compared to experiments performed at the NASA LaRC’s UPWT facility. This report concerns CFD simulation of the primary flow-path in the empty wind tunnel, including the settling chamber, nozzle, test section, and some of the tunnel downstream of the test section. The Mach number in the test section ranges from M~2.4 to M~4.6, and the required area ratio variation is achieved by translation of a nozzle block which constricts the area of a loosely S-shaped throat. Flow past protuberances in the settling chamber and into this S-bend throat are predicted by CFD to generate streamwise vorticity that may, or may not, persist through the throat and into the test section as coherent vortices. Some flow conditions are notably unsteady at frequencies well below those of turbulence, due to unsteady separated flow ahead of the nozzle block. The bulk flow moves at velocities ranging from 'walking speed' in the settling chamber to M~4.6 in the test section. Heat transfer to the settling chamber walls and buoyancy are significant at high Mach number. Subtle variations in surface curvature in the nozzle generate Mach waves that propagate into the test section. CFD of the empty tunnel serves two purposes. Firstly, the full-tunnel simulations are used to provide upstream boundary conditions for CFD of vehicle aerodynamics which are generally performed in a domain confined to the wind tunnel test section; these companion studies are addressed in other papers. Secondly, it is a challenging test for CFD to resolve all of the empty tunnel flow phenomena relevant to flow in the test section. It requires a more complete definition of geometry than was originally anticipated. In addition, it requires good spatial and temporal accuracy, and turbulence modeling that performs well on specific phenomena such as corner flows. The boundary conditions and solution algorithms must perform well from incompressible to almost hypersonic speeds. The final state of the pre-test CFD was a product of an iterative self-improvement process. The initial simulations of the empty tunnel were deficient in many respects, but hints to those deficiencies were recognized in the solutions, and remedies were implemented. Possible further improvements will be studied in the post-test phase when comparisons with experimental data are possible. The CFD was performed by five separate collaborative teams using four different flow solvers: FUN3D, Overflow, Star-CCM+ and USM3D. The level of effort of these teams varied significantly, but each made important contributions to the goals of the work. A concerted effort to use uncertainty quantification methods (UQ) in CFD is also a goal of this work. To this end, variations in CFD results due to grid refinement, turbulence modeling, and boundary conditions have been characterized. Code-to-code variation is another means of assessing CFD uncertainty. All CFD solvers predict similarity among the primary flow features; these include the variation of Mach number due to changes in Reynolds number, and the bulk flow angularity due to tunnel-wall curvature. All CFD solvers also predict similar trends in secondary flows, such as the downwash in the side-wall boundary layers. Three of the high-spatial resolution simulations give similar predictions of a complex secondary flow phenomena, streamwise vortices generated in the S-bend throat that persist into the side-wall boundary layers of the test section. Two of the highest-resolution simulations, run with the same turbulence model in different CFD solvers, gave encouragingly similar predictions of a complex tertiary flow phenomenon, small transient "sprites" of upwash flows, resulting from vortices that presumably originate in the separated flow near the leading edge of the nozzle block. While the CFD was done in a "blind pre-test" mode, requests from the experimental team for CFD results pertaining to unsteadiness and total temperature variations in the test section resulted in CFD runs that included a cooled wall in settling chamber. This then led to a change in the standard practice for running the Overflow results, which would not have occurred without this "release" of this experimental information. CFD was also used to guide some measurements. The paper will focus on establishing the consensus among CFD results and understanding differences among those results. Initial findings from the efforts to characterize CFD uncertainty have been done and will be included in the paper.

Robert Edward Childs↗

DOE OSTI2024

Ferroelectric Semimetals with α-Bi/SnSe van der Waals Heterostructures and Their Topological Currents

We show that proximity effects can be utilized to engineer van der Waals heterostructures (vdWHs) displaying semimetallic spin-ferroelectricity locking, where ferroelectricity and semimetallic spin states are confined to different layers, but are correlated by means of proximity effects. Our findings are supported by first principles calculations involving α-Bi/SnSe bilayers. We show that such systems support ferroelectrically switchable nonlinear anomalous Hall effect originating from large Berry curvature dipoles as well as direct and inverse spin Hall effects with giant bulk spin-charge interconversion efficiencies. In conclusion, the giant efficiencies are consequences of the proximity-induced semimetallic nature of low energy electron states, which are shown to behave as two-dimensional pseudo-Weyl fermions by means of symmetry analysis and first principles calculations as well as direct angle-resolved photoemission spectroscopy measurements.

2-dimensional systems↗

NASA NTRS1992

Observation of meteors by MST radar

The observation of meteor trails by a vertical mesosphere - stratosphere - troposphere (MST) radar beam has the advantage of good height resolution and an approximate knowledge of the zenith angle since the trails are horizontal or near-horizontal. An extension of the ablation theory of meteors was developed for near horizontal trails which takes into account the curvature of the earth. Observations of the Geminid meteor shower by MST radar reveal the 'diffusion heights' to be in fair agreement with the true height, but with some discrepancies that can amount to 4 km. The true heights are almost entirely confined to the range 87-91 km, although the upper limit is attributed to the coherent integration time of the existing MST radar processing.

Jones, William↗

NASA NTRS2001

Carbon Nanotubes in Water: MD Simulations of Internal and External Flow, Self Organization

We have developed computational tools, based on particle codes, for molecular dynamics (MD) simulation of carbon nanotubes (CNT) in aqueous environments. The interaction of CNTs with water is envisioned as a prototype for the design of engineering nano-devices, such as artificial sterocillia and molecular biosensors. Large scale simulations involving thousands of water molecules are possible due to our efficient parallel MD code that takes long range electrostatic interactions into account. Since CNTs can be considered as rolled up sheets of graphite, we expect the CNT-water interaction to be similar to the interaction of graphite with water. However, there are fundamental differences between considering graphite and CNTs, since the curvature of CNTs affects their chemical activity and also since capillary effects play an important role for both dynamic and static behaviour of materials inside CNTs. In recent studies Gordillo and Marti described the hydrogen bond structure as well as time dependent properties of water confined in CNTs. We are presenting results from the development of force fields describing the interaction of CNTs and water based on ab-initio quantum mechanical calculations. Furthermore, our results include both water flows external to CNTs and the behaviour of water nanodroplets inside heated CNTs. In the first case (external flows) the hydrophobic behaviour of CNTs is quantified and we analyze structural properties of water in the vicinity of CNTs with diagnostics such as hydrogen bond distribution, water dipole orientation and radial distribution functions. The presence of water leads to attractive forces between CNTs as a result of their hydrophobicity. Through extensive simulations we quantify these attractive forces in terms of the number and separation of the CNT. Results of our simulations involving arrays of CNTs indicate that these exhibit a hydrophobic behaviour that leads to self-organising structures capable of trapping water clusters. In the second case (internal flows) we study the behaviour of water droplets confined inside CNTs. Constant temperature simulations allow us to capture structural properties such as the contact angles and density profiles of the equilibrated drops. By heating and subsequently cooling of the CNT, we are able to measure the evaporation and the condensation rate of the entrapped water.

Jaffe, Richard L.↗

DOE OSTI2025

A curvature-based mechanism for the spontaneous emission of electromagnetic radiation at curved plasma–vacuum interfaces

In this work, we derive from first principles the dispersion relation for electromagnetic radiation that is spontaneously emitted at a 2D curved plasma–vacuum interface in the absence of both background and external electric and magnetic fields. This plasma geometry is motivated by plasmas generated through hypervelocity impact projectiles—an area that is of importance to satellite safety in the aerospace community, and may be of importance for the developing study of projectile impact-driven inertial confinement fusion. We find spontaneous radiation occurs when thermal surface waves propagate along the edge of a circular plasma, with a frequency lower than the plasma frequency and spatial decay proportional to 1/r. This is distinct from more well-known mechanisms of radiation, in which the radiation frequency is greater than or equal to the plasma frequency. This theory is validated against 2D, electromagnetic particle-in-cell simulations and shows reasonable agreement. Finally, we discuss the effects of plasma geometry and attribute the radiation generation mechanism to the acceleration of charged surface waves as they traverse azimuthally along a curved path.

Lau, Raymond (ORCID:0000000299318101)↗

DOE OSTI2025

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗

DOE OSTI2026

Detecting shearless phase-space transport barriers in global gyrokinetic turbulence simulations with test particle map models

In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.

fusion plasma↗