Monte Carlo method projective estimators for angular and temporal characteristics evaluation of polarized radiation
|
|
- Melina Smith
- 5 years ago
- Views:
Transcription
1 Monte Carlo method projective estimators for angular and temporal characteristics evaluation of polarized radiation Natalya Tracheva Sergey Ukhinov Institute of Computational Mathematics and Mathematical Geophysics SB RAS and Novosibirsk State University, Novosibirsk, Russia 1
2 Mathematical model 2 Let us consider the following stationary integro-differential radiative transfer equation with polarization: ω Φ(r, ω) + σφ(r, ω) = σ s P (ω, ω)φ(r, ω ) dω + f 0 (r, ω), or in the operator form Ω LΦ + σφ = SΦ + f 0, where Φ = (Φ 1, Φ 2, Φ 3, Φ 4 ) T is a vector function describing the flux density of particles (of vector photons ), or, in other words, the radiation intensity; Ω is the space of unit direction vectors, ω Ω, r D R 3 ; P (ω, ω) is the matrix scattering function, σ = σ(r) is the total cross section, σ = σ s + σ c, σ c is the absorption cross section, σ s is the scattering cross section; f 0 = (f (1) 0, f (2) 0, f (3) 0, f (4) 0 )T is the vector distribution density of particles.
3 Phase matrix definition 3 The matrix P (ω, ω) is defined by the relation P (ω, ω) = L(π i 2 )R(ω, ω)l( i 1 ), where L( ) is a special rotation matrix cos 2i sin 2i 0 L(i) = 0 sin 2i cos 2i 0, R( ) is the scattering matrix; i 1 is the angle between the plane ω, s and the scattering plane ω, ω ; i 2 is the angle between the scattering plane ω, ω and the plane ω, s; s = (0, 0, 1).
4 Scattering matrix 4 In the case of an isotropic medium, the scattering matrix is of the form: r 11 r R(µ, r) = 1 r 21 r π 0 0 r 33 r 34, 0 0 r 43 r 44 where µ = (ω, ω ) is cosine of the scattering angle, r 11 is the scattering function (or indicatrix), 1 1 r 11(µ) dµ = 1. If the scattering particles are homogeneous spheres, then r 11 = r 22, r 12 = r 21, r 33 = r 44, r 34 = r 43. In the case of molecular (Rayleigh) scattering, the matrix R is defined as follows: 3(1 + µ 2 )/8 3(1 µ 2 )/8 0 0 R(µ) = 1 2π 3(1 µ 2 )/8 3(1 + µ 2 )/ µ/ µ/4
5 Angular characteristics
6 Formal problem statement 6 Consider the radiation transport through a flat layer 0 < z < H of scattering and capturing media from the radiation source located on the border of z = 0 and directed in some direction ω 0. Vector-function Φ s (r, ω) of the angular distribution of the scattered radiation flux at r point on the surface z = h, 0 h H is Φ s (r, µ, ϕ) = µ Φ(r, µ, ϕ), where µ = (ω, s), ϕ is the azimuthal angle. Total radiation flux density (surface illumination) equal to P h = Φ (1) (r, ω) µ dω = 2π 1 Φ (1) (r, µ, ϕ) µ dµdϕ. Ω 0 0 In the case of isotropic radiation, the total flux density is P h Lambert = πφ (1) (r) and the corresponding angular distribution of the radiation flux on a surface with a normalized density Φ (1) s Lambert (r, µ, ϕ)/ph Lambert µ/π, µ (0, 1), ϕ (0, 2π) is well-known Lambert distribution.
7 Angular characteristics 7 Study angular distribution of backscattered and transmitted radiation intensity Φ (1) s (x, y)/p H, x (0, 1), y (0, π), degree of polarization p(r, ω) = Φ (2) (r, ω) 2 + Φ (3) (r, ω) 2 + Φ (4) (r, ω) 2. Φ (1) (r, ω)
8 Two types of the orthonormal function systems on (0, 1) (0, 2π) with weight x 8 1: {ψ i (x) φ j (y)}, i, j = 0,.., such that 2π x ψ i (x)φ k (y)ψ j (x)φ l (y) dx dy = { 1, i = j, k = l 0, else. 1 0 x ψ i (x)ψ j (x) dx = { 1, i = j 0, i j. 2π 0 φ k (y)φ l (y) dy = { 1, k = l 0, k l. 2: {H i,j (x, y), i = 0,..., j = i,.., i} such that 2π x H i,k (x, y)h j,l (x, y) dx dy = { 1, i = j, k = l 0, else.
9 The explicit form of the basis functions 1 9 N. V. Tracheva, S. A. Ukhinov (2017) Numerical statistical study of the angular distribution of the polarized radiation scattered by medium// Russian Journal of Numerical Analysis and Mathematical Modelling. Vol. 32(2), Tracheva, N. V. (2017) The use of two-dimensional projective Monte Carlo estimators for solution of number of problems of theory of radiation transfer// Proceedings of the Conference of young scientists of Institute of computational mathematics and mathematical geophysics of SB RAN. Novosibirsk, Functions ψ i (x) were obtained from Jacobi polynomials: ψ i (x) = 2i + 2 i k=0 ( 1) k (2 i + 1 k)! (i k)! k! (i + 1 k)! xi k. Functions φ j (x) were obtained from Legendre polynomials: φ j (y) = 2j + 1 ( 1) j 2π j k=0 (j + k)! (j k)! k! k! (2π) k ( x)k.
10 The explicit form of the basis functions 2 10 Mahotkin O. A. (1996) Analysis of radiative transfer between surfaces by hemispherical harmonics//journal of Quantitative Spectroscopy and Radiative Transfer. Vol. 56(6), The basis is constructed on the basis of hemispherical functions {H i,j (µ, ϕ)} of a form H 0,0 (µ, ϕ) = 1 i + 1 i ( 1) k (2 i + 1 k)! ; H i,0 (µ, ϕ) = π π (i k)! k! (i + 1 k)! µi k, i = 1, 2,...; k=0 2(i + 1) H i,j (µ, ϕ) = cos(jϕ) 2 j π j k=1 [(i + 1)2 k 2 ] i j i!(i + j + k + 1)!(i + j + 1)!(µ 1) k (1 (2µ 1) 2 ) j/2, i = 1, 2,..., j = 1,..., i; k!(l + j + 1)!(j + k)!(i j k)!(i + 1)! k=0 2(i + 1) H i, j (µ, ϕ) = sin(jϕ) 2 j π j k=1 [(i + 1)2 k 2 ] i j k=0 i!(i + j + k + 1)!(i + j + 1)!(µ 1) k (1 (2µ 1) 2 ) j/2, i = 1, 2,..., j = 1,..., i. k!(l + j + 1)!(j + k)!(i j k)!(i + 1)!
11 Series expansion 11 Vector-function of the angular distribution of the scattered radiation flux on the surface z = h, 0 h H : where Φ s (x, y) = x a i,j = 2π i i=0 j= i a i,j H i,j (x, y), Φ s (x, y)h i,j (x, y) dx dy.
12 Monte Carlo estimator 12 ξ i,j = N t 1 k=1 σ(h z k) q k Q k+1 H i,j (µ k+1, ϕ k+1 )e µ k+1 s,µ (z k, µ k+1 ), Eξ i,j = a i,j. ξ h = N t 1 k=1 σ(h z k) q k Q k+1 e µ k+1 s,µ (z k, µ k+1 ), Eξ h = P h. constructed on a Markov chain {x k }, where x k = (z k, ω k ), z k z-coordinate of the collision point, z 0 = 0; ω k the direction of the particle before the collision event, ω 1 = ω 0 ; µ k = (ω k, e z ); { 1, if zk < h and µ s,µ (z k, µ k+1 ) = k+1 > 0 0, else, i.e. s,µ the indicator of the surface z = h intersection by the particle direction after the scattering event in the point z k < h, N t random number of the last collision event. Direction of the scattering ω k in point z k 1 simulated according to the element R 11 of the scattering matrix R(ω k 1, ω k ). Vector weight Q k = ( Q (1) k, Q(2) k, Q(3) k, Q(4) k ) T defines by: Q 1 = ( I (1) 0, I(2) 0, I(3) 0, I(4) 0 ) T, Q k = P (ω k 1, ω k )Q k 1 /R 11 (ω k 1, ω k ), k 2.
13 Randomized projective estimator 13 Φ(x, y) where ξ (m) i,j n i i=0 j= i a i,j H i,j (x, y) = Φ n (x, y) Φ n (x, y) = α i,j = E N ξ i,j 1 N N m=0 ξ (m) i,j, random value realization ξ i,j on m-th trajectory, Eα i,j = Eξ i,j = a i,j, E Φ n (x, y) = Φ n (x, y). n i i=0 j= i α i,j H i,j (x, y),
14 Randomized projective estimator variation 14 Componentwise representation of variation of random vectorfunction Φ n (x) D Φ n (x, y) = n i n j i=0 k= i j=0 l= j cov(α i,k, α j,l )H i,k (x, y)h j,l (x, y) = 1 N n i n j i=0 k= i j=0 l= j cov(ξ i,k, ξ j,l )H i,k (x, y)h j,l (x, y).
15 Polarization degree estimator 15 p n (x, y) = ( Φ (2) n (x, y)) 2 + ( Φ (3) n (x, y)) 2 + ( Φ (4) n (x, y)) 2 Φ (1) n (x, y).
16 The variance of the asymptotic distribution of p n (x, y) estimator 16 D N p n (x) = 1 g 2 1 ( g D N g 1 g g2 2 D Ng 2 + g3 2 D Ng 3 + g4 2 D Ng 4 g 2 cov N(g 1, g 2 ) g 2 + cov N (g 1, g 3 ) g 3 + cov N (g 1, g 4 ) g 4 g cov N(g 2, g 3 ) g 2 g 3 + cov N (g 2, g 4 ) g 2 g 4 + cov N (g 3, g 4 ) g 3 g 4 g + o(n 1 ), where g m = Φ (m) n (x, y), g = g g2 3 + g2 4, ) + cov N (g m, g t ) = 1 N n i n j i=0 k= i j=0 l= j cov N (ξ (m) i,k, ξ(t) j,l ) H i,k(x, y)h j,l (x, y), m, t = 1,..., 4.
17 Temporal characteristics
18 Exponential asymptotics parameter 18 An exponential asymptotics parameter is the principal eigenvalue λ homogeneous stationary transfer equation: of the LΦ + (σ + λ/v)φ = SΦ, with standard boundary conditions (see Davison B. Neutron Transport Theory. Clarendon, Oxford ). In the work of Mikhailov G. A, Tracheva N. V., Uhinov S. A. (Monte Carlo study of time asymptotics of the polarized radiation intensity // Computational Mathematics and Mathematical Physics V. 47, N. 7.) this statement was extended to polarized radiation. In the spatially homogeneous case (i.e., when the entire space is filled with a homogeneous medium), this goal can be achieved rather easily through weighted simulation of radiative transfer. Moreover, it proves that λ = λ = σ cv irrespective of the polarization type. Here, σ c is the absorption cross section and v is the velocity of the particles. It is well known that λ = λ = σ cv for a half-space as well.
19 A problem of estimation of time asymptotic parameters of polarized radiation in scattering and capturing media 19 A number of publications, e.g. Romanova L. M. (1965) Limit cases of the free path distribution function of photons exiting a thick light-diffusing layer.// Izv. Akad. Nauk SSSR, Ser. Fiz. Atmos. Okeana., 1, No. 6, Zege E.P., Katsev I.L. (1973)Time asymptotical solutions of radiation transfer equation and their application. Preprint. Inst. Fiz. Akad. Nauk. BSSR, Minsk. show that for intensity functionals without polarization the following asymptotic relation exists: J(r, ω, t) C(r, ω) t α e σ cvt, t +. We are concerned with obtaining of the similar relation in pre-asymptotical times.
20 System of orthonormal functions 20 Let us consider a system ψ i (t) of orthonormal over (0, ) with weight e σ cvt functions, such that 0 e σ cvt ψ i (t)ψ j (t) dt = { 1, i = j, 0, else. where Φ(r, ω, t) = C(r, ω)e σ cvt a i ψ i (t), a i = C(r, ω) 0 i=0 Φ(t)ψ i (t) dt.
21 Monte Carlo estimator 21 ξ i = N t 1 k=1 q k Q k+1 ψ i (t k+1 ), Eξ i = a i built on Markov chain {x k }, x k = (r k, ω k ), r k coordinates of the collision point, r 0 = (0, 0, 0); ω k is the direction of the particle before the collision, ω 1 = ω 0 ; µ k = (ω k, e z ), N t random number of the last collision event, the scattering direction ω k at the point r k 1 is simulated according to an element R 11 of scattering matrix R(ω k 1, ω k ). Vector weigt Q k = ( Q (1) k is determined by Q 1 = ( I (1) 0, I(2) 0, I(3) 0, ) T I(4) 0,, Q(2) k, Q(3) k Q k = P (ω k 1, ω k )Q k 1 /R 11 (ω k 1, ω k ), k 2., Q(4) k ) T
22 Randomized projective Monte Carlo estimator 22 Φ(r, ω, t) C(r, ω)e σ cvt n i=0 a i ψ i (t) = Φ n (r, ω, t) where ξ (m) i Φ n (r, ω, t) = C(r, ω)e σ cvt α i = E N ξ i 1 N N m=0 n i=0 ξ (m) i, α i ψ i (t), random value ξ i realization on m-th trajectory, Eα i = Eξ i = a i, E Φ n (r, ω, t) = Φ n (r, ω, t).
23 Variation of projective Monte Carlo estimator 23 Variation of random function Φ n (t) is equal to D Φ n (t) = C 2 (r, ω)e 2σ cvt n n i=0 j=0 cov(α i, α j )ψ i (t)ψ j (t) = = C(r, ω) 2 e 2σ cvt 1 N n n i=0 j=0 cov(ξ i, ξ j )ψ i (t)ψ j (t).
24 Basis functions 24 As orthonormal basis let us consider a modyfied Laguerre polynomials L n (x), orthogonal with weight e x over (0, ). We can obtain the following explicit form for functions ψ i (t): ψ i (x) = σ c v i k=0 ( 1) k i! k! (i k)! k! σ c k v k t k.
25 Numerical results
26 Simulated medium characteristics 26 Deirmendjian, D. (1969)Electromagnetic Scattering on Spherical Polydispersions. Amer. Elsevier, New York. Aerosol scattering. Coefficient of refraction on particles n = i (water), particle size distribution is logarithmically normal with density f(r) = 1 r exp( 1 ln 2 ( r 2σg 2 r g )), r (0, 10mkm), r g = 0.12mkm, σ g = 0.5, wave-length of radiation equals to 0.65 mkm. Mean cosine of scattering angle - µ 0 = Marchuk, G. I., Mikhailov, G. A., Nazaraliev, M. A. et al. (1980)The Monte Carlo Methods in Atmospheric Optics. Amer. Springer-Verlag, Heidelberg. Molecular scattering. 3(1 + µ 2 )/8 3(1 µ 2 )/8 0 0 R(µ) = 1 2π 3(1 µ 2 )/8 3(1 + µ 2 )/ µ/ µ/4
27 Angular distribution function of the normalized flux of radiation Φ (1) s (cos(θ), ϕ)/p h, backscattered by a layer of optical thickness H = 5. Aerozol scattering. θ 0 = 45, ϕ 0 = 0. Hemispherical basis. 27 Lambertian approximation n = 2
28 Angular distribution function of the normalized flux of radiation Φ (1) s (cos(θ), ϕ)/p h, backscattered by a layer of optical thickness H = 5. Aerozol scattering. θ 0 = 45, ϕ 0 = 0. Hemispherical basis. 28 n = 3 n = 5
29 Angular distribution function of the normalized flux of radiation Φ (1) s (cos(θ), ϕ)/p h, transmitted by a layer of optical thickness H = 5. Aerozol scattering. θ 0 = 45, ϕ 0 = 0. Factorized basis, n i = n j = n. 29 n = 1 n = 3
30 Angular distribution function of the normalized flux of radiation Φ (1) s (cos(θ), ϕ)/p h, transmitted by a layer of optical thickness H = 5. Aerozol scattering. θ 0 = 45, ϕ 0 = 0. Factorized basis, n i = n j = n. 30 n = 5 n = 8
31 Angular distribution functions of the normalized flux of radiation Φ (1) s (cos(θ), ϕ)/p h and the degree of polarization p(cos(θ), ϕ), backscattered (1-st row) and transmitted (2-d row) by a layer of optical thickness H = 5. Molecular scattering. θ 0 = 45, ϕ 0 = 0. Hemispherical basis, n = Φ (1) s (cos(θ), ϕ)/p h p(cos(θ), ϕ)
32 Angular distribution functions of the normalized flux of radiation Φ (1) s (cos(θ), ϕ)/p h and the degree of polarization p(cos(θ), ϕ), backscattered (1-st row) and transmitted (2-d row) by a layer of optical thickness H = 5. Molecular scattering. θ 0 = 45, ϕ 0 = 0. Hemispherical basis, n = Φ (1) s (cos(θ), ϕ)/p h p(cos(θ), ϕ)
33 Angular characteristics of the radiation, backscattered by the layer of optical thickness H = 10. (a) (c) The normalized radiation flux Φ (1) s (cos(θ), ϕ)/p h. (d) (f) The degree of polarization p n (cos(θ), ϕ). ϕ 0 = 0. Hemispherical basis. 33 θ 0 = 30 θ 0 = 45 θ 0 = 60
34 Angular characteristics of the radiation, backscattered by the layer of optical thickness H = 10. (a) (c) The normalized radiation flux Φ (1) s (cos(θ), ϕ)/p h. (d) (f) The degree of polarization p n (cos(θ), ϕ). ϕ 0 = 0. Hemispherical basis. 34 θ 0 = 30 θ 0 = 45 θ 0 = 60
35 σ c = 0.1. Solid line is for an exponential curve, dashed line is for molecular scattering matrix curve, dash-and-dot line is for aerosol scattering matrix curve. (a) Time distribution of radiation intensity, normalized over total flux, in comparison to an exponential asymptotics e σ cvt. (b) Degree of polarization. 35 (a) (b)
36 σ c = Solid line is for an exponential curve, dashed line is for molecular scattering matrix curve, dash-and-dot line is for aerosol scattering matrix curve. (a) Time distribution of radiation intensity, normalized over total flux, in comparison to an exponential asymptotics e σ cvt. (b) Degree of polarization. 36 (a) (b)
37 Thank you for attention!
Light Scattering in Inhomogeneous Atmospheres
Edgard G. Yanovitskij Light Scattering in Inhomogeneous Atmospheres Translated by Sergij Ginsheimer and Oleg Yanovitskij With 37 Figures and 54 Tables Springer Contents Introduction 1 1. Basic Concepts,
More information6. LIGHT SCATTERING 6.1 The first Born approximation
6. LIGHT SCATTERING 6.1 The first Born approximation In many situations, light interacts with inhomogeneous systems, in which case the generic light-matter interaction process is referred to as scattering
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 informationOptical Imaging Chapter 5 Light Scattering
Optical Imaging Chapter 5 Light Scattering Gabriel Popescu University of Illinois at Urbana-Champaign Beckman Institute Quantitative Light Imaging Laboratory http://light.ece.uiuc.edu Principles of Optical
More informationPC4262 Remote Sensing Scattering and Absorption
PC46 Remote Sensing Scattering and Absorption Dr. S. C. Liew, Jan 003 crslsc@nus.edu.sg Scattering by a single particle I(θ, φ) dφ dω F γ A parallel beam of light with a flux density F along the incident
More informationSpectral analysis for a class of linear pencils arising in transport theory
Spectral analysis for a class of linear pencils arising in transport theory Petru A. Cojuhari AGH University of Science and Technology Kraków, Poland Vienna, December 17 20, 2016 Petru A. Cojuhari Spectral
More informationMCRT L10: Scattering and clarification of astronomy/medical terminology
MCRT L10: Scattering and clarification of astronomy/medical terminology What does the scattering? Shape of scattering Sampling from scattering phase functions Co-ordinate frames Refractive index changes
More informationSOLUTION OF PARAMETRIC INVERSE PROBLEM OF ATMOSPHERIC OPTICS BY MONTE CARLO METHODS
TWMS Jour. Pure Appl. Math. V.3, N.2, 2012, pp.220-230 SOLUTION OF PARAMETRIC INVERSE PROBLEM OF ATMOSPHERIC OPTICS BY MONTE CARLO METHODS KANAT SHAKENOV 1 Abstract. In the article parametric inverse problems
More informationLecture 06. Fundamentals of Lidar Remote Sensing (4) Physical Processes in Lidar
Lecture 06. Fundamentals of Lidar Remote Sensing (4) Physical Processes in Lidar Physical processes in lidar (continued) Doppler effect (Doppler shift and broadening) Boltzmann distribution Reflection
More informationAnalytic Assessment of Eigenvalues of the Neutron Transport Equation *
Analytic Assessment of igenvalues of the Neutron Transport quation * Sung Joong KIM Massachusetts Institute of Technology Nuclear Science & ngineering Department May 0, 2005 Articles for Current Study
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 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 informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More information2 The Radiative Transfer Equation
9 The Radiative Transfer Equation. Radiative transfer without absorption and scattering Free space or homogeneous space I (r,,) I (r,,) r -r d da da Figure.: Following a pencil of radiation in free space
More informationAnalysis of Scattering of Radiation in a Plane-Parallel Atmosphere. Stephanie M. Carney ES 299r May 23, 2007
Analysis of Scattering of Radiation in a Plane-Parallel Atmosphere Stephanie M. Carney ES 299r May 23, 27 TABLE OF CONTENTS. INTRODUCTION... 2. DEFINITION OF PHYSICAL QUANTITIES... 3. DERIVATION OF EQUATION
More informationMultiple scattering of light by water cloud droplets with external and internal mixing of black carbon aerosols
Chin. Phys. B Vol. 21, No. 5 (212) 5424 Multiple scattering of light by water cloud droplets with external and internal mixing of black carbon aerosols Wang Hai-Hua( 王海华 ) and Sun Xian-Ming( 孙贤明 ) School
More informationRegularities of Angular Distribution of Near-Horizon Sky Brightness in the Cloudless Atmosphere
Regularities of Angular Distribution of Near-Horizon Sky Brightness in the Cloudless Atmosphere S.M. Sakerin, T.B. Zhuravleva, and I.M. Nasrtdinov Institute of Atomospheric Optics SB RAS Tomsk, Russia
More informationScattering cross-section (µm 2 )
Supplementary Figures Scattering cross-section (µm 2 ).16.14.12.1.8.6.4.2 Total scattering Electric dipole, a E (1,1) Magnetic dipole, a M (1,1) Magnetic quardupole, a M (2,1). 44 48 52 56 Wavelength (nm)
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 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 informationPreface to the Second Edition. Preface to the First Edition
Contents Preface to the Second Edition Preface to the First Edition iii v 1 Introduction 1 1.1 Relevance for Climate and Weather........... 1 1.1.1 Solar Radiation.................. 2 1.1.2 Thermal Infrared
More informationSerco Assurance. Resonance Theory and Transport Theory in WIMSD J L Hutton
Serco Assurance Resonance Theory and Transport Theory in WIMSD J L Hutton 2 March 2004 Outline of Talk Resonance Treatment Outline of problem - pin cell geometry U 238 cross section Simple non-mathematical
More informationLesson 6: Diffusion Theory (cf. Transport), Applications
Lesson 6: Diffusion Theory (cf. Transport), Applications Transport Equation Diffusion Theory as Special Case Multi-zone Problems (Passive Media) Self-shielding Effects Diffusion Kernels Typical Values
More informationAn asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations
An asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations Institute of Applied Physics and Computational Mathematics, Beijing NUS, Singapore, March 2-6, 2015 (joint
More informationMonte Carlo Sampling
Monte Carlo Sampling Sampling from PDFs: given F(x) in analytic or tabulated form, generate a random number ξ in the range (0,1) and solve the equation to get the randomly sampled value X X ξ = F(x)dx
More informationUniversity of Illinois at Chicago Department of Physics. Electricity & Magnetism Qualifying Examination
University of Illinois at Chicago Department of Physics Electricity & Magnetism Qualifying Examination January 7, 28 9. am 12: pm Full credit can be achieved from completely correct answers to 4 questions.
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 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 informationPhys 622 Problems Chapter 6
1 Problem 1 Elastic scattering Phys 622 Problems Chapter 6 A heavy scatterer interacts with a fast electron with a potential V (r) = V e r/r. (a) Find the differential cross section dσ dω = f(θ) 2 in the
More informationANALYTICAL ONE-PARAMETER SEA WATER LIGHT SCATTERING PHASE FUNCTION WITH INTEGRAL PROPERTIES OF EXPERIMENTAL PHASE FUNCTIONS
ANALYTICAL ONE-PARAMETER SEA WATER LIGHT SCATTERING PHASE FUNCTION WITH INTEGRAL PROPERTIES OF EXPERIMENTAL PHASE FUNCTIONS Vladimir I. Haltrin Naval Research Laboratory, Ocean Optics Section, Code 7333,
More informationMore On Carbon Monoxide
More On Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results Jerry Gilfoyle The Configurations of CO 1 / 26 More On Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results Jerry Gilfoyle The Configurations
More informationThe collision probability method in 1D part 1
The collision probability method in 1D part 1 Alain Hébert alain.hebert@polymtl.ca Institut de génie nucléaire École Polytechnique de Montréal ENE6101: Week 8 The collision probability method in 1D part
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationLaser Beam Interactions with Solids In absorbing materials photons deposit energy hc λ. h λ. p =
Laser Beam Interactions with Solids In absorbing materials photons deposit energy E = hv = hc λ where h = Plank's constant = 6.63 x 10-34 J s c = speed of light Also photons also transfer momentum p p
More informationFundametals of Rendering - Radiometry / Photometry
Fundametals of Rendering - Radiometry / Photometry Physically Based Rendering by Pharr & Humphreys Chapter 5: Color and Radiometry Chapter 6: Camera Models - we won t cover this in class Realistic Rendering
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationCollection of formulae Quantum mechanics. Basic Formulas Division of Material Science Hans Weber. Operators
Basic Formulas 17-1-1 Division of Material Science Hans Weer The de Broglie wave length λ = h p The Schrödinger equation Hψr,t = i h t ψr,t Stationary states Hψr,t = Eψr,t Collection of formulae Quantum
More informationIntroduction to Polarization
Phone: Ext 659, E-mail: hcchui@mail.ncku.edu.tw Fall/007 Introduction to Polarization Text Book: A Yariv and P Yeh, Photonics, Oxford (007) 1.6 Polarization States and Representations (Stokes Parameters
More informationMATH 31BH Homework 5 Solutions
MATH 3BH Homework 5 Solutions February 4, 204 Problem.8.2 (a) Let x t f y = x 2 + y 2 + 2z 2 and g(t) = t 2. z t 3 Then by the chain rule a a a D(g f) b = Dg f b Df b c c c = [Dg(a 2 + b 2 + 2c 2 )] [
More informationPolarized light propagation and scattering in random media
Polarized light propagation and scattering in random media Arnold D. Kim a, Sermsak Jaruwatanadilok b, Akira Ishimaru b, and Yasuo Kuga b a Department of Mathematics, Stanford University, Stanford, CA
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationSupplementary Information
1 Supplementary Information 3 Supplementary Figures 4 5 6 7 8 9 10 11 Supplementary Figure 1. Absorbing material placed between two dielectric media The incident electromagnetic wave propagates in stratified
More informationFUNDAMENTALS OF POLARIZED LIGHT
FUNDAMENTALS OF POLARIZED LIGHT A STATISTICAL OPTICS APPROACH Christian Brosseau University of Brest, France A WILEY-INTERSCIENCE PUBLICATION JOHN WILEY & SONS, INC. New York - Chichester. Weinheim. Brisbane
More informationAy121 Problem Set 1. TA: Kathryn Plant. October Some parts are based on solutions by previous TAs Yuguang Chen and Rachel Theios.
Ay2 Problem Set TA: Kathryn Plant October 28 Some parts are based on solutions by previous TAs Yuguang Chen and Rachel Theios. Temperatures One point for each part. (a) Since we are observing at MHz and
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 informationPHYS 771, Quantum Mechanics, Final Exam, Fall 2011 Instructor: Dr. A. G. Petukhov. Solutions
PHYS 771, Quantum Mechanics, Final Exam, Fall 11 Instructor: Dr. A. G. Petukhov Solutions 1. Apply WKB approximation to a particle moving in a potential 1 V x) = mω x x > otherwise Find eigenfunctions,
More informationGeneralized Burgers equations and Miura Map in nonabelian ring. nonabelian rings as integrable systems.
Generalized Burgers equations and Miura Map in nonabelian rings as integrable systems. Sergey Leble Gdansk University of Technology 05.07.2015 Table of contents 1 Introduction: general remarks 2 Remainders
More informationUNIVERSITY OF DUBLIN
UNIVERSITY OF DUBLIN TRINITY COLLEGE JS & SS Mathematics SS Theoretical Physics SS TSM Mathematics Faculty of Engineering, Mathematics and Science school of mathematics Trinity Term 2015 Module MA3429
More informationFun With Carbon Monoxide. p. 1/2
Fun With Carbon Monoxide p. 1/2 p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results C V (J/K-mole) 35 30 25
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 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 informationWhat is it good for? RT is a key part of remote sensing and climate modeling.
Read Bohren and Clothiaux Ch.; Ch 4.-4. Thomas and Stamnes, Ch..-.6; 4.3.-4.3. Radiative Transfer Applications What is it good for? RT is a key part of remote sensing and climate modeling. Remote sensing:
More informationSeparation of Variables in Polar and Spherical Coordinates
Separation of Variables in Polar and Spherical Coordinates Polar Coordinates Suppose we are given the potential on the inside surface of an infinitely long cylindrical cavity, and we want to find the potential
More informationInteraction X-rays - Matter
Interaction X-rays - Matter Pair production hν > M ev Photoelectric absorption hν MATTER hν Transmission X-rays hν' < hν Scattering hν Decay processes hν f Compton Thomson Fluorescence Auger electrons
More informationElectrodynamics Exam Solutions
Electrodynamics Exam Solutions Name: FS 215 Prof. C. Anastasiou Student number: Exercise 1 2 3 4 Total Max. points 15 15 15 15 6 Points Visum 1 Visum 2 The exam lasts 18 minutes. Start every new exercise
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
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 information10. Atmospheric scattering. Extinction ( ε ) = absorption ( k ) + scattering ( m ): ε. = = single scattering albedo (SSA).
1. Atmospheric scattering Extinction ( ε ) = absorption ( k ) + scattering ( m ): ε σ = kσ + mσ, mσ mσ = = single scattering albedo (SSA). k + m ε σ σ σ As before, except for polarization (which is quite
More informationWave Propagation in Heterogeneous Media: Born and Rytov Approximations. Chris Sherman
Wave Propagation in Heterogeneous Media: Born and Rytov Approximations Chris Sherman Stochastic Scalar Wave Equation Wave Equation: &! " %& 1 # t V x ' ) () u(x,t) = 0 (*) Velocity Perturbation: = V o
More informationON THE CORRECT FORMULATION OF A MULTIDIMENSIONAL PROBLEM FOR STRICTLY HYPERBOLIC EQUATIONS OF HIGHER ORDER
Georgian Mathematical Journal 1(1994), No., 141-150 ON THE CORRECT FORMULATION OF A MULTIDIMENSIONAL PROBLEM FOR STRICTLY HYPERBOLIC EQUATIONS OF HIGHER ORDER S. KHARIBEGASHVILI Abstract. A theorem of
More informationParticle Interactions in Detectors
Particle Interactions in Detectors Dr Peter R Hobson C.Phys M.Inst.P. Department of Electronic and Computer Engineering Brunel University, Uxbridge Peter.Hobson@brunel.ac.uk http://www.brunel.ac.uk/~eestprh/
More informationReview Article Mathematical Methods in Biomedical Optics
ISRN Biomedical Engineering Volume 2013, Article ID 464293, 8 pages http://dx.doi.org/10.1155/2013/464293 Review Article Mathematical Methods in Biomedical Optics Macaveiu Gabriela Academic Center of Optical
More informationMathematical Tripos Part IB Michaelmas Term Example Sheet 1. Values of some physical constants are given on the supplementary sheet
Mathematical Tripos Part IB Michaelmas Term 2015 Quantum Mechanics Dr. J.M. Evans Example Sheet 1 Values of some physical constants are given on the supplementary sheet 1. Whenasampleofpotassiumisilluminatedwithlightofwavelength3
More informationLet b be the distance of closest approach between the trajectory of the center of the moving ball and the center of the stationary one.
Scattering Classical model As a model for the classical approach to collision, consider the case of a billiard ball colliding with a stationary one. The scattering direction quite clearly depends rather
More informationChapter 5 Fast Multipole Methods
Computational Electromagnetics; Chapter 1 1 Chapter 5 Fast Multipole Methods 5.1 Near-field and far-field expansions Like the panel clustering, the Fast Multipole Method (FMM) is a technique for the fast
More informationSummary: angular momentum derivation
Summary: angular momentum derivation L = r p L x = yp z zp y, etc. [x, p y ] = 0, etc. (-) (-) (-3) Angular momentum commutation relations [L x, L y ] = i hl z (-4) [L i, L j ] = i hɛ ijk L k (-5) Levi-Civita
More informationEfficiency Random Walks Algorithms for Solving BVP of meta Elliptic Equations
Efficiency Random Walks Algorithms for Solving BVP of meta Elliptic Equations Vitaliy Lukinov Russian Academy of Sciences ICM&MG SB RAS Novosibirsk, Russia Email: Vitaliy.Lukinov@gmail.com Abstract In
More informationInfluence of an Electric Field on the Propagation of a Photon in a Magnetic field
Journal of Physics: Conference Series PAPER OPEN ACCESS Influence of an Electric Field on the Propagation of a Photon in a Magnetic field To cite this article: V M Katkov 06 J. Phys.: Conf. Ser. 73 0003
More informationWAVE PROPAGATION AND SCATTERING IN RANDOM MEDIA
WAVE PROPAGATION AND SCATTERING IN RANDOM MEDIA AKIRA ISHIMARU UNIVERSITY of WASHINGTON IEEE Antennas & Propagation Society, Sponsor IEEE PRESS The Institute of Electrical and Electronics Engineers, Inc.
More informationLecture notes 5: Diffraction
Lecture notes 5: Diffraction Let us now consider how light reacts to being confined to a given aperture. The resolution of an aperture is restricted due to the wave nature of light: as light passes through
More information221B Lecture Notes Notes on Spherical Bessel Functions
Definitions B Lecture Notes Notes on Spherical Bessel Functions We would like to solve the free Schrödinger equation [ h d l(l + ) r R(r) = h k R(r). () m r dr r m R(r) is the radial wave function ψ( x)
More informationeff (r) which contains the influence of angular momentum. On the left is
1 Fig. 13.1. The radial eigenfunctions R nl (r) of bound states in a square-well potential for three angular-momentum values, l = 0, 1, 2, are shown as continuous lines in the left column. The form V (r)
More informationApplied Nuclear Physics (Fall 2006) Lecture 19 (11/22/06) Gamma Interactions: Compton Scattering
.101 Applied Nuclear Physics (Fall 006) Lecture 19 (11//06) Gamma Interactions: Compton Scattering References: R. D. Evans, Atomic Nucleus (McGraw-Hill New York, 1955), Chaps 3 5.. W. E. Meyerhof, Elements
More informationElectromagnetic optics!
1 EM theory Electromagnetic optics! EM waves Monochromatic light 2 Electromagnetic optics! Electromagnetic theory of light Electromagnetic waves in dielectric media Monochromatic light References: Fundamentals
More informationAstro 305 Lecture Notes Wayne Hu
Astro 305 Lecture Notes Wayne Hu Set 1: Radiative Transfer Radiation Observables From an empiricist s point of view there are 4 observables for radiation Energy Flux Direction Color Polarization Energy
More informationHigh-Order Similarity Relations in Radiative Transfer: Supplementary Document
High-Order Similarity Relations in Radiative Transfer: Supplementary Document Shuang hao Ravi Ramamoorthi Kavita Bala Cornell University University of California, Berkeley Cornell University Overview This
More informationModelling the Directionality of Light Scattered in Translucent Materials
Modelling the Directionality of Light Scattered in Translucent Materials Jeppe Revall Frisvad Joint work with Toshiya Hachisuka, Aarhus University Thomas Kjeldsen, The Alexandra Institute March 2014 Materials
More informationFinite element investigations for the construction of composite solutions of even-parity transport equation
Finite element investigations for the construction of composite solutions of even-parity transport equation Anwar M. Mirza and Shaukat Iqbal Abstract. Finite element investigations have been carried out
More informationLecture 06. Fundamentals of Lidar Remote Sensing (4) Physical Processes in Lidar
Lecture 06. Fundamentals of Lidar Remote Sensing (4) Physical Processes in Lidar Light interactions with objects (continued) Resonance fluorescence Laser induced fluorescence Doppler effect (Doppler shift
More informationNEUTRAL PARTICLE TRANSPORT IN PLASMA EDGE USING TRANSMISSION/ESCAPE PROBABILITY (TEP) METHOD
NEUTRAL PARTICLE TRANSPORT IN PLASMA EDGE USING TRANSMISSION/ESCAPE PROBABILITY (TEP) METHOD A Dissertation Presented to The Academic Faculty By Dingkang Zhang In Partial Fulfillment of the Requirements
More informationScattering in Cold- Cathode Discharges
Simulating Electron Scattering in Cold- Cathode Discharges Alexander Khrabrov, Igor Kaganovich*, Vladimir I. Demidov**, George Petrov*** *Princeton Plasma Physics Laboratory ** Wright-Patterson Air Force
More informationFORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked
More informationLecture 10. Lidar Effective Cross-Section vs. Convolution
Lecture 10. Lidar Effective Cross-Section vs. Convolution q Introduction q Convolution in Lineshape Determination -- Voigt Lineshape (Lorentzian Gaussian) q Effective Cross Section for Single Isotope --
More informationarxiv: v1 [hep-th] 10 Dec 2015
arxiv:52.03203v [hep-th] 0 Dec 205 Influence of an Electric Field on the Propagation of a Photon in a Magnetic field V.M. Katkov Budker Institute of Nuclear Physics, Novosibirsk, 630090, Russia e-mail:
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 informationPolynomial Chaos and Karhunen-Loeve Expansion
Polynomial Chaos and Karhunen-Loeve Expansion 1) Random Variables Consider a system that is modeled by R = M(x, t, X) where X is a random variable. We are interested in determining the probability of the
More informationStochastic modeling of Extinction coefficients for solar power applications
1(77) Stochastic modeling of Extinction coefficients for solar power applications Ingemar Mathiasson January 2007 Department for Energy and Environment Division of Electric Power Engineering Chalmers University
More informationSolar radiation / radiative transfer
Solar radiation / radiative transfer The sun as a source of energy The sun is the main source of energy for the climate system, exceeding the next importat source (geothermal energy) by 4 orders of magnitude!
More informationLecture 6 Scattering theory Partial Wave Analysis. SS2011: Introduction to Nuclear and Particle Physics, Part 2 2
Lecture 6 Scattering theory Partial Wave Analysis SS2011: Introduction to Nuclear and Particle Physics, Part 2 2 1 The Born approximation for the differential cross section is valid if the interaction
More informationStatistical mechanics of random billiard systems
Statistical mechanics of random billiard systems Renato Feres Washington University, St. Louis Banff, August 2014 1 / 39 Acknowledgements Collaborators: Timothy Chumley, U. of Iowa Scott Cook, Swarthmore
More informationGravitational lensing by gravitational wave pulse
G.S. Bisnovatyi-Kogan abc and ac a Space Research Institute of Russian Academy of Science, Profsoyuznaya 84/3, Moscow 7997, Russia b Joint Institute for Nuclear Research, Dubna, Russia c Moscow Engineering
More informationarxiv:physics/ v2 [physics.optics] 6 Apr 2007
Unified theory for Goos-Hänchen and Imbert-Fedorov effects Chun-Fang Li Department of Physics, Shanghai University, arxiv:physics/0611047v2 [physics.optics] 6 Apr 2007 Shanghai 200444, P. R. China and
More informationDesign of one-dimensional Lambertian diffusers of light
INSTITUTE OF PHYSICS PUBLISHING Waves Random Media 11 (2001) 529 533 WAVES IN RANDOM MEDIA PII: S0959-7174(01)18801-4 Design of one-dimensional Lambertian diffusers of light A A Maradudin 1, I Simonsen
More informationIl picco di Bragg. G. Battistoni INFN Milano. 08/06/2015 G. Battistoni
Il picco di Bragg G. Battistoni INFN Milano 08/06/015 G. Battistoni 1 Φ(z) The physics of Bragg Peak 180 MeV proton in water Longitudinal profile: Transversel profile: Φ(z,x) dominated by interaction with
More informationSection 22. Radiative Transfer
OPTI-01/0 Geometrical and Instrumental Optics Copyright 018 John E. Greivenkamp -1 Section Radiative Transfer Radiometry Radiometry characterizes the propagation of radiant energy through an optical system.
More informationIntroduction to Elementary Particle Physics I
Physics 56400 Introduction to Elementary Particle Physics I Lecture 2 Fall 2018 Semester Prof. Matthew Jones Cross Sections Reaction rate: R = L σ The cross section is proportional to the probability of
More informationSOLUTIONS OF THE LINEAR BOLTZMANN EQUATION AND SOME DIRICHLET SERIES
SOLUTIONS OF THE LINEAR BOLTZMANN EQUATION AND SOME DIRICHLET SERIES A.V. BOBYLEV( ) AND I.M. GAMBA( ) Abstract. It is shown that a broad class of generalized Dirichlet series (including the polylogarithm,
More informationRadiative heat transfer
Radiative heat transfer 22 mars 2017 Energy can be transported by the electromagnetic field radiated by an object at finite temperature. A very important example is the infrared radiation emitted towards
More informationIntroduction to Condensed Matter Physics
Introduction to Condensed Matter Physics Diffraction I Basic Physics M.P. Vaughan Diffraction Electromagnetic waves Geometric wavefront The Principle of Linear Superposition Diffraction regimes Single
More informationFokker-Planck Equation with Detailed Balance
Appendix E Fokker-Planck Equation with Detailed Balance A stochastic process is simply a function of two variables, one is the time, the other is a stochastic variable X, defined by specifying: a: the
More information