The collision probability method in 1D part 1
|
|
- Arleen Townsend
- 5 years ago
- Views:
Transcription
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 1 1/28
2 Content (week 8) 1 The collision probability method Scattering reduction and power iteration The Matlab scriptaleig Calculation of 1D collision probabilities slab geometry ENE6101: Week 8 The collision probability method in 1D part 1 2/28
3 The collision probability method 1 The collision probability (CP) method results from the spatial discretization of the integral transport equation in multigroup form, assuming isotropic particle sources. For a problem with I regions, produces a I I matrix in each energy group. Preferred for treating general unstructured meshes and few-region problems. The integral infinite-domain form of the transport equation is (1) φ g (r,ω) = 0 ds e τ g(s) Q g (r sω,ω) Integrating over the solid angles, we directly obtain the integrated flux φ g (r): (2) φ g (r) = 4π d 2 Ω φ g (r,ω) = 1 4π where the optical path τ g (s) is given by 4π d 2 Ω 0 ds e τ g(s) Q g (r sω) (3) τ g (s) = s 0 ds Σ g (r s Ω). ENE6101: Week 8 The collision probability method in 1D part 1 3/28
4 The collision probability method 2 r' s r Ω We now introduce the change of variable r = r sω with d 3 r = s 2 d 2 Ω ds. We obtain (4) φ g (r) = 1 4π d 3 r e τ g(s) s 2 Q g (r ) with s = r r. Consider an infinite lattice of unit cells, each of them represented as i V i V i represents the infinite set of regions V i belonging to all the cells in the lattice The sources of secondary neutrons are uniform and equal to Q i,g on each region V i. After multiplication by Σ g (r) and integration over each region V i, Eq. (4) can be written (5) V j d 3 rσ g (r)φ g (r) = 1 4π V j d 3 r Σ g (r) i Q i,g V i d 3 r e τ g(s) s 2 ENE6101: Week 8 The collision probability method in 1D part 1 4/28
5 The collision probability method 3 The neutron source in Eq. (5) is defined as (6) Q i,g = h Σ s0,i,g h φ i,h + 1 K eff Q fiss i,g. The fission source in Eq. (6) is defined as (7) fiss J Q fiss i,g = j=1 χ j,g νσ f,j,h φ i,h h where j is a fissionable isotope index, χ j,g is the fission spectrum of isotope j and Σ f,j,h is the macroscopic fission cross section of isotope j for neutrons in group h. Equation (5) simplifies to (8) V j Σ j,g φ j,g = i Q i,g V i P ij,g ENE6101: Week 8 The collision probability method in 1D part 1 5/28
6 The collision probability method 4 where (9) φ j,g = 1 V j V j d 3 rφ g (r), (10) Σ j,g = 1 V j φ j,g V j d 3 rσ g (r)φ g (r) and (11) P ij,g = 1 4πV i V i d 3 r V j d 3 rσ g (r) e τg(s) s 2. The collision probabilityp ij,g is the probability for a neutron born uniformally and isotropically in any of the regions V i of the lattice to undergo its first collision in region V j of a unit cell or assembly. If the total cross section Σ g (r) is constant and equal to Σ j,g in region V j, reduced CPs are (12) p ij,g = P ij,g Σ j,g = 1 4πV i V i d 3 r V j d 3 r e τg(s) s 2. ENE6101: Week 8 The collision probability method in 1D part 1 6/28
7 The collision probability method 5 Reduced CPs generally remain finite in the limit where Σ j,g tends to zero. This ensures the correct behaviour of the collision probability theory in cases where some regions of the lattice are voided. Other interesting properties of CPs are the reciprocity and conservation properties, which can be written (13) and (14) p ij,g V i = p ji,g V j p ij,g Σ j,g = 1 ; i. j Using the reciprocity property, Eq. (8) can be further simplified to (15) φ i,g = j Q j,g p ij,g. ENE6101: Week 8 The collision probability method in 1D part 1 7/28
8 The collision probability method 6 The CP method generally proceeds in three steps: 1. A tracking process is applied over the lattice geometry to span a sufficiently large number of neutron trajectories. In a 2D domain, the tracking parameters are the number of azimuthal angles and the number of parallel tracks per centimetre. The tracking is generally independent of the energy group. 2. A numerical integration process is required to compute the CPs, using tracking information and knowledge of the macroscopic total cross sections in each region. This integration should be repeated for each energy group. 3. Once the CPs are known, the integrated flux can be computed from Eqs. (6) and (15). Collision probability techniques can also be applied to the case of a domain D surrounded by a surface D. The free path lengths are restricted to finite lengths defined inside D. Infinite lattices can still be described by incorporating reflective or periodic boundary conditions over D. ENE6101: Week 8 The collision probability method in 1D part 1 8/28
9 Power iteration 1 The within-group scattering term is first included in the left-hand side of Eq. (15) to obtain (16) φ i,g j p ij,g Σ s0,j,g g φ j,g = j Q j,g p ij,g where Q i,g includes the fission sources and the diffusion sources from all groups except group g. It is obtained from Eq. (6) as (17) Q i,g = h g Σ s0,i,g h φ i,h + 1 K eff Q fiss i,g. Equation (16) can be written in matrix form as (18) Φ g = W g Q g where Φ g = {φ i,g ; i} and Q g = {Q i,g ; i}. ENE6101: Week 8 The collision probability method in 1D part 1 9/28
10 Power iteration 2 The W g matrix is the scattering-reduced collision probability matrix. It is defined as (19) W g = [I P g S s0,g g ] 1 P g where I is the identity matrix, P g = {p ij,g ; i and j} and S s0,g g = diag{σ s0,i,g g ; i}. Two iterative processes are generally superimposed on these 1 G flux solution methods: 1. The inner process iterates over the diffusion up-scattering sources until convergence of the thermal flux. This iteration process is accelerated using two different techniques. The rebalancing technique produces a group-dependent factor which restores the exact multigroup flux distribution homogenized over all regions. The variational acceleration technique. 2. The outer (or power) iteration process is over the eigenvalue and is not required for fixed source problems. It is generally not accelerated in a lattice code and consists in computing neutron flux at outer iteration k +1 from source at iteration k, using (20) Φ (k+1) g = W g Q (k) g The critical parameter (K eff ) is then adjusted at the end of each power iteration. ENE6101: Week 8 The collision probability method in 1D part 1 10/28
11 The Matlab scriptaleig 1 The numerical solution of an eigenvalue equation is possible with Matlab scriptaleig. This script is called as [iter,eval,evect]=aleig(a,b,eps) ; This script find the fondamental eigenvalue and corresponding eigenvector of equation using the inverse power method. Dummy variables are defined as follows: ( A 1 ) B Φ = 0 K eff variablesaandbare the input matrices eps is the convergence parameter of the inverse power method. The script returns a list containing the number of iterations the fundamental eigenvalue 1/K eff the fundamental eigenvector Φ. ENE6101: Week 8 The collision probability method in 1D part 1 11/28
12 Examples of 2D geometries insybilt: 1 Only compatible with the collision probability method. 2D Physical" cell examples 2D assembly example (use the interface current method) ENE6101: Week 8 The collision probability method in 1D part 1 12/28
13 Example of 2D geometries innxt: 1 Compatible with the CP, MOC and Monte-Carlo methods in 2D Can generate prismatic 3D geometries ENE6101: Week 8 The collision probability method in 1D part 1 13/28
14 Example of 3D geometries innxt: 1 Compatible with the CP, MOC and Monte-Carlo methods in 3D z x y ENE6101: Week 8 The collision probability method in 1D part 1 14/28
15 Example of 2D geometries in SALT: 1 mesh-splitting for flat-source approximation same geometry without mesh-splitting ENE6101: Week 8 The collision probability method in 1D part 1 15/28
16 Example of 2D geometries in SALT: 1 mesh-splitting for flat-source approximation same geometry without mesh-splitting ENE6101: Week 8 The collision probability method in 1D part 1 16/28
17 Example of 2D tracking 1 tracking with two values of DELX Dir= 0 Plane=0 IU= 0 IV= 0 Dir= 0 Plane=0 IU= 0 IV= y 2 y x x DELX = 0.2 DELX = 0.01 ENE6101: Week 8 The collision probability method in 1D part 1 17/28
18 Slab geometry 1 The CP formulation in one-dimensional (1D) slab geometry makes possible the analytical integration of Eq. (12) with respect to some dependent variables. We will first consider the case of a geometry made of an infinite lattice of identical unit cells. Each unit cell is made of a succession of I regions, each of volume V i with 1 i I. The reference unit cell is defined in x 1/2 x x I+1/2 and the reference region i is defined in x i 1/2 x x i+1/2. C i is the infinite set of all indices i, representing instances of a region repeating itself by translation of the unit cell ρ x i-1/2 r' ε x i+1/2 s x i+3/2 θ r Ω x x i region i x i+1 region i+1 ENE6101: Week 8 The collision probability method in 1D part 1 18/28
19 Slab geometry 2 Equation (12) is the expression of a CP component in the infinite domain case. It can be rewritten in a convenient form. r = r +sω and d 3 r = s 2 d 2 Ωds, so that (21) p ij = 1 4πV i V i d 3 r 4π d 2 Ω dse τ(s) I j where we have omitted the energy group index in order to simplify the equations. The quantity I j is the set of points belonging simultaneously to to the half straight line of origin r and direction Ω; to a single instance of volume V j. We also define the direction Ω of the particle in term of the colatitude θ and azimuth ǫ using Ω = cosθi+sinθ cosǫj +sinθ sinǫk. With this definition and with the help of figure, we can write (22) d 2 Ω = dθdǫ sinθ l(ρ) = τ(s) cosθ dx = ds cosθ where ρ is a positive number representing the projection of s on the x axis. The material properties are independent of ǫ so that the integration in this variable can be performed analytically. ENE6101: Week 8 The collision probability method in 1D part 1 19/28
20 Slab geometry 3 Equation (21) can be simplified as (23) where p ij = 1 2 x i i C i xi +1/2 π/2 xj+1/2 l(x,x) dx dθ tanθ dx e cosθ x i 1/2 0 x j 1/2 (24) x i = V i = x i+1/2 x i 1/2 andl(x,x) = l(x,x ) is the projected optical path l(ρ) defined between pointsx andxas (25) l(x,x) = x x dx Σ(x ) so that l(x,x j+1/2 ) = l(x,x j 1/2 )+ x j Σ j. Summation over Vi includes instances of V i belonging to all unit cells of the lattice. The presence of neighboring cells will be taken into account in a specular way. ENE6101: Week 8 The collision probability method in 1D part 1 20/28
21 Slab geometry 4 Equation (23) can be rewritten so as to explicitly represent all contributions from other cells: (26) p ij = xi+1/2 π/2 1 dx dθ tanθ 2 x i m=0 x i 1/2 0 dx e ml cell +l(x,x) cosθ +e (m+1)l cell l(x,x) cosθ x j 1/2 xj+1/2 where 1 i < j I and where the unit-cell optical path l cell is defined as (27) l cell = xi+1/2 x 1/2 dx Σ(x ). ENE6101: Week 8 The collision probability method in 1D part 1 21/28
22 Slab geometry 5 The self-region collision probability p ii is (28) p ii = x i x m=0 x i 1/2 dx xi+1/2 x dx xi+1/2 dx x i 1/2 π/2 0 dθ tanθ e ml cell +l(x,x ) cosθ +e (m+1)l cell l(x,x ) cosθ e ml cell +l(x,x) cosθ +e (m+1)l cell l(x,x) cosθ where 1 i I. ENE6101: Week 8 The collision probability method in 1D part 1 22/28
23 Exponential functions 1 Exponential functions are defined as (29) E n (x) = 1 0 duu n 2 e x/u = 1 duu n e xu with n The exponential functions satisfy relations 1 E 1 (x) x x due n (u) = E n+1 (x) E n+1 (x ) 0.5 E 2 (x) E 3 (x) and E n (x) = d dx E n+1(u) E 4 (x) x They can be evaluated efficiently using the matlab scripttaben. ENE6101: Week 8 The collision probability method in 1D part 1 23/28
24 Slab geometry cont n 1 We will now present the Kavenoky method for obtaining the numerical values of the CP components in slab geometry. The principle of this method is to keep the sum over n in Eqs. (26) and (28) and to introduce exponential functions E n (x), with n 1. Integration of Eq. (26) in θ leads to (30) p ij = xi+1/2 1 dx 2 x i m=0 x i 1/2 xj+1/2 [ [ ]} dx {E 1 ml cell +l(x,x) ]+E 1 (m+1)l cell l(x,x) x j 1/2 if i j and integration of Eq. (28) in θ leads to (31) p ii = x i { x m=0 x i 1/2 dx xi+1/2 x dx xi+1/2 x i 1/2 dx [ [ ]} {E 1 ml cell +l(x,x ) ]+E 1 (m+1)l cell l(x,x ) [ [ {E 1 ml cell +l(x,x) ]+E 1 (m+1)l cell l(x,x)]} }. ENE6101: Week 8 The collision probability method in 1D part 1 24/28
25 Slab geometry cont n 2 Equations (30) and (31) involve two types on spatial integration which can be written in functional form as (32) and (33) R ij,m {±l} = 1 2 xi+1/2 x i 1/2 dx R ii,m {±l} = 1 dx 2 x i 1/2 { x xi+1/2 x i 1/2 dxe 1 xj+1/2 x j 1/2 dxe 1 [ ] ml cell ±l(x,x) where i < j [ ] xi+1/2 [ ] } ml cell ±l(x,x ) + dxe 1 ml cell ±l(x,x) x so that (34) p ij = 1 x i m=0 R ij,m {l}+r ij,m+1 { l} where i j. ENE6101: Week 8 The collision probability method in 1D part 1 25/28
26 Slab geometry cont n 3 The collision probabilities defined in Eq. (34) are representing an infinite lattice and are normalized as I (35) p ij Σ j = 1, i. j=1 We will now find analytical reductions for functionals R ij,m {±l} and use them to evaluate the CPs. Integration in x and x is performed, leading to the following relations: a)i < j,σ i 0 andσ j 0 : 1 [ ] R ij,m {±l} = {E 3 ml cell ±l(x i+1/2,x j 1/2 ) 2Σ i Σ j ] E 3 [ml cell ±l(x i+1/2,x j+1/2 ) [ ] + E 3 [ml cell ±l(x i 1/2,x j+1/2 ) ] E } 3 ml cell ±l(x i 1/2,x j 1/2 ) b)i < j,σ i = 0 andσ j 0 : R ij,m {±l} = ± x i 2Σ j [ [ ] {E 2 ml cell ±l(x i 1/2,x j 1/2 ) ] E } 2 ml cell ±l(x i 1/2,x j+1/2 ) ENE6101: Week 8 The collision probability method in 1D part 1 26/28
27 Slab geometry cont n 4 c)i < j,σ i 0 andσ j = 0 : R ij,m {±l} = ± x j 2Σ i d)i < j andσ i = Σ j = 0 : [ [ ] {E 2 ml cell ±l(x i+1/2,x j 1/2 ) ] E } 2 ml cell ±l(x i 1/2,x j 1/2 ) R ij,m {±l} = x i x j 2 ] E 1 [ml cell ±l(x i 1/2,x j 1/2 ) e)i = j andσ i 0 : R ii,m {±l} = ± x i E 2 [ ml cell ] Σ i 1 Σ 2 i [ ] {E 3 [ml cell ] E } 3 ml cell ±l(x i 1/2,x i+1/2 ) f)i = j andσ i = 0 : R ii,m {±l} = E 1 [ ml cell ] x 2 i 2. ENE6101: Week 8 The collision probability method in 1D part 1 27/28
28 Matlab scripts 1 function f=rij_f(sg,tau0,sigi,sigj,segi,segj) if sigi = 0 && sigj = 0 f=0.5*(taben(3,tau0)-taben(3,tau0+sg*sigi*segi)-... taben(3,tau0+sg*sigj*segj)+... taben(3,tau0+sg*sigi*segi+sg*sigj*segj))/(sigi*sigj) ; elseif sigi == 0 && sigj = 0 f=sg*0.5*segi*(taben(2,tau0)-taben(2,tau0+sg*sigj*segj))/sigj ; elseif sigi = 0 && sigj == 0 f=sg*0.5*segj*(taben(2,tau0)-taben(2,tau0+sg*sigi*segi))/sigi ; else f=0.5*segi*segj*taben(1,tau0) ; end function f=rii_f(sg,tau0,sigi,segi) if sigi = 0 f=sg*segi*taben(2,tau0)/sigi-(taben(3,tau0)-... taben(3,tau0+sg*sigi*segi))/sigiˆ2 ; else f=0.5*segiˆ2*taben(1,tau0) ; end ENE6101: Week 8 The collision probability method in 1D part 1 28/28
The collision probability method in 1D part 2
The collision probability method in 1D part 2 Alain Hébert alain.hebert@polymtl.ca Institut de génie nucléaire École Polytechnique de Montréal ENE611: Week 9 The collision probability method in 1D part
More informationSpace-time kinetics. Alain Hébert. Institut de génie nucléaire École Polytechnique de Montréal.
Space-time kinetics Alain Hébert alain.hebert@polymtl.ca Institut de génie nucléaire École Polytechnique de Montréal ENE613: Week 1 Space-time kinetics 1/24 Content (week 1) 1 Space-time kinetics equations
More informationNumerical methods part 2
Numerical methods part 2 Alain Hébert alain.hebert@polymtl.ca Institut de génie nucléaire École Polytechnique de Montréal ENE6103: Week 6 Numerical methods part 2 1/33 Content (week 6) 1 Solution of an
More informationA PWR ASSEMBLY COMPUTATIONAL SCHEME BASED ON THE DRAGON V4 LATTICE CODE
Supercomputing in Nuclear Applications (M&C + SNA 2007) Monterey, California, April 15-19, 2007, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2007) A PWR ASSEMBLY COMPUTATIONAL SCHEME BASED
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 informationA Hybrid Stochastic Deterministic Approach for Full Core Neutronics Seyed Rida Housseiny Milany, Guy Marleau
A Hybrid Stochastic Deterministic Approach for Full Core Neutronics Seyed Rida Housseiny Milany, Guy Marleau Institute of Nuclear Engineering, Ecole Polytechnique de Montreal, C.P. 6079 succ Centre-Ville,
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 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 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 informationLecture 18 Neutron Kinetics Equations
24.505 Lecture 18 Neutron Kinetics Equations Prof. Dean Wang For a nuclear reactor to operate at a constant power level, the rate of neutron production via fission reactions should be exactly balanced
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 informationLesson 8: Slowing Down Spectra, p, Fermi Age
Lesson 8: Slowing Down Spectra, p, Fermi Age Slowing Down Spectra in Infinite Homogeneous Media Resonance Escape Probability ( p ) Resonance Integral ( I, I eff ) p, for a Reactor Lattice Semi-empirical
More informationA Reformulation of the Transport-Transport SPH Equivalence Technique
Ottawa Marriott Hotel, Ottawa, Ontario, Canada, October 8-, 05 A Reformulation of the Transport-Transport SPH Equivalence Technique A Hébert École Polytechnique de Montréal, Montréal, QC, Canada alainhebert@polymtlca
More informationUSE OF LATTICE CODE DRAGON IN REACTOR CALUCLATIONS
USE OF LATTICE CODE DRAGON IN REACTOR CALUCLATIONS ABSTRACT Dušan Ćalić ZEL-EN razvojni center Hočevarjev trg 1 Slovenia-SI8270, Krško, Slovenia dusan.calic@zel-en.si Andrej Trkov, Marjan Kromar J. Stefan
More informationLesson 9: Multiplying Media (Reactors)
Lesson 9: Multiplying Media (Reactors) Laboratory for Reactor Physics and Systems Behaviour Multiplication Factors Reactor Equation for a Bare, Homogeneous Reactor Geometrical, Material Buckling Spherical,
More informationENHANCEMENT OF COMPUTER SYSTEMS FOR CANDU REACTOR PHYSICS SIMULATIONS
ENHANCEMENT OF COMPUTER SYSTEMS FOR CANDU REACTOR PHYSICS SIMULATIONS E. Varin, M. Dahmani, W. Shen, B. Phelps, A. Zkiek, E-L. Pelletier, T. Sissaoui Candu Energy Inc. WORKSHOP ON ADVANCED CODE SUITE FOR
More informationA Cumulative migration method for computing rigorous transport cross sections and diffusion coefficients for LWR lattices with Monte Carlo
A Cumulative migration method for computing rigorous transport cross sections and diffusion coefficients for LWR lattices with Monte Carlo The MIT Faculty has made this article openly available. Please
More informationChapter 22. Dr. Armen Kocharian. Gauss s Law Lecture 4
Chapter 22 Dr. Armen Kocharian Gauss s Law Lecture 4 Field Due to a Plane of Charge E must be perpendicular to the plane and must have the same magnitude at all points equidistant from the plane Choose
More informationEnergy Dependence of Neutron Flux
Energy Dependence of Neutron Flux B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec. 2015 September 1 Contents We start the discussion of the energy
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 informationOn the use of SERPENT code for few-group XS generation for Sodium Fast Reactors
On the use of SERPENT code for few-group XS generation for Sodium Fast Reactors Raquel Ochoa Nuclear Engineering Department UPM CONTENTS: 1. Introduction 2. Comparison with ERANOS 3. Parameters required
More informationUncertainty quantification using SCALE 6.2 package and GPT techniques implemented in Serpent 2
6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy September 26 th -30 th, 2016 Uncertainty quantification using SCALE 6.2 package and GPT techniques implemented in Serpent
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More information1. Introductory Examples
1. Introductory Examples We introduce the concept of the deterministic and stochastic simulation methods. Two problems are provided to explain the methods: the percolation problem, providing an example
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 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 informationNuclear Reactor Physics I Final Exam Solutions
.11 Nuclear Reactor Physics I Final Exam Solutions Author: Lulu Li Professor: Kord Smith May 5, 01 Prof. Smith wants to stress a couple of concepts that get people confused: Square cylinder means a cylinder
More information8: Source-Sink Problems in 1 Energy Group
8: Source-Sink Problems in 1 Energy Group B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 015 Sept.-Dec. 015 September 1 Contents Solving the 1-group diffusion
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational
More informationCross Section Generation Strategy for High Conversion Light Water Reactors Bryan Herman and Eugene Shwageraus
Cross Section Generation Strategy for High Conversion Light Water Reactors Bryan Herman and Eugene Shwageraus 1 Department of Nuclear Science and Engineering Massachusetts Institute of Technology 77 Massachusetts
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 informationLecture 20 Reactor Theory-V
Objectives In this lecture you will learn the following We will discuss the criticality condition and then introduce the concept of k eff.. We then will introduce the four factor formula and two group
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More information51. General Surface Integrals
51. General urface Integrals The area of a surface in defined parametrically by r(u, v) = x(u, v), y(u, v), z(u, v) over a region of integration in the input-variable plane is given by d = r u r v da.
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 informationRobot Control Basics CS 685
Robot Control Basics CS 685 Control basics Use some concepts from control theory to understand and learn how to control robots Control Theory general field studies control and understanding of behavior
More informationABSTRACT 1 INTRODUCTION
A NODAL SP 3 APPROACH FOR REACTORS WITH HEXAGONAL FUEL ASSEMBLIES S. Duerigen, U. Grundmann, S. Mittag, B. Merk, S. Kliem Forschungszentrum Dresden-Rossendorf e.v. Institute of Safety Research P.O. Box
More informationBackground. Background. C. T. Kelley NC State University tim C. T. Kelley Background NCSU, Spring / 58
Background C. T. Kelley NC State University tim kelley@ncsu.edu C. T. Kelley Background NCSU, Spring 2012 1 / 58 Notation vectors, matrices, norms l 1 : max col sum... spectral radius scaled integral norms
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 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 informationEigenvalues and Eigenvectors
Eigenvalues and Eigenvectors Definition 0 Let A R n n be an n n real matrix A number λ R is a real eigenvalue of A if there exists a nonzero vector v R n such that A v = λ v The vector v is called an eigenvector
More informationBeer-Lambert (cont.)
The Beer-Lambert Law: Optical Depth Consider the following process: F(x) Absorbed flux df abs F(x + dx) Scattered flux df scat x x + dx The absorption or scattering of radiation by an optically active
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 informationPoint Load Generated Surface Waves in a Transversely Isotropic Half-Space using Elastodynamic Reciprocity
Advances in Copyright 4 Tech Science Press 5 Point Load Generated Surface Waves in a Transversely Isotropic Half-Space using Elastodynamic eciprocity A. C. Wijeyewickrema 1, M. Hattori 1, S. Leungvichcharoen
More informationLeakage-Corrected Discontinuity Factors for a Second-Generation Th-Pu Pressure-Tube SCWR. Katarzyna Carisse
Leaage-Corrected Discontinuity Factors for a Second-Generation Th-Pu Pressure-Tube SCWR by Katarzyna Carisse A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree of Master of Applied
More informationStandard Finite Elements and Weighted Regularization
Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the
More informationExtension of the MCBEND Monte Carlo Code to Perform Adjoint Calculations using Point Energy Data
Extension of the MCBEND Monte Carlo Code to Perform Adjoint Calculations using Point Energy Data Malcolm Grimstone Abstract In radiation transport calculations there are many situations where the adjoint
More informationContinuous Energy Neutron Transport
Continuous Energy Neutron Transport Kevin Clarno Mark Williams, Mark DeHart, and Zhaopeng Zhong A&M Labfest - Workshop IV on Parallel Transport May 10-11, 2005 College Station, TX clarnokt@ornl.gov (865)
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 informationPHY492: Nuclear & Particle Physics. Lecture 3 Homework 1 Nuclear Phenomenology
PHY49: Nuclear & Particle Physics Lecture 3 Homework 1 Nuclear Phenomenology Measuring cross sections in thin targets beam particles/s n beam m T = ρts mass of target n moles = m T A n nuclei = n moles
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 informationOn the Chord Length Sampling in 1-D Binary Stochastic Media
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 informationQuantum Dynamics. March 10, 2017
Quantum Dynamics March 0, 07 As in classical mechanics, time is a parameter in quantum mechanics. It is distinct from space in the sense that, while we have Hermitian operators, X, for position and therefore
More informationAppendix A. Longwave radiation in a canopy stand
213 Appendix A. ongwave radiation in a canopy stand The longwave radiation transfer equations can be written as (Ross, 1981): di ld (,r) d cosθ = G(,r)i ld (,r)+η,r ( ) (A.1) di lu ( ),r d cosφ = G,r (
More informationCASMO-5/5M Code and Library Status. J. Rhodes, K. Smith, D. Lee, Z. Xu, & N. Gheorghiu Arizona 2008
CASMO-5/5M Code and Library Status J. Rhodes, K. Smith, D. Lee, Z. Xu, & N. Gheorghiu Arizona 2008 CASMO Methodolgy Evolution CASMO-3 Homo. transmission probability/external Gd depletion CASMO-4 up to
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 informationGG612 Lecture 3. Outline
GG61 Lecture 3 Strain and Stress Should complete infinitesimal strain by adding rota>on. Outline Matrix Opera+ons Strain 1 General concepts Homogeneous strain 3 Matrix representa>ons 4 Squares of line
More information11. Point Detector Tally
and suppose that we want half the source to be in the range 005 < E < 1 and the other half to be in the range 1 < E < 20 Then the input is SPn -5 a SIn 005 1 20 SBn C 0 5 1 MCNP breaks up the function
More informationBenchmark Calculation of KRITZ-2 by DRAGON/PARCS. M. Choi, H. Choi, R. Hon
Benchmark Calculation of KRITZ-2 by DRAGON/PARCS M. Choi, H. Choi, R. Hon General Atomics: 3550 General Atomics Court, San Diego, CA 92121, USA, Hangbok.Choi@ga.com Abstract - Benchmark calculations have
More informationAC : TEACHING THE SN METHOD: ZERO TO INTERNATIONAL BENCHMARK IN SIX WEEKS
AC 2008-657: TEACHING THE SN METHOD: ZERO TO INTERNATIONAL BENCHMARK IN SIX WEEKS Erich Schneider, University of Texas at Austin Dr. Schneider is an Assistant Professor of Nuclear and Radiation Engineering
More information1.1. Rotational Kinematics Description Of Motion Of A Rotating Body
PHY 19- PHYSICS III 1. Moment Of Inertia 1.1. Rotational Kinematics Description Of Motion Of A Rotating Body 1.1.1. Linear Kinematics Consider the case of linear kinematics; it concerns the description
More informationThe XY-Model. David-Alexander Robinson Sch th January 2012
The XY-Model David-Alexander Robinson Sch. 08332461 17th January 2012 Contents 1 Introduction & Theory 2 1.1 The XY-Model............................... 2 1.2 Markov Chains...............................
More informationDemonstration of Full PWR Core Coupled Monte Carlo Neutron Transport and Thermal-Hydraulic Simulations Using Serpent 2/ SUBCHANFLOW
Demonstration of Full PWR Core Coupled Monte Carlo Neutron Transport and Thermal-Hydraulic Simulations Using Serpent 2/ SUBCHANFLOW M. Daeubler Institute for Neutron Physics and Reactor Technology (INR)
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 informationHydrodynamic Limit with Geometric Correction in Kinetic Equations
Hydrodynamic Limit with Geometric Correction in Kinetic Equations Lei Wu and Yan Guo KI-Net Workshop, CSCAMM University of Maryland, College Park 2015-11-10 1 Simple Model - Neutron Transport Equation
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 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 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 informationChapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set
Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless
More informationIterative Solution methods
p. 1/28 TDB NLA Parallel Algorithms for Scientific Computing Iterative Solution methods p. 2/28 TDB NLA Parallel Algorithms for Scientific Computing Basic Iterative Solution methods The ideas to use iterative
More informationADVANCED PROGRAMME MATHEMATICS MARKING GUIDELINES
GRADE EXAMINATION NOVEMBER 0 ADVANCED PROGRAMME MATHEMATICS MARKING GUIDELINES Time: hours 00 marks These marking guidelines are prepared for use by examiners and sub-examiners, all of whom are required
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 informationInstructor: Sudipta Mukherji, Institute of Physics, Bhubaneswar
Chapter 1 Lorentz and Poincare Instructor: Sudipta Mukherji, Institute of Physics, Bhubaneswar 1.1 Lorentz Transformation Consider two inertial frames S and S, where S moves with a velocity v with respect
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 informationOrder of Magnitude Astrophysics - a.k.a. Astronomy 111. Photon Opacities in Matter
1 Order of Magnitude Astrophysics - a.k.a. Astronomy 111 Photon Opacities in Matter If the cross section for the relevant process that scatters or absorbs radiation given by σ and the number density of
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 informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More information2 Eigenvectors and Eigenvalues in abstract spaces.
MA322 Sathaye Notes on Eigenvalues Spring 27 Introduction In these notes, we start with the definition of eigenvectors in abstract vector spaces and follow with the more common definition of eigenvectors
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 information7. FREE ELECTRON THEORY.
7. FREE ELECTRON THEORY. Aim: To introduce the free electron model for the physical properties of metals. It is the simplest theory for these materials, but still gives a very good description of many
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 informationContents. 1 Basic Equations 1. Acknowledgment. 1.1 The Maxwell Equations Constitutive Relations 11
Preface Foreword Acknowledgment xvi xviii xix 1 Basic Equations 1 1.1 The Maxwell Equations 1 1.1.1 Boundary Conditions at Interfaces 4 1.1.2 Energy Conservation and Poynting s Theorem 9 1.2 Constitutive
More informationAppendix B. Solutions to Chapter 2 Problems
Appendix B Solutions to Chapter Problems Problem Problem 4 Problem 3 5 Problem 4 6 Problem 5 6 Problem 6 7 Problem 7 8 Problem 8 8 Problem 9 Problem Problem 3 Problem et f ( ) be a function of one variable
More informationCost-accuracy analysis of a variational nodal 2D/1D approach to pin resolved neutron transport
Cost-accuracy analysis of a variational nodal 2D/1D approach to pin resolved neutron transport ZHANG Tengfei 1, WU Hongchun 1, CAO Liangzhi 1, LEWIS Elmer-E. 2, SMITH Micheal-A. 3, and YANG Won-sik 4 1.
More informationAppendix: SU(2) spin angular momentum and single spin dynamics
Phys 7 Topics in Particles & Fields Spring 03 Lecture v0 Appendix: SU spin angular momentum and single spin dynamics Jeffrey Yepez Department of Physics and Astronomy University of Hawai i at Manoa Watanabe
More informationNeutron Diffusion Theory: One Velocity Model
22.05 Reactor Physics - Part Ten 1. Background: Neutron Diffusion Theory: One Velocity Model We now have sufficient tools to begin a study of the second method for the determination of neutron flux as
More informationBohr & Wheeler Fission Theory Calculation 4 March 2009
Bohr & Wheeler Fission Theory Calculation 4 March 9 () Introduction The goal here is to reproduce the calculation of the limiting Z /A against spontaneous fission Z A lim a S. (.) a C as first done by
More informationEfficient Reduced Order Modeling of Low- to Mid-Frequency Vibration and Power Flow in Complex Structures
Efficient Reduced Order Modeling of Low- to Mid-Frequency Vibration and Power Flow in Complex Structures Yung-Chang Tan Graduate Student Research Assistant Matthew P. Castanier Assistant Research Scientist
More informationPROBLEM L. (3) Noting that since the aperture emits diffusely, I e = E/π (see Eq ), and hence
PROBLEM 1.004 KNOWN: Furnace with prescribed aperture and emissive power. FIND: (a) Position of gauge such that irradiation is G = 1000 W/m, (b) Irradiation when gauge is tilted θ d = 0 o, and (c) Compute
More informationEigenvalues and Eigenvectors A =
Eigenvalues and Eigenvectors Definition 0 Let A R n n be an n n real matrix A number λ R is a real eigenvalue of A if there exists a nonzero vector v R n such that A v = λ v The vector v is called an eigenvector
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
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 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 informationBefore you begin read these instructions carefully:
NATURAL SCIENCES TRIPOS Part IB & II (General) Friday, 27 May, 2016 9:00 am to 12:00 pm MATHEMATICS (2) Before you begin read these instructions carefully: You may submit answers to no more than six questions.
More informationSolutions: Homework 7
Solutions: Homework 7 Ex. 7.1: Frustrated Total Internal Reflection a) Consider light propagating from a prism, with refraction index n, into air, with refraction index 1. We fix the angle of incidence
More informationChapter 10. Rotation of a Rigid Object about a Fixed Axis
Chapter 10 Rotation of a Rigid Object about a Fixed Axis Angular Position Axis of rotation is the center of the disc Choose a fixed reference line. Point P is at a fixed distance r from the origin. A small
More informationProblem 1: Spin 1 2. particles (10 points)
Problem 1: Spin 1 particles 1 points 1 Consider a system made up of spin 1/ particles. If one measures the spin of the particles, one can only measure spin up or spin down. The general spin state of a
More informationDocumentation of the Solutions to the SFPE Heat Transfer Verification Cases
Documentation of the Solutions to the SFPE Heat Transfer Verification Cases Prepared by a Task Group of the SFPE Standards Making Committee on Predicting the Thermal Performance of Fire Resistive Assemblies
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
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 information