Skip to content

Add neutral equation term diagnostics #542

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
sun51027 wants to merge 12 commits into
base: neutrals-flux-limmiters
Choose a base branch
from
Open

Add neutral equation term diagnostics #542

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
12 commits
Select commit Hold shift + click to select a range
664bd23
add vth diagnose
sun51027 Apr 1, 2026
9669d56
add physcis terms
sun51027 Apr 1, 2026
161f8bf
add terms
sun51027 Apr 1, 2026
60bdbc6
add all diagnostic for each evolution term
sun51027 Apr 1, 2026
873c934
fix ddtNVn_par_advection += to =
sun51027 Apr 9, 2026
0210aa3
add viscosity term
sun51027 Apr 10, 2026
d5fc9f3
[bot] Apply format changes
sun51027 Apr 10, 2026
96fbef6
[bot] Add last format changes commit to ignore file
sun51027 Apr 10, 2026
75314e5
add viscosity and anomalous_D term
sun51027 Apr 10, 2026
5eeba36
initialise terms = 0.0 if boolean variable is required
sun51027 Apr 10, 2026
acf4625
[bot] Apply format changes
sun51027 Apr 20, 2026
27c0bf8
[bot] Add last format changes commit to ignore file
sun51027 Apr 20, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 2 additions & 0 deletions .git-blame-ignore-revs
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2,3 +2,5 @@ f5de620399e88f139e236e40675facfd0cfd6b13
07616298519738fad0c1bb259e4d021cd6f51d58
37efa12abe0961add8715b797d3926664977d453
8e49d5933131ddcb2f97bb2ac7095a1be80ad7b7
d5fc9f3f5f0a75b8de0f35aaafe541114a241951
acf4625dd0940720ccfdc38af61493e19cde21ed
26 changes: 26 additions & 0 deletions include/neutral_mixed.hxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,6 +45,9 @@ private:
std::vector<std::string> collision_names; ///< Collisions used for collisionality
std::string
diffusion_collisions_mode; ///< Collision selection, either afn or multispecies
Field3D Vnth_pf = 0.0;
Field3D Vnth_hf = 0.0;

Field3D nu; ///< Collisionality to use for diffusion
Field3D nu_pseudo_mfp; ///< Pseudo-collision frequency based on mean free path
Field3D nu_total; ///< Total collision frequency used for diffusion, including
Expand Down Expand Up @@ -105,6 +108,29 @@ private:
bool output_ddt; ///< Save time derivatives?
bool diagnose; ///< Save additional diagnostics?

// Physic terms
Field3D ddtN_par_advection = 0.0;
Field3D ddtN_perp_diffusion = 0.0;

Field3D ddtPn_par_advection = 0.0;
Field3D ddtPn_work_done = 0.0;
Field3D ddtPn_perp_advection = 0.0;
Field3D ddtPn_par_conduction = 0.0;
Field3D ddtPn_perp_conduction = 0.0;

Field3D ddtNVn_par_advection = 0.0;
Field3D ddtNVn_pressure_gradient = 0.0;
Field3D ddtNVn_perp_advection = 0.0;

Field3D par_viscosity_source = 0.0;
Field3D perp_viscosity_source = 0.0;
Field3D ddtNVn_viscosity = 0.0; // par_viscosity_source + perp_viscosity_source
Field3D ddtPn_viscosity = 0.0;

Field3D ddtN_anomalous_transport = 0.0;
Field3D ddtNVn_anomalous_transport = 0.0;
Field3D ddtPn_anomalous_transport = 0.0;

// Flow diagnostics
Field3D pf_adv_perp_xlow, pf_adv_perp_ylow, pf_adv_par_ylow;
Field3D mf_adv_perp_xlow, mf_adv_perp_ylow, mf_adv_par_ylow;
Expand Down
Loading