Space-time kinetics. Alain Hébert. Institut de génie nucléaire École Polytechnique de Montréal.
|
|
- Erica Fields
- 5 years ago
- Views:
Transcription
1 Space-time kinetics Alain Hébert Institut de génie nucléaire École Polytechnique de Montréal ENE613: Week 1 Space-time kinetics 1/24
2 Content (week 1) 1 Space-time kinetics equations Point kinetics equations Equivalent point kinetics parameters Point kinetics equations (classical form) Inhour equation with precursors Calculation of an exponential matrix The theta finite-difference scheme The theta finite-difference scheme Theta scheme with point kinetics Stability condition Space-time kinetics with Θ = 1 Space-time kinetics with Θ = ENE613: Week 1 Space-time kinetics 2/24
3 Space-time kinetics equations 1 1 V n,g = t φ g(r, t) D g (r) φ g (r, t) + Σ rg (r) φ g (r, t) + X l GX h=1 h g Σ g h (r) φ h (r, t) + χ pr g (r) (1 β) G X χ del l,g (r) λ l c l (r, t) + q g (r, t) h=1 νσ fh (r) φ h (r, t) c l (r, t) t X G = β l νσ fh (r) φ h (r, t) λ l c l (r, t) ; l = 1, N d h=1 β = X l β l. ENE613: Week 1 Space-time kinetics 3/24
4 Point kinetics equations 1 1 V n (E) φ(e, t) + Σ(E) φ(e, t) = t + χ pr (E)(1 β) + X l Z Z χ del l (E) λ l c l (t) + q(e, t) de Σ s, (E E ) φ(e, t) de νσ f (E ) φ(e, t) together with the set of N d precursor equations: c l (t) t = β l Z de νσ f (E) φ(e, t) λ l c l (t) ; l = 1, N d. ENE613: Week 1 Space-time kinetics 4/24
5 Equivalent point kinetics parameters 1 The time-dependent reactivity of the reactor is ρ(t) = 1 R de φ (E)»Σ(E) φ(e, t) R de φ (E) χss (E) R R de Σ s, (E E ) φ(e, t) de νσ f (E ) φ(e, t), the steady-state fission spectrum is χ ss (E) = (1 β) χ pr (E) + X l β l χ del l (E), ENE613: Week 1 Space-time kinetics 5/24
6 Equivalent point kinetics parameters 2 the mean neutron generation time of the reactor is Λ(t) = R R de V 1 n (E) φ (E) φ(e, t) R de φ (E) χss (E) de νσ f (E ) φ(e, t), the average delayed fractions are β l (t) = β l R R de φ (E) χdel l (E) de φ (E) χss (E) ENE613: Week 1 Space-time kinetics 6/24
7 Equivalent point kinetics parameters 3 and β(t) = X l β l (t) = 1 (1 β) R R de φ (E) χpr (E) de φ (E) χss (E) and where the neutron population, precursor concentration and external source are n(t) = Z de 1 V n (E) φ (E) φ(e, t), c l (t) = c l (t) Z de φ (E) χ del l (E) and q(t) = Z de φ (E) q(e, t). ENE613: Week 1 Space-time kinetics 7/24
8 Point kinetics equations (classical form) 1 d dt n(t) = ρ(t) β(t) Λ(t) n(t) + X l λ l c l (t) + q(t) d dt c l(t) = β l (t) Λ(t) n(t) λ l c l (t) ; l = 1, N d ENE613: Week 1 Space-time kinetics 8/24
9 Point kinetics equations (classical form) 2 Point kinetics equations are written in matrix form as x() = x d x(t) = A x(t) + S dt x(t) = [ n(t) c 1 (t) c 2 (t)... c Nd (t)], S = [ q... ] A = ρ β Λ β 1 Λ β 2 Λ. β Nd Λ λ 1 λ 2... λ Nd λ 1... λ λ Nd ENE613: Week 1 Space-time kinetics 9/24
10 Point kinetics equations (classical form) 3 Steady-state equation: A x = c l () c l, = β l Λ λ l n ; l = 1, N d. Each side of the point kinetics equation is multiplied by the integrating factor e At and integrated in time. We obtain the analytical solution of the point-kinetics equations as x(t) = e At x + (e At I) A 1 S where we made use of the identity Z t» d dt e At dt x(t ) A x(t ) = e At x(t) x(). ENE613: Week 1 Space-time kinetics 1/24
11 Calculation of an exponential matrix 1 A straightforward technique consists to obtain the N d + 1 eigenvalues and eigenvectors of linear system (A ω i I) p i = where 1 i N d + 1 The eigenvalues are the root of the characteristic polynomial det (A ω I) = which can be obtained as the roots of the inhour equation : ρ β Λ ω + X l β l λ l λ l + ω =. ENE613: Week 1 Space-time kinetics 11/24
12 Inhour equation with precursors 1 ω = ρ β Λ ρ = β ω = ω = λ 6 ω = β Λ ρ = ENE613: Week 1 Space-time kinetics 12/24
13 Inhour equation with precursors 2 We reproduce here a simple Matlab script for solving the inhour equation ρ β Λ ω + X l β l λ l λ l + ω = function y = inhour(lambda, L, beta, rho) % Find the roots of the inhour equation % function y = inhour(lambda, L, beta, rho) % (c) 27 Alain Hebert, Ecole Polytechnique de Montreal [nn, dd]=residue(lambda.*beta,-lambda,[]) ; y=sort(roots(conv([-l, rho-sum(beta)],dd)+[ nn])) ; The script parameters are defined as lambda= array of dimension N d containing the decay constants λ l for the delayed neutrons L= neutron generation time Λ beta= array of dimension N d containing the delayed fractions β l rho= reactivity ρ y= array of dimension N d + 1 containing the roots ω i of the inhour equation. ENE613: Week 1 Space-time kinetics 13/24
14 Calculation of an exponential matrix 1 The analytical solution also requires the knowledge of the matrix P whose columns are the eigenvectors of square matrix A. This matrix is equal to P = {p i,j ; i = 1, N d + 1 and j = 1, N d + 1} where p 1,j = 1 and p i+1,j = β i Λ (λ i + ω j ). ENE613: Week 1 Space-time kinetics 14/24
15 Calculation of an exponential matrix 1 Similarly, we can show that P 1 = {y i,j ; i = 1, N d + 1 and j = 1, N d + 1} where y i,1 = " 1 + X l β l λ l Λ (λ l + ω i ) 2 # 1 and y i,j+1 = y i,1 λ j λ j + ω i. Theorem If P is the matrix whose columns are the eigenvectors of square matrix A, and if ω i are the corresponding eigenvalues, then e At = P diag(e ω it ) P 1. ENE613: Week 1 Space-time kinetics 15/24
16 The theta finite-difference scheme 1 A finite difference relation is used for time discretization of neutron flux and precursor equations, written as (1 Θ) t φ g(r, t) + Θ tn 1 t φ g(r, t) = φ g(r, t n ) φ g (r, t n 1 ) tn t n (1 Θ) t c l(r, t) + Θ tn 1 t c l(r, t) = c l(r, t n ) c l (r, t n 1 ) tn t n where l = 1, N d. Setting Θ =, Θ = 1 and Θ = 1/2 yield the explicit scheme, the fully implicit scheme and the Crank-Nicholson scheme, respectively. ENE613: Week 1 Space-time kinetics 16/24
17 Theta scheme with point kinetics 1 x (n+1) x (n) t n = Θ A x (n+1) + (1 Θ) A x (n) + S where x (n) x(t n ). The theta scheme can be applied to the point-kinetics equations and used to study its stability characteristics. It can also be applied to the space-time kinetics equations for solving full-core problems. Above equation with S = may be written as x (n+1) = R x (n) where I is the identity matrix and where the iterative matrix is defined as R =» I 1» I Θ A + (1 Θ) A t n t n. ENE613: Week 1 Space-time kinetics 17/24
18 Theta scheme with point kinetics 2 The eigenvalue problem associated to the above iterative scheme is written γ i p i R p i = where γ i is the i th eigenvalue of matrix R and p i is the corresponding eigenvector. We can easily show that these eigenvectors are identical to those of the inhour equation and that the corresponding eigenvalues are related using γ i = 1 + (1 Θ) ω i t n 1 Θ ω i t n. The eigenvectors of matrix R are linearly independent, so that the initial condition x can be expressed as a linear combination: x x () = X i c i p i. ENE613: Week 1 Space-time kinetics 18/24
19 Stability condition 1 After a progression of n time-steps, we have x (n) = R x (n 1) = R X i c i γ n 1 i p i = X i c i γ n i p i. We apply a negative reactivity ρ < to a point kinetics problem and verify that the numerical solution given by previous equation vanishes as t increases. The stability will be guaranteed if γ i < 1 for all values of i. If ρ <, we observe from inhour figure that all eigenvalues ω i are negative. γ i equation leads to condition γ i < 1 if (1 + Θ ω i t n ) < 1 (1 Θ) ω i t n < 1 + Θ ω i t n which is satisfied if 1 2Θ 2 ω i t n < 1 ; i. ENE613: Week 1 Space-time kinetics 19/24
20 Stability condition 2 We conclude that the explicit scheme (Θ = ) is stable if the time-step size is set below a threshold value written as 2 t n < max i ( ω i ). the fully implicit and Crank-Nicholson schemes are unconditionally stable. The practical choice of Θ for production calculations is selected in interval 1/2 Θ 1. ENE613: Week 1 Space-time kinetics 2/24
21 Space-time kinetics with Θ = V n,g t n φ g (r, t n ) D g (r) φ g (r, t n ) + Σ rg (r) φ g (r, t n ) = S imp g (r, t n ) + + ( χ ss g (r) X l GX h=1 h g χ del l,g (r) Σ g h (r) φ h (r, t n ) ) β l X G νσ fh (r) φ h (r, t n ) 1 + λ l t n h=1 c l (r, t n ) = λ l t n " G # X c l (r, t n 1 ) + β l t n νσ fh (r) φ h (r, t n ) h=1 ENE613: Week 1 Space-time kinetics 21/24
22 Space-time kinetics with Θ = 1 2 with the fixed source defined as S imp g (r, t n ) = 1 V n,g t n φ g (r, t n 1 ) + X l χ del l,g (r) λ l 1 + λ l t n c l (r, t n 1 ) + q g (r, t n ) ENE613: Week 1 Space-time kinetics 22/24
23 Space-time kinetics with Θ = 1 1 V n,g t n φ g (r, t n ) = S exp g (r, t n ) X G c l (r, t n ) = (1 λ l t n ) c l (r, t n 1 ) + β l t n νσ fh (r) φ h (r, t n 1 ) h=1 ENE613: Week 1 Space-time kinetics 23/24
24 Space-time kinetics with Θ = 2 with the fixed source defined as S exp g (r, t n ) = 1 V n,g t n φ g (r, t n 1 ) + D g (r) φ g (r, t n 1 ) Σ rg (r) φ g (r, t n 1 ) + + " + X l χ ss g (r) X l GX h=1 h g β l χ del l,g (r) # G X Σ g h (r) φ h (r, t n 1 ) h=1 νσ fh (r) φ h (r, t n 1 ) χ del l,g (r) λ l c l (r, t n 1 ) + q g (r, t n 1 ). ENE613: Week 1 Space-time kinetics 24/24
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 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 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 informationThe 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 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 information"Control Rod Calibration"
TECHNICAL UNIVERSITY DRESDEN Institute of Power Engineering Training Reactor Reactor Training Course Experiment "Control Rod Calibration" Instruction for Experiment Control Rod Calibration Content: 1...
More informationCOMPARISON OF DIRECT AND QUASI-STATIC METHODS FOR NEUTRON KINETIC CALCULATIONS WITH THE EDF R&D COCAGNE CODE
PHYSOR 2012 Advances in Reactor Physics Linking Research, Industry, and Education, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2012) COMPARISON OF DIRECT AND QUASI-STATIC METHODS FOR NEUTRON
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationSection 8.2 : Homogeneous Linear Systems
Section 8.2 : Homogeneous Linear Systems Review: Eigenvalues and Eigenvectors Let A be an n n matrix with constant real components a ij. An eigenvector of A is a nonzero n 1 column vector v such that Av
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 informationCayley-Hamilton Theorem
Cayley-Hamilton Theorem Massoud Malek In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n Let A be an n n matrix Although det (λ I n A
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 informationDelayed neutrons in nuclear fission chain reaction
Delayed neutrons in nuclear fission chain reaction 1 Critical state Temporal flow Loss by leakage Loss by Absorption If the number of neutrons (the number of fission reactions) is practically constant
More informationMatrix Solutions to Linear Systems of ODEs
Matrix Solutions to Linear Systems of ODEs James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 3, 216 Outline 1 Symmetric Systems of
More informationNumerical simulation of non-steady state neutron kinetics of the TRIGA Mark II reactor Vienna
arxiv:137.767v1 [physics.ins-det] 29 Jul 213 Numerical simulation of non-steady state neutron kinetics of the TRIGA Mark II reactor Vienna Julia Riede, Helmuth Boeck TU Wien, Atominstitut, A-12 Wien, Stadionallee
More informationCalculus C (ordinary differential equations)
Calculus C (ordinary differential equations) Lesson 9: Matrix exponential of a symmetric matrix Coefficient matrices with a full set of eigenvectors Solving linear ODE s by power series Solutions to linear
More informationReactor Operation Without Feedback Effects
22.05 Reactor Physics - Part Twenty-Six Reactor Operation Without Feedback Effects 1. Reference Material: See pp. 363-368 of the article, Light Water Reactor Control Systems, in Wiley Encyclopedia of Electrical
More informationSTABILITY FOR PARABOLIC SOLVERS
Review STABILITY FOR PARABOLIC SOLVERS School of Mathematics Semester 1 2008 OUTLINE Review 1 REVIEW 2 STABILITY: EXPLICIT METHOD Explicit Method as a Matrix Equation Growing Errors Stability Constraint
More informationJordan normal form notes (version date: 11/21/07)
Jordan normal form notes (version date: /2/7) If A has an eigenbasis {u,, u n }, ie a basis made up of eigenvectors, so that Au j = λ j u j, then A is diagonal with respect to that basis To see this, let
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationA = 3 1. We conclude that the algebraic multiplicity of the eigenvalues are both one, that is,
65 Diagonalizable Matrices It is useful to introduce few more concepts, that are common in the literature Definition 65 The characteristic polynomial of an n n matrix A is the function p(λ) det(a λi) Example
More informationOperational Reactor Safety
Operational Reactor Safety 22.091/22.903 Professor Andrew C. Kadak Professor of the Practice Lecture 3 Reactor Kinetics and Control Page 1 Topics to Be Covered Time Dependent Diffusion Equation Prompt
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
More informationELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky
ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky Equilibrium points and linearization Eigenvalue decomposition and modal form State transition matrix and matrix exponential Stability ELEC 3035 (Part
More informationMath 216 Final Exam 24 April, 2017
Math 216 Final Exam 24 April, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 6: Numerical solution of the heat equation with FD method: method of lines, Euler
More informationStudy guide - Math 220
Study guide - Math 220 November 28, 2012 1 Exam I 1.1 Linear Equations An equation is linear, if in the form y + p(t)y = q(t). Introducing the integrating factor µ(t) = e p(t)dt the solutions is then in
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More information6 Linear Equation. 6.1 Equation with constant coefficients
6 Linear Equation 6.1 Equation with constant coefficients Consider the equation ẋ = Ax, x R n. This equating has n independent solutions. If the eigenvalues are distinct then the solutions are c k e λ
More informationMATH 5640: Functions of Diagonalizable Matrices
MATH 5640: Functions of Diagonalizable Matrices Hung Phan, UMass Lowell November 27, 208 Spectral theorem for diagonalizable matrices Definition Let V = X Y Every v V is uniquely decomposed as u = x +
More informationDifferential equations
Differential equations Math 7 Spring Practice problems for April Exam Problem Use the method of elimination to find the x-component of the general solution of x y = 6x 9x + y = x 6y 9y Soln: The system
More informationReview of Linear Algebra
Review of Linear Algebra Definitions An m n (read "m by n") matrix, is a rectangular array of entries, where m is the number of rows and n the number of columns. 2 Definitions (Con t) A is square if m=
More informationMathematics 22: Lecture 19
Mathematics 22: Lecture 19 Legendre s Equation Dan Sloughter Furman University February 5, 2008 Dan Sloughter (Furman University) Mathematics 22: Lecture 19 February 5, 2008 1 / 11 Example: Legendre s
More informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
More informationXenon Effects. B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec.
enon Effects B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec. 2015 September 1 Contents We study the importance of e-135 in the operation of
More informationDiagonalization. P. Danziger. u B = A 1. B u S.
7., 8., 8.2 Diagonalization P. Danziger Change of Basis Given a basis of R n, B {v,..., v n }, we have seen that the matrix whose columns consist of these vectors can be thought of as a change of basis
More informationbut mostly as the result of the beta decay of its precursor 135 I (which has a half-life of hours).
8. Effects of 135Xe The xenon isotope 135 Xe plays an important role in any power reactor. It has a very large absorption cross section for thermal neutrons and represents therefore a considerable load
More informationSolutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 28, 2011
Solutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 8, Solutions to Problems 5 are omitted since they involve topics no longer covered on the Comprehensive
More informationLinear Algebra Exercises
9. 8.03 Linear Algebra Exercises 9A. Matrix Multiplication, Rank, Echelon Form 9A-. Which of the following matrices is in row-echelon form? 2 0 0 5 0 (i) (ii) (iii) (iv) 0 0 0 (v) [ 0 ] 0 0 0 0 0 0 0 9A-2.
More informationComputational Methods CMSC/AMSC/MAPL 460. Eigenvalues and Eigenvectors. Ramani Duraiswami, Dept. of Computer Science
Computational Methods CMSC/AMSC/MAPL 460 Eigenvalues and Eigenvectors Ramani Duraiswami, Dept. of Computer Science Eigen Values of a Matrix Recap: A N N matrix A has an eigenvector x (non-zero) with corresponding
More informationChapter 9 Part II Maintainability
Chapter 9 Part II Maintainability 9.4 System Repair Time 9.5 Reliability Under Preventive Maintenance 9.6 State-Dependent Systems with Repair C. Ebeling, Intro to Reliability & Maintainability Chapter
More informationReactor Physics: Basic Definitions and Perspectives. Table of Contents
Reactor Physics - Basic Definitions and Perspectives Reactor Physics: Basic Definitions and Perspectives prepared by Wm. J. Garland, Professor, Department of Engineering Physics, McMaster University, Hamilton,
More informationMath 54. Selected Solutions for Week 10
Math 54. Selected Solutions for Week 10 Section 4.1 (Page 399) 9. Find a synchronous solution of the form A cos Ωt+B sin Ωt to the given forced oscillator equation using the method of Example 4 to solve
More informationNumerical Techniques for the Neutron Diffusion Equations in the Nuclear Reactors
Adv. Studies Theor. Phys., Vol. 6, 01, no. 14, 649-664 Numerical Techniques for the Neutron Diffusion Equations in the Nuclear Reactors Abdallah A. Nahla Department of Mathematics, Faculty of Science Taif
More informationMCE 366 System Dynamics, Spring Problem Set 2. Solutions to Set 2
MCE 366 System Dynamics, Spring 2012 Problem Set 2 Reading: Chapter 2, Sections 2.3 and 2.4, Chapter 3, Sections 3.1 and 3.2 Problems: 2.22, 2.24, 2.26, 2.31, 3.4(a, b, d), 3.5 Solutions to Set 2 2.22
More informationSolution: In standard form (i.e. y + P (t)y = Q(t)) we have y t y = cos(t)
Math 380 Practice Final Solutions This is longer than the actual exam, which will be 8 to 0 questions (some might be multiple choice). You are allowed up to two sheets of notes (both sides) and a calculator,
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationGoals: Introduction to Systems of Differential Equations Solving systems with distinct real eigenvalues and eigenvectors
Week #10 : Systems of DEs Goals: Introduction to Systems of Differential Equations Solving systems with distinct real eigenvalues and eigenvectors 1 Systems of Differential Equations - Introdution - 1
More informationMATH 320 INHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS
MATH 2 INHOMOGENEOUS LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS W To find a particular solution for a linear inhomogeneous system of differential equations x Ax = ft) or of a mechanical system with external
More informationMTH739U/P: Topics in Scientific Computing Autumn 2016 Week 6
MTH739U/P: Topics in Scientific Computing Autumn 16 Week 6 4.5 Generic algorithms for non-uniform variates We have seen that sampling from a uniform distribution in [, 1] is a relatively straightforward
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationApplicability of Point Reactor Kinetic Equations for Subcriticality Monitoring
Proceedings of Korean Nuclear Society Spring Meeting Cheji, Korea, May 21 Applicability of Point Reactor Kinetic Equations for Subcriticality Monitoring Yoichiro Shimazu, Masashi Tsuji Hokkaido University
More informationWe shall finally briefly discuss the generalization of the solution methods to a system of n first order differential equations.
George Alogoskoufis, Dynamic Macroeconomic Theory, 2015 Mathematical Annex 1 Ordinary Differential Equations In this mathematical annex, we define and analyze the solution of first and second order linear
More informationMa 227 Review for Systems of DEs
Ma 7 Review for Systems of DEs Matrices Basic Properties Addition and subtraction: Let A a ij mn and B b ij mn.then A B a ij b ij mn 3 A 6 B 6 4 7 6 A B 6 4 3 7 6 6 7 3 Scaler Multiplication: Let k be
More informationAM 034 Applied Mathematics - II Brown University Spring 2018 Solutions to Homework, Set 3 Due February 28
AM 034 Applied Mathematics - II Brown University Spring 18 Solutions to Homework, Set 3 Due February 28 31 ( points) Find all square roots of the following matrices: [ ] 2 3 9 31 18 (a), (b) 3 2 9 22 27
More informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More informationGeneral Exam Part II, Fall 1998 Quantum Mechanics Solutions
General Exam Part II, Fall 1998 Quantum Mechanics Solutions Leo C. Stein Problem 1 Consider a particle of charge q and mass m confined to the x-y plane and subject to a harmonic oscillator potential V
More informationEntropy-dissipation methods I: Fokker-Planck equations
1 Entropy-dissipation methods I: Fokker-Planck equations Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Introduction Boltzmann equation Fokker-Planck equations Degenerate parabolic
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 1. Linear systems Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Linear systems 1.1 Differential Equations 1.2 Linear flows 1.3 Linear maps
More informationNEW FUEL IN MARIA RESEARCH REACTOR, PROVIDING BETTER CONDITIONS FOR IRRADIATION IN THE FAST NEUTRON SPECTRUM.
NEW FUEL IN MARIA RESEARCH REACTOR, PROVIDING BETTER CONDITIONS FOR IRRADIATION IN THE FAST NEUTRON SPECTRUM. M. LIPKA National Centre for Nuclear Research Andrzeja Sołtana 7, 05-400 Otwock-Świerk, Poland
More informationNORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009
NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and
More information14 : Theory of Variational Inference: Inner and Outer Approximation
10-708: Probabilistic Graphical Models 10-708, Spring 2014 14 : Theory of Variational Inference: Inner and Outer Approximation Lecturer: Eric P. Xing Scribes: Yu-Hsin Kuo, Amos Ng 1 Introduction Last lecture
More informationApplications of Diagonalization
Applications of Diagonalization Hsiu-Hau Lin hsiuhau@phys.nthu.edu.tw Apr 2, 200 The notes cover applications of matrix diagonalization Boas 3.2. Quadratic curves Consider the quadratic curve, 5x 2 4xy
More informationAPPM 2360: Midterm exam 3 April 19, 2017
APPM 36: Midterm exam 3 April 19, 17 On the front of your Bluebook write: (1) your name, () your instructor s name, (3) your lecture section number and (4) a grading table. Text books, class notes, cell
More informationOctober 4, 2017 EIGENVALUES AND EIGENVECTORS. APPLICATIONS
October 4, 207 EIGENVALUES AND EIGENVECTORS. APPLICATIONS RODICA D. COSTIN Contents 4. Eigenvalues and Eigenvectors 3 4.. Motivation 3 4.2. Diagonal matrices 3 4.3. Example: solving linear differential
More informationPh.D. Katarína Bellová Page 1 Mathematics 2 (10-PHY-BIPMA2) EXAM - Solutions, 20 July 2017, 10:00 12:00 All answers to be justified.
PhD Katarína Bellová Page 1 Mathematics 2 (10-PHY-BIPMA2 EXAM - Solutions, 20 July 2017, 10:00 12:00 All answers to be justified Problem 1 [ points]: For which parameters λ R does the following system
More informationADAPTIVE EXTREMUM SEEKING CONTROL OF CONTINUOUS STIRRED TANK BIOREACTORS 1
ADAPTIVE EXTREMUM SEEKING CONTROL OF CONTINUOUS STIRRED TANK BIOREACTORS M. Guay, D. Dochain M. Perrier Department of Chemical Engineering, Queen s University, Kingston, Ontario, Canada K7L 3N6 CESAME,
More informationSolutions to Math 53 Math 53 Practice Final
Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points
More informationSecond Order Systems
Second Order Systems independent energy storage elements => Resonance: inertance & capacitance trade energy, kinetic to potential Example: Automobile Suspension x z vertical motions suspension spring shock
More informationPeriodic Linear Systems
Periodic Linear Systems Lecture 17 Math 634 10/8/99 We now consider ẋ = A(t)x (1) when A is a continuous periodic n n matrix function of t; i.e., when there is a constant T>0 such that A(t + T )=A(t) for
More informationSeptember 26, 2017 EIGENVALUES AND EIGENVECTORS. APPLICATIONS
September 26, 207 EIGENVALUES AND EIGENVECTORS. APPLICATIONS RODICA D. COSTIN Contents 4. Eigenvalues and Eigenvectors 3 4.. Motivation 3 4.2. Diagonal matrices 3 4.3. Example: solving linear differential
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationChemical Engineering 412
Chemical Engineering 412 Introductory Nuclear Engineering Lecture 18 Nuclear Reactor Theory IV Reactivity Insertions 1 Spiritual Thought 2 Mosiah 2:33 33 For behold, there is a wo pronounced upon him who
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationDiagonalization of Matrices
LECTURE 4 Diagonalization of Matrices Recall that a diagonal matrix is a square n n matrix with non-zero entries only along the diagonal from the upper left to the lower right (the main diagonal) Diagonal
More informationFind the Fourier series of the odd-periodic extension of the function f (x) = 1 for x ( 1, 0). Solution: The Fourier series is.
Review for Final Exam. Monday /09, :45-:45pm in CC-403. Exam is cumulative, -4 problems. 5 grading attempts per problem. Problems similar to homeworks. Integration and LT tables provided. No notes, no
More informationj=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.
Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 0
MATH 56A: STOCHASTIC PROCESSES CHAPTER 0 0. Chapter 0 I reviewed basic properties of linear differential equations in one variable. I still need to do the theory for several variables. 0.1. linear differential
More informationChapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors
Chapter 7 Canonical Forms 7.1 Eigenvalues and Eigenvectors Definition 7.1.1. Let V be a vector space over the field F and let T be a linear operator on V. An eigenvalue of T is a scalar λ F such that there
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationDynamic interpretation of eigenvectors
EE263 Autumn 2015 S. Boyd and S. Lall Dynamic interpretation of eigenvectors invariant sets complex eigenvectors & invariant planes left eigenvectors modal form discrete-time stability 1 Dynamic interpretation
More informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 State Estimator In previous section, we have discussed the state feedback, based on the assumption that all state variables are
More informationStability of generalized Runge Kutta methods for stiff kinetics coupled differential equations
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J Phys A: Math Gen 39 (2006) 859 876 doi:0088/0305-4470/39/8/006 Stability of generalized Runge Kutta methods for stiff kinetics
More informationRadioactivity. Radioactivity
The Law of Radioactive Decay. 72 The law of radioactive decay. It turns out that the probability per unit time for any radioactive nucleus to decay is a constant, called the decay constant, lambda, ".
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationPower Changes in a Critical Reactor. The Critical Reactor
Chapter 8 Power Changes in a Critical Reactor n For very small reactivity increases n For small reactivity increases n For large reactivity increases/decreases The Critical Reactor < k = hfpel f L t =
More informationChapter III. Stability of Linear Systems
1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,
More information8 Ecosystem stability
8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single
More informationLocalization I: General considerations, one-parameter scaling
PHYS598PTD A.J.Leggett 2013 Lecture 4 Localization I: General considerations 1 Localization I: General considerations, one-parameter scaling Traditionally, two mechanisms for localization of electron states
More informationDynamical Systems & Scientic Computing: Homework Assignments
Fakultäten für Informatik & Mathematik Technische Universität München Dr. Tobias Neckel Summer Term Dr. Florian Rupp Exercise Sheet 3 Dynamical Systems & Scientic Computing: Homework Assignments 3. [ ]
More informationLecture 6. Eigen-analysis
Lecture 6 Eigen-analysis University of British Columbia, Vancouver Yue-Xian Li March 7 6 Definition of eigenvectors and eigenvalues Def: Any n n matrix A defines a LT, A : R n R n A vector v and a scalar
More informationTranscendence in Positive Characteristic
Transcendence in Positive Characteristic Galois Group Examples and Applications W. Dale Brownawell Matthew Papanikolas Penn State University Texas A&M University Arizona Winter School 2008 March 18, 2008
More informationDesigning Information Devices and Systems I Discussion 4B
Last Updated: 29-2-2 9:56 EECS 6A Spring 29 Designing Information Devices and Systems I Discussion 4B Reference Definitions: Matrices and Linear (In)Dependence We ve seen that the following statements
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More informationVectors, matrices, eigenvalues and eigenvectors
Vectors, matrices, eigenvalues and eigenvectors 1 ( ) ( ) ( ) Scaling a vector: 0.5V 2 0.5 2 1 = 0.5 = = 1 0.5 1 0.5 ( ) ( ) ( ) ( ) Adding two vectors: V + W 2 1 2 + 1 3 = + = = 1 3 1 + 3 4 ( ) ( ) a
More informationProcess Modelling. Table of Contents
Process Modelling 1 Process Modelling prepared by Wm. J. Garland, Professor, Department of Engineering Physics, McMaster University, Hamilton, Ontario, Canada More about this document Summary: The general
More informationLecture Notes for Math 251: ODE and PDE. Lecture 27: 7.8 Repeated Eigenvalues
Lecture Notes for Math 25: ODE and PDE. Lecture 27: 7.8 Repeated Eigenvalues Shawn D. Ryan Spring 22 Repeated Eigenvalues Last Time: We studied phase portraits and systems of differential equations with
More information22 APPENDIX 1: MATH FACTS
22 APPENDIX : MATH FACTS 22. Vectors 22.. Definition A vector has a dual definition: It is a segment of a a line with direction, or it consists of its projection on a reference system xyz, usually orthogonal
More information18.06 Problem Set 8 - Solutions Due Wednesday, 14 November 2007 at 4 pm in
806 Problem Set 8 - Solutions Due Wednesday, 4 November 2007 at 4 pm in 2-06 08 03 Problem : 205+5+5+5 Consider the matrix A 02 07 a Check that A is a positive Markov matrix, and find its steady state
More information