Skip to content

Commit

Permalink
Bugfix: make electron conduction kappa Braginskii consistent
Browse files Browse the repository at this point in the history
  • Loading branch information
mikekryjak committed Jul 8, 2024
1 parent 8ad24b8 commit baf4606
Showing 1 changed file with 4 additions and 2 deletions.
6 changes: 4 additions & 2 deletions src/evolve_pressure.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -76,9 +76,11 @@ EvolvePressure::EvolvePressure(std::string name, Options& alloptions, Solver* so
.doc("Include parallel heat conduction?")
.withDefault<bool>(true);

// This is consistent with Bragkinskii definition of kappa and tau_e but with the sqrt(2)
// moved into the kappa so that the collision time is consistent with Fitzpatrick.
kappa_coefficient = options["kappa_coefficient"]
.doc("Numerical coefficient in parallel heat conduction. Default is 3.16 for electrons, 3.9 otherwise")
.withDefault((name == "e") ? 3.16 : 3.9);
.doc("Numerical coefficient in parallel heat conduction. Default is 3.16/sqrt(2) for electrons, 3.9 otherwise")
.withDefault((name == "e") ? 3.16/sqrt(2) : 3.9);

kappa_limit_alpha = options["kappa_limit_alpha"]
.doc("Flux limiter factor. < 0 means no limit. Typical is 0.2 for electrons, 1 for ions.")
Expand Down

0 comments on commit baf4606

Please sign in to comment.