From ba4856c67b519f301cbae7048bc5d58ee82eea88 Mon Sep 17 00:00:00 2001 From: Jorge Gonzalez Date: Fri, 29 Nov 2024 10:01:00 +0100 Subject: [PATCH] New boundary files. Now they follow Z = 22.5 T^0.6 --- bc_50ns_T10Z6.csv | 5 +++++ bc_50ns_T30Z11.csv | 5 +++++ bc_50ns_T60Z16.csv | 5 +++++ bc_50ns_T6Z4.csv | 5 +++++ bc_5ns.csv | 5 ----- bc_fa_T30Z11.csv | 3 +++ bc_fullAblation.csv | 3 --- plasmaExpansion.f90 | 2 +- 8 files changed, 24 insertions(+), 9 deletions(-) create mode 100644 bc_50ns_T10Z6.csv create mode 100644 bc_50ns_T30Z11.csv create mode 100644 bc_50ns_T60Z16.csv create mode 100644 bc_50ns_T6Z4.csv delete mode 100644 bc_5ns.csv create mode 100644 bc_fa_T30Z11.csv delete mode 100644 bc_fullAblation.csv diff --git a/bc_50ns_T10Z6.csv b/bc_50ns_T10Z6.csv new file mode 100644 index 0000000..98fb235 --- /dev/null +++ b/bc_50ns_T10Z6.csv @@ -0,0 +1,5 @@ +t (s),n (m^-3),u (m s^-1),T (eV),Z +0.000000E-000,1.111111E+025,0.000000E+000,1.000000E+001,6.000000E+000 +5.000000E-008,1.111111E+025,0.000000E+000,1.000000E+001,6.000000E+000 +5.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,6.000000E+000 +2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,6.000000E+000 diff --git a/bc_50ns_T30Z11.csv b/bc_50ns_T30Z11.csv new file mode 100644 index 0000000..856d5d6 --- /dev/null +++ b/bc_50ns_T30Z11.csv @@ -0,0 +1,5 @@ +t (s),n (m^-3),u (m s^-1),T (eV),Z +0.000000E-000,1.111111E+025,0.000000E+000,3.000000E+001,1.100000E+001 +5.000000E-008,1.111111E+025,0.000000E+000,3.000000E+001,1.100000E+001 +5.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,1.100000E+001 +2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,1.100000E+001 diff --git a/bc_50ns_T60Z16.csv b/bc_50ns_T60Z16.csv new file mode 100644 index 0000000..a3eff24 --- /dev/null +++ b/bc_50ns_T60Z16.csv @@ -0,0 +1,5 @@ +t (s),n (m^-3),u (m s^-1),T (eV),Z +0.000000E-000,1.111111E+025,0.000000E+000,6.000000E+001,1.600000E+001 +5.000000E-008,1.111111E+025,0.000000E+000,6.000000E+001,1.600000E+001 +5.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,1.600000E+001 +2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,1.600000E+001 diff --git a/bc_50ns_T6Z4.csv b/bc_50ns_T6Z4.csv new file mode 100644 index 0000000..4d623c7 --- /dev/null +++ b/bc_50ns_T6Z4.csv @@ -0,0 +1,5 @@ +t (s),n (m^-3),u (m s^-1),T (eV),Z +0.000000E-000,1.111111E+025,0.000000E+000,6.000000E+000,4.000000E+000 +5.000000E-008,1.111111E+025,0.000000E+000,6.000000E+000,4.000000E+000 +5.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,4.000000E+000 +2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,4.000000E+000 diff --git a/bc_5ns.csv b/bc_5ns.csv deleted file mode 100644 index 306b9b5..0000000 --- a/bc_5ns.csv +++ /dev/null @@ -1,5 +0,0 @@ -t (s),n (m^-3),u (m s^-1),T (eV),Z -0.000000E-000,1.111111E+025,0.000000E+000,3.000000E+001,9.000000E+000 -5.000000E-008,1.111111E+025,0.000000E+000,3.000000E+001,9.000000E+000 -5.500000E-008,1.000000E+020,0.000000E+000,5.000000E+000,9.000000E+000 -2.000000E-006,1.000000E+020,0.000000E+000,5.000000E+000,9.000000E+000 diff --git a/bc_fa_T30Z11.csv b/bc_fa_T30Z11.csv new file mode 100644 index 0000000..2c0c946 --- /dev/null +++ b/bc_fa_T30Z11.csv @@ -0,0 +1,3 @@ +t (s),n (m^-3),u (m s^-1),T (eV),Z +0.000000E-000,1.111111E+025,0.000000E+000,3.000000E+001,1.100000E+001 +1.000000E-006,1.111111E+025,0.000000E+000,3.000000E+001,1.100000E+001 diff --git a/bc_fullAblation.csv b/bc_fullAblation.csv deleted file mode 100644 index 487aee3..0000000 --- a/bc_fullAblation.csv +++ /dev/null @@ -1,3 +0,0 @@ -t (s),n (m^-3),u (m s^-1),T (eV),Z -0.000000E-000,1.111111E+025,0.000000E+000,3.000000E+001,9.000000E+000 -1.000000E-006,1.111111E+025,0.000000E+000,3.000000E+001,9.000000E+000 diff --git a/plasmaExpansion.f90 b/plasmaExpansion.f90 index fd24795..f715510 100644 --- a/plasmaExpansion.f90 +++ b/plasmaExpansion.f90 @@ -93,7 +93,7 @@ program plasmaExpansion phi_ref = kb * Temp_ref / qe ! Set input parameters (remember these have to be in non-dimensional units) - c_s = sqrt(9.0_dp * gamma_i * 1.0_dp) + c_s = sqrt(11.0_dp * gamma_i * 1.0_dp) bc_file = 'bc.csv' call boundaryConditions%init(bc_file)