CFD MODELOVÁNÍ POŽÁRU V MÍSTNOSTI Autoři Doc. Dr. Ing. Milan JAHODA,VŠCHT Praha,
[email protected] Mgr. Jan ANGELIS, TÚPO, MV-GŘ HZS ČR,
[email protected] Ing. Otto DVOŘÁK, Ph.D., TÚPO, MV-GŘ HZS ČR,
[email protected]
Anotace Článek se zabývá predikcí teplotního a rychlostního pole v místnosti s otevřeným ohněm metodou počítačové dynamiky tekutin. Zdrojem ohně byl plynový hořák o výkonu 62,9 kW, který byl umístěný v místnosti s otevřeným dveřním otvorem. Experimenty byly provedeny Stecklerem a kol. a publikovány v literatuře. K simulačním výpočtům byl užit komerční řešič Fluent 6.3. Oheň byl modelován jednak jako zdroj tepla bez chemické reakce, jednak s probíhající jednoduchou reakcí spalování metanu. Dále byl testován vliv hustoty sítě, velikost časového kroku při nestacionární simulaci a vliv modelu turbulentního proudění.
Annotation The purpose of this study is to predict the flow and temperature fields due to a fire in small compartment using CFD simulations and code Fluent 6.3. The non-spreading fire was created by a centrally located 62.9 kW methane burner. The experiments were conducted by Steckler et al. and the results were presented in the literature. The fire has been simulated both a thermal source method and non-premixed combustion method. The Reynolds Averaging Navier-Stokes equations method together with standard k-epsilon turbulence model was compared with the Large Eddy Simulation technique. The domain and grid size dependence are documented. The results from CFD simulations were compared with experimentally obtained ones.
Úvod Matematické modelování požárů a jevů souvisejících s požáry patří mezi rychle se rozvíjející oblasti počítačové dynamiky tekutin. Je to především spojeno s rozvojem teorií spojených s hořením, numerických metod pro řešení soustav řídicích rovnic a vývojem počítačové techniky. V požárně bezpečnostním inženýrství můžeme pomocí fungujícího matematického modelu předpovídat rizika a důsledky požárů nebo odhadovat možné příčiny vzniku a šíření požárů při posuzování jejich možných hypotéz . Nedílnou součástí modelování je také validace získaných výsledků pomocí vhodného fyzikálního modelu. Matematické modely požárů lze rozdělovat do dvou základních typů: pravděpodobnostní a deterministické, přičemž každý typ užívá jiných soustav rovnic. Pravděpodobnostní (statistické) modely neřeší přímo fyzikální a chemické procesy spojené s požáry, ale vycházejí z historicky získaných dat o chování požárů. Výsledkem jsou jednoduché modely, které předpovídají jaký průběh požáru lze za daných podmínek očekávat. Deterministické modely využívají popisu fyzikálních a chemických jevů doprovázející požár pomocí matematických rovnic. Prakticky se užívají dva typy deterministických modelů: zónové modely a modely typu pole. Principem zónových modelů je rozdělení modelovaného prostoru na malý počet zón, pro které software řeší soustavu základních diferenciálních rovnic. Nejčastěji se užívá rozdělení prostoru na dvě zóny, tyto modely pak označujeme jako dvou-zónové. Při hoření se prostor rozdělí na horní teplou vrstvu, kterou tvoří horké spaliny z požáru, a dolní studenou vrstvu tvořenou přítomným/přicházejícím vzduchem. Do horní teplé vrstvy se počítá také ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 1
vzestupný proud spalin. Přehled užitých rovnic pro řešení zachování hmoty, energie a hybnosti uvádí např. Dudáček (2000). Pomocí zónových modelů dostaneme rychlou předpověď rozvoje požárů v uzavřených prostorech. Základní a velkou nevýhodou jsou omezené informace, neboť předpokládáme homogenní výstupní parametry (teplota, hustota kouře, koncentrace plynů) v každé zóně. Tuto nevýhodu odstraňují metody typu pole založené na počítačové dynamice tekutin (anglicky: Computational Fluid Dynamics, CFD). Prostřednictvím počítačové dynamiky tekutin získáme daleko přesnější výsledky. V současné době existují specializované programy pro simulace požárů, např. FDS (Fire Dynamics Simulator), SmartFire nebo FLACS (zaměřený navíc na modelování explozí). Vedle těchto produktů lze užít programy určené na pokročilé simulace proudění, které navíc umí řešit procesy spojené s chemickými reakcemi a sdílením tepla. Mezi tyto produkty patří také program Fluent od firmy Fluent Inc., který byl užit pro simulace v tomto příspěvku. Předkládaná práce se zabývá CFD predikcí teplotního a rychlostního pole v místnosti, které se vyvinuly v důsledku působení otevřeného ohně. Jde o počáteční simulace, které měly hlavní cíl odhalit problematiky v řešení a vymezit směr dalšího studia. Dlouhodobým cílem je užití CFD simulací průběhu požáru pro určení příčin vzniku požáru nebo pro požární prevenci.
Fyzikální model Experimentální hodnoty byly získány z požárních testů publikovaných Stecklerem a kol. (1982). Pro porovnání s daty ze CFD simulací bylo vybráno uspořádání prezentované na obr. 1. V místnosti o rozměrech 2,8x2,8x2,18 m byl centrálně umístěný hořák ve výšce 0,3 m. Místnost měla dveřní otvor o rozměrech 0,74x1,83 m otevřený do volného prostoru. Ve střední rovině dveří byly umístěny teplotní senzory a dvousměrové snímače rychlosti. Stěny a strop místnosti byly pokryty izolačními keramickými vlákny za účelem získání ustálených podmínek během 30 minut hoření. Palivem byl metan průmyslového složení s konstantním průtokem. Tepelný výkon hořáku o průměru 0,3 m byl 62,9 kW. Okolní teplota vzduchu byla 13 °C.
0,74 2,18 1,83
0,30 0,10 0,30 0,48
0,42
1,40
2,80
1,40 2,80
Obrázek 1: Geometrie experimentálního prostoru (rozměry v metrech) ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 2
Matematický model Vytvoření výpočetní sítě a následné řešení rovnic bylo děláno v prostředí programového balíku Fluent. Nejdříve byla definována geometrie sledované oblasti v preprocesoru Gambit. K modelové místnosti byla připojena vnější oblast pro modelování toku vzduchu/spalin vně místnosti a ve dveřním prostoru, obr. 2. DVEŘNÍ PROSTOR
pevné stěny volné stěny
vstup metanu
VNĚJŠÍ PROSTOR
MÍSTNOST
Obrázek 2: Ukázka definic hraničních oblastí a vnitřních prostor U modelů typu pole lze obecně simulace požárů resp. výpočet teplotních a rychlostních polí provádět dvěma různými způsoby: bez hoření (oheň je nahrazen zdrojem tepla), nebo s hořením, kdy je simulován reálný plamen včetně chemických reakcí. Simulace bez hoření dávají rychlé, ale pouze orientační hodnoty o polích teploty a rychlostí, zatímco u modelů s hořením získáváme navíc informace o koncentracích složek přítomných látek. Oba způsoby byly testovány. Dále byly testován vliv výpočetní sítě a vliv modelu turbulence. Jednotlivá nastavení modelových simulací je uvedeno níže. I. Vliv zdroje tepla (modelování ohně), kdy byl oheň buď nahrazen pouze zdrojem tepla (výpočet I.a, I.b), nebo byl oheň reprezentován spalováním metanu (výpočet I.c). Výpočet I.a: Vnější prostor před místností měl rozměry 1,4x2,8x2,8 m (délka x šířka x výška). Výpočetní síť byla tvořena šestistěny (kvádry) s celkovým počtem buněk cca 203 tisíc. Úloha byla počítána jako časově závislá s časovým krokem 0,5 s, s 20 iteracemi mezi každým časovým krokem a s celkovým počtem časových kroků 1800 což odpovídá době 15 minut reálného času, tedy době kdy by proudění mělo být celkově ustálené. K počítání turbulence byl vybrán k-epsilon standardní model. Oheň byl simulován zdrojem tepla o průměru 0,3 m a výšce 0,3 m, v němž se generovala energie 2966981 W/m3, což odpovídá výkonu 62,9 kW. Model P1 byl užit pro modelování radiace, jedná se o standardní model, který poskytuje reálné výsledky a přitom podstatně nezvyšuje výpočetní čas. Výpočet I.b: Velikost vnější oblasti byla oproti výpočtu I.a zvětšena na délku 4,2 m se zachováním původní šířky a výšky. Výpočetní síť byla tvořena mnohostěny s celkovým počtem buněk 248 tisíc. Jedná se o nový typ výpočetní sítě, která snižuje počet buněk (a také výpočetní čas) při zachování kvality výpočtu. Úloha byla počítána jako časově závislá s časovým krokem stejným jako při výpočtu I.a. K simulaci turbulentního proudění byl ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 3
vybrán k-epsilon RNG model. Oheň byl nahrazen zdrojem tepla o stejném měrném výkonu jako ve výpočtu I.a. Výpočet I.c: Geometrie vnější oblasti, výpočetní síť a nastavení výpočtu včetně časového kroku byly shodné s výpočtem I.a. Oheň byl reprezentován spalováním metanu s jednou chemickou reakcí (2O2 + CH4 → CO2+2H2O), kdy z válcového hořáku vstupoval metan do místnosti byl s konstantním hmotnostním průtokem 0,00113 kg/s, což přibližně odpovídá výkonu 62,9 kW. K modelování toku složek (hoření) byl užit model „non-premixed combustion“, kdy palivo (např. metan) a oxidant (např. vzdušný kyslík) přicházejí do prostoru hoření (reakční oblasti) odděleně. Předpokládá se, že směšovací doba složek je daleko delší než doba reakční. Dále se předpokládají shodné difúzní koeficienty reakčních složek. Termochemický stav tekutiny je zjednodušený do jedné skalární veličiny: směsného zlomku (anglicky: Mixture Fraction). Směsný zlomek je podíl hmotnosti paliva (např. metan) k sumě hmotností paliva a oxidačního činidla (např. vzdušný kyslík). Transportní rovnice jsou zjednodušeny na transportní rovnici směsi bez zdrojového členu, hoření je tedy zjednodušeno pouze na problematiku směšování. Výhodou směsného zlomku je, že další skalární veličiny, které jsou funkcí zlomku (např. hustota nebo hmotnostní zlomek složky) může být vypočítán z lokální hodnoty směsného zlomku a jeho čtverci směrodatné odchylky. Termochemické výpočty jsou předpřipravené a tabelované ve Fluentu. Vztah mezi turbulencí a chemií je vyjádřen předpokládaným tvarem pravděpodobnostně hustotní funkcí (anglicky: Probability Density Function, PDF). Metoda je výpočetně nenáročná, není třeba řešit velkou soustavu transportních rovnic složek. Prostorová a časová diskretizace rovnic byla protiproudá, druhého řádu. Turbulence byla modelována pomocí standardního k-epsilon modelu rozšířeného o vztlakový člen. V tomto typu výpočtu byly definovány čtyři objemové zóny (viz obr. 2): místnost, dveřní prostor, vnější prostor a vrstva tekutiny o tloušťce 2 cm a průměru 30 cm nad vstupem plynného paliva do místnosti. Tato vrstva tekutiny sloužila k zážehu paliva. Horní rovina „zážehové“ vrstvy byla na shodné úrovni jako horní stěna hořáku. Okrajové podmínky. Reálné stěny byly nastaveny jako pevné stěny se standardní stěnovou funkcí. Stěny přidružené oblasti byly definovány jako tlakový vstup (anglicky: pressure inlet), strop přidružené oblasti jako tlakový výstup (anglicky: pressure outlet). Počáteční podmínky. Nulové složky rychlosti, nulový přetlak vzduchu, objemové zóny neobsahují metan, teplota zážehové vrstvy 1 000 K a teplota ostatních zón 286 K. II. Vliv hustoty a typu výpočetní sítě, byl testován na příkladu spalování metanu. Ve všech případech byla velikost vnějšího prostoru totožná s velikostní místnosti. Ostatní nastavení byly shodné s výpočtem I.c, časový krok při iteracích byl 1 s a dosažený reálný čas byl 200 s. Výpočet II.a: Výpočetní síť byla tvořena čtyřstěny (nestrukturovaný typ sítě) a obsahovala 757 tisíc buněk. Výpočet II.b: Výpočetní síť byla tvořena opět čtyřstěny s obsahem cca polovičního počtu buněk, tj. 314 tisíc buněk. Výpočet II.c: Výpočetní síť pro výpočet II.c vznikla ze sítě výpočtu II.b převedením na vícestěny (polyhedra). Na vázali jsme tak na výpočty I.a a I.b, kdy zdrojem tepla nebyla chemická reakce. Počet buněk ve výpočetní síti klesl na 62 tisíc. Souhrnné informace o počtech buněk v jednotlivých oblastech jsou uvedeny v tabulce 1. Ukázky jednotlivých typů sítí jsou znázorněny na obr. 3 a 4.
ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 4
Tabulka 1: Počty buněk výpočetní sítě ve výpočtu II Síť Oblast II.a (757k) II.b (314k) II.c (62k) Místnost 395 501 154 403 31 565 Dveřní prostor 6 061 6 061 1 392 Vnější prostor 355 026 155 403 29 235 Zážehová vrstva 345 345 145 756 933 314 018 62 337 Celkem buněk
Obrázek 3: Ukázka výpočetní sítě tvořené z čtyřstěnů (výpočet II.b)
Obrázek 4: Ukázka výpočetní sítě tvořené z mnohostěnů (výpočet II.c) III. Vliv modelu turbulence a hmotnostního toku metanu, byl opět testován na příkladu spalování metanu. Vnější prostor byl zvětšen na 1,5 násobek délky místnosti a 1,33 násobek výšky při zachování shodné šířky s místností. Výpočetní prostor byl nejprve rozdělen nestrukturovanou sítí tvořenou čtyřstěny sestávající z 460 tisíc buněk v preprocesoru Gambit. Tato síť byla poté načtena do řešiče Fluent, kde byly čtyřstěny převedeny na mnohostěny. Počet buněk sítě tak klesl na 90 tisíc, viz. obr. 5. ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 5
Obrázek 5: Výpočetní síť pro simulace č. III, 90 tisíc buněk (mnohostěny) V části III bylo provedeno celkem sedm typů výpočtů. Šest bylo časově závislých s časovým krokem 0,2 s. Časově závislé výpočty se od sebe lišily hmotnostním průtokem metanu a použitým modelem turbulence. Steckler a kol. (1982) ve své zprávě uvádějí pouze výkon hořáku a informaci, že použili metan o průmyslovém složení. Přesné procento obsahu metanu v palivu však uvedeno není. Za předpokladu, že se jednalo o zemní plyn, tak hmotnostní % metanu se mohou pohybovat asi v rozmezí 83 – 98,5. Z tohoto důvodu se spalné teplo zemního plynu mění dle složení. Po přepočtu na hmotnostní průtok byly následně testovány tři hodnoty: 0,00113, 0,00126 a 0,00136 kg/s. Vedle standardního modelu k-epsilon modelu turbulence a časově zprůměrněné NavierStokesovy rovnice (anglicky Reynolds-Averaged-Navier-Stokes, RANS) byly k výpočtu použity filtrované Navier-Stokesovy rovnice a metoda velkých vírů (anglicky Large Eddy Simulation, LES). Při kterých jsou velké a střední víry počítány přímo a víry menší než je velikost buněk výpočetní sítě jsou modelovány pomocí podsíťových (anglicky Sub-grid Scale, SGS) modelů, např. dynamického Smagorinského-Lillyho modelu, který byl užit při našich výpočtech. Tabulka 2 souhrnně znázorňuje všechny užité kombinace ve výpočtu č. III Tabulka 2: Přehled provedených simulací ve výpočtu č. III Výpočet Časová Průtok CH4 Model -1 závislost [kg s ] turbulence III.a Ano 0,00113 LES III.b Ano 0,00113 k-epsilon III.c Ano 0,00126 LES III.d Ano 0,00126 k-epsilon III.e Ano 0,00136 LES III.f Ano 0,00136 k-epsilon III.g Ne 0,00136 k-epsilon
Reálný čas simulace (s) 130 100 170 150 315 73 -
Reálný čas simulace odpovídá skutečné době, při kterém již bylo dosaženo pseudoustálenému stavu polí teplot/rychlostí. Pro kontrolu byl výpočet III.e následně prodloužen až ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 6
na 30 min reálného času, tj. na dobu shodnou s dobou trvání fyzikálního modelu, přičemž získané výsledky je již prakticky neměnily. Okrajové a počáteční podmínky byly jinak shodné s nastavením ve výpočtech č. II. včetně nastavení řešiče.
Výsledky a diskuse Na obrázku 6 je uvedena ukázka rozvinutého teplotního a rychlostního pole v případě metody LES s chemickou reakcí. U všech ostatních uspořádání jsou výsledné profily podobné. Na obrázku jsou patrné teplotní vrstvy – horní horká a spodní "studená", přičemž stoupající teplotní proud plamene od hořáku se počítá do horní horké vrstvy.
Obrázek 6: Ukázka rozvinutého teplotního a rychlostního (x-ová složka rychlosti) pole Pro lepší srovnání vypočtených hodnot s experimentem byly porovnány výsledky teplot a rychlostí ve střední linii dveřního prostoru. Na obrázku 7 je znázorněno porovnání teplotních profilů v linii dveřního prostoru. Je zde patrná dobrá shoda s experimentem pro výpočet I.a, ve kterém byl oheň modelován jako zdroj tepla (bez chemické reakce). Také výpočet I.b probíhal bez chemické reakce, patrná odlišnost v horní teplé zóně bude pravděpodobně více způsobený použitým modelem turbulence (RNG k-epsilon), než změnou výpočetní sítě na mnohostěny, která při jinak stejných podmínek dává shodné výsledky při zachování hustoty sítě. Jakmile se však použije model s chemickou reakcí (I.c), dochází ke zhoršení shody s experimentem. To je dáno složitějším výpočtem a vstupem dalších proměnných ovlivňující predikci teplotního a rychlostního pole. Proto jsme se dále zaměřili na systémy s chemickou reakcí.
ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 7
2 1,8
Experiment I.a I.b I.c
h (m)
1,6 1,4 1,2 1 0,8 0,6 0,4 0,2 0 0
50
100
150
200
250
t (°C)
Obrázek 7: Porovnání teplotních profilů pro výpočet I (zliv tepelného zdroje) Ukázka vlivu hustoty výpočetní sítě (výpočet II) je uvedena na obrázku 8. Ve všech třech případech došlo k velké odchylce vypočtených teplotních polí od hodnot z fyzikálního modelu. Porovnáme-li změny teplotních profilů s hustotou sítě, dostaneme, že zvýšením počtu buněk z 315 tisíc na víc než dvojnásobek, tj. 757 tisíc, nedochází v daném nastavení řešení ke zlepšení. S porovnáním s výsledky výpočtu I.c lze usuzovat, že pravděpodobnou příčinou velké odchylky bude dlouhý časový krok při řešení - jedna sekunda. Snížením jeho délky dostaneme zlepšení řešení, viz. obr. 9. Při řešení výpočtu III byl časový krok 0,2 s. Na obrázku je také znázorněn vliv množství přiváděného paliva. Přestože došlo ke zlepšení shody teplot v horní části dveřního prostoru, kde odcházejí produkty hoření do vnějšího prostoru, nebyla tak dobrá jako ve výpočtu I.a. Proto byl testován jiný model turbulence – metoda velkých vírů (LES). Tato metoda je úspěšně používána ve speciálním řešiči FDS pro simulace požárů. 2 1,8
h (m)
1,6 1,4 1,2 1 0,8 0,6
Experiment II.a: 757k II.b: 314k II.c: 62k
0,4 0,2 0 0
50
100
150
200
250
300
350
400
450
t (°C)
Obrázek 8: Porovnání teplotních profilů pro výpočet II (vliv hustoty výpočetní sítě)
ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 8
2 1,8
h (m)
1,6 1,4 1,2 1 0,8
Experiment III.b: SKE; 0,00113 kg/s III.d: SKE; 0,00126 kg/s III.f: SKE; 0,00136 kg/s III.g: SKE-stac; 0,00136 kg/s
0,6 0,4 0,2 0 0
50
100
150
200
250
t (°C)
Obrázek 9: Porovnání teplotních profilů pro výpočet III (vliv hm. průtoku metanu), metoda SKE 2 1,8
h (m)
1,6 1,4 1,2 1 0,8 0,6
Experiment III.a: LES; 0,00113 kg/s III.c: LES; 0,00126 kg/s III.e: LES; 0,00136 kg/s
0,4 0,2 0 0
20
40
60
80
100
120
140
160
180
200
t (°C)
Obrázek 10: Porovnání teplotních profilů pro výpočet III (vliv hm. průtoku metanu), metoda LES Na obrázku 10 jsou znázorněny výsledné hodnoty teplotních profilů v oblasti dveří při použití metody LES. Shoda experimentálními výsledky je velmi dobrá. Vyšší objemové průtoky dávají vyšší hodnoty teploty. Díky nestacionární simulaci velkých vírů může docházet k časově-prostorové fluktuaci výstupního proudu horkých spalin, které mohou ovlivňovat teplotu vstupního proudu vzduchu do oblasti dveří. Je zajímavé, že ačkoliv výsledné teplotní profily v oblasti dveří se mnohdy značně odlišují, získané rychlostní profily jsou v dobré shodě, jak pro výpočty bez chemické reakce, obr. 11, tak s chemickou reakcí pro různé modely turbulentního proudění, obr. 12. Odchylky patrné na obr. 12 ve spodní části dveřního prostoru jsou způsobené časovou proměnlivostí fluktuace ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 9
rychlosti díky metodě velkých vírů, neboť uvedené profily jsou okamžité hodnoty rychlostí ve stejný časový okamžik. Na základě tohoto poznatku je třeba se dále zaměřit na modelování spalování. 2 Experiment Výpočet I.a Výpočet I.b Výpočet I.c
1,8 1,6 1,4 h (m)
1,2 1 0,8 0,6 0,4 0,2 0 2,5
2
1,5
1
0,5
0
-0,5
-1
-1,5
v (m/s)
Obrázek 11: Výškové profily x-ové složky rychlosti v linii dveří pro výpočet I
2 1,8 1,6
h (m)
1,4 1,2 1
Experiment III.a: LES; 0,00113 kg/s III.b: SKE; 0,00113 kg/s III.c: LES; 0,00126 kg/s III.d: SKE; 0,00126 kg/s III.e: LES; 0,00136 kg/s III.f: SKE; 0,00136 kg/s III.g: SKE-stac; 0,00136 kg/s
0,8 0,6 0,4 0,2 0
2,5
2
1,5
1
0,5
0
-0,5
-1
-1,5
v (m/s)
Obrázek 12: Výškové profily x-ové složky rychlosti v linii dveří pro výpočet III
ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 10
Závěry Základním cílem této studie bylo získání zkušeností s modelováním požárů v místnostech pomocí programu Fluent. Dále definovat problematické části při simulacích a vytipovat další směry výzkumu. Pro výpočet teplotního pole, kdy není třeba zjišťovat koncentrace látek spojených s procesem hoření, je výhodnější nahradit otevřený oheň pouze zdrojem tepla bez chemické reakce. Toto vede k uspokojivým výsledků a ke snížení výpočetní náročnosti. Ze získaných zkušeností lze nastínit následující možné směry pro další studium: Pro vybranou výpočetní síť testovat vliv relaxačních faktorů pro úspěšnou konvergenci výpočtů v závislosti na časovém kroku. Vyzkoušet jiné modely pro radiaci. Testovat další modely pro transport složek s chemickou reakcí (hoření). Na základě vhodného fyzikálního modelu vypočítat koncentrace složek spalin v závislosti na čase a prostoru.
Poděkování Tato práce byla finančně podporována MV-GŘ HZS ČR, Technickým ústavem PO v rámci řešení jeho dílčího výzkumného úkolu z výzkumného projektu Programu bezpečnostního výzkumu č. VD20062010A07 a Ministerstvem školství, mládeže a tělovýchovy ČR projektem č. MSM 6046137306.
Literatura [1] Dudáček, A., 2000. Úvod do software pro matematické modelování požárů, elektronická učebnice pro VŠB-TU Ostrava [on-line], naposledny dostupné z http://homel.vsb.cz/~dud20/SMP/T1/soft_mmp.htm [2] Steckler, K. D., Quintiere, J. G., Rinkinen, W. J., 1982. Flow induced by fire in a compartment, NBSIR 82-2520, National Bureau of Standards, Center for Fire Research, Washington, USA. [3] Fluent 6.3 User’s Guide, Fluent Inc., Libanon, USA, 2006 [4] Dvořák, O. a kol. Vývoj a validace požárních modelů pro stanovení vývinu/šíření tepla a kouře, toxických plynů, tlakových vln pro simulaci/interpretaci scénářů požárů/výbuchů a jejich ničivých účinků. Závěrečná výzkumná zpráva výzkumného projektu č. VD20062007A01. Praha: MV-GŘ HZS ČR-TÚPO, 2008. [5] Angelis, J., Bursíková, P., Dvořák,O. Praktické využití matematických modelů Smartfire a Fluent k rekonstrukci scénářů požáru pro potřebu PTE. Sborník přednášek XVI. ročníku mezinárodní konference Požární ochrana 2007,.Ostrava: VŠB-TUO, 2007, 17-25.
ANSYS konference 2008 16. ANSYS FEM Users‘ Meeting & 14. ANSYS CFD Users’ Meeting Luhačovice 5. - 7. listopadu 2008 11