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 829 records · Page 46

Effective viscosities in a hydrodynamically expanding boost-invariant QCD plasma

Background: The near-equilibrium properties of a QCD plasma can be encoded into transport coefficients such as bulk and shear viscosity. In QCD, the ratio of these transport coefficients to entropy density, ζ/$\textit{s}$ and η/$\textit{s}$, depends nontrivially on the plasma's temperature. This is unlike in conformal systems where they take constant values such as η/$\textit{s}$ = 1/(4π). Purpose: In this work, we show that in a (0 + 1)D boost-invariant fluid with no transverse expansion, a temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ can be described by an equivalent effective viscosity $\langleζ/s\rangle_{\text{eff}}$ or $\langleη/s\rangle_{\text{eff}}$. This effective viscosity combines the actual temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ with the temperature profile of the fluid. Further, we extend the concept of effective viscosity in systems with transverse expansion and discuss how effective viscosities can be used to identify families of ζ/$\textit{s(T)}$ and η /$\textit{s(T)}$ that lead to similar hydrodynamic evolution. Methods: The Navier-Stokes relativistic hydrodynamic equations are used to provide a first definition of effective viscosity, in (0 + 1)D and (1 + 1)D. In the (0 + 1)D case, the analysis is extended to Israel-Stewart-type second-order hydrodynamics to clarify the effect of higher-order hydrodynamics corrections on the effective viscosity. Results: In a boost-invariant fluid with no transverse expansion [(0 + 1)D], the effective viscosity is expressed as a simple integral of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ over temperature, with a weight determined by the speed of sound of the fluid. The result is general for any equation of state with a moderate temperature dependence of the speed of sound, including the QCD equation of state. This definition of effective viscosity can be used to identify infinite families of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ that produce essentially indistinguishable temperature profiles. In a boost-invariant cylindrical system [(1 + 1)D], a similar definition of effective viscosity is obtained in terms of characteristic trajectories in time and transverse direction. This leads to an infinite number of constraints on an infinite functional space for ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$. Realistic examples are presented by using a finite number of constraints on a finite functional space. Conclusions: The definition of effective viscosity in a (0 + 1)D system clarifies how infinite families of ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$ can result in nearly identical hydrodynamic temperature profiles. By extending the study to a boost-invariant cylindrical [(1 + 1)D] fluid, we identify an approximate but more general definition of effective viscosity that highlights the potential and limits of the concept of effective viscosity in fluids with limited symmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum simulations of hydrodynamics via the Madelung transformation

Developing numerical methods to simulate efficiently nonlinear fluid dynamics on universal quantum computers is a challenging problem. In this paper, a generalization of the Madelung transform is defined to solve quantum relativistic charged fluid equations interacting with external electromagnetic forces via the Dirac equation. The Dirac equation is discretized into discrete-time quantum walks which can be efficiently implemented on universal quantum computers. A variant of this algorithm is proposed to implement simulations using current noisy intermediate scale quantum (NISQ) devices in the case of homogeneous external forces. High resolution (up to N=2 17 grid points) numerical simulations of relativistic and nonrelativistic hydrodynamical shocks on current IBM NISQs are performed with this algorithm. Here, this paper demonstrates that fluid dynamics can be simulated on NISQs, and opens the door to simulating other fluids, including plasmas, with more general quantum walks and quantum automata.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Areal density evolution of isolated surface perturbations at the onset of x-ray ablation Richtmyer-Meshkov growth

Isolated defects on inertial confinement fusion ignition capsules are a concern as defects taller than a few hundred nanometers are calculated to form jets of high-Z material, which enter the main fuel. If this mixing of high-Z material is not controlled, a serious degradation in thermonuclear burn can occur. A path towards controlling the growth of defects on the outer surface of plastic capsules is currently under development, but requires accurate predictions of defect evolution driven by the early time ablative Richtmyer-Meshkov (RM) effect. The chief uncertainty is the Equation of State (EOS) for polystyrene and its effect on ablative RM. We report on measurements of the growth of isolated defects made at the onset of ablative RM oscillations driven by x-ray ablation to differentiate between EOS models used in design calculations. Experiments at the OMEGA laser [T. R. Boehly et al., Opt. Commun. 133, 495 (1997)] used on-axis area backlighting radiography and x-ray framing cameras to determine bump areal densities at discrete times. Bumps 12 and 14 {mu}m tall and 33 {mu}m FWHM were found to grow to 2 x their initial areal density by 3 ns after the start of the drive laser pulse. Shock speed measurements established target conditions resulting from the ablation process. The tabular LEOS 5310 [D. Young and E. Corey, J. Appl. Phys. 78, 3748 (1995)] model shows good agreement with measured shock speeds and bump growth whereas the QEOS model [R. More et al., Phys. Fluids 31, 3059 (1988)] over predicts shock speed and under predicts bump growth by 6x. Differences in ablative RM behavior were also found for x-ray ablation compared to laser ablation, which result in an overestimation (or non-existence) of oscillation frequency for x-ray ablation as predicted by theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatio-temporal evolution and breaking of double layers: A description using Lagrangian hydrodynamics

The nonlinear development and collapse (breaking) of double layers in the long scale length limit is well described by equations for the cold ion fluid with quasineutrality. It is shown that electron dynamics is responsible for giving an 'equation of state' with negative ratio of specific heats to this fluid. Introducing a transformation for the density variable, the governing equation for the transformed quantity in terms of Lagrange variables turns out exactly to be a linear partial differential equation. This equation has been analyzed in various limits of interest. Nonlinear development of double layers with a sinusoidal initial disturbance and collapse of double layers with an initial perturbation in the form of a density void are analytically investigated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comet 1P/Halley multifluid MHD model for the Giotto fly-by

The interaction of comets with the solar wind has been the focus of many studies including numerical modeling. We compare the results of our multifluid MHD simulation of comet 1P/Halley to data obtained during the flyby of the European Space Agency's Giotto spacecraft in 1986. The model solves the full set of MHD equations for the individual fluids representing the solar wind protons, the cometary light and heavy ions, and the electrons. The mass loading, charge-exchange, dissociative ion-electron recombination, and collisional interactions between the fluids are taken into account. The computational domain spans over several million kilometers, and the close vicinity of the comet is resolved to the details of the magnetic cavity. The model is validated by comparison to the corresponding Giotto observations obtained by the Ion Mass Spectrometer, the Neutral Mass Spectrometer, the Giotto magnetometer experiment, and the Johnstone Plasma Analyzer instrument. The model shows the formation of the bow shock, the ion pile-up, and the diamagnetic cavity and is able to reproduce the observed temperature differences between the pick-up ion populations and the solar wind protons. We give an overview of the global interaction of the comet with the solar wind and then show the effects of the Lorentz force interaction between the different plasma populations.

79 ASTRONOMY AND ASTROPHYSICS↗

Four-dimensional QCD equation of state with multiple chemical potentials

Here, we construct a four-dimensional version of the equation of state (EoS) model neos, neos-4d, as a function of the temperature and chemical potentials of baryon, electric charge, and strangeness for the hot and dense quantum chromodynamics (QCD) matter created in relativistic nuclear collisions. This EoS enables multiple conserved charge current evolution in a relativistic fluid. Input from lattice QCD simulations and a hadron resonance gas model is considered for constructing the equation of state. We investigate its applicability to the relativistic hydrodynamic description of nuclear collisions and present a method for efficient numerical implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multifluid Block-Adaptive-Tree Solar Wind Roe-Type Upwind Scheme: Magnetospheric Composition and Dynamics During Geomagnetic Storms-Initial Results

The magnetosphere contains a significant amount of ionospheric O+, particularly during geomagnetically active times. The presence of ionospheric plasma in the magnetosphere has a notable impact on magnetospheric composition and processes. We present a new multifluid MHD version of the Block-Adaptive-Tree Solar wind Roe-type Upwind Scheme model of the magnetosphere to track the fate and consequences of ionospheric outflow. The multifluid MHD equations are presented as are the novel techniques for overcoming the formidable challenges associated with solving them. Our new model is then applied to the May 4, 1998 and March 31, 2001 geomagnetic storms. The results are juxtaposed with traditional single-fluid MHD and multispecies MHD simulations from a previous study, thereby allowing us to assess the benefits of using a more complex model with additional physics. We find that our multifluid MHD model (with outflow) gives comparable results to the multispecies MHD model (with outflow), including a more strongly negative Dst, reduced CPCP, and a drastically improved magnetic field at geosynchronous orbit, as compared to single-fluid MHD with no outflow. Significant differences in composition and magnetic field are found between the multispecies and multifluid approach further away from the Earth. We further demonstrate the ability to explore pressure and bulk velocity differences between H+ and O+, which is not possible when utilizing the other techniques considered

Glocer, A.↗

Comparison of kinetic and extended magnetohydrodynamics computational models for the linear ion temperature gradient instability in slab geometry

We perform linear stability studies of the ion temperature gradient (ITG) instability in unsheared slab geometry using kinetic and extended magnetohydrodynamics (MHD) models, in the regime k{sub ∥}/k{sub ⊥}≪1. The ITG is a parallel (to B) sound wave that may be destabilized by finite ion Larmor radius (FLR) effects in the presence of a gradient in the equilibrium ion temperature. The ITG is stable in both ideal and resistive MHD; for a given temperature scale length L{sub Ti0}, instability requires that either k{sub ⊥}ρ{sub i} or ρ{sub i}/L{sub Ti0} be sufficiently large. Kinetic models capture FLR effects to all orders in either parameter. In the extended MHD model, these effects are captured only to lowest order by means of the Braginskii ion gyro-viscous stress tensor and the ion diamagnetic heat flux. We present the linear electrostatic dispersion relations for the ITG for both kinetic Vlasov and extended MHD (two-fluid) models in the local approximation. In the low frequency fluid regime, these reduce to the same cubic equation for the complex eigenvalue ω=ω{sub r}+iγ. An explicit solution is derived for the growth rate and real frequency in this regime. These are found to depend on a single non-dimensional parameter. We also compute the eigenvalues and the eigenfunctions with the extended MHD code NIMROD, and a hybrid kinetic δf code that assumes six-dimensional Vlasov ions and isothermal fluid electrons, as functions of k{sub ⊥}ρ{sub i} and ρ{sub i}/L{sub Ti0} using a spatially dependent equilibrium. These solutions are compared with each other, and with the predictions of the local kinetic and fluid dispersion relations. Kinetic and fluid calculations agree well at and near the marginal stability point, but diverge as k{sub ⊥}ρ{sub i} or ρ{sub i}/L{sub Ti0} increases. There is good qualitative agreement between the models for the shape of the unstable global eigenfunction for L{sub Ti0}/ρ{sub i}=30 and 20. The results quantify how far fluid calculations can be extended accurately into the kinetic regime. We conclude that for the linear ITG problem in slab geometry with unsheared magnetic field when k{sub ∥}/k{sub ⊥}≪1, the extended MHD model may be a reliable physical model for this problem when ρ{sub i}/L{sub Ti0}<10{sup −2} and k{sub ⊥}ρ{sub i}<0.2.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microtearing modes in tokamak discharges

Microtearing modes (MTMs) have been identified as a source of significant electron thermal transport in tokamak discharges. In order to describe the evolution of these discharges, it is necessary to improve the prediction of electron thermal transport. This can be accomplished by utilizing a model for transport driven by MTMs in whole device predictive modeling codes. The objective of this paper is to develop the dispersion relation that governs the MTM driven transport. A unified fluid/kinetic approach is used in the development of a nonlinear dispersion relation for MTMs. The derivation includes the effects of electrostatic and magnetic fluctuations, arbitrary electron-ion collisionality, electron temperature and density gradients, magnetic curvature, and the effects associated with the parallel propagation vector. An iterative nonlinear approach is used to calculate the distribution function employed in obtaining the nonlinear parallel current and the nonlinear dispersion relation. The third order nonlinear effects in magnetic fluctuations are included, and the influence of third order effects on a multi-wave system is considered. An envelope equation for the nonlinear microtearing modes in the collision dominant limit is introduced in order to obtain the saturation level. In the limit that the mode amplitude does not vary along the field line, slab geometry, and strong collisionality, the fluid dispersion relation for nonlinear microtearing modes is found to agree with the kinetic dispersion relation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications

Here, the analysis of many fusion applications such as liquid-metal blankets requires application of computational fluid dynamics (CFD) methods for electrically conductive liquids in geometrically complex regions and in the presence of a strong magnetic field. A current state of the art general purpose CFD code allows modeling of the flow in complex geometric regions, with simultaneous conjugated heat transfer analysis in liquid and surrounding solid parts. Together with a magnetohydrodynamics (MHD) capability, the general purpose CFD code will be a valuable tool for the design and optimization of fusion devices. This paper describes an introduction of MHD capability into the general purpose CFD code CFX, part of the ANSYS Workbench. The code was adapted for MHD problems using a magnetic induction approach. CFX allows introduction of user-defined variables using transport or Poisson equations. For MHD adaptation of the code three additional transport equations were introduced for the components of the magnetic field, in addition to the Poisson equation for electric potential. The Lorentz force is included in the momentum transport equation as a source term. Fusion applications usually involve very strong magnetic fields, with values of the Hartmann number of up to tens of thousands. In this situation a system of MHD equations become very rigid with very large source terms and very strong variable gradients. To increase system robustness, special measures were introduced during the iterative convergence process, such as linearization using source coefficient for momentum equations. The MHD implementation in general purpose CFD code was tested against benchmarks, specifically selected for liquid-metal blanket applications. Results of numerical simulations using present implementation closely match analytical solutions for a Hartmann number of up to 1500 for a 2-D laminar flow in the duct of square cross section, with conducting and nonconducting walls. Results for a 3-D test case are also included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The kinetic effects of Alfven wave pressure in the solar wind

The results of a kinetic model for the radial evolution of the proton component of the solar wind in the presence of Alfven waves are presented. The calculation is based on general quasi-linear equations developed to describe the temporal and spatial evolution of the ion distribution functions of a multispecies plasma in presence of waves, using a short wavelength expansion. These equations include new wave-particle interaction terms arising from temporal and spatial inhomogeneities in the plasma. Numerical solutions are obtained of these equations specialized to the case of Alfven waves in a spherically symmetric solar wind. The Alfven wave effects on the proton distribution function vary strongly in velocity space. Protons with small transverse velocities are primarily decelerated with respect to the wave rest frame. This deceleration becomes less important with increasing transverse velocity, as wave induced diffusion to larger transverse velocity becomes the dominant effect. The competition of these effects results in interesting distortions of evolving proton distribution functions which give rise the wave acceleration well known from fluid theory.

Goodrich, C. C.↗

Self-Consistent Implementation of a Zero-Equation Transport Model Into a Predictive Model for a Hall Effect Thruster

The performance of an axisymmetric multi-fluid Hall thruster code that incorporates a self-consistent, data-driven closure model for the anomalous electron transport is investigated. Five different operating conditions of the H9 magnetically shielded Hall thruster are simulated with the Jet Propulsion Laboratory’s Hall2De. In order to capture the inherent uncertainty associated with the closure model, O(100) simulations are run for each condition, each of using a coefficient set sampled randomly from a probability distribution. The results of these simulations provide probabilistic predictions of thruster performance quantities including thrust, and discharge current, as well as several component efficiencies and centerline plasma properties. The model is found to yield converged solutions at all conditions, with large 10 kHzrange oscillations and performance trends with voltage and flow rate similar to experiment. The model under-predicts the thrust by 15-25% and over-predicts the discharge current by 20% on average compared to experiments at the same discharge voltage and mass flow rate. This performance discrepancy is due to lower beam utilization, mass utilization, and divergence efficiency than experiment, resulting from high Hall parameters in the acceleration region, which lead to a protracted ion acceleration region. The physical processes underlying this result are discussed in the context of future data-driven modeling efforts.

Jorns, Benjamin A.↗

Kinetic fluid moments closure for a magnetized plasma with collisions

A novel method aimed at a kinetic moments closure for a magnetized plasma with collisions is proposed. The velocity distribution function for each species is expanded in 8 Gaussian Radial Basis Functions (GRBFs) which are essentially shifted Maxwellians at eight representative 3D-velocity points of drift. The vector of 8 fluid moments (for particle density, 3 particle fluxes, total energy density, and 3 energy fluxes) has an 8x8 analytic linear matrix relation to the vector of 8 GRBF density weighs in 3D-real space. The 8 fluid moments with sources for each species are advanced in time while the 8 GRBF weighs are determined from the 8x8 inverse matrix. The two closure moments (for the stress tensor and the energy weighted stress tensor) are linearly determined from the GRBF weighs. Most importantly the velocity moments of the nonlinear Coulomb Fokker-Planck collision operator [Rosenbluth et al, Phys. Rev 107, 1957] are evaluated from the GRBF weights. Generalization from 8 to 12, 16, 20 .., in an energy weighted moment hierarchy is straightforward. The electric field follows from a generalized vorticity (quasi-neutral charge conservation) equation. A strong drift ordering approximation can be applied to eliminate any spuriously unstable high frequency cyclotron motions. A novel weak drift ordering two-time step scheme avoids the vorticity equation by following ion cyclotron motion in time to get the electric field with ion gyroaveraging. Inclusion of low-beta magnetic perturbations is straightforward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-dependent simulation of oblique MHD cosmic-ray shocks using the two-fluid model

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.

Frank, Adam↗

Using the maximum entropy distribution to describe electrons in reconnecting current sheets

Particle distributions in weakly collisional environments such as the magnetosphere have been observed to show deviations from the Maxwellian distribution. These can often be reproduced in kinetic simulations, but fluid models, which are used in global simulations of the magnetosphere, do not necessarily capture any of this. We apply the maximum entropy fluid closure of Levermore, which leads to well posed moment equations, to reconstruct particle distributions from a kinetic simulation in a reconnection region. Our results show that without information other than the moments, the model can reproduce the general structure of the distributions but not all of the finer details. In conclusion, the advantages of the closure over the traditional Grad closure are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A quasilinear operator retaining magnetic drift effects in tokamak geometry

The interaction of radio frequency waves with charged particles in a magnetized plasma is usually described by the quasilinear operator that was originally formulated by Kennel & Engelmann ( Phys. Fluids , vol. 9, 1966, pp. 2377–2388). In their formulation the plasma is assumed to be homogenous and embedded in a uniform magnetic field. In tokamak plasmas the Kennel–Engelmann operator does not capture the magnetic drifts of the particles that are inherent to the non-uniform magnetic field. To overcome this deficiency a combined drift and gyrokinetic derivation is employed to derive the quasilinear operator for radio frequency heating and current drive in a tokamak with magnetic drifts retained. The derivation requires retaining the magnetic moment to higher order in both the unperturbed and perturbed kinetic equations. The formal prescription for determining the perturbed distribution function then follows a novel procedure in which two non-resonant terms must be evaluated explicitly. The systematic analysis leads to a diffusion equation that is compact and completely expressed in terms of the drift kinetic variables. The equation is not transit averaged, and satisfies the entropy principle, while retaining the full poloidal angle variation without resorting to Fourier decomposition. As the diffusion equation is in physical variables, it can be implemented in any computational code. In the Kennel–Engelmann formalism, the wave–particle resonant delta function is either for the Landau resonance or the Doppler shifted cyclotron resonance. In the combined gyro and drift kinetic approach, a term related to the magnetic drift modifies the resonance condition.

Physics↗

A quasilinear operator retaining magnetic drift effects in tokamak geometry

The interaction of radio frequency waves with charged particles in a magnetized plasma is usually described by the quasilinear operator that was originally formulated by Kennel and Engelmann (Phys. Fluids vol. 9, 1966, pp. 2377-2388). In their formulation the plasma is assumed to be homogenous and embedded in a uniform magnetic field. In tokamak plasmas the Kennel-Engelmann operator does not capture the magnetic drifts of the particles that are inherent to the non-uniform magnetic field. To overcome this deficiency a combined drift and gyrokinetic derivation is employed to derive the quasilinear operator for radio frequency heating and current drive in a tokamak with magnetic drifts retained. The derivation requires retaining the magnetic moment to higher order in both the unperturbed and perturbed kinetic equations. The formal prescription for determining the perturbed distribution function then follows a novel procedure in which two non-resonant terms must be evaluated explicitly. The systematic analysis leads to a diffusion equation that is compact and completely expressed in terms of the drift kinetic variables. The equation is not transit averaged, and satisfies the entropy principle, while retaining the full poloidal angle variation without resorting to Fourier decomposition. As the diffusion equation is in physical variables, it can be implemented in any computational code. In the Kennel-Engelmann formalism, the wave-particle resonant delta function is either for the Landau resonance or the Doppler shifted cyclotron resonance. In the combined gyro and drift kinetic approach, a term related to the magnetic drift modifies the resonance condition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons

A conservative scheme has been formulated and confirmed for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from the drift kinetic equation is used to time advance the electron density perturbation by using the perturbed mechanical flow calculated from the parallel vector potential, and the parallel vector potential is solved by using the perturbed canonical flow from the perturbed distribution function. In gyrokinetic particle simulations using this new scheme, the shear Alfvén wave dispersion relation in the shearless slab and continuum damping in the sheared cylinder have been recovered. The new scheme overcomes the stringent requirement in the conventional perturbative simulation method that perpendicular grid size needs to be as small as electron collisionless skin depth even for the long wavelength Alfvén waves. The new scheme also circumvents the problem in the typical approach that an unphysically large parallel electric field arises due to the inconsistency between electrostatic potential calculated from the perturbed density and vector potential calculated from the perturbed canonical flow. Lastly, the gyrokinetic particle simulations of the Alfvén waves in sheared cylinder have superior numerical properties compared with the fluid simulations, which suffer from numerical difficulties associated with singular mode structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗