Repository navigation
InnerLayer.SLAYER - Add D_perp, tau_E and other options for P_perp #509
Description
Summary
The InnerLayer (IL) SLAYER model has a single way to build the perpendicular Prandtl number P_⊥, from the electron perpendicular heat diffusivity χ⊥,e. Add a [SLAYER] P_perp_model option so P_⊥ can be built from the closures in the layer-theory literature, and an optional ion perpendicular heat diffusivity chi_i profile in the kinetic HDF5 schema to feed them.
Implemented in #501.
Options
P_perp_model |
P_⊥ | Reference |
|---|---|---|
chi_perp_e |
τ_R χ⊥,e / r_s² (current behaviour) | — |
chi_perp_i |
τ_R χ⊥,i / r_s² | — |
D_perp |
τ_R D⊥ / r_s², with D⊥ from c_β, η, χ⊥,e, χ⊥,i and the n, T gradients | Fitzpatrick, Phys. Plasmas 29, 032507 (2022), Eq. 16 |
c_beta |
c_β², the χ⊥ → 0 limit of D_perp (a floor on P_⊥) |
J.-K. Park, Phys. Plasmas 29, 072506 (2022), Eq. 10, thermal-conduction term K dropped |
P_phi |
P_⊥ = P_φ | Fitzpatrick, Phys. Plasmas 30, 092512 (2023), Eq. 143 |
tau_E |
P_⊥ = P_φ = τ_R / τ_E, from a user-supplied tau_E [s] |
Fitzpatrick, Phys. Plasmas 30, 092512 (2023), Eq. 131 |
Interface
- New
[SLAYER] P_perp_modelandtau_Ekeys ingpec.toml. chi_perprenamed tochi_perp_ein[SLAYER]and in thebuild_slayer_inputs/slayer_parameterskeywords.- Optional
chi_idataset in the kinetic HDF5 file, read only by the models that need it.
The default model is tracked separately in #508.
Activity
logan-nc commentedon Oct 27, 2026
All six models, the
chi_perp_erename,tau_Eand the optionalchi_iprofile are in #501 at its current head, including d-burg's corrections to theD_perpdenominator and the local τ_E conversion. Close with #501.