From 0b3c9a44bcded961f0d1b672ed240f34d11a51b9 Mon Sep 17 00:00:00 2001 From: Antoine Hoffmann <antoine.hoffmann@epfl.ch> Date: Wed, 20 Sep 2023 09:55:34 +0200 Subject: [PATCH] new test case for ExB shearing --- testcases/CBC_ExBshear/fort_00.90 | 97 +++++++++++++++++++ testcases/ExB_shear_atomic_tests/fort_00.90 | 101 ++++++++++++++++++++ testcases/ExB_shear_atomic_tests/fort_01.90 | 98 +++++++++++++++++++ testcases/ExB_shear_atomic_tests/fort_02.90 | 101 ++++++++++++++++++++ testcases/ExB_shear_atomic_tests/fort_03.90 | 101 ++++++++++++++++++++ testcases/ExB_shear_atomic_tests/fort_04.90 | 101 ++++++++++++++++++++ testcases/ExB_shear_atomic_tests/fort_05.90 | 100 +++++++++++++++++++ 7 files changed, 699 insertions(+) create mode 100644 testcases/CBC_ExBshear/fort_00.90 create mode 100644 testcases/ExB_shear_atomic_tests/fort_00.90 create mode 100644 testcases/ExB_shear_atomic_tests/fort_01.90 create mode 100644 testcases/ExB_shear_atomic_tests/fort_02.90 create mode 100644 testcases/ExB_shear_atomic_tests/fort_03.90 create mode 100644 testcases/ExB_shear_atomic_tests/fort_04.90 create mode 100644 testcases/ExB_shear_atomic_tests/fort_05.90 diff --git a/testcases/CBC_ExBshear/fort_00.90 b/testcases/CBC_ExBshear/fort_00.90 new file mode 100644 index 00000000..a65bd044 --- /dev/null +++ b/testcases/CBC_ExBshear/fort_00.90 @@ -0,0 +1,97 @@ +&BASIC + nrun = 99999999 + dt = 0.01 + tmax = 50 + maxruntime = 72000 + job2load = -1 +/ +&GRID + pmax = 2 + jmax = 1 + Nx = 64 + Lx = 120 + Ny = 48 + Ly = 120 + Nz = 16 + SG = .f. + Nexc = 0 +/ +&GEOMETRY + geom = 's-alpha' + !geom = 'miller' + q0 = 1.4 + shear = 0.8 + eps = 0.18 + kappa = 1.0 + s_kappa= 0.0 + delta = 0.0 + s_delta= 0.0 + zeta = 0.0 + s_zeta = 0.0 + parallel_bc = 'dirichlet' + shift_y= 0.0 +/ +&DIAGNOSTICS + dtsave_0d = 1 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 1 + dtsave_5d = 20 + write_doubleprecision = .f. + write_gamma = .t. + write_hf = .t. + write_phi = .t. + write_Na00 = .t. + write_Napj = .t. + write_dens = .t. + write_fvel = .t. + write_temp = .t. +/ +&MODEL + LINEARITY = 'nonlinear' + Na = 1 ! number of species + mu_x = 0.0 + mu_y = 0.0 + N_HD = 4 + mu_z = 1.0 + HYP_V = 'hypcoll' + mu_p = 0.0 + mu_j = 0.0 + nu = 0.05 + beta = 0.0 + ADIAB_E = .t. + ExBrate = 0.01 +/ +&CLOSURE + hierarchy_closure='truncation' + dmax = -1 + nonlinear_closure='truncation' + nmax = 0 +/ +&SPECIES + ! ions + name_ = 'ions' + tau_ = 1.0 + sigma_= 1.0 + q_ = 1.0 + k_N_ = 2.22 + k_T_ = 6.96 +/ + +&COLLISION + collision_model = 'DG' !DG/SG/PA/LD (dougherty, sugama, pitch angle, landau) + GK_CO = .f. + INTERSPECIES = .true. + mat_file = 'gk_sugama_P_20_J_10_N_150_kpm_8.0.h5' +/ +&INITIAL + INIT_OPT = 'blob' + ACT_ON_MODES = 'donothing' + init_background = 0.0 + init_noiselvl = 0.005 + iseed = 42 +/ +&TIME_INTEGRATION + !numerical_scheme = 'RK4' + numerical_scheme = 'SSP_RK3' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_00.90 b/testcases/ExB_shear_atomic_tests/fort_00.90 new file mode 100644 index 00000000..29ec2972 --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_00.90 @@ -0,0 +1,101 @@ +&BASIC + nrun = 100000000 + dt = 0.00075 + tmax = 2.0 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 186 + Lx = 100 + Ny = 186 + Ly = 100 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.01 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.01 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'nonlinear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 1 + mu_y = 1 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 1 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 0.0 + ikxZF = 2 + ZFamp = 400 + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =0 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'ricci' + init_amp = 0.1 +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_01.90 b/testcases/ExB_shear_atomic_tests/fort_01.90 new file mode 100644 index 00000000..2aa6249f --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_01.90 @@ -0,0 +1,98 @@ +&BASIC + nrun = 100000000 + dt = 0.00075 + tmax = 0.5 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 186 + Lx = 100 + Ny = 186 + Ly = 100 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.01 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.01 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'linear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 1 + mu_y = 1 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 1 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 5.0 + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =0 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'ricci' +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_02.90 b/testcases/ExB_shear_atomic_tests/fort_02.90 new file mode 100644 index 00000000..898315bb --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_02.90 @@ -0,0 +1,101 @@ +&BASIC + nrun = 100000000 + dt = 0.001 + tmax = 5.0 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 96 + Lx = 100 + Ny = 64 + Ly = 100 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.05 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.05 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'linear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 0 + mu_y = 0 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 0 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 3.0 + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =-1 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'mom00_single_mode' + init_amp = 100 + ikx_init = 3 + iky_init = 3 +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_03.90 b/testcases/ExB_shear_atomic_tests/fort_03.90 new file mode 100644 index 00000000..d1cbbdfa --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_03.90 @@ -0,0 +1,101 @@ +&BASIC + nrun = 100000000 + dt = 0.001 + tmax = 5.0 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 96 + Lx = 100 + Ny = 64 + Ly = 100 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.05 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.05 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'nonlinear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 0 + mu_y = 0 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 0 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 0.0 + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =-1 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'mom00_single_mode' + init_amp = 100 + ikx_init = 2 + iky_init = 2 +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_04.90 b/testcases/ExB_shear_atomic_tests/fort_04.90 new file mode 100644 index 00000000..05e8113b --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_04.90 @@ -0,0 +1,101 @@ +&BASIC + nrun = 100000000 + dt = 0.001 + tmax = 5.0 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 96 + Lx = 100 + Ny = 64 + Ly = 100 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.05 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.05 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'nonlinear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 0 + mu_y = 0 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 0 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 0.0 + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =-1 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'mom00_single_mode' + init_amp = 100 + ikx_init = 3 + iky_init = 3 +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ diff --git a/testcases/ExB_shear_atomic_tests/fort_05.90 b/testcases/ExB_shear_atomic_tests/fort_05.90 new file mode 100644 index 00000000..9fcce56d --- /dev/null +++ b/testcases/ExB_shear_atomic_tests/fort_05.90 @@ -0,0 +1,100 @@ +&BASIC + nrun = 100000000 + dt = 0.001 + tmax = 5.0 + maxruntime = 356400 + job2load = -1 +/ +&GRID + pmax = 0 + jmax = 0 + Nx = 96 + Lx = 120 + Ny = 64 + Ly = 90 + Nz = 1 + SG = .false. + Nexc = 1 +/ +&GEOMETRY + geom = 'zpinch'!'miller' + q0 = 1.0 + shear = 0.0 + eps = 0.0 + kappa = 1 + s_kappa = 0 + delta = 0 + s_delta = 0 + zeta = 0 + s_zeta = 0 + parallel_bc = 'dirichlet' + shift_y = 0 + Npol = 1 + PB_PHASE= .false. +/ +&DIAGNOSTICS + dtsave_0d = 0.05 + dtsave_1d = -1 + dtsave_2d = -1 + dtsave_3d = 0.05 + dtsave_5d = 100 + write_doubleprecision = .true. + write_gamma = .true. + write_hf = .true. + write_phi = .true. + write_Na00 = .true. + write_Napj = .true. + write_dens = .true. + write_temp = .true. +/ +&MODEL +LINEARITY = 'nonlinear' +RM_LD_T_EQ= .false. + Na = 1 + mu_x = 0 + mu_y = 0 + N_HD = 4 + mu_z = 0 + HYP_V = 'hypcoll' + mu_p = 0 + mu_j = 0 + nu = 0 + k_gB = 0 + k_cB = 0 + lambdaD = 0 + beta = 0.0 + ExBrate = 3.0 + ExB_NL_CORRECTION = .true. + ADIAB_E = .true. + ADIAB_I = .false. + tau_i = 1 + MHD_PD = .false. +/ +&CLOSURE + hierarchy_closure='truncation' + dmax =-1 + nonlinear_closure='truncation' + nmax =-1 +/ +&SPECIES + name_ = 'ions' + tau_ = 1 + sigma_ = 1 + q_ = 1 + K_N_ = 0 + K_T_ = 0 +/ +&COLLISION + collision_model = 'DG' + GK_CO = .false. + INTERSPECIES = .true. + mat_file = '/home/ahoffman/gyacomo/iCa/null' + collision_kcut = 1 +/ +&INITIAL + INIT_OPT = 'mom00_single_mode' + init_amp = 100 +/ +&TIME_INTEGRATION + numerical_scheme = 'RK4' +/ -- GitLab