Skip to content

Commit 4397928

Browse files
johnomotanijohnomotani
authored
Merge pull request #429 from mabarnes/2D1V-inputs
Inputs for 2D instability demonstrations
2 parents 710afcc + 92c49eb commit 4397928

109 files changed

Lines changed: 11286 additions & 0 deletions

File tree

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

‎publication_inputs/2D-instability/input_files/2D1V-instability-periodic-test/2D1V-instability-periodic-Lz0.02-test-Fourier-r-AT0.1-An-0.005.toml‎

Lines changed: 99 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,99 @@
1+
[composition]
2+
T_wall = 0.1
3+
T_e = 0.2
4+
n_ion_species = 1
5+
n_neutral_species = 0
6+
7+
[evolve_moments]
8+
density = false
9+
parallel_flow = false
10+
pressure = false
11+
12+
[krook_collisions]
13+
use_krook = false
14+
15+
[geometry]
16+
option = "constant-helical"
17+
pitch = 0.1
18+
19+
[manufactured_solns]
20+
use_for_init = true
21+
use_for_advance = true
22+
type = "2D-instability"
23+
24+
[ion_species_1]
25+
initial_temperature = 0.2
26+
initial_density = 1.0
27+
28+
[z_IC_ion_species_1]
29+
initialization_option = "sinusoid"
30+
density_amplitude = -0.005
31+
temperature_amplitude = 0.1
32+
density_phase = 0.0
33+
upar_amplitude = 0.0
34+
temperature_phase = 0.0
35+
upar_phase = 0.0
36+
37+
[vpa_IC_ion_species_1]
38+
initialization_option = "gaussian"
39+
density_amplitude = 1.0
40+
temperature_amplitude = 0.0
41+
density_phase = 0.0
42+
upar_amplitude = 0.0
43+
temperature_phase = 0.0
44+
upar_phase = 0.0
45+
46+
[r]
47+
ngrid = 4
48+
discretization = "fourier_pseudospectral"
49+
nelement = 1
50+
bc = "periodic"
51+
L = 1.0e-4 # equivalent to 1mm
52+
53+
[z]
54+
ngrid = 5
55+
discretization = "gausslegendre_pseudospectral"
56+
nelement = 32
57+
#nelement_local = 8
58+
bc = "periodic"
59+
L = 0.02 # 1m vertical height with pitch=0.1 corresponds to 10m of parallel connection length.
60+
61+
[vpa]
62+
ngrid = 6
63+
discretization = "gausslegendre_pseudospectral"
64+
nelement = 31
65+
L = 10.0
66+
element_spacing_option = "coarse_tails8"
67+
bc = "zero"
68+
69+
[timestepping]
70+
type = "Fekete4(3)"
71+
nstep = 3000000
72+
dt = 1.0e-5
73+
minimum_dt = 1.0e-5
74+
CFL_prefactor = 2.0
75+
rtol = 1.0e-5
76+
atol = 1.0e-12
77+
add_random_noise = 1.0e-9
78+
random_noise_seed = 1234
79+
nwrite = 10000
80+
nwrite_dfns = 10000
81+
steady_state_residual = true
82+
converged_residual_value = 1.0e-12
83+
84+
[ion_numerical_dissipation]
85+
r_dissipation_coefficient = 0.0e-10
86+
#vpa_dissipation_coefficient = 1.0e-1
87+
vpa_dissipation_coefficient = 1.0e-2
88+
#vpa_dissipation_coefficient = 1.0e-3
89+
#force_minimum_pdf_value = 0.0
90+
91+
[electron_numerical_dissipation]
92+
#vpa_dissipation_coefficient = 2.0
93+
#force_minimum_pdf_value = 0.0
94+
95+
[neutral_numerical_dissipation]
96+
#vz_dissipation_coefficient = 1.0e-1
97+
#vz_dissipation_coefficient = 1.0e-2
98+
#vz_dissipation_coefficient = 1.0e-3
99+
#force_minimum_pdf_value = 0.0

‎publication_inputs/2D-instability/input_files/2D1V-instability-periodic-test/2D1V-instability-periodic-Lz0.02-test-Fourier-r-AT0.1-An-0.01.toml‎

Lines changed: 99 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,99 @@
1+
[composition]
2+
T_wall = 0.1
3+
T_e = 0.2
4+
n_ion_species = 1
5+
n_neutral_species = 0
6+
7+
[evolve_moments]
8+
density = false
9+
parallel_flow = false
10+
pressure = false
11+
12+
[krook_collisions]
13+
use_krook = false
14+
15+
[geometry]
16+
option = "constant-helical"
17+
pitch = 0.1
18+
19+
[manufactured_solns]
20+
use_for_init = true
21+
use_for_advance = true
22+
type = "2D-instability"
23+
24+
[ion_species_1]
25+
initial_temperature = 0.2
26+
initial_density = 1.0
27+
28+
[z_IC_ion_species_1]
29+
initialization_option = "sinusoid"
30+
density_amplitude = -0.01
31+
temperature_amplitude = 0.1
32+
density_phase = 0.0
33+
upar_amplitude = 0.0
34+
temperature_phase = 0.0
35+
upar_phase = 0.0
36+
37+
[vpa_IC_ion_species_1]
38+
initialization_option = "gaussian"
39+
density_amplitude = 1.0
40+
temperature_amplitude = 0.0
41+
density_phase = 0.0
42+
upar_amplitude = 0.0
43+
temperature_phase = 0.0
44+
upar_phase = 0.0
45+
46+
[r]
47+
ngrid = 4
48+
discretization = "fourier_pseudospectral"
49+
nelement = 1
50+
bc = "periodic"
51+
L = 1.0e-4 # equivalent to 1mm
52+
53+
[z]
54+
ngrid = 5
55+
discretization = "gausslegendre_pseudospectral"
56+
nelement = 32
57+
#nelement_local = 8
58+
bc = "periodic"
59+
L = 0.02 # 1m vertical height with pitch=0.1 corresponds to 10m of parallel connection length.
60+
61+
[vpa]
62+
ngrid = 6
63+
discretization = "gausslegendre_pseudospectral"
64+
nelement = 31
65+
L = 10.0
66+
element_spacing_option = "coarse_tails8"
67+
bc = "zero"
68+
69+
[timestepping]
70+
type = "Fekete4(3)"
71+
nstep = 3000000
72+
dt = 1.0e-5
73+
minimum_dt = 1.0e-5
74+
CFL_prefactor = 2.0
75+
rtol = 1.0e-5
76+
atol = 1.0e-12
77+
add_random_noise = 1.0e-9
78+
random_noise_seed = 1234
79+
nwrite = 10000
80+
nwrite_dfns = 10000
81+
steady_state_residual = true
82+
converged_residual_value = 1.0e-12
83+
84+
[ion_numerical_dissipation]
85+
r_dissipation_coefficient = 0.0e-10
86+
#vpa_dissipation_coefficient = 1.0e-1
87+
vpa_dissipation_coefficient = 1.0e-2
88+
#vpa_dissipation_coefficient = 1.0e-3
89+
#force_minimum_pdf_value = 0.0
90+
91+
[electron_numerical_dissipation]
92+
#vpa_dissipation_coefficient = 2.0
93+
#force_minimum_pdf_value = 0.0
94+
95+
[neutral_numerical_dissipation]
96+
#vz_dissipation_coefficient = 1.0e-1
97+
#vz_dissipation_coefficient = 1.0e-2
98+
#vz_dissipation_coefficient = 1.0e-3
99+
#force_minimum_pdf_value = 0.0

‎publication_inputs/2D-instability/input_files/2D1V-instability-periodic-test/2D1V-instability-periodic-Lz0.02-test-Fourier-r-AT0.1-An-0.02.toml‎

Lines changed: 99 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,99 @@
1+
[composition]
2+
T_wall = 0.1
3+
T_e = 0.2
4+
n_ion_species = 1
5+
n_neutral_species = 0
6+
7+
[evolve_moments]
8+
density = false
9+
parallel_flow = false
10+
pressure = false
11+
12+
[krook_collisions]
13+
use_krook = false
14+
15+
[geometry]
16+
option = "constant-helical"
17+
pitch = 0.1
18+
19+
[manufactured_solns]
20+
use_for_init = true
21+
use_for_advance = true
22+
type = "2D-instability"
23+
24+
[ion_species_1]
25+
initial_temperature = 0.2
26+
initial_density = 1.0
27+
28+
[z_IC_ion_species_1]
29+
initialization_option = "sinusoid"
30+
density_amplitude = -0.02
31+
temperature_amplitude = 0.1
32+
density_phase = 0.0
33+
upar_amplitude = 0.0
34+
temperature_phase = 0.0
35+
upar_phase = 0.0
36+
37+
[vpa_IC_ion_species_1]
38+
initialization_option = "gaussian"
39+
density_amplitude = 1.0
40+
temperature_amplitude = 0.0
41+
density_phase = 0.0
42+
upar_amplitude = 0.0
43+
temperature_phase = 0.0
44+
upar_phase = 0.0
45+
46+
[r]
47+
ngrid = 4
48+
discretization = "fourier_pseudospectral"
49+
nelement = 1
50+
bc = "periodic"
51+
L = 1.0e-4 # equivalent to 1mm
52+
53+
[z]
54+
ngrid = 5
55+
discretization = "gausslegendre_pseudospectral"
56+
nelement = 32
57+
#nelement_local = 8
58+
bc = "periodic"
59+
L = 0.02 # 1m vertical height with pitch=0.1 corresponds to 10m of parallel connection length.
60+
61+
[vpa]
62+
ngrid = 6
63+
discretization = "gausslegendre_pseudospectral"
64+
nelement = 31
65+
L = 10.0
66+
element_spacing_option = "coarse_tails8"
67+
bc = "zero"
68+
69+
[timestepping]
70+
type = "Fekete4(3)"
71+
nstep = 3000000
72+
dt = 1.0e-5
73+
minimum_dt = 1.0e-5
74+
CFL_prefactor = 2.0
75+
rtol = 1.0e-5
76+
atol = 1.0e-12
77+
add_random_noise = 1.0e-9
78+
random_noise_seed = 1234
79+
nwrite = 10000
80+
nwrite_dfns = 10000
81+
steady_state_residual = true
82+
converged_residual_value = 1.0e-12
83+
84+
[ion_numerical_dissipation]
85+
r_dissipation_coefficient = 0.0e-10
86+
#vpa_dissipation_coefficient = 1.0e-1
87+
vpa_dissipation_coefficient = 1.0e-2
88+
#vpa_dissipation_coefficient = 1.0e-3
89+
#force_minimum_pdf_value = 0.0
90+
91+
[electron_numerical_dissipation]
92+
#vpa_dissipation_coefficient = 2.0
93+
#force_minimum_pdf_value = 0.0
94+
95+
[neutral_numerical_dissipation]
96+
#vz_dissipation_coefficient = 1.0e-1
97+
#vz_dissipation_coefficient = 1.0e-2
98+
#vz_dissipation_coefficient = 1.0e-3
99+
#force_minimum_pdf_value = 0.0

‎publication_inputs/2D-instability/input_files/2D1V-instability-periodic-test/2D1V-instability-periodic-Lz0.02-test-Fourier-r-AT0.1-An-0.03.toml‎

Lines changed: 99 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,99 @@
1+
[composition]
2+
T_wall = 0.1
3+
T_e = 0.2
4+
n_ion_species = 1
5+
n_neutral_species = 0
6+
7+
[evolve_moments]
8+
density = false
9+
parallel_flow = false
10+
pressure = false
11+
12+
[krook_collisions]
13+
use_krook = false
14+
15+
[geometry]
16+
option = "constant-helical"
17+
pitch = 0.1
18+
19+
[manufactured_solns]
20+
use_for_init = true
21+
use_for_advance = true
22+
type = "2D-instability"
23+
24+
[ion_species_1]
25+
initial_temperature = 0.2
26+
initial_density = 1.0
27+
28+
[z_IC_ion_species_1]
29+
initialization_option = "sinusoid"
30+
density_amplitude = -0.03
31+
temperature_amplitude = 0.1
32+
density_phase = 0.0
33+
upar_amplitude = 0.0
34+
temperature_phase = 0.0
35+
upar_phase = 0.0
36+
37+
[vpa_IC_ion_species_1]
38+
initialization_option = "gaussian"
39+
density_amplitude = 1.0
40+
temperature_amplitude = 0.0
41+
density_phase = 0.0
42+
upar_amplitude = 0.0
43+
temperature_phase = 0.0
44+
upar_phase = 0.0
45+
46+
[r]
47+
ngrid = 4
48+
discretization = "fourier_pseudospectral"
49+
nelement = 1
50+
bc = "periodic"
51+
L = 1.0e-4 # equivalent to 1mm
52+
53+
[z]
54+
ngrid = 5
55+
discretization = "gausslegendre_pseudospectral"
56+
nelement = 32
57+
#nelement_local = 8
58+
bc = "periodic"
59+
L = 0.02 # 1m vertical height with pitch=0.1 corresponds to 10m of parallel connection length.
60+
61+
[vpa]
62+
ngrid = 6
63+
discretization = "gausslegendre_pseudospectral"
64+
nelement = 31
65+
L = 10.0
66+
element_spacing_option = "coarse_tails8"
67+
bc = "zero"
68+
69+
[timestepping]
70+
type = "Fekete4(3)"
71+
nstep = 3000000
72+
dt = 1.0e-5
73+
minimum_dt = 1.0e-5
74+
CFL_prefactor = 2.0
75+
rtol = 1.0e-5
76+
atol = 1.0e-12
77+
add_random_noise = 1.0e-9
78+
random_noise_seed = 1234
79+
nwrite = 10000
80+
nwrite_dfns = 10000
81+
steady_state_residual = true
82+
converged_residual_value = 1.0e-12
83+
84+
[ion_numerical_dissipation]
85+
r_dissipation_coefficient = 0.0e-10
86+
#vpa_dissipation_coefficient = 1.0e-1
87+
vpa_dissipation_coefficient = 1.0e-2
88+
#vpa_dissipation_coefficient = 1.0e-3
89+
#force_minimum_pdf_value = 0.0
90+
91+
[electron_numerical_dissipation]
92+
#vpa_dissipation_coefficient = 2.0
93+
#force_minimum_pdf_value = 0.0
94+
95+
[neutral_numerical_dissipation]
96+
#vz_dissipation_coefficient = 1.0e-1
97+
#vz_dissipation_coefficient = 1.0e-2
98+
#vz_dissipation_coefficient = 1.0e-3
99+
#force_minimum_pdf_value = 0.0

0 commit comments

Comments
 (0)