Merge pull request #207 from mabarnes/pre-electron-updates

Preparation for adding electrons
mabarnes · Apr 29, 2024 · c0880a5 · c0880a5
2 parents 1800f80 + 7d3a3f7
commit c0880a5
Show file tree

Hide file tree

Showing 105 changed files with 5,642 additions and 4,178 deletions.
diff --git a/.gitignore b/.gitignore
@@ -26,3 +26,4 @@ makie_postproc.so
 plots_postproc.so
 precompile-temp
 post_processing_input.toml
+.DS_Store
diff --git a/docs/src/external_sources_notes.md b/docs/src/external_sources_notes.md
@@ -14,8 +14,8 @@ S_n &= A_n(r,z) \frac{1}{(\pi)^{3/2} (2 T_{\mathrm{source},n} / m_n)^{3/2}} \exp
 or in 1V simulations that do not include $v_\perp$, $v_\zeta$, $v_r$ dimensions
 ```math
 \begin{align}
-S_i &= A_i(r,z) \frac{1}{sqrt{\pi} \sqrt{2 T_{\mathrm{source},i} / m_i}} \exp\left( -\frac{v_\perp^2}{T_{\mathrm{source},i}} \right) \\
-S_n &= A_n(r,z) \frac{1}{sqrt{\pi} \sqrt{2 T_{\mathrm{source},n} / m_n}} \exp\left( -\frac{v_z^2}{T_{\mathrm{source},n}} \right)
+S_i &= A_i(r,z) \frac{1}{\sqrt{\pi} \sqrt{2 T_{\mathrm{source},i} / m_i}} \exp\left( -\frac{v_\perp^2}{T_{\mathrm{source},i}} \right) \\
+S_n &= A_n(r,z) \frac{1}{\sqrt{\pi} \sqrt{2 T_{\mathrm{source},n} / m_n}} \exp\left( -\frac{v_z^2}{T_{\mathrm{source},n}} \right)
 \end{align}
 ```
 

diff --git a/docs/src/moment_kinetic_equations.md b/docs/src/moment_kinetic_equations.md
@@ -1018,7 +1018,7 @@ and for neutrals where several of the ionization terms cancel
   & \qquad-\frac{w_{\|,n}}{2}\frac{1}{p_{\|,n}}\left(-\frac{\partial q_{\|,n}}{\partial z}
   - R_{in}\left(n_{i}p_{\|,n} - n_{n}p_{\|,i}
   - n_{n}n_{i}\left(u_{n} - u_{i}\right)^{2}\right)
-  + \int dv_\parallel S_{n} + u_{n}^2\int dv_\parallel v_\parallel^2 S_{n}
+  + \int dv_\parallel v_\parallel^2 S_{n} + u_{n}^2\int dv_\parallel S_{n}
   + v_{\mathrm{th},n}w_{\|,n}\frac{\partial p_{\|,n}}{\partial z}\right) \\
   & \qquad\left. + \frac{w_{\parallel,n}}{2}\frac{1}{n_{n}}\int dv_\parallel S_{n}
   + \frac{w_{\|,n}^{2}}{2}\frac{v_{\mathrm{th},n}}{n_{n}}\frac{\partial n_{n}}{\partial z}\right]\frac{\partial g_{n}}{\partial w_{\|,n}} \\

diff --git a/examples/fokker-planck-1D2V/fokker-planck-1D2V-even_nz-shorttest-nstep200.toml b/examples/fokker-planck-1D2V/fokker-planck-1D2V-even_nz-shorttest-nstep200.toml
@@ -86,7 +86,7 @@ vperp_discretization = "gausslegendre_pseudospectral"
 #vzeta_bc = "periodic"
 #vzeta_discretization = "chebyshev_pseudospectral"
 
-#[numerical_dissipation]
+#[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 0.0
 #vperp_dissipation_coefficient = 0.0
 #z_dissipation_coefficient = 0.1

diff --git a/examples/geometry/1D-mirror.toml b/examples/geometry/1D-mirror.toml
@@ -71,7 +71,7 @@ vz_L = 18.0
 vz_bc = "both_zero"
 vz_discretization = "chebyshev_pseudospectral"
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 vpa_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
 vperp_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
 force_minimum_pdf_value = 0.0

diff --git a/examples/gk-ions/2D-periodic-gk.toml b/examples/gk-ions/2D-periodic-gk.toml
@@ -73,15 +73,19 @@ vz_L = 18.0
 vz_bc = "both_zero"
 vz_discretization = "chebyshev_pseudospectral"
 
-[numerical_dissipation]
-vpa_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
-vperp_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
-#r_disspipation_coefficient = 1.0e-3
-#force_minimum_pdf_value = 0.0
-
 [geometry]
 #option="1D-mirror"
 DeltaB=0.0
 option="constant-helical"
 pitch=0.1
 rhostar= 0.1
+
+[ion_numerical_dissipation]
+vpa_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
+vperp_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
+#r_disspipation_coefficient = 1.0e-3
+#force_minimum_pdf_value = 0.0
+
+[neutral_numerical_dissipation]
+moment_dissipation_coefficient = 0.0001
+vz_dissipation_coefficient = 0.01
diff --git a/examples/numerical-dissipation/num-diss-relaxation.toml b/examples/numerical-dissipation/num-diss-relaxation.toml
@@ -58,7 +58,7 @@ nwrite = 2000
 nwrite_dfns = 2000
 split_operators = false
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 vpa_dissipation_coefficient = 0.1
 vperp_dissipation_coefficient = 0.1
 z_dissipation_coefficient = -1.0

diff --git a/examples/recycling-fraction/wall-bc_recyclefraction0.5_split3-init.toml b/examples/recycling-fraction/wall-bc_recyclefraction0.5_split3-init.toml
@@ -81,8 +81,14 @@ z_width = 0.125
 source_strength = 2.0
 source_T = 2.0
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3
 force_minimum_pdf_value = 0.0
+
+[neutral_numerical_dissipation]
+vz_dissipation_coefficient = 1.0e-1
+#vz_dissipation_coefficient = 1.0e-2
+#vz_dissipation_coefficient = 1.0e-3
+force_minimum_pdf_value = 0.0
diff --git a/examples/recycling-fraction/wall-bc_recyclefraction0.5_split3.toml b/examples/recycling-fraction/wall-bc_recyclefraction0.5_split3.toml
@@ -81,8 +81,14 @@ z_width = 0.125
 source_strength = 2.0
 source_T = 2.0
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3
 force_minimum_pdf_value = 0.0
+
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-1
+#vz_dissipation_coefficient = 1.0e-2
+#vz_dissipation_coefficient = 1.0e-3
+force_minimum_pdf_value = 0.0
diff --git a/examples/recycling-fraction/wall-bc_recyclefraction0.5_split3_fekete104.toml b/examples/recycling-fraction/wall-bc_recyclefraction0.5_split3_fekete104.toml
@@ -87,8 +87,14 @@ z_width = 0.125
 source_strength = 2.0
 source_T = 2.0
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3
 force_minimum_pdf_value = 0.0
+
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-1
+#vz_dissipation_coefficient = 1.0e-2
+#vz_dissipation_coefficient = 1.0e-3
+force_minimum_pdf_value = 0.0
diff --git a/examples/recycling-fraction/wall-bc_recyclefraction0.5_split3_fekete43.toml b/examples/recycling-fraction/wall-bc_recyclefraction0.5_split3_fekete43.toml
@@ -87,8 +87,14 @@ z_width = 0.125
 source_strength = 2.0
 source_T = 2.0
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3
 force_minimum_pdf_value = 0.0
+
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-1
+#vz_dissipation_coefficient = 1.0e-2
+#vz_dissipation_coefficient = 1.0e-3
+force_minimum_pdf_value = 0.0
diff --git a/examples/recycling-fraction/wall-bc_recyclefraction0.5_split3_fekete64.toml b/examples/recycling-fraction/wall-bc_recyclefraction0.5_split3_fekete64.toml
@@ -88,7 +88,7 @@ source_strength = 2.0
 source_T = 2.0
 
 [numerical_dissipation]
-#vpa_dissipation_coefficient = 1.0e-1
-#vpa_dissipation_coefficient = 1.0e-2
-#vpa_dissipation_coefficient = 1.0e-3
+#vz_dissipation_coefficient = 1.0e-1
+#vz_dissipation_coefficient = 1.0e-2
+#vz_dissipation_coefficient = 1.0e-3
 force_minimum_pdf_value = 0.0
diff --git a/examples/wall-bc/wall-bc_cheb.toml b/examples/wall-bc/wall-bc_cheb.toml
@@ -65,6 +65,10 @@ dt = 1.0e-5
 nwrite = 1000
 split_operators = false
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
 force_minimum_pdf_value = 0.0
+
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
+force_minimum_pdf_value = 0.0
diff --git a/examples/wall-bc/wall-bc_cheb_split1.toml b/examples/wall-bc/wall-bc_cheb_split1.toml
@@ -66,6 +66,10 @@ dt = 1.0e-5
 nwrite = 1000
 split_operators = false
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
 force_minimum_pdf_value = 0.0
+
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
+force_minimum_pdf_value = 0.0
diff --git a/examples/wall-bc/wall-bc_cheb_split2.toml b/examples/wall-bc/wall-bc_cheb_split2.toml
@@ -66,6 +66,10 @@ dt = 1.0e-5
 nwrite = 1000
 split_operators = false
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
 force_minimum_pdf_value = 0.0
+
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
+force_minimum_pdf_value = 0.0
diff --git a/examples/wall-bc/wall-bc_cheb_split3.toml b/examples/wall-bc/wall-bc_cheb_split3.toml
@@ -66,6 +66,10 @@ dt = 1.0e-5
 nwrite = 1000
 split_operators = false
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
 force_minimum_pdf_value = 0.0
+
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-3 #1.0e-2 #1.0e-1
+force_minimum_pdf_value = 0.0
diff --git a/examples/wall-bc/wall-bc_no-neutrals.toml b/examples/wall-bc/wall-bc_no-neutrals.toml
@@ -77,7 +77,7 @@ z_width = 0.125
 source_strength = 8.0
 source_T = 1.0
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3

diff --git a/examples/wall-bc/wall-bc_no-neutrals_split1.toml b/examples/wall-bc/wall-bc_no-neutrals_split1.toml
@@ -77,7 +77,7 @@ z_width = 0.125
 source_strength = 8.0
 source_T = 1.0
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3

diff --git a/examples/wall-bc/wall-bc_no-neutrals_split2.toml b/examples/wall-bc/wall-bc_no-neutrals_split2.toml
@@ -77,7 +77,7 @@ z_width = 0.125
 source_strength = 8.0
 source_T = 1.0
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3

diff --git a/examples/wall-bc/wall-bc_no-neutrals_split3.toml b/examples/wall-bc/wall-bc_no-neutrals_split3.toml
@@ -77,7 +77,7 @@ z_width = 0.125
 source_strength = 8.0
 source_T = 1.0
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3

diff --git a/examples/wall-bc/wall-bc_volumerecycle.toml b/examples/wall-bc/wall-bc_volumerecycle.toml
@@ -84,12 +84,18 @@ active = true
 source_type = "recycling"
 recycling_controller_fraction = 0.25
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3
 #force_minimum_pdf_value = 0.0
 
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-1
+#vz_dissipation_coefficient = 1.0e-2
+#vz_dissipation_coefficient = 1.0e-3
+force_minimum_pdf_value = 0.0
+
 [krook_collisions]
 use_krook = true
 frequency_option = "reference_parameters"
diff --git a/examples/wall-bc/wall-bc_volumerecycle_split1.toml b/examples/wall-bc/wall-bc_volumerecycle_split1.toml
@@ -84,12 +84,18 @@ active = true
 source_type = "recycling"
 recycling_controller_fraction = 0.25
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3
 #force_minimum_pdf_value = 0.0
 
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-1
+#vz_dissipation_coefficient = 1.0e-2
+#vz_dissipation_coefficient = 1.0e-3
+#force_minimum_pdf_value = 0.0
+
 [krook_collisions]
 use_krook = true
 frequency_option = "reference_parameters"
diff --git a/examples/wall-bc/wall-bc_volumerecycle_split2.toml b/examples/wall-bc/wall-bc_volumerecycle_split2.toml
@@ -84,12 +84,18 @@ active = true
 source_type = "recycling"
 recycling_controller_fraction = 0.25
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3
 #force_minimum_pdf_value = 0.0
 
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-1
+#vz_dissipation_coefficient = 1.0e-2
+#vz_dissipation_coefficient = 1.0e-3
+#force_minimum_pdf_value = 0.0
+
 [krook_collisions]
 use_krook = true
 frequency_option = "reference_parameters"
diff --git a/examples/wall-bc/wall-bc_volumerecycle_split3.toml b/examples/wall-bc/wall-bc_volumerecycle_split3.toml
@@ -84,12 +84,18 @@ active = true
 source_type = "recycling"
 recycling_controller_fraction = 0.25
 
-[numerical_dissipation]
+[ion_numerical_dissipation]
 #vpa_dissipation_coefficient = 1.0e-1
 #vpa_dissipation_coefficient = 1.0e-2
 #vpa_dissipation_coefficient = 1.0e-3
 force_minimum_pdf_value = 0.0
 
+[neutral_numerical_dissipation]
+#vz_dissipation_coefficient = 1.0e-1
+#vz_dissipation_coefficient = 1.0e-2
+#vz_dissipation_coefficient = 1.0e-3
+force_minimum_pdf_value = 0.0
+
 [krook_collisions]
 use_krook = true
 frequency_option = "reference_parameters"
diff --git a/machines/marconi/julia.env b/machines/marconi/julia.env
@@ -2,7 +2,7 @@
 
 module purge
 module load env-skl profile/base profile/advanced
-module load gnu/7.3.0 openmpi/3.1.4--gnu--7.3.0 intel/pe-xe-2018--binary python/3.9.4
+module load gnu/7.3.0 openmpi/3.1.4--gnu--7.3.0 intel/pe-xe-2018--binary python/3.9.4 git/2.17
 # Needed because Julia's Cairo library complains about libz<1.2.9
 module load zlib/1.2.11--intel--pe-xe-2018--binary