SEARCH · Engineering Papers
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.
The barotropic instability of zonal flows due to infinitesimal normal-mode perturbations is considered, with the zonal flow assumed to be continuous and either monotonic or nonmonotonic and allowed to have one or more inflection points. A sufficient condition for instability is derived for this general flow profile that complements that of Arnol'd (1965), which restricted itself to the prediction of stability only.
Limit-cycle oscillations are studied for ion temperature gradient turbulence, which, in the absence of large diamagnetic (mean) shear flows, saturates through energy transfer from unstable modes to large-scale stable modes via zonal-flow intermediary modes. Oscillations of zonal flow and turbulence levels are strongly constrained by the reactive, largely non-dissipative character of the zonal flows. Since existing predator–prey models for observed oscillations in experiments do not include energy transfer through zonal flows to stable modes, low-order fluid models with this physics are constructed and investigated. A simple three-wave truncation produces low-amplitude zonal flows that slowly oscillate around a zero mean, with turbulence oscillations between coupled wavenumbers that exceed linear frequencies by orders of magnitude. This inconsistency with experimental observations is caused by the weak non-linear drive of zonal flows in three-wave systems and the lack of multiple-wavenumber turbulent interactions. A more comprehensive model that preserves multiple wavenumber interactions within the context of conservative zonal-flow-mediated energy transfer to stable modes accurately reflects observed dynamics when the phase between stable and unstable modes is occasionally randomized.
Gyrokinetic and kinetic-MHD simulations are performed for the fishbone instability in the DIII-D discharge #178631, chosen for validation of first-principles simulations to predict the energetic particle (EP) transport in an ITER prefusion baseline scenario. Fishbone modes are found to generate zonal flows, which dominate the fishbone saturation. The underlying mechanisms of the two-way fishbone-zonal flows nonlinear interplay are discussed in details. Numerical and analytical analyses identify the fishbone-induced EP redistribution as the dominant generation mechanism for zonal flows. The zonal flows modify the nonlinear dynamics of phase space zonal structures, which reduces the amount of EPs able to resonate with the mode, leading to an early fishbone saturation. Simulation results including zonal flows agree quantitatively with DIII-D experimental measurements of the fishbone saturation amplitude and EP transport, supporting this novel saturation mechanism by self- generated zonal flows. Moreover, the wave-particle mode-locking mechanism is shown to determine quantitatively the fishbone frequency down-chirping, as evident in GTC simulation results in agreement with predictions from analytical theory. Finally, the
fishbone-induced zonal flows are possibly responsible for the formation of an ion-ITB in the DIII-D discharge. Based on the low EP transport and the large zonal flow shearing rates associated with the fishbone instability in gyrokinetic simulations of the ITER scenario, it is conjectured that high performance scenarios could be designed in ITER burning plasmas through fishbone-induced ITBs.
The role of zonal flow (ZF) in the turbulence saturation of trapped electron modes (TEM) in magnetically confined plasmas is revisited. Here, in this study, we examine ZF excitation and saturation mechanisms in TEM turbulence using detailed free-energy transfer diagnostics from nonlinear gyrokinetic simulations. When the ion channel becomes subdominant, electrons shift to transferring energy to the zonal component, which results in zonal flows always playing a significant though subdominant role, in the case of the temperature gradient driven TEM. The velocity-space structure of energy transfer reveals that trapped electrons directly couple to ZF. Regarding the saturation physics, the zonal flow advection, drift-wave–drift-wave interactions, and stable modes at the pump waves are compared. When ZF saturation is weak, drift-wave–drift-wave interactions play a larger role in maintaining turbulence saturation. The findings elucidate the roles of nonzonal to zonal flow coupling, stable modes, and nonzonal wave-wave interactions, and provide a basis for improving reduced models of turbulent transport.
This study addresses key aspects of momentum transport in hydrodynamic disks, which is critical for understanding zonal flow generation and turbulence in compressible hydrodynamic disks. We find that nonlinear momentum/density transport leads to the formation of zonal flows from the Rossby wave instability in disks. We analytically derive the generation and location of zonal flows and describe a modified Taylor identity applicable to compressible disk flows. We further present a self-regulation model, revealing a dynamic interplay between zonal flow and fluctuations driven by Rossby wave instability that regulates the nonlinear saturation state. This theoretical framework contributes insights into the dynamics of disks such as protoplanetary disks, shedding light on the intricate processes governing momentum/density transport and the emergence of zonal flows in the saturation of protoplanetary disks.
Strongly driven ion-scale turbulence in tokamak plasmas is shown to be regulated by a new propagating zonal flow mode, the toroidal secondary mode, which is nonlinearly supported by the turbulence. The mode grows and propagates due to the combined effects of zonal flow shearing and advection by the magnetic drift. Above a threshold in the turbulence level, small-scale toroidal secondary modes become unstable and shear apart turbulent eddies, forcing the turbulence level to remain near the threshold. This threshold condition is used to derive scaling laws for the turbulent heat flux, fluctuation spectra, and zonal flow amplitude, which are validated in nonlinear gyrokinetic simulations and explain previous experimental observations.
Dataset and input files to generate figures in "Theory of zonal flow generation and propagation in toroidal geometry."
The linear theory of deep zonal flows developed by Busse (1976) for the origins of deep motions of the atmospheres of Jupiter and Saturn is extended into the nonlinear regime. Relationships for the relative magnitudes of convective heat and momentum transports are formulated. A perturbation approach is taken to the problem, with the amplitude of the convection serving as the small parameter, and the basic equations being expanded in terms of the Prandtl number. The Boussinesq approximation is employed, together with an assumption of a low Rossby number for the Jovian and Saturn atmospheres. Differences in the amplitude of the Jovian equatorial jet relative to that of Saturn are explored in terms of a low equatorial convective heat flux on Jupiter.
Integrations of the linear stability problem, in the present numerical investigation of weakly nonparallel zonal flows barotropic instability having localized intense shear regions, reveal the existence of unstable localized wave packets. The spatial structure and eigenfrequencies of these packets depend on two parameters measuring the degree of supercriticality and the zonal length scale of the shear region. It is found that instability structure is defined by conditions ensuring the decay of the wave packet at infinity, and the transition from long to short waves across a turning point region that is controlled by nonparallel effects.
Gierasch et al. (1979) have shown that the question of the stability of zonal flows on Jupiter is influenced strongly by the lower boundary condition. In the case of Jupiter, the troposphere rests above a deep layer of fluid which is probably extremely close to adiabatic. The present study extends the earlier investigations to a wider range of conditions, taking into account the nongeostrophic limit and the presence of horizontal shear. The major difference between the behavior of the model derived in the study and the behavior of a model with a rigid lower boundary, corresponding to the terrestrial conditions, is that baroclinic instability by the Eady mechanism is suppressed. The influence of beta (Charney-type instability) operates only at the top of the troposphere, not at its base. The suppression of the Eady modes seems to apply all the way down to the nongeostrophic limit where the Rossby number is order unity.
The toroidal geometry of tokamaks and stellarators is known to play a crucial role in the linear physics of zonal flows (ZFs), leading to e.g. the Rosenbluth–Hinton residual and geodesic acoustic modes. However, descriptions of the nonlinear ZF growth from a turbulent background typically resort to simplified models of the geometry. We present a generalised theory of the secondary instability to model the ZF growth from turbulent fluctuations in toroidal geometry, demonstrating that the radial magnetic drift substantially affects the nonlinear ZF dynamics. In particular, the toroidicity gives rise to a new branch of propagating ZFs, the toroidal secondary mode, which is nonlinearly supported by the turbulence. We present a theory of this mode and compare the theory against gyrokinetic simulations of the secondary mode. The connection with other secondary modes—the ion-temperature-gradient and Rogers–Dorland–Kotschenreuther secondary modes—is also examined.
In order to reexamine the magnitude differences of the Jovian atmosphere's jets, as determined by Voyager 1 and 2 images, a novel approach is used to ascertain the zonal mean east-west component of motion. This technique is based on digital pattern matching, and is applied on pairs of mapped images to yield a profile of the mean zonal component that reproduces the exact locations of the easterly and westerly jets between + and - 60 deg latitude. Results have been obtained for all of the Voyager 1 and 2 cylindrical mosaics; the correlation coefficient between Voyagers 1 and 2 in mean zonal flow between + and - 60 deg latitude, determined from violet filter mosaics, is 0.998.
In order to reexamine the magnitude differences of the Jovian atmosphere's jets, as determined by Voyager 1 and 2 images, a novel approach is used to ascertain the zonal mean east-west component of motion. This technique is based on digital pattern matching, and is applied on pairs of mapped images to yield a profile of the mean zonal component that reproduces the exact locations of the easterly and westerly jets between + and 60 deg latitude. Results were obtained for all of the Voyager 1 and 2 cylindrical mosaics; the correlation coefficient between Voyagers 1 and 2 in mean zonal flow between + and - 60 deg latitude, determined from violet filter mosaics, is 0.998.
Voyager imaging, infrared, and radio observations for Saturn have been recently interpreted by Smith et al. (1982) as an indication that the jet streams observed at the cloud tops extend to depths greater than the 10,000-bar level. This analysis assumes a maximum latitudinal temperature contrast of a few percent, a mean atmospheric rotation rate at depth given by Saturn's ratio period, and no variation with latitude of the bottom pressure level for the zonal flow system. These assumptions are not, however, firmly constrained by observation. The diagnostic analysis of plausible alternative configurations for Saturn's atmospheric structure demonstrates that a thin weather layer system (confined at mid to high latitudes to levels above 200 bar) cannot be excluded by any of the available observations. A quantitative estimate of the effects of moisture condensation (including the differentiation of mean molecular weight) suggests that these might provide the buoyancy contrasts necessary to support a thin-layer flow provided that Saturn's outer envelope is enriched approximately 10 times in water abundance relative to a solar composition atmosphere and strongly differentiated with latitude at the condensation level.
A simple model of the equatorial zonal jet in the giant planets is studied in which the flow is assumed uniform on cylinders concentric with the spin axis, and viscous and magnetic torques on each cylinder are balanced. This 'Taylor constraint' is solved simultaneously with the dynamo equation to obtain the velocity and magnetic field in the equatorial plane. The model is used to reproduce the widely differing jet widths of Jupiter and Saturn eddy viscosity of 2500 sq cm/s. For Saturn, the model has a large magnetic Reynolds number where the Chandrasekhar number Q = 1, and hence exhibits substantial axisymmetrization of the field in the equatorial plane.
The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.
The response of linear planetary wave 2 to changes in the isentropic zonally symmetric distribution of potential vorticity (PV) is investigated numerically. Wave 2 is sensitive to the width and position of a region where the meridional derivative of PV is weak, denoted the 'surf zone', in the middle stratosphere. A narrow surf zone leads to an amplification of wave 2, and confines some of its Eliassen-Palm wave activity in high latitudes. For wider surf zones the wave activity is concentrated near the associated critical surface, and the amplitude of the wave decreases. Large changes in the wave amplitude and in the distribution of its activity are associated with the subtle changes in the zonal winds produced by modest modifications in the distribution of PV. Basic states that include regions of reversed meridional gradients of PV lead to wave overreflection and strong poleward focusing of wave activity. The amplitude of wave 2 is enhanced in the presence of negative gradients, with large responses occurring for eastward traveling waves.
A closed-loop, overlapped, velocity-coupling procedure has been utilized to combine a two-dimensional potential flow panel code and a Navier-Stokes code. The fully coupled, two-zone code has been used to compute the flow past a NACA 0012 airfoil. For this case, the zonal method has shown that the grid domain size can be reduced to 0.14 chord lengths with less than 1 percent loss in accuracy. Further, the required computation time is reduced by a factor of approximately four.