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

Adding nlte solver #2171

Merged
merged 38 commits into from
Jan 19, 2023
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
38 commits
Select commit Hold shift + click to select a range
eb7f5cb
adding solver part in the calculate method
sonachitchyan Dec 7, 2022
3e10cc0
adding the objective function
sonachitchyan Dec 7, 2022
87caec5
added the solution vector method
sonachitchyan Dec 7, 2022
ddc3f11
final touches for the solver
sonachitchyan Dec 7, 2022
82bf704
tests, attempt 1
sonachitchyan Dec 8, 2022
73043db
ran black
sonachitchyan Dec 8, 2022
e0ac106
fixed the issue with tests
sonachitchyan Jan 9, 2023
92444bd
ran black
sonachitchyan Jan 9, 2023
3ac7629
added missing doctrings in the nlte_rate_equation_solver.py
sonachitchyan Jan 9, 2023
e67976a
added doctrings to tests
sonachitchyan Jan 9, 2023
572914a
changed the test to explicitly calculate the lte solution ion number …
sonachitchyan Jan 9, 2023
33553b0
ran black
sonachitchyan Jan 9, 2023
58ce29e
added the atom data file to download_reference_data.sh
sonachitchyan Jan 9, 2023
42b136d
fixed download_reference_data.sh
sonachitchyan Jan 9, 2023
2765722
Revert "added the atom data file to download_reference_data.sh"
sonachitchyan Jan 9, 2023
59f4db8
got rid of index i, kept only shell
sonachitchyan Jan 12, 2023
dceb551
switched DataFrame to pandas.DataFrame in the docstrings
sonachitchyan Jan 12, 2023
e170c63
docstring bug fix
sonachitchyan Jan 12, 2023
b6d90d3
got rid of the deepcopy of nlte atomic data file
sonachitchyan Jan 12, 2023
f6169b6
changed the number of shells to 5
sonachitchyan Jan 12, 2023
789fd00
Merge branch 'master' of github.com:tardis-sn/tardis into adding_nlte…
sonachitchyan Jan 12, 2023
b129081
Update tardis/io/tests/data/tardis_configv1_nlte.yml
sonachitchyan Jan 12, 2023
16d3d0f
Update tardis/plasma/properties/nlte_rate_equation_solver.py
sonachitchyan Jan 12, 2023
33a6562
Update tardis/io/tests/data/tardis_configv1_nlte.yml
sonachitchyan Jan 12, 2023
63d222f
Update tardis/io/tests/data/tardis_configv1_nlte.yml
sonachitchyan Jan 12, 2023
7076c06
ran black
sonachitchyan Jan 12, 2023
fe3ac0c
Merge branch 'adding_nlte_solver' of github.com:tardis-sn/tardis into…
sonachitchyan Jan 12, 2023
c63fa0c
Update tardis/plasma/properties/nlte_rate_equation_solver.py
sonachitchyan Jan 12, 2023
b9a5eec
Update tardis/plasma/tests/test_nlte_solver.py
sonachitchyan Jan 12, 2023
4f28528
restructured a summary in a docstring for the solution_vector_block
sonachitchyan Jan 12, 2023
884a67a
ran black
sonachitchyan Jan 12, 2023
67ad1bd
added the test for w=0 case
sonachitchyan Jan 17, 2023
3f5b6b9
Update tardis/plasma/properties/nlte_rate_equation_solver.py
sonachitchyan Jan 18, 2023
7b5f504
Update tardis/plasma/properties/nlte_rate_equation_solver.py
sonachitchyan Jan 18, 2023
34b3554
Update tardis/plasma/tests/test_nlte_solver.py
sonachitchyan Jan 18, 2023
07a6faa
fixed an issue in a docsting
sonachitchyan Jan 18, 2023
f36aa84
fixed an issue in a docsting
sonachitchyan Jan 18, 2023
352e3c9
removed a TODO comment
sonachitchyan Jan 18, 2023
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
100 changes: 58 additions & 42 deletions tardis/io/tests/data/tardis_configv1_nlte.yml
Original file line number Diff line number Diff line change
@@ -1,54 +1,70 @@
tardis_config_version: v1.0

supernova:
luminosity_requested: 2.8e9 solLum
time_explosion: 13 day
luminosity_requested: 9.05 log_lsun
time_explosion: 16 day

atom_data: kurucz_atom_pure_simple.h5
atom_data: TestNLTE_Ti.h5

model:
structure:
type: specific
velocity:
start: 1.1e4 km/s
stop: 2.0e4 km/s
num: 20
density:
type: branch85_w7
abundances:
type: uniform
H: 0.1
O: 0.09
Mg: 0.03
Si: 0.52
S: 0.19
Ar: 0.04
Ca: 0.03
structure:
type: specific
velocity:
start: 6395 km/s
stop: 12500 km/s
num: 5
density:
type : power_law
time_0: 16.0 day
rho_0: 1.948e-14 g/cm^3
v_0: 8000 km/s
exponent: -10

abundances:
type: uniform
H: 0.9
He: 0.099
Ti: 1e-3

plasma:
ionization: lte
excitation: lte
radiative_rates_type: dilute-blackbody
line_interaction_type: macroatom
continuum_interaction:
species:
- H I
nlte_ionization_species: [H I]
disable_electron_scattering: no
ionization: lte
excitation: lte
radiative_rates_type: dilute-blackbody
line_interaction_type: macroatom
initial_t_inner: 12000 K
link_t_rad_t_electron: 1.0
nlte_ionization_species: [H I, He II, Ti II]

continuum_interaction:
species:
- H I
- Ti II
- He II

enable_adiabatic_cooling: True



montecarlo:
seed: 23111963
no_of_packets : 2.0e+5
iterations: 5
last_no_of_packets: 5.0e+5
no_of_virtual_packets: 5
convergence_strategy:
type: damped
damping_constant: 0.5
threshold: 0.05
lock_t_inner_cycles: 1
t_inner_update_exponent: -0.5
seed: 23111963
no_of_packets: 100000
iterations: 1
nthreads: 1

last_no_of_packets: 1000
no_of_virtual_packets: 0

convergence_strategy:
type: damped
damping_constant: 0.5
threshold: 0.05
fraction: 0.8
hold_iterations: 3

spectrum:
start: 500 angstrom
stop: 20000 angstrom
num: 10000
start: 250 angstrom
stop: 10000 angstrom
num: 10000
method: integrated
compute: Automatic
Loading