MODITIC - report on source term modelling

The European Defence Agency (EDA) project B-1097-ESM4-GP “MOdelling the DIspersion of Toxic Industrial Chemicals in urban environments” (MODITIC) has in the period from 2012 to 2016 studied the release and transport of neutral and non-neutral chemicals in complex urban environments, in order to enhance the understanding of the dominating physical processes involved, and to support improvements in modelling techniques. WP3000 concerns agent characterization and source modelling, this includes modelling the flashing, expansion, evaporation and rainout processes. Four studies have been conducted within WP3000 by four research institutes: • DGA tested Reynolds-averaged Navier-Stokes models for simulating the release of a pure gas or two-phase ammonia jet impinging an obstacle. The small time step required makes this method inappropriate for large geometrical domains. In addition, unphysical results for some parameters were obtained. • FFI used large eddy simulations approaches for the same two cases. The demands for the time step size and the grid resolution is even stricter in this case, precluding simulation of the entire experimental duration, as well as for large domains. Furthermore, in order to calculate the flashing of the liquid jet and the evaporation of the liquid droplets, a two-way coupling between the gas phase solver and the Lagrangian liquid droplet solver is needed, but the model used does not include such a two-way coupling at the present time. • FOI tested an approach where the conditions in the jet are calculated by an analytical approach and given as input to a subsequent computational fluid dynamics (CFD) simulation of the dispersion of gas and droplets. They validated the approach with field tests of the release and dispersion of liquified ammonia, and the results are encouraging. • Finally, INERIS used a large eddy simulations approach where the inlet conditions is specified from real (experimental) energy spectrum profiles. The vapour source term is specified according to an analytical calculation. The simulations are compared with field experiments of the release of ammonia, and the results are promising regarding the complexity for describing both near and far field dispersion. None of the institutes involved in MODITIC have a CFD capability that can fully model all the physical processes associated with a two-phase release of ammonia. The CFD methodologies produced somewhat unphysical results. In addition, the computer resources needed for a full CFD simulation makes such models inappropriate for large domains. On the other hand, a decoupled approach, where the rapid processes (flashing and expansion of the jet) are calculated with an empirical or analytical approach, and the slow processes (gas dispersion and air entrainment) are simulated with CFD, yields encouraging results. It is recommended to use the decoupled approach when modelling such demanding processes.

European Defence Agency-prosjektet B-1097-ESM4-GP MODITIC (MOdelling the Dispersion of Toxic Industrial Chemicals in urban environments) har i perioden 2012 – 2016 studert utslipp og spredning av nøytrale og ikke-nøytrale industrikjemikalier i urbane miljø. Målene er både å øke kunnskapen om de dominerende fysikalske prosessene som er involvert, og å støtte opp om forbedring av modelleringsteknikker. Arbeidspakke WP3000 dreier seg om beskrivelse av utslippskilden. Dette inkluderer beskrivelse av såkalt “flashing”, utvidelse av væsken og gassen etter utslippet, fordamping og væskedamdannelse. Fire studier er gjenomført innenfor WP3000 av fire forskningsinstitutt: • DGA brukte Reynolds-averaged Navier-Stokes-modeller for å studere utslipp av ammoniakk som gass og som tofase der utslippsstrømmen treffer en vegg. Tidsstegene som var nødvendige for simuleringene, var imidlertid så små at metoden er uegnet for store geometriske domener. I tillegg resulterte metoden i ufysikalske verdier på en del variabler. • FFI brukte large eddy-simuleringer for de samme to scenarioene. Disse metodene har imidlertid enda sterkere krav til tidssteget og størrelsen på beregningsgrid, slik at det ble uhåndterlig å simulere hele eksperimentet. Metodene blir også for tidkrevende for store geometriske domener. Modellen som ble testet, gir heller ikke mulighet til å simulere fullt ut de fysikalske prosessene som er involvert i flashing av væskestrømmen og fordamping av væskedråpene. • FOI testet en metodikk der forholdene til jetstrømmen blir beregnet analytisk og gitt som inngangsparametre til en computational fluid dynamics-beregning (CFD) av den videre spredningen av gass. Simuleringene er sammenlignet med feltforsøk av utslipp og spredning av ammoniakk, og resultatene er lovende. • INERIS har brukt en large eddy-simuleringsmetode der ingangsparametrene er spesifisert fra reelle (eksperimentelle) målinger av energispektrumsprofiler. Gasskilden er spesifisert fra en analytisk beregning. Simuleringene er sammenlignet med feltforsøk med utslipp av ammoniakk, og resultatene er lovende med tanke på å beskrive spredning både i nærfeltet og langt fra kilden. Ingen av instituttene har en CFD-kapabilitet som fullt ut kan beskrive alle fysikalske prosesser med et tofaseutslipp av ammoniakk. CFD-modellene som er brukt, ga til dels ufysikalske resultater, i tillegg til at de er for tidkrevende å brukes på store geometriske domener. I stedet anbefales det å benytte en totrinnsmodell, der de raske prosessene (flashing og utvidelse at strømmen) beregnes analytisk, og de langsomme prosessene (transport av gass og innblanding av luft) simuleres med CFD. En slik framgangsmåte har vist lovende resultater.