---
title: Trapped Gyro-Landau Fluid (TGLF) Model
url: https://www.emergentmind.com/topics/trapped-gyro-landau-fluid-tglf
type: topic
---

# Trapped Gyro-Landau Fluid (TGLF) Model

The Trapped Gyro-Landau Fluid (TGLF) model is a quasilinear, local gyro-fluid transport model for tokamak turbulence. It projects the five-dimensional gyrokinetic system onto a closed set of fluid-moment equations for each species, incorporates trapped-particle bounce averaging and Landau-fluid closures, and converts linear eigenmode spectra into turbulent particle, heat, and momentum fluxes through empirically calibrated saturation rules. In practical integrated modeling, TGLF is used to predict local anomalous transport from local gradients, geometry, collisionality, and finite-\(\beta\) parameters; a typical TGLF call takes \(\mathcal O(1\text{--}10)\,\mathrm{ms}\), but whole-device simulations requiring thousands of evaluations remain computationally expensive [2509.07024].

## 1. Foundational formulation and ordering assumptions

TGLF is described as a quasilinear, local gyro-fluid transport model whose basic assumptions include drift ordering, small-amplitude perturbations, a slab-like ballooning representation of toroidal geometry, and moment truncation with Landau-fluid closure. In the ST40 study, these assumptions are listed as
\[
\omega \lesssim \omega_*, \qquad k_\parallel \ll k_\perp, \qquad \rho\,k_\perp \lesssim 1,
\]
together with \(|\delta f| \ll f_0\) [2502.04993].

A complementary derivation writes the perturbed distribution for species \(s\) as
\[
f_s = F_{0s}(v) + \delta f_s(\mathbf R,v_\|,\mu),
\]
with \(F_{0s}\) a Maxwellian background depending on \(n_s(r)\) and \(T_s(r)\), and introduces the non-adiabatic response
\[
h_s=\delta f_s + (q_s\varphi/T_s)\,F_{0s}.
\]
Under the standard drift ordering \(\omega\sim k_\|\vth\ll\Omega_s\), with \(\Omega_s\) the cyclotron frequency and \(\vth=\sqrt{2T_s/m_s}\), the gyrokinetic dynamics are reduced to a fluid-moment hierarchy in which trapped particles are bounce averaged and passing populations are closed by Padé-type operators [2509.07024].

The significance of these approximations is operational rather than merely formal. TGLF is not a nonlinear gyrokinetic solver; it is a reduced-order model that retains the microinstability drives and transport-channel structure needed for transport prediction while remaining fast enough for repeated calls inside profile-evolution or flux-matching workflows.

## 2. Moment hierarchy, trapped-particle physics, and Landau-fluid closure

The model’s central reduction is from gyrokinetic phase space to a finite set of fluid moments. In the ST40 summary, TGLF is described as truncating at density, parallel momentum, and temperature moments, closed by Landau-fluid (“Hammett–Perkins”) operators. In the NSTX application summary, the retained state is described more explicitly as including density \(\delta n_s\), parallel velocity \(u_{\|s}\), parallel and perpendicular temperatures \(\delta T_{\|,s}\), \(\delta T_{\perp,s}\), and the corresponding parallel and perpendicular heat fluxes \(q_{\|,s}\), \(q_{\perp,s}\), giving up to eight moments for each species [2509.04359].

A compact expression of the linear system is
\[
\partial_t M + i \omega_\star M + L_{\rm geom}(M) + L_{\rm EM}(\delta A_\|,\delta B_\perp) = C(M) + S,
\]
where \(M\equiv (\delta n, u_\|, \delta T_\|, \delta T_\perp, q_\|, q_\perp)^{\mathsf T}\), \(\omega_\star\) contains diamagnetic drives proportional to \(R/L_T\) and \(R/L_n\), \(L_{\rm geom}\) encodes curvature and geometric couplings such as \(q\), \(\hat s\), and \(\epsilon\), and \(L_{\rm EM}\) incorporates finite-\(\beta\) electromagnetic effects [2509.04359].

Trapped-particle physics enters through bounce averaging and an effective Landau-pole closure. In the ST40 description, trapped electrons are handled by replacing the true \(k_\parallel\) in the Landau closure with an effective, bounce-averaged parallel wavenumber \(k_{\parallel,\rm eff}\), so that the trapped response captures the essence of trapped-electron Landau damping and bounce-resonant drive of trapped-particle modes such as TEM and UM. The NSTX summary likewise emphasizes the split between passing and trapped electrons, noting that trapped electrons carry their own fluid moments and respond non-adiabatically to curvature and \(\nabla B\) drifts [2502.04993].

The reduced hierarchy admits a quadratic free-energy description. One expression given for the decomposed free energy is
\[
\mathcal E \;=\;\sum_s \Bigl(\delta n_s^2 + m_s\,\delta u_{\|s}^2 + \tfrac{3}{2}n_s\,\delta T_s^2\Bigr)\;+\,\cdots,
\]
which identifies the density, parallel-flow, and temperature perturbations as the basic energetic content of the fluid reduction [2509.07024].

## 3. Quasilinear transport construction and computational workflow

TGLF computes transport from a linear eigenmode spectrum evaluated over binormal wavenumber. At each \(k_y\), linear analysis produces mode eigenfunctions \(\phi_{k_y}(x)\) and growth rate \(\gamma_{k_y}\). The quasilinear particle and heat fluxes are then assembled as spectral sums,
\[
\Gamma_s = \sum_{j=1}^{N_y}\Gamma_s^j(k_y^j), \qquad
Q_s = \sum_{j=1}^{N_y}Q_s^j(k_y^j),
\]
with per-mode contributions
\[
\Gamma_s^j = \frac{k_y^j}{B}\,\mathrm{Im}\bigl\langle\widetilde\phi^*_{k_y^j}\,\delta n_s^{k_y^j}\bigr\rangle,
\qquad
Q_s^j = \frac{k_y^j}{B}\,\mathrm{Im}\bigl\langle\widetilde\phi^*_{k_y^j}\,\delta T_s^{k_y^j}\bigr\rangle.
\]
A diffusive–pinch form is also used,
\[
\Gamma_s=-D_s\,\nabla n_s + V_s\,n_s, \qquad
Q_s=-\chi_s\,\nabla T_s + V_{T,s}\,T_s,
\]
where \(D_s\) and \(\chi_s\) are diffusivities and \(V_s\) and \(V_{T,s}\) are convective pinches obtained by projection of the quasilinear spectrum onto gradient directions [2509.07024].

In operational terms, a single native TGLF evaluation proceeds in three stages:

- **Inputs**: 31 variables, including local normalized gradient scale lengths \(a/L_{n_s}=-a\partial_r\ln n_s\) and \(a/L_{T_s}\), safety factor \(q\), magnetic shear \(\hat s\), effective charge \(Z_{\rm eff}\), elongation \(\kappa\), triangularity \(\delta\), normalized collisionality \(\nu^*\), beta \(\beta\), geometry factors \((r/R_0)\), and Shafranov shift.
- **Linear analysis**: solution of the bounce-averaged fluid eigenvalue problem for discrete \(k_y\)-values to obtain \(\gamma(k_y)\), real frequencies \(\omega_r(k_y)\), and quasilinear weight functions \(W_s(k_y)\).
- **Saturation and spectral summation**: application of empirically calibrated saturation amplitudes \(A(k_y)\), specifically the “SAT2” rule in the surrogate study, followed by
  \[
  \Gamma_s^j = A(k_y^j)\,W_s(k_y^j),\qquad
  Q_s^j = A(k_y^j)\,W_{T,s}(k_y^j),
  \]
  and summation over \(j\) to return \(\Gamma_e,\Pi_i,Q_e,Q_i\) [2509.07024].

The NSTX study emphasizes that TGLF uses physics-based saturation rules denoted SAT0, SAT1, and related variants, while the ST40 and surrogate studies explicitly discuss SAT2. This saturation dependence is not a secondary detail: it is one of the principal knobs governing transport stiffness and channel balance in applications.

## 4. Coupling to integrated transport solvers

TGLF is routinely embedded in transport frameworks rather than used in isolation. In the ST40 predictive workflow, ASTRA is coupled to SPIDER, NUBEAM, NCLASS, and a reduced scrape-off-layer model for last-closed-flux-surface boundary conditions. TGLF supplies the local anomalous fluxes, NCLASS supplies neoclassical coefficients and pinch velocities, NUBEAM provides particle, heat, and torque source terms, and SPIDER fixes the magnetic boundary via EFIT [2502.04993].

The transport equations advanced by ASTRA are written in flux coordinates \(\psi\) as
\[
\frac{\partial n_s}{\partial t}
+\frac1{V'}\frac{\partial}{\partial\psi}\Bigl[V'\,\Gamma_s(\psi)\Bigr]
= S_{n,s}(\psi),
\]
\[
\frac{\partial (n_s T_s)}{\partial t}
+\frac1{V'}\frac{\partial}{\partial\psi}\Bigl[V'\,Q_s(\psi)\Bigr]
= S_{T,s}(\psi),
\]
\[
\frac{\partial (n_i m_iR^2\Omega)}{\partial t}
+\frac1{V'}\frac{\partial}{\partial\psi}\Bigl[V'\,\Pi_\phi(\psi)\Bigr]
= T_\phi(\psi),
\]
with total fluxes split into neoclassical and TGLF contributions. In the same workflow, the reduced SOL “two-point” model sets LCFS boundary conditions through relations for \(n_e(1)\) and \(T_e(1)\) in terms of the outgoing core fluxes [2502.04993].

Quantitatively, the ST40 study reports typical mid-radius TGLF-predicted fluxes at \(\psi_N\approx0.5\) of \(Q_i \simeq 4.3\,\mathrm{GB}\), \(Q_e \simeq 3.8\,\mathrm{GB}\), and \(\Gamma_i \simeq 1.1\times 10^{20}\,\mathrm{s}^{-1}\). At the core, \(\psi_N<0.3\), strong \(\gamma_{E\times B}\) shear \((\sim0.08\,\omega_s)\) and fast-ion dilution yield near-neoclassical transport levels of \(\sim0.05\text{--}0.1\,\mathrm{GB}\), whereas toward the edge, \(\psi_N>0.8\), the anomalous \(\chi_i\) and \(\chi_e\) rise to a few \(\mathrm{m}^2/\mathrm{s}\) \((\sim1\text{--}3\,\mathrm{GB})\). The same study reports agreement in global quantities and kinetic profiles between predictive and interpretative modeling as well as experimental measurements, and summarizes the fully predictive coupled simulations as achieving global parameters and profiles within \(\sim20\%\) of the reference [2502.04993].

The role of TGLF in these workflows is therefore not limited to linear microstability assessment. It is a local transport closure used inside time-evolving, source-driven, geometry-constrained simulations.

## 5. Electromagnetic treatment, validation against higher-fidelity models, and reported shortcomings

A recurring issue in TGLF applications is the treatment of electromagnetic physics. The NSTX study distinguishes electrostatic (ES) and electromagnetic (EM) formulations: ES TGLF retains only \(\delta\phi\), whereas EM TGLF also includes parallel vector potential \(\delta A_\|\) via Ampère’s law and perpendicular magnetic fluctuations \(\delta B_\perp\) in the curvature drift. The paper states that at NSTX’s high \(\beta\), the ES model grossly overpredicts \(\chi_e\) and \(\chi_i\) by \(\approx2\times\), while the EM model saturates to a finite \(\beta\)-stabilized level [2509.04359].

The device-specific assessments reported in the ST40 and NSTX studies are summarized below.

| Regime | Reported TGLF behavior | Reported issue |
|---|---|---|
| ST40, \(\psi_N\ge 0.5\) | reproduces GS2 linear growth rates to \(\lesssim 20\%\) and nonlinear GS2 fluxes to \(\lesssim 30\%\) | satisfactory for TEM/UM/ITG-dominated transport |
| ST40, \(\psi_N\le 0.3\) | misses the KBM branch seen by GS2 | \(\delta B_\|\) effects omitted |
| NSTX, ES TGLF | \(T_e\) overpredicted by \(\approx 93\%\pm27\%\); \(T_i\) overpredicted by \(47\%\pm12\%\) | high-\(\beta\) electromagnetic physics absent |
| NSTX, EM TGLF | \(T_e\) overpredicted by \(46\%\pm30\%\); \(T_i\) underpredicted by \(25\%\pm5\%\) | persistent channel imbalance and saturation sensitivity |

In the ST40 analysis, TGLF is described as reproducing GS2 within \(\lesssim20\%\) at \(\psi_N\ge0.5\), but discrepancies grow in the core where fast-ion \(\beta\sim1\text{--}2\%\) and electromagnetic stabilization matter. The same analysis reports that TGLF is somewhat less stiff in its sensitivity to \(\gamma_{E\times B}\) shear and beam-ion dilution than GS2, and that mid-edge transport can be overpredicted by up to \(50\%\) unless tuned [2502.04993].

In the NSTX database of 37 high-performance discharges, the time-dependent TRANSP study finds that MMM more consistently agrees with the observations than TGLF, despite TGLF requiring orders of magnitude greater computational cost. EM-TGLF predicts \(T_e/T_i\approx1.8\), compared with experimental \(\approx0.9\), and the paper concludes that at NSTX’s high \(\beta\) the electrostatic TGLF model is unusable and electromagnetic fluctuations must be included. Even then, the study emphasizes the need for careful calibration of the spectral-shift model, specifically SAT1geo and SAT2, to avoid systematic over- or under-prediction [2509.04359].

These results delimit the present operating envelope of TGLF. The model is effective in many electrostatic or moderately electromagnetic regimes, but finite-\(\beta\), KBM-sensitive, MTM-sensitive, and fast-ion-dominated cases remain strongly dependent on closure choices and saturation calibration.

## 6. Surrogate models, differentiability, and accelerated TGLF workflows

Because native TGLF remains costly in integrated simulations, recent work has focused on neural surrogates. The TGLF-SINN study replaces the spectral summation with a fully differentiable neural-network surrogate and identifies three key innovations:

- **Feature engineering**: an element-wise asinh transform followed by standardization applied to both the total fluxes and each per-wavenumber contribution,
  \[
  \tilde{\mathbf F} = \mathcal S\bigl(\sinh^{-1}(\mathbf F)\bigr), \qquad
  \tilde{\mathbf F}^j = \mathcal S\bigl(\sinh^{-1}(\mathbf F^j)\bigr),
  \]
  which reduces target prediction range and permits negative predictions without clipping.
- **Spectra-informed architecture**: \(N_y=24\) parallel branches, one per \(k_y^j\), each taking the 31 physical features plus the scalar \(k_y^j\), processed by a shared Encoder–ResNet–Decoder to predict
  \[
  f_\theta(x,k_y^j)\approx \tilde{\mathbf F}^j,
  \]
  whose sum is inverted to recover the total flux.
- **Physics-guided spectral regularization and Bayesian active learning**: a combined loss
  \[
  \mathcal L = \mathcal L_f + \mathcal L_s
  \]
  with
  \[
  \mathcal L_f =\bigl\|\sum_{j=1}^{N_y}f_\theta(x,k_y^j)-\tilde{\mathbf F}\bigr\|^2,\qquad
  \mathcal L_s =\sum_{j=1}^{N_y}\bigl\|f_\theta(x,k_y^j)-\tilde{\mathbf F}^j\bigr\|^2,
  \]
  together with Expected Information Gain acquisition
  \[
  g_{\rm EIG}(x) =H\bigl[p(\hat{\mathbf F}\!\mid\!x)\bigr]
  - \mathbb E_{\theta^*}\!\bigl[H\bigl[p(\hat{\mathbf F}\!\mid\!x,\theta^*)\bigr]\bigr],
  \]
  estimated via Monte Carlo dropout and ensembling [2509.07024].

On the Major + Minor perturbation set of \(3\,\mathrm{M}\) outlier-filtered points, TGLF-SINN achieves an average LRMSE of \(0.0583\) versus \(0.0666\) for the SOTA TGLF-NN, a \(12.5\%\) relative reduction. In the BAL setting, using only \(25\%\) of the data raises LRMSE by only \(0.0165\) above the full-data TGLF-NN and by \(0.0248\) above the full-data TGLF-SINN. In downstream flux matching within the FUSE suite for DIII-D L-mode and H-mode, replacing native TGLF with TGLF-SINN reduces end-to-end runtime from 15 min on 128 threads to 20 s, a \(45\times\) speedup, while TGLF-NN in Julia takes 1.5 s, or \(600\times\) faster than full TGLF. Both surrogates yield experimentally validated profiles with sub-\(10\%\) relative errors [2509.07024].

The NSTX study reports a parallel surrogate effort in which TGLF-NN models are trained on \(O(10^6)\) TGLF runs from NSTX plasmas over \(\rho=0.1\text{--}0.9\), reproducing TGLF transport coefficients with \(\lesssim10\%\) local flux error. Embedded in the FUSE flux-matching solver rather than TRANSP/PT, TGLF-NN is reported as \(\approx600\times\) faster than full TGLF per slice, and a transfer-learned GKNN surrogate to QLGYRO is reported as \(\approx10^6\times\) faster than QLGYRO, with \(T_e\) RMSE \(\approx13\%\) and \(T_i\approx9\%\) in a 700-case validation [2509.04359].

A defining feature of TGLF-SINN is full differentiability with respect to the 31 transport inputs. The study explicitly connects this to gradient-based coupling to neoclassical, equilibrium, and control-system modules in whole-device simulation, adjoint-based optimization of magnetic configurations or heating profiles, and model-predictive schemes for real-time control of plasma turbulence. In that sense, the recent surrogate literature does not replace TGLF’s transport role; it re-expresses the TGLF mapping in a form that is faster, differentiable, and more data-efficient for integrated fusion workflows.

Source: https://www.emergentmind.com/topics/trapped-gyro-landau-fluid-tglf