Monte Carlo Methods & Virtual Photonics
|
|
- Herbert Hudson
- 6 years ago
- Views:
Transcription
1 Monte Carlo Methods & Virtual Photonics Tuesday-Thursday 9 10:30 Room 3201, Nat. Sci. II Jerry Spanier P220 BLI jspanier@uci.edu Please help yourselves to refreshments at the back of the room and leave your name, affiliation and contact information on the sign-in sheet. L A M M P
2 Monte Carlo Methods & Virtual Photonics Purposes Provide basic instruction in MC methods for BLI, CCBS and MCSB Integrate M&C across LAMMP projects Advance plans for LAMMP Workshops on VTS/ATKs L A M M P
3 Specific Goals Provide overview of conventional + more advanced MC methods Outline conceptual framework for MC algorithm development from user perspective (VTS/ATK) Describe features of newer adaptive MC methods that converge geometrically L A M M P
4 Detecting Cancer Non-invasively At BLI we use a number of optical technologies to detect, diagnose and treat various subtle transformations in tissue, including those caused by cancer, wounds/burns, etc. The effective use of these methods requires accurate models of the interactions of light with tissue. L A M M P
5 L A M M P Tissue Pathology
6 L A M M P Histology: Squamous Epithelium Note irregular pattern
7 L A M M P Histology: Breast Carcinoma More irregularities, dark features
8 L A M M P BLI Optical Technologies MI
9 L A M M P Mathematical Models
10 Background Radiative transport theory provides a widely accepted gold standard for validation in biomedical optics For general geometries, fully heterogeneous descriptions of turbid media, such as tissue, Monte Carlo has unique capabilities Naive simulation of RTE is rarely effective (more sophisticated, not necessarily intuitive methods are needed) L A M M P
11 Handouts Brief reading list Outline of topics covered L A M M P
12 L A M M P Reading List 1. The Monte Carlo Method in Atmospheric Optics by G.I. Marchuk et. al., Springer Series in Optical Sciences, Springer- Verlag, Neutron Transport Theory by B. Davison, Oxford Univ. Press, Nuclear Reactor Theory by G.I. Bell & S. Glasstone, Krieger Pub. Co., Wave Propagation and Scattering in Random Media by A. Ishimaru, Academic Press 1978, IEEE Press, Monte Carlo Principles and Neutron Transport Problems by J. Spanier and E.M. Gelbard, original ed. 1969, reprint ed. Dover, 2008 (to appear) 6. Monte Carlo Particle Transport Methods: Neutron and Photon Calculations, I. Lux and L. Koblinger, CRC Press, 1991
13 Short Course Topics Outline Monte Carlo Basics Equations of Radiation Transport Diffusion & Higher Order Approximations Keys to Efficient Use of Monte Carlo The Virtual Tissue Model Adaptive Monte Carlo Methods L A M M P
14 L A M M P Monte Carlo Modeling Natural Transport Applications Atmospheric Optics Nuclear Reactor Design/Engineering Biomedical Optics Radiation Therapy Planning Financial Modeling Important Numerical Applications Estimation of Integrals Solution of Matrix Equations (e.g., linear systems arising from discretization of PDEs and/or integral equations Solution of Integral Equations The MC method is very robust and can be used for phenomenological descriptions that incorporate new knowledge quite easily. Monte Carlo simulations account for well over 50% of all computer time worldwide. Optimizing their use is important.
15 The Virtual Tissue System: First Look Visualization Optimization Analysis Solver Problem Definition L A M M P
16 L A M M P ATK/VTS Relationships
17 Tissue Primitives : Towards a Virtual Tissue Model L A M M P
18 Why Monte Carlo? MC is a versatile method for producing answers by computer using random sampling Estimate π: 4 (14/18) = 3.11 L A M M P
19 L A M M P Estimate Integrals The same principle can be used to estimate I b = a f( x) dx and higher order integrals f(x) a b
20 L A M M P Simple Random Walk What is average number of steps (each uniform on [0,1]) to pass through unit slab?
21 L A M M P The Exact Solution Let X = # steps to escape. X = n means (u 1,,u n-1 )εt n-1 &(u 1,,u n ) ε T n so Pr{X = n} = 1/(n-1)! 1/n! = 1/n(n-2)! = f(n) E[X] = nf ( n) = e = n= 2 and σ = T n
22 Model Problem Characteristics Solution is an infinite sum. The MC solution does not truncate this sum! Instead, it picks out every term with correct probability by sampling. Thus, the only error is statistical; i.e., because # of samples is finite. These qualities make the MC method quite special. L A M M P
23 L A M M P Monte Carlo Estimate Tchebycheff s inequality: Pr{ estimate(n) e > kσ(n)} 1/k 2 n 3.84(100) 2 (.77/2.72) 2 = 3100 to achieve relative error ~ 1% with probability.95. n = 1000 gives est(1000) = 2.737
24 Doing Better: Importance Sampling Exaggerate large value of X (sample from g(n) rather than f(n) based on uniform steps but count nf(n)/g(n) instead of n for X. Try h(u) = const. exp(-ku) instead of h(u) = 1. With n steps chosen this way, correct weight factor is [h(u 1 )h(u 2 ) h(u n )] -1 which is 4 times as efficient as before. n = 1000 gave est(1000) = L A M M P
25 Random Walk Processes Consider the general (operator) equation (1) φ = Ҝφ + s where Ҝ < 1 (or even Ҝ m < 1 for some integer m). Then under rather general circumstances, (1) has a unique solution φ = (I Ҝ) -1 s = (I + Ҝ + Ҝ 2 + )s. These simple ideas are basic and relevant to what we do later. L A M M P
26 L A M M P Two Important Special Cases If Ҝ is an nxn matrix, H, and φ, s are n- dimensional vectors, (1) becomes x = H x + s. MC methods can be applied to many important matrix problems. Okten has recently shown that such methods compete very favorably with the best deterministic matrix solvers; e.g., based on conjugate gradient methods.
27 L A M M P Okten s Findings From G. Okten, SIAM J. Sci. Comp. 27, pp : For 0.9 hij = ; bi = i/ N N+ i+ j time to converge is
28 More from Okten The MC method uses sequential correlated sampling (discussed in final class) to achieve geometric convergence in a few iterations. Table 7.4 shows the norm of the rel. residual and # adaptive stages to converge. L A M M P
29 L A M M P Table 7.4
30 L A M M P More About Matrix Problems The study of MC methods for solving matrix problems introduces most of the important ideas and methods used in solving continuous problems. The main difference is that the phase space for matrix problems consists of the integers {1, N}, N = order of matrix, while for continous problems, the phase space is continuous and multidimensional, in general.
31 Linear Integral Equations When Ҝ is an integral operator, K, and ψ, S are continuous functions, (1) becomes Ψ(P) = (K ψ)(p) + S(P) = K(P,Q) Ψ(Q) dq + S(P) which provides a model for all of our RTE applications. L A M M P
32 Duality and Reciprocity When solving (1) φ = Ҝφ + s, often we need only an accurate value of one or more functionals <s*, φ>, where (2) φ* = Ҝ*φ* + s*. It follows easily from the definitions of the operator adjoint operation * that (3) <s*, φ> = < φ*, s>. L A M M P
33 L A M M P Important Connections When we estimate either <s*, φ> or < φ *, s> by MC, we do so in the context of a probability space Р and an estimator (random variable) ξ: Р R so that E[ξ] = ξ dμ = <s*, φ> and Var [ξ] is as small as possible. Here, Р consists of all random walks on the phase space Γ (which is {1,..., N} for matrix problems or part of R d for d-dimensional integral equations).
34 L A M M P Further Connections Then φ = Ҝφ + s gives rise to a natural random walk process {s, Ҝ} where s (suitably normalized) is used to create initial collision points in Γ, and Ҝ (suitably normalized) is used to move from state to state in Γ as well as to terminate the random walks (in a finite number of steps with probability 1 if Ҝ m < 1).
35 Still More Connections Now define P n (x) = Pr density of making collision n at x for n = 1,2, Then P n (x) = [Ҝ P n-1 ](x), n > 1 and P 1 (x) = Pr{x 1 = x} = s(x). Define a collision estimator on the space Р of random walks on Γ by X(x) = number density of collisions at x εγ L A M M P
36 Completing the Connections Then if an abstract collision density is defined by P(x) = E[X(x)], we find that (*) P(x) = s(x) + [Ҝ P](x). The solution of (*) is always accomplished in practice through the Neumann series representation (I Ҝ) -1 s = (I + Ҝ + Ҝ 2 + )s based on the assumption that Ҝ m < 1 for some integer m. L A M M P
37 L A M M P The Final Connections Finally, the entire theory can be carried through for importance sampling; i.e., based on replacing the analog source s by ŝ and the operator Ҝ by ˆK under quite general conditions. The validity of importance sampling is established through. E[ξ] = ξ dμ = [ξdμ/dν] dν and the generalized derivative dμ/dν serves to correct the otherwise biased results.
38 Monte Carlo Methods and Virtual Photonics Acknowledgements: For co-sponsorship: Bruce Tromberg (BLI and LAMMP Director) and Arthur Lander (CCSB Director & Chair, Dev. & Cell Biology) NIH P-41-RR01192 (LAMMP), NSF/DMS , UCOP (LANL) The core Virtual Photonics group : (Vasan Venugopalan, Carole Hayakawa, Katya Bhan, David Cuccia, Albert Cerussi, Tony Durkin) Rong Kong, Martin Ambrose and the entire Friday morning M&C gang L A M M P
NIH Public Access Author Manuscript J Comput Phys. Author manuscript; available in PMC 2009 September 22.
NIH Public Access Author Manuscript Published in final edited form as: J Comput Phys. 2008 December 1; 227(23): 9762 9777. doi:10.1016/j.jcp.2008.07.016. A new proof of geometric convergence for general
More informationSTOCHASTIC & DETERMINISTIC SOLVERS
STOCHASTIC & DETERMINISTIC SOLVERS Outline Spatial Scales of Optical Technologies & Mathematical Models Prototype RTE Problems 1-D Transport in Slab Geometry: Exact Solution Stochastic Models: Monte Carlo
More informationAdaptive Monte Carlo Algorithms Applied to Heterogeneous Transport Problems
Adaptive Monte Carlo Algorithms Applied to Heterogeneous Transport Problems Katherine Bhan, Rong Kong, and Jerome Spanier Abstract We apply three generations of geometrically convergent adaptive Monte
More informationMonte Carlo Methods:
Short Course on Computational Monte Carlo Methods: Fundamentals Shuang Zhao Assistant Professor Computer Science Department University of California, Irvine Shuang Zhao 1 Teaching Objective Introducing
More informationDiffuse Optical Spectroscopy: Analysis of BiophotonicsSignals Short Course in Computational Biophotonics Albert Cerussi Day 4/Lecture #2
Diffuse Optical Spectroscopy: Analysis of BiophotonicsSignals 2013 Short Course in Computational Biophotonics Albert Cerussi Day 4/Lecture #2 Disclosures No financial conflicts to disclose Learning Objectives
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 1 (2/3/04) Overview -- Interactions, Distributions, Cross Sections, Applications
.54 Neutron Interactions and Applications (Spring 004) Chapter 1 (/3/04) Overview -- Interactions, Distributions, Cross Sections, Applications There are many references in the vast literature on nuclear
More informationModel Counting for Logical Theories
Model Counting for Logical Theories Wednesday Dmitry Chistikov Rayna Dimitrova Department of Computer Science University of Oxford, UK Max Planck Institute for Software Systems (MPI-SWS) Kaiserslautern
More informationMCRT: L4 A Monte Carlo Scattering Code
MCRT: L4 A Monte Carlo Scattering Code Plane parallel scattering slab Optical depths & physical distances Emergent flux & intensity Internal intensity moments Constant density slab, vertical optical depth
More informationMonte Carlo method projective estimators for angular and temporal characteristics evaluation of polarized radiation
Monte Carlo method projective estimators for angular and temporal characteristics evaluation of polarized radiation Natalya Tracheva Sergey Ukhinov Institute of Computational Mathematics and Mathematical
More informationRadiative transfer equations with varying refractive index: a mathematical perspective
Radiative transfer equations with varying refractive index: a mathematical perspective Guillaume Bal November 23, 2005 Abstract This paper reviews established mathematical techniques to model the energy
More informationDetectors in Nuclear Physics: Monte Carlo Methods. Dr. Andrea Mairani. Lectures I-II
Detectors in Nuclear Physics: Monte Carlo Methods Dr. Andrea Mairani Lectures I-II INTRODUCTION Sampling from a probability distribution Sampling from a probability distribution X λ Sampling from a probability
More informationNATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS. (Semester I: AY ) Time Allowed: 2 Hours
NATIONAL UNIVERSITY OF SINGAPORE PC55 NUMERICAL RECIPES WITH APPLICATIONS (Semester I: AY 0-) Time Allowed: Hours INSTRUCTIONS TO CANDIDATES. This examination paper contains FIVE questions and comprises
More informationReduction of Variance. Importance Sampling
Reduction of Variance As we discussed earlier, the statistical error goes as: error = sqrt(variance/computer time). DEFINE: Efficiency = = 1/vT v = error of mean and T = total CPU time How can you make
More informationMultiple-source optical diffusion approximation for a multilayer scattering medium
Multiple-source optical diffusion approximation for a multilayer scattering medium Joseph. Hollmann 1 and ihong V. Wang 1,2, * 1 Optical Imaging aboratory, Department of Biomedical Engineering, Texas A&M
More informationMonte Carlo Radiation Transport Kenny Wood
Monte Carlo Radiation Transport Kenny Wood kw25@st-andrews.ac.uk A practical approach to the numerical simulation of radiation transport Develop programs for the random walks of photons and neutrons using
More informationReduction of Random Variables in Structural Reliability Analysis
Reduction of Random Variables in Structural Reliability Analysis S. Adhikari and R. S. Langley Department of Engineering University of Cambridge Trumpington Street Cambridge CB2 1PZ (U.K.) February 21,
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationMCRT L9: Neutron Transport
MCRT L9: Neutron Transport Outline of code for mono-energetic particles, if change particle energy must change cross sections Bank locations of fission neutrons for subsequent generations Criticality calculations,
More informationOutline. Indrin J. Chetty, AAPM 2006 Monte Carlo CE course. Indrin J. Chetty Henry Ford Hospital. David W. O. Rogers Carleton University
AAPM Task Group Report No. 105: Issues associated with clinical implementation of Monte Carlo-based photon and electron external beam treatment planning Indrin J. Chetty Henry Ford Hospital David W. O.
More informationTechnical Meeting on Modern Neutron Detection
Technical Meeting (F1-TM-55243) Vienna International Center Room M6 (Please note change of room compared to Note Verbale) Technical Meeting on Modern Neutron Detection Monday 4 September 12:30 Welcome
More informationInferring from data. Theory of estimators
Inferring from data Theory of estimators 1 Estimators Estimator is any function of the data e(x) used to provide an estimate ( a measurement ) of an unknown parameter. Because estimators are functions
More informationMonte Carlo Radiation Transfer I
Monte Carlo Radiation Transfer I Monte Carlo Photons and interactions Sampling from probability distributions Optical depths, isotropic emission, scattering Monte Carlo Basics Emit energy packet, hereafter
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationVARIATONAL STATIONARY FUNCTIONALS IN MONTE CARLO COMPUTATIONS
Chapter 25 VARIATONAL STATIONARY FUNCTIONALS IN MONTE CARLO COMPUTATIONS M. Ragheb 2/0/203 INTRODUCTION Use of stationary functionals from variational theory in Monte Carlo computations is discussed. Minimum
More informationA Backward Particle Interpretation of Feynman-Kac Formulae
A Backward Particle Interpretation of Feynman-Kac Formulae P. Del Moral Centre INRIA de Bordeaux - Sud Ouest Workshop on Filtering, Cambridge Univ., June 14-15th 2010 Preprints (with hyperlinks), joint
More informationL p Functions. Given a measure space (X, µ) and a real number p [1, ), recall that the L p -norm of a measurable function f : X R is defined by
L p Functions Given a measure space (, µ) and a real number p [, ), recall that the L p -norm of a measurable function f : R is defined by f p = ( ) /p f p dµ Note that the L p -norm of a function f may
More informationSpring 2014: Computational and Variational Methods for Inverse Problems CSE 397/GEO 391/ME 397/ORI 397 Assignment 4 (due 14 April 2014)
Spring 2014: Computational and Variational Methods for Inverse Problems CSE 397/GEO 391/ME 397/ORI 397 Assignment 4 (due 14 April 2014) The first problem in this assignment is a paper-and-pencil exercise
More informationPHYS 5020 Computation and Image Processing
PHYS 5020 and Image Processing : Monte Carlo Thursday 2 August 2012 Monte Carlo (MC) is a numerical method that uses random sampling of probability distributions to simulate stochastic processes in nature,
More informationImportance Sampling in Monte Carlo Simulation of Rare Transition Events
Importance Sampling in Monte Carlo Simulation of Rare Transition Events Wei Cai Lecture 1. August 1, 25 1 Motivation: time scale limit and rare events Atomistic simulations such as Molecular Dynamics (MD)
More informationEstimation of probability density functions by the Maximum Entropy Method
Estimation of probability density functions by the Maximum Entropy Method Claudio Bierig, Alexey Chernov Institute for Mathematics Carl von Ossietzky University Oldenburg alexey.chernov@uni-oldenburg.de
More informationMonte Carlo Radiation Transport Kenny Wood
MCRT: L0 Some background, what previous courses students should look over, gentle introduction/recap of probabilities Get an idea of computer programming experience of the class Overview of course structure
More informationExpectation of geometric distribution. Variance and Standard Deviation. Variance: Examples
Expectation of geometric distribution Variance and Standard Deviation What is the probability that X is finite? Can now compute E(X): Σ k=f X (k) = Σ k=( p) k p = pσ j=0( p) j = p ( p) = E(X) = Σ k=k (
More informationSequential Monte Carlo Samplers for Applications in High Dimensions
Sequential Monte Carlo Samplers for Applications in High Dimensions Alexandros Beskos National University of Singapore KAUST, 26th February 2014 Joint work with: Dan Crisan, Ajay Jasra, Nik Kantas, Alex
More informationNonclassical Particle Transport in Heterogeneous Materials
Nonclassical Particle Transport in Heterogeneous Materials Thomas Camminady, Martin Frank and Edward W. Larsen Center for Computational Engineering Science, RWTH Aachen University, Schinkelstrasse 2, 5262
More informationThe multigroup Monte Carlo method part 1
The multigroup Monte Carlo method part 1 Alain Hébert alain.hebert@polymtl.ca Institut de génie nucléaire École Polytechnique de Montréal ENE6103: Week 11 The multigroup Monte Carlo method part 1 1/23
More informationExpectation of geometric distribution
Expectation of geometric distribution What is the probability that X is finite? Can now compute E(X): Σ k=1f X (k) = Σ k=1(1 p) k 1 p = pσ j=0(1 p) j = p 1 1 (1 p) = 1 E(X) = Σ k=1k (1 p) k 1 p = p [ Σ
More informationProbing the covariance matrix
Probing the covariance matrix Kenneth M. Hanson Los Alamos National Laboratory (ret.) BIE Users Group Meeting, September 24, 2013 This presentation available at http://kmh-lanl.hansonhub.com/ LA-UR-06-5241
More informationEstimating functional uncertainty using polynomial chaos and adjoint equations
0. Estimating functional uncertainty using polynomial chaos and adjoint equations February 24, 2011 1 Florida State University, Tallahassee, Florida, Usa 2 Moscow Institute of Physics and Technology, Moscow,
More informationSimulation. Version 1.1 c 2010, 2009, 2001 José Fernando Oliveira Maria Antónia Carravilla FEUP
Simulation Version 1.1 c 2010, 2009, 2001 José Fernando Oliveira Maria Antónia Carravilla FEUP Systems - why simulate? System any object or entity on which a particular study is run. Model representation
More informationSolving the steady state diffusion equation with uncertainty Final Presentation
Solving the steady state diffusion equation with uncertainty Final Presentation Virginia Forstall vhfors@gmail.com Advisor: Howard Elman elman@cs.umd.edu Department of Computer Science May 6, 2012 Problem
More informationSolution of Stochastic Nonlinear PDEs Using Wiener-Hermite Expansion of High Orders
Solution of Stochastic Nonlinear PDEs Using Wiener-Hermite Expansion of High Orders Dr. Mohamed El-Beltagy 1,2 Joint Wor with Late Prof. Magdy El-Tawil 2 1 Effat University, Engineering College, Electrical
More informationUse of probability gradients in hybrid MCMC and a new convergence test
Use of probability gradients in hybrid MCMC and a new convergence test Ken Hanson Methods for Advanced Scientific Simulations Group This presentation available under http://www.lanl.gov/home/kmh/ June
More informationOptimal Prediction for Radiative Transfer: A New Perspective on Moment Closure
Optimal Prediction for Radiative Transfer: A New Perspective on Moment Closure Benjamin Seibold MIT Applied Mathematics Mar 02 nd, 2009 Collaborator Martin Frank (TU Kaiserslautern) Partial Support NSF
More informationWave Turbulence and Condensation in an Optical Experiment
Wave Turbulence and Condensation in an Optical Experiment S. Residori, U. Bortolozzo Institut Non Linéaire de Nice, CNRS, France S. Nazarenko, J. Laurie Mathematics Institute, University of Warwick, UK
More informationMATH36061 Convex Optimization
M\cr NA Manchester Numerical Analysis MATH36061 Convex Optimization Martin Lotz School of Mathematics The University of Manchester Manchester, September 26, 2017 Outline General information What is optimization?
More informationSeismic Waves Propagation in Complex Media
H4.SMR/1586-1 "7th Workshop on Three-Dimensional Modelling of Seismic Waves Generation and their Propagation" 5 October - 5 November 004 Seismic Waves Propagation in Comple Media Fabio ROMANELLI Dept.
More informationLecture 4 Lebesgue spaces and inequalities
Lecture 4: Lebesgue spaces and inequalities 1 of 10 Course: Theory of Probability I Term: Fall 2013 Instructor: Gordan Zitkovic Lecture 4 Lebesgue spaces and inequalities Lebesgue spaces We have seen how
More informationAn introduction to PDE-constrained optimization
An introduction to PDE-constrained optimization Wolfgang Bangerth Department of Mathematics Texas A&M University 1 Overview Why partial differential equations? Why optimization? Examples of PDE optimization
More informationMarkov Chain Monte Carlo Methods
Markov Chain Monte Carlo Methods Sargur Srihari srihari@cedar.buffalo.edu 1 Topics Limitations of Likelihood Weighting Gibbs Sampling Algorithm Markov Chains Gibbs Sampling Revisited A broader class of
More informationAuxiliary Particle Methods
Auxiliary Particle Methods Perspectives & Applications Adam M. Johansen 1 adam.johansen@bristol.ac.uk Oxford University Man Institute 29th May 2008 1 Collaborators include: Arnaud Doucet, Nick Whiteley
More informationLecture 3: From Random Walks to Continuum Diffusion
Lecture 3: From Random Walks to Continuum Diffusion Martin Z. Bazant Department of Mathematics, MIT February 3, 6 Overview In the previous lecture (by Prof. Yip), we discussed how individual particles
More informationThe Larmor Formula (Chapters 18-19)
2017-02-28 Dispersive Media, Lecture 12 - Thomas Johnson 1 The Larmor Formula (Chapters 18-19) T. Johnson Outline Brief repetition of emission formula The emission from a single free particle - the Larmor
More informationPLATFORM DEPENDENT EFFICIENCY OF A MONTE CARLO CODE FOR TISSUE NEUTRON DOSE ASSESSMENT
RAD Conference Proceedings, vol., pp. 2-25, 206 www.rad-proceedings.org PLATFORM DEPENDENT EFFICIENCY OF A MONTE CARLO CODE FOR TISSUE NEUTRON DOSE ASSESSMENT Slobodan Milutinović, Filip Jeremić, Marko
More information2 Chance constrained programming
2 Chance constrained programming In this Chapter we give a brief introduction to chance constrained programming. The goals are to motivate the subject and to give the reader an idea of the related difficulties.
More informationNonlinear Computational Methods for Simulating Interactions of Radiation with Matter in Physical Systems
Nonlinear Computational Methods for Simulating Interactions of Radiation with Matter in Physical Systems Dmitriy Y. Anistratov Department of Nuclear Engineering North Carolina State University NE Seminar,
More information1 Random walks and data
Inference, Models and Simulation for Complex Systems CSCI 7-1 Lecture 7 15 September 11 Prof. Aaron Clauset 1 Random walks and data Supposeyou have some time-series data x 1,x,x 3,...,x T and you want
More informationAsymptotics and Simulation of Heavy-Tailed Processes
Asymptotics and Simulation of Heavy-Tailed Processes Department of Mathematics Stockholm, Sweden Workshop on Heavy-tailed Distributions and Extreme Value Theory ISI Kolkata January 14-17, 2013 Outline
More informationMONTE CARLO POWER ITERATION: ENTROPY AND SPATIAL CORRELATIONS
MONTE CARLO POWER ITERATION: ENTROPY AND SPATIAL CORRELATIONS ANDREA ZOIA, M. NOWAK (CEA/SACLAY) E. DUMONTEIL, A. ONILLON (IRSN) J. MIAO, B. FORGET, K. S. SMITH (MIT) NEA EGAMCT meeting Andrea ZOIA DEN/DANS/DM2S/SERMA/LTSD
More informationReduction of Systematic Error in Radiopharmaceutical Activity by Entropy Based Mutual Information
World Journal of Nuclear Science and Technology, 01,, 1-5 http://dx.doi.org/10.436/wjnst.01.1001 Published Online January 01 (http://www.scirp.org/journal/wjnst) Reduction of Systematic Error in Radiopharmaceutical
More informationA high accuracy approximation for half - space problems with anisotropic scattering
IOSR Journal of Applied Physics (IOSR-JAP e-iss: 78-486. Volume 4, Issue (May. - Jun., PP 7-46 A high accuracy approximation for half - space problems with anisotropic scattering M. H. Haggag,, H.F. Mashali
More informationMCRT L8: Neutron Transport
MCRT L8: Neutron Transport Recap fission, absorption, scattering, cross sections Fission products and secondary neutrons Slow and fast neutrons Energy spectrum of fission neutrons Nuclear reactor safety
More informationPartial Differential Equations and Homogenization for Simulating Aberrant Crypt Foci Dynamics
Partial Differential Equations and Homogenization for Simulating Aberrant Crypt Foci Dynamics Isabel Narra Figueiredo 1, Carlos Leal 1, Giuseppe Romanazzi 1, Björn Engquist 2 1 University of Coimbra, Portugal
More informationEnrico Fermi and the FERMIAC. Mechanical device that plots 2D random walks of slow and fast neutrons in fissile material
Monte Carlo History Statistical sampling Buffon s needles and estimates of π 1940s: neutron transport in fissile material Origin of name People: Ulam, von Neuman, Metropolis, Teller Other areas of use:
More informationExploring Monte Carlo Methods
Exploring Monte Carlo Methods William L Dunn J. Kenneth Shultis AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO ELSEVIER Academic Press Is an imprint
More informationKEYWORDS: chord length sampling, random media, radiation transport, Monte Carlo method, chord length probability distribution function.
On the Chord Length Sampling in 1-D Binary Stochastic Media Chao Liang and Wei Ji * Department of Mechanical, Aerospace, and Nuclear Engineering Rensselaer Polytechnic Institute, Troy, NY 12180-3590 *
More informationOptimal Control for Radiative Heat Transfer Model with Monotonic Cost Functionals
Optimal Control for Radiative Heat Transfer Model with Monotonic Cost Functionals Gleb Grenkin 1,2 and Alexander Chebotarev 1,2 1 Far Eastern Federal University, Sukhanova st. 8, 6995 Vladivostok, Russia,
More informationRandom Variable. Pr(X = a) = Pr(s)
Random Variable Definition A random variable X on a sample space Ω is a real-valued function on Ω; that is, X : Ω R. A discrete random variable is a random variable that takes on only a finite or countably
More informationSymmetries and reduction techniques for dissipative models
Symmetries and reduction techniques for dissipative models M. Ruggieri and A. Valenti Dipartimento di Matematica e Informatica Università di Catania viale A. Doria 6, 95125 Catania, Italy Fourth Workshop
More informationIII. Iterative Monte Carlo Methods for Linear Equations
III. Iterative Monte Carlo Methods for Linear Equations In general, Monte Carlo numerical algorithms may be divided into two classes direct algorithms and iterative algorithms. The direct algorithms provide
More informationUniversity of Cyprus. Reflectance and Diffuse Spectroscopy
University of Cyprus Biomedical Imaging and Applied Optics Reflectance and Diffuse Spectroscopy Spectroscopy What is it? from the Greek: spectro = color + scope = look at or observe = measuring/recording
More informationEECS 126 Probability and Random Processes University of California, Berkeley: Fall 2014 Kannan Ramchandran November 13, 2014.
EECS 126 Probability and Random Processes University of California, Berkeley: Fall 2014 Kannan Ramchandran November 13, 2014 Midterm Exam 2 Last name First name SID Rules. DO NOT open the exam until instructed
More informationSimple benchmark for evaluating self-shielding models
Simple benchmark for evaluating self-shielding models The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published Publisher
More informationA new class of interacting Markov Chain Monte Carlo methods
A new class of interacting Marov Chain Monte Carlo methods P Del Moral, A Doucet INRIA Bordeaux & UBC Vancouver Worshop on Numerics and Stochastics, Helsini, August 2008 Outline 1 Introduction Stochastic
More informationSimultaneous Robust Design and Tolerancing of Compressor Blades
Simultaneous Robust Design and Tolerancing of Compressor Blades Eric Dow Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics Massachusetts Institute of Technology GTL Seminar
More informationZero variance-based sampling schemes (a.k.a. path guiding) Jaroslav Křivánek Charles University, Prague Render Legion/Chaos Group
Zero variance-based sampling schemes (a.k.a. path guiding) Jaroslav Křivánek Charles University, Prague Render Legion/Chaos Group Path guiding The idea Zero-variance path sampling in volumes A theoretical
More informationNumerical Synthesis of Pontryagin Optimal Control Minimizers Using Sampling-Based Methods
Numerical Synthesis of Pontryagin Optimal Control Minimizers Using Sampling-Based Methods Humberto Gonzalez Department of Electrical & Systems Engineering Washington University in St. Louis April 7th,
More informationNonparametric density estimation for elliptic problems with random perturbations
Nonparametric density estimation for elliptic problems with random perturbations, DqF Workshop, Stockholm, Sweden, 28--2 p. /2 Nonparametric density estimation for elliptic problems with random perturbations
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 1. Title: Neutron Life Cycle
Lectures on Nuclear Power Safety Lecture No 1 Title: Neutron Life Cycle Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture Infinite Multiplication Factor, k Four Factor Formula
More informationComputational statistics
Computational statistics Markov Chain Monte Carlo methods Thierry Denœux March 2017 Thierry Denœux Computational statistics March 2017 1 / 71 Contents of this chapter When a target density f can be evaluated
More informationExpectation Propagation Algorithm
Expectation Propagation Algorithm 1 Shuang Wang School of Electrical and Computer Engineering University of Oklahoma, Tulsa, OK, 74135 Email: {shuangwang}@ou.edu This note contains three parts. First,
More informationMean Intensity. Same units as I ν : J/m 2 /s/hz/sr (ergs/cm 2 /s/hz/sr) Function of position (and time), but not direction
MCRT: L5 MCRT estimators for mean intensity, absorbed radiation, radiation pressure, etc Path length sampling: mean intensity, fluence rate, number of absorptions Random number generators J Mean Intensity
More informationRadiation Transport Calculations by a Monte Carlo Method. Hideo Hirayama and Yoshihito Namito KEK, High Energy Accelerator Research Organization
Radiation Transport Calculations by a Monte Carlo Method Hideo Hirayama and Yoshihito Namito KEK, High Energy Accelerator Research Organization Monte Carlo Method A method used to solve a problem with
More informationA Probability Primer. A random walk down a probabilistic path leading to some stochastic thoughts on chance events and uncertain outcomes.
A Probability Primer A random walk down a probabilistic path leading to some stochastic thoughts on chance events and uncertain outcomes. Are you holding all the cards?? Random Events A random event, E,
More informationON THE EXISTENCE OF TRANSMISSION EIGENVALUES. Andreas Kirsch1
Manuscript submitted to AIMS Journals Volume 3, Number 2, May 29 Website: http://aimsciences.org pp. 1 XX ON THE EXISTENCE OF TRANSMISSION EIGENVALUES Andreas Kirsch1 University of Karlsruhe epartment
More informationRobert Collins CSE586, PSU Intro to Sampling Methods
Intro to Sampling Methods CSE586 Computer Vision II Penn State Univ Topics to be Covered Monte Carlo Integration Sampling and Expected Values Inverse Transform Sampling (CDF) Ancestral Sampling Rejection
More informationThe Monte Carlo Method in Medical Radiation Physics
Monte Carlo in Medical Physics The Monte Carlo Method in Medical Radiation Physics P Andreo, Professor of Medical Radiation Physics Stockholm University @ Karolinska Univ Hospital, Stockholm, Sweden ICRM
More informationGeneral Lower Bounds for Resonances in One Dimension*
Commun. Math. Phys. 86, 221-225 (1982) Communications in Mathematical Physics Springer-Verlag 1982 General Lower Bounds for Resonances in One Dimension* Evans M. Harrell II Department of Mathematics, The
More informationPHAS3135 The Physics of Stars
PHAS3135 The Physics of Stars Exam 2013 (Zane/Howarth) Answer ALL SIX questions from Section A, and ANY TWO questions from Section B The numbers in square brackets in the right-hand margin indicate the
More informationA Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions
mathematics of computation volume 39, number 159 july 1982, pages 201-206 A Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions By John P. Boyd Abstract. The theorem proved here extends
More informationCS5314 Randomized Algorithms. Lecture 18: Probabilistic Method (De-randomization, Sample-and-Modify)
CS5314 Randomized Algorithms Lecture 18: Probabilistic Method (De-randomization, Sample-and-Modify) 1 Introduce two topics: De-randomize by conditional expectation provides a deterministic way to construct
More informationScattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion.
Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion. David M. Ambrose Jay Gopalakrishnan Shari Moskow Scott Rome June
More informationMCSHAPE: A MONTE CARLO CODE FOR SIMULATION OF POLARIZED PHOTON TRANSPORT
Copyright JCPDS - International Centre for Diffraction Data 2003, Advances in X-ray Analysis, Volume 46. 363 MCSHAPE: A MONTE CARLO CODE FOR SIMULATION OF POLARIZED PHOTON TRANSPORT J.E. Fernández, V.
More informationD. Guirado Instituto de Astrofísica de Andalucía CSIC, Granda, Spain Polarization school, Aussois, 06/06/2013
MONTE CARLO METHODS D. Guirado Instituto de Astrofísica de Andalucía CSIC, Granda, Spain Polarization school, Aussois, 06/06/2013 STOCHASTIC ALGORITHMS PROBLEM = GET OUTPUT FROM INPUT THROUGH AN ALGORITHM
More informationThe Radiative Transfer Equation
The Radiative Transfer Equation R. Wordsworth April 11, 215 1 Objectives Derive the general RTE equation Derive the atmospheric 1D horizontally homogenous RTE equation Look at heating/cooling rates in
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationBrief Introduction to: Transport Theory/Monte Carlo Techniques/MCNP
22:54 Neutron Interactions and Applications Brief Introduction to: Transport Theory/Monte Carlo Techniques/MCNP The behavior of individual neutrons and nuclei cannot be predicted. However, the average
More informationA LeVeque-type lower bound for discrepancy
reprinted from Monte Carlo and Quasi-Monte Carlo Methods 998, H. Niederreiter and J. Spanier, eds., Springer-Verlag, 000, pp. 448-458. A LeVeque-type lower bound for discrepancy Francis Edward Su Department
More informationLecture 16: October 29
CS294 Markov Chain Monte Carlo: Foundations & Applications Fall 2009 Lecture 16: October 29 Lecturer: Alistair Sinclair Scribes: Disclaimer: These notes have not been subjected to the usual scrutiny reserved
More informationNewton-Raphson. Relies on the Taylor expansion f(x + δ) = f(x) + δ f (x) + ½ δ 2 f (x) +..
2008 Lecture 7 starts here Newton-Raphson When the derivative of f(x) is known, and when f(x) is well behaved, the celebrated (and ancient) Newton- Raphson method gives the fastest convergence of all (
More informationBasics on quantum information
Basics on quantum information Mika Hirvensalo Department of Mathematics and Statistics University of Turku mikhirve@utu.fi Thessaloniki, May 2014 Mika Hirvensalo Basics on quantum information 1 of 49 Brief
More information