Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA FLUID EQUATIONS”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Kelvin-Helmholtz instability is the result of parity-time symmetry breaking

Parity-Time (PT)-symmetry is being actively investigated as a fundamental property of observables in quantum physics. We show that the governing equations of the classical two-fluid interaction and the incompressible fluid system are PT-symmetric, and the well-known Kelvin-Helmholtz instability is the result of spontaneous PT-symmetry breaking. It is expected that all classical conservative systems governed by Newton's law admit PT-symmetry, and the spontaneous breaking thereof is a generic mechanism for classical instabilities. Here discovering the PT-symmetry of systems in fluid dynamics and plasma physics and identifying the PT-symmetry breaking responsible for instabilities enable new techniques to classical physics and enrich the physics of PT-symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

The Role of Extracellular Fluid in Biokinetic Modeling

Here, the pharmacokinetic equations of Pierson et al. describing the behavior of bromide in rat provide a general approach to the modeling of extracellular fluid (ECF). The movement of material into ECF spaces is rapid and is completely characterized by tissue volumes and vascular flow rates to and from a tissue, the volumes of the tissue, and the ECF associated with the tissue. Early-time measurements are needed to characterize ECF. Measurements of DTPA disappearance from plasma by Wedeking et al. are discussed as an example of such measurements. In any biokinetic model, the fastest transfer rates are not determinable with the usual datasets, and if determined empirically, these rates will have very large and highly correlated uncertainties, so particular values of these rates, even though the model fits the available data, are not significant. A pharmacokinetic front-end provides values for these fast rates. An example of such a front-end for a 200–g rat is given.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1977

Implications of Navier-Stokes turbulence theory for plasma turbulence

The methodology of Navier-Stokes fluid turbulence theory is reviewed, with emphasis placed on the relevance of the Navier-Stokes concepts for understanding plasma turbulence. After a brief consideration of the three-dimensional case, two-dimensional problems are discussed. In addition, MHD turbulence and turbulence in Vlasov plasmas are treated. In particular, the direct interaction approximation developed by Kraichnan (1959) is generalized from Navier-Stokes turbulence theory to produce a computable set of differentio-integral equations for a Vlasov plasma.

Montgomery, D.↗
DOE OSTI2016

Experimental investigation of flow induced dust acoustic shock waves in a complex plasma

We report on experimental observations of flow induced large amplitude dust-acoustic shock waves in a complex plasma. The experiments have been carried out in a Π shaped direct current glow discharge experimental device using kaolin particles as the dust component in a background of Argon plasma. A strong supersonic flow of the dust fluid is induced by adjusting the pumping speed and neutral gas flow into the device. An isolated copper wire mounted on the cathode acts as a potential barrier to the flow of dust particles. A sudden change in the gas flow rate is used to trigger the onset of high velocity dust acoustic shocks whose dynamics are captured by fast video pictures of the evolving structures. The physical characteristics of these shocks are delineated through a parametric scan of their dynamical properties over a range of flow speeds and potential hill heights. The observed evolution of the shock waves and their propagation characteristics are found to compare well with model numerical results based on a modified Korteweg-de-Vries-Burgers type equation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Symmetries of a reduced fluid-gyrokinetic system

Symmetries of a fluid-gyrokinetic model are investigated using Lie group techniques. Specifically the nonlinear system constructed by Zocco and Schekochihin (Zocco & Schekochihin 2011), which combines nonlinear fluid equations with a drift-kinetic description of parallel electron dynamics, is studied. Significantly, this model is fully gyrokinetic, allowing for arbitrary $k_{\bot \rho I}$ where $k_\bot$ is the perpendicular wave vector of the fluctuations and $\rho I}$ is the ion gyroradius. The model includes integral operators corresponding to gyroaveraging as well as the moment equations relating fluid variables to the kinetic distribution function. A large variety of exact symmetries is uncovered, some of which have unexpected form. Using these results, new nonlinear solutions are constructed, including a helical generalization of the Chapman-Kendall solution for a collapsing current sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1997

Predictive Modeling in Plasma Reactor and Process Design

Research continues toward the improvement and increased understanding of high-density plasma tools. Such reactor systems are lauded for their independent control of ion flux and energy enabling high etch rates with low ion damage and for their improved ion velocity anisotropy resulting from thin collisionless sheaths and low neutral pressures. Still, with the transition to 300 mm processing, achieving etch uniformity and high etch rates concurrently may be a formidable task for such large diameter wafers for which computational modeling can play an important role in successful reactor and process design. The inductively coupled plasma (ICP) reactor is the focus of the present investigation. The present work attempts to understand the fundamental physical phenomena of such systems through computational modeling. Simulations will be presented using both computational fluid dynamics (CFD) techniques and the direct simulation Monte Carlo (DSMC) method for argon and chlorine discharges. ICP reactors generally operate at pressures on the order of 1 to 10 mTorr. At such low pressures, rarefaction can be significant to the degree that the constitutive relations used in typical CFD techniques become invalid and a particle simulation must be employed. This work will assess the extent to which CFD can be applied and evaluate the degree to which accuracy is lost in prediction of the phenomenon of interest; i.e., etch rate. If the CFD approach is found reasonably accurate and bench-marked with DSMC and experimental results, it has the potential to serve as a design tool due to the rapid time relative to DSMC. The continuum CFD simulation solves the governing equations for plasma flow using a finite difference technique with an implicit Gauss-Seidel Line Relaxation method for time marching toward a converged solution. The equation set consists of mass conservation for each species, separate energy equations for the electrons and heavy species, and momentum equations for the gas. The sheath is modeled by imposing the Bohm velocity to the ions near the walls. The DSMC method simulates each constituent of the gas as a separate species which would be analogous in CFD to employing separate species mass, momentum, and energy equations. All particles including electrons are moved and allowed to collide with one another with the stipulation that the electrons remain tied to the ions consistent with the concept of ambipolar diffusion. The velocities of the electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. These benefits come at a price in terms of computational time and memory. The DSMC and CFD are made as consistent as possible by using similar chemistry and power deposition models. Although the comparison of CFD and DSMC is interesting, the main goal of this work is the increased understanding of high-density plasma flowfields that can then direct improvements in both techniques. This work is unique in the level of the physical models employed in both the DSMC and CFD for high-density plasma reactor applications. For example, the electrons are simulated in the present DSMC work which has not been done before for low temperature plasma processing problems. In the CFD approach, for the first time, the charged particle transport (discharge physics) has been self-consistently coupled to the gas flow and heat transfer.

Hash, D. B.↗
DOE OSTI2019

Heating and ablation of high- Z cryogenic pellets in high temperature plasmas

A kinetic model is developed here to determine the power deposition from energetic electrons into the neutral gas shield of an ablating high-Z pellet. For high-Z, the velocity distribution of the hot electrons is nearly isotropic, and we use this feature to develop solutions to the kinetic equation. In contrast to pre-existing models, we consider the effect of gyro-motion as well as elastic scattering of the hot electrons. Two limits are considered: when the gyro-frequency is much greater than the elastic collision frequency, the hot electrons diffuse longitudinally along the field lines as they slow down; but if the gyro-frequency is much less than the elastic collision frequency, the hot electrons diffuse radially. In both limits, it is possible to express the kinetic equation describing the hot electrons independently of the gas density profile, and therefore, the power deposition model is universal in this respect. The emphasis on elastic scattering yields an ablation rate which scales as Z -7/6 which is different than Z -2/3 shown in Sergeev et al (2006 Plasma Phys. Rep. 32 5). It is also demonstrated that the sheath potential required to maintain ambipolarity in the cloud scales as Z -1/3 . Fluid simulations yield ablation rates that scale with the four-thirds power of the pellet radius in agreement with [2].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Sub- and super-luminar propagation of structures satisfying Poynting-like theorem for incompressible generalized hydrodynamic fluid model depicting strongly coupled dusty plasma medium

The strongly coupled dusty plasma has often been modelled by the Generalized Hydrodynamic (GHD) model used for representing visco-elastic fluid systems. The incompressible limit of the model which supports transverse shear wave mode is studied in detail. In particular, dipole structures are observed to emit transverse shear waves in both the limits of sub- and super-luminar propagation, where the structures move slower and faster than the phase velocity of the shear waves, respectively. In the sub-luminar limit the dipole gets engulfed within the shear waves emitted by itself, which then backreacts on it and ultimately the identity of the structure is lost. However, in the super-luminar limit the emission appears like a wake from the tail region of the dipole. The dipole, however, keeps propagating forward with little damping but minimal distortion in its form. A Poynting-like conservation law with radiative, convective, and dissipative terms being responsible for the evolution of W, which is similar to “enstrophy” like quantity in normal hydrodynamic fluid systems, has also been constructed for the incompressible GHD equations. The conservation law is shown to be satisfied in all the cases of evolution and collision amidst the nonlinear structures to a great accuracy. It is shown that monopole structures which do not move at all but merely radiate shear waves, the radiative term, and dissipative losses solely contribute to the evolution of W. The dipolar structures, on the other hand, propagate in the medium and hence convection also plays an important role in the evolution of W.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A nonlinear eigenvalue problem for self-similar spherical force-free magnetic fields

An axisymmetric force-free magnetic field B(r, θ) in spherical coordinates is defined by a function r sin θB{sub φ}=Q(A) relating its azimuthal component to its poloidal flux-function A. The power law r sin θB{sub φ}=aA|A|{sup 1/n}, n a positive constant, admits separable fields with A=(A{sub n}(θ))/(r{sup n}) , posing a nonlinear boundary-value problem for the constant parameter a as an eigenvalue and A{sub n}(θ) as its eigenfunction [B. C. Low and Y. Q Lou, Astrophys. J. 352, 343 (1990)]. A complete analysis is presented of the eigenvalue spectrum for a given n, providing a unified understanding of the eigenfunctions and the physical relationship between the field's degree of multi-polarity and rate of radial decay via the parameter n. These force-free fields, self-similar on spheres of constant r, have basic astrophysical applications. As explicit solutions they have, over the years, served as standard benchmarks for testing 3D numerical codes developed to compute general force-free fields in the solar corona. The study presented includes a set of illustrative multipolar field solutions to address the magnetohydrodynamics (MHD) issues underlying the observation that the solar corona has a statistical preference for negative and positive magnetic helicities in its northern and southern hemispheres, respectively; a hemispherical effect, unchanging as the Sun's global field reverses polarity in successive eleven-year cycles. Generalizing these force-free fields to the separable form B=(H(θ,φ))/(r{sup n+2}) promises field solutions of even richer topological varieties but allowing for φ-dependence greatly complicates the governing equations that have remained intractable. The axisymmetric results obtained are discussed in relation to this generalization and the Parker Magnetostatic Theorem. The axisymmetric solutions are mathematically related to a family of 3D time-dependent ideal MHD solutions for a polytropic fluid of index γ = 4/3 as discussed in the Appendix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Dependence on ion temperature of shallow-angle magnetic presheaths with adiabatic electrons

The magnetic presheath is a boundary layer occurring when magnetized plasma is in contact with a wall and the angle α between the wall and the magnetic field B is oblique. Here, we consider the fusion-relevant case of a shallow-angle, α << 1, electron-repelling sheath, with the electron density given by a Boltzmann distribution, valid for α/$\sqrt {τ + 1}$ >> $\sqrt {m_e/m_i}$, where m e is the electron mass, m i is the ion mass, τ = T i /ZT e , T e is the electron temperature, T i is the ion temperature, and Z is the ionic charge state. The thickness of the magnetic presheath is of the order of a few ion sound Larmor radii ρ s = p m i (ZT e + T i )/ZeB, where e is the proton charge and B = |B| is the magnitude of the magnetic field. We study the dependence on τ of the electrostatic potential and ion distribution function in the magnetic presheath by using a set of prescribed ion distribution functions at the magnetic presheath entrance, parameterized by τ. The kinetic model is shown to be asymptotically equivalent to Chodura’s fluid model at small ion temperature, τ << 1, for | ln α| > 3| ln τ| >> 1. In this limit, despite the fact that fluid equations give a reasonable approximation to the potential, ion gyro-orbits acquire a spatial extent that occupies a large portion of the magnetic presheath. At large ion temperature, τ>> 1, relevant because Ti is measured to be a few times larger than Te near divertor targets of fusion devices, ions reach the Debye sheath entrance (and subsequently the wall) at a shallow angle whose size is given by $\sqrt α$ or 1/ $\sqrt τ$ , depending on which is largest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Low-pressure hydrogen discharge maintenance in a large-size plasma source with localized high radio-frequency power deposition

The development of the two-dimensional fluid-plasma model of a low-pressure hydrogen discharge, presented in the study, is regarding description of the plasma maintenance in a discharge vessel with the configuration of the SPIDER source. The SPIDER source, planned for the neutral-beam-injection plasma-heating system of ITER, is with localized high RF power deposition to its eight drivers (cylindrical-coil inductive discharges) and a large-area second chamber, common for all the drivers. The continuity equations for the charged particles (electrons and the three types of positive ions) and for the neutral species (atoms and molecules), their momentum equations, the energy balance equations for electrons, atoms and molecules and the Poisson equations are involved in the discharge description. In addition to the local processes in the plasma volume, the surface processes of particle reflection and conversion on the walls as well as for a heat exchange with the walls are included in the model. The analysis of the results stresses on the role of the fluxes (particle and energy fluxes) in the formation of the discharge structure. The conclusion is that the discharge behavior is completely obeyed to non-locality. The latter is displayed by: (i) maximum values of plasma parameters (charged particle densities and temperatures of the neutral species) outside the region of the RF power deposition, (ii) shifted maxima of the electron density and temperature, of the plasma potential and of the electron production, (iii) an electron flux, with a vortex structure, strongly exceeding the total ion flux which gives evidence of a discharge regime of non-ambipolarity and (iv) a spatial distribution of the densities of the neutral species resulting from their fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

TWO-DIMENSIONAL RADIATIVE MAGNETOHYDRODYNAMIC SIMULATIONS OF THE IMPORTANCE OF PARTIAL IONIZATION IN THE CHROMOSPHERE

The bulk of the solar chromosphere is weakly ionized and interactions between ionized particles and neutral particles likely have significant consequences for the thermodynamics of the chromospheric plasma. We investigate the importance of introducing neutral particles into the MHD equations using numerical 2.5D radiative MHD simulations obtained with the Bifrost code. The models span the solar atmosphere from the upper layers of the convection zone to the low corona, and solve the full MHD equations with non-gray and non-LTE radiative transfer, and thermal conduction along the magnetic field. The effects of partial ionization are implemented using the generalized Ohm's law, i.e., we consider the effects of the Hall term and ambipolar diffusion in the induction equation. The approximations required in going from three fluids to the generalized Ohm's law are tested in our simulations. The Ohmic diffusion, Hall term, and ambipolar diffusion show strong variations in the chromosphere. These strong variations of the various magnetic diffusivities are absent or significantly underestimated when, as has been common for these types of studies, using the semi-empirical VAL-C model as a basis for estimates. In addition, we find that differences in estimating the magnitude of ambipolar diffusion arise depending on which method is used to calculate the ion-neutral collision frequency. These differences cause uncertainties in the different magnetic diffusivity terms. In the chromosphere, we find that the ambipolar diffusion is of the same order of magnitude or even larger than the numerical diffusion used to stabilize our code. As a consequence, ambipolar diffusion produces a strong impact on the modeled atmosphere. Perhaps more importantly, it suggests that at least in the chromospheric domain, self-consistent simulations of the solar atmosphere driven by magnetoconvection can accurately describe the impact of the dominant form of resistivity, i.e., ambipolar diffusion. This suggests that such simulations may be more realistic in their approach to the lower solar atmosphere (which directly drives the coronal volume) than previously assumed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2012

Magnetic Field Effects on Plasma Plumes

Here, we will discuss our numerical studies of plasma jets and loops, of basic interest for plasma propulsion and plasma astrophysics. Space plasma propulsion systems require strong guiding magnetic fields known as magnetic nozzles to control plasma flow and produce thrust. Propulsion methods currently being developed that require magnetic nozzles include the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR) [1] and magnetoplasmadynamic thrusters. Magnetic nozzles are functionally similar to de Laval nozzles, but are inherently more complex due to electromagnetic field interactions. The two crucial physical phenomenon are thrust production and plasma detachment. Thrust production encompasses the energy conversion within the nozzle and momentum transfer to a spacecraft. Plasma detachment through magnetic reconnection addresses the problem of the fluid separating efficiently from the magnetic field lines to produce maximum thrust. Plasma jets similar to those of VASIMR will be studied with particular interest in dual jet configurations, which begin as a plasma loops between two nozzles. This research strives to fulfill a need for computational study of these systems and should culminate with a greater understanding of the crucial physics of magnetic nozzles with dual jet plasma thrusters, as well as astrophysics problems such as magnetic reconnection and dynamics of coronal loops.[2] To study this problem a novel, hybrid kinetic theory and single fluid magnetohydrodynamic (MHD) solver known as the Magneto-Gas Kinetic Method is used.[3] The solver is comprised of a "hydrodynamic" portion based on the Gas Kinetic Method and a "magnetic" portion that accounts for the electromagnetic behaviour of the fluid through source terms based on the resistive MHD equations. This method is being further developed to include additional physics such as the Hall effect. Here, we will discuss the current level of code development, as well as numerical simulation results

Ebersohn, F.↗
DOE OSTI2019

Ellipticity conditions for the extended MHD Grad-Shafranov-Bernoulli equilibrium equations

In this study, we find the points of transition between elliptic and hyperbolic regimes for the axisymmetric extended magnetohydrodynamic (MHD) equilibrium equations. The ellipticity condition is expressed via a single inequality but is more involved than the corresponding two-fluid ones due to the imposition of the quasineutrality condition and is also more complicated than the Hall MHD ellipticity condition, due to electron inertia. In fact, the inclusion of electron inertia is responsible for peculiar results; namely, even the static equilibrium equations can become hyperbolic.

Speed of sound↗
DOE OSTI2023

Improved melt model for power flow

Accelerators that drive z-pinch experiments transport current densities in excess of 1 MA/cm 2 in order to melt or ionize the target and implode it on axis. These high current densities stress the transmission lines upstream from the target, where rapid electrode heating causes plasma formation, melt, and possibly vaporization. These plasmas negatively impact accelerator efficiency by diverting some portion of the current away from the target, referred to as “current loss”. Simulations that are able to reproduce this behavior may be applied to improving the efficiency of existing accelerators and to designing systems operating at ever higher current densities. The relativistic particle-in-cell code CHICAGO ® is the primary code for modeling power flow on Sandia National Laboratories’ Z accelerator. We report here on new algorithms that incorporate vaporization and melt into the standard power-flow simulation framework. Taking a hybrid approach, the CHICAGO® kinetic/multi-fluid treatment has been expanded to include vaporization while the quasi-neutral equation-of-motion has been updated for melt at high current-densities. For vaporization, a new one-dimensional substrate model provides a more accurate calculation of electrode thermal, mass, and magnetic field diffusion as well as a means of emitting absorbed contaminants and vaporized metal ions. A quasi-fluid model has been implemented expressly to mimic the motion of imploding liners for accurate inductance histories. For melt, a multi-ion Hall-MHD option has been implemented and benchmarked against Alegra MHD. This new model is described with sufficient detail to reproduce these algorithms in any hybrid kinetic code. Physics results from the new code are also presented. A CHICAGO ® Hall-MHD simulation of a radial transmission line demonstrates that Hall physics, not included in Alegra, has no significant impact on the diffusion of electrode material. When surface contaminant desorption is mocked in as a hydrogen surface plasma, both the surface and bulk-material plasmas largely compress under the influence of the j × B force. Similar results are seen in Alegra, which also shows magnetic and material diffusion scaling with peak current. Test vaporization simulations using MagLIF and a power-flow experimental geometry show Fe + ions diffuse only a few hundred µm from the electrodes, so present models of Z power flow remain valid.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

A tutorial introduction to the statistical theory of turbulent plasmas, a half-century after Kadomtsev’s Plasma Turbulence and the resonance-broadening theory of Dupree and Weinstock

In honour of the 50th anniversary of the influential review/monograph on plasma turbulence by B. B. Kadomtsev as well as the seminal works of T. H. Dupree and J. Weinstock on resonance-broadening theory, an introductory tutorial is given about some highlights of the statistical–dynamical description of turbulent plasmas and fluids, including the ideas of nonlinear incoherent noise, coherent damping, and self-consistent dielectric response. The statistical closure problem is introduced. Incoherent noise and coherent damping are illustrated with a solvable model of passive advection. Self-consistency introduces turbulent polarization effects that are described by the dielectric function ${\mathcal{D}}$ . Dupree’s method of using ${\mathcal{D}}$ to estimate the saturation level of turbulence is described; then it is explained why a more complete theory that includes nonlinear noise is required. The general theory is best formulated in terms of Dyson equations for the covariance $C$ and an infinitesimal response function $R$ , which subsumes ${\mathcal{D}}$ . An important example is the direct-interaction approximation (DIA). It is shown how to use Novikov’s theorem to develop an $\boldsymbol{x}$ -space approach to the DIA that is complementary to the original $\boldsymbol{k}$ -space approach of Kraichnan. A dielectric function is defined for arbitrary quadratically nonlinear systems, including the Navier–Stokes equation, and an algorithm for determining the form of ${\mathcal{D}}$ in the DIA is sketched. The independent insights of Kadomtsev and Kraichnan about the problem of the DIA with random Galilean invariance are described. The mixing-length formula for drift-wave saturation is discussed in the context of closures that include nonlinear noise (shielded by ${\mathcal{D}}$ ). The role of $R$ in the calculation of the symmetry-breaking (zonostrophic) instability of homogeneous turbulence to the generation of inhomogeneous mean flows is addressed. The second-order cumulant expansion and the stochastic structural stability theory are also discussed in that context. In conclusion, various historical research threads are mentioned and representative entry points to the literature are given. In addition, some outstanding conceptual issues are enumerated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗