Simulation of Kamenna RDD field tests with the RDD_MCC computer code - - PowerPoint PPT Presentation

simulation of kamenna rdd field tests with the rdd mcc
SMART_READER_LITE
LIVE PREVIEW

Simulation of Kamenna RDD field tests with the RDD_MCC computer code - - PowerPoint PPT Presentation

Presentation to the 5th Meeting of the EMRAS II Working Group Urban Areas Simulation of Kamenna RDD field tests with the RDD_MCC computer code Juraj ran, Department of Accident Management and Risk Assessment VUJE Inc., Trnava,


slide-1
SLIDE 1

1/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

Presentation to the 5th Meeting of the EMRAS II Working Group “Urban Areas”

Simulation of Kamenna RDD field tests with the RDD_MCC computer code

Juraj Ďúran, Department of Accident Management and Risk Assessment VUJE Inc., Trnava, Slovak Republic 24-28 January 2011, Vienna, Austria

slide-2
SLIDE 2

2/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

CONTENTS

  • 1. Introduction
  • 2. Mathematical model
  • 3. Uncertainties
  • 4. Input data
  • 5. Preliminary results of comparison
  • 6. Conclusion
slide-3
SLIDE 3

3/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

  • 1. INTRODUCTION
  • Simulations for Test1 and Test2
  • Calculations were carry out with RDD_MMC

computer code (Method Monte Carlo)

  • Main task for present time – to find an appropriate

model for volume source term and distribution of particles

  • Assumptions – model and its validation will be

completed to end of June

slide-4
SLIDE 4

4/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

  • 2. MATHEMATICAL MODEL

Model approach

  • Lagrangian particle model
  • Simulation of fluctuation by the Monte Carlo

method

  • Fortran 77 code

Meteorological wind field (time independent)

  • homogenous wind field
  • computer code CALMET
  • computer code MATHEW
slide-5
SLIDE 5

5/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

  • 2. MATHEMATICAL MODEL

Source term

  • Initial volume of source term can be composed

from up to 8 independent volumes

  • Up to 8 diameter of particles (aerosols)
  • Each diameter of particle have own density
  • Initial distribution of particles in sources (random,

hat, gaussian) Removal processes

  • Gravitational settlement, Dry and wet deposition,

Radioactive decay

slide-6
SLIDE 6

6/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

  • 2. MATHEMATICAL MODEL

Dose rate

  • Calculation for 1 mono-energetic photon emission
  • Dose rate [Gy/h] in air from ground exposure
  • Numerical integration of contribution from all area

(step 1 x 1 m)

  • Berger form of a build-up factor
  • Energy of gamma rays, E = 0.140510 MeV
slide-7
SLIDE 7

7/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

  • 3. UNCERTAINTIES
  • Meteorological uncertainties – wind speed, wind

direction, category of stability

  • Geometric shape and dimension of initial clouds
  • Distribution of volume activity in initial clouds
  • Distribution of the activity as function of the

aerosol (particle) diameters

  • Distribution of particles with different diameter as

function of the height

  • Density for each diameter of aerosols (particles)
slide-8
SLIDE 8

8/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

  • 4. INPUT DATA

We suppose that geometrical shape of initial source term is composed from 2 independent volumes Table: Simulation of source term volumes for Test2

No. Center of volume X,Y,Z [m] Length, Wide, Thickness [m] Fraction of emission [-] 1 0.0, 0.0, 1.0 13.0, 6.0, 2.0 0.20 2 0.0, 0.0, 3.5 4.0, 2.0, 3.0 0.80

slide-9
SLIDE 9

9/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

slide-10
SLIDE 10

10/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

Measured and used meteorological input data

Date Explosion time Wind speed Wind direction * Stability category 06-12-2007 12:45 pm 4.00 m/s SW C (ETE) 15-05-2008 11:30 am 0.59 m/s S A (ETE) 05-05-2009 12:22 pm 1.30 m/s WNW B (ETE) 14-07-2010 12:42 pm 0.20 m/s S A (estimate)

05-05-2009 – data from 1st meteorological station, * measured value, for calculation was used the direction South + 10 degrees ETE – meteorological station near to NPP Temelin

slide-11
SLIDE 11

11/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

List of simulated particles parameters

No. Diameter [µm] Density [g/cm3] Dry deposition velocity [m/s] Fraction of emission [-] 1 0.2 1.0 5.0E-03 0.20 2 1.0 1.0 1.5E-04 0.15 3 8.0 1.5 1.0E-03 0.50 4 20.0 2.0 8.0E-03 0.15

slide-12
SLIDE 12

12/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

Calculated deposition for TEST 1

0.00E+00 5.00E+06 1.00E+07 1.50E+07 2.00E+07 2.50E+07 3.00E+07

  • 5

5 10 15 20 25 30 35 40 45

slide-13
SLIDE 13

13/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

Calculated deposition for TEST 2

0.00E+00 2.00E+07 4.00E+07 6.00E+07 8.00E+07 1.00E+08 1.20E+08

  • 5

5 10 15 20 25 30 35 40 45

slide-14
SLIDE 14

14/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

  • 5. PRELIMINARY RESULTS OF COMPARISON
  • Test 1 – result for deposition are in approximately

agreement with measurements (max. value ≈ 107 Bq/m2)

  • Test 2 – result for deposition are not in agreement

with measurements (max. value ≈ 105 Bq/m2)

  • Main task for present time – to find an appropriate

model for volume source term and distribution of particles

  • Assumptions – model and its validation will be

completed to end of June

slide-15
SLIDE 15

15/28

WG 9 “Urban Areas”

VUJE, Inc. Okružná 5, 918 64 Trnava, Slovakia

Thank you for your attention