tardis.energy_input.gamma_ray_transport module

tardis.energy_input.gamma_ray_transport.calculate_average_energies(raw_isotope_abundance, gamma_ray_lines)[source]

Function to calculate average energies of positrons and gamma rays from a list of gamma ray lines from nndc.

Parameters:
raw_isotope_abundancepd.DataFrame

isotope abundance in mass fractions

gamma_ray_linespd.DataFrame

decay data

Returns:
average_energies_listList

list of gamma ray energies

average_positron_energies_listList

list of positron energies

gamma_ray_line_array_listList

list of gamma ray lines

tardis.energy_input.gamma_ray_transport.calculate_average_power_per_mass(energy_per_mass, time_delta)[source]
tardis.energy_input.gamma_ray_transport.calculate_ejecta_velocity_volume(model)[source]
tardis.energy_input.gamma_ray_transport.calculate_energy_per_mass(decay_energy, raw_isotope_abundance, cell_masses)[source]

Function to calculate decay energy per mass for each isotope chain.

Parameters:
decay_energyDict

dictionary of decay chain energies for each isotope in each shell

raw_isotope_abundancepd.DataFrame

isotope abundance in mass fractions.

cell_massesnumpy.ndarray

shell masses in units of g

Returns:
energy_per_masspd.DataFrame

decay energy per mass for each isotope chain For e.g Ni56 has 2 decay chains: Ni56 -> Co56 -> Fe56. It will calculate the decay energy per mass for each chain and store as a dataframe.

tardis.energy_input.gamma_ray_transport.calculate_total_decays_old(inventories, time_delta)[source]

Function to create inventories of isotope

Parameters:
modeltardis.Radial1DModel

The tardis model to calculate gamma ray propagation through

time_endfloat

End time of simulation in days

Returns:
Total decay listList

list of total decays for x g of isotope for time ‘t’

tardis.energy_input.gamma_ray_transport.create_inventories_dict_old(isotope_dict)[source]

Function to create dictionary of inventories for each shell

Parameters:
isotope_dictDict

dictionary of isotopes for each shell with their masses.

Returns:
invDict

dictionary of inventories for each shell For eg: {0: {‘Ni56’: <radioactivedecay.Inventory at 0x7f7d2c0d0e50>,

‘Co56’: <radioactivedecay.Inventory at 0x7f7d2c0d0e50>},

{1: {‘Ni56’: <radioactivedecay.Inventory at 0x7f7d2c0d0e50>,

‘Co56’: <radioactivedecay.Inventory at 0x7f7d2c0d0e50>}} etc

tardis.energy_input.gamma_ray_transport.create_isotope_dicts_old(raw_isotope_abundance, cell_masses)[source]

Function to create a dictionary of isotopes for each shell with their masses.

Parameters:
raw_isotope_abundancepd.DataFrame

isotope abundance in mass fractions.

cell_massesnumpy.ndarray

shell masses in units of g

Returns:
isotope_dictsDict

dictionary of isotopes for each shell with their masses. For eg: {0: {(28, 56): {‘Ni56’: 0.0001}, (27, 57): {‘Co56’: 0.0001}}

{1: {(28, 56): {‘Ni56’: 0.0001}, (27, 57): {‘Co56’: 0.0001}}} etc

tardis.energy_input.gamma_ray_transport.decay_chain_energies(average_energies, total_decays)[source]

Function to calculate decay chain energies.

Parameters:
raw_isotope_abundancepd.DataFrame

isotope abundance in mass fractions

average_energies_listList

list of gamma ray energies

average_positron_energies_listList

list of positron energies

gamma_ray_line_array_listList

list of gamma ray lines

total_decaysDict

dictionary of total decays for each isotope in each shell

Returns:
decay_energyDict

dictionary of decay chain energies for each isotope in each shell

tardis.energy_input.gamma_ray_transport.distribute_packets(decay_energy, total_energy, num_packets)[source]
tardis.energy_input.gamma_ray_transport.fractional_decay_energy(decay_energy)[source]

Function to calculate fractional decay energy

Parameters:
decay_energyDict

dictionary of decay chain energies for each isotope in each shell

Returns:
fractional_decay_energyDict

dictionary of fractional decay chain energies for each isotope in each shell

tardis.energy_input.gamma_ray_transport.get_chain_decay_power_per_ejectamass(inventory, time, time_start, average_energies, average_positron_energies, tau_end)[source]
tardis.energy_input.gamma_ray_transport.get_nuclide_atomic_number(nuclide)[source]
tardis.energy_input.gamma_ray_transport.get_taus(raw_isotope_abundance)[source]

Function to calculate taus for each isotope

Parameters:
raw_isotope_abundancepd.DataFrame

isotope abundance in mass fractions

Returns:
tausDict

dictionary of taus for each isotope

parentsDict

dictionary of parents for each isotope

tardis.energy_input.gamma_ray_transport.iron_group_fraction_per_shell(model)[source]
tardis.energy_input.gamma_ray_transport.packets_per_isotope(fractional_decay_energy, decayed_packet_count_dict)[source]