|
|
The QLKNN outputs fluxes in
|
|
|
[QuaLiKiz GB units](https://gitlab.com/qualikiz-group/QuaLiKiz/wikis/Input-and-output-variables#gyrobohm).
|
|
|
The inputs are QuaLiKiz convenient units.
|
|
|
|
|
|
| variable | definition | calculation | TCI field |
|
|
|
| -------- | ---------- | ----------- | --------- |
|
|
|
$`A_{ns}`$ | Normalized logarithmic density gradient | $`-\frac{R_0}{n_{s, FSA}} * \frac{dn_{s, FSA}}{dr_{minor}}`$ | - |
|
|
|
$`A_{ts}`$ | Normalized logarithmic temperature gradient | $`-\frac{R_0}{T_s} \frac{dT_s}{dr_{minor}}`$ | - |
|
|
|
$`n_{s,LFS}`$ | Density in [m**-3] on LFS midplane | - | `ni_LFS`/`ne_LFS` |
|
|
|
$`n_{s,FSA}`$ | Density in [m**-3] Flux Surface Avaraged | - | `ni`/`ne` |
|
|
|
$`T_e`$ | temperature in [J] | - | `Te` |
|
|
|
$`q_s`$ | Heat transport SI| $`(n_{s,LFS} * T_s * \chi_{GB})/a_{ref} * q_{s,NN}`$ | `heat_fluxes` |
|
|
|
$`\Gamma_e`$ | Particle transport SI| $`(n_{s,LFS} * \chi_{GB})/a_{ref} * \Gamma_{e,NN}`$ | `particle_fluxes` |
|
|
|
$`\chi_{s, eff}`$ | Effective heat diffusivity | $`q_s / (n_{s,FSA} * A_{ts} * \frac{T_s}{R_0})`$ | `chieff` |
|
|
|
$`V_{s, t, eff}`$ | Effective heat pinch | $`q_s / (n_{s,FSA} T_s)`$ | - |
|
|
|
$`D_{s, eff}`$ | Effective diffusivity | $`\Gamma_{FSA} / (A_{ns} \frac{n_{s, FSA}}{R_0})`$ | `deff` |
|
|
|
$`V_{s, p, eff}`$ | Effective particle pinch | $`\Gamma_{FSA} / n_{s, FSA}`$ | - |
|
|
|
|
|
|
Note that QuaLiKiz $`R_0 == R_{0,TCI}`$, and we don't have heavy impurity 2D
|
|
|
effects, so `$\Gamma_{FSA} == \Gamma_{asym} == \Gamma$`. `$aref=Rmin_in(1)$`.
|
|
|
|
|
|
Now there are multiple options to get the $`\chi`$, $`V_t`$, $`D`$ and $`V`$
|
|
|
terms for the transport equation solver (TCI `chi`, `Vt`, `d` and `Vp` respectively)
|
|
|
|
|
|
## 1. effective flux
|
|
|
If `heat_fluxes` * $`A_{ts}`$ > 0 and $`|A_{ts}|`$ >=0.1:
|
|
|
|
|
|
`chi` = $`\chi_{s, eff}`$
|
|
|
|
|
|
`Vt` = 0
|
|
|
|
|
|
else (as RAPTOR doesn't have `Vt`, this is never true in RAPTOR)
|
|
|
|
|
|
`chi` = 0
|
|
|
|
|
|
`Vt` = $`V_{s, t, eff}`$
|
|
|
|
|
|
Similar for the particle fluxes, if `particle_fluxes` * $`A_{ns}`$ > 0 and
|
|
|
$`|A_{ns}|`$ >=0.1:
|
|
|
|
|
|
|
|
|
`D` = $`D_{s, eff}`$
|
|
|
|
|
|
`Vp` = 0
|
|
|
|
|
|
else (as RAPTOR doesn't have `Vt`, this is never true in RAPTOR)
|
|
|
|
|
|
`D` = 0
|
|
|
|
|
|
`Vp` = $`V_{s, p, eff}`$
|
|
|
|
|
|
2. Using the electron D's/V's
|
|
|
|
|
|
# Transport in transport codes
|
|
|
In general it depends on the transport equation being solved which transport
|
|
|
coefficients QLKNN has to provide. In jetto this is for _all_ ions and electrons:
|
|
|
|
|
|
``` fortran
|
|
|
TRANSPORT%chi
|
|
|
TRANSPORT%Vt
|
|
|
TRANSPORT%d
|
|
|
TRANSPORT%Vp
|
|
|
TRANSPORT%chimom
|
|
|
``` |