Nonlinear Computational Methods for Simulating Interactions of Radiation with Matter in Physical Systems
|
|
- Harold Ray
- 5 years ago
- Views:
Transcription
1 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, Oregon State University October 14,
2 Introduction The simulation of interaction of particles with matter is a challenging problem. The transport equation is a basis for mathematical models of this physical phenomenon. Transport problems are difficult to solve high dimensionality, an integro-differential equation, coefficients depend on the state of matter the state of matter is affected by fluxes of particles. Applications: reactor physics, astrophysics (stars), plasma physics (laser fusion), atmospheric sciences. 2
3 Introduction There exists a family of efficient nonlinear methods for solving the transport equation: Quasidiffusion method (V. Gol din, 1964), Nonlinear Diffusion Acceleration (K. Smith, 2002), Flux methods (T. Germogenova, V. Gol din, 1969). These methods are defined by a system of nonlinearly coupled high-order and low-order problems that is equivalent to the original linear transport problem. The Nonlinear Projective Iteration (NPI) methods possess certain advantages for their use in multiphysics applications. The low-order equations of NPI methods can be used to formulate approximate particle transport models. These low-order equations can also be utilized as a basis for development of hybrid Monte Carlo computational methods. The NPI methods are distinct from each other by the definition of the low-order equations which gives rise to differences in features of these methods. 3
4 Introduction Outline A Tutorial on nonlinear methods for solving the transport equation The Quasidiffusion (QD) method for transport problems in 2D Cartesian geometry on grids composed of arbitrary quadrilaterals Nonlinear Weighted Flux (NWF) methods for particle transport problems in 2D Cartesian geometry on orthogonal grid This work was performed in collaboration with my Ph.D. students at NCSU William A. Wieselquist (to graduate in 2008 Paul Scherrer Institute, Switzerland) Loren Roberts (graduated in 2008 Baker Hughes Inc.) 4
5 Transport Problem Let us consider the following single-group slab geometry transport problem with isotropic scattering and source: μ x ψ(x, μ)+σ t(x)ψ(x, μ) = 1 1 Σs (x) 2 μ 1, 0 x L, ψ(x, μ )dμ + Q(x), φ(x) = 1 ψ(x, μ)dμ, J(x) = 1 μψ(x, μ)dμ, with reflective left boundary ψ(0,μ)=ψ(0, μ), 0 <μ 1, and vacuum boundary condition on the right ψ(l, μ) =0, μ<0. 5
6 Low-Order Equations The transport equation is integrated over μ 1 with weights 1 and μ dj dx +(Σ t Σ s )φ = Q, d 1 μ 2 ψdμ +Σt J =0, dx where φ(x) = 1 ψ(x, μ)dμ, J(x) = 1 μψ(x, μ)dμ. What should we do with the extra moment? 1 μ 2 ψdμ =? 6
7 Approximate Closure: Diffusion Theory Diffusion theory (approximation) ψ(x, μ) = 1 2 (φ(x)+3μj(x)), 1 μ 2 ψdμ = 1 3 φ(x). P 1 equations Diffusion equation dj dx +(Σ t Σ s )φ = Q, 1 dφ 3 dx +Σ tj =0, d 1 dφ dx 3Σ t dx +(Σ t Σ s )φ = Q, 7
8 Approximate Closure: Variable Eddington Factor Method 1 μ 2 ψdμ F(x)φ(x) Minerbo closure (approximation) ψ(x, μ) =α(x)e μβ(x) F(x) =1 2 J(x) Z(x) φ(x), Kershaw closure (approximation) ( F(x) = 1 ( ) ) 2 J(x) φ(x) Levermore-Pomraning closure (approximation) F(x) = J(x) φ(x) coth(z(x)), coth(z) 1 Z = J φ P 1 -like nonlinear low-order equations dj dx +(Σ t Σ s )φ = Q, d dx (Fφ)+Σ tj =0, F = F(φ, J) 8
9 Quasidiffusion (QD) Method: Exact Closure The Quasidiffusion method (exact) 1 1 μ 2 ψdμ μ 2 ψdμ = 1 ψdμ The quasidiffusion factor. E(x) = μ 2 ψdμ ψdμ ψdμ = E(x)φ(x), =< μ 2 > Quasidiffusion low-order equations dj dx +(Σ t Σ s )φ = Q, deφ dx +Σ tj =0, E = E[ψ]. 9
10 System of Equations of the Quasidiffusion (QD) Method The transport (high-order) equation μ x ψ(x, μ)+σ tψ(x, μ) = 1 ( Σs φ(x)+q ). 2 QD (variable Eddington) factors E(x) = 1 1 μ 2 ψ(x, μ)dμ ψ(x, μ)dμ, C L = the low-order QD equations d dx J(x)+(Σ t Σ s )φ(x) =Q, ( ) d E(x)φ(x) +Σ t J(x) =0. dx 0 0 μψ(l, μ)dμ ψ(l, μ)dμ. J(0) = 0, J(L) =C L φ(l). 10
11 QD Method: Iteration Process High-order problem (transport sweep) μ x ψ(s+1/2) +Σ t ψ (s+1/2) = 1 2 (Σ sφ (s) + Q), Calculation of QD factors E (s+1/2) (x) = 1 μ2 ψ (s+1/2) (x, μ)dμ 1 ψ(s+1/2) (x, μ)dμ, C L = Low-order QD problem d dx J (s+1) +(Σ t Σ s )φ (s+1) = Q, d ) (E (s+1/2) φ (s+1) +Σ t J (s+1) =0, dx J (s+1) (0) = 0, 0 0 μψ (s+1/2) (L, μ)dμ ψ (s+1/2) (L, μ)dμ. J (s+1) (L) =C (s+1/2) L φ (s+1) (L). 11
12 Nonlinear Diffusion Acceleration (NDA) High-order problem (transport sweep) Calculation of D (s+1/2) J (s+1/2) (x) = μ x ψ(s+1/2) +Σ t ψ (s+1/2) = 1 2 (Σ sφ (s) + Q), 1 μψ (s+1/2) (x, μ)dμ, φ (s+1/2) (x) = D (s+1/2) = 1 φ (s+1/2) ( 1 J (s+1/2) + 1 ) dφ (s+1/2) 3Σ t dx Low-order NDA problem d dx J (s+1) +(Σ t Σ s )φ (s+1) = Q, 1 dφ (s+1) 3Σ t dx + D (s+1/2) φ (s+1) + J (s+1) =0, (Kord Smith) ψ (s+1/2) (x, μ)dμ,. 12
13 QD Method for the Multidimensional Transport Equation The transport (high-order) equation Ω ψ (k+1/2) ( r, Ω) + σ t ( r)ψ (k+1/2) ( r, Ω) = 1 4π σ s( r)φ (k) ( r)+ 1 4π q( r). The QD (Eddington) tensor E (k+1/2) αβ ( r) = 4π Ē = ( ) Exx E xy E xy E yy Ω α Ω β ψ (k+1/2) ( r, Ω)d Ω / 4π E αβ ( r) = average value of Ω α Ω β. The low-order QD (LOQD) equations J (k+1) ( r)+σ a ( r)φ (k+1) ( r) =q( r), J (k+1) ( r) = 1 ) (Ē (k+1/2) ( r)φ k+1 ( r). σ t ( r), ψ (k+1/2) ( r, Ω)d Ω, α,β = x, y, 13
14 Discretization of the LOQD Equations The unknowns: φ i, φ iω,andj iω = J iω n iω, The balance equation in ith cell J iω A iω + σ a,i φ i V i = q i V i, where ω cell face. ω Jim Morel (2002) proposed a cell-centered discretization for the diffusion equation on meshes of arbitrary polyhedrons. { φ} iω =(φ iω φ i ) α iω + 1 (φ iω ) β iωω, V i ω where α iω and β iωω are based on geometry. We apply Morel s method J iω = 1 { } (Ēφ). σ t,i iω { } ({ } { } ) ({ } Exx φ Exy φ Exy φ (Ēφ) = e x + + e y iω x iω y iω x iω + { } ) Eyy φ y iω. 14
15 Discretization of the LOQD Equations Conditions at interfaces between cells Strong current and weak scalar flux continuity conditions J R = J 2L J R = J 3L φ R A R = φ 2L A 2L + φ 3L A 3L 15
16 Analytic Test Problem for the LOQD Equations 0 x 1, 0 y 1. σ t =1, σ a =0.5. Analytic solution [ ( φ(x, y) =5 tanh 100 x 2) 1 2 ( y 1 ) ] 2 2 The QD tensor (E αβ (x, y)) is given in analytic form φ(x, y) E xx (x, y) E xy (x, y) 16
17 Results on 3-Level Grids Relative error in the cell-average scalar flux Grid 1: 40 cells orthogonal grid randomized (20%) grid 17
18 Results on 3-Level Grids Relative error in the cell-average scalar flux Grid 2: 160 cells orthogonal grid randomized (20%) grid 18
19 Results on 3-Level Grids Relative error in the cell-average scalar flux Grid 3: 640 cells orthogonal grid randomized (20%) grid 19
20 Results on 3-Level Grids Relative error in the cell-average scalar flux Grid 4: 2560 cells orthogonal grid randomized (20%) grid 20
21 Transport Test Problem 21
22 Numerical Results Results of Single Level Grids Orthogonal Grid Randomized Grid Relative Difference F φ L φ R F φ L φ R F φ L φ R e e e-2-4.0e-3 2.2e e e e-3-2.2e-3 1.4e e e e-3-2.4e-4 5.0e e e e-3 1.2e-4 8.4e-4 Numerically Estimated Spatial Order of Convergence Orthogonal Grid Randomized Grid F φ L φ R F φ L φ R single-level two-level (with refinement on the left) two-level (with refinement on the right)
23 Nonlinear Weighted Flux (NWF) Methods The transport (high-order) equation The factors G m ( r) =Γ m F α m ( r) =Γ m Ω ψ( r, Ω) + σ t ( r)ψ( r, Ω) = 1 4π σ s( r)φ( r)+ 1 4π q( r). ω m ω m w(ω x, Ω y )ψ( r, Ω)d Ω / Ω α w(ω x, Ω y )ψ( r, Ω)d Ω ψ( r, Ω)d Ω, Γ m = / ω m dω ω m w(ω x, Ω y )dω. ω m ψ( r, Ω)d Ω, α = x, y, The low-order NWF equations φ m ( r) = ω m ψ( r, Ω)d Ω, m =1,...,4 ν x m ( ) F x m φ m + ν y m x y ( F y m φ m ) + σt G m φ m = 1 4 (σ sφ + q), φ = 4 m=1 φ m., ν x 1 = ν y 1 =1, νx 2 =, ν y 2 =1,νx 3 = ν y 3 =, νx 4 =1, ν y 4 =, 23
24 Asymptotic Diffusion Limit of the Transport Eqation We define a small parameter ε, scaled cross sections and source: Consider that ε 0. σ t = ˆσ t ε, σ a = εˆσ a, q = εˆq, We now assume that the solution can be expanded in power series of ε ψ = ε n ψ [n], n=0 The leading-order solutio of the transport the equation meets to the diffusion equation 1 ψ [0] x3σ t x 1 ψ [0] y 3σ t y + σ aψ [0] = q. 24
25 Asymptotic Diffusion Limit of the Low-Order NWF Equations We consider various weights w(ω x, Ω y )=1 w(ω x, Ω y )= Ω x + Ω y w(ω x, Ω y )=1+ Ω x + Ω y w(ω x, Ω y )=1+β( Ω x + Ω y ) Values of the diffusion coefficients (D) for specific NWF methods Weight 1 Ω x + Ω y 1+ Ω x + Ω y 1+β( Ω x + Ω y ) ( 1 D π+2 ) 2 1 ( 4σ t 3π σ t 1 4+5π ) σ t 12π σ t σ t 3σ t 25
26 Numerical Results Test problem A unit square domain σ t = 1 ɛ, σ a = ɛ, andq = ɛ ɛ =10 2, 10 3, 10 4, and 10 5 Boundary conditions are vacuum. A uniform spatial mesh consists of equal cells. The compatible product quadruple-range quadrature set uses 3 polar and 3 azimuthal angles per octant. Relative errors of the cell-average scalar flux in the cell located at the center Weight ɛ 1 Ω x + Ω y 1+ Ω x + Ω y 1+β( Ω x + Ω y ) E E E E E E E E E E E E E E E E-3 26
27 Asymptotic Boundary Condition There exist a boundary layer with the width of order of ε. The scalar flux at the boundary of the diffusive domain φ [0] (X, y) =2 W ( n Ω )ψ in (X, y, Ω)d Ω, n Ω<0 It is possible to improve the performance of the NWF method by using w(ω x, Ω y )=1+λ( Ω x + Ω y )+κ Ω x Ω y Relative errors of the asymptotic boundary conditions for the NWF methods BLD Transport Weight τ (mfp) Method 1 Ω x + Ω y 1+ Ω x + Ω y w β w λ,κ % 0.033% 0.033% 0.033% 0.033% % % % % % % % % % % % % 0.697% % % % % % 1.36% % % % % % 1.25% % % % % % 0.755% 27
28 Problem with an Unresolved Boundary Layer A domain 0 x, y 11 with two subregions: 1. 0 x 1, σ t =2,σ s =0,q =0,Δx =0.1, and Δy =1, 2. 1 x 11, σ t = σ s = 100, q =0,andΔx =Δy =1. There is an isotropic incoming angular flux with magnitude 1 boundary 2π on the left Other boundary condition are vacuum. This problem enables one to test the ability of a method to reproduce an accurate diffusion solution in the interior of the diffusive region with the spatially unresolved boundary layer. 28
29 Problem with an Unresolved Boundary Layer 29
30 Conclusions We have developed a new method for solving the LOQD equations on arbitrary quadrilateral and AMR-like grids The resulting QD method for solving transport problems demonstrated good performance on both randomized and AMR-like grids. Observe second order spatial convergence in numerical tests We defined a new parameterized family of nonlinear flux methods for solving the 2D transport equation. The asymptotic diffusion analysis enabled us to find a particular method of this family the solution of which satisfies a good approximation of the diffusion equation in diffusive regions. As a result, we developed a 2D flux method with an important feature for solving transport problems with optically thick regions. Note that it is possible to formulate a linear version of the proposed methods. 30
ABSTRACT. flux methods for the accurate and efficient solution of multidimensional transport problems in
ABSTRACT ROBERTS, LOREN P. Nonlinear Weighted Flux Methods for Solving Multidimensional Transport Problems. (Under the direction of Dmitriy Y. Anistratov.) The goal of the research presented in this dissertation
More informationAn Angular Multigrid Acceleration Method for S n Equations with Highly Forward-Peaked Scattering
An Angular Multigrid Acceleration Method for S n Equations with Highly Forward-Peaked Scattering Bruno Turcksin & Jean C. Ragusa & Jim E. Morel Texas A&M University, Dept. of Nuclear Engineering Bruno
More informationHomogenization Methods for Full Core Solution of the Pn Transport Equations with 3-D Cross Sections. Andrew Hall October 16, 2015
Homogenization Methods for Full Core Solution of the Pn Transport Equations with 3-D Cross Sections Andrew Hall October 16, 2015 Outline Resource-Renewable Boiling Water Reactor (RBWR) Current Neutron
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 informationNonlinear Iterative Solution of the Neutron Transport Equation
Nonlinear Iterative Solution of the Neutron Transport Equation Emiliano Masiello Commissariat à l Energie Atomique de Saclay /DANS//SERMA/LTSD emiliano.masiello@cea.fr 1/37 Outline - motivations and framework
More informationRECXWH2 W/o s3-1
.,, Nuclear Energy Research Initiative (NERI) DE-FG03-99SF21922 Program Tasks 1 & 3 Quarterly Technical Progress Report July 1- September 30, 2000 An Innovative Reactor Analysis Methodology Based on a
More informationKey words. finite elements, discontinuous Galerkin, neutron transport, diffusion limit
DISCONTINUOUS GALERKIN METHODS FOR THE RADIATIVE TRANSPORT EQUATION. JEAN-LUC GUERMOND, G. KANSCHAT, AND J.-C. RAGUSA Abstract. This notes presents some recent results regarding the approximation of the
More informationSummary of time independent electrodynamics
hapter 10 Summary of time independent electrodynamics 10.1 Electrostatics Physical law oulomb s law charges as origin of electric field Superposition principle ector of the electric field E(x) in vacuum
More informationLecture 2 Solutions to the Transport Equation
Lecture 2 Solutions to the Transport Equation Equation along a ray I In general we can solve the static transfer equation along a ray in some particular direction. Since photons move in straight lines
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 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 informationTHREE-DIMENSIONAL INTEGRAL NEUTRON TRANSPORT CELL CALCULATIONS FOR THE DETERMINATION OF MEAN CELL CROSS SECTIONS
THREE-DIMENSIONAL INTEGRAL NEUTRON TRANSPORT CELL CALCULATIONS FOR THE DETERMINATION OF MEAN CELL CROSS SECTIONS Carsten Beckert 1. Introduction To calculate the neutron transport in a reactor, it is often
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 information3D Cartesian Transport Sweep for Massively Parallel Architectures on top of PaRSEC
3D Cartesian Transport Sweep for Massively Parallel Architectures on top of PaRSEC 9th Scheduling for Large Scale Systems Workshop, Lyon S. Moustafa, M. Faverge, L. Plagne, and P. Ramet S. Moustafa, M.
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 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 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 informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationPhD Qualifying Exam Nuclear Engineering Program. Part 1 Core Courses
PhD Qualifying Exam Nuclear Engineering Program Part 1 Core Courses 9:00 am 12:00 noon, November 19, 2016 (1) Nuclear Reactor Analysis During the startup of a one-region, homogeneous slab reactor of size
More informationEFFECTS OF THE TEMPERATURE DISCRETIZATION ON NUMERICAL METHODS FOR THERMAL RADIATION TRANSPORT
International Conference on Mathematics, Computational Methods & Reactor Physics M&C 2009 Saratoga Springs, New York, May 3 7, 2009, on CD-ROM, American Nuclear Society, LaGrange Park, IL 2009 EFFECTS
More informationStars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11)
Stars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11) Kenneth Wood, Room 316 kw25@st-andrews.ac.uk http://www-star.st-and.ac.uk/~kw25 What is a Stellar Atmosphere? Transition from dense
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 informationStellar Atmospheres: Basic Processes and Equations
Stellar Atmospheres: Basic Processes and Equations Giovanni Catanzaro Abstract The content of this chapter is a very quick summary of key concepts that concern the interaction between photons created in
More informationIterative Properties of Parallel Block Jacobi-Integral Transport Matrix Method with Source Iteration Preconditioning
M&C 017 - International Conference on Mathematics & Computational Methods Applied to Nuclear Science & Engineering, Jeju, Korea, April 16-0, 017, on USB 017) Iterative Properties of Parallel Block Jacobi-Integral
More informationLecture 3 Numerical Solutions to the Transport Equation
Lecture 3 Numerical Solutions to the Transport Equation Introduction I There are numerous methods for solving the transport problem numerically. First we must recognize that we need to solve two problems:
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 informationFundamental Stellar Parameters
Fundamental Stellar Parameters Radiative Transfer Specific Intensity, Radiative Flux and Stellar Luminosity Observed Flux, Emission and Absorption of Radiation Radiative Transfer Equation, Solution and
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 informationExtensions of the TEP Neutral Transport Methodology. Dingkang Zhang, John Mandrekas, Weston M. Stacey
Extensions of the TEP Neutral Transport Methodology Dingkang Zhang, John Mandrekas, Weston M. Stacey Fusion Research Center, Georgia Institute of Technology, Atlanta, GA 30332-0425, USA Abstract Recent
More information1. (a) (5 points) Find the unit tangent and unit normal vectors T and N to the curve. r (t) = 3 cos t, 0, 3 sin t, r ( 3π
1. a) 5 points) Find the unit tangent and unit normal vectors T and N to the curve at the point P 3, 3π, r t) 3 cos t, 4t, 3 sin t 3 ). b) 5 points) Find curvature of the curve at the point P. olution:
More informationconditional cdf, conditional pdf, total probability theorem?
6 Multiple Random Variables 6.0 INTRODUCTION scalar vs. random variable cdf, pdf transformation of a random variable conditional cdf, conditional pdf, total probability theorem expectation of a random
More informationA Model for Radiation Interactions with Matter
Eastern Michigan University DigitalCommons@EMU Senior Honors Theses Honors College 2004 A Model for Radiation Interactions with Matter Carrie Beck Follow this and additional works at: http://commons.emich.edu/honors
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 informationTo remain at 0 K heat absorbed by the medium must be removed in the amount of. dq dx = q(l) q(0) dx. Q = [1 2E 3 (τ L )] σ(t T 4 2 ).
294 RADIATIVE HEAT TRANSFER 13.7 Two infinite, isothermal plates at temperatures T 1 and T 2 are separated by a cold, gray medium of optical thickness τ L = L (no scattering). (a) Calculate the radiative
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 informationReview of Fundamental Equations Supplementary notes on Section 1.2 and 1.3
Review of Fundamental Equations Supplementary notes on Section. and.3 Introduction of the velocity potential: irrotational motion: ω = u = identity in the vector analysis: ϕ u = ϕ Basic conservation principles:
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationStability and dispersion analysis of high order FDTD methods for Maxwell s equations in dispersive media
Contemporary Mathematics Volume 586 013 http://dx.doi.org/.90/conm/586/11666 Stability and dispersion analysis of high order FDTD methods for Maxwell s equations in dispersive media V. A. Bokil and N.
More informationDIFFUSION PRECONDITIONER FOR DISCONTINUOUS GALERKIN TRANSPORT PROBLEMS. A Thesis ANTHONY PETRU BARBU
DIFFUSION PRECONDITIONER FOR DISCONTINUOUS GALERKIN TRANSPORT PROBLEMS A Thesis by ANTHONY PETRU BARBU Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationRadiative transfer equation in spherically symmetric NLTE model stellar atmospheres
Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres Jiří Kubát Astronomický ústav AV ČR Ondřejov Zářivě (magneto)hydrodynamický seminář Ondřejov 20.03.2008 p. Outline 1.
More information2. The Steady State and the Diffusion Equation
2. The Steady State and the Diffusion Equation The Neutron Field Basic field quantity in reactor physics is the neutron angular flux density distribution: Φ( r r, E, r Ω,t) = v(e)n( r r, E, r Ω,t) -- distribution
More informationAND ITERATIVE EFFICIENCY. A Dissertation YUNHUANG ZHANG
EVEN-PARITY S N ADJOINT METHOD INCLUDING SP N MODEL ERROR AND ITERATIVE EFFICIENCY A Dissertation by YUNHUANG ZHANG Submitted to the Office of Graduate and Professional Studies of Texas A&M University
More informationA DISCRETE ORDINATES APPROXIMATION TO THE NEUTRON TRANSPORT EQUATION APPLIED TO GENERALIZED GEOMETRIES. A Dissertation.
A DISCRETE ORDINATES APPROXIMATION TO THE NEUTRON TRANSPORT EQUATION APPLIED TO GENERALIZED GEOMETRIES A Dissertation by MARK DAVID DeHART Submitted to the Office of Graduate Studies of Texas A&M University
More informationSimulating experiments for ultra-intense laser-vacuum interaction
Simulating experiments for ultra-intense laser-vacuum interaction Nina Elkina LMU München, Germany March 11, 2011 Simulating experiments for ultra-intense laser-vacuum interaction March 11, 2011 1 / 23
More informationTwo Posts to Fill On School Board
Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83
More informationMixed Black Hole - Neutron Star Simulations with Whisky
Mixed Black Hole - Neutron Star Simulations with Whisky Roland Haas (Georgia Tech) Center for Relativistic Astrophysics Collaborators: Tanja Bode, Jim Healy, Pablo Laguna, Deirdre Shoemaker Oct 06, 2009
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 informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More informationStatus of the MAGIX Spectrometer Design. Julian Müller MAGIX collaboration meeting 2017
Status of the MAGIX Spectrometer Design Julian Müller MAGIX collaboration meeting 2017 Magneto Optic Design Requirements relative momentum resolution Δp p < 10 4 resolution of the scattering angle Δθ
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
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 informationMultivariate Random Variable
Multivariate Random Variable Author: Author: Andrés Hincapié and Linyi Cao This Version: August 7, 2016 Multivariate Random Variable 3 Now we consider models with more than one r.v. These are called multivariate
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationPlasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University
Plasma Astrophysics Chapter 1: Basic Concepts of Plasma Yosuke Mizuno Institute of Astronomy National Tsing-Hua University What is a Plasma? A plasma is a quasi-neutral gas consisting of positive and negative
More informationNonlinear wave-wave interactions involving gravitational waves
Nonlinear wave-wave interactions involving gravitational waves ANDREAS KÄLLBERG Department of Physics, Umeå University, Umeå, Sweden Thessaloniki, 30/8-5/9 2004 p. 1/38 Outline Orthonormal frames. Thessaloniki,
More informationDirect Sum. McKelvey, and Madie Wilkin Adviser: Dr. Andrew Christlieb Co-advisers: Eric Wolf and Justin Droba. July 23, Michigan St.
and Adviser: Dr. Andrew Christlieb Co-advisers: Eric Wolf and Justin Droba Michigan St. University July 23, 2014 The N-Body Problem Solar systems The N-Body Problem Solar systems Interacting charges, gases
More informationLecture 3: Specific Intensity, Flux and Optical Depth
Lecture 3: Specific Intensity, Flux and Optical Depth We begin a more detailed look at stellar atmospheres by defining the fundamental variable, which is called the Specific Intensity. It may be specified
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationA Tutorial on Wavelets and their Applications. Martin J. Mohlenkamp
A Tutorial on Wavelets and their Applications Martin J. Mohlenkamp University of Colorado at Boulder Department of Applied Mathematics mjm@colorado.edu This tutorial is designed for people with little
More informationTypical anisotropies introduced by geometry (not everything is spherically symmetric) temperature gradients magnetic fields electrical fields
Lecture 6: Polarimetry 1 Outline 1 Polarized Light in the Universe 2 Fundamentals of Polarized Light 3 Descriptions of Polarized Light Polarized Light in the Universe Polarization indicates anisotropy
More informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More informationPhysics 236a assignment, Week 2:
Physics 236a assignment, Week 2: (October 8, 2015. Due on October 15, 2015) 1. Equation of motion for a spin in a magnetic field. [10 points] We will obtain the relativistic generalization of the nonrelativistic
More informationThe Scattering of Electromagnetic Waves from Two-Dimensional Randomly Rough Surfaces
The Scattering of Electromagnetic Waves from Two-Dimensional Randomly Rough Surfaces Ingve Simonsen 1, Alexei A. Maradudin and Tamara A. Leskova 1 Department of Physics Norwegian University of Science
More informationLecture10: Plasma Physics 1. APPH E6101x Columbia University
Lecture10: Plasma Physics 1 APPH E6101x Columbia University Last Lecture - Conservation principles in magnetized plasma frozen-in and conservation of particles/flux tubes) - Alfvén waves without plasma
More informationHomogeneous Bianisotropic Medium, Dissipation and the Non-constancy of Speed of Light in Vacuum for Different Galilean Reference Systems
Progress In Electromagnetics Research Symposium Proceedings, Cambridge, USA, July 5 8, 2010 489 Homogeneous Bianisotropic Medium, Dissipation and the Non-constancy of Speed of Light in Vacuum for Different
More informationElectromagnetically Induced Flows in Water
Electromagnetically Induced Flows in Water Michiel de Reus 8 maart 213 () Electromagnetically Induced Flows 1 / 56 Outline 1 Introduction 2 Maxwell equations Complex Maxwell equations 3 Gaussian sources
More informationFundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres
Fundamental Stellar Parameters Radiative Transfer Stellar Atmospheres Equations of Stellar Structure Basic Principles Equations of Hydrostatic Equilibrium and Mass Conservation Central Pressure, Virial
More informationHigh performance computing for neutron diffusion and transport equations
High performance computing for neutron diffusion and transport equations Horizon Maths 2012 Fondation Science Mathématiques de Paris A.-M. Baudron, C. Calvin, J. Dubois, E. Jamelot, J.-J. Lautard, O. Mula-Hernandez
More informationIn this method, one defines
1 Feautrier s Method for Radiative Transfer This project is to use Feautrier s method to solve the Gray Atmosphere problem for a plane-parallel atmosphere Although we have developed alternate techniques
More informationChemistry 431. NC State University. Lecture 17. Vibrational Spectroscopy
Chemistry 43 Lecture 7 Vibrational Spectroscopy NC State University The Dipole Moment Expansion The permanent dipole moment of a molecule oscillates about an equilibrium value as the molecule vibrates.
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 9
WiSe 202 20.2.202 Prof. Dr. A-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg
More informationLong-wave Instability in Anisotropic Double-Diffusion
Long-wave Instability in Anisotropic Double-Diffusion Jean-Luc Thiffeault Institute for Fusion Studies and Department of Physics University of Texas at Austin and Neil J. Balmforth Department of Theoretical
More informationSolitons and point particles
Solitons and point particles Martin Speight University of Leeds July 28, 2009 What are topological solitons? Smooth, spatially localized solutions of nonlinear relativistic classical field theories Stable;
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationEMAFF: MAgnetic Equations with FreeFem++
EMAFF: MAgnetic Equations with FreeFem++ The Grad-Shafranov equation & The Current Hole Erwan DERIAZ Bruno DESPRÉS Gloria FACCANONI š Kirill Pichon GOSTAF Lise-Marie IMBERT-GÉRARD Georges SADAKA Remy
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
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 informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationVorticity Equation Marine Hydrodynamics Lecture 9. Return to viscous incompressible flow. N-S equation: v. Now: v = v + = 0 incompressible
13.01 Marine Hydrodynamics, Fall 004 Lecture 9 Copyright c 004 MIT - Department of Ocean Engineering, All rights reserved. Vorticity Equation 13.01 - Marine Hydrodynamics Lecture 9 Return to viscous incompressible
More informationRadiation hydrodynamics of tin targets for laser-plasma EUV sources
Radiation hydrodynamics of tin targets for laser-plasma EUV sources M. M. Basko, V. G. Novikov, A. S. Grushin Keldysh Institute of Applied Mathematics, Moscow, Russia RnD-ISAN, Troitsk, Moscow, Russia
More information' Liberty and Umou Ono and Inseparablo "
3 5? #< q 8 2 / / ) 9 ) 2 ) > < _ / ] > ) 2 ) ) 5 > x > [ < > < ) > _ ] ]? <
More informationComparison of Neutrino Transport Approximations in Core-Collapse Supernova Simulations
Comparison of Neutrino Transport Approximations in Core-Collapse Supernova Simulations Sherwood Richers 1 Evan O Connor 2 Christian Ott 1 1 TAPIR, California Institute of Technology 2 CITA, University
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationLagrangian Formulation of Elastic Wave Equation on Riemannian Manifolds
RWE-C3-EAFIT Lagrangian Formulation of Elastic Wave Equation on Riemannian Manifolds Hector Roman Quiceno E. Advisors Ph.D Jairo Alberto Villegas G. Ph.D Diego Alberto Gutierrez I. Centro de Ciencias de
More informationNumerical Simulations. Duncan Christie
Numerical Simulations Duncan Christie Motivation There isn t enough time to derive the necessary methods to do numerical simulations, but there is enough time to survey what methods and codes are available
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationThe MRI in a Collisionless Plasma
The MRI in a Collisionless Plasma Eliot Quataert (UC Berkeley) Collaborators: Prateek Sharma, Greg Hammett, Jim Stone Modes of Accretion thin disk: energy radiated away (relevant to star & planet formation,
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
More information3.5 Vorticity Equation
.0 - Marine Hydrodynamics, Spring 005 Lecture 9.0 - Marine Hydrodynamics Lecture 9 Lecture 9 is structured as follows: In paragraph 3.5 we return to the full Navier-Stokes equations (unsteady, viscous
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 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 informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationTime dependent perturbation theory 1 D. E. Soper 2 University of Oregon 11 May 2012
Time dependent perturbation theory D. E. Soper University of Oregon May 0 offer here some background for Chapter 5 of J. J. Sakurai, Modern Quantum Mechanics. The problem Let the hamiltonian for a system
More informationModeling of ELM Dynamics for ITER
Modeling of ELM Dynamics for ITER A.Y. PANKIN 1, G. BATEMAN 1, D.P. BRENNAN 2, A.H. KRITZ 1, S. KRUGER 3, P.B. SNYDER 4 and the NIMROD team 1 Lehigh University, 16 Memorial Drive East, Bethlehem, PA 18015
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.323: Relativistic Quantum Field Theory I PROBLEM SET 2
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.323: Relativistic Quantum Field Theory I PROBLEM SET 2 REFERENCES: Peskin and Schroeder, Chapter 2 Problem 1: Complex scalar fields Peskin and
More information