Landau Fluid Research Articles

Study of the evolution of density and particle transport during density relaxation phase after pellet fueling on the EAST tokamak

Abstract The impact of pellet injection depth on particle transport characteristics and density profile evolution during the post-pellet density relaxation phase is investigated on the experimental advanced superconducting tokamak. Experiments were conducted under similar background plasma conditions and the results showed that variations in pellet deposition depths led to differences in density profiles despite comparable core fueling efficiencies. By computing the effective particle transport coefficients immediately after the pellet injection, it is found that the differences in density profiles are mainly due to significantly different transport velocities in the inner region of ρ < 0.4. The potential influence of pellets on electrostatic turbulence in this region is investigated using the trapped gyro Landau fluid model and CO 2 laser collective scattering diagnostics. The simulation suggests that the increase in transport coefficients following pellet injection is attributed to the destabilization of lower-k modes rather than the higher-k modes within the diagnostic measurement range.

Ti/Te effects on transport in EAST low q95 plasmas

Enhanced confinement is observed in EAST low q95 ( q95<3.5 ) plasmas with the temperature ratio ( Ti/Te ) increasing. A statistically significant positive correlation between H98y2 and Ti/Te has been established. Increased Ti triggerd by Neutral Beam Injection in experiments maintaining constant Te leads initially to an improved Ti/Te ratio, which subsequently mitigates heat transport in the ion channel. The Trapped Gyro–Landau Fluid transport model was applied to analyse the dependence of the turbulent transport on Ti/Te . Simulation results demonstrate that the Ion Temperature Gradient mode is stabilized in high Ti/Te plasma, with ion flux results from predictive modelling further validating this effect. These findings corroborate the Weiland model and complement one another. This study contributes to a deeper understanding of how the dimensionless parameter Ti/Te influences plasma transport, and helps demonstrate the feasibility of operations under the low q95 condition required for the ITER-inductive scenario in the new Ti/Te regime.

Validation of D–T fusion power prediction capability against 2021 JET D–T experiments

JET experiments using the fuel mixture envisaged for fusion power plants, deuterium and tritium (D–T), provide a unique opportunity to validate existing D–T fusion power prediction capabilities in support of future device design and operation preparation. The 2021 JET D–T experimental campaign has achieved D–T fusion powers sustained over 5 s in ITER-relevant conditions i.e. operation with the baseline or hybrid scenario in the full metallic wall. In preparation of the 2021 JET D–T experimental campaign, extensive D–T predictive modelling was carried out with several assumptions based on D discharges. To improve the validity of ITER D–T predictive modelling in the future, it is important to use the input data measured from 2021 JET D–T discharges in the present core predictive modelling, and to specify the accuracy of the D–T fusion power prediction in comparison with the experiments. This paper reports on the validation of the core integrated modelling with TRANSP, JINTRAC, and ETS coupled with a quasilinear turbulent transport model (Trapped Gyro Landau Fluid or QualLiKiz) against the measured data in 2021 JET D–T discharges. Detailed simulation settings and the heating and transport models used are described. The D–T fusion power calculated with the interpretive TRANSP runs for 38 D–T discharges (12 baseline and 26 hybrid discharges) reproduced the measured values within 20. This indicates the additional uncertainties, that could result from the measurement error bars in kinetic profiles, impurity contents and neutron rates, and also from the beam-thermal fusion reaction modelling, are less than in total. The good statistical agreement confirms that we have the capability to accurately calculate the D–T fusion power if correct kinetic profiles are predicted, and indicates that any larger deviation of the D–T fusion power prediction from the measured fusion power could be attributed to the deviation of the predicted kinetic profiles from the measured kinetic profiles in these plasma scenarios. Without any posterior adjustment of the simulation settings, the ratio of predicted D–T fusion power to the measured fusion power was found as 65%–96% for the D–T baseline and 81%–97% for D–T hybrid discharge. Possible reasons for the lower D–T prediction are discussed and future works to improve the fusion power prediction capability are suggested. The D–T predictive modelling results have also been compared to the predictive modelling of the counterpart D discharges, where the key engineering parameters are similar. Features in the predicted kinetic profiles of D–T discharges such as underprediction of ne are also found in the prediction results of the counterpart D discharges, and it leads to similar levels of the normalized neutron rate prediction between the modelling results of D–T and the counterpart D discharges. This implies that the credibility of D–T fusion power prediction could be a priori estimated by the prediction quality of the preparatory D discharges, which will be attempted before actual D–T experiments.

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes’ contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes’ contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and ρ⩽0.7 in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

Global Effects on Drift Wave Microturbulence in Tokamak Plasmas

Microturbulence excited by ion temperature gradient (ITG)-dominant and trapped electron mode (TEM)-dominant instabilities is investigated by employing an extended fluid code (ExFC) based on the so-called Landau fluid model, which includes the trapped electron dynamics. Firstly, the global effect is emphasized through direct comparison of ITG and TEM instability domains based on local and global simulations. The global effect makes differences in both linear instability and nonlinear transport, including the fluxes and the structure of zonal flow. The transitions among ITG, TEM, and ITG & TEM (ITG & TEM represents that ITG and TEM coexist with different wavelengths) instabilities/turbulence depend not only on the three key drive forces (R/L n, R/L Te, R/L Ti) but also on their global (profile) effects. Secondly, a lot of electrostatic linear gyro-fluid simulations are concluded to obtain a distribution of the instability.

Electromagnetic drift wave instability in tokamak plasmas with strong pedestal gradient

The linear eigenmode characterizations and the nonlinear turbulence energy spreading of the drift waves in a tokamak plasma with strong pedestal gradient are numerically investigated based on an electromagnetic Landau fluid model. By the linear eigenmode analysis, it is found that the dominant instability in the low β regime is the ion-temperature-gradient (ITGc) mode and the electron drift wave instability (eDWI p ) in the core and edge region with strong density gradient, respectively. Multiple eigenstates of the eDWI p with different peak locations in the poloidal direction can be obtained by the eigenvalue problem solver. The dominant one is the high order eDWI p corresponding to the unconventional ballooning mode structure with multiple peaks in the poloidal position, in contrast to the conventional modes that peak at the outboard mid-plane, and has been verified through initial value simulation. In the high β regime, the dominant eigenmodes in the core and edge region are the conventional and unconventional kinetic ballooning modes respectively. In the nonlinear simulation, an inward turbulence spreading phenomenon during the quasi-saturation phase of the edge turbulence is clearly observed. The inward speed of the turbulence energy front in the high β regime is much faster than that in the low β regime. It is interestingly found that the speed of the turbulence energy front increases with the increase of the plasma β in the low β regime, while it is almost unchanged in the high β regime. It is identified that the turbulence spreading in the low and high β regimes are determined by the nonlinear dynamics and the linear toroidal coupling respectively.

Global simulations of energetic electron excitation of beta-induced Alfvén eigenmodes

The energetic electron (EE) excitation of beta-induced Alfvén eigenmodes is investigated by using the newly developed global eigenvalue code MAS, which is based on a hybrid model that consists of Landau fluid bulk plasma and drift kinetic EE. Specifically, the bulk plasma kinetic effects such as finite Larmor radius, diamagnetic drifts and Landau dampings, and the EE adiabatic fluid response of convection and non-adiabatic kinetic response of precessional drift resonance are incorporated in the simulations. The global eigenmode equation is solved for e-BAE mode structure and linear dispersion relation in tokamak non-perturbatively. The radial width of e-BAE mode structure becomes narrower as the toroidal mode number increases, which can be explained by the change of Alfvén continuous spectra that interact with kinetic Alfvén waves for corresponding eigenmode formation. The e-BAE growth rate exhibits a non-monotonic variation with toroidal mode number for precessional drift resonance destabilization, while the e-BAE real frequency is close to the continuum accumulation point that almost remains the same. The parametric dependence of e-BAE stability on EE density and that on temperature are analyzed by MAS non-perturbative simulations, which shows that the EE density can affect e-BAE real frequency and thus changes the resonance condition, resulting in e-BAE stabilization in the strong EE drive regime. Further, the EE non-perturbative effect on the symmetry breaking of e-BAE mode structure is reported. The poloidal symmetry breaking characterized by the ‘boomerang’ shape two-dimensional (2D) structure can be greatly enhanced by increasing EE temperature, together with the large radial variation of the poloidal phase angle of dominant principal poloidal harmonic. The radial symmetry breaking of e-BAE mode structure arises when EE density/temperature drive is not symmetric with respect to corresponding rational surface, which can lead to a net volume-averaged value of e-BAE parallel wave number which drives plasma intrinsic rotation. These results are helpful in understanding the e-BAE dynamics observed in recent experiments.

Linear simulation of kinetic Peeling-Ballooning mode with bootstrap current under the BOUT++ gyro landau fluid code

The kinetic peeling-ballooning mode (KPBM) plays a crucial role in the edge turbulence and transport in a tokamak plasma. However, the impact of the bootstrap current on KPBM is still unclear. Simulations of KPBM using the BOUT++ gyro-Landau-fluid code are presented in this study. To investigate the KPBM in a real tokamak equilibrium, the global equilibrium solver CORSICA was employed to generate a set of realistic equilibria of shifted circular geometry, including the Shafranov shift, elongation effects and bootstrap current. The linear instability property of the KPBM is observed in a wide range of pressure gradient and parallel current density in the pedestal region. Our results indicate that the bootstrap current has a stabilizing effect on the high-n KPBM over the entire pedestal region. Here, n is the toroidal mode number. In comparison with the ideal magnetohydrodynamic peeling-ballooning mode stability diagram, the unstable KPBM region shrinks and shifts to low ( where denotes the pressure at the peak pressure gradient position) with the increase of the bootstrap current. Moreover, we find that the low-n kink modes are driven unstable by the bootstrap current in the second stable region. With an additional external current drive on top of the bootstrap current, the low-n KPBMs can be stabilized when the total edge current is sufficiently large, although the low-n kink mode is still unstable in the region with .

Explaining the lack of power degradation of energy confinement in wide pedestal quiescent H-modes via transport modeling

Wide pedestal quiescent H (WPQH)-mode is an attractive scenario for future burning plasmas as they operate without ELMs. WPQH is characterized by formation of a wider and higher pedestal (than quiescent H-mode), and broadband fluctuations in the pedestal. Unlike conventional H-modes, where the energy confinement time reduces with increasing heating power, the WPQH plasmas reported in this paper do not show power degradation of the energy confinement. As the injected neutral beam power was increased, reduced core (ρ ⩽ 0.45) transport calculated by transp, as well as increased core temperatures, pressure gradient and diamagnetic E × B shear rate were observed. The reduction in the heat transport and rapid decrease in the ion temperature gradient scale length suggest the formation of an ion internal transport barrier (ITB) that was accompanied by increased stored energy in the core. Quasilinear turbulent transport modeling using the trapped gyro Landau fluid (tglf) code was used to predict the ITB and its turbulence stability properties. By using profiles and equilibria produced by matching the transp transport fluxes with the tglf transport model within the tgyro transport solver, the energy confinement time captures the experimentally observed insensitivity to the increased P NBI. Linear stability analysis reveals that drift-wave instabilities in the core are stabilized by E × B shear, T i/T e ratio and Shafranov shift; the latter was found to have the strongest effect on the turbulence suppression at the highest heating level.

Simulation prediction of micro-instability transition and associated particle transport in tokamak plasmas

Two reduced simulation approaches are exploited to predict the parametric boundary of dominant instability regime with global effects and the characteristics of corresponding turbulent particle fluxes in tokamak plasmas. One is usual numerical simulation of coexisting ion temperature gradient (ITG) mode and trapped electron mode (TEM) turbulence employing an extended fluid code (ExFC) based on the so-called Landau–Fluid model including the trapped electron dynamics. Here the density gradient (i.e. R/L n ) driven TEM (∇n-TEM) is emphasized. The other one is a surrogate turbulence transport model, taking a neural network (NN) based approach with speeding calculation. It is shown that the turbulent particle flux, particularly their directions depend on the type of micro-instability as ITG and/or TEM. On the other hand, the density gradient may govern the direction of the turbulent particle fluxes in general circumstances. Specifically, in the parameter regime explored here, the ITG and the electron temperature gradient driven TEM (∇T e-TEM) are destabilized for flat density profile, generally causing an inward particle flux, i.e. particle pinch. Contrarily, for steep density profile, the ∇n-TEM or coexisting ITG and TEM turbulence are dominant so that the particle always diffuses outwards. An empirical criterion is obtained to predict the dominant instability and the direction of particle flux for medium density gradients, involving the gradients of both ion and electron temperature as well as the density. These two transport models are applied to analyze the spontaneous excitation of a quasi-coherent mode in the turbulence modulation discharge by MHD magnetic island observed on tokamak HL-2A, clearly showing a dynamic transition from ITG to TEM. Furthermore, the ExFC-NN model can predict and speed up the analysis of the turbulence transport in tokamak experiments.

Drift reduced Landau fluid model for magnetized plasma turbulence simulations in BOUT++ framework

Recently the drift-reduced Landau fluid six-field turbulence model within the BOUT++ framework [1] has been upgraded. In particular, this new model employs a new normalization, adds a volumetric flux-driven source option, the Landau fluid closure for parallel heat flux and a Laplacian inversion solver which is able to capture n=0 axisymmetric mode evolution in realistic tokamak configurations. These improvements substantially extended model's capability to study a wider range of tokamak edge phenomena, and are essential to build a fully self-consistent edge turbulence model capable of both transient (e.g., ELM, disruption) and transport time-scale simulations.

A non-local fluid closure for modeling cyclotron resonance in collisionless magnetized plasmas

A fluid description for collisionless magnetized plasmas that takes into account the effect of cyclotron resonance has been developed. Following the same approach as the Landau fluid closure, the heat flux components associated with transverse electromagnetic fluctuations are approximated by a linear combination of lower-order moments in wavenumber space. The closure successfully reproduces the linear cyclotron resonance for electromagnetic waves propagating parallel to the ambient magnetic field. In the presence of finite temperature anisotropy, the model gives approximately correct prediction for an instability destabilized via the cyclotron resonance. A nonlinear simulation demonstrates the wave growth consistent with the linear theory followed by the reduction of initial anisotropy, and finally, the saturation of the instability. The isotropization may be understood in terms of quasilinear theory, which is developed within the framework of the fluid model but very similar to its fully kinetic counterpart. The result indicates that both linear and nonlinear collisionless plasma responses are approximately incorporated in the fluid model.

Local and global properties of energy transfer in models of plasma turbulence

The nature of the turbulent energy transfer rate is studied using direct numerical simulations of weakly collisional space plasmas. This is done comparing results obtained from hybrid Vlasov–Maxwell simulations of collisionless plasmas, Hall magnetohydrodynamics and Landau fluid models reproducing low-frequency kinetic effects, such as the Landau damping. In this turbulent scenario, estimates of the local and global scaling properties of different energy channels are obtained using a proxy of the local energy transfer. This approach provides information on the structure of energy fluxes, under the assumption that the turbulent cascade transfers most of the energy that is then dissipated at small scales by various kinetic processes in these kinds of plasmas.

Deep learning surrogate model for kinetic Landau-fluid closure with collision

In this work, the kinetic Landau-fluid (LF) closure with collision and periodic boundary condition is used in the development of the deep learning (DL) surrogate model. A classical neural network, namely, feedforward neural network or sometimes termed multilayer perceptron, is constructed and trained to learn the kinetic LF closure in the static limit and arbitrary mean free path in configuration space. The preliminary relation between best hyperparameters and critical parameters for data generation is found. Compared with the numerical approach (non-Fourier method) of the LF closure, the deep learning surrogate model shows an order of magnitude of improvement in terms of accuracy. Perhaps most importantly, the surrogate model closure has been integrated for the first time with fluid simulations. Our DL-enabled fluid simulations, for the first time, give the correct Landau damping rate for a wide range of wave vectors, while the Hammett–Perkins closure cannot produce the correct damping rate. We correctly connect the collisionless Hammett–Perkins closure and collisional Braginskii closure to reproduce the intrinsic nonlocal feature of the heat flux with DL techniques. We address the most concerning error accumulation problem and find that simulations with the deep learning surrogate model are as good as, if not better than, simulations with the analytic closure in terms of long-term numerical stability in the linear Landau damping test.

Machine learning surrogate models for Landau fluid closure

The first result of applying the machine/deep learning technique to the fluid closure problem is presented in this paper. As a start, three different types of neural networks [multilayer perceptron (MLP), convolutional neural network (CNN), and two-layer discrete Fourier transform (DFT) network] were constructed and trained to learn the well-known Hammett–Perkins Landau fluid closure in configuration space. We find that in order to train a well-preformed network, a minimum size of the training data set is needed; MLP also requires a minimum number of neurons in the hidden layers that equals the degrees of freedom in Fourier space, despite the fact that training data are being fed into the configuration space. Out of the three models, DFT performs the best for the clean data, most likely due to the existence of the simple Fourier expression for the Hammett–Perkins closure, but it is the least robust with respect to input noise. Overall, with appropriate tuning and optimization, all three neural networks are able to accurately predict the Hammett–Perkins closure and reproduce the intrinsic nonlocal feature, suggesting a promising path to calculating more sophisticated closures with the machine/deep learning technique.

An introductory guide to fluid models with anisotropic temperatures. Part 1. CGL description and collisionless fluid hierarchy

We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modelling the turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non-Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the fluid model of Chew–Goldberger–Low (CGL) in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the fourth-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the ‘normal’ closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the ‘normal’ closure, and beyond the fourth-order moment, Landau fluid closures are required.

An introductory guide to fluid models with anisotropic temperatures. Part 2. Kinetic theory, Padé approximants and Landau fluid closures

In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(\unicode[STIX]{x1D701})$ , and the associated plasma response function $R(\unicode[STIX]{x1D701})=1+\unicode[STIX]{x1D701}Z(\unicode[STIX]{x1D701})=-Z^{\prime }(\unicode[STIX]{x1D701})/2$ . We then consider a one-dimensional (1-D) (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the fourth-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1-D closures at higher-order moments than the fourth order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3-D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the fourth-order moment level. In appendix A, we provide comprehensive tables with Padé approximants of $R(\unicode[STIX]{x1D701})$ up to the eighth-pole order, with many given in an analytic form.

Basic features of the multiscale interaction between tearing modes and slab ion-temperature-gradient modes**Project supported by the National Key R&D Program of China (Grant Nos. 2017YFE0301100 and 2017YFE0300500), the National Natural Science Foundation of China (Grant Nos. 11675038, 11775069, and 11305027), and the Fundamental Research Funds for the Central Universities of China (Grant No. DUT17RC(4)54).

Nonlinear interaction between tearing modes (TM) and slab ion-temperature-gradient (ITG) modes is numerically investigated by using a Landau fluid model. It is observed that the energy spectra with respect to wavenumbers become broader during the transition phase from the ITG-dominated stage to TM-dominated stage. Accompanied with the fast growth of the magnetic island, the frequency of TM/ITG with long/short wavelength fluctuations in the electron/ion diamagnetic direction decreases/increases respectively. The decrease of TM frequency is identified to result from the effect of the profile flattening in the vicinity of the magnetic island, while the increase of the frequencies of ITG fluctuations is due to the eigenmode transition of ITG induced by the large scale zonal flow and zonal current related to TM. Roles of zonal current induced by the ITG fluctuations in the instability of TM are also analyzed. Finally, the electromagnetic transport features in the vicinity of the magnetic island are discussed.

Perturbative transport modeling of cold-pulse dynamics in Alcator C-Mod Ohmic plasmas

Perturbative transport experiments in magnetically confined plasmas have shown, for more than 20 years, that the injection of cold pulses at the plasma edge can trigger the increase of core temperature. Predictive heat transport simulations with the trapped gyro Landau fluid (TGLF) quasilinear transport model demonstrate that the increase of core temperature in some regimes, and lack thereof in other regimes, can be explained by a change in dominant linear micro-instability in Alcator C-Mod. The effect of density and plasma current on the cold pulse are well captured by TGLF, including the relative change in position of the temperature flex point as current density changes. Linear stability analysis of simulated density and current scans reveals a competition between trapped electron and ion temperature gradient modes as the main driver of the core transient response. These results further demonstrate that cold-pulse propagation and associated phenomenology in the cases studied are well explained within the local transport paradigm, without resorting to non-local effects.

Drift Instabilities in Thin Current Sheets Using a Two‐Fluid Model With Pressure Tensor Effects

AbstractThe integration of kinetic effects in fluid models is important for global simulations of the Earth's magnetosphere. We use a two‐fluid 10‐moment model, which includes the pressure tensor and has been used to study reconnection, to study the drift kink and lower hybrid drift instabilities. Using a nonlocal linear eigenmode analysis, we find that for the kink mode, the 10‐moment model shows good agreement with kinetic calculations with the same closure model used in reconnection simulations, while the electromagnetic and electrostatic lower hybrid instabilities require modeling the effects of the ion resonance using a Landau fluid closure. Comparisons with kinetic simulations and the implications of the results for global magnetospheric simulations are discussed.