The collision probability method in 1D part 2
|
|
- Easter Barnett
- 5 years ago
- Views:
Transcription
1 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 2 1/25
2 Content (week 9) 1 Calculation of 1D collision probabilities plane 2D geometry cylindrical 1D geometry spherical 1D geometry ENE611: Week 9 The collision probability method in 1D part 2 2/25
3 The collision probability method 1 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 on each region V i. 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 Σ(r) is constant and equal to Σ j in region V j, reduced CPs are (1) p ij = P ij Σ j = 1 4πV i V i d 3 r d 3 r e τ(s) V j s 2 where the optical path τ(s) is given by (2) τ(s) = s ds Σ(r s Ω). ENE611: Week 9 The collision probability method in 1D part 2 3/25
4 Plane 2D geometry 1 We will now study the case of a 2D geometry defined in the x y plane and homogeneous along the z axis. A cylinder or tube of infinite length is an example of such a geometry. In this case, it is possible to integrate analytically the integral transport equation along the axial angle θ. This angle is use to define the direction cosine of the particle relative to the z axis. z Ω r = (x,y,z) The direction of the particle is Ω = sinθ cosǫi+sinθ sinǫj +cosθk. Integration over θ is possible by taking the projection of each particle free path on the x y plane. We write θ y d 2 Ω = dθdǫ sinθ x ε (x,y) (3) τ(ρ) = τ(s) sinθ dρ = ds sinθ where ρ is the projection of s on the x y plane. ENE611: Week 9 The collision probability method in 1D part 2 4/25
5 Plane 2D geometry 2 We will limit ourself on the integration of the integral transport equation, defined on the domain of a unique unit cell. We use the relations r = r +sω and d 3 r = s 2 d 2 Ωds and obtain (4) p ij = 1 4πV i d 3 r V i 4π d 2 Ω dse τ(s) I j where the quantity I j is the set of points belonging simultaneously to the half straight line of origin r and direction Ω; to volume V j of the unit cell. With the help of figure, Eq. (4) can be rewritten as (5) p ij = 1 4πV i 2π dǫ d 2 r V i π dθ dρe I j τ(ρ) sinθ where V i now represents a surface in the x y plane and I j is the set of points belonging simultaneously to the half straight line of origin r and direction ǫ; to surface V j of the unit cell. ENE611: Week 9 The collision probability method in 1D part 2 5/25
6 Bickley functions 1 Bickley functions are defined as (6) with n 1. Ki n (x) = π/2 dθ sin n 1 θe x sinθ = π/2 dθ cos n 1 θe x cosθ Ki1(x) Ki2(x) Ki3(x) The following equations are relating Bickley functions of different orders: x x du Ki n (u) = Ki n+1 (x) Ki n+1 (x ),.5 Ki4(x) Ki5(x) and d dx Ki n(x) = Ki n 1 (x) x They can be evaluated efficiently using the matlab scriptakin. ENE611: Week 9 The collision probability method in 1D part 2 6/25
7 Plane 2D geometry 1 Equation (5) can be simplified by introducing the Bickley functions. We obtain (7) p ij = 1 2πV i 2π dǫ d 2 r dρki 1 [τ(ρ)]. V i I j Sets of integration lines, referred as tracks, are draw over the complete domain. Each set is characterized by a given angle ǫ and contains parallel tracks covering the domain. Tracks in a set can be separated by a constant distance h or can be placed at optimal locations, as depicted in sub-figures (a) and (b), respectively. Each track is used forward and backward, corresponding to angles ǫ and ǫ. Three-point formula h Three-point formula Three-point formula (a) Rectangular quadrature (b) Gaussian quadrature ENE611: Week 9 The collision probability method in 1D part 2 7/25
8 Plane 2D geometry 2 In the case of convex volumes i and j, Eq. (7) can be rewritten as (8) p ij = 1 2πV i 2π dǫ hmax h min dh lj dl dl Ki 1 (τ ij +Σ i l +Σ j l) if i j where τ ij is the optical path of the materials between regions i and j and (9) p ii = 1 2πV i 2π dǫ hmax h min dh dl l dl Ki 1 [Σ i (l l )]. y V i l l i l ij l j l' r V j r' h ε x ENE611: Week 9 The collision probability method in 1D part 2 8/25
9 Cylindrical 1D geometry 1 In the case of 1D cylindrical geometry, the tracks are identical for any value of the angle ǫ. track 1 track 2 The tubular volumes are concave, making extra terms to appear. The corresponding geometry is depicted in figure. Two tracks are represented, both contributing to collision probability components p ij with i < j and p ii. Track 1 corresponds to the integration domain where region i is concave and track 2 corresponds to the integration domain where it is convex. In this figure, τ ij and τ ii are the optical paths of the materials located between regions i and j or between two instances of region i. Vj Vi l j l j τij l i τii r j-1/2 r l i-1/2 i r i+1/2 l i l j l j r j+1/2 ENE611: Week 9 The collision probability method in 1D part 2 9/25
10 Cylindrical 1D geometry 2 In this case, Eq. (8) can be rewritten as (1) p ij { = 2 ri 1/2 dh V i + lj dl dl ] + Ki 1 (τ ij +Σ i l +Σ j l) ri+1/2 r i 1/2 dh [ Ki 1 [(τ ij +τ ii +Σ i l i )+Σ i l +Σ j l] } lj dl dl Ki 1 (τ ij +Σ i l +Σ j l) if i < j where r i±1/2 are the radii bounding region i and Eq. (9) can be rewritten as { p ii = 2 ri 1/2 dh V i (11) + ri+1/2 r i 1/2 dh dl dl [2 dl Ki 1 [Σ i (l l )]+ l } l dl Ki 1 [Σ i (l l )] dl Ki 1 [τ ii +Σ i (l +l)] where the first and second terms are the contributions from tracks 1 and 2, respectively. ] ENE611: Week 9 The collision probability method in 1D part 2 1/25
11 Cylindrical 1D geometry 3 These equations can be simplified by introducing the following definitions: (12) C ij (τ ) = lj dl dl Ki 1 (τ +Σ i l +Σ j l) and (13) so that (14) D i = dl dl Ki 1 [Σ i (l l )] l { p ij = 2 ri 1/2 dh[c ij (τ ij +τ ii +Σ i l i )+C ij (τ ij )]+ V i if i < j, and ri+1/2 r i 1/2 dhc ij (τ ij ) } (15) { p ii = 2 ri 1/2 dh[2d i +C ii (τ ii )]+ V i ri+1/2 r i 1/2 dhd i }. ENE611: Week 9 The collision probability method in 1D part 2 11/25
12 Cylindrical 1D geometry 4 Integration in l and l is next performed, leading to the following relations: a)σ i andσ j : 1 [ C ij (τ ) = Ki 3 (τ ) Ki 3 (τ +Σ i l i ) Ki 3 (τ +Σ j l j ) Σ i Σ j ] + Ki 3 (τ +Σ i l i +Σ j l j ) b)σ i = andσ j : c)σ i andσ j = : C ij (τ ) = l i Σ j [Ki 2 (τ ) Ki 2 (τ +Σ j l j )] C ij (τ ) = l j Σ i [Ki 2 (τ ) Ki 2 (τ +Σ i l i )] d)σ i = Σ j = : C ij (τ ) = l i l j Ki 1 (τ ) ENE611: Week 9 The collision probability method in 1D part 2 12/25
13 Cylindrical 1D geometry 5 e)σ i : f)σ i = : D i = l i Σ i 1 Σ 2 i [Ki 3 () Ki 3 (Σ i l i )] D i = πl2 i 4. ENE611: Week 9 The collision probability method in 1D part 2 13/25
14 Cylindrical 1D geometry 6 The above relations are implemented in Matlab using the following two scripts: function f=cij_f(tau,sigi,sigj,segmenti,segmentj) if sigi = && sigj = f=(akin(3,tau)-akin(3,tau+sigi*segmenti)-... akin(3,tau+sigj*segmentj)+... akin(3,tau+sigi*segmenti+sigj*segmentj))/(sigi*sigj) ; elseif sigi == && sigj = f=(akin(2,tau)-akin(2,tau+sigj*segmentj))*segmenti/sigj ; elseif sigi = && sigj == f=(akin(2,tau)-akin(2,tau+sigi*segmenti))*segmentj/sigi ; else f=akin(1,tau)*segmenti*segmentj ; end function f=di_f(sig,segment) if sig = f=segment/sig-(akin(3,)-akin(3,sig*segment))/sigˆ2 ; else f=pi*segmentˆ2/4 ; end ENE611: Week 9 The collision probability method in 1D part 2 14/25
15 Tracking 1 The tracking of 1D cylindrical and spherical geometries is similar and can be generated with the same algorithm. A matlab scriptf=sybt1d(rad,lgsph,ngauss) is presented to produce a tracking object in these cases. The figure presents an example with I = 2 regions and K = 2 tracks per regions. Tracks 1 and 2 have two segments, identified as l 1 and l 2. Tracks 3 and 4 have only one segment, identified as l φ l 2 3 l 2 l 2 l 1 l 1 l 2 h 1 h 3 h 2 h 4 h ENE611: Week 9 The collision probability method in 1D part 2 15/25
16 Tracking 2 A K point Gauss-Jacobi quadrature is used The four tracks are obtained using rad=[.75.15] ; track=sybt1d(rad,false,2) ; track(i,j) is a 2D cell array containing I K elements, one per track. The i index refers to the tracks located in r i 1/2 < h < r i+1/2. The (i,j) th cell is also a cell array containing three elements: the weight w k of the track, the product w k cosφ, a 1D array of dimension I i+1 containing the track segments. The l 1 segment length of track 2 is l2=track{1,2}{3}(1) leading tol2 = Any lengthtrack{ik,il}{3}(1) corresponding to a segment that cross the h axis must be multiplied by 2. 2 l 2 l 1 1 l 2 l 1 4 l 2 φ 3 l 2 h 1 h 3 h 2 h 4 h ENE611: Week 9 The collision probability method in 1D part 2 16/25
17 Tracking 3 The following script will perform a numerical integration of p 11 in 1D cylindrical tube of radius rad and macroscopic total cross sectionsigt: function f=p11_cyl(rad,sigt) % compute the p11 component in 1D cylindrical geometry track=sybt1d(rad,false,5) ; val= ; for i=1:5 segment=2*track{1,i}{3}(1) ; if sigt = d=segment/sigt-(akin(3,)-akin(3,sigt*segment))/sigtˆ2 ; else d=pi*segmentˆ2/4 ; end val=val+track{1,i}{1}*d ; end f=val/(pi*radˆ2) ; ENE611: Week 9 The collision probability method in 1D part 2 17/25
18 Boundary conditions 1 The external boundary condition of a 1D cylindrical geometry is introduced by choosing an albedo β + set to zero for representing a voided boundary, or set to one for representing reflection of particles. The voided boundary condition is similar to the vacuum condition used for 1D slab geometry. The reflecting boundary condition is implemented in the context of the Wigner-Seitz approximation for lattice calculations in reactor physics. The Wigner-Seitz approximation consists of replacing the exact boundary with an equivalent cylindrical boundary, taking care to conserve the amount of moderator present in the cell. The correct condition in this case is the white boundary condition in which the particles are reflected with an isotropic angular distribution. Fuel Fuel ENE611: Week 9 The collision probability method in 1D part 2 18/25
19 Boundary conditions 2 The application of a white boundary condition is based on the transformation of the reduced collision probability matrix P = {p ij, i = 1,I and j = 1,I} corresponding to a vacuum boundary condition. We first compute the escape probability vector P is = {P is, i = 1,I}, a column vector whose components are obtained from (16) I P is = 1 p ij Σ j j=1 representing the probability for a neutron born in region i to escape from the unit cell without collision. The probability for a neutron entering in the unit cell by its boundary, with an isotropic angular distribution, to first collide in region i is represented by the component P Si. It is possible to show that (17) where S is the surface of the boundary. p Si = P Si Σ i = 4V i S P is ENE611: Week 9 The collision probability method in 1D part 2 19/25
20 Boundary conditions 3 The transmission probability P SS is therefore given as (18) P SS = 1 I p Si Σ i. i=1 The reduced collision probability matrix P = { p ij, i = 1,I and j = 1,I} corresponding to a white boundary condition can now be obtained using P = P+P is [ β + +(β + ) 2 P SS +(β + ) 3 P 2 SS +(β+ ) 4 P 3 SS +...] p Sj or, using a geometric progression, as (19) P = P+ β + 1 β + P SS P is p Sj. The collision probability matrix P can be scattering-reduced and used in the matrix flux equation of the CP method. ENE611: Week 9 The collision probability method in 1D part 2 2/25
21 Spherical 1D geometry 1 Computing the tracking and collision probabilities in 1D spherical geometry is similar to the cylindrical geometry case. The tracking is obtained by thesybt1d Matlab script The expression of reduced collision probabilities is written (2) p ij = 1 d 3 r V i V i I j dse τ(s) ENE611: Week 9 The collision probability method in 1D part 2 21/25
22 Spherical 1D geometry 2 which can be rewritten as (21) p ij = 2π V i { ri 1/2 dhh ] + e (τ ij+σ i l +Σ j l) + ri+1/2 r i 1/2 dhh lj dl dl } lj dl dl e (τ ij+σ i l +Σ j l) [ e [(τ ij+τ ii +Σ i l i )+Σ i l +Σ j l] if i < j where r i±1/2 are the radii bounding region i and (22) p ii = 2π V i + { ri 1/2 dhh ri+1/2 r i 1/2 dhh dl dl [2 l dl e [Σ i (l l l dl e [Σ i (l l )] } )] + dl e [τ ii+σ i (l +l)] where the first and second terms are the contributions from tracks 1 and 2, respectively. ] ENE611: Week 9 The collision probability method in 1D part 2 22/25
23 Spherical 1D geometry 3 These equations can be simplified by introducing the following definitions: (23) C ij (τ ) = lj dl dl e (τ +Σ i l +Σ j l) and (24) so that p ij = 2π V i (25) D i = dl l dl e Σ i (l l { ri 1/2 dhh[c ij (τ ij +τ ii +Σ i l i )+C ij (τ ij )]+ ) ri+1/2 r i 1/2 dhhc ij (τ ij ) } if i < j and (26) p ii = 2π V i { ri 1/2 dhh[2d i +C ii (τ ii )]+ ri+1/2 r i 1/2 dhhd i }. ENE611: Week 9 The collision probability method in 1D part 2 23/25
24 Spherical 1D geometry 4 Integration in l and l is next performed, leading to the following relations: a)σ i andσ j : C ij (τ ) = 1 Σ i Σ j [ e τ e (τ +Σ i l i ) e (τ +Σ j l j ) +e (τ +Σ i l i +Σ j l j ) ] b)σ i = andσ j : c)σ i andσ j = : d)σ i = Σ j = : C ij (τ ) = l i Σ j [ e τ e (τ +Σ j l j ) ] C ij (τ ) = l j Σ i [ e τ e (τ +Σ i l i ) ] C ij (τ ) = l i l j e τ ENE611: Week 9 The collision probability method in 1D part 2 24/25
25 Spherical 1D geometry 5 e)σ i : f)σ i = : D i = l i Σ i 1 Σ 2 i [1 e Σ i l i ] D i = l2 i 2. ENE611: Week 9 The collision probability method in 1D part 2 25/25
The 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 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 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 informationDiscrete Euclidean Curvature Flows
Discrete Euclidean Curvature Flows 1 1 Department of Computer Science SUNY at Stony Brook Tsinghua University 2010 Isothermal Coordinates Relation between conformal structure and Riemannian metric Isothermal
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 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 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 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 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 informationSEAFLOOR MAPPING MODELLING UNDERWATER PROPAGATION RAY ACOUSTICS
3 Underwater propagation 3. Ray acoustics 3.. Relevant mathematics We first consider a plane wave as depicted in figure. As shown in the figure wave fronts are planes. The arrow perpendicular to the wave
More informationxy 2 e 2z dx dy dz = 8 3 (1 e 4 ) = 2.62 mc. 12 x2 y 3 e 2z 2 m 2 m 2 m Figure P4.1: Cube of Problem 4.1.
Problem 4.1 A cube m on a side is located in the first octant in a Cartesian coordinate system, with one of its corners at the origin. Find the total charge contained in the cube if the charge density
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 informationElectromagnetic Field Theory (EMT)
Electromagnetic Field Theory (EMT) Lecture # 9 1) Coulomb s Law and Field Intensity 2) Electric Fields Due to Continuous Charge Distributions Line Charge Surface Charge Volume Charge Coulomb's Law Coulomb's
More information1. ELECTRIC CHARGES AND FIELDS
1. ELECTRIC CHARGES AND FIELDS 1. What are point charges? One mark questions with answers A: Charges whose sizes are very small compared to the distance between them are called point charges 2. The net
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 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 informationIntegrals in cylindrical, spherical coordinates (Sect. 15.7)
Integrals in clindrical, spherical coordinates (Sect. 15.7 Integration in spherical coordinates. Review: Clindrical coordinates. Spherical coordinates in space. Triple integral in spherical coordinates.
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 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 informationProblem Solving 1: Line Integrals and Surface Integrals
A. Line Integrals MASSACHUSETTS INSTITUTE OF TECHNOLOY Department of Physics Problem Solving 1: Line Integrals and Surface Integrals The line integral of a scalar function f ( xyz),, along a path C is
More informationChapter 21: Gauss law Tuesday September 13 th. Gauss law and conductors Electrostatic potential energy (more likely on Thu.)
Chapter 21: Gauss law Tuesday September 13 th LABS START THIS WEEK Quick review of Gauss law The flux of a vector field The shell theorem Gauss law for other symmetries A uniformly charged sheet A uniformly
More informationChapter 1 The Electric Force
Chapter 1 The Electric Force 1. Properties of the Electric Charges 1- There are two kinds of the electric charges in the nature, which are positive and negative charges. - The charges of opposite sign
More informationMaxwell s equations for electrostatics
Maxwell s equations for electrostatics October 6, 5 The differential form of Gauss s law Starting from the integral form of Gauss s law, we treat the charge as a continuous distribution, ρ x. Then, letting
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 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 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 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 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 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 informationThe Dancoff correction. in various geometries. /. Carlvik, B. Pershagen
AE-16 The Dancoff correction in various geometries /. Carlvik, B. Pershagen AKTIEBOLAGET ATOMENERGI STOCKHOLM SWEDEN 1959 AE-16 The Dancoff correction in various geometries I. Carhik, B. Persbagen Summary;
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 informationNote: Each problem is worth 14 points except numbers 5 and 6 which are 15 points. = 3 2
Math Prelim II Solutions Spring Note: Each problem is worth points except numbers 5 and 6 which are 5 points. x. Compute x da where is the region in the second quadrant between the + y circles x + y and
More informationPHYS 390 Lecture 23 - Photon gas 23-1
PHYS 39 Lecture 23 - Photon gas 23-1 Lecture 23 - Photon gas What's Important: radiative intensity and pressure stellar opacity Text: Carroll and Ostlie, Secs. 9.1 and 9.2 The temperature required to drive
More informationSENIOR_ 2017_CLASS_12_PHYSICS_ RAPID REVISION_1_ DERIVATIONS IN FIRST FIVE LESSONS Page 1
INDIAN SCHOOL MUSCAT Department of Physics Class XII Rapid Revision -1 DERIVATIONS IN FIRST FIVE LESSONS 1) Field due to an infinite long straight charged wire Consider an uniformly charged wire of infinite
More information3 Chapter. Gauss s Law
3 Chapter Gauss s Law 3.1 Electric Flux... 3-2 3.2 Gauss s Law (see also Gauss s Law Simulation in Section 3.10)... 3-4 Example 3.1: Infinitely Long Rod of Uniform Charge Density... 3-9 Example 3.2: Infinite
More information14.1. Multiple Integration. Iterated Integrals and Area in the Plane. Iterated Integrals. Iterated Integrals. MAC2313 Calculus III - Chapter 14
14 Multiple Integration 14.1 Iterated Integrals and Area in the Plane Objectives Evaluate an iterated integral. Use an iterated integral to find the area of a plane region. Copyright Cengage Learning.
More informationMath 3435 Homework Set 11 Solutions 10 Points. x= 1,, is in the disk of radius 1 centered at origin
Math 45 Homework et olutions Points. ( pts) The integral is, x + z y d = x + + z da 8 6 6 where is = x + z 8 x + z = 4 o, is the disk of radius centered on the origin. onverting to polar coordinates then
More informationInstructions: No books. No notes. Non-graphing calculators only. You are encouraged, although not required, to show your work.
Exam 3 Math 850-007 Fall 04 Odenthal Name: Instructions: No books. No notes. Non-graphing calculators only. You are encouraged, although not required, to show your work.. Evaluate the iterated integral
More informationTutorial on Surface Ricci Flow, Theory, Algorithm and Application
Tutorial on Surface Ricci Flow, Theory, Algorithm and Application 1 1 Department of Computer Science University of New York at Stony Brook Graduate Summer School: Computer Vision, IPAM, UCLA Collaborators
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 informationGraduate Written Examination Spring 2014 Part I Thursday, January 16th, :00am to 1:00pm
Graduate Written Examination Spring 2014 Part I Thursday, January 16th, 2014 9:00am to 1:00pm University of Minnesota School of Physics and Astronomy Examination Instructions Part 1 of this exam consists
More informationIntroduction and Vectors Lecture 1
1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum
More informationDownloaded from
Question 1.1: What is the force between two small charged spheres having charges of 2 10 7 C and 3 10 7 C placed 30 cm apart in air? Repulsive force of magnitude 6 10 3 N Charge on the first sphere, q
More informationAnswer sheet: Final exam for Math 2339, Dec 10, 2010
Answer sheet: Final exam for Math 9, ec, Problem. Let the surface be z f(x,y) ln(y + cos(πxy) + e ). (a) Find the gradient vector of f f(x,y) y + cos(πxy) + e πy sin(πxy), y πx sin(πxy) (b) Evaluate f(,
More informationMultiple Choice. Compute the Jacobian, (u, v), of the coordinate transformation x = u2 v 4, y = uv. (a) 2u 2 + 4v 4 (b) xu yv (c) 3u 2 + 7v 6
.(5pts) y = uv. ompute the Jacobian, Multiple hoice (x, y) (u, v), of the coordinate transformation x = u v 4, (a) u + 4v 4 (b) xu yv (c) u + 7v 6 (d) u (e) u v uv 4 Solution. u v 4v u = u + 4v 4..(5pts)
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 informationf dr. (6.1) f(x i, y i, z i ) r i. (6.2) N i=1
hapter 6 Integrals In this chapter we will look at integrals in more detail. We will look at integrals along a curve, and multi-dimensional integrals in 2 or more dimensions. In physics we use these integrals
More informationGeometrical models for spheroidal cosmological voids
Geometrical models for spheroidal cosmological voids talk by: Osvaldo M. Moreschi collaborator: Ezequiel Boero FaMAF, Universidad Nacional de Córdoba, Instituto de Física Enrique Gaviola (IFEG), CONICET,
More informationUNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 110A. Homework #6. Benjamin Stahl. February 17, 2015
UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS A Homework #6 Benjamin Stahl February 7, 5 GRIFFITHS, 5.9 The magnetic field at a point, P, will be found for each of the steady current
More informationANTENNAS. Vector and Scalar Potentials. Maxwell's Equations. E = jωb. H = J + jωd. D = ρ (M3) B = 0 (M4) D = εe
ANTENNAS Vector and Scalar Potentials Maxwell's Equations E = jωb H = J + jωd D = ρ B = (M) (M) (M3) (M4) D = εe B= µh For a linear, homogeneous, isotropic medium µ and ε are contant. Since B =, there
More informationMath Exam IV - Fall 2011
Math 233 - Exam IV - Fall 2011 December 15, 2011 - Renato Feres NAME: STUDENT ID NUMBER: General instructions: This exam has 16 questions, each worth the same amount. Check that no pages are missing and
More informationOne side of each sheet is blank and may be used as scratch paper.
Math 244 Spring 2017 (Practice) Final 5/11/2017 Time Limit: 2 hours Name: No calculators or notes are allowed. One side of each sheet is blank and may be used as scratch paper. heck your answers whenever
More informationNeutron reproduction. factor ε. k eff = Neutron Life Cycle. x η
Neutron reproduction factor k eff = 1.000 What is: Migration length? Critical size? How does the geometry affect the reproduction factor? x 0.9 Thermal utilization factor f x 0.9 Resonance escape probability
More informationCANONICAL EQUATIONS. Application to the study of the equilibrium of flexible filaments and brachistochrone curves. By A.
Équations canoniques. Application a la recherche de l équilibre des fils flexibles et des courbes brachystochrones, Mem. Acad. Sci de Toulouse (8) 7 (885), 545-570. CANONICAL EQUATIONS Application to the
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 informationPhysics 208, Spring 2016 Exam #3
Physics 208, Spring 206 Exam #3 A Name (Last, First): ID #: Section #: You have 75 minutes to complete the exam. Formulae are provided on an attached sheet. You may NOT use any other formula sheet. You
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 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 informationROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring
ROTATIONAL STRESSES INTRODUCTION High centrifugal forces are developed in machine components rotating at a high angular speed of the order of 100 to 500 revolutions per second (rps). High centrifugal force
More informationSections minutes. 5 to 10 problems, similar to homework problems. No calculators, no notes, no books, no phones. No green book needed.
MTH 34 Review for Exam 4 ections 16.1-16.8. 5 minutes. 5 to 1 problems, similar to homework problems. No calculators, no notes, no books, no phones. No green book needed. Review for Exam 4 (16.1) Line
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 informationMultiple Integrals and Vector Calculus: Synopsis
Multiple Integrals and Vector Calculus: Synopsis Hilary Term 28: 14 lectures. Steve Rawlings. 1. Vectors - recap of basic principles. Things which are (and are not) vectors. Differentiation and integration
More informationMATH 52 FINAL EXAM SOLUTIONS
MAH 5 FINAL EXAM OLUION. (a) ketch the region R of integration in the following double integral. x xe y5 dy dx R = {(x, y) x, x y }. (b) Express the region R as an x-simple region. R = {(x, y) y, x y }
More informationFLUX OF VECTOR FIELD INTRODUCTION
Chapter 3 GAUSS LAW ntroduction Flux of vector field Solid angle Gauss s Law Symmetry Spherical symmetry Cylindrical symmetry Plane symmetry Superposition of symmetric geometries Motion of point charges
More information第 1 頁, 共 8 頁 Chap32&Chap33 1. Test Bank, Question 2 Gauss' law for magnetism tells us: the net charge in any given volume that the line integral of a magnetic around any closed loop must vanish the magnetic
More informationPhysics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017
Physics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
More informationMath 3c Solutions: Exam 1 Fall Graph by eliiminating the parameter; be sure to write the equation you get when you eliminate the parameter.
Math c Solutions: Exam 1 Fall 16 1. Graph by eliiminating the parameter; be sure to write the equation you get when you eliminate the parameter. x tan t x tan t y sec t y sec t t π 4 To eliminate the parameter,
More informationPhys 4322 Final Exam - Solution May 12, 2015
Phys 4322 Final Exam - Solution May 12, 2015 You may NOT use any book or notes other than that supplied with this test. You will have 3 hours to finish. DO YOUR OWN WORK. Express your answers clearly and
More information(You may need to make a sin / cos-type trigonometric substitution.) Solution.
MTHE 7 Problem Set Solutions. As a reminder, a torus with radii a and b is the surface of revolution of the circle (x b) + z = a in the xz-plane about the z-axis (a and b are positive real numbers, with
More informationCompton Storage Rings
Compton Polarimetry @ Storage Rings Wolfgang Hillert ELectron Stretcher Accelerator Physics Institute of Bonn University Møller-Polarimeter Compton-Polarimeter Mott-Polarimeter Compton Scattering Differential
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 informationX. Assembling the Pieces
X. Assembling the Pieces 179 Introduction Our goal all along has been to gain an understanding of nuclear reactors. As we ve noted many times, this requires knowledge of how neutrons are produced and lost.
More informationθ N e θ 0 θ 0 = λ A 0 (θ,φ) 2 sinθdθdφ = 1 (3) A i(θ,φ)a j (θ,φ)sinθdθdφ (4)
1 Introduction 1.1 Angular Spectrum of Plane Waves A Gaussian beam with a beam waist of w 0 has a corresponding angular representation A 0 (θ,φ) = 1 ( ) 2 θ N e θ 0 (1) where and the normalization N is
More informatione x3 dx dy. 0 y x 2, 0 x 1.
Problem 1. Evaluate by changing the order of integration y e x3 dx dy. Solution:We change the order of integration over the region y x 1. We find and x e x3 dy dx = y x, x 1. x e x3 dx = 1 x=1 3 ex3 x=
More information10. Scattering from Central Force Potential
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 215 1. Scattering from Central Force Potential Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative
More informationSolutions to Sample Questions for Final Exam
olutions to ample Questions for Final Exam Find the points on the surface xy z 3 that are closest to the origin. We use the method of Lagrange Multipliers, with f(x, y, z) x + y + z for the square of the
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 informationMATH 241 Practice Second Midterm Exam - Fall 2012
MATH 41 Practice Second Midterm Exam - Fall 1 1. Let f(x = { 1 x for x 1 for 1 x (a Compute the Fourier sine series of f(x. The Fourier sine series is b n sin where b n = f(x sin dx = 1 = (1 x cos = 4
More informationf(p i )Area(T i ) F ( r(u, w) ) (r u r w ) da
MAH 55 Flux integrals Fall 16 1. Review 1.1. Surface integrals. Let be a surface in R. Let f : R be a function defined on. efine f ds = f(p i Area( i lim mesh(p as a limit of Riemann sums over sampled-partitions.
More informationPhysics 221. Exam III Spring f S While the cylinder is rolling up, the frictional force is and the cylinder is rotating
Physics 1. Exam III Spring 003 The situation below refers to the next three questions: A solid cylinder of radius R and mass M with initial velocity v 0 rolls without slipping up the inclined plane. N
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 informationProblem Set #4: 4.1,4.7,4.9 (Due Monday, March 25th)
Chapter 4 Multipoles, Dielectrics Problem Set #4: 4.,4.7,4.9 (Due Monday, March 5th 4. Multipole expansion Consider a localized distribution of charges described by ρ(x contained entirely in a sphere of
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 informationProblem Set 5 Math 213, Fall 2016
Problem Set 5 Math 213, Fall 216 Directions: Name: Show all your work. You are welcome and encouraged to use Mathematica, or similar software, to check your answers and aid in your understanding of the
More informationMoments of Inertia (7 pages; 23/3/18)
Moments of Inertia (7 pages; 3/3/8) () Suppose that an object rotates about a fixed axis AB with angular velocity θ. Considering the object to be made up of particles, suppose that particle i (with mass
More informationHeriot-Watt University
Heriot-Watt University Distinctly Global www.hw.ac.uk Thermodynamics By Peter Cumber Prerequisites Interest in thermodynamics Some ability in calculus (multiple integrals) Good understanding of conduction
More informationPractice Final Solutions
Practice Final Solutions Math 1, Fall 17 Problem 1. Find a parameterization for the given curve, including bounds on the parameter t. Part a) The ellipse in R whose major axis has endpoints, ) and 6, )
More informationSAMPLING SPECIAL DISTRIBUTIONS
SAMPLING SPECIAL DISTRIBUTIONS M. Ragheb 4/9/4 INTRODUCTION Monte Carlo simulations require the sampling of various distributions at different steps in the simulation process. There exist a multitude of
More informationTopic 7. Electric flux Gauss s Law Divergence of E Application of Gauss Law Curl of E
Topic 7 Electric flux Gauss s Law Divergence of E Application of Gauss Law Curl of E urface enclosing an electric dipole. urface enclosing charges 2q and q. Electric flux Flux density : The number of field
More informationElastic Scattering. R = m 1r 1 + m 2 r 2 m 1 + m 2. is the center of mass which is known to move with a constant velocity (see previous lectures):
Elastic Scattering In this section we will consider a problem of scattering of two particles obeying Newtonian mechanics. The problem of scattering can be viewed as a truncated version of dynamic problem
More informationNotes 19 Gradient and Laplacian
ECE 3318 Applied Electricity and Magnetism Spring 218 Prof. David R. Jackson Dept. of ECE Notes 19 Gradient and Laplacian 1 Gradient Φ ( x, y, z) =scalar function Φ Φ Φ grad Φ xˆ + yˆ + zˆ x y z We can
More informationCalculus III 2004 Summer Practice Final 8/3/2004
.. Calculus III 4 ummer Practice Final 8/3/4. Compute the following limits if they exist: (a) lim (x,y) (,) e xy x+. cos x (b) lim x. (x,y) (,) x 4 +y 4 (a) ince lim (x,y) (,) exy and lim x + 6 in a (x,y)
More informationCreated by T. Madas SURFACE INTEGRALS. Created by T. Madas
SURFACE INTEGRALS Question 1 Find the area of the plane with equation x + 3y + 6z = 60, 0 x 4, 0 y 6. 8 Question A surface has Cartesian equation y z x + + = 1. 4 5 Determine the area of the surface which
More informationApplications of Gauss Law
Applications of Gauss Law The Gauss Law of electrostatics relates the net electric field flux through a complete surface S of some volume V to the net electric charge inside that volume, Φ E [S] def =
More information1. (30 points) In the x-y plane, find and classify all local maxima, local minima, and saddle points of the function. f(x, y) = 3y 2 2y 3 3x 2 + 6xy.
APPM 35 FINAL EXAM FALL 13 INSTUTIONS: Electronic devices, books, and crib sheets are not permitted. Write your name and your instructor s name on the front of your bluebook. Work all problems. Show your
More informationVolume in n Dimensions
Volume in n Dimensions MA 305 Kurt Bryan Introduction You ve seen that if we have two vectors v and w in two dimensions then the area spanned by these vectors can be computed as v w = v 1 w 2 v 2 w 1 (where
More informationChapter 5 Trigonometric Functions of Angles
Chapter 5 Trigonometric Functions of Angles Section 3 Points on Circles Using Sine and Cosine Signs Signs I Signs (+, +) I Signs II (+, +) I Signs II (, +) (+, +) I Signs II (, +) (+, +) I III Signs II
More information1. (3) Write Gauss Law in differential form. Explain the physical meaning.
Electrodynamics I Midterm Exam - Part A - Closed Book KSU 204/0/23 Name Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try to tell about the physics involved,
More information1 CHAPTER 5 ABSORPTION, SCATTERING, EXTINCTION AND THE EQUATION OF TRANSFER
1 CHAPTER 5 ABSORPTION, SCATTERING, EXTINCTION AND THE EQUATION OF TRANSFER 5.1 Introduction As radiation struggles to make its way upwards through a stellar atmosphere, it may be weakened by absorption
More informationWeek 7: Integration: Special Coordinates
Week 7: Integration: Special Coordinates Introduction Many problems naturally involve symmetry. One should exploit it where possible and this often means using coordinate systems other than Cartesian coordinates.
More information