230 likes | 265 Views
NUMERICKÉ SCHÉMA TVORBY DCEŘINNÝCH NUKLIDŮ BĚHEM POSTUPU RADIOAKTIVNÍHO MRAKU NAD TERÉNEM Petr Pecha Ústav teorie informace a automatizace AV ČR, odd. AS pecha@utia.cas.cz. parent (p) daughter (d) stable
E N D
NUMERICKÉ SCHÉMA TVORBY DCEŘINNÝCH NUKLIDŮ BĚHEM POSTUPU RADIOAKTIVNÍHO MRAKU NAD TERÉNEM Petr Pecha Ústav teorie informace a automatizace AV ČR, odd. AS pecha@utia.cas.cz XXIX. DNY RADIAČNÍ OCHRANY
parent (p) daughter (d) stable • Kr88(2,8 hod)Rb88(17,8 min)Sr88 Te132(78,2 hod)I132(2,3 hod)Xe132 • B)izomerní přechod : • d mY91 mTc99 m • • p d ……Sr91 Y91 …… Mo99 Tc99 …… • C)atd atd atd …… Metastable isomers of a particular isotope are usually designated with an "m„- excitace 1 nebo více nukleonů XXIX. DNY RADIAČNÍ OCHRANY
Jednoduchá genetická souvislost RODIČ DCERA Formulace: (izolovaná soustava) Řešení: Různá asymptotická přiblížení: Tp Td …….. Tp Td ……. …. atd. …. XXIX. DNY RADIAČNÍ OCHRANY
Transport aktivity v proudících vzdušných masách Momentové rovnice + difúze (teplota, vodní pára) + difúze aktivity + kontinuita + stavová rov. 8 rovnic pro 8 neznámých: vx, vy, vz, p, , , s, a • zdroje (bodové či lineární úniky do atm., příspěvky rodiče) a propady (rad. rozpad, suché vypadávání, vymývání) ------------------------------------------------------------------------------- PRAXE ŘEŠENÍ: generovaná meteopředpověď vstupuje jako počáteční podmínka do rovnice bilance aktivity XXIX. DNY RADIAČNÍ OCHRANY
MATCH Výpočetní oblast modelu ALADIN XXIX. DNY RADIAČNÍ OCHRANY
Speciální řešení difúzní rovnice ve tvaruGaussovy rovnice šíření má tvar: GPM – Gaussian Plume Model faktor ochuzení XXIX. DNY RADIAČNÍ OCHRANY
Příklad 1: ETE: Přímočarý únik směr SV – aktivita v přízemní vrstvě vzduchu Kategorie počasí = F; Rychlost větru (m/s) = 1.31; Trvání a výška úniku 3600 sec; 99.0 m Vzdalenost nadm. vysky drsnosti srazky landuse (m) (m) (m) (mm/h) -- 1000.0 498.00 .64 3.00E-02 3 2000.0 455.00 .35 3.00E-02 3 3000.0 452.00 .30 3.00E-02 3 4000.0 453.00 .35 3.00E-02 3 5000.0 446.00 .84 3.00E-02 3 ………. …….. …… . …….. …… ……… …….. ……. ………. ….. 30000.0 503.00 .57 3.00E-02 3 35000.0 501.00 .83 3.00E-02 3 40000.0 511.00 .59 3.00E-02 3 45000.0 696.00 .52 3.00E-02 3 50000.0 717.00 .50 3.00E-02 3 55000.0 662.00 .53 3.00E-02 3 60000.0 718.00 .54 3.00E-02 3 65000.0 732.00 .55 3.00E-02 3 Te132(78,2 hod) I132 (2,3 hod) Xe132 Varianty úniku: 1) pouze I132 = 3.92E+18 Bq 2) pouze Te132 = 1.48E+18 Bq XXIX. DNY RADIAČNÍ OCHRANY
faktor ochuzení (dcery) XXIX. DNY RADIAČNÍ OCHRANY
Příklad 2: ETE: Segmentovaný únik(podle STEP II b–case II) depozice aktivity na povrchu Kr88(2,8 hod)Rb88(17,8 min)Sr88 Hodinový únik: A) pouze Rb88 = 1.48E+18 Bq Hodinové meteopředpovědi: 28. června 2002, 00.00 UTM XXIX. DNY RADIAČNÍ OCHRANY
2. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] Počáteční únik: pouze 1.48E+18 Bq of Rb88) Situace: přesně po 2. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
3. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] Počáteční únik: pouze 1.48E+18 Bq of Rb88) Situace: přesně po 3. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
4. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] Počáteční únik: pouze 1.48E+18 Bq of Rb88) Situace: přesně po 4. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
5. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] Počáteční únik: pouze 1.48E+18 Bq of Rb88) Situace: přesně po 5. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
Příklad 2b: ETE: Segmentovaný únik(podle STEP II b–case II) depozice aktivity na povrchu Kr88(2,8 hod)Rb88(17,8 min)Sr88 Hodinový únik: B) pouze Kr88 = 1.48E+18 Bq Hodinové meteopředpovědi: 28. června 2002, 00.00 UTM XXIX. DNY RADIAČNÍ OCHRANY
2. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] (dcera Kr88) Počáteční únik: pouze 1.48E+18 Bq of Kr88) Situace: přesně po 2. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
3. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] (dcera Kr88) Počáteční únik: pouze 1.48E+18 Bq of Kr88) Situace: přesně po 3. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
4. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] (dcera Kr88) Počáteční únik: pouze 1.48E+18 Bq of Kr88) Situace: přesně po 4. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
5. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] (dcera Kr88) Počáteční únik: pouze 1.48E+18 Bq of Kr88) Situace: přesně po 5. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
6. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] (dcera Kr88) Počáteční únik: pouze 1.48E+18 Bq of Kr88) Situace: přesně po 6. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
7. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] (dcera Kr88) Počáteční únik: pouze 1.48E+18 Bq of Kr88) Situace: přesně po 7. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
8. hodina úniku [Bq/m2] Měrná depozice aktivity Rb88 [Bq/m2] (dcera Kr88) Počáteční únik: pouze 1.48E+18 Bq of Kr88) Situace: přesně po 8. hodině úniku XXIX. DNY RADIAČNÍ OCHRANY
HLAVNÍ OBLASTI POUŽITÍ • Podklady pro zavádění protiopatření v časné fázi • Identifikace možného mechanismu přenosu aktivity • do větších vzdáleností • Nutnost zohlednění při konstrukci operátoru pozorování • pro techniky asimilace modelu s měřeními XXIX. DNY RADIAČNÍ OCHRANY