Turbulent strømning og fordampning av væske i sjakt

Abstract

Formålet med denne studien er å bidra til bedre forståelse og grunnleggende kunnskap om hvordan kjemiske trusselstoffer oppfører seg ved utslipp og hvordan de spres i luft. Slik kunnskap er et avgjørende grunnlag for konsekvensvurderinger av mulige utslipp og for planlegging av beskyttelsesog mottiltak, og for beredskapstrening og øvelser for Forsvaret og for sivile beredskapsmyndigheter. Det er utviklet en rekke fareprediksjonsverktøy og spredningsmodeller for å vurdere konsekvenser av angrep, sabotasje eller ulykker. En komplett modell inneholder en matematisk og fysisk beskrivelser av en rekke ulike fenomener som til sammen utgjør en spredningshendelse, og kan grovt deles inn i hhv kildemodellering, transportmodellering, og effektmodellering. Forbedring av denne type modeller er et felt der mye forskning gjenstår, da det per i dag er til dels stor usikkerhet knyttet til de tre leddene i modellkjeden. Denne studien fokuserer på kildemodellering. FFI har bidratt med et scenario for innendørs spredning av sarin i NATO-studien RTO-SAS-061 ’Defence against CBRN-attacks in the changing NATO strategic environment’. Scenarioet ble også simulert av andre deltakere, og det viste seg at de ulike deltakerne brukte en fordampningsrate som varierte med inntil en faktor 10. Dette bidrar til tilsvarende store variasjoner for beregninger av konsekvenser og skadeomfang, og er bakgrunnen for dette arbeidet. En nøyaktig beskrivelse av fordampningen er avgjørende for å modellere spredningen realistisk, både når det gjelder tidsforløp og konsentrasjon. Egenskaper til luftstrømmen over den avdampende væskeflaten har stor innvirkning på fordampningsraten. Spesielt viktig er det turbulente grensesjiktet rett over dammen. Hensikten med denne studien er å studere fordampning av en væske, i dette tilfellet nervegassen sarin, i en sjakt, for å bidra til en bedre beskrivelse av fordampning av væsker generelt. Dette kan brukes både i kombinasjon med enkle fareprediksjonsverktøy, enkle spredningsberegninger og mer komplekse modeller. I dette arbeidet er turbulent strømning i en kvadratisk sjakt simulert ved hjelp av ulike turbulensmodeller i ’Computational Fluid Dynamics’-koden Fluent. Resultatene er sammenlignet med direkte numeriske simuleringer for å verifisere den valgte metoden. Basert på de turbulente hastighetsfeltene fra simuleringene er fordampning fra en dam i sjakten, og videre spredning av dampen gjennom sjakten, simulert. Dette er spesielt relevant for innendørs spredning der ventilasjonsanlegg blir brukt, men også viktig i andre scenarier der fordampning av kjemiske trusselstoffer fra en dam inngår. Resultatene viser at luftstrømmen påvirker væskeflaten, og dette har stor innvirkning på fordampningsraten. Fordampningsraten beregnet med simuleringene er omtrent fire ganger høyere enn fordampningsraten beregnet med en enkel modell som baserer seg på molekulær diffusjon i en laminær strøm. Resultatene viser også at en del fordampet sarin fester seg på veggene av sjakten, og dette vil gi opphav til sekundær fordampning. Dette arbeidet vil bli videreført ved å se på andre trusselstoffer og variere egenskaper ved strømningen. Planen er å lage en database med startbetingelser for spredningsberegninger, for senere å kunne gjennomføre fordampningsberegninger raskt og nøyaktig

The aim of this study is to contribute to a better understanding and fundamental knowledge on the behaviour of harmful chemical substances when released, and how they are dispersed in air. Such knowledge is a crucial foundation for estimating the consequences of possible chemical incidents, and for the planning of protective actions and countermeasures, and also for training and exercises for the Defence and for civilian emergency authorities. There are numerous models for calculating the dispersion and estimating hazard areas, and possible consequences of attacks, sabotage or accidents. A complete model includes a mathematical and physical description of a number of phenomena that constitute a dispersion incident, and can roughly be divided in source modelling, transport modelling, and effect modelling respectively. Improvements of such models is a field that still requires a lot of research, as there are large uncertainties in the three parts of the modelling chain. This study focuses on source modelling. FFI has contributed to the NATO study RTO-SAS-061 ’Defence against CBRN-attacks in the changing NATO strategic environment’ with a scenario of indoor dispersion of sarin. The scenario was also simulated by other participants, and it turned out that different participants used evaporation rates that were set apart with a factor of up to 10. This gives correspondingly large variations on the calculations of the consequences and the extent of damages, and is the background for this study. An accurate description of the evaporation is vital for a realistic modelling of the dispersion, both for time and concentration. The properties of the air flow above the evaporating liquid surface have a great influence on the evaporation rate. The turbulent boundary layer directly above the pool is particularly important. The object of this work is to study the evaporation from a liquid, in this case the nerve agent sarin, in a shaft, in order to contribute to a better description of the evaporation of liquids in general. This can be used both in combination with simple tools for estimating the hazard areas, simple dispersion models, and more complex models. In this work, the turbulent flow through a quadratic shaft has been simulated by means of various turbulence models in the Computational Fluid Dynamics code Fluent. The results are compared with direct numerical simulations in order to verify the chosen method. Based on the turbulent velocity fields from the simulations, the evaporation from a pool in the shaft and the subsequent dispersion of the vapour are simulated. This is especially important for indoor dispersion where ventilation systems are utilised, but also important in other scenarios that include evaporation of chemicals from a pool. The results show that the air flow perturbs the liquid surface, and this has a great effect on the evaporation rate. The evaporation rate calculated from the simulations is about four times greater than the corresponding evaporation rate calculated with a simple formula based on molecular diffusion through a laminar flow. They also show that some of the sarin binds to the walls of the shaft, and this causes a secondary evaporation. This work will be continued by examining other chemical substances and by varying the properties of the flow. The plan is to construct a data base with initial conditions for dispersion calculations, so that later calculations of evaporation can be executed quickly and accurate.