Comparison between test field data and Gaussian plume model
|
|
- Colleen Davis
- 6 years ago
- Views:
Transcription
1 Environmental Modelling for RAdiation Safety II Working group 9 Comparison between test field data and Gaussian plume model Laura Urso Helmholtz Zentrum München Institut für Strahlenschutz AG Radioecological Modelling and Retrospective Dosimetry(REM) Wien, January
2 Simulation program for determination of population exposure to high doses after the explosion of an RDD device ) Gaussian model with known metereological parameters ) With at least 3 TLD measurements the free parameters can be inversely determined 3) Mathematical approach for inverse modelling: Levenberg-Marquardt Algorithm Two calculated examples: A) Synthetic data produced with HOTSPOT.7 National project: Retrospective dosimetry for the population in emergency situations Contract No 367S456 Bundesamt für Strahlenschutz (BfS) Federal Ministry for the environment, Nature Conservation and Nuclear Safety (BMBF)
3 Equations from SSK report No. 37 pp Gaussian dispersion model χ(x, y, z; H) = Dispersion coefficients a,b,c depend on stability class (from HOTSPOT guide.7) Depletion factor (from HOTSPOT guide.7)) Wet deposition σ y,z (x) = DF(x) = W (x) = πσ y (x)σ z (x)u(x) }{{} ax ( + bx) c [ exp( x exp( Λ πσy (x)u(x) e exp y σy (x) } {{ } ( e (z H) σz (x) H σ σz(x ) z(x )) dx y σy (x) ) + e (z+h) σz (x) } {{ } 3 ] v du π Ground deposition B r (x, y) =Q r (v d DF(x) χ(x, y, ) + W (x))e λ r t Dose conversion factors: submersion gw,r, inhalation gh,r, deposition gb,r submersion dose inhalation dose deposition dose External dose H wr (x, y, z) =Q r χ(x, y, z)g wr H hr (x, y, ) = Q r χ(x, y, z)g hr H br (x, y, ) = Q r (χ(x, y, )v d DF(x)+W (x))b g br K br H tot (x, y, z) =H wr (x, y, z)+h br (x, y, )
4 Dose conversion factors (Zähringer-Sempau BfS-IAR-/97) submersion gw,r (Gy s / Bq m 3 ) TABLE A.3, deposition gb,r (Gy s / Bq m ) TABLE A. Tc-99m Xe-33 Te-3 I-3 Cs-37 SZ Cloud Emitted photon energy (kev) DCF deposition log((gy m )/(Bq s)) Te-3 Tc-99m I-3 Cs-37 6 ZS DRY ZS WET Emitted photon energy (kev) Source homogeneously distributed in the air Source exponentially distributed in the soil with relaxation mass per unit area β DRY =. g/cm WET = g/cm
5 CODE: OPTLMDOSE.f9 MAIN PROGRAM SUBROUTINE FCN.f9 calculates objective function as log(dosedata) - log(dose) SUBROUTINE LMDIF from MINPACK runs optimisation SUBROUTINE COVAR calculates covariance matrix for error estimation cartesian axis: wind direction is x-axis INPUT data namelist: &global_para rnuclide='tc-99m' wind_ref=3.3d theta=.d stability_class = 'A' H=.5d vd=.d h_ref=.d dep_model= DRY I_rain =.d eq_model= EXPONENTIALX xdata=.d Dt_plot=6.d Qr = 5.8D8/ filename_read: x (km), y(km), Dose(Sv), Surface activity (kbq/m ), Dt(s) Radionuclide implemented are: Cs-37, I-3, Xe-33, Te-3, Tc-99m OUTPUT data filename_save: info M 3 N opt_value NORM.6464 unbiased sigmax other output files to produce plots
6 CODE: LEVENBERG-MARQUARDT ALGORITHM in MINPACK F = log(dosedata)-log(dose) If xsol is a solution of a non-linear least square problem then x solves: and orthogonality condition is valid F (x sol ) T F (x sol ) = The algorithm looks for a correction p such that F(x+p) F(x) To find appropriate p, the algorithm solves the problem: min{ f=j p : D p Δ} where D is diagonal scaling matrix and Δ is a step bound m i= f i (x) f i (x) = LMDIF runs various convergency tests between approximation x and the solution xsol INFO : if the final norm of the residual has K significant decimal digits compared to initial one (the assumed tolerance -K is set to square root of machine precision) INFO : the larger components of (D x) have K significant digits compared to initial ones INFO 3: if both and are fulfilled INFO 4: if the norm of the residuals is orthogonal to the Jacobian matrix.this should be examined further: could be F(x)=, some local minimum and accuracy is not implicit
7 Test data from HOTSPOT Surface Activity 5 8 y (m) Total Dose 5 e y (m) e 9 5 e
8 Test data from HOTSPOT - In principle with identical input values and initial guess, initial value for surface activity and dose should be the same as test data. BUT there is a difference of about 4-5% between the two DEPLETION FACTOR x [m] - HOTSPOT CODE and OP_LM_BfS almost identical: the only difference is integration for Depletion factor! - In OPT_LM_BfS GAUSS integration is used to increase the number of steps during integration. HOSPOT uses trapezoidal rule but no possibility to check it
9 Test data from HOTSPOT: result info M 3 N opt_value NORM.6464 unbiased sigmax (sigmax divided by M-N) convergence achieved after 9 iterations 7.5 Exp Dose INI Dose OPT Dose 8 Dose log(sv) x log(m)
10 Test data from HOTSPOT: cloudshine and groundshine Dose cloudshine Dose groundshine 5 e 5 y (m) e 9 y (m) e e Groundshine/Cloudshine
11 Test data from HOTSPOT: results.5. RESIDUAL PLOT OPT RESIDUAL INI RESIDUAL During optimisation, the residuals decrease and mean value goes to zero (.7 - ) The norm of the residual decreases from.6 to.6 (y exp y model ) There is a clear trend in the residual plot - residual is not random! Uncertainty on source term decreases with increasing the number of points Small uncertainty in the result of the fit has to be expected as by fixing the meteorological data the shape of the curve is fixed Optimised source term x points 5 points points 7 points 3 points points point Number of experimental points
12 Experimental data from Prouza et al. TEST : measurements of surface activity &global_para rnuclide='tc-99m' wind_ref=.d theta=.d stability_class ='B' H= 5.d vd=. h_ref=.d I_rain =.d Dt=45.d Qr = 9.D4/ - large uncertainty on deposition velocity vd - Initial source term (measured) is 9 MBq - Along x-direction, experimental profile of the plume is NOT an exponential - --> slow wind and clear Gaussian profile suggest diffusive process also in X direction: possibility for users to choose! - No source partitioning is included - objective function which is minimised is log(bdata+) - log(br+) 8 6 Experiment Gauss fit Surface activity (Bq/m ) diffusivex ( ) ( ) ( ) x x / ó x x = e Fitted Gaussian profile x=5 m, σx=
13 Experimental data from Prouza et al. TEST : results info M N opt_value unbiased sigmax Result strongly depends on deposition velocity for v_d =.8 m/s result is not physical anymore? CONTOUR PLOT log (B r +) INI CONTOUR PLOT log (B +) INI r CONTOUR PLOT log (B +) OPT r Optimised source term log(bq) 9 v d =.5 m/s.5 m/s. m/s.7 m/s 8.8 m/s. m/s.5 m/s. m/s.8 m/s Deposition velocity log(m/s)
14 Experimental data from Prouza et al. TEST EXPERIMENTAL INITIAL OPTIMAL INITIAL FINAL Surface activity log (B r + ) Surface activity log (B r + ) RESIDUAL y (m) The norm of the residual decreases from 4 to The residual clearly follows a trend and is not random but around zero The standard deviation is very small...again by fixing meteorological data the form of the curve underlying the fit is fixed!
Comparison between test field data and Gaussian plume model
Environmental Modelling for RAdiation Safety II Working group 9 Comparison between test field data and Gaussian plume model Dr. Laura Urso Helmholtz Zentrum München Institut für Strahlenschutz AG Risikoanalyse
More informationSimulation of Kamenna RDD field tests with the RDD_MCC computer code
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
More informationSIMULATION OF KAMENNA-EXPERIMENTS TEST 1 4 WITH THE DECISION SUPPORT MODEL LASAIR
SIMULATION OF KAMENNA-EXPERIMENTS TEST 1 4 WITH THE DECISION SUPPORT MODEL LASAIR IAEA, EMRAS II Working Group 9 Urban Areas Sevilla, Spain, 08. - 10.06. 2010 Federal Office for Radiation Protection Ingolstaedter
More informationCTBT: Science and Technology 2017 Conference June 2017, Vienna, Austria
ROSATOM STATE CORPORATION FOR ATOMIC ENERGY Express-Method for Identification of Underground Nuclear Explosion-Generated Radionuclide Seepage Area during On-Site Inspection Igor V. Markov Gennadi V. Orlov
More informationINTERNAL RADIATION DOSIMETRY
INTERNAL RADIATION DOSIMETRY Introduction to Internal Dosimetry Importance of External and Internal Exposures by Radiation Type Charged particle radiation (α, p, β) Generally important only for internal
More informationTHE DECISION SUPPORT SYSTEM LASAIR: NEW FEATURES FOR EVALUATING DIRTY BOMB SCENARIOS
THE DECISION SUPPORT SYSTEM LASAIR: NEW FEATURES FOR EVALUATING DIRTY BOMB SCENARIOS 16th HARMO 2014 08-11 September 2014 Varna, Bulgaria Federal Office for Radiation Protection Ingolstaedter Landstrasse
More informationAP1000 European 15. Accident Analyses Design Control Document EVALUATION MODELS AND PARAMETERS FOR ANALYSIS OF RADIOLOGICAL CONSEQUENCES OF ACCIDENTS
APPENDIX 15A EVALUATION MODELS AND PARAMETERS FOR ANALYSIS OF RADIOLOGICAL CONSEQUENCES OF ACCIDENTS This appendix contains the parameters and models that form the basis of the radiological consequences
More informationCFD calculations of the test 2-4 experiments. Author: G. de With
CFD calculations of the test 2-4 experiments Author: G. de With 34. Model setup and boundary conditions Dimensions CFD model: x=1000m / y=100m / z=2000m. CFD Model: Transient simulation, with steady-state
More information13. Nonlinear least squares
L. Vandenberghe ECE133A (Fall 2018) 13. Nonlinear least squares definition and examples derivatives and optimality condition Gauss Newton method Levenberg Marquardt method 13.1 Nonlinear least squares
More informationDose Reconstruction Methods and Source Term Assessment using Data from Monitoring Networks and Mobile Teams A German Approach
Dose Reconstruction Methods and Source Term Assessment using Data from Monitoring Networks and Mobile Teams A German Approach M. Bleher, U. Stöhlker, F. Gering Federal Office for Radiation Protection (BfS),
More informationStudy of the Effect of Weather, Topography and Radionuclide on the TEDE in a Fire Scenario Involving a Dispersion of a Plume in the Atmosphere
Middle-East Journal of Scientific Research 18 (8): 1192-1198, 2013 ISSN 1990-9233 IDOSI Publications, 2013 DOI: 10.5829/idosi.mejsr.2013.18.8.12405 Study of the Effect of Weather, Topography and Radionuclide
More informationDose Calculations in ARGOS
Dose Calculations in ARGOS Table of Contents Introduction... 2 Dose Pathways... 2 Dose Types... 2 Sheltering... 3 Age Groups... 3 Integration Time... 4 Database Tables used by ARGOS Dose Calculations...
More informationMODELING FOR ENVIRONMENTAL RADIATION DOSE RECONSTRUCTION. Bruce Napier 23 May 2011
MODELING FOR ENVIRONMENTAL RADIATION DOSE RECONSTRUCTION Bruce Napier 23 May 2011 1 Topics NCRP Report No. 163 Atmospheric dispersion modeling Reconstruction of dose from releases of iodines and noble
More informationAssessment of atmospheric dispersion and radiological consequences for the Fukushima Dai-ichi Nuclear Power Plant accident
Assessment of atmospheric dispersion and radiological consequences for the Fukushima Dai-ichi Nuclear Power Plant accident IRPA 13 - Fukushima session A.Mathieu, I.Korsakissok, D.Quélo, J.Groëll, M.Tombette,
More informationNon-linear least squares
Non-linear least squares Concept of non-linear least squares We have extensively studied linear least squares or linear regression. We see that there is a unique regression line that can be determined
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run
More informationModeling the Physical Processes that Impact the Fate and Fallout of Radioactive Materials
Modeling the Physical Processes that Impact the Fate and Fallout of Radioactive Materials J.V. Ramsdell, Jr. Radiological Science and Engineering Group Pacific Northwest National Laboratory Richland, Washington
More informationSURVEY OF THE 137 Cs CONTAMINATION IN BELGIUM AFTER THE CHERNOBYL ACCIDENT BY IN-SITU GAMMA SPECTROMETRY
SURVEY OF THE 137 Cs CONTAMINATION IN BELGIUM AFTER THE CHERNOBYL ACCIDENT BY IN-SITU GAMMA SPECTROMETRY J. Uyttenhove 1, S. Pommé 2, B. Van Waeyenberge 3 1 University Gent, Physics Lab, Krijgslaan 281
More informationADAPTATION OF A LAGRANGIAN PARTICLE-DISPERSION MODEL FOR THE USE IN RADIOPROTECTION AND RADIOECOLOGY
ADAPTATION OF A LAGRANGIAN PARTICLE-DISPERSION MODEL FOR THE USE IN RADIOPROTECTION AND RADIOECOLOGY Christoph Haustein TÜV Industrie Service GmbH TÜV Süd Gruppe, Munich, Germany INTRODUCTION The assessment
More informationKarlsruhe Chart of the Nuclides Nuclear Data. Zsolt Sóti Institute for Transuranium Elements - Karlsruhe Joint Research Centre
Karlsruhe Chart of the Nuclides Nuclear Data Zsolt Sóti Institute for Transuranium Elements - Karlsruhe Joint Research Centre Overview Nuclide Charts Karlsruhe Chart of Nuclides Electronic Nuclide Charts
More informationNumerical methods, midterm test I (2018/19 autumn, group A) Solutions
Numerical methods, midterm test I (2018/19 autumn, group A Solutions x Problem 1 (6p We are going to approximate the limit 3/2 x lim x 1 x 1 by substituting x = 099 into the fraction in the present form
More informationDr. Ryohji Ohba (Nuclear Safety Research Association )
Survey of estimation method for amount of radioactive materials emitted from nuclear power station during severe accident Dr. Ryohji Ohba (Nuclear Safety Research Association ) Funded by Ministry of Education
More informationICONE DRAFT - EVALUATION OF THE CONSERVATISM OF NPP SAFETY ANALYSIS DOSE CALCULATIONS AS TYPICAL FOR LICENSING PURPOSES
Proceedings of the 18th International Conference on Nuclear Engineering ICONE18 May 17-21, 21, Xi'an, China ICONE18-29465 DRAFT - EVALUATION OF THE CONSERVATISM OF NPP SAFETY ANALYSIS DOSE CALCULATIONS
More informationRoot Finding (and Optimisation)
Root Finding (and Optimisation) M.Sc. in Mathematical Modelling & Scientific Computing, Practical Numerical Analysis Michaelmas Term 2018, Lecture 4 Root Finding The idea of root finding is simple we want
More informationAssessment of atmospheric dispersion and radiological consequences for the Fukushima Dai-ichi Nuclear Power Plant accident
Assessment of atmospheric dispersion and radiological consequences for the Fukushima Dai-ichi Nuclear Power Plant accident IRPA 13 - Fukushima session A.Mathieu, I.Korsakissok, D.Quélo, J.Groëll, M.Tombette,
More informationAnalysis on the formation process of high dose rate zone in the northwest direction of the Fukushima Daiichi nuclear power plant
Analysis on the formation process of high dose rate zone in the northwest direction of the Fukushima Daiichi nuclear power plant Jun. 13, 2011 Research group of Japan Atomic Energy Agency (JAEA) has analyzed
More informationMeteorological Data Collection, X/Q and D/Q, Critical Receptors
Meteorological Data Collection, X/Q and D/Q, Critical Receptors Ken Sejkora Entergy Nuclear Northeast Pilgrim Station Presented at the 23rd Annual RETS-REMP Workshop Training Session Westminster, CO /
More informationComputational Methods. Least Squares Approximation/Optimization
Computational Methods Least Squares Approximation/Optimization Manfred Huber 2011 1 Least Squares Least squares methods are aimed at finding approximate solutions when no precise solution exists Find the
More informationFixed Points and Contractive Transformations. Ron Goldman Department of Computer Science Rice University
Fixed Points and Contractive Transformations Ron Goldman Department of Computer Science Rice University Applications Computer Graphics Fractals Bezier and B-Spline Curves and Surfaces Root Finding Newton
More informationIntroduction to Probability and Stocastic Processes - Part I
Introduction to Probability and Stocastic Processes - Part I Lecture 2 Henrik Vie Christensen vie@control.auc.dk Department of Control Engineering Institute of Electronic Systems Aalborg University Denmark
More informationO R D E R OF THE HEAD OF THE STATE NUCLEAR POWER SAFETY INSPECTORATE
O R D E R OF THE HEAD OF THE STATE NUCLEAR POWER SAFETY INSPECTORATE ON THE APPROVAL OF NUCLEAR SAFETY REQUIREMENTS BSR-1.9.1-2011 STANDARDS OF RELEASE OF RADIONUCLIDES FROM NUCLEAR INSTALLATIONS AND REQUIREMENTS
More informationRapid Extraction of Plutonium from Urine by Pyrosulfate Fusion and
Rapid Extraction of Plutonium from Urine by Pyrosulfate Fusion and PERALS Spectroscopy R.L. Metzger, P.H. Pouquette, and G.W. Klingler Radiation Safety Engineering, Inc. Chandler, AZ Abstract To effectively
More informationCALIBRATION OF INSTRUMENTS MEASURING RADON OVER A LARGE ACTIVITY RANGE. N. Michielsen, V. Voisin
Radon in the Living Environment, 011 CALIBRATION OF INSTRUMENTS MEASURING RADON OVER A LARGE ACTIVITY RANGE N. Michielsen, V. Voisin Institut de Protection et de Sûreté Nucléaire, Département de Prévention
More informationMath and Numerical Methods Review
Math and Numerical Methods Review Michael Caracotsios, Ph.D. Clinical Associate Professor Chemical Engineering Department University of Illinois at Chicago Introduction In the study of chemical engineering
More informationNumerical solution of Least Squares Problems 1/32
Numerical solution of Least Squares Problems 1/32 Linear Least Squares Problems Suppose that we have a matrix A of the size m n and the vector b of the size m 1. The linear least square problem is to find
More informationNuclear Spectroscopy: Radioactivity and Half Life
Particle and Spectroscopy: and Half Life 02/08/2018 My Office Hours: Thursday 1:00-3:00 PM 212 Keen Building Outline 1 2 3 4 5 Some nuclei are unstable and decay spontaneously into two or more particles.
More informationJoint Distributions. (a) Scalar multiplication: k = c d. (b) Product of two matrices: c d. (c) The transpose of a matrix:
Joint Distributions Joint Distributions A bivariate normal distribution generalizes the concept of normal distribution to bivariate random variables It requires a matrix formulation of quadratic forms,
More informationData Fitting and Uncertainty
TiloStrutz Data Fitting and Uncertainty A practical introduction to weighted least squares and beyond With 124 figures, 23 tables and 71 test questions and examples VIEWEG+ TEUBNER IX Contents I Framework
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationRADIATION MONITORING NETWORKS and some results from the reactor accident in Fukushima
RADIATION MONITORING NETWORKS and some results from the reactor accident in Fukushima Sixth Annual Workshop on Radiation Measurement Cross Calibration Project June 13-15, 2011 Middle East Scientific Institute
More informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
More informationIn-situ measurement of Cs distribution in the soil. Sabina Markelj
In-situ measurement of Cs distribution in the soil Sabina Markelj Supervisor: Matej Lipoglavšek April 12, 2004 Abstract In-situ gamma spectrometry is a new method for measuring radionuclides in soil. It
More informationLaboratoire National Henri Becquerel (CEA/LIST/LNHB), France (2) ENEA-Radiation Protection Institute, Bologna, Italy (3)
PROPOSAL FOR EYE-LENS DOSEMETER CALIBRATION AND TYPE TESTING ORAMED WP2 J.-M. Bordy (1), J. Daures (1), M. Denozière (1), G. Gualdrini (2), M. Guijaume (3), E. Carinou (4),F. Vanhavere (5) (1) Laboratoire
More informationRSMC MONTRÉAL USERS' INTERPRETATION GUIDELINES ATMOSPHERIC TRANSPORT MODEL OUTPUTS. Version 9
RSMC MONTRÉAL USERS' INTERPRETATION GUIDELINES ATMOSPHERIC TRANSPORT MODEL OUTPUTS Version 9 Environmental Emergency Response Section RSMC Montréal Canadian Meteorological Centre Meteorological Service
More informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
More informationminimize x x2 2 x 1x 2 x 1 subject to x 1 +2x 2 u 1 x 1 4x 2 u 2, 5x 1 +76x 2 1,
4 Duality 4.1 Numerical perturbation analysis example. Consider the quadratic program with variables x 1, x 2, and parameters u 1, u 2. minimize x 2 1 +2x2 2 x 1x 2 x 1 subject to x 1 +2x 2 u 1 x 1 4x
More informationTritium Environmental Assessment for Fusion Reactor Using TAS3.0
Tritium Environmental Assessment for Fusion Reactor Using TAS3.0 Baojie Nie, Muyi Ni, Minghuang Wang, Jing Song Contributed by FDS Team Key laboratory of Neutronics and Radiation Safety Institute of Nuclear
More informationMeteorology and Dispersion Forecast in Nuclear Emergency in Argentina
Meteorology and Dispersion Forecast in Nuclear Emergency in Argentina Kunst, J.J.; Boutet, L.I.; Jordán, O.D.; Hernández, D.G.; Guichandut, M.E. and Chiappesoni, H. Presentado en: 12 th International Congress
More informationRadiological impact of rutile covered welding electrodes
LABORATORIO DE MEDIDAS DE BAJA ACTIVIDAD Universidad del País Vasco / Euskal Herriko Unibertsitatea Radiological impact of rutile covered welding electrodes Margarita Herranz 1, Saroa Rozas 1, Carmen Pérez
More informationExpert judgement and uncertainty quantification for climate change
Expert judgement and uncertainty quantification for climate change 1. Probabilistic inversion as applied to atmospheric dispersion Atmospheric dispersion is mathematically simple, yet involves many features
More informationName: Date: Block: Algebra 1 Function/Solving for a Variable STUDY GUIDE
Algebra Functions Test STUDY GUIDE Name: Date: Block: SOLs: A.7, A. Algebra Function/Solving for a Variable STUDY GUIDE Know how to Plot points on a Cartesian coordinate plane. Find the (x, y) coordinates
More informationRadioactivity in the Risø District January-June Risø-R-Report. DTU Nutech Report
Radioactivity in the Risø District January-June 2012 Risø-R-Report DTU Nutech Report Sven P. Nielsen, Kasper G. Andersson and Arne Miller DTU-Nutech-R-0003(EN) December 2012 Author: Sven P. Nielsen, Kasper
More informationHence a root lies between 1 and 2. Since f a is negative and f(x 0 ) is positive The root lies between a and x 0 i.e. 1 and 1.
The Bisection method or BOLZANO s method or Interval halving method: Find the positive root of x 3 x = 1 correct to four decimal places by bisection method Let f x = x 3 x 1 Here f 0 = 1 = ve, f 1 = ve,
More informationEvaluation of detection abilities of monitoring networks using multiple assessment criteria. Radek Hofman, Petr Pecha* and Václav Šmídl
Int. J. Environment and Pollution, Vol. 55, Nos. 1/2/3/4, 2014 129 Evaluation of detection abilities of monitoring networks using multiple assessment criteria Radek Hofman, Petr Pecha* and Václav Šmídl
More informationUser Guide for the Terrestrial Food Chain and Dose Module FDMT of
Version 1.1 User Guide for the Terrestrial Food Chain and Dose Module FDMT of RODOS PV6.0 User Guide for the Terrestrial Food Chain and Dose Module FDMT of RODOS PV6.0 Florian Gering, Heinz Müller USER
More informationLATE PHASE OF NUCLEAR ACCIDENT IN AGRICULTURE WEAKNESSES IN DECISION MAKING DUE TO LACK OF DATA, KNOWLEDGE AND SKILL AND POSSIBLE INNOVATION
IAEA expert meeting on Radiation Protection after Fukushima Accident, Vienna 17-21 Feb 2014 LATE PHASE OF NUCLEAR ACCIDENT IN AGRICULTURE WEAKNESSES IN DECISION MAKING DUE TO LACK OF DATA, KNOWLEDGE AND
More informationNuclear risks from atmospheric dispersion in Northern Europe
NKS NordRisk: Nuclear risks from atmospheric dispersion in Northern Europe Bent Lauritzen Risø DTU Presented at the NKS meeting 26-27 March, Stockholm NKS NordRisk NordRisk 2005-2006 NordRisk II 2008-2009
More informationCOURSE Iterative methods for solving linear systems
COURSE 0 4.3. Iterative methods for solving linear systems Because of round-off errors, direct methods become less efficient than iterative methods for large systems (>00 000 variables). An iterative scheme
More information8.) Our radiological environment
8.) Our radiological environment Te average effective annual dose Indoor radon Rn-222 is te biggest problem, because of a relatively large abundance of its element of origin, U- 238, and a relatively long
More information5 Operations on Multiple Random Variables
EE360 Random Signal analysis Chapter 5: Operations on Multiple Random Variables 5 Operations on Multiple Random Variables Expected value of a function of r.v. s Two r.v. s: ḡ = E[g(X, Y )] = g(x, y)f X,Y
More informationMAPPING OF DOSE-RATE OF GAMMA-IRRADIATION BY 137 Cs OVER THE EARTH SURFACE FOR THE TERRITORY OF RUSSIA
MAPPING OF DOSE-RATE OF GAMMA-IRRADIATION BY 137 Cs OVER THE EARTH SURFACE FOR THE TERRITORY OF RUSSIA Yu.A.Izrael, E.V.Kvasnikova, I.M.Nazarov, A.N.Pegoev, E.D.Stukin Institute of Global Climate and Ecology,
More informationJUST THE MATHS UNIT NUMBER INTEGRATION APPLICATIONS 10 (Second moments of an arc) A.J.Hobson
JUST THE MATHS UNIT NUMBER 13.1 INTEGRATION APPLICATIONS 1 (Second moments of an arc) by A.J.Hobson 13.1.1 Introduction 13.1. The second moment of an arc about the y-axis 13.1.3 The second moment of an
More informationFeb 21 and 25: Local weighted least squares: Quadratic loess smoother
Feb 1 and 5: Local weighted least squares: Quadratic loess smoother An example of weighted least squares fitting of data to a simple model for the purposes of simultaneous smoothing and interpolation is
More informationRegression. Oscar García
Regression Oscar García Regression methods are fundamental in Forest Mensuration For a more concise and general presentation, we shall first review some matrix concepts 1 Matrices An order n m matrix is
More informationAlgebra. Pang-Cheng, Wu. January 22, 2016
Algebra Pang-Cheng, Wu January 22, 2016 Abstract For preparing competitions, one should focus on some techniques and important theorems. This time, I want to talk about a method for solving inequality
More informationM. Rogozina, M. Zhukovsky, A. Ekidin, M. Vasyanovich. Institute of Industrial Ecology, Ural Branch Russian Academy of Sciences
THORON PROGENY SIZE DISTRIBUTION IN MONAZITE STORAGE FACILITY M. Rogozina, M. Zhukovsky, A. Ekidin, M. Vasyanovich Institute of Industrial Ecology, Ural Branch Russian Academy of Sciences Thoron ( 220
More informationAssessment of Radiological Dose around a 3-MW TRIGA Mark-II Research Reactor
International Letters of Chemistry, Physics and Astronomy Online: 2013-06-29 ISSN: 2299-3843, Vol. 15, pp 183-200 doi:10.18052/www.scipress.com/ilcpa.15.183 2013 SciPress Ltd., Switzerland Assessment of
More informationEnvironmental radiation measurements immediately after the accident and dose evaluations based on soil deposition
Environmental radiation measurements immediately after the accident and dose evaluations based on soil deposition Masahiro YOSHIDA* and Tominori SUZUKI Nuclear Safety Technology Center, 5-1-3-101 Hakusan
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationBACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES
No. of Printed Pages : 5 BCS-054 BACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 058b9 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES Time : 3 hours Maximum Marks
More informationMinimum Detectable and Maximum Missable Activities
and Maximum Missable The HPA Approach David Williams 18 May 2005 Centre for Radiation, Chemical and Environmental Hazards Radiation Protection Division formerly the National Radiological Protection Board
More informationCOMPARISON OF PERSONNEL RADIATION MONITORING DOSIMETERS DESIGNED FOR MEDICAL FIELD
COMPARISON OF PERSONNEL RADIATION MONITORING DOSIMETERS DESIGNED FOR MEDICAL FIELD Kirill SKOVORODKO, Birute GRICIENE, Milda PETKELYTE Radiation Protection Division, Vilnius University Hospital Santaros
More informationINTRODUCTION TO COMPUTATIONAL MATHEMATICS
INTRODUCTION TO COMPUTATIONAL MATHEMATICS Course Notes for CM 271 / AMATH 341 / CS 371 Fall 2007 Instructor: Prof. Justin Wan School of Computer Science University of Waterloo Course notes by Prof. Hans
More informationOptimization: Nonlinear Optimization without Constraints. Nonlinear Optimization without Constraints 1 / 23
Optimization: Nonlinear Optimization without Constraints Nonlinear Optimization without Constraints 1 / 23 Nonlinear optimization without constraints Unconstrained minimization min x f(x) where f(x) is
More informationNORM Monitor-IS. Performance data
NORM Monitor-IS Performance data Contents 1. Introduction... 3 Scintillator probe... 3 GM probe... 3 2. Energy response... 4 Scintillator probe gamma energy response... 4 Scintillator probe alpha/beta
More informationIterative Methods for Ax=b
1 FUNDAMENTALS 1 Iterative Methods for Ax=b 1 Fundamentals consider the solution of the set of simultaneous equations Ax = b where A is a square matrix, n n and b is a right hand vector. We write the iterative
More informationCounting Statistics and Error Prediction
Counting Statistics and Error Prediction Outline Characterization of data Statistical models Error propagation Optimization of counting experiments Energy resolution Detection efficiency Dead time Characterization
More informationSTUDIES ON THE AEROSOL RADIOACTIVITY LEVEL AND AIR QUALITY AROUND NUCLEAR RESEARCH INSTITUTE AREA
STUDIES ON THE AEROSOL RADIOACTIVITY LEVEL AND AIR QUALITY AROUND NUCLEAR RESEARCH INSTITUTE AREA A. STOCHIOIU 1, F. MIHAI 1, C. STOCHIOIU 2 1 Horia Hulubei National Institute for Physics and Nuclear Engineering,
More informationResults of the (i) Fifth Airborne Monitoring Survey and (ii) Airborne Monitoring Survey Outside 80km from the Fukushima Dai-ichi NPP
September 28, 2012 Results of the (i) Fifth Airborne Monitoring Survey and (ii) Airborne Monitoring Survey Outside 80km from the The results of the airborne monitoring survey within a 80km radius of the
More informationNuclear Medicine Intro & Physics from Medical Imaging Signals and Systems, Chapter 7, by Prince and Links
Nuclear Medicine Intro & Physics from Medical Imaging Signals and Systems, Chapter 7, by Prince and Links NM - introduction Relies on EMISSION of photons from body (versus transmission of photons through
More informationQuantifying Electrostatic Resuspension of Radionuclides from Surface Contamination
Quantifying Electrostatic Resuspension of Radionuclides from Surface Contamination Shaun Marshall 1, Charles Potter 2, David Medich 1 1 Worcester Polytechnic Institute, Worcester, MA 01609 2 Sandia National
More informationSome nuclei are unstable Become stable by ejecting excess energy and often a particle in the process Types of radiation particle - particle
Radioactivity George Starkschall, Ph.D. Lecture Objectives Identify methods for making radioactive isotopes Recognize the various types of radioactive decay Interpret an energy level diagram for radioactive
More informationImportance of uncertainties in dose assessment to prove compliance with radiation protection standards
Importance of uncertainties in dose assessment to prove compliance with radiation protection standards Manfred TSCHURLOVITS, Atominstitute of Austrian Universities, Vienna University of Technology, Stadionallee
More information3. Gas Detectors General introduction
3. Gas Detectors 3.1. General introduction principle ionizing particle creates primary and secondary charges via energy loss by ionization (Bethe Bloch, chapter 2) N0 electrons and ions charges drift in
More informationContinuous Optimization
Continuous Optimization Sanzheng Qiao Department of Computing and Software McMaster University March, 2009 Outline 1 Introduction 2 Golden Section Search 3 Multivariate Functions Steepest Descent Method
More informationIterative Methods. Splitting Methods
Iterative Methods Splitting Methods 1 Direct Methods Solving Ax = b using direct methods. Gaussian elimination (using LU decomposition) Variants of LU, including Crout and Doolittle Other decomposition
More informationUncertainty in Measurements
Uncertainty in Measurements Joshua Russell January 4, 010 1 Introduction Error analysis is an important part of laboratory work and research in general. We will be using probability density functions PDF)
More informationPresenting Uncertain Information in Radiological Emergencies. University of Warwick, Public Health England, Met Office
Presenting Uncertain Information in Radiological Emergencies University of Warwick, Public Health England, Met Office 20 th September 2016, 9.30am 12.30pm Accident scenario example Presentation of theoretical
More informationPolytechnic Institute of NYU MA 2132 Final Practice Answers Fall 2012
Polytechnic Institute of NYU MA Final Practice Answers Fall Studying from past or sample exams is NOT recommended. If you do, it should be only AFTER you know how to do all of the homework and worksheet
More informationAssignment 1, SOLUTIONS
MAT 10: Calculus for the Life Sciences I 01.10.2008 Pawel Lorek Assignment 1, SOLUTIONS University of Ottawa Problem 1: [4 points] Suppose that every morning a patient receives the same dose of drug. From
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationA MATLAB Implementation of the. Minimum Relative Entropy Method. for Linear Inverse Problems
A MATLAB Implementation of the Minimum Relative Entropy Method for Linear Inverse Problems Roseanna M. Neupauer 1 and Brian Borchers 2 1 Department of Earth and Environmental Science 2 Department of Mathematics
More informationApplied Mathematics 205. Unit I: Data Fitting. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit I: Data Fitting Lecturer: Dr. David Knezevic Unit I: Data Fitting Chapter I.4: Nonlinear Least Squares 2 / 25 Nonlinear Least Squares So far we have looked at finding a best
More informationINTRODUCTION, FOUNDATIONS
1 INTRODUCTION, FOUNDATIONS ELM1222 Numerical Analysis Some of the contents are adopted from Laurene V. Fausett, Applied Numerical Analysis using MATLAB. Prentice Hall Inc., 1999 2 Today s lecture Information
More informationMeteorological monitoring system for NPP Kozloduy. Hristomir Branzov
Meteorological monitoring system for NPP Kozloduy Hristomir Branzov National Institute of Meteorology and Hydrology 66 Tzarigradsko schaussee, Sofia 784, Bulgaria Tel. (+359 2) 975 35 9, E-mail: Hristomir.Branzov@meteo.bg
More informationQuestion. 1. Which natural source of background radiation do you consider as dominant?
Question 1. Which natural source of background radiation do you consider as dominant? 2. Is the radiation background constant or does it change with time and location? 3. What is the level of anthropogenic
More informationGate simulation of Compton Ar-Xe gamma-camera for radionuclide imaging in nuclear medicine
Journal of Physics: Conference Series PAPER OPEN ACCESS Gate simulation of Compton Ar-Xe gamma-camera for radionuclide imaging in nuclear medicine To cite this article: L Yu Dubov et al 2017 J. Phys.:
More informationSLOVENIAN ALARM AND MONITORING NETWORKS
SLOVENIAN ALARM AND MONITORING NETWORKS Jasmina Kožar Logar, Denis Glavič-Cindro "Jožef Stefan" Institute,, Slovenia CHERNOBYL 25 YEARS LATER: STUDIES, REMARKS AND RECENT FINDINGS June 21-22-23, 2011 UDINE
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationScientific Highlight February 2011
Scientific Highlight February 2011 co-ordinated with the Director of the Institute / Research Unit Institute/ Research Unit / Clinical Co-operation Group / Junior Research Group: Institute of Radiation
More information