Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Multiple source functionality #250

Merged
merged 24 commits into from
Sep 13, 2024
Merged
Changes from 1 commit
Commits
Show all changes
24 commits
Select commit Hold shift + click to select a range
af49432
Split get_settings function for neutrals and ions, in preparation for…
LucasMontoya4 Sep 6, 2024
0c67d26
Convert form of ion_settings, electron_settings and neutral_settings …
LucasMontoya4 Sep 6, 2024
372988a
Add options to ion and neutral source profile type error
LucasMontoya4 Sep 6, 2024
4a74181
Change names of structs from 'profile' to 'data'
LucasMontoya4 Sep 6, 2024
a381240
Change create_moments_ion and create_moments_neutral functions to be …
LucasMontoya4 Sep 7, 2024
fb0b082
Restructure main initialisation files and update loops to view each s…
LucasMontoya4 Sep 8, 2024
2e5527a
Update file_io.jl to write out external_source_amplitudes for each mo…
LucasMontoya4 Sep 8, 2024
0cf1c36
Update load_data.jl for multiple sources
LucasMontoya4 Sep 8, 2024
0362f21
Update source functionality for all moment kinetic advances of densit…
LucasMontoya4 Sep 8, 2024
41759d9
clean up external_sources.jl
LucasMontoya4 Sep 8, 2024
2622a24
Clean up small typos
LucasMontoya4 Sep 8, 2024
6123063
Update all example input files and test files
LucasMontoya4 Sep 8, 2024
dddad29
Fix bug in velocity_moments.jl PI_density_controller loop
LucasMontoya4 Sep 9, 2024
da0f366
Add exponential wall decay profile option to simulate neutral ionisat…
LucasMontoya4 Sep 9, 2024
9bfb631
Temporary fix for file_io.jl, where only the first source profile of …
LucasMontoya4 Sep 9, 2024
25bdbfe
Add collisionality_scan option to krook operator for ions
LucasMontoya4 Sep 9, 2024
8dd9025
Allow NamedTuple 'coords' in create_dynamic_variable!()
johnomotani Sep 10, 2024
69b80ba
Merge remote-tracking branch 'origin/simplify-io-variable-creation' i…
LucasMontoya4 Sep 10, 2024
8d61d1e
Allow output of all sources as one variable (after John's update to a…
LucasMontoya4 Sep 11, 2024
d4dd565
bug fixes for plotting sources in makie_post_processing
LucasMontoya4 Sep 11, 2024
e8f2926
Use ion_sources directly instead of `counter` to set up electron sources
johnomotani Sep 12, 2024
150d05a
Fix import of external source function in electron_kinetic_equation
johnomotani Sep 12, 2024
fe922fa
Fix typo in external_electron_source!()
johnomotani Sep 12, 2024
5ecabac
Fix output of neutral source terms
johnomotani Sep 12, 2024
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
Update all example input files and test files
LucasMontoya4 committed Sep 8, 2024
commit 6123063d5ed901c9ce516d746f3980b748a5b663
Original file line number Diff line number Diff line change
@@ -65,7 +65,7 @@ sd_mi = 0.5 # mD/malpha
sd_me = 0.000013616 # 0.25/1836.0 me/malpha
sd_q = 1.0

[ion_source]
[ion_source_1]
active=false
source_strength=0.0
source_T=0.005
Original file line number Diff line number Diff line change
@@ -64,7 +64,7 @@ sd_mi = 0.5 # mD/malpha
sd_me = 0.000013616 # 0.25/1836.0 me/malpha
sd_q = 1.0

[ion_source]
[ion_source_1]
active=true
source_strength=1.0
source_T=0.005
Original file line number Diff line number Diff line change
@@ -64,7 +64,7 @@ sd_mi = 0.5 # mD/malpha
sd_me = 0.000013616 # 0.25/1836.0 me/malpha
sd_q = 1.0

[ion_source]
[ion_source_1]
active=true
source_strength=1.0
source_T=0.005
Original file line number Diff line number Diff line change
@@ -74,7 +74,7 @@ nwrite_dfns = 100000
steady_state_residual = true
converged_residual_value = 1.0e-3

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -74,7 +74,7 @@ nwrite_dfns = 100000
steady_state_residual = true
converged_residual_value = 1.0e-3

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -84,7 +84,7 @@ steady_state_residual = true
converged_residual_value = 1.0e-3
#write_after_fixed_step_count = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -90,7 +90,7 @@ minimum_dt = 1.0e-9
initialization_residual_value = 2.5
converged_residual_value = 0.1 #1.0e-3

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -75,7 +75,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -80,7 +80,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -80,7 +80,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -80,7 +80,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -75,7 +75,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -80,7 +80,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -74,7 +74,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -82,7 +82,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -82,7 +82,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -82,7 +82,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -82,7 +82,7 @@ converged_residual_value = 1.0e-3
[krook_collisions]
use_krook = true

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
2 changes: 1 addition & 1 deletion examples/wall-bc/wall-bc_no-neutrals.toml
Original file line number Diff line number Diff line change
@@ -70,7 +70,7 @@ nwrite = 5000
nwrite_dfns = 5000
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
2 changes: 1 addition & 1 deletion examples/wall-bc/wall-bc_no-neutrals_split1.toml
Original file line number Diff line number Diff line change
@@ -70,7 +70,7 @@ nwrite = 5000
nwrite_dfns = 5000
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
2 changes: 1 addition & 1 deletion examples/wall-bc/wall-bc_no-neutrals_split2.toml
Original file line number Diff line number Diff line change
@@ -70,7 +70,7 @@ nwrite = 5000
nwrite_dfns = 5000
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
2 changes: 1 addition & 1 deletion examples/wall-bc/wall-bc_no-neutrals_split3.toml
Original file line number Diff line number Diff line change
@@ -70,7 +70,7 @@ nwrite = 5000
nwrite_dfns = 5000
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
4 changes: 2 additions & 2 deletions examples/wall-bc/wall-bc_volumerecycle.toml
Original file line number Diff line number Diff line change
@@ -72,14 +72,14 @@ nwrite = 10000
nwrite_dfns = 10000
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
source_strength = 2.0
source_T = 2.0

[neutral_source]
[neutral_source_1]
active = true
source_type = "recycling"
recycling_controller_fraction = 0.25
4 changes: 2 additions & 2 deletions examples/wall-bc/wall-bc_volumerecycle_split1.toml
Original file line number Diff line number Diff line change
@@ -72,14 +72,14 @@ nwrite = 10000
nwrite_dfns = 10000
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
source_strength = 2.0
source_T = 2.0

[neutral_source]
[neutral_source_1]
active = true
source_type = "recycling"
recycling_controller_fraction = 0.25
4 changes: 2 additions & 2 deletions examples/wall-bc/wall-bc_volumerecycle_split2.toml
Original file line number Diff line number Diff line change
@@ -72,14 +72,14 @@ nwrite = 10000
nwrite_dfns = 10000
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
source_strength = 2.0
source_T = 2.0

[neutral_source]
[neutral_source_1]
active = true
source_type = "recycling"
recycling_controller_fraction = 0.25
4 changes: 2 additions & 2 deletions examples/wall-bc/wall-bc_volumerecycle_split3.toml
Original file line number Diff line number Diff line change
@@ -72,14 +72,14 @@ nwrite = 10000
nwrite_dfns = 10000
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
source_strength = 2.0
source_T = 2.0

[neutral_source]
[neutral_source_1]
active = true
source_type = "recycling"
recycling_controller_fraction = 0.25
2 changes: 1 addition & 1 deletion moment_kinetics/debug_test/kinetic_electron_inputs.jl
Original file line number Diff line number Diff line change
@@ -81,7 +81,7 @@ test_input = Dict("n_ion_species" => 1,
"vz_L" => 6.0,
"vz_bc" => "zero",
"vz_discretization" => "chebyshev_pseudospectral",
"ion_source" => Dict{String,Any}("active" => true,
"ion_source_1" => Dict{String,Any}("active" => true,
"z_profile" => "gaussian",
"z_width" => 0.125,
"source_strength" => 2.0,
2 changes: 1 addition & 1 deletion moment_kinetics/debug_test/recycling_fraction_inputs.jl
Original file line number Diff line number Diff line change
@@ -73,7 +73,7 @@ test_input = Dict("n_ion_species" => 1,
"vz_L" => 6.0,
"vz_bc" => "zero",
"vz_discretization" => "chebyshev_pseudospectral",
"ion_source" => Dict("active" => true,
"ion_source_1" => Dict("active" => true,
"z_profile" => "gaussian",
"z_width" => 0.125,
"source_strength" => 2.0,
2 changes: 1 addition & 1 deletion moment_kinetics/test/harrisonthompson.jl
Original file line number Diff line number Diff line change
@@ -112,7 +112,7 @@ test_input_finite_difference = Dict("n_ion_species" => 1,
"vz_L" => 8.0,
"vz_bc" => "zero",
"vz_discretization" => "finite_difference",
"ion_source" => Dict("active" => true,
"ion_source_1" => Dict("active" => true,
"source_strength" => ionization_frequency,
"source_T" => 0.25,
"z_profile" => "constant",
2 changes: 1 addition & 1 deletion moment_kinetics/test/recycling_fraction_tests.jl
Original file line number Diff line number Diff line change
@@ -83,7 +83,7 @@ test_input = Dict("n_ion_species" => 1,
"vz_L" => 18.0,
"vz_bc" => "zero",
"vz_discretization" => "chebyshev_pseudospectral",
"ion_source" => Dict("active" => true,
"ion_source_1" => Dict("active" => true,
"z_profile" => "gaussian",
"z_width" => 0.125,
"source_strength" => 2.0,
Original file line number Diff line number Diff line change
@@ -75,7 +75,7 @@ nwrite_dfns = 10000
n_rk_stages = 4
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -75,7 +75,7 @@ nwrite_dfns = 10000
n_rk_stages = 4
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -75,7 +75,7 @@ nwrite_dfns = 10000
n_rk_stages = 4
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125
Original file line number Diff line number Diff line change
@@ -75,7 +75,7 @@ nwrite_dfns = 10000
n_rk_stages = 4
split_operators = false

[ion_source]
[ion_source_1]
active = true
z_profile = "gaussian"
z_width = 0.125