Introduction Numerical methods Title Parameters for computation of the short-scale scenario Introduction Results Short scale scenario Vladim´ ır Fuka, Josef Brechler Dept. of Meteorology and Environment Protection, Charles University in Prague Sevilla, 8.–10. June 2010 Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Numerical methods Title Parameters for computation of the short-scale scenario Introduction Results Introduction Our model is general-purpose high resolution code for Computational Fluid Dynamics. We have very limited experience in modelling of radioactive substances. The model uses a Large Eddy Simulation framework, which requires high resolution and has large demands on computational time. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Large Eddy Simulation Numerical methods Subgrid models Parameters for computation of the short-scale scenario Other numerical methods Results Large Eddy Simulation Filtered momentum equations ∂ t + ∂ u i u j + ∂ 2 ν S ij + ∂τ ij ∂ u i = − ∂ p (1) ∂ x j ∂ x i ∂ x j ∂ x j τ ij = u i u j − u i u j (2) Continuity equation ∂ u i = 0 (3) ∂ x i This filtering is performed implicitly in our model. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Large Eddy Simulation Numerical methods Subgrid models Parameters for computation of the short-scale scenario Other numerical methods Results Smagorinsky model Turbulent viscosity τ ij = 2 ν t S ij + 1 3 τ ll δ ij (4) Smagorinsky model ν t = ( C S ¯ ∆ ) 2 | S | (5) � | S | = S ij S ij (6) Too diffusive for most flows. Theoretical value C S ≈ 0 . 2 is too high in many cases. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Large Eddy Simulation Numerical methods Subgrid models Parameters for computation of the short-scale scenario Other numerical methods Results Vreman model � B β ν t = c (7) α ij α ij where (8) α ij = ∂ u j (9) ∂ x i � ∆ 2 β ij = m α mi α mj (10) m B β = β 11 β 22 − β 2 12 + β 22 β 33 − β 2 23 + β 33 β 11 − β 2 (11) 31 Gives zero ν t for laminar flows with large strain rates unlike the Smagorinsky model. The constant c is also not a universal constant, but it’s optimal value changes for different problems less than C S . Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Large Eddy Simulation Numerical methods Subgrid models Parameters for computation of the short-scale scenario Other numerical methods Results Other numerical methods Projection (fractional step) method with second order Runge-Kutta method for time discretization. Finite volume spatial discretization on a staggered grid. Second order central discretization for advective and diffusive fluxes. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Large Eddy Simulation Numerical methods Subgrid models Parameters for computation of the short-scale scenario Other numerical methods Results Wall model Wall model Because of poor wall resolution is expected in atmospheric applications and we use aerodynamically rough walls a wall model must be used. ν + ν t = ( u τ ) 2 y / u (12) y + < y + � y + if u + = 0 (13) 0 , 41 ln y + + 5 , 2 y + > y + 1 if 0 or (14) 1 u + = 0 , 41 ln ( y / z 0 ) (15) Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Large Eddy Simulation Numerical methods Subgrid models Parameters for computation of the short-scale scenario Other numerical methods Results Dry deposition Dry deposition For dry deposition of particles we employ a semi-empirical scheme by Kharchenko (1997). v d = f ( ν, d , D , z 0 , u ∗ , τ + ) (16) where (17) d = particle diameter (18) D = coefficient of Brownian diffusion (19) u 2 τ + = u sed ∗ (20) g ν Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Large Eddy Simulation Numerical methods Subgrid models Parameters for computation of the short-scale scenario Other numerical methods Results Turbulent inflow On of difficult problems in LES is imposing inflow boundary conditions, that are turbulent and as realistic, as possible. One possibility is a separate computations of the incoming flow (e.g. a periodic boundary layer or channel), but it is very resource demanding. The other possibility is generation of the incoming field using random numbers. The generated flow has to have propper length scale and spectrum to contain the lower wavenumbers. White noise would be damped very quickly and turns out to be useless. Also correct correlations between velocity components are important. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Large Eddy Simulation Numerical methods Subgrid models Parameters for computation of the short-scale scenario Other numerical methods Results Turbulent inflow We employ a method proposed by Klein, Sadiki and Janicka (2003) in a more effective form by Xie and Castro (2008). First a random field U i is generated with correct point to point and time statistics using filtering. The filter is constructed, so that the resulting field has prescribed autocorrelation function. Then this field is transformed using u i = ¯ u i + a ij U j , (21) where a ij is a function of Reynolds tensor. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Large Eddy Simulation Numerical methods Subgrid models Parameters for computation of the short-scale scenario Other numerical methods Results Turbulent inflow At this time our implementation requires uniform grid in y and z . We use the Reynolds tensor components reported in literature for neutral boundary layer scaled by the friction velocity (i.e. σ u / u ∗ , σ v / u ∗ , σ w / u ∗ ). The mean wind and stress profiles can be prescribed by the user. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Numerical methods Parameters for computation of the short-scale scenario Results Parameters I. The grid is always uniform, although generally the code is capable of nonuniform grid, some of the methods can be used only for uniform grid (for now at least). Two domains for computation: up to 50 m and up to 2000 m (both with some additional buffer zone) These were 2 separate computations and therefore there is a discontinuity in the results. None of the two tests had favourable wind speed and directions. Values used for computations were avereaged from the measurement: 2.3 m/s with 30 degrees to the right and 0.4 m/s in the direction of the main axis respectively. The reference height for the wind was 2 m. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Numerical methods Parameters for computation of the short-scale scenario Results Parameters II. We used 4 separate classes of particle sizes that were treated independently each with fixed diameter and initial percentage. 0 . 2 µ m 10% 1 µ m 13% 6 µ m 64% 20 µ m 13% Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Numerical methods Parameters for computation of the short-scale scenario Results Parameters III. The incoming flow had logaritmic profile. The turbulent stresses were computed using roughness length 10 cm. The roughness length for the ground surface was 1 cm. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Numerical methods Parameters for computation of the short-scale scenario Results Results The main feature that is common to all our computations is much weaker decay of deposite activity with distance than measured (in other scenarios with released data). This could be changed by allowing more larger particles in the initial cloud but that changes the maximum values. We would like to know opinions and results of other modellers. Also the 2 km run was at rather coarse resolution of 12.5 m. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Numerical methods Parameters for computation of the short-scale scenario Results Example result Deposition of activity T3 (grid lines every 100 m, aligned to the wind direction) Vladim´ ır Fuka, Josef Brechler Short scale scenario
Introduction Numerical methods Parameters for computation of the short-scale scenario Results Percentiles of contaminated area The percentiles of the contaminated area (only up to 2 km): test3 test4 50 % 3.95 ha 3.25 ha 75 % 12.2 ha 11.1 ha 95 % 25.9 ha 25.8 ha Our model does not compute γ dose rates and we would appreciate any help how to do that. Vladim´ ır Fuka, Josef Brechler Short scale scenario
Recommend
More recommend
