Excising the polar fluxes has proven unstable, disable by default

AFD-Illinois · Oct 4, 2024 · b8deac0 · b8deac0
1 parent 7a0d7c8
commit b8deac0
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 5 deletions.
diff --git a/pars/bondi/bondi_b_vertical.par b/pars/bondi/bondi_b_vertical.par
@@ -55,7 +55,9 @@ on = true
 # Transmitting boundaries, since there will be material at the pole
 inner_x2 = transmitting
 outer_x2 = transmitting
-excise_polar_flux = true
+# This might additionally help, but the option is unstable
+#excise_polar_flux = true
+
 # Allow in material from boundary
 check_inflow_outer_x1 = false
 

diff --git a/pars/tori_3d/mad.par b/pars/tori_3d/mad.par
@@ -66,8 +66,9 @@ two_sync = false
 <boundaries>
 inner_x2 = transmitting
 outer_x2 = transmitting
-# This reduces the polar "wake," but may shrink timestep
-excise_polar_flux = true
+# This reduces the polar "wake," but shrinks the timestep,
+# and has caused instabilities. YMMV.
+#excise_polar_flux = true
 
 <floors>
 frame        = drift

diff --git a/pars/tori_3d/sane.par b/pars/tori_3d/sane.par
@@ -64,8 +64,9 @@ two_sync = false
 <boundaries>
 inner_x2 = transmitting
 outer_x2 = transmitting
-# This reduces the polar "wake," but may shrink timestep
-excise_polar_flux = true
+# This reduces the polar "wake," but shrinks the timestep,
+# and has caused instabilities. YMMV.
+#excise_polar_flux = true
 
 <floors>
 # Now with FOFC it might work to set "frame = normal",