Implicit schemes of Lax-Friedrichs type for systems with source terms
|
|
- Irene Osborne
- 6 years ago
- Views:
Transcription
1 Implicit schemes of Lax-Friedrichs type for systems with source terms A. J. Forestier 1, P. González-Rodelas 2 1 DEN/DSNI/Réacteurs C.E.A. Saclay (France 2 Department of Applied Mathematics, Universidad de Granada (Spain Castro Urdiales, 7-11 Sept., 2009 (Spain
2 Index 1 Introduction 2 Implicit Lax-Friedrichs and local Lax-Friedrichs schemes in the scalar case 3 Implicit Lax-Friedrichs and local Lax-Friedrichs schemes for systems
3 Implicit schemes Are taking more and more importance in the CFD and other areas. Better stability properties. Allows us to employ larger CFL numbers. Reaching the specified instant time with fewer iterations. But with possibly loosing some transient information. So specially appropriate for stationary solutions. Leads to solve large non linear systems of equations. It is convenient to linearize them. Need of efficient linear algebra routines (specially for tri-diagonal and block-diagonal matrices.
4 Index 1 Introduction 2 Implicit Lax-Friedrichs and local Lax-Friedrichs schemes in the scalar case 3 Implicit Lax-Friedrichs and local Lax-Friedrichs schemes for systems
5 Formulation of the Lax-Friedrichs scheme As a variant of the unconditionally unstable FTCS scheme in the approximation of the time derivative, done now as t u ( un 1 +un +1 2 t n, x t = 2un+1 u n 1 un +1 2 t we obtain the so called (explicit Lax-Friedrichs scheme where σ t x = u n n σ ( f f n and independently of n N {0} we have LF LF f f == f +1 f 1 2 u +1 2u + u 1 2σ
6 Lax-Friedrichs numerical flux This results in the following classical L-F numerical flux f LF (u, v = 1 (f (u + f (v 1σ 2 (v u (1 where f 1 2 with u n f LF ( u 1, u, f+ 1 2 f± 1 2 = 1 2 u ( t n, x and f n ( f ±1 + f 1 σ f ( u n f LF ( u, u +1. in such a way that it has good stability properties. ( u±1 u (2
7 Rusanov or Local Lax-Friedrichs flux Is ust a variation of the above L-F one, replacing the constant viscosity coefficient ( 1 σ = x t by a new one, calculated locally for each one of the Riemann problems ζ = max f (q q [u,v] (where q [u, v] indicates that q can be lengthwise the segment [u, v] := {q = t u + (1 t v / t [0, 1]} For example, in the case of a convex or concave flux ζ = max { f (u, f (v }.
8 Formulation of the so-called θ-schemes In a similar way, depending on the specific numerical flux f n+1 taken, writing f := f ( 1 1, un+1, we can consider the 2 ( explicit scheme = u n σ f n f n ( implicit scheme n+1 n+1 + σ f f = u n 2 2 or even taking any convex combination of these two, the θ-schemes with θ [0, 1]: ( n+1 n+1 + θ σ f f = u n n (1 θ σ ( f f n (3
9 Formulation and study of the implicit schemes They would lead in general to systems of nonlinear equations that would have to be solved for each time step. Nevertheless, we can develop the expressions in the implicit parts, if we can easily find ṽ n+1 v n+1 in such a way that ( f ṽ n+1 when f ( ṽ n+1 ( f f v n+1 = +1 un+1 1 ( +1 un+1 1 ( +1 ( f 1 +1 un+1 1 f n+1 +1 f n un+1 1 ; but in any case we could also write ( ( f v n+1 +1 un+1 1 = f n+1 n+1 +1 f 1
10 Implicit θ-schemes In the case of considering periodic boundary conditions, the coefficients matrix appearing in each iteration step will have the following form β n 1 γ n α n 1 α n 2 β n 2 γ n α n 1 β n 1 γ n α n β n γ n 0 0 α n +1 β n +1 γ n α n J 1 β n J 1 γ n J 1 γ n J α n J β n J J 1 J = ( δ n 1 δ n 2 δ n 1 δ n δ n +1 δ n J 1 δ n J
11 Implicit θ-schemes So we have, δ n 0 = αn 0 un 1 + βn 0 un 0 + γn 0 un 1 = α n J un J 1 + βn J un J + γn J un J+1 = δn J whereas if the boundary conditions would be of Dirichlet type, also the first and last equations above would not have more sense, and we would have to restrict ourselves to the indices = 1,..., J 1. But both in the case of periodic or Dirichlet boundary conditions, we can apply the same type of reasoning.
12 FTCS and BTCS θ-schemes In the case of a particular θ-scheme, with the centered flux: + σ ( ( 2 θ f v n+1 +1 un+1 1 = u n (1 θ σ 2 f ( v n ( u+1 1 n un δ n (θ the linear system would be clearly tridiagonal α n 1 +βn +γ n +1 = δn (θ, n = 0, 1,... ; = 0, 1,..., J (4 considering periodic boundary conditions ( 0 = J, 1 and un+1 1, n N {0} with un+1 J 1 J+1 α n α n (θ := σ 2 θ f ( v n+1 β n β n (θ := 1 γ n γ n (θ := σ 2 θ f ( v n+1 σ 2 θ f ( ṽ n+1 σ ( 2 θ f ṽ n+1
13 Inversibility and positive definiteness of the matrices We can apply the following Proposition If the square matrix of the linear system is strictly diagonally dominant ( α n + γ n < β n,, and also it is its transpose matrix ( + <,, and the diagonal elements α n +1 γ n 1 β n are positive (β n > 0, then the original matrix not only is invertible, but also definite positive. obtaining a restriction over the time step of the following type: t < x f
14 Lax-Friedrichs θ-schemes In this case the tridiagonal system reads α n α n (θ := θ ( ( σ f v n θ ( 2 β n β n (θ := 1 + θ γ n γ n (θ := θ ( ( σ f v n θ ( 2 σ f ( ṽ n+1 σ f ( ṽ n We also can see that the system (for θ 0 is not symmetric, since γ n := θ ( ( σ f v n α+1 n := θ ( ( σ f v n γ n α+1 n = σ ( ( ( 2 θ f v n+1 + f v n+1 n+1 +1 = σ θ f ( v +1/2 [ ] for some v n+1 +1/2 v n+1, v n+1 +1 (if we suppose f continuous.
15 Lax-Friedrichs θ-schemes But always that α n and γ n have the same sign: α n, γn 0 (because 0 α n, γn is not possible since it would imply the ( ( following contradiction σ f v n+1 1 < 1 σ f v n+1 ; ( so 1 σ f x 1 (equivalent to t it will v n+1 also be verified that 0 α n + γ n = θ < 1 + θ = β n f ( v n+1 β n whereas in the case that α n and γ n have distinct sign, it would be necessary to study the situation (suposing in any case θ 0, because the case θ = 0 is trivial, since then the scheme is reduces to the explicit case.
16 Lax-Friedrichs θ-schemes Proposition In order that the matrix of coefficients of the linear system that appears in the application of this schemes, would be strictly diagonally dominant it would suffice that: As 1+θ θ or t x f or always that t > t < 1+θ θ f ( x v n+1 ( x f v n+1 for some, then > 1 for every θ ]0, 1], a sufficient condition of invertibility of these matrices would be the following bound t < 1 + θ θ x f (5
17 Lax-Friedrichs θ-schemes Using the Proposition 1, we can also study the diagonal dominance of the corresponding transpose matrix, Case Always that α+1 n and γn 1 will not have the same sign and ( ( ( f v n+1 + f = 2 f will not vanish then we will 1 v n+1 +1 v n+1 have a restriction over the time step than can be used in order to maintain the strict diagonal dominance of the transpose matrix of the implicit scheme of the form t < ( θ f ( x v n+1
18 Lax-Friedrichs θ-schemes Case But if α+1 n and γn 1 may have both the same sign then this restriction will take the following forms ( t < x (, if α n θ f v n+1 +1, γn 1 0 t < 1 θ x f ( v n+1, if α n +1, γn 1 0 It would be also interesting to investigate now in which practical situations of the application of this implicit scheme to a concrete problem we would have to deal with each of these cases, but it would be the subect of other part of the study, when we will apply this scheme to concrete examples.
19 Rusanov or Local Lax-Friedrichs θ-schemes We can rewrite now this implicit scheme σ ( ( 2 θ f v n+1 ( ζ n+1 1/ σ ( 2 θ ζ n+1 + σ ( ( 2 θ f v n+1 +1/2 + ζn+1 1/2 ζ n+1 +1/2 +1 = δ n (θ f ( v n ( u+1 u n (1 θ σ 1 n un ( 2 +ζ 1/2 n u n u 1 n ( u+1 n un ζ n +1/2
20 Rusanov or Local Lax-Friedrichs θ-schemes Or in a tridiagonal way, α n 1 +βn +γ n +1 = δn (θ, n = 0, 1,... ; = 0, 1,..., J identifying appropriately the coefficients α n α n (θ := σ ( ( 2 θ f v n+1 ζ n+1 1/2 σ ( ( 2 θ f ṽ n+1 β n β n (θ := 1 + σ ( 2 θ ζ n+1 +1/2 + ζn+1 1/2 γ n γ n (θ := σ ( ( 2 θ f v n+1 σ ( 2 θ such that, for all and n ζ n+1 +1/2 f ( ṽ n+1 0 < 1 = α n (θ + β n (θ + γ n (θ β n (θ, θ [0, 1] ζ n+1 1/2 ζ n+1 +1/2
21 Rusanov or Local Lax-Friedrichs θ-schemes We obtain easily the relations, for all and n 0 < 1 = α n (θ + β n (θ + γ n (θ β n (θ, θ [0, 1] 0 < β n (θ = β n (θ and α n (θ + γ n (θ 0; in such a way that, if we know the no possitivity of both coefficients, α n, γn 0 0 α n + γ n = α n γ n = σ ( 2 θ < 1 + σ 2 θ ( ζ n+1 +1/2 + ζn+1 1/2 = β n = ζ n+1 1/2 + ζn+1 +1/2 we also could assure the strict diagonally dominance of the matrices involved, so that we also have their invertibility and the good setting out of the implicit scheme and a positivity result (positive definite matrix. β n
22 Rusanov or Local Lax-Friedrichs θ-schemes If we know that α (θ and γ (θ have the same sign, we would have that Case In the case that both are non positive 0 α n (θ + γ n (θ = σ ( 2 θ ζ n+1 < 1 + σ ( 2 θ ζ n+1 1/2 + ζn+1 +1/2 = β = β 1/2 + ζn+1 +1/2 having assured in this way the good setting out of the corresponding implicit scheme (for θ ]0, 1], whereas for θ = 0 we would obtain the identity matrix and we would recover the usual explicit scheme.
23 Rusanov or Local Lax-Friedrichs θ-schemes Case In the case that one of then were non negative (for example α n (θ 0 we can see easily that this would imply the non positiveness of the other because ( α n (θ 0 f v n+1 ( f v n+1 ζ n+1 1/2 0 0 ζ n+1 +1/2 γn (θ 0 (and the same can be easily proved for the implication γ n (θ 0 = α n (θ 0 so that if we assume the non negativeness of both coefficients, this will lead to a contradiction, unless both vanish, that also imply that 1 = un+1 = +1 would be a sonic point of the flux function (where f (q vanish.
24 Rusanov or Local Lax-Friedrichs θ-schemes So, we can also ask to ourselves the following question, In which situations we will have that α n and γ n take values of distinct sign? We have done an exhaustive study of this question, arriving to the following bounds for the time steps to be taken in these cases (having assured previously the not vanishing of the denominator: case 0 < γ n then t < ( θ (f x v n+1 ζ n+1 +1/2 case 0 < α n we must have then t < ( θ ( f x v n+1 ζ n+1 1/2
25 Index 1 Introduction 2 Implicit Lax-Friedrichs and local Lax-Friedrichs schemes in the scalar case 3 Implicit Lax-Friedrichs and local Lax-Friedrichs schemes for systems
26 Implicit Lax-Friedrichs schemes In the case of hyperbolic systems of conservation laws with source term t U + x f (U = g (U, x U R m the expression of the L-F and local L-F (Rusanov schemes would have the same form that in the scalar case and ĝ (U n 1, Un, Un +1 ;σ should relate adequately with the original function g in the source term. U n+1 = U n σ ( fn f n + ĝ (U n , Un, Un +1 ;σ 2 2 using the corresponding numerical fluxes.
27 Implicit Lax-Friedrichs schemes The well-known Lax-Friedrichs numerical flux reads flf (u, v = 1 (f (u + f (v 1σ 2 (v u (6 Whereas the Rusanov s numerical flux will read as follows fru (u, v = 1 (f (u + f (v ζ (v u, (7 2 where ζ := max (ρ (Df (q, for q [u, v], the segment oining both states. ρ (Df (q is the corresponding espectral radius of the Jacobian matrix of f evaluated in q, that we will also notates Df (q A (q.
28 Implicit Lax-Friedrichs schemes The implicit version for systems of the Lax-Friedichs type schemes would write as usual ( U n+1 + σ fn+1 f n+1 = U n ĝ (U n 1, Un, Un +1 ;σ 2 2 (8 or, if we also evaluate the source terms ĝ at time n + 1 and not at time n ( ( U n+1 + σ fn+1 f n+1 + ĝ U n , Un+1, U n+1 +1 ;σ = U n (9 2 2 This last option (9 will leads to the resolution of large systems of non linear equations in each time step, and will not be considered for the moment.
29 Implicit Lax-Friedrichs schemes On the other hand, the first option (8 can be developed using the expression (6 or (7, of the corresponding numerical flux = U n+1 + σ ( f n+1 2 = ( σ ( 2 A taking V n+1 ( f U n fn+1 1 V n+1 U n + ĝ (U n 1, Un, Un +1 ;σ σ ( 2 ζ U n Un+1 + U n+1 1 σ 2 ζ I ( σ ( + 2 A U n+1 1 V n+1 + (1 + σ ζ Un+1 σ 2 ζ I the corresponding Roe s state, satisfying ( ( ( f U n+1 1 = A V n+1 U n+1 +1 Un+1 1 U n+1 +1
30 Implicit Lax-Friedrichs type schemes So, in order to try to show the inversibility of the corresponding matrices associated to the application of these implicit schemes, following the ideas of partitioned matrices explained in the book of Varga: Geršgorin and His Circles (see also the classical articles of Feingold&Varga-1962 and Varah-1972, as partitioned as follows (not considering into account for the moment boundary conditions A A, 1 A, A,
31 Implicit Lax-Friedrichs schemes where, these m m matrix blocks are the following (for adequate subindexes A, := (1 + θ σ ζ I ( A, 1 := σ 2 (A θ ( A,+1 := σ 2 (A θ V n+1 V n+1 ζ I + ζ I so that, to have the strictly block diagonally dominance of the matrix A (according to the definition in Varga it will suffice to take any of the matrix-norms usually taken, that come from the corresponding vector norms for p = 1, 2 or respectively.
32 Partitioned matrices The interesting thing is that this prominent role of the diagonal entries of a matrix can be generalized also for block matrices, through the use of partitions of C n (or R n in the case of real matrices. By a such partition we mean a direct decomposition in pairwise disoint linear subspaces, each having at least dimension one, whose direct sum is still C n (or R n. We will consider a general partition of C n = W 1... W l with p { p } l =0 such that the nonnegative integers (p N satisfy that p 0 := 0 < p 1 < p 2 <... < p l := n assuming that, for all L {1, 2,..., l}, W = span { e k : p k p } (where e k ( δ k,1,..., δ k,n denotes the usual column basis vectors in C n or R n ; and the Kronecker symbol δ k, := 1 only if k = while vanish otherwise.
33 Partitioned matrices So, given any matrix A ( a i, C n n and a partition p { } l p =0 of Cn = W 1... W l, the matrix results partitioned according to p as A 1,1 A 1,2... A 1,l A 2,1 A 2,2... A 2,l A =.. ( A i, A l,1 A l,2... A l,l i, L (10 Each block submatrix A i, C (p i p i 1 (p p 1 still can represents a linear transformation from W to W i for every i, L (with dim W = p p 1 1, L.
34 Partitioned matrices We can also consider a norm l-tuple φ (φ 1,..., φ l, according to the given partition p { } l p =0, where each φ ( L is a norm on the corresponding subspace W. We will denote the collection of all such norm l-tuples associated with p simply by Φ p, so that we will also have the usual matrix norms associated, for any, k L ( A,k φ φ A,k x ( = sup = sup φ A,k x (11 x W k φ k (x x W k x 0 φ k (x=1 Anyway, we always can consider the more used matrix norms, that come from the well-known corresponding vector norms for p = 1, 2 or, respectively.
35 Partitioned matrices For each block diagonal submatrix A i,i (i L we have m ( ( φ i Ai,i x A i,i := inf (12 x W i φ i (x x 0 so that if A i,i is non-singular then m ( A A i,i = 1/ 1 and φ m ( A i,i = 0 otherwise. This quantity (12 can be also called the reciprocal norm of A i,i. We set for any i L (with the convention that r φ (A := 0 if l = 1 r φ i,p (A := L\{i} i,p i,i A i, φ. (13
36 Partitioned matrices We can now give the main generalizations by Feingol&Varga in the case of partitioned matrices Definition Given a partition p of C n and a norm l-tuple φ Φ p then a matrix A = ( A i, partitioned following p will be said strictly i, L block diagonally dominant with respect to φ if Theorem m ( A i,i > r φ i,p (A, i L. (14 Given a matrix A = ( A i, partitioned following the partition i, L p of C n and a norm l-tuple φ Φ p then, if A is strictly block diagonally dominant with respect to φ we can also assure that A is nonsingular.
37 Implicit Lax-Friedrichs schemes In the L-F case (using now that m ( A, = 1 + σ θ ζ = 1 + θ, and, for p { } l p := {0, m,..., m} and choosing the l-tuple =0 norm φ (φ 1,..., φ l, l copies of any of the usual vector and matrix norms (for p = 1, 2 or then = θ σ 2 ( A ( V n+1 r φ,p (A := k L\{} A,k φ ( + ζ I + A V n+1 ζ I p p θ σ ( ( 2 A V n+1 p + 2ζ I 2 p ( ( p = θ 1 + σ A V n+1
38 Implicit local Lax-Friedrichs θ-schemes Now the m m matrix blocks are the following (for adequate subindexes ( A, := (1 + σ 2 θ ζ n+1 +1/2 + ζn+1 1/2 I ( A, 1 := σ 2 (A θ V n+1 + ζ n+1 ( 1/2 I A,+1 := σ 2 (A θ V n+1 ζ n+1 +1/2 I Now m ( σ ( A, = θ ζ n+1 +1/2 + ζn+1 1/2, obtaining ust the same restriction as in the case of the Lax-Friedrichs θ-scheme for systems t < 1 θ max x ( A V n+1 p.
39 Conclusions An exhaustive theoretical study of well-posed implicit θ-schemes of Lax-Friedrichs type. Invertibility and positive definiteness of the matrices in the scalar case. With explicit time step restrictions. That can be also determined for concrete problems. Use of partitioned matrices and generalizations of the Geršgorin circles theorem in the system case. The source terms do not pose special problems a priori. We are doing also many practical calculations to see other subtleties in concrete problems. All this altogether will complete the present theoretical study.
40 Bibliography D. G Feingold, R. S. Varga, Block Diagonally Dominant Matrices and Generalizations of the Gerschgorin Circle Theorem, Pacific Journal of Mathematics, Vol. 12 (1962, pp E. Godlewski, P. A. Raviart, Numerical Approximation of Hyperbolic Systems of Conservation Laws. Applied Mathematical Sciences-Springer (1996. R. J. Leveque, Finite Volume Methods for Hyperbolic Problems. Cambridge Texts in Applied Mathematics-Cambridge Univ. Press (2003. R. Varga, Gerschgorin and His Circles, Springer Series in Computational Mathematics, 36 (2004. J. M. Varah, On the Solution of Block-Tridiagonal Systems Arising from Certain Finite-Difference Equations, Mathematics of Computation, Vol. 26 Number 120 (1972, pp
41 End THANK YOU VERY MUCH FOR YOUR ATTENTION!
Foundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationChapter 2: Matrix Algebra
Chapter 2: Matrix Algebra (Last Updated: October 12, 2016) These notes are derived primarily from Linear Algebra and its applications by David Lay (4ed). Write A = 1. Matrix operations [a 1 a n. Then entry
More informationFUNCTIONAL ANALYSIS LECTURE NOTES: COMPACT SETS AND FINITE-DIMENSIONAL SPACES. 1. Compact Sets
FUNCTIONAL ANALYSIS LECTURE NOTES: COMPACT SETS AND FINITE-DIMENSIONAL SPACES CHRISTOPHER HEIL 1. Compact Sets Definition 1.1 (Compact and Totally Bounded Sets). Let X be a metric space, and let E X be
More informationREAL RENORMINGS ON COMPLEX BANACH SPACES
REAL RENORMINGS ON COMPLEX BANACH SPACES F. J. GARCÍA PACHECO AND A. MIRALLES Abstract. In this paper we provide two ways of obtaining real Banach spaces that cannot come from complex spaces. In concrete
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More informationMatrices and Matrix Algebra.
Matrices and Matrix Algebra 3.1. Operations on Matrices Matrix Notation and Terminology Matrix: a rectangular array of numbers, called entries. A matrix with m rows and n columns m n A n n matrix : a square
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationMatrices and Vectors. Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A =
30 MATHEMATICS REVIEW G A.1.1 Matrices and Vectors Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A = a 11 a 12... a 1N a 21 a 22... a 2N...... a M1 a M2... a MN A matrix can
More informationThere are two things that are particularly nice about the first basis
Orthogonality and the Gram-Schmidt Process In Chapter 4, we spent a great deal of time studying the problem of finding a basis for a vector space We know that a basis for a vector space can potentially
More information2. Linear algebra. matrices and vectors. linear equations. range and nullspace of matrices. function of vectors, gradient and Hessian
FE661 - Statistical Methods for Financial Engineering 2. Linear algebra Jitkomut Songsiri matrices and vectors linear equations range and nullspace of matrices function of vectors, gradient and Hessian
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationMAT2342 : Introduction to Applied Linear Algebra Mike Newman, fall Projections. introduction
MAT4 : Introduction to Applied Linear Algebra Mike Newman fall 7 9. Projections introduction One reason to consider projections is to understand approximate solutions to linear systems. A common example
More informationVector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition)
Vector Space Basics (Remark: these notes are highly formal and may be a useful reference to some students however I am also posting Ray Heitmann's notes to Canvas for students interested in a direct computational
More information(Linear Programming program, Linear, Theorem on Alternative, Linear Programming duality)
Lecture 2 Theory of Linear Programming (Linear Programming program, Linear, Theorem on Alternative, Linear Programming duality) 2.1 Linear Programming: basic notions A Linear Programming (LP) program is
More informationGaussian elimination
Gaussian elimination October 14, 2013 Contents 1 Introduction 1 2 Some definitions and examples 2 3 Elementary row operations 7 4 Gaussian elimination 11 5 Rank and row reduction 16 6 Some computational
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
More informationFirst, we review some important facts on the location of eigenvalues of matrices.
BLOCK NORMAL MATRICES AND GERSHGORIN-TYPE DISCS JAKUB KIERZKOWSKI AND ALICJA SMOKTUNOWICZ Abstract The block analogues of the theorems on inclusion regions for the eigenvalues of normal matrices are given
More informationReview Let A, B, and C be matrices of the same size, and let r and s be scalars. Then
1 Sec 21 Matrix Operations Review Let A, B, and C be matrices of the same size, and let r and s be scalars Then (i) A + B = B + A (iv) r(a + B) = ra + rb (ii) (A + B) + C = A + (B + C) (v) (r + s)a = ra
More informationLinear Algebra. The analysis of many models in the social sciences reduces to the study of systems of equations.
POLI 7 - Mathematical and Statistical Foundations Prof S Saiegh Fall Lecture Notes - Class 4 October 4, Linear Algebra The analysis of many models in the social sciences reduces to the study of systems
More informationImplicit Scheme for the Heat Equation
Implicit Scheme for the Heat Equation Implicit scheme for the one-dimensional heat equation Once again we consider the one-dimensional heat equation where we seek a u(x, t) satisfying u t = νu xx + f(x,
More informationMATH 426, TOPOLOGY. p 1.
MATH 426, TOPOLOGY THE p-norms In this document we assume an extended real line, where is an element greater than all real numbers; the interval notation [1, ] will be used to mean [1, ) { }. 1. THE p
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space 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 An inner product
More informationThe extreme rays of the 5 5 copositive cone
The extreme rays of the copositive cone Roland Hildebrand March 8, 0 Abstract We give an explicit characterization of all extreme rays of the cone C of copositive matrices. The results are based on the
More informationElectromagnetic Modeling and Simulation
Electromagnetic Modeling and Simulation Erin Bela and Erik Hortsch Department of Mathematics Research Experiences for Undergraduates April 7, 2011 Bela and Hortsch (OSU) EM REU 2010 1 / 45 Maxwell s Equations
More informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationSome Remarks on the Discrete Uncertainty Principle
Highly Composite: Papers in Number Theory, RMS-Lecture Notes Series No. 23, 2016, pp. 77 85. Some Remarks on the Discrete Uncertainty Principle M. Ram Murty Department of Mathematics, Queen s University,
More informationMATH 2331 Linear Algebra. Section 2.1 Matrix Operations. Definition: A : m n, B : n p. Example: Compute AB, if possible.
MATH 2331 Linear Algebra Section 2.1 Matrix Operations Definition: A : m n, B : n p ( 1 2 p ) ( 1 2 p ) AB = A b b b = Ab Ab Ab Example: Compute AB, if possible. 1 Row-column rule: i-j-th entry of AB:
More informationMATH2210 Notebook 2 Spring 2018
MATH2210 Notebook 2 Spring 2018 prepared by Professor Jenny Baglivo c Copyright 2009 2018 by Jenny A. Baglivo. All Rights Reserved. 2 MATH2210 Notebook 2 3 2.1 Matrices and Their Operations................................
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationIn particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with
Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient
More informationLINEAR ALGEBRA: THEORY. Version: August 12,
LINEAR ALGEBRA: THEORY. Version: August 12, 2000 13 2 Basic concepts We will assume that the following concepts are known: Vector, column vector, row vector, transpose. Recall that x is a column vector,
More informationOptimization Theory. A Concise Introduction. Jiongmin Yong
October 11, 017 16:5 ws-book9x6 Book Title Optimization Theory 017-08-Lecture Notes page 1 1 Optimization Theory A Concise Introduction Jiongmin Yong Optimization Theory 017-08-Lecture Notes page Optimization
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationChapter 2 Notes, Linear Algebra 5e Lay
Contents.1 Operations with Matrices..................................1.1 Addition and Subtraction.............................1. Multiplication by a scalar............................ 3.1.3 Multiplication
More informationCopositive Plus Matrices
Copositive Plus Matrices Willemieke van Vliet Master Thesis in Applied Mathematics October 2011 Copositive Plus Matrices Summary In this report we discuss the set of copositive plus matrices and their
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More information1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )
Direct Methods for Linear Systems Chapter Direct Methods for Solving Linear Systems Per-Olof Persson persson@berkeleyedu Department of Mathematics University of California, Berkeley Math 18A Numerical
More information7.6 The Inverse of a Square Matrix
7.6 The Inverse of a Square Matrix Copyright Cengage Learning. All rights reserved. What You Should Learn Verify that two matrices are inverses of each other. Use Gauss-Jordan elimination to find inverses
More informationA Multi Dimensional Stochastic Differential Equation Model
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn:2319-765x. Volume 10, Issue 2 Ver. V (Mar-Apr. 2014), PP 64-69 N. Hema 1 & Dr. A. Jayalakshmi 2 1 Research Scholar, SCSVMV, Kanchipuram,
More informationProblem List MATH 5143 Fall, 2013
Problem List MATH 5143 Fall, 2013 On any problem you may use the result of any previous problem (even if you were not able to do it) and any information given in class up to the moment the problem was
More information3.2 Iterative Solution Methods for Solving Linear
22 CHAPTER 3. NUMERICAL LINEAR ALGEBRA 3.2 Iterative Solution Methods for Solving Linear Systems 3.2.1 Introduction We continue looking how to solve linear systems of the form Ax = b where A = (a ij is
More informationIrrationality exponent and rational approximations with prescribed growth
Irrationality exponent and rational approximations with prescribed growth Stéphane Fischler and Tanguy Rivoal June 0, 2009 Introduction In 978, Apéry [2] proved the irrationality of ζ(3) by constructing
More informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
More informationCSL361 Problem set 4: Basic linear algebra
CSL361 Problem set 4: Basic linear algebra February 21, 2017 [Note:] If the numerical matrix computations turn out to be tedious, you may use the function rref in Matlab. 1 Row-reduced echelon matrices
More informationCopyrighted Material. Type A Weyl Group Multiple Dirichlet Series
Chapter One Type A Weyl Group Multiple Dirichlet Series We begin by defining the basic shape of the class of Weyl group multiple Dirichlet series. To do so, we choose the following parameters. Φ, a reduced
More informationStability of a Class of Singular Difference Equations
International Journal of Difference Equations. ISSN 0973-6069 Volume 1 Number 2 2006), pp. 181 193 c Research India Publications http://www.ripublication.com/ijde.htm Stability of a Class of Singular Difference
More informationSolving Linear Systems
Solving Linear Systems Iterative Solutions Methods Philippe B. Laval KSU Fall 207 Philippe B. Laval (KSU) Linear Systems Fall 207 / 2 Introduction We continue looking how to solve linear systems of the
More informationMath 408 Advanced Linear Algebra
Math 408 Advanced Linear Algebra Chi-Kwong Li Chapter 4 Hermitian and symmetric matrices Basic properties Theorem Let A M n. The following are equivalent. Remark (a) A is Hermitian, i.e., A = A. (b) x
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationMath Linear Algebra
Math 220 - Linear Algebra (Summer 208) Solutions to Homework #7 Exercise 6..20 (a) TRUE. u v v u = 0 is equivalent to u v = v u. The latter identity is true due to the commutative property of the inner
More informationMatrix Mathematics. Theory, Facts, and Formulas with Application to Linear Systems Theory. Dennis S. Bernstein
Matrix Mathematics Theory, Facts, and Formulas with Application to Linear Systems Theory Dennis S. Bernstein PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Contents Special Symbols xv Conventions, Notation,
More informationMATH SOLUTIONS TO PRACTICE MIDTERM LECTURE 1, SUMMER Given vector spaces V and W, V W is the vector space given by
MATH 110 - SOLUTIONS TO PRACTICE MIDTERM LECTURE 1, SUMMER 2009 GSI: SANTIAGO CAÑEZ 1. Given vector spaces V and W, V W is the vector space given by V W = {(v, w) v V and w W }, with addition and scalar
More informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationMATH 115, SUMMER 2012 LECTURE 12
MATH 115, SUMMER 2012 LECTURE 12 JAMES MCIVOR - last time - we used hensel s lemma to go from roots of polynomial equations mod p to roots mod p 2, mod p 3, etc. - from there we can use CRT to construct
More informationACI-matrices all of whose completions have the same rank
ACI-matrices all of whose completions have the same rank Zejun Huang, Xingzhi Zhan Department of Mathematics East China Normal University Shanghai 200241, China Abstract We characterize the ACI-matrices
More informationThe following definition is fundamental.
1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic
More informationProblems in Linear Algebra and Representation Theory
Problems in Linear Algebra and Representation Theory (Most of these were provided by Victor Ginzburg) The problems appearing below have varying level of difficulty. They are not listed in any specific
More informationChapter 1. Vectors, Matrices, and Linear Spaces
1.6 Homogeneous Systems, Subspaces and Bases 1 Chapter 1. Vectors, Matrices, and Linear Spaces 1.6. Homogeneous Systems, Subspaces and Bases Note. In this section we explore the structure of the solution
More informationMODULE 8 Topics: Null space, range, column space, row space and rank of a matrix
MODULE 8 Topics: Null space, range, column space, row space and rank of a matrix Definition: Let L : V 1 V 2 be a linear operator. The null space N (L) of L is the subspace of V 1 defined by N (L) = {x
More informationFinite Differences: Consistency, Stability and Convergence
Finite Differences: Consistency, Stability and Convergence Varun Shankar March, 06 Introduction Now that we have tackled our first space-time PDE, we will take a quick detour from presenting new FD methods,
More information1 Determinants. 1.1 Determinant
1 Determinants [SB], Chapter 9, p.188-196. [SB], Chapter 26, p.719-739. Bellow w ll study the central question: which additional conditions must satisfy a quadratic matrix A to be invertible, that is to
More informationLecture 6: Geometry of OLS Estimation of Linear Regession
Lecture 6: Geometry of OLS Estimation of Linear Regession Xuexin Wang WISE Oct 2013 1 / 22 Matrix Algebra An n m matrix A is a rectangular array that consists of nm elements arranged in n rows and m columns
More information10.2 ITERATIVE METHODS FOR SOLVING LINEAR SYSTEMS. The Jacobi Method
54 CHAPTER 10 NUMERICAL METHODS 10. ITERATIVE METHODS FOR SOLVING LINEAR SYSTEMS As a numerical technique, Gaussian elimination is rather unusual because it is direct. That is, a solution is obtained after
More informationNumerical Linear Algebra
Numerical Linear Algebra The two principal problems in linear algebra are: Linear system Given an n n matrix A and an n-vector b, determine x IR n such that A x = b Eigenvalue problem Given an n n matrix
More informationMATH 581D FINAL EXAM Autumn December 12, 2016
MATH 58D FINAL EXAM Autumn 206 December 2, 206 NAME: SIGNATURE: Instructions: there are 6 problems on the final. Aim for solving 4 problems, but do as much as you can. Partial credit will be given on all
More informationIterative Methods. Splitting Methods
Iterative Methods Splitting Methods 1 Direct Methods Solving Ax = b using direct methods. Gaussian elimination (using LU decomposition) Variants of LU, including Crout and Doolittle Other decomposition
More informationNotes on Linear Algebra and Matrix Theory
Massimo Franceschet featuring Enrico Bozzo Scalar product The scalar product (a.k.a. dot product or inner product) of two real vectors x = (x 1,..., x n ) and y = (y 1,..., y n ) is not a vector but a
More informationA 3 3 DILATION COUNTEREXAMPLE
A 3 3 DILATION COUNTEREXAMPLE MAN DUEN CHOI AND KENNETH R DAVIDSON Dedicated to the memory of William B Arveson Abstract We define four 3 3 commuting contractions which do not dilate to commuting isometries
More informationLecture Notes for Inf-Mat 3350/4350, Tom Lyche
Lecture Notes for Inf-Mat 3350/4350, 2007 Tom Lyche August 5, 2007 2 Contents Preface vii I A Review of Linear Algebra 1 1 Introduction 3 1.1 Notation............................... 3 2 Vectors 5 2.1 Vector
More informationHarmonic Polynomials and Dirichlet-Type Problems. 1. Derivatives of x 2 n
Harmonic Polynomials and Dirichlet-Type Problems Sheldon Axler and Wade Ramey 30 May 1995 Abstract. We take a new approach to harmonic polynomials via differentiation. Surprisingly powerful results about
More informationCourse Summary Math 211
Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.
More informationMATH 221: SOLUTIONS TO SELECTED HOMEWORK PROBLEMS
MATH 221: SOLUTIONS TO SELECTED HOMEWORK PROBLEMS 1. HW 1: Due September 4 1.1.21. Suppose v, w R n and c is a scalar. Prove that Span(v + cw, w) = Span(v, w). We must prove two things: that every element
More informationFoundations of Analysis. Joseph L. Taylor. University of Utah
Foundations of Analysis Joseph L. Taylor University of Utah Contents Preface vii Chapter 1. The Real Numbers 1 1.1. Sets and Functions 2 1.2. The Natural Numbers 8 1.3. Integers and Rational Numbers 16
More informationTopics in linear algebra
Chapter 6 Topics in linear algebra 6.1 Change of basis I want to remind you of one of the basic ideas in linear algebra: change of basis. Let F be a field, V and W be finite dimensional vector spaces over
More informationM. Matrices and Linear Algebra
M. Matrices and Linear Algebra. Matrix algebra. In section D we calculated the determinants of square arrays of numbers. Such arrays are important in mathematics and its applications; they are called matrices.
More information1 Number Systems and Errors 1
Contents 1 Number Systems and Errors 1 1.1 Introduction................................ 1 1.2 Number Representation and Base of Numbers............. 1 1.2.1 Normalized Floating-point Representation...........
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationChapter 4 - MATRIX ALGEBRA. ... a 2j... a 2n. a i1 a i2... a ij... a in
Chapter 4 - MATRIX ALGEBRA 4.1. Matrix Operations A a 11 a 12... a 1j... a 1n a 21. a 22.... a 2j... a 2n. a i1 a i2... a ij... a in... a m1 a m2... a mj... a mn The entry in the ith row and the jth column
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationAnalysis II: The Implicit and Inverse Function Theorems
Analysis II: The Implicit and Inverse Function Theorems Jesse Ratzkin November 17, 2009 Let f : R n R m be C 1. When is the zero set Z = {x R n : f(x) = 0} the graph of another function? When is Z nicely
More informationJim Lambers MAT 610 Summer Session Lecture 1 Notes
Jim Lambers MAT 60 Summer Session 2009-0 Lecture Notes Introduction This course is about numerical linear algebra, which is the study of the approximate solution of fundamental problems from linear algebra
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationDIFFERENTIAL GEOMETRY Multivariable Calculus, a refresher
DIFFERENTIAL GEOMETRY Multivariable Calculus, a refresher 1 Preliminaries In this section I collect some results from multivariable calculus that play an important role in the course. If you don t know
More informationContents. 1 Vectors, Lines and Planes 1. 2 Gaussian Elimination Matrices Vector Spaces and Subspaces 124
Matrices Math 220 Copyright 2016 Pinaki Das This document is freely redistributable under the terms of the GNU Free Documentation License For more information, visit http://wwwgnuorg/copyleft/fdlhtml Contents
More informationIowa State University. Instructor: Alex Roitershtein Summer Exam #1. Solutions. x u = 2 x v
Math 501 Iowa State University Introduction to Real Analysis Department of Mathematics Instructor: Alex Roitershtein Summer 015 Exam #1 Solutions This is a take-home examination. The exam includes 8 questions.
More informationSolution Set 7, Fall '12
Solution Set 7, 18.06 Fall '12 1. Do Problem 26 from 5.1. (It might take a while but when you see it, it's easy) Solution. Let n 3, and let A be an n n matrix whose i, j entry is i + j. To show that det
More information1 Finite difference example: 1D implicit heat equation
1 Finite difference example: 1D implicit heat equation 1.1 Boundary conditions Neumann and Dirichlet We solve the transient heat equation ρc p t = ( k ) (1) on the domain L/2 x L/2 subject to the following
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationLinearizing Symmetric Matrix Polynomials via Fiedler pencils with Repetition
Linearizing Symmetric Matrix Polynomials via Fiedler pencils with Repetition Kyle Curlett Maribel Bueno Cachadina, Advisor March, 2012 Department of Mathematics Abstract Strong linearizations of a matrix
More informationIntroduction to Matrices
POLS 704 Introduction to Matrices Introduction to Matrices. The Cast of Characters A matrix is a rectangular array (i.e., a table) of numbers. For example, 2 3 X 4 5 6 (4 3) 7 8 9 0 0 0 Thismatrix,with4rowsand3columns,isoforder
More informationA classification of sharp tridiagonal pairs. Tatsuro Ito, Kazumasa Nomura, Paul Terwilliger
Tatsuro Ito Kazumasa Nomura Paul Terwilliger Overview This talk concerns a linear algebraic object called a tridiagonal pair. We will describe its features such as the eigenvalues, dual eigenvalues, shape,
More informationSAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra
SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to 1.1. Introduction Linear algebra is a specific branch of mathematics dealing with the study of vectors, vector spaces with functions that
More informationLinear Algebra. Solving Linear Systems. Copyright 2005, W.R. Winfrey
Copyright 2005, W.R. Winfrey Topics Preliminaries Echelon Form of a Matrix Elementary Matrices; Finding A -1 Equivalent Matrices LU-Factorization Topics Preliminaries Echelon Form of a Matrix Elementary
More information