Numerical Solutions to Partial Differential Equations
|
|
- Neil Potter
- 5 years ago
- Views:
Transcription
1 Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University
2 A Model Problem in a 2D Box Region Let us consider a model problem of parabolic equation: u t = a u xx + u yy, x, y Ω, t > 0, ux, y, 0 = u 0 x, y, 0, x, y Ω, ux, y, t = 0, x, y Ω, t > 0, where a > 0 is a constant, Ω = 0, X 0, Y R 2. 1 For integers N x 1 and N y 1, let h x = x = XNx 1 and h y = y = YNy 1 be the grid sizes in the x and y directions; 2 a uniform parallelopiped grid with the set of grid nodes J Ω R+ = {x j, y k, t m : 0 j N x, 0 k N y, m 0}, where x j = j h x, y k = k h y, t m = mτ τ > 0 time step size. 3 the space of grid functions U = {U m = Ux j, y k, t m : 0 j N x, 0 k N y, m 0}.
3 Explicit and Implicit Schemes in Box Regions The Explicit Scheme in a 2D Box Region The Forward Explicit Scheme in a 2D Box Region The forward explicit scheme and its error equation U m+1 U m [ U m j+1,k 2U m = a + Um j 1,k τ hx 2 + Um +1 2Um + ] Um 1 hy 2. e m+1 = [1 2µ x + µ y ] e+µ m x e m j+1,k + ej 1,k m +µy e m +1 + e 1 m T m where µ x = a τ, µ hx 2 y = a τ hy 2 are the grid ratios in x and y directions. The truncation error is Tux, y, t = 1 2 u ttx, y, tτ a 4 12 x ux, y, thx 2 + y 4 ux, y, thy 2 + Oτ 2 + hx 4 + hy 4. τ, 3 / 35
4 Explicit and Implicit Schemes in Box Regions The Explicit Scheme in a 2D Box Region The Forward Explicit Scheme in a 2D Box Region 1 the condition for the maximum principle: µ x + µ y 1 2 ; 2 Fourier modes and amplification factors λ l l = l x, l y : U m = λm l e iαx x j +α y y k = λ m l e ilx jπn 1 x +l y kπny 1, α x = lx π X, N x + 1 l x N x, α y = ly π Y, N y + 1 l y N y ; l x, l y represent the frequency or wave number in x, y direction. [ ] 3 2 αx hx amplification factor λ l = 1 4 µ x sin 2 + µ y sin 2 α y h y 2 [ = 1 4 µ x sin 2 lx π 2N x + µ y sin 2 l y π 2N x ]; 4 L 2 stable if and only if µ x + µ y 1/2; 5 Convergence rate is Oτ + h 2 x + h 2 y. 4 / 35
5 The θ-scheme in a 2D Box Region [ ] [ ] U m+1 U m δx 2 = 1 θa τ hx 2 + δ2 y hy 2 U m δx 2 + θa hx 2 + δ2 y hy 2 [ U m j+1,k 2U m = 1 θa + Um j 1,k + Um +1 2Um + ] Um 1 + θa h 2 x [ U m+1 j+1,k 2Um+1 + U m+1 j 1,k hx 2 + Um+1 The error equation: [1 + 2θµ x + µ y ] e m+1 h 2 y +1 2Um+1 + U m+1 h 2 y 1 U m+1 ] = θ [µ x e m+1 j+1,k + em+1 j 1,k + µ y e m em [1 21 θµ x + µ y ] e m + 1 θ [ µ x e m j+1,k + e m j 1,k + µy e m +1 + e m 1] τt m+ 1 2, ].
6 Explicit and Implicit Schemes in Box Regions The θ-scheme in a 2D Box Region The θ-scheme in a 2D Box Region The truncation error T m+ j = { Oτ 2 + h 2 x + h 2 y, if θ = 1 2, Oτ + h 2 x + h 2 y, if θ 1 2, 1 the condition for the maximum principle: 2µ x + µ y 1 θ 1; 2 Fourier modes: U m = λm l e iαx x j +α y y k = λ m l e l = l x, l y, α x = lx π X, αy = ly π Y, αx x j = 3 amplification factor 1 41 θ λ l = 1 + 4θ lx jπ Nx, αy y ly kπ k = Ny ; [ µ x sin 2 lx π 2N x [ µ x sin 2 lx π 2N x lx jπ ly kπ i + Nx Ny, ] + µ y sin 2 l y π 2N x + µ y sin 2 l y π 2N x ] ; 6 / 35
7 Explicit and Implicit Schemes in Box Regions The θ-scheme in a 2D Box Region The θ-scheme in a 2D Box Region 4 for θ 1/2, unconditionally L 2 stable; 5 for 0 θ < 1/2, L 2 stable iff 21 2θµ x + µ y 1; 6 the matrix of the linear system is still symmetric positive definite and diagonal dominant, however, each row has now up to 5 nonzero elements with a band width of the order Oh 1 ; 7 if solved by the Thompson method, the cost is Oh 1 times of that of the explicit scheme; 8 in 3D, µ x + µ y µ x + µ y + µ z, and Oh 1 Oh 2, if solved by the Thompson method, the cost is Oh 2 times of that of the explicit scheme. 7 / 35
8 2D and 3D ADI and LOD Schemes Alternative Approaches for Solving n-d Parabolic Equations To reduce the computational cost, we may consider to apply highly efficient iterative methods to solve the linear algebraic equations, for example, the preconditioned conjugate gradient method; the multi-grid method; etc.. Alternatively, to avoid the shortcoming of the implicit difference schemes for high space dimensions, we may develop the alternating direction implicit ADI schemes; the locally one dimensional LOD schemes. 8 / 35
9 2D and 3D ADI and LOD Schemes A 2D Alternating Direction Implicit ADI Scheme A Fractional Steps 2D ADI Scheme by Peaceman and Rachford 1 in the odd fractional steps implicit in x and explicit in y, and in even fractional steps implicit in y and explicit in x: µ xδx 2 U m+ 1 2 = µ y δy 2 U, m µ y δ 2 y U m+1 = µ xδ 2 x U m+ 1 2, 2 numerical boundary conditions are easily imposed directly by those of the original problem, since U m+ 1 2 u m+ 1 2 ; 9 / 35
10 2D and 3D ADI and LOD Schemes A 2D Alternating Direction Implicit ADI Scheme A Fractional Steps 2D ADI Scheme by Peaceman and Rachford 3 one step equivalent scheme µ xδx µ y δy 2 U m+1 = µ xδx µ y δy 2 U. m 4 Crank-Nicolson scheme µ xδx µ y δy 2 U m+1 = µ xδx µ y δy 2 U. m 5 Since µ x µ y δ 2 xδ 2 y δ t u m+ 1 2 = a 2 τ 3 [u xxyyt ] m Oτ 5 + τ 3 h 2 x + h 2 y, the truncation error of the scheme is Oτ 2 + h 2 x + h 2 y. 10 / 35
11 2D and 3D ADI and LOD Schemes A 2D Alternating Direction Implicit ADI Scheme Stability of the 2D ADI Scheme by Peaceman and Rachford 1 For the Fourier mode U m = λm l e ilx jπn 1 x +l y kπn 1 y, λ l = 1 2 µ x sin 2 lx π 2N x 1 2 µ y sin 2 ly π 2N y µ x sin 2 lx π 2N x µ y sin 2 ly π 2N y, the scheme is unconditionally L 2 stable; 2 the two fractional steps can be equivalently written as 1 + µ xu m+ 1 2 = 1 µ y U m + µy U m 1 + U m µ x µ y U m+1 = 1 µ xu m µx U m+ 1 2 j 1,k Um+ 2 j+1,k + µy 2 thus, the maximum principle holds if max{µ x, µ y } 1; U m+ 1 2 j 1,k + Um+ 1 2 j+1,k, U m U m+1 +1, 11 / 35
12 2D and 3D ADI and LOD Schemes A 2D Alternating Direction Implicit ADI Scheme Cost of the 2D ADI Scheme by Peaceman and Rachford 1 order the linear system by N x 1 k + j and N y 1 j + k in odd and even steps, the corresponding matrixes are tridiagonal; 2 the computational cost: 3 times of that of the explicit scheme. 12 / 35
13 2D and 3D ADI and LOD Schemes A 2D Alternating Direction Implicit ADI Scheme The Idea of Peaceman and Rachford Doesn t Work for 3D In 3D, we can still construct fractional step scheme by 1 dividing each time step into 3 fractional steps; 2 introducing U m+ 1 3 and U m+ 2 3 at t m+ 1 3 and t m+ 2 ; 3 3 in the 3 fractional time steps, applying schemes which are in turn implicit in x-, y- and z-direction and explicit in the other 2 directions respectively. However, the equivalent one step scheme is, up to a higher order term, the same as the θ-scheme with θ = 1/3, which has a local truncation error Oτ + h 2 x + h 2 y instead of what we expect to have Oτ 2 + h 2 x + h 2 y for an ADI method. 13 / 35
14 The Key Properties that Make an ADI Scheme Successful 1 Implicit only in one dimension in each fractional time step, thus, if properly ordered, the matrix of the linear algebraic equations is tridiagonal as well as diagonally dominant, hence the computational cost is significantly reduced. 2 In the equivalent one step scheme, the implicit part is the product of 1D implicit difference operators, the explicit part is the product of 1D explicit difference operators in a half time step, + certain h.o.t. small perturbations, which guarantees the scheme is unconditionally L 2 stable; 3 The difference between the equivalent one step scheme and the Crank-Nicolson scheme is a higher order term, which guarantees that the local truncation error is Oτ 2 + h 2.
15 2D and 3D ADI and LOD Schemes 3D Alternate Direction Implicit ADI Schemes Target Equivalent One Step Schemes of ADI Methods For the 3D model problem, we aim to develop ADI finite difference schemes which have a one time step equivalent scheme of the form µ xδx µ y δy µ zδz 2 U m+1 = µ xδx µ y δy µ zδz 2 U m, or µ xδx µ y δy µ zδz 2 = µ xδx µ y δy 2 U m µ zδ 2 z U m + h.o.t. 15 / 35
16 2D and 3D ADI and LOD Schemes 3D Alternate Direction Implicit ADI Schemes An Extendable ADI Scheme Proposed by D yakonov For the 3D model problem, the ADI Scheme of D yakonov is 1 µ x 2 δ2 x U m+ = 1 + µ x 2 δ2 x 1 + µ y 2 δ2 y 1 + µ z 2 δ2 z U m, 1 µ y 2 δ2 y U m+ = U m+, 1 µ z 2 δ2 z U m+1 = U m+. The scheme obviously has all the key properties. However, there is one thing need to be taken care of: U m+, U m+ are not approximate solutions, their boundary conditions may not be directly derived from that of the original problem. 16 / 35
17 2D and 3D ADI and LOD Schemes 3D Alternate Direction Implicit ADI Schemes Numerical Boundary Conditions for U m+, U m+ Since U m+1 u m+1, Um+1 = u m+1 on the boundary j = 0, N x, or k = 0, N y, or l = 0, N z. Thus, the boundary condition { U m+ } : j = 0, N x, 0 k N y, 0 < l < N z, or k = 0, N y ; 0 j N x, 0 < l < N z can be obtained from U m+ = 1 µz 2 δ2 z U m+1. The boundary { condition } U m+ : j = 0, N x, 0 < k < N y, 0 < l < N z. can be obtained from U m+ = 1 µy 2 δ2 y U m+ and the boundary values for U m+ on the nodes j = 0, N x. 17 / 35
18 2D and 3D ADI and LOD Schemes 3D Alternate Direction Implicit ADI Schemes An Extendable ADI Scheme Proposed by Douglas and Rachford For the 3D model problem, the Scheme is of the form U m+1 = U m µ xδx 2 U m+1 + U m + µ y δy 2 U m + µ zδz 2 U m, U m+1 = U m µ y δy 2 U m+1 + U m µ y δy 2 U m, U m+1 = U m µ zδz 2 U m+1 +U m µ z δz 2 U m. Since U m+1 and U m+1 are approximate solutions at t m+1, their boundary conditions can be directly derived from that of the original problem.the scheme has all the key properties with the h.o.t.= 1 4 µ xµ y µ z δ 2 xδ 2 y δ 2 z U m. 18 / 35
19 2D and 3D ADI and LOD Schemes 3D Alternate Direction Implicit ADI Schemes The Idea of the ADI Scheme of Douglas and Rachford The construction of the scheme may be viewed as a predictioncorrection process, which may also be explained as follows: In the 1st fractional step, Crank-Nicolson to x, explicit to y, z. In the 2nd fractional step, to improve stability and accuracy in y, the y-direction is corrected by the Crank-Nicolson scheme: U m+1 = U m µ xδx 2 U m+1 + U m µ y δy 2 U m+1 + U m + µ z δz 2 U m. Note, this minus the 1st equation gives the second in the scheme. 19 / 35
20 2D and 3D ADI and LOD Schemes 3D Alternate Direction Implicit ADI Schemes The Idea of the ADI Scheme of Douglas and Rachford In the 3rd, z-direction is corrected by the Crank-Nicolson scheme: U m+1 = U m µ xδ 2 x µ y δ 2 y U m+1 U m+1 + U m + U m µ zδz 2 U m+1 + U m. Note, the 3rd equation in the scheme is simply the difference of the above two equations. Since U m+1 and U m+1 are approximations of U m+1, numerical boundary conditions can be directly derived from those of the original problem. 20 / 35
21 2D and 3D ADI and LOD Schemes 3D Locally One Dimensional LOD Schemes An Extendable Locally One Dimensional LOD Scheme 1 In the ith fractional step, the problem is treated as a one dimensional problem of the ith dimension. 2 1D Crank-Nicolson scheme is applied to each fractional step. 3 For the 3D model problem, the LOD scheme is given as µ xδx 2 U m+ = µ xδx 2 U, m µ y δy 2 U m+ = µ y δy 2 U m+, µ zδ 2 z U m+1 = µ zδ 2 z 4 The scheme has all of the key properties. U m+. 21 / 35
22 2D and 3D ADI and LOD Schemes 3D Locally One Dimensional LOD Schemes Impose Boundary Conditions for the LOD Scheme However, U m+, U m+ are nonphysical, their boundary conditions must be handled with special care. 1 multiplying µ zδz 2 on the 3rd equation leads to µ2 zδz 4 U m+ = 1 µ z δz µ2 zδz 4 U m+1, Hence, by omitting the Oτ 2 order terms, the boundary condition for U m+ can be given by: U m+ = 1 µ z δz 2 U m+1, k = 0, N y, 0 j N x, 0 < l < N z, j = 0, N x, 0 k N y, 0 < l < N z. 22 / 35
23 2D and 3D ADI and LOD Schemes 3D Locally One Dimensional LOD Schemes Impose Boundary Conditions for the LOD Scheme 2 Similarly, multiplying µ y δ 2 y on the 2nd equation, omitting the Oτ 2 and higher order terms, the boundary condition for U m+ can be given by: U m+ = 1 µ y δy 2 U m+, j = 0, N x, 0 < k < N y, 0 < l < N z. 23 / 35
24 FDM for General Parabolic Problems in High Dimensions General Domain and Boundary Conditions General Domain and Boundary Conditions 1 Approximate boundary conditions can be established with the similar methods used in for elliptic problems. 2 Approximate boundary conditions can further restrict the stability condition for explicit schemes, thus implicit schemes are even more preferable. 3 ADI and LOD schemes are in more complicated forms. 4 Some special regions with non-planar boundaries can be transformed into box regions by curvilinear coordinate systems, say, the polar coordinate system can be used for a circular region; and the cylindrical coordinate system can be used for cylindrical regions, etc / 35
25 FDM for General Parabolic Problems in High Dimensions Variable-coefficient Equations Variable-coefficient and nonlinear Equations 1 ADI and LOD schemes can also be extended to high dimensional variable-coefficient linear, and even certain nonlinear problems. 2 In such cases, the difference operators 1 ± 1 2 µ xδx 2 and 1 ± 1 2 µ y δy 2 are generally noncommutative, i.e. 1 ± 1 2 µ xδx 2 1 ± 1 2 µ y δy 2 1 ± 1 2 µ y δy 2 1 ± 1 2 µ xδx 2, and this will introduce the so called splitting error. 3 Boundary conditions are more difficult to handle. 25 / 35
26 FDM for General Parabolic Problems in High Dimensions Variable-coefficient Equations Implicit, Semi-implicit Schemes For nonlinear problems, nonlinear algebraic equations derived from implicit schemes need to be solved with iterative methods, such as semi-implicit schemes. The idea is to apply an implicit scheme only to the principal linear part of the nonlinear equation; approximate the residual nonlinear part by an explicit scheme. 26 / 35
27 FDM for General Parabolic Problems in High Dimensions Variable-coefficient Equations Implicit, Semi-implicit Schemes and Elliptic Solver To apply the implicit schemes to linear parabolic problems and the semi-implicit schemes to nonlinear parabolic problems, it is a sequence of linear algebraic systems corresponding to certain elliptic problems we actually need to solve in the end. So, Fast solvers for elliptic problems play important roles in solving parabolic problems. 27 / 35
28 FDM for General Parabolic Problems in High Dimensions Variable-coefficient Equations Asymptotic Analysis and Extrapolation Methods The asymptotic analysis and extrapolation methods introduced in 1.5 for elliptic problems can also be extended to analyze the finite difference approximation error for parabolic problems. In particular, the error bounds can be estimated by the numerical solutions obtained on grids with different mesh sizes. For example, for the explicit scheme of the heat equation, let τ = µh 2, µ 1/2, then, we have U m hj = um hj + Oh2 + Oh 4, U 1m j 4U4m h/22j U m hj 3 = u m hj + Oh4. 28 / 35
29 Finite Difference Method for Hyperbolic equations Introduction to Hyperbolic Equations n-d 1st Order Linear Hyperbolic Partial Differential Equation: n 1 Scalar case u R 1, u t + a i u xi + b u = ψ 0, i=1 where a i, b and ψ 0 are real functions x = x 1,..., x n and t. 2 Vector case u = u 1,, u p T R p, n u t + A i u xi + B u = ψ 0, i=1 where A i, B R p p, ψ 0 R p are real functions of t and x = x 1,..., x n, and α R n, Ax, t = n i=1 α ia i x, t is real diagonalizable, i.e. Ax, t has p linearly independent eigenvectors corresponding to real eigenvalues. 3 If Ax, t = n i=1 α ia i x, t has p mutually different real eigenvalues, the system is called strictly hyperbolic.
30 Finite Difference Method for Hyperbolic equations Introduction to Hyperbolic Equations n-d 2nd Order Linear Hyperbolic Partial Differential Equation A general 2nd order scalar equation u R 1, n n n u tt + 2 a i u xi t + b 0 u t a ij u xi x j + b i u xi + cu = ψ 0, i=1 i,j=1 where a i, a ij, b i, c and ψ 0 are real functions x = x 1,..., x n and t, a ij is real symmetric positive definite. Define v = u, v 0 = u t, v i = u xi, then the above 2nd order scalar equation transforms into a first order linear system of partial differential equations for v = v, v 0, v 1,, v n T n Av t + A i v xi + Bv = ψ 0, i=1 i=1
31 Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 2nd Order Scalar Transforms to n-d 1st Order System p = n a i a 1i a ni A = 0 0 a 11 a 1n, A i = 0 a 1i 0 0, a n1 a nn 0 a ni c b 0 b 1 b n ψ 0 B = , ψ 0 = / 35
32 Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 2nd Order Scalar Transforms to n-d 1st Order System p = n Let R T AR = I A is real symmetric positive definite; 2 Introduce a new variable w = R 1 v; 3 Denote Âi = R T A i R, ˆB = R T BR, ˆψ 0 = R T ψ; 4 Âx, t = n i=1 α iâix, t is a real symmetric matrix and thus is real diagonizable for all real α i, i = 1,..., n; 5 The 2nd order scalar equation now transforms into a 1st order linear hyperbolic system of partial differential equations for w R n+2 : n w t + Â i w xi + ˆBw = ˆψ 0. i=1 32 / 35
33 Introduction to Hyperbolic Equations The Hyperbolic Equations Standard Form of n-d 1st Order Linear Hyperbolic Equations 1 The standard form of 1st order linear hyperbolic equation: n u t + a i u xi = ψ, ψ = ψ 0 b u. i=1 2 The standard form of 1st order linear hyperbolic system: n u t + A i u xi = ψ, ψ = ψ 0 B u; i=1 3 The equation system is said to be homogeneous, if ψ = 0; 33 / 35
34 Introduction to Hyperbolic Equations The Hyperbolic Equations Standard Form of n-d 1st Order Linear Hyperbolic Equations 4 In general, a higher order linear hyperbolic equation system of equations can always be transformed into a first order linear hyperbolic system of equations. 5 An equation system is said to be nonlinear, if at least one of the coefficients depends on the unknown or the right hand side term is a nonlinear function of the unknown. 34 / 35
35 Thank You! SK 2µ20, 21
Numerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
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 informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 1: Introduction to Multigrid 2000 Eric de Sturler 1 12/02/09 MG01.prz Basic Iterative Methods (1) Nonlinear equation: f(x) = 0 Rewrite as x = F(x), and iterate x i+1
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
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 informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationFinite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More informationNumerical Algorithms for Visual Computing II 2010/11 Example Solutions for Assignment 6
Numerical Algorithms for Visual Computing II 00/ Example Solutions for Assignment 6 Problem (Matrix Stability Infusion). The matrix A of the arising matrix notation U n+ = AU n takes the following form,
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 informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationA CCD-ADI method for unsteady convection-diffusion equations
A CCD-ADI method for unsteady convection-diffusion equations Hai-Wei Sun, Leonard Z. Li Department of Mathematics, University of Macau, Macao Abstract With a combined compact difference scheme for the
More informationNumerical Solution of partial differential equations
G. D. SMITH Brunei University Numerical Solution of partial differential equations FINITE DIFFERENCE METHODS THIRD EDITION CLARENDON PRESS OXFORD Contents NOTATION 1. INTRODUCTION AND FINITE-DIFFERENCE
More informationMulti-Factor Finite Differences
February 17, 2017 Aims and outline Finite differences for more than one direction The θ-method, explicit, implicit, Crank-Nicolson Iterative solution of discretised equations Alternating directions implicit
More informationMaths for Signals and Systems Linear Algebra in Engineering
Maths for Signals and Systems Linear Algebra in Engineering Lectures 13 15, Tuesday 8 th and Friday 11 th November 016 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More information1 Separation of Variables
Jim ambers ENERGY 281 Spring Quarter 27-8 ecture 2 Notes 1 Separation of Variables In the previous lecture, we learned how to derive a PDE that describes fluid flow. Now, we will learn a number of analytical
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 informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationCHAPTER 5. Basic Iterative Methods
Basic Iterative Methods CHAPTER 5 Solve Ax = f where A is large and sparse (and nonsingular. Let A be split as A = M N in which M is nonsingular, and solving systems of the form Mz = r is much easier than
More informationECE539 - Advanced Theory of Semiconductors and Semiconductor Devices. Numerical Methods and Simulation / Umberto Ravaioli
ECE539 - Advanced Theory of Semiconductors and Semiconductor Devices 1 General concepts Numerical Methods and Simulation / Umberto Ravaioli Introduction to the Numerical Solution of Partial Differential
More informationand let s calculate the image of some vectors under the transformation T.
Chapter 5 Eigenvalues and Eigenvectors 5. Eigenvalues and Eigenvectors Let T : R n R n be a linear transformation. Then T can be represented by a matrix (the standard matrix), and we can write T ( v) =
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE
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 informationAn Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems
An Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems Yongbin Ge, 1 Zhen F. Tian, 2 Jun Zhang 3 1 Institute of Applied Mathematics and Mechanics, Ningxia University,
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 informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationCoordinate Update Algorithm Short Course Operator Splitting
Coordinate Update Algorithm Short Course Operator Splitting Instructor: Wotao Yin (UCLA Math) Summer 2016 1 / 25 Operator splitting pipeline 1. Formulate a problem as 0 A(x) + B(x) with monotone operators
More informationFundamentals of Engineering Analysis (650163)
Philadelphia University Faculty of Engineering Communications and Electronics Engineering Fundamentals of Engineering Analysis (6563) Part Dr. Omar R Daoud Matrices: Introduction DEFINITION A matrix is
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationEcon Slides from Lecture 7
Econ 205 Sobel Econ 205 - Slides from Lecture 7 Joel Sobel August 31, 2010 Linear Algebra: Main Theory A linear combination of a collection of vectors {x 1,..., x k } is a vector of the form k λ ix i for
More informationEigenvalues and Eigenvectors
/88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix
More informationAn Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations.
An Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations by Tong Chen A thesis submitted in conformity with the requirements
More informationHomework For each of the following matrices, find the minimal polynomial and determine whether the matrix is diagonalizable.
Math 5327 Fall 2018 Homework 7 1. For each of the following matrices, find the minimal polynomial and determine whether the matrix is diagonalizable. 3 1 0 (a) A = 1 2 0 1 1 0 x 3 1 0 Solution: 1 x 2 0
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 informationLecture 1. Finite difference and finite element methods. Partial differential equations (PDEs) Solving the heat equation numerically
Finite difference and finite element methods Lecture 1 Scope of the course Analysis and implementation of numerical methods for pricing options. Models: Black-Scholes, stochastic volatility, exponential
More informationComputational Methods. Systems of Linear Equations
Computational Methods Systems of Linear Equations Manfred Huber 2010 1 Systems of Equations Often a system model contains multiple variables (parameters) and contains multiple equations Multiple equations
More informationq is large, to ensure accuracy of the spatial discretization 1 In practice, q = 6 is a good choice.
Time-stepping beyond CFL: a locally one-dimensional scheme for acoustic wave propagation Leonardo Zepeda-Núñez 1*), Russell J. Hewett 1, Minghua Michel Rao 2, and Laurent Demanet 1 1 Dept. of Mathematics
More informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More informationIntroduction to Matrix Algebra
Introduction to Matrix Algebra August 18, 2010 1 Vectors 1.1 Notations A p-dimensional vector is p numbers put together. Written as x 1 x =. x p. When p = 1, this represents a point in the line. When p
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationLecture 4.5 Schemes for Parabolic Type Equations
Lecture 4.5 Schemes for Parabolic Type Equations 1 Difference Schemes for Parabolic Equations One-dimensional problems: Consider the unsteady diffusion problem (parabolic in nature) in a thin wire governed
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationChapter 4 Euclid Space
Chapter 4 Euclid Space Inner Product Spaces Definition.. Let V be a real vector space over IR. A real inner product on V is a real valued function on V V, denoted by (, ), which satisfies () (x, y) = (y,
More informationEIGENVALUE PROBLEMS. EIGENVALUE PROBLEMS p. 1/4
EIGENVALUE PROBLEMS EIGENVALUE PROBLEMS p. 1/4 EIGENVALUE PROBLEMS p. 2/4 Eigenvalues and eigenvectors Let A C n n. Suppose Ax = λx, x 0, then x is a (right) eigenvector of A, corresponding to the eigenvalue
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationChapter 5. Formulation of FEM for Unsteady Problems
Chapter 5 Formulation of FEM for Unsteady Problems Two alternatives for formulating time dependent problems are called coupled space-time formulation and semi-discrete formulation. The first one treats
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 informationBackground. Background. C. T. Kelley NC State University tim C. T. Kelley Background NCSU, Spring / 58
Background C. T. Kelley NC State University tim kelley@ncsu.edu C. T. Kelley Background NCSU, Spring 2012 1 / 58 Notation vectors, matrices, norms l 1 : max col sum... spectral radius scaled integral norms
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differential Equations Steffen Börm Compiled July 12, 2018, 12:01 All rights reserved. Contents 1. Introduction 5 2. Finite difference methods 7 2.1. Potential equation.............................
More informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
More informationLecture 4 Eigenvalue problems
Lecture 4 Eigenvalue problems Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationClassical Mechanics. Luis Anchordoqui
1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body
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 informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
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 informationSolving Linear Systems of Equations
November 6, 2013 Introduction The type of problems that we have to solve are: Solve the system: A x = B, where a 11 a 1N a 12 a 2N A =.. a 1N a NN x = x 1 x 2. x N B = b 1 b 2. b N To find A 1 (inverse
More informationTMA4212 Numerical solution of partial differential equations with finite difference methods
TMA4 Numerical solution of partial differential equations with finite difference methods Brynjulf Owren January 3, 07 Translated and amended by E. Celledoni Preface The preparation of this note started
More informationLecture 1 INF-MAT3350/ : Some Tridiagonal Matrix Problems
Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems Tom Lyche University of Oslo Norway Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems p.1/33 Plan for the day 1. Notation
More informationBasics of Discretization Methods
Basics of Discretization Methods In the finite difference approach, the continuous problem domain is discretized, so that the dependent variables are considered to exist only at discrete points. Derivatives
More informationMatrix Algebra: Summary
May, 27 Appendix E Matrix Algebra: Summary ontents E. Vectors and Matrtices.......................... 2 E.. Notation.................................. 2 E..2 Special Types of Vectors.........................
More informationarxiv: v1 [math.na] 23 May 2013
ADI SPLITTING SCHEMES FOR A FOURTH-ORDER NONLINEAR PARTIAL DIFFERENTIAL EQUATION FROM IMAGE PROCESSING Luca Calatroni Cambridge Centre for Analysis, University of Cambridge Wilberforce Road, CB3 0WA, Cambridge,
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationWe wish to solve a system of N simultaneous linear algebraic equations for the N unknowns x 1, x 2,...,x N, that are expressed in the general form
Linear algebra This chapter discusses the solution of sets of linear algebraic equations and defines basic vector/matrix operations The focus is upon elimination methods such as Gaussian elimination, and
More informationMath 126 Final Exam Solutions
Math 126 Final Exam Solutions 1. (a) Give an example of a linear homogeneous PE, a linear inhomogeneous PE, and a nonlinear PE. [3 points] Solution. Poisson s equation u = f is linear homogeneous when
More informationFinite Difference Methods (FDMs) 2
Finite Difference Methods (FDMs) 2 Time- dependent PDEs A partial differential equation of the form (15.1) where A, B, and C are constants, is called quasilinear. There are three types of quasilinear equations:
More information3.5 Finite Differences and Fast Poisson Solvers
3.5 Finite Differences and Fast Poisson Solvers It is extremely unusual to use eigenvectors to solve a linear system KU = F. You need to know all the eigenvectors of K, and (much more than that) the eigenvector
More informationMath Bootcamp An p-dimensional vector is p numbers put together. Written as. x 1 x =. x p
Math Bootcamp 2012 1 Review of matrix algebra 1.1 Vectors and rules of operations An p-dimensional vector is p numbers put together. Written as x 1 x =. x p. When p = 1, this represents a point in the
More informationApril 26, Applied mathematics PhD candidate, physics MA UC Berkeley. Lecture 4/26/2013. Jed Duersch. Spd matrices. Cholesky decomposition
Applied mathematics PhD candidate, physics MA UC Berkeley April 26, 2013 UCB 1/19 Symmetric positive-definite I Definition A symmetric matrix A R n n is positive definite iff x T Ax > 0 holds x 0 R n.
More informationMathematical Foundations
Chapter 1 Mathematical Foundations 1.1 Big-O Notations In the description of algorithmic complexity, we often have to use the order notations, often in terms of big O and small o. Loosely speaking, for
More informationModule 3: BASICS OF CFD. Part A: Finite Difference Methods
Module 3: BASICS OF CFD Part A: Finite Difference Methods THE CFD APPROACH Assembling the governing equations Identifying flow domain and boundary conditions Geometrical discretization of flow domain Discretization
More informationAssignment on iterative solution methods and preconditioning
Division of Scientific Computing, Department of Information Technology, Uppsala University Numerical Linear Algebra October-November, 2018 Assignment on iterative solution methods and preconditioning 1.
More informationHigher-Order Difference and Higher-Order Splitting Methods for 2D Parabolic Problems with Mixed Derivatives
International Mathematical Forum, 2, 2007, no. 67, 3339-3350 Higher-Order Difference and Higher-Order Splitting Methods for 2D Parabolic Problems with Mixed Derivatives Jürgen Geiser Department of Mathematics
More informationAlgebra Workshops 10 and 11
Algebra Workshops 1 and 11 Suggestion: For Workshop 1 please do questions 2,3 and 14. For the other questions, it s best to wait till the material is covered in lectures. Bilinear and Quadratic Forms on
More informationAppendix A: Matrices
Appendix A: Matrices A matrix is a rectangular array of numbers Such arrays have rows and columns The numbers of rows and columns are referred to as the dimensions of a matrix A matrix with, say, 5 rows
More information88 CHAPTER 3. SYMMETRIES
88 CHAPTER 3 SYMMETRIES 31 Linear Algebra Start with a field F (this will be the field of scalars) Definition: A vector space over F is a set V with a vector addition and scalar multiplication ( scalars
More information