Particle depletion in a differentially heated cavity

Particle depletion in a differentially heated
cavity
Jarmo Kalilainen
SAFIR2014 FINAL SEMINAR, 20.3.2015
Introduction
• Aerosol behaviour in a containment building during an accident
scenario has an importance on the mitigation of the FP release to the
environment.
• Large scale SA experiments (e.g. Phebus FPT2-3) observed relatively
large deposition to vertical containment walls
Possible effects of turbulent natural convective flow to the particle wall
deposition?
• The aim of this study is to study fluid and heat transfer and particle
retention in a Differentially Heated Cavity (DHC). The experimental
results are used in validation of a large eddy simulation (LES) of the
DHC and Lagrangian particle tracking (LPT) simulations.
DIANA Facility
•
DIfferentially heated cavity with Aerosol in turbulent NAtural
convection (DIANA) facility.
•
The facility has two vertical isothermal aluminium walls and
four adiabatic glass walls.
•
The facility must allow optical access for laser-based
measurement devices used to determine the flow properties as
well as particle deposition rates
•
In order for the Boussinesq approximation, used in the
simulation work to be valid the heat difference between the
walls must remain below 50 K.
•
Fluid: air, T = 57 °C - 18 ° C = 39 ° C, Rayleigh number Ra ~ 109.
Turbulent flow.
PIV
Cavity
Two part study
• Experimental study:
Measurement of flow field and gas temperature in the
DIANA cavity
Measurement of particle deposition rates using
monodisperse SiO2 particles with diameters 1 µm and 2.5
µm.
• Computational study:
Validation of the LES model using experimental temp.
boundary conditions (BC).
Comparison of measurement to Lagrangian particle tracking
data obtained using the validated LES model.
Particle image velocimetry (PIV)
Mean and fluctuating velocity x-y components of the
flow field measured near the cavity lateral center
plane using a PIV.
•Imaging measurement of velocity fields of particles in a transparent fluid
•Flow visualization
•Double pulse laser and double frame CCD-camera
•Analysis of flow fields from two images taken at short interval
Fluid velocity
A series of 120 PIV measurements at the lower
left corner, 330 mm from the front adiabatic wall.
An average velocity, calculated from 1200 PIV
images at the lower left corner, 330 mm from the
front adiabatic wall.
Fluid velocity
Vmag=
+
=
•
Mean velocity magnitudes and the turbulence intensity at the cavity center plane.
•
Turbulent flow encircling a stagnant core next to the isothermal and horizontal walls.
LES of the cavity
•
Two LES with different BCs at the
horizontal wall were used in the particle
tracking simulations.
•
In the first LES (LES-WT), the thermal BCs
for isothermal and horizontal walls were
obtained from the wall temperature
measurement data.
•
The second simulations using idealized
adiabatic BCs on horizontal walls (LESIDEAL), was validated against the DNS
data by Puragliesi (PhD thesis, EPFL, 2010) with
similar conditions.
•
The LES-WT validated using the
measurement data.
•
Flow and temp. field differed with LESIDEAL – Flow geometry different, less
turbulent.
Horizontal and vertical temperature profiles measured with K-type
thermocouple, compared against LES-WT and LES-IDEAL simulation data.
Both show stratified temperature distribution at cavity core.
LES of the cavity
Comparison between PIV measurement of mean vertical and horizontal
velocity components and LES-WT simulation.
Rms of velocity along horizontal and vertical profiles from PIV
measurements and LES-WT simulations .
LES of the cavity
Velocity magnitude
LES-IDEAL
PIV
Turbulence kinetic
energy
LES-IDEAL
LES-WT
Particle deposition measurements
•
Monodisperse SiO2 (1 and 2.5 µm)
particles seeded from the bottom of
the cavity.
•
The change of particle concentration
was investigated by introducing a
laser sheet to the cavity through the
top or front glass wall and measuring
the intensity of the reflected light
from the particles by a CCD camera.
•
Tapered Element Oscillating
Microbalance (TEOM) was used in
addition to laser intensity
measurements.
Particle deposition measurements
Experiment matrix.
•
Exp.
Particle
diam.
[µm]
Laser
window
placem.
Exp.
Particle
diam.
[µm]
Laser
window
placem.
1
1.0
xz hot
6
2.5
xz hot
2
1.0
xy hot
7
2.5
xy hot
3
1.0
xy
centre
8
2.5
xy
centre
4
1.0
xy cold
9
2.5
xy cold
5
1.0 (dry
atm.)
xz hot
Particle distribution in the cavity direction
(experiments 1 & 6).
•
Particles approx. uniformly distributed to
the cavity at the lateral direction.
•
The results indicated uniform deposition
rates throughout the cavity atmosphere.
•
No effect of dry atmosphere on particle
depletion rate was observed in exp. 5.
Particle tracking simulations
Simulation matrix.
1 µm particles
LES-IDEAL-CRW tp ON
LES-IDEAL-CRW tp OFF
LES-WT-CRW tp ON
LES-WT-CRW tp OFF
LES-WT tp ON
2.5 µm particle
LES-IDEAL-CRW tp ON
LES-IDEAL-CRW tp OFF
LES-WT-CRW tp ON
LES-WT-CRW tp OFF
LES-WT tp OFF
number
of particle
thermophoretic
force considered
computation
ended [s]
10000
10000
10000
10000
10000
Yes
No
Yes
No
Yes
~6300
~6300
~6300
~6300
~1000
10000
10000
10000
10000
100000
Yes
No
Yes
No
No
~2300
~2300
~2300
~2300
~1700
•
In Continuous Random Walk (CRW) simulations the fluctuating
fluid velocity is modelled using a Markov chain based on the
normalized Langevin equation which takes into account the
inhomogeneities of the turbulence (Dehbi, 2008). The mean flow and
temperature fields, along with the average Reynolds stresses were
extracted from LES and used to calculate the fluctuating velocity
component.
•
CRW compared against pure LES particle tracking. Good match.
•
Particle tracking simulations made
using the LES-WT with realistic BCs
and LES-IDEAL with ideal adiabatic
BCs at horizontal walls.
•
Thermophoretic force (tp)
considered in some simulations.
Particle tracking simulations
•
In the stirred settling case, particles are kept uniformly
distributed in a cubic volume with side length L, and are
deposited only through gravitational settling to the
enclosure floor (Hinds, 1999).
•
The inclusion of thermophoretic force to the simulation
has only relatively small effect on the deposition speed.
Only effects on distribution of deposit on different
surfaces.
Average decay constants from the laser intensity and mass
concentration measurements, particle tracking simulations and
form the stirred settling calculation.
Comparison of particle decay from the cavity
atmosphere with dp = 1 µm particles.
dp = 1 µm
dp = 2.5 µm
Average , TEOM [s]
5220 ± 190 s
1100 ± 90 s
Average , laser intensity [s]
4970 ± 60 s
1800 ± 80 s
LES-WT-CRW tp ON
4890 s
1510 s
LES-WT-CRW tp OFF
5060 s
1520 s
LES-WT tp (1 µm ON) (2.5 µm OFF)
5530 s
1610 s
LES-IDEAL-CRW tp ON
6360 s
1860 s
LES-IDEAL-CRW tp OFF
6920 s
1840 s
Theoretical Stirred Settling (Hinds, 1999)
10210 s
1780 s
Summary
•
Flow and temperature measurements were performed in DIANA cavity with
turbulent natural convective flow.
LES model using DIANA cavity BCs produced valid representation of the flow and temp.
fields.
•
Aerosol depletion in DIANA cavity was investigated experimentally using
monodisperse SiO2 particles with dp = 1 µm and dp = 2.5 µm.
Depletion rates from Lagrangian particle tracking with realistic BCs at the cavity walls
agreed well with the measurement data.
•
Comparison to theoretical stirred settling indicated large discrepancy with 1 µm
particle.
LPT with ideal boundary conditions resulted on slower particle depletion due to altered
flow geometry and decreased level of turbulence in the cavity.
Stirred settling not applicable in DHC geometry even in the cavity core where the flow is
almost stagnant.