evolve_momentum: Change fix_boundary_momentum_flux to false

bendudson · bendudson · commit 7bff6997a6d4 · 2023-03-07T14:48:28.000-08:00
This change introduces a parallel viscosity that results in
smoother profiles into the sheath, and better energy conservation
in transport simulations.
diff --git a/src/evolve_momentum.cxx b/src/evolve_momentum.cxx
@@ -46,7 +46,7 @@ EvolveMomentum::EvolveMomentum(std::string name, Options &alloptions, Solver *so
 
   fix_momentum_boundary_flux = options["fix_momentum_boundary_flux"]
     .doc("Fix Y boundary momentum flux to boundary midpoint value?")
-    .withDefault<bool>(true);
+    .withDefault<bool>(false);
 }
 
 void EvolveMomentum::transform(Options &state) {

Original file line number	Diff line number	Diff line change
`@@ -46,7 +46,7 @@ EvolveMomentum::EvolveMomentum(std::string name, Options &alloptions, Solver *so`
`46`	`46`
`47`	`47`	`fix_momentum_boundary_flux = options["fix_momentum_boundary_flux"]`
`48`	`48`	`.doc("Fix Y boundary momentum flux to boundary midpoint value?")`
`49`		`- .withDefault<bool>(true);`
	`49`	`+ .withDefault<bool>(false);`
`50`	`50`	`}`
`51`	`51`
`52`	`52`	`void EvolveMomentum::transform(Options &state) {`