Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “oscillating flow”

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

Effect of the divertor configuration on the JET edge radial electric field

The role of the divertor configuration and divertor plasma physics on the L–H transition is poorly understood, leading to large uncertainties in predicting the L–H power threshold in future devices. This contribution reports on edge perpendicular plasma flow measurements by Doppler backscattering in JET L–H transition experiments with the outer divertor strike-point at different positions: horizontal target (HT), vertical target (VT), and in the corner configuration (between the HT and VTs). The edge perpendicular flow was found to be significantly affected by changes in the divertor configuration in the region inside the separatrix. Our results do not show evidence for the existence of a critical edge flow shear needed to achieve H-mode for different divertor configurations, with a larger shear observed for the VT configuration. No significant change in the shear flow and in the density fluctuation level is measured preceding the L–H transition in the region just inside the separatrix. The dynamics of the L–H transitions are also influenced by the divertor configuration with divertor oscillations (DOs) observed only in the HT configuration. Interestingly, DOs are associated with marked changes in the edge perpendicular flow around the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Convective Momentum Transport and Its Impact on the Madden‐Julian Oscillation in E3SM‐MMF

Abstract Convective momentum transport (CMT) is the process of vertical redistribution of horizontal momentum by small‐scale turbulent flows from moist convection. Traditional general circulation models (GCMs) and their multiscale modeling framework (MMF) versions poorly represent CMT due to insufficient information of subgrid‐scale flows at each GCM grid. Here the explicit scalar momentum transport (ESMT) scheme for representing CMT is implemented in the Energy Exascale Earth System Model‐Multiscale Modeling Framework (E3SM‐MMF) with embedded 2‐D cloud‐resolving models (CRMs), and verified against E3SM‐MMF simulations with 3‐D CRMs and observations. The results show that representing CMT by ESMT helps reduce climatological mean precipitation model bias over the western Pacific and the ITCZ regions, which is attributed to the weakened mean easterlies over the Pacific. Also, CMT from simulations with 2‐D and 3‐D CRMs impose a similar impact on Kelvin waves by reducing their variability and slowing down their phase speed, but opposite impacts on the Madden‐Julian Oscillation (MJO) variability. The ESMT scheme readily captures the climatological mean spatial patterns of the zonal and meridional components of CMT and their variability across multiple time scales, but shows some differences in estimating its magnitude. CMT mainly affects the MJO by decelerating its winds in the free troposphere, but accelerates its near‐surface winds. This study serves as a prototype for implementing CMT scheme in the MMF simulations, highlighting its crucial role in reducing model bias in mean state and spatiotemporal variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Theoretical modeling of a bottom-raised oscillating surge wave energy converter structural loadings and power performances

Here, this study presents theoretical formulations to evaluate the fundamental parameters and performance characteristics of a bottom-raised oscillating surge wave energy converter (OSWEC) device. Employing a flat plate assumption and potential flow formulation in elliptical coordinates, closed-form equations for the added mass, radiation damping, and excitation forces/torques in the relevant pitch-pitch and surge-pitch directions of motion are developed and used to calculate the system's response amplitude operator and the forces and moments acting on the foundation. The model is benchmarked against numerical simulations using WAMIT and WEC-Sim, showcasing excellent agreement. The sensitivity of plate thickness on the analytical hydrodynamic solutions is investigated over several thickness-to-width ratios ranging from 1:80 to 1:10. The results show that as the thickness of the benchmark OSWEC increases, the deviation of the analytical hydrodynamic coefficients from the numerical solutions grows from 3% to 25%. Differences in the excitation forces and torques, however, are contained within 12%. While the flat plate assumption is a limitation of the proposed analytical model, the error is within a reasonable margin for use in the design space exploration phase before a higher-fidelity (and thus more computationally expensive) model is employed. A parametric study demonstrates the ability of the analytical model to quickly sweep over a domain of OSWEC dimensions, illustrating the analytical model's utility in the early phases of design.

13 HYDRO ENERGY↗
DOE OSTI2024

Mechanisms of Low-Level Jet Formation in the U.S. Mid-Atlantic Offshore

Abstract Low-level jets (LLJs), in which the wind speed attains a local maximum at low altitudes, have been found to occur in the U.S. mid-Atlantic offshore, a region of active wind energy deployment as of 2023. In contrast to widely studied regions such as the U.S. southern Great Plains and the California coastline, the mechanisms that underlie LLJs in the U.S. mid-Atlantic are poorly understood. This work analyzes floating lidar data from buoys deployed in the New York Bight to understand the characteristics and causes of coastal LLJs in the region. Application of the Hilbert–Huang transform, a frequency analysis technique, to LLJ case studies reveals that mid-Atlantic jets frequently occur during times of adjustment in synoptic-scale motions, such as large-scale temperature and pressure gradients or frontal passages, and that they do not coincide with motions at the native inertial oscillation frequency. Subsequent analysis with theoretical models of inertial oscillation and thermal winds further reveals that these jets can form in the stationary geostrophic wind profile from horizontal temperature gradients alone—in contrast to canonical LLJs, which arise from low-level inertial motions. Here, inertial oscillation can further modulate the intensity and altitude of the wind speed maximum. Statistical evidence indicates that these oscillations arise from stable stratification and the associated frictional decoupling due to warmer air flowing over a cold sea surface during the springtime land–sea breeze. These results improve our conceptual understanding of mid-Atlantic jets and may be used to better predict low-level wind speed maxima. Significance Statement The purpose of this work is to identify and characterize the atmospheric mechanisms that result in an occasional low-level maximum in the wind speed off the U.S. mid-Atlantic coastline. Our findings show that these low-level jets form due to horizontal temperature gradients arising from fronts and synoptic systems, as well as from the land–sea breeze that forces warmer air over the cold ocean surface. This work aids predictability of such jets, improves our understanding of this coastal environment, and has implications for future deployment of offshore wind energy in this region.

17 WIND ENERGY↗
DOE OSTI2023

Turbulence and E × B flow correlations across the L-H transition in DIII-D deuterium and hydrogen plasmas

The isotope dependence of the low- to high confinement-mode (L- to H-mode) transition power threshold P LH presents significant challenges for the initial (non-nuclear) hydrogen operations phase of the international thermonuclear experimental reactor (ITER). Here, we examine the isotope dependence of turbulence and E × B flow correlation properties in the L-mode edge plasma, leading up to the L-H transition. At marginal auxiliary power (near P LH ), turbulence is initially suppressed periodically during limit cycle oscillations (LCO) that precede the transition to sustained H-mode confinement. Here, we present evidence that the long-range (toroidal) correlation of the E × B edge plasma flow across the LCO phase is much weaker in hydrogen than in deuterium in the DIII-D tokamak, congruent with the higher threshold power P LH required to access LCO and H-mode in hydrogen. Concomitantly, the time required to initially quench edge turbulence via localized edge E × B flow shear is significantly longer in hydrogen (1–1.5 ms) than in deuterium (approx. 100 µs). No toroidal long-range correlation of the turbulence amplitude is observed, in agreement with expectations based on the relatively short poloidal turbulence correlation length. Radial edge turbulence and flow correlation lengths are longer in deuterium than in hydrogen plasmas as one would expect from ‘naïve’ gyro-Bohm isotope transport scaling, despite the substantially higher thermal flux across the last closed flux surface in hydrogen before the L-H transition. This article is part of a discussion meeting issue 'H-mode transition and pedestal studies in fusion plasmas'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Code Demonstration: Approximation of nearly-periodic symplectic maps via structure-preserving neural networks

In this notebook, we use the symplectic gyroceptron architecture from [Duruisseaux et al., 2022] to learn a surrogate map for the nearly-periodic symplectic flow map associated to a nearly-periodic Hamiltonian system composed of two nonlinearly coupled oscillators, where one of them oscillates significantly faster than the other: $\Big\lbrace$ $ {\dot{q}_1 = p_1 \hspace{5mm} \dot{p}_1 = -q_1 - \epsilon \partial_{q1}U(q_1 , q_2 )} \atop {\dot{q}_2 = \epsilon p_2 \hspace{4mm} \dot{p}_2 = -\epsilon q_2 - \epsilon \partial_{q2}U(q_1 , q_2 )}$. These equations of motion are the Hamilton's equations associated to the Hamiltonian $H_{\epsilon}(q_1,q_2,p_1,p_2)=\frac{1}{2}(q^2_1+p^2_1)+\frac{1}{2}\epsilon(q^2_2+p^2_2)+\epsilon$$U(q_1,q_2)$.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Hydrokinetic Power Conversion Using Vortex-Induced Oscillation with Cubic Restoring Force

A cubic-spring restoring function with high-deformation stiffening is introduced to passively improve the harnessed marine hydrokinetic power by using flow-induced oscillations/vibrations (FIO/V) of a cylinder. In these FIO/V experiments, a smooth, rigid, single-cylinder on elastic end-supports is tested at Reynolds numbers ranging from 24,000 < Re < 120,000. The parameters of the tested current energy converter (CEC) are cubic stiffness and linear damping. Using the second generation of digital virtual spring-damping (Vck) controller developed by the Marine Renewable Energy Laboratory (MRELab), the cubic modeling of the oscillator stiffness is tested. Experimental results show the influence of the parameter variation on the amplitude, frequency, energy conversion, energy efficiency, and power of the converter. All experiments are conducted in the low turbulence-free surface water (LTFSW) channel of the MRELab of the University of Michigan. The main conclusions are: (1) The nonlinearity in the cubic oscillator is an effective way to extend the vortex-induced vibration (VIV) upper branch, which results in higher harnessing power and efficiency compared to the linear stiffness cylinder converter. (2) Compared to the linear converter, the overall power increase is substantial. The nonlinear power optimum, occurring at the end of the VIV upper branch, is 63% higher than its linear counterpart. (3) The cubic stiffness converter with low harnessing damping achieves consistently good performance in all the VIV regions because of the hardening restoring force, especially at higher flow velocity.

13 HYDRO ENERGY↗
DOE OSTI2021

Cooling and instabilities in colliding flows

ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Flowing plasma rearrangement in the presence of static perturbing fields

Charged particles interacting with electromagnetic waves have a portion of their energy tied up in wave-driven oscillations. When these waves are localized to the exhaust of linear magnetic confinement systems, this ponderomotive effect can be utilized to enhance particle confinement. The same effect can be derived for particles moving via an E×B drift into a region of a static perturbation to the electromagnetic fields which has a large wave vector component in the direction of the motion. In this work, we use a simplified slab model to self-consistently solve for the electromagnetic fields within the fluid flowing plasma of a static flute-like (k∥=0) perturbation and evaluate the resulting ponderomotive potential. We find that two types of perturbations can exist within the flowing plasma, which are an O wave and an X wave in the frame moving with the fluid. In the case of tenuous plasma, these perturbations are magnetostatic or electrostatic multipole-analog perpendicular to the guiding magnetic field in the lab frame, respectfully. For denser plasmas, the O wave-like perturbation is screened at the electron skin depth scale, and the X wave-like perturbation is a combination of a similar perpendicular electric perturbation and parallel magnetic perturbation. The ponderomotive potential generated in the X wave-like case is gyrofrequency-dependent and can be used as either potential barriers or potential wells, depending on the direction of the flow velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Computational Characterization and Boundary Condition Models of a Jet Interaction Fluidic Oscillator

Fluidically oscillating jet actuators that generate sweeping jets when supplied with a pressurized fluid have been used to mitigate separation and reduce drag in a range of flow control applications. Their implementation in future flight vehicles would require fundamental understanding and accurate predictive techniques of the physics of their internal flow and jet formation. The present investigations focus on high-fidelity, time-accurate simulations to characterize the flow physics of the actuator in quiescent conditions. An important element of the present simulations is to demonstrate the ability of the computational fluid dynamics (CFD) solver to predict the jet characteristics and provide a basis for the development of improved boundary conditions (BC) without entirely resolving the geometrical features of the fluidic device. The CFD-predicted oscillation frequencies of the engendered jets were found to be in excellent agreement with experiments, even on two-dimensional meshes. The study revealed that three-dimensional simulations are required to capture some of the flow features of the sweeping jet such as the double peak in time-averaged velocity distributions downstream of the actuator’s orifice that were measured in experiments. Furthermore, several approaches for modeling the actuator were implemented and assessed in quiescent conditions. The evaluation of a boundary condition at the device throat, based on the phase-averaged flow variables, provides the basis for devising surface-based boundary conditions. The influence and necessity of including turbulent characteristics as part of the boundary conditions have also been identified.

42 ENGINEERING↗
DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

INL POSTER 2025- KRYSTAL FOSTER

The Advanced Test Reactor (ATR) employs the N-16 system to measure reactor power by detecting beta particles emitted from nitrogen-16 isotopes. During startup, the N-16 pumps cause rapid oscillations in the discharge piping due to pressure transients, which cease when the valve supplying N-16 flow to the reactor is opened. To address this, a soft-starter is proposed to gradually increase the motor voltage, reducing mechanical stress and preventing pressure surges. This solution enhances system stability, lowers maintenance costs, extends equipment lifespan, and is more cost-effective than Variable Frequency Drives (VFDs), making it ideal for the infrequently activated N-16 pumps.

42 - ENGINEERING↗
DOE OSTI2023

A numerical extension of the spatially-filtered Euler equations for contact discontinuities

Solving the spatially-filtered Euler equations with a kinetic energy-preserving centered discretization is an efficient framework for simulating shock-dominated flows. Here, the present work describes a numerical treatment of the spatially-filtered Euler equations to minimize oscillations when simulating contact discontinuities. To counteract dispersive errors inherent to centered schemes, a WENO-like correction term is applied to the enthalpy transport in the energy equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A new form of the Saint-Venant equations for variable topography

The solution stability of river models using the one-dimensional (1D) Saint-Venant equations can be easily undermined when source terms in the discrete equations do not satisfy the Lipschitz smoothness condition for partial differential equations. Although instability issues have been previously noted, they are typically treated as model implementation issues rather than as underlying problems associated with the form of the governing equations. This study proposes a new reference slope form of the Saint-Venant equations to ensure smooth slope source terms and eliminate one source of potential numerical oscillations. It is shown that a simple algebraic transformation of channel geometry provides a smooth reference slope while preserving the correct cross-section flow area and the total Piezometric pressure gradient that drives the flow. The reference slope method ensures the slope source term in the governing equations is Lipschitz continuous while maintaining all the underlying complexity of the real-world geometry. The validity of the mathematical concept is demonstrated with the open-source Simulation Program for River Networks (SPRNT) model in a series of artificial test cases and a simulation of a small urban creek. Validation comparisons are made with analytical solutions and the Hydrologic Engineering Center's River Analysis System (HEC-RAS) model. The new method reduces numerical oscillations and instabilities without requiring ad hoc smoothing algorithms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Geometric control of hyperbolic exciton-polariton condensate dimers

Coupled many-body quantum systems exhibit rich emergent physics with diverse stationary and dynamical behaviours. By engineering platforms with tunable and distinct coupling mechanisms, new insights emerge into the collective behaviour of coupled many body systems. Particles can be exchanged via evanescent or ballistic coupling: the former, based on proximity, yields large spectral splitting, while the latter requires strict phase-matching, analogous to phase-coupled harmonic oscillators and has a smaller impact on the energy landscape. We demonstrate an all-optically tunable quantum fluid dimer based on exciton-polariton condensates in a photonic crystal waveguide with hyperbolic (saddle-like) dispersion. Varying the dimer’s angle relative to the grating tunes the coupling from evanescent to ballistic. We directly observe spectral features and mass flow shaped by the saddle dispersion. This work highlights photonic crystals as powerful platforms to explore condensed matter phenomena lying at the interface between delay-coupled nonlinear oscillators and tight binding physics.

Georgakilas, Ioannis [IBM Research-Zurich, Rüschli↗
DOE OSTI2021

Chemical pumps and flexible sheets spontaneously form self-regulating oscillators in solution

The synchronization of self-oscillating systems is vital to various biological functions, from the coordinated contraction of heart muscle to the self-organization of slime molds. Through modeling, we design bioinspired materials systems that spontaneously form shape-changing self-oscillators, which communicate to synchronize both their temporal and spatial behavior. Here, catalytic reactions at the bottom of a fluid-filled chamber and on mobile, flexible sheets generate the energy to “pump” the surrounding fluid, which also transports the immersed sheets. The sheets exert a force on the fluid that modifies the flow, which in turn affects the shape and movement of the flexible sheets. This feedback enables a single coated (active) and even an uncoated (passive) sheet to undergo self-oscillation, displaying different oscillatory modes with increases in the catalytic reaction rate. Two sheets (active or passive) introduce excluded volume, steric interactions. This distinctive combination of the hydrodynamic, fluid–structure, and steric interactions causes the sheets to form coupled oscillators, whose motion is synchronized in time and space. We develop a heuristic model that rationalizes this behavior. These coupled self-oscillators exhibit rich and tunable phase dynamics, which depends on the sheets’ initial placement, coverage by catalyst and relative size. Moreover, through variations in the reactant concentration, the system can switch between the different oscillatory modes. This breadth of dynamic behavior expands the functionality of the coupled oscillators, enabling soft robots to display a variety of self-sustained, self-regulating moves.

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

Guiding center and gyrokinetic orbit theory for large electric field gradients and strong shear flows

The guiding center and gyrokinetic theory of magnetized particle motion is extended to the regime of large electric field gradients perpendicular to the magnetic field. A gradient in the electric field directly modifies the oscillation frequency and causes the Larmor orbits to deform from circular to elliptical trajectories. In order to retain a good adiabatic invariant, there can only be strong dependence on a single coordinate at lowest order, so that resonances do not generate chaotic motion that destroys the invariant. When the gradient across magnetic flux surfaces is dominant, the guiding center drift velocity becomes anisotropic in response to external forces and additional curvature drifts must be included. The electric polarization density remains gyrotropic, but both the polarization and magnetization are modified by the change in gyrofrequency. The theory can be applied to shear flows that are even stronger than those observed in the edge transport barrier of a high-performance tokamak (H-mode) pedestal, even if the toroidal field is as small as or even smaller than the poloidal field. Yet, the theory retains a mathematical form that is similar to the standard case and can readily be implemented within existing simulation tools.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗