Rapport om brandsimulering. Niras Safety - 2005 - Ejendomsviden

Projekt:
Mini-elevator
BILAG
B 1.3
Metodens orden
Metodens orden, i.e. målet for beregningsmodellens konvergenshastighed mod en
griduafhængig løsning fremgår af , og er bestemt iht.:
f
3
− f
ln
f
2
− f
p =
ln
( 2)
2
1
Benævnelserne f
i
svarer til løsningerne beregnet på de respektive gridkonfigurationer, se Tabel
15.
Konfiguration, a [ m]
Orden, p [ 1]
0,8 3,9
1,4 1,7
2,0 0,7
Tabel 17: Metodens orden.
FDS anvender finite difference schemes af anden orden, i.e. er den forventede teoretiske orden
p = 2 . Ovenstående afvigelser kan blandt andet tilskrives stærk ikke-linearitet i forbindelse
med løsningen af de strømningsmekaniske bevarelsesligninger samt anvendelsen af en
turbulensmodel. På det sted, hvor GCI-analysen er foretaget er turbulencen meget høj pga.
window plumens passage af vinduesoverliggeren. Sammen med det faktum, at der sker
forbrænding ud af vinduet, er løsningen på det evaluerede sted stærkt ikke-lineær og det må
derfor forventes, at den teoretiske orden ikke opnås.
B 1.3.1
GCI for de tre anvendte gridkonfigurationer
GCI
ij
imellem de anvendte gurationerne evalueres som:
(4)
GCI
f
i
− f
f
i
ij
= Fs
p
j
2 −1
Her er F = 1, 25 en sikkerhedsfaktor foreslået anbefalet i Roache, P. J. (1997).
s
(5)
GCI
ij
for de enkelte gridkonfigurationer fremgår af Tabel 18.
Konfiguration, a [ m]
GCI [%]
GCI [%]
Asymp. mål, A [ 1]
23
0,8 1,9 0,1 1,01
1,4 1,5 0,5 1,00
2,0 26,5 18,2 0,92
Tabel 18: GCI for de enkelte gridkonfigurationer.
12
Brandteknisk analyse af 01-10-05:
me_bta_udg_MBM_051001_1.doc
- 49 -
Safety