HPC Particle Transport Methodologies for Simulation of Nuclear Systems
|
|
- Anthony O’Brien’
- 6 years ago
- Views:
Transcription
1 HPC Particle Transport Methodologies for Simulation of Nuclear Systems Prof. Alireza Haghighat Virginia Tech Virginia Tech Transport Theory Group (VT 3 G) Director of Nuclear Engineering and Science Lab (NSEL) at Arlington Nuclear Engineering Program, Mechanical Engineering Department VT High Performance Computing Day, Blacksburg, VA, April 11,
2 Objective Particle Transport Theory Determine the expected number of particles in a phase space (d 3 rded) at time t: n( r, 3 E, ˆ, t) d rded z d 3 r dω Ω de r y x 7 independent variables: space (x, y, z), energy (E), direction (,, time (t) 2
3 Simulation Approaches Deterministic Methods Solve the Linear Boltzmann Equation (LBE) to obtain the expected flux in a phase space Statistical Monte Carlo Methods Perform particle transport experiments using random numbers (RN s) on a computer to estimate average properties of a particle in phase space 3
4 Deterministic Linear Boltzmann Equation Integro-differential form ) ˆ,, ( ) ˆ ',, ( ), ( ' ' 4 ) ( ) ˆ ',, ( ) ˆ ' ˆ, ', ( ' ' ) ˆ,, ( ), ( ) ˆ,, (. ˆ E r S E r E r d de E E r E E r d de E r E r E r f s streaming collision scattering fission Independent source 4 Ψ,, Ω,, Ω) Where, angular flux is defined by Particle speed Particle density Source depends on the unknown; requires iteration
5 Integro-differential - Solution Method Spatial variable Integrated over fine meshes using FD or FE methods Angular variable: Discrete Ordinates (Sn) method: A discrete set of directions { ˆ } and associated weights m {w m } are selected ˆ ˆ. (,, ) (, ) (,, ˆ r E r E r E ) q( r, E, ˆ ) m m m, g, A V ijk d 3 r V m, g ijk ( r ) 0 0 m m X V V Energy variable Integrate over energy intervals to prepare multigroup cross sections, σ g 1 1 Example: a small size shield: 100x100x100 (space)x 80 (directions) x 50 (groups) x 8 byte/word = 32 GB 5
6 Monte Carlo (MC) Methods Perform an experiment on a computer; exact simulation of a physical process, e.g., a shield Path length r ln t Type of collision s Scattering angle 2 1 t (isotropic scattering) 0 Fixed source S (r, E, Ω) absorbed absorbed Tally (count) absorbed Issue Precise average values; i.e., small relative uncertainty, x R, require large x computation time x Therefore, Variance Reduction techniques are needed for real world problems! 6
7 Deterministic (Det.) vs. Monte Carlo (MC) Item Det. MC Geometry Discrete/ Exact Exact Energy treatment cross section Direction Discrete Discrete/ Truncated series Exact Exact Input preparation Difficult simple Computer memory Large Small Computer time Small Large Numerical issues Convergence Statistical uncertainty Amount of information Large Limited Variance Reduction ---- Unknown parameters Parallel computing Complex Trivial 7
8 Why not? Hybrid methods? MC Det. (for automated variance reduction) Det. Det. (use of region dependent technique) Parallel processing? Det. memory partitioning, and domain decomposition MC embarrassingly parallel 8
9 VT 3 G Vision Development of accurate particle transport methods for real time simulation of nuclear systems (power, security & medicine) Year Methodology Computer code system Wall clock time 2015 MRT TITAN IR 2014 MRT RAPID 2013 MRT AIMS MRT Hybrid MC det. (AVR) Hybrid det. det. INSPCT s ADIES ( ) TITAN (n, ϒ) 1997 Hybrid MC det. (automated VR AVR) A 3 MCNP (n, ϒ) 1996 Parallel (3 D) PENTRAN (n, ϒ) 1992 Vector & parallel (2 D) 1989 Parallel processing (1 D) 1986 Vector processing (1 D) 9
10 PENTRAN (Parallel Environment Neutral-particle TRANsport) Code System (G. Sjoden and A. Haghighat, 1996) Pre-processing PENMSH-XP (prepares mesh, source, and material distributions) Transport Calculation Domain Decomposition FORTRAN 90 with MPI library Hybrid spatial, angular, energy decomposition algorithms Partition memory; parallel I/O Post-processing PENPRL (3-D linear interpolation for comparison with experimental data) 10 10
11 PENTRAN performance for the Box in the Box test problem water Neutron source Run Scalability Results using PENTRAN on a Constant Size Problem (BP4) All results converged to a 1.0E-5 relative tolerance. PENTRAN Decomposition Strategy # Procs Angle # Procs Group # Procs Space #Total Procs P Sp= Speedup Factor Ep= Efficiency (Sp /P, %) Comm Overhead: % w allclock 1 Angular * Group Spatial Angular Angular-Group Angular Spatial Angular-Group Angular-Spatial Angular Angular-Spatial Angular-Group-Spatial Angular-Group Angular-Spatial Spatial * Average of 2-batch and interactive runs at CNSF 11
12 PENTRAN Performance Amdahl s law Pf =
13 Benchmarking Application of PENTRAN TM Kobayashi 3-D Benchmarks VENUS-3 Benchmark facility in SCK.CEN, Belgium C5G7 criticality benchmark Real-world problems BWR core-shroud and internals Pulsed Gamma Neutron Activation Analysis (PGNAA) device for assaying of waste X-Ray room CT Scan Time-of-Flight (TOF) for cross section measurement Simulation of spent fuel storage CASK 13
14 Venus-3 experimental benchmark (performed in ) ~85,000 cells 26 Groups P3-S8 Discrete Ordinates 32 IBM SP2 Procs, 84 min Parallel Performance on IBM SP2 95% C/E values +/-10%; 5% within +/-15% Case Number of processors Decomposition algorithms (A/G/S) 1 Wall-clock Speedup Efficiency time (min) 2 (%) 1 4 4/1/ /1/ /1/ /1/ (A/G/S) refers to the number of angular, group, and spatial sub-domains. 2 Time is obtained in a BATCH mode. 14
15 PENTRAN BWR Core shroud Modeling (performed in 1999) 67 G, S8/P3; 265,264 mesh 48 Procs (8 A/6 S) on IBM SP2: Within 5-15% A 3 MCNP (continuous energy), Wall-clock time = 12 hours Speedup = 10 (compared to a standard code) Scalability Case No. of Directions No. of Processors Decomposition (A/G/S) 1 Wallclock/iteration (s) /1/ /1/ /1/ /1/ (A/G/S) refers to the number of Angular, Group, and Spatial subdomains. 15
16 PENTRAN Simulation Storage Cask CASK library (22n, 18g) 17 Materials 318,426 fine meshes (1000 coarse meshes) (40 z-levels) P 3, S 12 (168 directions) 1.48 GB per processor 8 processors (~12 GB Total) Model # CPU Dose Ratio Run Time (hrs) PENTRAN Large Model
17 TITAN Hybrid Sn & CM Algorithm (C. Yi, A. Haghighat, 2004) Partition the problem into coarse meshes and allow for the use of different methods in different coarse meshes; Discrete ordinates (Sn) Solver Characteristic Method (CM) solver TITAN is written in Fortran 90 (with some features in Fortran 2003 standard, such as dynamic memory allocation and object oriented23) and MPI library 17
18 TITAN Sn CM Algorithm Sn out gink gink CM in gis Q ink gin gisink ginke (1 e ) Q gi in out gin gink gink s gi ink gi A B E.g., LD Scheme; 2 x o u t x i n 2 y o u t y i n 2 z o u t z i n ( n ) g ijk ( n ) g ijk ( n ) g ijk gin g ink ( Aink sink ) Aink gink Q k gin 1 k ( A s ) ( A s ) k ink ink gi gi ink ink k P B AB A 18
19 TITAN Benchmarking and Application OECD/NEA Benchmarks C5G7 MOX Kobayashi 3 D parameter space VENUS 2 Applications Adjoint calculation for the AIMS active interrogation simulation tool mpower reactor core and external modeling Modeling of a penetration duct in a nuclear reactor Benchmarking the multigroup SDM (subgroup decomposition method) algorithm (developed by Georgia Tech) Medical applications Modeling a CT machine Developed an image reconstruction algorithm, TITAN IR 19
20 TITAN simulation of CT scan 2.500E E-03 Flux (#/cm 2 s) 1.500E E E-04 Case 1: MCNP ref Case 4: Sn Pn-Tn S20 t11.2 t11.2 Case 5: hybrid Pn-Tn S20 t11.2 t E Detector Mesh (cm) CM algorithm z (cm) Method Time Speedup MCNP5 (σ<1%) 3510 sec 1.0 TITAN (P N T N S 100 ) 14 sec 249 y (cm) x (cm) Sn algorithm 20
21 A 3 MCNP (Automated Adjoint Accelerated MCNP) (John Wagner and Alireza Haghighat, 1997) A 3 MCNP performs variance reduction technique based on the CADIS (Consistent Adjoint Driven Importance Sampling) methodology by using an approximate deterministic importance function A 3 MCNP input (MCNP input + A 3 MCNP cards) TORT input A 3 MCNP GIP input Step 1 mesh distribution material composition input files TORT GIP multi-group cross sections S N adjoint function A 3 MCNP Step 2 VR parameters A 3 MCNP Step 3 non-analog MC Calculation 21
22 CADIS (Consistent Adjoint Driven Importance Sampling) methodology (John Wagner and Alireza Haghighat, 1997) Uses a 3-D S N importance function distribution for source biasing transport biasing in a consistent manner*, within the weight-window ( r, E) technique. Source biasing Biased source Transport biasing Splitting/rouletting qˆ( p) ( p) q( p) ( p) q( p) dp p ( p) q( p) R ( p) If ( p') <1, particles are processed through the Russian roulette, Otherwise, particles are split 22
23 PWR Cavity dosimetry Application of A 3 MCNP For determination of neutron interaction rates with dosimetry materials placed at the reactor cavity, and estimation of fluence at the reactor pressure vessel BWR Core Shroud Storage cask Determination of neutron and gamma fields at the reactor pressure vessel Determination of neutron and gamma fields at the cask s outside surface 23
24 A 3 MCNP simulation of Storage Cask Objective: Estimation of Neutron and gamma dose on the cask surface Speedup of A 3 MCNP vs. MCNP Problem Size: 180 x 180 x 840 cm 3 A Deep penetration problem 24
25 Remarks Even Fast particle transport codes, with parallel and hybrid algorithms, are not able to obtain solutions in real time! 25
26 Development of Transport Formulations for Real Time Applications Physics Based transport methodologies are needed: Developed Based on problem physics partition a problem into stages (subproblems), For each stage employ response method, adjoint function methodology, or develop/identify a fast algorithm Pre calculate response function or adjoint function using an accurate and fast transport code Solve a linear system of equations to couple all the stages 26
27 Examples for Nondestructive testing: Optimization of the Westinghouse s PGNNA active interrogation system for detection of RCRA (Resource Conversation and Recovery Act) (e.g., lead, mercury, cadmium) in waste drums (partial implementation of MRT; 1999) Nuclear Safeguards: Monitoring of spent fuel pools for detection of fuel diversion (2009) (funded by LLNL) Nuclear nonproliferation: Active interrogation of cargo containers for simulation of special nuclear materials (SNMs) (2013) (in collaboration with GaTech) Spent fuel safety and security: Real time simulation of spent fuel pools for determination of eigenvalue, subcritical multiplication, and material identification (partly funded by I 2 S project, led by GaTech) (Filing for a patent2014) Image reconstruction for SPECT (Single Photon Emission Computed Tomography): Real time simulation of an SPECT device for generation of project images using an MRT methodology and Maximum Likelihood Estimation Maximization (MLEM) (filed for a patent, 2015) 27
28 INSPCT s (William Walters & Alireza Haghighat, 2009) Inspection of a spent fuel pool : Online Calculation of detector response (R): IAEA inspection of Atucha 1 pool A A R n S n n Neutron source Adjoint (Importance) function Source (S = S intrinsic + S subcritical Multiplication ) Stage 1 Intrinsic Source Spontaneous fission & (, n) from fuel burnup calculation (ORIGEN ARP) (Created a database) Stage 2 Subcritical Multiplication (Hybrid method) Simplified fission matrix (FM) method Use MCNP Monte Carlo to obtain a i,j for each pool type (Created a database for coef. a ij ) Multiplication Source Strength (#/source) F i Multiplication Source by Assembly (Fission Matrix) S4 S5S S1 S2S Assembly Position (x,y) N i1 a i, j ( F j S int. j ) Safeguards Adjoint function Stage 3 Is obtained using the PENTRAN transport code (Created a database for multigroup adjoint for different lattice sizes)
29 INSPCT S (Inspection of Nuclear Spent fuel Pool Computing Tool Spreadsheet) INSPCT S checks for fuel diversion; timing <1 s INPUT OUTPUT src file C:\Users\ali\Documents\haghD\ufttg\LLNL\INSPCT-s\se.dsrc COLUMNS 8 fm file C:\Users\ali\Documents\haghD\ufttg\LLNL\INSPCT-s\seResponse Tolerance Detector Normalization ROWS 6 imp file C:\Users\ali\Documents\haghD\ufttg\LLNL\INSPCT-s\se 15.00% 5.28E-10 run Burnup Independent Source (x,y) (x,y) E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E Cooling time Fission Source (x,y) (x,y) E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E Response (experimental) Response(Calculated) (x,y) (x,y) Response Difference (x,y) % % % % %
30 Active Interrogation for Monitoring Special nuclear materials (AIMS) (K. Royston, W. Walters & A. Haghighat, 2013) Detection of SNM in a cargo container Developed a multi stage responsefunction algorithm Gamma Detector HEU g n (x, y, z) F g n (x x 0 ) 2 (y y 0 ) 2 (z z 0 ) 2 ( r, E, ˆ ) S d g n Detector modeling: IFLEX by Georgia Tech Sample Results Computation Time for 5 Scanning Locations AIMS for all HEU positions (serial) MCNP5* (8 processors) HEU position 0or 5 HEU position 1 HEU position 2, 3, or hours 495 hours 254 hours 63 hours 30
31 RAPID (Real time Analysis of spent nuclear fuel Pool for In Situ Detection) (filing for a patent) (W. Walters, A. Haghighat, N. Roskoff, 2014) RAPID uses a MRT methodology based on the Fission matrix technique; it calculates eigenvalue, subcritical multiplication, and pinwise, axially dependent fission density throughout a pool Sample Cases Assembly wise fission density All Case Results *total pre computation time for all FM : 4600 minutes 31 31
32 RAPID Results: 3 D Fission Density Visualization Y LEVEL ANIMATION Z LEVEL ANIMATION 32
33 TITAN-IR (TITAN deterministic Image Reconstruction for SPECT) (Katherine Royston and Alireza Haghighat, 2015) Single Photon Emission Computed Tomography (SPECT): radiopharmaceutical injected into patient 2 D projection images acquired at different angles around the patient Projection images reconstructed to form 3 D radionuclide distribution Accuracy of flux distribution, back of collimator TITAN IR TITAN includes a hybrid formulation with Sn method with fictitious quadrature and ray tracing for modeling the collimator An iterative technique based on Maximumlikelihood expectation maximization (ML EM) for image reconstruction using TITAN for forward projection step Filed for patent (VTIP disclosure number) Reference Solution Projection Angle Detector Element Sinogram Simulation of two collimators Collimator MCNP* TITAN IR* Speed up LEGP 46.9 hrs 132 sec 1353 LEHS 21.4 hrs 14 sec 5503 *on 8 processors Reconstructed Images: Iteration 0 Iteration 1 Iteration 5 Iteration
34 DEMONSTRATION Jaszczak Phantom Reconstruction with TITAN IR (disclosure#: VTIP ) Reconstruction of noiseless SIMIND data with no collimator blur Computation time for 40 iterations Parallel calculation with 16 processors on dedicated computer cluster: 23.6 sec Serial calculation on MacBook Pro:175.6 sec 34
35 Conclusions The MRT methodology (that is a physics based approach) enables real time simulation: This requires careful understanding of the problem physics, and the use of response methods, adjoint function methodology, or development problem dependent formulations; Multiple MRT calculations can be run in parallel for sensitivity studies and uncertainty quantification using parallel computers For pre calculation of response and adjoint functions can be performed in parallel, Advanced codes such as parallel and hybrid Monte Carlo with automated variance reduction (e.g., A 3 MCNP) and parallel and hybrid deterministic (e.g., PENTRAN, TITAN) should be used 35
36 Thanks! Questions? 36
Dosimetry Simulations with the UF-B Series Phantoms using the PENTRAN-MP Code System
Dosimetry Simulations with the UF-B Series Phantoms using the PENTRAN-MP Code System A. Al-Basheer, M. Ghita, G. Sjoden, W. Bolch, C. Lee, and the ALRADS Group Computational Medical Physics Team Nuclear
More informationParallel PENTRAN Applications. G. E. Sjoden and A. Haghighat Nuclear and Radiological Engineering University of Florida
Parallel PENTRAN Applications G. E. Sjoden and A. Haghighat Nuclear and Radiological Engineering University of Florida Overview Introduction Parallel Computing & MPI Boltzmann & Transport PENTRAN TM Code
More informationEvaluation of RAPID for a UNF cask benchmark problem
Evaluation of RAPID for a UNF cask benchmark problem Valerio Mascolino 1,a, Alireza Haghighat 1,b, and Nathan J. Roskoff 1,c 1 Nuclear Science & Engineering Lab (NSEL), Virginia Tech, 900 N Glebe Rd.,
More informationELECTRON DOSE KERNELS TO ACCOUNT FOR SECONDARY PARTICLE TRANSPORT IN DETERMINISTIC SIMULATIONS
Computational Medical Physics Working Group Workshop II, Sep 30 Oct 3, 2007 University of Florida (UF), Gainesville, Florida USA on CD-ROM, American Nuclear Society, LaGrange Park, IL (2007) ELECTRON DOSE
More informationClick to edit Master title style
Fun stuff with the built-in response matrix solver 7th International Serpent UGM, Gainesville, FL, Nov. 6 9, 2017 Jaakko Leppänen VTT Technical Research Center of Finland Click to edit Master title Outline
More informationPSG2 / Serpent a Monte Carlo Reactor Physics Burnup Calculation Code. Jaakko Leppänen
PSG2 / Serpent a Monte Carlo Reactor Physics Burnup Calculation Code Jaakko Leppänen Outline Background History The Serpent code: Neutron tracking Physics and interaction data Burnup calculation Output
More informationDevelopment of a Radiation Shielding Monte Carlo Code: RShieldMC
Development of a Radiation Shielding Monte Carlo Code: RShieldMC Shenshen GAO 1,2, Zhen WU 1,3, Xin WANG 1,2, Rui QIU 1,2, Chunyan LI 1,3, Wei LU 1,2, Junli LI 1,2*, 1.Department of Physics Engineering,
More informationIMPROVEMENTS TO MONK & MCBEND ENABLING COUPLING & THE USE OF MONK CALCULATED ISOTOPIC COMPOSITIONS IN SHIELDING & CRITICALITY
IMPROVEMENTS TO MONK & MCBEND ENABLING COUPLING & THE USE OF MONK CALCULATED ISOTOPIC COMPOSITIONS IN SHIELDING & CRITICALITY N. Davies, M.J. Armishaw, S.D. Richards and G.P.Dobson Serco Technical Consulting
More informationAdaptation of the Nagra Activation Analysis Methodology to Serpent
Adaptation of the Nagra Activation Analysis Methodology to Serpent Valentyn Bykov May 31, 2018, Serpent UGM 2018 Nuclear Power in Switzerland Reactor Type Net First MWe power Beznau 1 PWR 365 1969 Beznau
More informationClick to edit Master title style
New features in Serpent 2 for fusion neutronics 5th International Serpent UGM, Knoxville, TN, Oct. 13-16, 2015 Jaakko Leppänen VTT Technical Research Center of Finland Click to edit Master title Outline
More informationState of the art of Monte Carlo technics for reliable activated waste evaluations
State of the art of Monte Carlo technics for reliable activated waste evaluations Matthieu CULIOLI a*, Nicolas CHAPOUTIER a, Samuel BARBIER a, Sylvain JANSKI b a AREVA NP, 10-12 rue Juliette Récamier,
More informationA premilinary study of the OECD/NEA 3D transport problem using the lattice code DRAGON
A premilinary study of the OECD/NEA 3D transport problem using the lattice code DRAGON Nicolas Martin, Guy Marleau, Alain Hébert Institut de Génie Nucléaire École Polytechnique de Montréal 28 CNS Symposium
More informationRadiological Characterization and Decommissioning of Research and Power Reactors 15602
Radiological Characterization and Decommissioning of Research and Power Reactors 15602 INTRODUCTION Faezeh Abbasi *, Bruno Thomauske *, Rahim Nabbi * RWTH University Aachen The production of the detailed
More informationBreaking Through the Barriers to GPU Accelerated Monte Carlo Particle Transport
Breaking Through the Barriers to GPU Accelerated Monte Carlo Particle Transport GTC 2018 Jeremy Sweezy Scientist Monte Carlo Methods, Codes and Applications Group 3/28/2018 Operated by Los Alamos National
More informationApplication of MCNP Code in Shielding Design for Radioactive Sources
Application of MCNP Code in Shielding Design for Radioactive Sources Ibrahim A. Alrammah Abstract This paper presents three tasks: Task 1 explores: the detected number of as a function of polythene moderator
More informationAutomatic Mesh Adaptivity for Hybrid Monte Carlo/Deterministic Neutronics Modeling of Difficult Shielding Problems. Ahmad Ibrahim
Automatic Mesh Adaptivity for Hybrid Monte Carlo/Deterministic Neutronics Modeling of Difficult Shielding Problems by Ahmad Ibrahim A dissertation submitted in partial fulfillment of the requirement for
More informationABSTRACT. W. T. Urban', L. A. Crotzerl, K. B. Spinney', L. S. Waters', D. K. Parsons', R. J. Cacciapouti2, and R. E. Alcouffel. 1.
COMPARISON OF' THREE-DIMENSIONAL NEUTRON FLUX CALCULATIONS FOR MAINE YANKEE W. T. Urban', L. A. Crotzerl, K. B. Spinney', L. S. Waters', D. K. Parsons', R. J. Cacciapouti2, and R. E. Alcouffel ABSTRACT
More informationDETERMINISTIC 3D RADIATION TRANSPORT SIMULATION FOR DOSE DISTRIBUTION AND ORGAN DOSE EVALUATION IN DIAGNOSTIC CT
DETERMINISTIC 3D RADIATION TRANSPORT SIMULATION FOR DOSE DISTRIBUTION AND ORGAN DOSE EVALUATION IN DIAGNOSTIC CT Monica Ghita,, Glenn Sjoden, Manuel Arreola, Ahmad Al-Basheer Basheer, Choonsik Lee, Wesley
More informationOF THE, AMERICAN NUCLEAR SOCIiZTY
OF THE, AMERICAN NUCLEAR SOCIiZTY June 4-8, 2000 Volume 82 San Diego Town and Country Hotel and Convention Center TANSAO 82 l-300 (2000) San Diego, California ISSN: 0003-O 18X Raymond H. Gabaldon III Technical
More informationElectron Dose Kernels (EDK) for Secondary Particle Transport in Deterministic Simulations
Electron Dose Kernels (EDK) for Secondary Particle Transport in Deterministic Simulations A. Al-Basheer, G. Sjoden, M. Ghita Computational Medical Physics Team Nuclear & Radiological Engineering University
More informationGPU implementation for rapid iterative image reconstruction algorithm
GPU implementation for rapid iterative image reconstruction algorithm and its applications in nuclear medicine Jakub Pietrzak Krzysztof Kacperski Department of Medical Physics, Maria Skłodowska-Curie Memorial
More information1 st International Serpent User Group Meeting in Dresden, Germany, September 15 16, 2011
1 st International Serpent User Group Meeting in Dresden, Germany, September 15 16, 2011 Discussion notes The first international Serpent user group meeting was held at the Helmholtz Zentrum Dresden Rossendorf
More informationDevelopment of a Variance Reduction Scheme in the Serpent 2 Monte Carlo Code Jaakko Leppänen, Tuomas Viitanen, Olli Hyvönen
Development of a Variance Reduction Scheme in the Serpent 2 Monte Carlo Code Jaakko Leppänen, Tuomas Viitanen, Olli Hyvönen VTT Technical Research Centre of Finland, Ltd., P.O Box 1000, FI-02044 VTT, Finland
More informationDRAGON SOLUTIONS FOR BENCHMARK BWR LATTICE CELL PROBLEMS
DRAGON SOLUTIONS FOR BENCHMARK BWR LATTICE CELL PROBLEMS R. Roy and G. Marleau Institut de Génie Nucléaire École Polytechnique de Montréal P.O.Box 6079, Station CV, Montreal, Canada roy@meca.polymtl.ca
More informationMCNP Monte Carlo & Advanced Reactor Simulations. Forrest Brown. NEAMS Reactor Simulation Workshop ANL, 19 May Title: Author(s): Intended for:
LA-UR- 09-03055 Approved for public release; distribution is unlimited. Title: MCNP Monte Carlo & Advanced Reactor Simulations Author(s): Forrest Brown Intended for: NEAMS Reactor Simulation Workshop ANL,
More informationSERPENT Cross Section Generation for the RBWR
SERPENT Cross Section Generation for the RBWR Andrew Hall Thomas Downar 9/19/2012 Outline RBWR Motivation and Design Why use Serpent Cross Sections? Modeling the RBWR Generating an Equilibrium Cycle RBWR
More information2-D Reflector Modelling for VENUS-2 MOX Core Benchmark
2-D Reflector Modelling for VENUS-2 MOX Core Benchmark Dušan Ćalić ZEL-EN d.o.o. Vrbina 18 8270, Krsko, Slovenia dusan.calic@zel-en.si ABSTRACT The choice of the reflector model is an important issue in
More informationNUC E 521. Chapter 6: METHOD OF CHARACTERISTICS
NUC E 521 Chapter 6: METHOD OF CHARACTERISTICS K. Ivanov 206 Reber, 865-0040, kni1@psu.edu Introduction o Spatial three-dimensional (3D) and energy dependent modeling of neutron population in a reactor
More informationMethodology for spatial homogenization in Serpent 2
Methodology for spatial homogenization in erpent 2 Jaakko Leppänen Memo 204/05/26 Background patial homogenization has been one of the main motivations for developing erpent since the beginning of the
More informationVerification of the Hexagonal Ray Tracing Module and the CMFD Acceleration in ntracer
KNS 2017 Autumn Gyeongju Verification of the Hexagonal Ray Tracing Module and the CMFD Acceleration in ntracer October 27, 2017 Seongchan Kim, Changhyun Lim, Young Suk Ban and Han Gyu Joo * Reactor Physics
More informationMichael Speiser, Ph.D.
IMPROVED CT-BASED VOXEL PHANTOM GENERATION FOR MCNP MONTE CARLO Michael Speiser, Ph.D. Department of Radiation Oncology UT Southwestern Medical Center Dallas, TX September 1 st, 2012 CMPWG Workshop Medical
More informationDEVELOPMENT OF METHODS FOR HIGH PERFORMANCE COMPUTING APPLICATIONS OF THE DETERMINISTIC STAGE OF COMET CALCULATIONS
DEVELOPMENT OF METHODS FOR HIGH PERFORMANCE COMPUTING APPLICATIONS OF THE DETERMINISTIC STAGE OF COMET CALCULATIONS A Dissertation Presented to The Academic Faculty by Kyle Eugene Remley In Partial Fulfillment
More informationA PRACTICAL LOOK AT MONTE CARLO VARIANCE REDUCTION METHODS IN RADIATION SHIELDING
A PRACTICAL LOOK AT MONTE CARLO VARIANCE REDUCTION METHODS IN RADIATION SHIELDING RICHARD H. OLSHER Health Physics Measurements Group, Los Alamos National Laboratory MS J573, P.O. Box 1663, Los Alamos,
More informationOutline. Monte Carlo Radiation Transport Modeling Overview (MCNP5/6) Monte Carlo technique: Example. Monte Carlo technique: Introduction
Monte Carlo Radiation Transport Modeling Overview () Lecture 7 Special Topics: Device Modeling Outline Principles of Monte Carlo modeling Radiation transport modeling with Utilizing Visual Editor (VisEd)
More informationMCNP Variance Reduction technique application for the Development of the Citrusdal Irradiation Facility
IYNC 2008 Interlaken, Switzerland, 20-26 September 2008 Paper No. 376 MCNP Variance Reduction technique application for the Development of the Citrusdal Irradiation Facility R Makgae Pebble Bed Modular
More informationModeling the ORTEC EX-100 Detector using MCNP
Modeling the ORTEC EX-100 Detector using MCNP MCNP is a general-purpose Monte Carlo radiation transport code for modeling the interaction of radiation with materials based on composition and density. MCNP
More informationAttila4MC. Software for Simplifying Monte Carlo. For more info contact or
Attila4MC Software for Simplifying Monte Carlo For more info contact attila@varian.com or Gregory.Failla@varian.com MCNP and MCNP6 are trademarks of Los Alamos National Security, LLC, Los Alamos National
More informationModeling the White Sands Missile Range Fast Burst Reactor Using a Discrete Ordinates Code, PENTRAN
Air Force Institute of Technology AFIT Scholar Theses and Dissertations 3-23-207 Modeling the White Sands Missile Range Fast Burst Reactor Using a Discrete Ordinates Code, PENTRAN Taylor R. Schulmeister
More informationA Method for Estimating Criticality Lower Limit Multiplication Factor. Yoshitaka NAITO NAIS Co., Ltd.
A Method for Estimating Criticality Lower Limit Multiplication Factor Yoshitaka NAITO NAIS Co., Ltd. Progress Grade up of computer performance Sub-criticality becomes to be decided using computed results
More informationStatus of the Serpent criticality safety validation package
VTT TECHNICAL RESEARCH CENTRE OF FINLAND LTD Status of the Serpent criticality safety validation package Serpent UGM 2017 Riku Tuominen and Ville Valtavirta, VTT Outline Criticality Safety Evaluation What
More informationDaedeok-daero, Yuseong-gu, Daejeon , Republic of Korea b Argonne National Laboratory (ANL)
MC 2-3/TWODANT/DIF3D Analysis for the ZPPR-15 10 B(n, α) Reaction Rate Measurement Min Jae Lee a*, Donny Hartanto a, Sang Ji Kim a, and Changho Lee b a Korea Atomic Energy Research Institute (KAERI) 989-111
More informationValidation of GEANT4 for Accurate Modeling of 111 In SPECT Acquisition
Validation of GEANT4 for Accurate Modeling of 111 In SPECT Acquisition Bernd Schweizer, Andreas Goedicke Philips Technology Research Laboratories, Aachen, Germany bernd.schweizer@philips.com Abstract.
More informationMonte-Carlo-Based Scatter Correction for Quantitative SPECT Reconstruction
Monte-Carlo-Based Scatter Correction for Quantitative SPECT Reconstruction Realization and Evaluation Rolf Bippus 1, Andreas Goedicke 1, Henrik Botterweck 2 1 Philips Research Laboratories, Aachen 2 Fachhochschule
More informationIMPROVING COMPUTATIONAL EFFICIENCY OF MONTE-CARLO SIMULATIONS WITH VARIANCE REDUCTION
International Conference on Mathematics and Computational Methods Applied to Nuclear Science & Engineering (M&C 2013), Sun Valley, Idaho, USA, May 5-9, 2013, on CD-ROM, American Nuclear Society, LaGrange
More informationMesh Human Phantoms with MCNP
LAUR-12-01659 Mesh Human Phantoms with MCNP Casey Anderson (casey_a@lanl.gov) Karen Kelley, Tim Goorley Los Alamos National Laboratory U N C L A S S I F I E D Slide 1 Summary Monte Carlo for Radiation
More informationParticle track plotting in Visual MCNP6 Randy Schwarz 1,*
Particle track plotting in Visual MCNP6 Randy Schwarz 1,* 1 Visual Editor Consultants, PO Box 1308, Richland, WA 99352, USA Abstract. A visual interface for MCNP6 has been created to allow the plotting
More informationTRANSX-2005 New Structure and Features R.E.MacFarlane Los Alamos National Laboratory
TRANSX-2005 New Structure and Features R.E.MacFarlane Los Alamos National Laboratory TRANSX-2005 is a translation of TRANSX to Fortran- 90/95 style with an extended code-management scheme. The new features
More informationEvaluation of RayXpert for shielding design of medical facilities
Evaluation of Raypert for shielding design of medical facilities Sylvie Derreumaux 1,*, Sophie Vecchiola 1, Thomas Geoffray 2, and Cécile Etard 1 1 Institut for radiation protection and nuclear safety,
More information3D DETERMINISTIC RADIATION TRANSPORT FOR DOSE COMPUTATIONS IN CLINICAL PROCEDURES
3D DETERMINISTIC RADIATION TRANSPORT FOR DOSE COMPUTATIONS IN CLINICAL PROCEDURES By AHMAD AL-BASHEER A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT
More informationVerification of the 3D Method of characteristics solver in OpenMOC
Verification of the 3D Method of characteristics solver in OpenMOC The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published
More informationDEVELOPMENT OF METHODS AND TOOLS FOR ON-THE-FLY RESPONSE FUNCTION GENERATION FOR CRITICALITY CALCULATIONS
DEVELOPMENT OF METHODS AND TOOLS FOR ON-THE-FLY RESPONSE FUNCTION GENERATION FOR CRITICALITY CALCULATIONS A Dissertation Presented to The Academic Faculty by Ryan Paul Hon In Partial Fulfillment of the
More informationIntroduction to Positron Emission Tomography
Planar and SPECT Cameras Summary Introduction to Positron Emission Tomography, Ph.D. Nuclear Medicine Basic Science Lectures srbowen@uw.edu System components: Collimator Detector Electronics Collimator
More informationStatus and development of multi-physics capabilities in Serpent 2
Status and development of multi-physics capabilities in Serpent 2 V. Valtavirta VTT Technical Research Centre of Finland ville.valtavirta@vtt.fi 2014 Serpent User Group Meeting Structure Click to of edit
More informationMultiphysics simulations of nuclear reactors and more
Multiphysics simulations of nuclear reactors and more Gothenburg Region OpenFOAM User Group Meeting Klas Jareteg klasjareteg@chalmersse Division of Nuclear Engineering Department of Applied Physics Chalmers
More informationA COARSE MESH RADIATION TRANSPORT METHOD FOR PRISMATIC BLOCK THERMAL REACTORS IN TWO DIMENSIONS
A COARSE MESH RADIATION TRANSPORT METHOD FOR PRISMATIC BLOCK THERMAL REACTORS IN TWO DIMENSIONS A Thesis Presented to The Academic Faculty By Kevin John Connolly In Partial Fulfillment Of the Requirements
More informationBasics of treatment planning II
Basics of treatment planning II Sastry Vedam PhD DABR Introduction to Medical Physics III: Therapy Spring 2015 Monte Carlo Methods 1 Monte Carlo! Most accurate at predicting dose distributions! Based on
More informationTeresa S. Bailey (LLNL) Marvin L. Adams (Texas A&M University)
A Piecewise i Linear Discontinuous Finite it Element Spatial Discretization of the S N Transport Equation for Polyhedral Grids in 3D Cartesian Geometry International Conference on Transport Theory September
More informationClick to edit Master title style
Introduction to Serpent Code Fusion neutronics workshop, Cambridge, UK, June 11-12, 2015 Jaakko Leppänen VTT Technical Research Center of Finland Click to edit Master title Outline style Serpent overview
More informationDose Calculations: Where and How to Calculate Dose. Allen Holder Trinity University.
Dose Calculations: Where and How to Calculate Dose Trinity University www.trinity.edu/aholder R. Acosta, W. Brick, A. Hanna, D. Lara, G. McQuilen, D. Nevin, P. Uhlig and B. Slater Dose Calculations - Why
More informationIdentification of Shielding Material Configurations Using NMIS Imaging
Identification of Shielding Material Configurations Using NMIS Imaging B. R. Grogan, J. T. Mihalczo, S. M. McConchie, and J. A. Mullens Oak Ridge National Laboratory, P.O. Box 2008, MS-6010, Oak Ridge,
More informationIndrin Chetty Henry Ford Hospital Detroit, MI. AAPM Annual Meeting Houston 7:30-8:25 Mon 08/07/28 1/30
Review of TG105: Issues associated with clinical implementation of Monte Carlo-based photon and electron external beam treatment planning D. W. O. Rogers, Carleton Laboratory for Radiotherapy Physics,
More informationA dedicated tool for PET scanner simulations using FLUKA
A dedicated tool for PET scanner simulations using FLUKA P. G. Ortega FLUKA meeting June 2013 1 Need for in-vivo treatment monitoring Particles: The good thing is that they stop... Tumour Normal tissue/organ
More informationImage Acquisition Systems
Image Acquisition Systems Goals and Terminology Conventional Radiography Axial Tomography Computer Axial Tomography (CAT) Magnetic Resonance Imaging (MRI) PET, SPECT Ultrasound Microscopy Imaging ITCS
More informationGraphical User Interface for High Energy Multi-Particle Transport
Graphical User Interface for High Energy Multi-Particle Transport Phase I Final Report PREPARED BY: P.O. Box 1308 Richland, WA 99352-1308 PHONE: (509) 539-8621 FAX: (509) 946-2001 Email: randyschwarz@mcnpvised.com
More informationDetermination of Three-Dimensional Voxel Sensitivity for Two- and Three-Headed Coincidence Imaging
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 50, NO. 3, JUNE 2003 405 Determination of Three-Dimensional Voxel Sensitivity for Two- and Three-Headed Coincidence Imaging Edward J. Soares, Kevin W. Germino,
More informationChallenges and developments in fusion neutronics a CCFE perspective
Challenges and developments in fusion neutronics a CCFE perspective A. Turner Applied Radiation Physics group SERPENT fusion neutronics workshop, Cambridge, June 2015 CCFE is the fusion research arm of
More informationAutomated ADVANTG Variance Reduction in a Proton Driven System. Kenneth A. Van Riper1 and Robert L. Metzger2
Automated ADVANTG Variance Reduction in a Proton Driven System Kenneth A. Van Riper1 and Robert L. Metzger2 1 White Rock Science, P. O. Box 4729, White Rock, NM 87547, kvr@rt66.com Radiation Safety Engineering,
More informationReconstruction from Projections
Reconstruction from Projections M.C. Villa Uriol Computational Imaging Lab email: cruz.villa@upf.edu web: http://www.cilab.upf.edu Based on SPECT reconstruction Martin Šámal Charles University Prague,
More informationMonte Carlo methods in proton beam radiation therapy. Harald Paganetti
Monte Carlo methods in proton beam radiation therapy Harald Paganetti Introduction: Proton Physics Electromagnetic energy loss of protons Distal distribution Dose [%] 120 100 80 60 40 p e p Ionization
More informationCHAPTER 10: TALLYING IN MCNP
_or_.e_sa_m_h_us_se_in 6:..:...7 ---"M.=o:.c.;;nte.:-C-"-=arlo Particle Transport with MCNP CHAPTER 10: TALLYING IN MCNP Tallying is the process of scoring the parameters of interest, Le. providing the
More informationTheoretical Investigations of Tomographic Methods used for Determination of the Integrity of Spent BWR Nuclear Fuel
a UPPSALA UNIVERSITY Department of Radiation Sciences Box 535, S-751 1 Uppsala, Sweden http://www.tsl.uu.se/ Internal report ISV-6/97 August 1996 Theoretical Investigations of Tomographic Methods used
More informationReview of PET Physics. Timothy Turkington, Ph.D. Radiology and Medical Physics Duke University Durham, North Carolina, USA
Review of PET Physics Timothy Turkington, Ph.D. Radiology and Medical Physics Duke University Durham, North Carolina, USA Chart of Nuclides Z (protons) N (number of neutrons) Nuclear Data Evaluation Lab.
More informationShielding factors for traditional safety glasses
Shielding factors for traditional safety glasses Malcolm McEwen, Hong Shen and Ernesto Mainegra-Hing Ionizing Radiation Standards, National Research Council Canada Alan DuSautoy, Radiation and Health Sciences
More informationThe quest for a high performance Boltzmann transport solver
The quest for a high performance Boltzmann transport solver P. N. Brown, B. Chang, U. R. Hanebutte & C. S. Woodward Center/or Applied Scientific Computing, Lawrence Livermore National Laboratory, USA.
More informationNON-COLLIMATED SCATTERED RADIATION TOMOGRAPHY
NON-COLLIMATED SCATTERED RADIATION TOMOGRAPHY Gorshkov V.A., Space Research Institute, Moscow, Russia Yumashev V.M., State corporation "Rosatom", Centre "Atom-innovation", Moscow, Russia Kirilenko K.V.,
More information2017 Summer Course on Optical Oceanography and Ocean Color Remote Sensing. Monte Carlo Simulation
2017 Summer Course on Optical Oceanography and Ocean Color Remote Sensing Curtis Mobley Monte Carlo Simulation Delivered at the Darling Marine Center, University of Maine July 2017 Copyright 2017 by Curtis
More informationContinuation Format Page
C.1 PET with submillimeter spatial resolution Figure 2 shows two views of the high resolution PET experimental setup used to acquire preliminary data [92]. The mechanics of the proposed system are similar
More informationSimulations in emission tomography using GATE
Simulations in emission tomography using GATE Irène Buvat buvat@imed.jussieu.fr Laboratory of Functional Imaging, U678 INSERM, Paris, France Outline Emission tomography and need for simulations GATE short
More informationNeutronics Analysis of TRIGA Mark II Research Reactor. R. Khan, S. Karimzadeh, H. Böck Vienna University of Technology Atominstitute
Neutronics Analysis of TRIGA Mark II Research Reactor R. Khan, S. Karimzadeh, H. Böck Vienna University of Technology Atominstitute 23-03-2010 TRIGA Mark II reactor MCNP radiation transport code MCNP model
More informationMathematical methods and simulations tools useful in medical radiation physics
Mathematical methods and simulations tools useful in medical radiation physics Michael Ljungberg, professor Department of Medical Radiation Physics Lund University SE-221 85 Lund, Sweden Major topic 1:
More informationI. INTRODUCTION. Figure 1. Radiation room model at Dongnai General Hospital
International Journal of Computational Engineering Research Vol, 04 Issue, 4 Simulation of Photon and Electron dose distributions 5 code for the treatment area using the linear electron accelerator (LINAC)
More informationVirtual EM Inc. Ann Arbor, Michigan, USA
Functional Description of the Architecture of a Special Purpose Processor for Orders of Magnitude Reduction in Run Time in Computational Electromagnetics Tayfun Özdemir Virtual EM Inc. Ann Arbor, Michigan,
More informationThe Monte Carlo simulation of a Package formed by the combination of three scintillators: Brillance380, Brillance350, and Prelude420.
EURONS I3 506065 JRA9 RHIB Report made during stay IEM-CSIC Madrid december 2006 MINISTERIO DE ASUNTOS EXTERIORES Y DE COOPERACIÓN AECI VICESECRETARÍA GENERAL The Monte Carlo simulation of a Package formed
More information- PARALLEL PROCESSING AND VISUALIZATION
RECENT DEVELOPMENT FOR THE ITS CODE SYSTEM - PARALLEL PROCESSING AND VISUALIZATION Wesley C. Fan C. David Turner John A. Halbleib, Sr. and Ronald P. Kensek Sandia National Laboratories P. 0. BOX 5800 Albuquerque,
More informationParallel 3D Sweep Kernel with PaRSEC
Parallel 3D Sweep Kernel with PaRSEC Salli Moustafa Mathieu Faverge Laurent Plagne Pierre Ramet 1 st International Workshop on HPC-CFD in Energy/Transport Domains August 22, 2014 Overview 1. Cartesian
More informationCPM-3 BENCHMARKING to the DOE/B&W CRITICAL EXPERIMENTS
CPM-3 BENCHMARKING to the DOE/B&W CRITICAL EXPERIMENTS Kenneth M. Smolinske and Rodney L. Grow Utility Resource Associates Corporation 1901 Research Boulevard, Suite 405 Rockville, Maryland 20850 ABSTRACT
More informationCALCULATION OF THE ACTIVITY INVENTORY FOR THE TRIGA REACTOR AT THE MEDICAL UNIVERSITY OF HANNOVER (MHH) IN PREPARATION FOR DISMANTLING THE FACILITY
CALCULATION OF THE ACTIVITY INVENTORY FOR THE TRIGA REACTOR AT THE MEDICAL UNIVERSITY OF HANNOVER (MHH) IN PREPARATION FOR DISMANTLING THE FACILITY Gabriele Hampel, Friedemann Scheller, Medical University
More informationBEAVRS benchmark calculations with Serpent-ARES code sequence
BEAVRS benchmark calculations with Serpent-ARES code sequence Jaakko Leppänen rd International Serpent User Group Meeting Berkeley, CA, Nov. 6-8, Outline Goal of the study The ARES nodal diffusion code
More informationCustomizable and Advanced Software for Tomographic Reconstruction
Customizable and Advanced Software for Tomographic Reconstruction 1 What is CASToR? Open source toolkit for 4D emission (PET/SPECT) and transmission (CT) tomographic reconstruction Focus on generic, modular
More informationReducing 3D MOC Storage Requirements with Axial Onthe-fly
Reducing 3D MOC Storage Requirements with Axial Onthe-fly Ray Tracing The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As
More informationOPTIMIZATION OF MONTE CARLO TRANSPORT SIMULATIONS IN STOCHASTIC MEDIA
PHYSOR 2012 Advances in Reactor Physics Linking Research, Industry, and Education Knoxville, Tennessee, USA, April 15-20, 2012, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2010) OPTIMIZATION
More information3/27/2012 WHY SPECT / CT? SPECT / CT Basic Principles. Advantages of SPECT. Advantages of CT. Dr John C. Dickson, Principal Physicist UCLH
3/27/212 Advantages of SPECT SPECT / CT Basic Principles Dr John C. Dickson, Principal Physicist UCLH Institute of Nuclear Medicine, University College London Hospitals and University College London john.dickson@uclh.nhs.uk
More informationarxiv: v2 [cond-mat.mtrl-sci] 5 Jan 2010
Gamma scattering scanning of concrete block for detection of voids. Shivaramu 1, Arijit Bose 2 and M. Margret 1 1 Radiological Safety Division, Safety Group, IGCAR, Kalpakaam - 63 12 (India) 2 Chennai
More informationA New Method for Coupling 2D and 3D Deterministic and Stochastic Radiation Transport Calculations
University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Masters Theses Graduate School 8-2011 A New Method for Coupling 2D and 3D Deterministic and Stochastic Radiation Transport
More informationComparison of Shutdown Dose Rate Results using MCNP6 Activation Capability and MCR2S
APPLIED RADIATION PHYSICS GROUP TECHNICAL NOTE ARP-097 July 2014 Comparison of Shutdown Dose Rate Results using MCNP6 Activation Capability and MCR2S A. Turner 1, Z. Ghani 1, J. Shimwell 2 1: CCFE, Culham
More informationThe Spherical Harmonics Discrete Ordinate Method for Atmospheric Radiative Transfer
The Spherical Harmonics Discrete Ordinate Method for Atmospheric Radiative Transfer K. Franklin Evans Program in Atmospheric and Oceanic Sciences University of Colorado, Boulder Computational Methods in
More informationConstructing System Matrices for SPECT Simulations and Reconstructions
Constructing System Matrices for SPECT Simulations and Reconstructions Nirantha Balagopal April 28th, 2017 M.S. Report The University of Arizona College of Optical Sciences 1 Acknowledgement I would like
More informationANALYSIS OF CT AND PET/SPECT IMAGES FOR DOSIMETRY CALCULATION
2009 International Nuclear Atlantic Conference - INAC 2009 Rio de Janeiro,RJ, Brazil, September27 to October 2, 2009 ASSOCIAÇÃO BRASILEIRA DE ENERGIA NUCLEAR - ABEN ISBN: 978-85-99141-03-8 ANALYSIS OF
More informationParallel computation performances of Serpent and Serpent 2 on KTH Parallel Dator Centrum
KTH ROYAL INSTITUTE OF TECHNOLOGY, SH2704, 9 MAY 2018 1 Parallel computation performances of Serpent and Serpent 2 on KTH Parallel Dator Centrum Belle Andrea, Pourcelot Gregoire Abstract The aim of this
More informationRobin P. Gardner* and. Lianyan Liu. Computalog U.S.A., Inc., 501 Winscott Road, Fort Worth, Texas 76126
NUCLEAR SCIENCE AND ENGINEERING: 133, 80 91 ~1999! Monte Carlo Simulation of Neutron Porosity Oil Well Logging Tools: Combining the Geometry-Independent Fine-Mesh Importance Map and One-Dimensional Diffusion
More information