Obsah
Konfigurační soubor
Specifikace
Každá simulační úloha je definována v jediném konfiguračním souboru ve formátu YAML.
Vstupní jednotky jsou základní jednotky SI, ale čas je v letech.
simulation_params
Základní nastavení úlohy
dt
Maximální časový krok simulace [roky]
Příklad: dt: 5
time
Celkový čas simulace [roky]
Příklad: time: 2000000
max_iter
Maximální počet iterací. Pokud je dosaženo tohoto limitu, simulace se ukončí (pravděpodobně numerická nestabilita).
Příklad: max_iter: 24
include_processes
Seznam fyzikálních procesů zahrnutých v simulaci
Možnosti položek seznamu:
conversion– Radioaktivní přeměny radionuklidů.max_solubility– Zohlednění maximální rozpustnosti jednotlivých prvků (jinak by byla neomezená).diffusion– Difúze radionuklidů v bentonitové bariéře.sorption– Sorpce radionuklidů v bentonitové bariéře (lineární Henryho izoterma).
output
Definuje formát výstupních dat
step_time
Perioda časových výstupů [roky]
Příklad: step_time: 500
quantities
Veličiny (sloupečky ve výstupním CSV)
quantities:
fuel_m_kg– Celkové VJP (v tyči, tj. nezdegradované, nerozpuštěné) [kg]fuel_n_mol– Celkové VJP [mol]fuel_A_Bq– Celkové VJP [Bq]deg_m_kg– Zdegradované VJP „písek“ [kg]deg_n_mol– Zdegradované VJP [mol]fuel+deg_m_kg– Celkové VJP (zdegradované i nezdegradované) [kg]cont.water_c_kg/m3– Kontejner: Koncentrace ve vodě [kg/m3]cont.water_sat_%– Kontejner: Koncentrace ve vodě [% rozpustnosti]cont.water_c_mol/m3– Kontejner: Koncentrace ve vodě [mol/m3]cont.water_A_Bq– Kontejner: Aktivita ve vodě [Bq]cont.water_AV_Bq/m3– Kontejner: Aktivita ve vodě [Bq/m3]bent.inner.bound.water_c_kg/m3– Vnitřní okraj bentonitu: Koncentrace ve vodě [kg/m3]bent.inner.bound.water_c_mol/m3– Vnitřní okraj bentonitu: Koncentrace ve vodě [mol/m3]bent.1/4.water_c_kg/m3– 1/4 bentonitu (width*0.25) od vnitřního okraje: Koncentrace ve vodě [kg/m3]bent.1/4.water_c_mol/m3– 1/4 bentonitu (width*0.25) od vnitřního okraje: Koncentrace ve vodě [mol/m3]bent.1/2.water_c_kg/m3– 1/2 bentonitu (width*0.50) od vnitřního okraje: Koncentrace ve vodě [kg/m3]bent.1/2.water_c_mol/m3– 1/2 bentonitu (width*0.50) od vnitřního okraje: Koncentrace ve vodě [mol/m3]bent.3/4.water_c_kg/m3– 3/4 bentonitu (width*0.75) od vnitřního okraje: Koncentrace ve vodě [kg/m3]bent.3/4.water_c_mol/m3– 3/4 bentonitu (width*0.75) od vnitřního okraje: Koncentrace ve vodě [mol/m3]bent.outer.bound.water_c_kg/m3– Vnější okraj bentonitu: Koncentrace ve vodě [kg/m3]bent.outer.bound.water_c_mol/m3– Vnější okraj bentonitu: Koncentrace ve vodě [mol/m3]bent.layers.water_c_mol/m3– Koncentrace ve vodě v jednotlivých vrstvách bentonitu [mol/m3] (defaultně po 5 cm)bent.layers.water_c_kg/m3– Koncentrace ve vodě v jednotlivých vrstvách bentonitu [kg/m3] (defaultně po 5 cm)bent.water_m_kg– Hmotnost rozpuštěné látky v bentonitu [kg]bent.rock_m_kg– Hmotnost sorbované látky v bentonitu [kg]bent_m_kg– Celkové množství radionuklidů v bentonitu [kg]granit.water_m_kg– Kumulativní výtok z bentonitu [kg]granit.water_n_mol– Kumulativní výtok z bentonitu [mol]granit.water_m(cumsum)_kg– Kumulativní výtok bez radioaktivního rozpadu [kg]granit.water_n(cumsum)_mol– Kumulativní výtok bez radioaktivního rozpadu [mol]
isotopes
Nuklidy (sloupečky ve výstupním CSV)
Příklad: isotopes: [all]
elements
Prvky (sloupečky ve výstupním CSV)
Příklad: elements: [all]
summary_quantity
Ukládá součty všech prvků/izotopů
Příklad: summary_quantity: no
degradation_rate
Ukládá výstup fuel_degradation_function
Příklad: degradation_rate: yes
columns_ordering
Určuje pořadí sloupců ve výstupním CSV souboru.
Mount je amount (množství).
Možnosti:
quantity– Sloupce jsou seřazeny podle veličin.quantity/name– Nejprve veličina, poté název izotopu.quantity/initial_mount– Nejprve veličina, poté počáteční množství izotopu.name/quantity– Nejprve název izotopu, poté veličina.initial_mount/quantity– Nejprve počáteční množství izotopu, poté veličina.
Příklad: columns_ordering: quantity
separators
Oddělovač sloupců.
Možnosti:
dec._col,– Desetinný oddělovač je., oddělovač sloupců je,.dec._col;– Desetinný oddělovač je., oddělovač sloupců je;.dec,_col;– Desetinný oddělovač je,, oddělovač sloupců je;.dec._colt– Desetinný oddělovač je., oddělovač sloupců je tabelátor (tab).
Příklad: separators: dec,_col;
number_format
Formát čísel
Možnosti:
G– general (obecný formát, např.0.0)E– exponential (exponenciální formát, např.0.0E+0)4E– dtto se 4 desetinnými místy
F– fixed (formát s pevním počtem desetinných míst)4F– dtto se 4 desetinnými místy
Příklad: number_format: G
file_name
Cesta k výstupnímu souboru CSV.
Příklad: file_name: .\edu.csv
statistics
Výběr statistických výstupů simulace.
Možnosti:
cont_max_c– Maximální koncentrace radionuklidů v kontejneru.bent_bound_max_c– Maximální koncentrace radionuklidů na vnějším okraji bentonitové bariéry.granit_inflow– Kumulativní výtok radionuklidů do horninového prostředí.
container
Parametry úložného kontejneru.
radius
Vnější poloměr kontejneru (jen kvůli geometrii bentonitu)
Příklad: radius: 0.4
height
Výška kontejneru (jen kvůli geometrii bentonitu)
Příklad: height: 4
free_volume
Vnitřní volný objem kontejneru
Příklad: free_volume: 0.2
T0_water_inflow_time
Čas od začátku simulace do průniku vody do kontejneru (začátek mokré fáze) [roky]
Příklad: T0_water_inflow_time: 2e5
fuel_degradation_function
Definice kinetiky degradace VJP.
model
Přepínač
Možnosti:
user_def– Ručně definovaná křivka degradace (sekcepowerse ignoruje).power– Mocninná funkce degradace (sekceuser_defse ignoruje).
user_def
Lineární lomená funkce
times
Čas od T0_water_inflow_time. Musí obsahovat aspoň dvě hodnoty; první musí být nula.
Příklad: times: [0, 1000000]
values_%
Procenta VJP, které je zdegradováno.
- První hodnota musí být
irf(indikuje počáteční hodnotu spočítanou jinde), - Poslední hodnota musí být 100 (vše bylo zdegradováno).
Příklad: values_%: [irf, 100]
power
Mocninná funkce
degradation_time_span
V tomto čase funkce nabývá 100%
Příklad: degradation_time_span: 8e5
power_factor
Mocnina
Příklad: power_factor: 0.3
time_step
Časová diskretizace (výpočet podle lomenné čáry)
Příklad: time_step: 100
dissolving
Rozpouštění
controlled
Možnosti:
no– všechny prvky se rozpouští do úrovně max. rozpustnosti pro každý prvekyes– navíc je omezeno poměrem kecontrolling_element
Příklad: controlled: yes
controlling_element
Prvek řídící rozpustnost ostatních prvků
Příklad: controlling_element: U
bentonite_contact_area_function
Funkce časového průběhu velikosti plochy odstraněného pláště kontejneru — lineární lomená funkce
times
Čas – první hodnota musí být nula
Příklad: times: [0]
values_%
Hodnota libovolné, poslední hodnota zůstává po zbytek simulace
Příklad: values_%: [100]
direct_outflow_into_granite_%/year
Tok prvků rozp. v kontejnerové vodě [% v kont. vodě / rok]
Příklad: direct_outflow_into_granite_%/year: 0
bentonite
width
Tloušťka bentonitové bariéry [m]
Příklad: width: 0.7
dr
Diskretizace tloušťky [m]
Příklad: dr: 0.05
density
Hustota suchého bentonitu [kg/m3]
Příklad: density: 1500
zones_outer_bound
Vnější hranice zón (od vnitřního okraje)
Příklad: zones_outer_bound: [0.7]
Následující 3 položky definují vlastnosti v zónách (lineární lomená funkce):
porosity,
De_b_modification a
Kd_b_modification.
porosity
Pórovitost
times
Příklad: times: [0]
zone1_%
Takových může být víc
Příklad: zone1_%: [30]
De_b_modification
times
Příklad: times: [0]
zone1_%
Příklad: zone1_%: [100]
Kd_b_modification
times
Příklad: times: [0]
zone1_%
Příklad: zone1_%: [100]
bound_model
Okrajová podmínka pro transport na vnějším okraji bentonitu.
Možnosti:
granite_layer– Přechod do horninové vrstvy.zero_conc– Koncentrace na vnější hranici je nulová (tj. okamžité odstranění látek z okraje bentonitu)
Příklad: bound_model: zero_conc
granite_layer
Jen pokud: bound_model: granite_layer
width
Tloušťka horninové vrstvy
Příklad: width: 0.002
De_source
Kde je zadán efektivní difuzní koeficient
Možnosti:
De_b_%: definován na řádku nížeDe_g: definován v bloku isotopes (ještě níže)
Příklad: De_source: De_b_%
De_b_%
Efektivní difuzní koeficient De v granitu zadán jako procento De v bentonitu.
Příklad: De_b_%: 100
conc_by_bentonite_bound_conc_%
Koncentrace na vnějším okraji granitu (zadáno jako procenta koncentrace na vnitřním okraji)
Příklad: conc_by_bentonite_bound_conc_%: 0
isotopes
Popis jednotlivých izotopů v simulaci.
name
Název izotopu (musí být: prvek a nukleové číslo)
Příklad: name: C14
M_kg/mol
Molární hmotnost izotopu [kg/mol]
Příklad: M_kg/mol: 0.014
T_1/2_year
Poločas rozpadu izotopu [roky]
Příklad: T_1/2_year: 5700
conversion_to
Produkt rozpadu (musí být v tomto seznamu izotopů)
Příklad: conversion_to: odpad
init_m_kg
Počáteční hmotnost izotopu ve VJP
Možnosti:
init_a_bq– Počáteční aktivita izotopu ve vyhořelém palivu [Bq]. Udává se jako počet radioaktivních přeměn za sekundu.init_n_mol– Počáteční množství izotopu ve vyhořelém palivu [mol]. Udává se v látkovém množství.
Příklad: init_m_kg: 0.000882
irf_%
Počáteční okamžitě uvolnitelná frakce (procenta z init_m_kg)
Příklad: irf_%: 10
c_max_mol/m3
Rozpustnost [mol/m3]
Příklad: c_max_mol/m3: inf
De_b_m2/s
Efektivní difuzní koeficient De v bentonitu [m2/s]
Příklad: De_b_m2/s: 1.2e-10
Kd_b_m3/kg
Distribuční koeficient Kd v bentonitu [m3/kg]
Příklad: Kd_b_m3/kg: 0
De_g_m2/s
Efektivní difuzní koeficient De v granitu [m2/s]
Příklad: De_g_m2/s: 5e-15