From f9688909c70cc16da71ee867f7638fae6b338fd9 Mon Sep 17 00:00:00 2001 From: d-burg Date: Wed, 28 Oct 2026 11:58:39 -0400 Subject: [PATCH] InnerLayer.SLAYER - DOCS - Describe chi_tor as the toroidal momentum diffusivity, not a heat diffusivity Co-Authored-By: Claude Opus 5.5 --- src/InnerLayer/SLAYER/LayerInputs.jl | 5 +++-- src/InnerLayer/SLAYER/LayerParameters.jl | 2 +- 2 files changed, 4 insertions(+), 3 deletions(-) diff --git a/src/InnerLayer/SLAYER/LayerInputs.jl b/src/InnerLayer/SLAYER/LayerInputs.jl index ed9bc0aa13..f1e37610f2 100644 --- a/src/InnerLayer/SLAYER/LayerInputs.jl +++ b/src/InnerLayer/SLAYER/LayerInputs.jl @@ -195,8 +195,9 @@ profiles, without an intermediate file round-trip. `:flux` label, whatever `rs_method` is. - `P_perp_model` -- how `P_perp` is built; see [`slayer_parameters`](@ref) (default `:D_perp`, as in the `[SLAYER]` deck). - - `chi_tor` -- toroidal heat diffusivity [m²/s]. Scalar or a callable - of `psi` (default `1.0`). + - `chi_tor` -- toroidal momentum diffusivity χ_φ [m²/s], setting the momentum + Prandtl number `P_tor` = τ_R/τ_φ with τ_φ = r_s²/χ_φ (Fitzpatrick, Phys. Plasmas + 30, 092512 (2023)). Scalar or a callable of `psi` (default `1.0`). - `dr_val` -- resistive interchange index `D_R = E + F + H²` (Glasser-Greene-Johnson 1975) feeding the critical-Δ formulas diff --git a/src/InnerLayer/SLAYER/LayerParameters.jl b/src/InnerLayer/SLAYER/LayerParameters.jl index a8d34110bf..b990f5d22a 100644 --- a/src/InnerLayer/SLAYER/LayerParameters.jl +++ b/src/InnerLayer/SLAYER/LayerParameters.jl @@ -29,7 +29,7 @@ de-normalization. The parametrization uses `P_perp`, `P_tor`, and | `c_beta` | Compressibility √(β_local / (1 + β_local)) | | `D_norm` | (d_β/r_s) · S^(1/3) · √ι_e (Fitzpatrick normalized scale) | | `P_perp` | Perpendicular Prandtl number τ_R / τ_⊥ (set by `P_perp_model`) | -| `P_tor` | Toroidal-direction Prandtl number τ_R / τ_‖tor | +| `P_tor` | Momentum Prandtl number τ_R / τ_φ, τ_φ = r_s²/χ_φ | | `Q_e` | Normalized electron diamagnetic: −tauk · ω_*e | | `Q_i` | Normalized ion diamagnetic: −tauk · ω_*i | | `iota_e` | Q_e / (Q_e − Q_i) |