Numerical Solution to Parabolic PDE Using Implicit Finite Difference Approach
|
|
- Kathlyn Lyons
- 5 years ago
- Views:
Transcription
1 Numerical Solution to arabolic DE Using Implicit Finite Difference Approac Jon Amoa-Mensa, Francis Oene Boateng, Kwame Bonsu Department of Matematics and Statistics, Sunyani Tecnical University, Sunyani, Gana Department of Interdisciplinary Studies, University of Education, Winneba,. O. Box, Kumasi, Gana Abstract Tis paper examines an implicit Finite Difference approac for solving te parabolic partial differential equation (DE) in one dimension. We consider te Cran Nicolson sceme wic offers a better truncation error for bot time and spatial dimensions as compared wit te explicit Finite Difference metod. In addition te sceme is consistent and unconditionally stable. One downside of implicit metods is te relatively ig computational cost involved in te solution process, owever tis is compensated by te ig level of accuracy of te approximate solution and efficiency of te numerical sceme. A pysical problem modelled by te eat equation wit Neumann boundary condition is solved using te Cran Nicolson sceme. Comparing te numerical solution wit te analytical solution, we observe tat te relative error increases sarply at te rigt boundary, owever it diminises as te spatial step size approaces zero. Keywords: artial Differential Equation, Implicit Finite Difference, Cran Nicolson Sceme. Introduction Myriads of pysical penomenon lie diffusion, wave propagation, laser beam models, financial models, etc., are modelled by partial differential equations. Most often te equations under consideration are so complicated tat finding teir solutions analytically (e.g. by Laplace and Fourier transform metods, or in te form of a power series) is eiter impossible or impracticable (Süli, ). Tis may be due to several reasons including nonlinearity of problems, variable coefficients, difficulty in evaluating integrals analytically, approximation based on truncation of infinite series required, inappropriate solution space, etc. And, even problems tat could be solved analytically, we depend on numerical procedure to plot te solution. Consequently, numerical metods are employed as alternative approac for te approximation of te unnown analytical solution (Bruaset, and Tveito, 6). Te main idea of any numerical metod for a differential equation is to discretize te given continuous problem wit infinitely many degrees of freedom to obtain a discrete problem of system of equations wit only finitely many unnowns tat may be solved using numerical algoritms (Claes, 9). Altoug tere are several different numerical metods available for solving te parabolic equations, we focus on finite difference scemes. Tis is motivated by te fact tat tey are very simple to understand and tey are easy to generalize to more complex boundary value problems, and also are easily implemented on a computer. Te finite difference metod comes in two forms namely; te explicit metods and te implicit metods. Explicit metods generally are consistent, owever teir stability is restricted (LeVeque, ). On te oter and te implicit metods are consistent as well as unconditionally stable, owever tey are computationally costly compared to te explicit metods (Douglas and Kim, ). Tis is compensated by te ig level of accuracy of approximate solution and efficiency of te numerical sceme. In tis paper we use te Cran Nicolson metod wic is an implicit finite difference metod to solve a pysical problem modelled by te parabolic equation. We coose tis metod over te explicit metod because of its attractive properties mentioned.. Finite Difference Metod In tis section, we discuss te concept of te finite difference metod by providing useful terminologies, definitions and teorems. Te differencing metod solves a DE numerically by setting up a regular grid in space and time to compute approximate solutions at space or time points of tis grid (Strauss, 8). Te metod consist of approximating te differential operator by replacing te derivatives in te equation using differential 4
2 quotients. Tis process is commonly termed as discretization. Te error between te numerical solution and te exact solution is determine by te error tat is committed by going from a differential operator to a difference operator. Tis error is called te discretization error or truncation error. Definition.. Suppose te function u is C continuous in te neigbourood of x. For any > we ave u( x ) u( x) u( x) u ( x ) () were is a number between and : x ] x, x [. Truncating () at te first order gives u( x ) u( x) u( x) O( ) () were te term O ( ) indicates tat error of te approximation is proportional to We refer to equation () as te forward difference approximation of u. We deduce from equation () tat tere exists a constant C suc tat, for > sufficiently small, we ave. were u( x ) u( x) u( x) C, C sup [ x, x ] u( ) for, given. Te approximation of u at te point x is said to be consistent at te first order. Te order of a first derivative can be defined in a more general sense as follows; Definition.. Te approximation of te derivative u at point x is of order p (p > ) if tere exists a constant C p >, independent of, suc tat te error between te derivative and its approximation is bounded by C, i.e., p is exactly O ( ). Similarly, te first order bacward difference approximation of u at point x can be defined as u x u x u x O ( ) ( ) ( ) ( ) (3) In order to improve te accuracy of te approximation of u we define a consistent approximation, called te central difference approximation, also referred to as symmetric quotient, by combining te forward difference and te bacward difference (Lapidus, and inder,999). Tat is, taing te points x and x into consideration. Definition.3. Suppose tat te function uc 3 ( ) were ] x, x [, ten we write, 3 u( x ) u( x) u( x) u ( x) u ( ) (6) 6 3 u( x ) u( x) u( x) u ( x) u ( ) () 6 were ] x, x [ and ] x, x[. Subtracting () from (6) to we obtain, 3 u( x ) u( x ) u( x) u ( ) (8) 6 were ] x, x [. It follows from te definition tat, for every ], [, we ave a bound on te approximation error as: 5
3 u( x ) u( x ) u x ( ) C, C sup [ x, x] u( ) 6 If u is C continuous, ten te approximation is consistent at te order one only. Te order of te approximation is related to te regularity of te function u (Alan, ). Te approximation of te second derivative and its associated truncation error is given by te following lemma. Lemma.. Suppose u is a C 4 continuous function on an interval[ x, x ],. Ten, tere exists a constant C > suc tat for every ], [, we ave u( x ) u( x) u( x ) u( x) C. (9) Te differential quotient ( u( x ) u( x) u( x )) is consistent second-order approximation of te second derivative u of u at point x. roof. We begin by using te Taylor expansion up to te fourt order to obtain te following: u x u x u x u x u x (4) ( ) ( ) ( ) ( ) ( ) u ( ) () u x u x u x u x u u 6 4 were ] x, x [ and ] x, x[. By te intermediate value teorem, we can write: 3 4 (4) ( ) ( ) ( ) ( ) ( ) ( ) () u( x ) u( x) u( x ) u x u for ] x, x [. Hence, from equation (9) we deduce te constant (4) u ( ) C sup. [ x, x] (4) ( ) ( ), Remar.. Te error estimates (9) depends on te regularity of te function u. If u is C 3 continuous, ten te error is of order only. 3. arabolic roblem in -dimension arabolic DE describes time dependent pysical penomenon and is mostly used in modeling diffusion problems suc as eat transfer (Morton and Mayer, 5). In tis section, we consider te one-dimensional omogenous eat equation posed in te bounded domain [, L] : u u, for ( xt, ) t t u( x,) u( x), for x, () BC : u(, t), u( L, t), for, were is a given scalar value and u C (, ). Te boundary condition (BC) specified in () is of te Diriclet boundary condition. Te boundary condition can be replaced by Neumann boundary condition given as, u u (, t), ( L, t),, x x 3.. Analytic Solution Te general analytic solution of te eat equation () wit omogeneous Diriclet boundary condition is as follows (Kersale, 3): 6
4 n nx t L n n L u( x, t) b sin e were L n bn u ( x)sin dx, n,,3, L L Tis solution can be obtained by using te metod of separation of variables and Fourier transform. Wen te boundary conditions are of te Neumann ind ten te solution is written as: were n t L nx n n L u( x, t) a a e cos L n an u ( x)cos dx, n,,, L L 3.. Implicit Finite Difference Sceme To begin wit, we introduce equi-distributed grid points ( x ) given by i i N xi i, wit L ( N ) to discretize te domain, in space. Furtermore, we consider time step, to discretize te time domain, ence defining a regular grid, t,. We denote by U( x, i t ) te value of te numerical solution at point ( x, i t ) and u( x, t ) te exact solution of te equation (). Te initial data u is discretized by: Ui u ( xi) for i {,, n }. Te boundary conditions are also andled in similar vein. Diriclet boundary conditions is represented in a discrete form in a straigt forward manner as: U, UN. In te case of Neumann boundary conditions, in order to improve te accuracy of te solution we apply te central difference approximation to discretize te derivatives involved. Tis is given by N N U U, U U, Using te finite difference approximations for te derivative in te eat equation we obtain te following: u u( x, t ) u( x, t) ( x, t) O( ) (3) t u u( x, t) u( x, t) u( x, t) ( x, t) O( ) (4) x Considering equations (3) and (4), te finite difference scemes can be written in a more general form as a weigted average of te explicit and implicit scemes, producing te teta sceme, also written as te -sceme : U i Ui Ui Ui ui Ui Ui Ui ( ) (5) We simplify equation (5) furter as: Ui Ui U i Ui U i ( ) U i Ui u i (6) were denotes a positive parameter, [,] and Te forward and bacward Euler metods, and te Cran Nicolson metod are considered to be special cases of te -sceme defined by (6), wit for te forward Euler, for bacward Euler, and for te Cran Nicolson metod. Te -sceme is an implicit sceme for.
5 Te difference equation of te -sceme (6) can be written in a vector form as: were I is a unit matrix, and I A U I ( ) A U () T,,, N U U U U... A Tus, to find U, we need to solve a tridiagonal linear system, wic can be done by employing Tomas algoritm. Lemma 3.. Let T be an N N tridiagonal matrix of te form: b c.. a b c. T a b c.. a b Te eigenvalues and te corresponding eigenvectors of T are: b a sin c N a sin ac cos, v c N, N N a N sin c N,, N (Kicaid and Ceney, ) Lemma 3.. Consider te -sceme (4), for [,]. (i) If, te -sceme is unconditionally stable in L - norm (and convergent) (ii) If, te -sceme is stable in L - norm (and convergent of tis norm) under te condition:. ( ) roof. To establis te stability of te -sceme, we represent te sceme () in a matrix form as U U (8) and we determine weter te magnitude of any eigenvalue of coefficient matrix exceeds. From equation (), we can write, It follows from (9) tat, U I A I ( ) A U (9) I A I ( ) A 8
6 Since te eigenvalues of I, are easily found be ( ) A and I A ( ) are te same, te eigenvalues appearing in te matrix were are te eigenvalues of A. Te values of can directly be deduced from lemma 3.. as, cos 4sin, N N,, N () For stability of te sceme, it is necessary tat, We denote ( ) N,,, N () () Substituting equations () and () into () gives, 4 ( ) sin 4 sin. (3) Since, by assumption and by definition, ten te rigt and side of (3) is positive, and terefore te inequality can be written as ( 4 sin ) 4 ( ) sin 4 sin (4) Te rigt part of te inequality is automatically satisfied for all. Te left part is satisfied wen Te strongest restriction on (and ence on te step size in time) occurs wen 4 ( ) sin (5) sin value, tat is. In tat case, (5) yields (6) assumes its largest Considering (6), we note tat te -sceme is unconditionally stable if, for any arbitrary. If, ten te -sceme is stable provided tat, 3... Cran Nicolson Sceme We consider Cran Nicolson metod for solving te eat equation in tis paper due to te fact tat te metod as second order accuracy in bot time and spatial dimensions, and is unconditionally stable (Arivis, et. al., 6). Tis means tat te stability of te sceme is not affected by any step size cosen for te time dimension as seen in lemma 3.. Te metod as mentioned earlier can be extracted from te -sceme by setting. We can write te sceme as follows: U i Ui Ui U i ui Ui U i Ui Rearranging () gives, () i i i i Ui U U i U Ui U Ui U (8) 9
7 were Now, before we generate te system of linear equations out of (8), we recall te discrete form of te boundary conditions discussed earlier. Let us consider for instance te discrete omogeneous Neumann boundary condition given by, N N U U U U For i,, N, we generate a system of linear equations out of equation (8), taing into consideration te boundary conditions. We obtain te following system: U U U U U U U U U U U U U U U U U U U U U U U U U U 3 3 Te above system can be written in a matrix form as U U U U U ( ) U U ( U U ) ( ) U U ( U U U ) ( ) U U ( U 3 U U ) (9) ( ) U U9 ( U U9 ) Eac time step, requires tat we solve te tridiagonal system (9). Tis maes it computationally costly compared to te explicit metods. However, for one dimensional eat problem, tere is not muc difference between te explicit metod and te implicit metod in terms of computational cost. 4. Test Case Consider te equation describing one-dimensional, single-pase, sligtly compressible flow in a producing petroleum reservoir wic is given, for x 5 and t, by;, px, t C p x t K t x, (3) te porous medium and te reservoir are omogeneous, tat te liquid is ideal, and tat gravitational effects are negligible. Te distance is represented by x, time by t, p te pressure, te dimensionless constant porosity of te medium, te viscosity te permeability of te medium and C te compressibility (in {pounds per square inc} - ). Assuming tat C K.4 days/ft, and tat te following conditions old: p x,.5, x 5,, 5, p t p t, t. Kx x 8
8 Determine te pressure p at t = 5, using te Cran Nicolson metod wit t.5 and x 5 (Burden et. al., 98). Te analytical solution for te eat problem described above is given by n.4 5 t 5 nxn p( xn, t).5 sin ne cos n n 5. Te pressure p, at time t = 5 given at te spatial step, x 5 is provided as follows x px (,5) x x x x x x x x x x x 4.. Numerical Solution We solve te eat problem (3) numerically using te Cran-Nicolson metod wit t.5 and x 5. We represent te approximate solution of te eat problem (3) at te next time step by ( x, t ), wic for simplicity can be written as. Expressing (3) by te Cran Nicolson sceme (8), we i obtain i i i i i i i i i (3) To incorporate te Neumann boundary condition,, 5, t t, Kx x we discretize te boundary data using central difference to obtain, i i. Te spatial domain is discretized using x = i. At te left boundary i = and so x =. Te boundary reduces to. We andle te rigt boundary in similar manner by setting i =, wic implies x = 5. Te boundary can ten be written as 9 Now, we generate a system of equations wic will enable us solve for te interior nodes. For i = in equation (3) we ave, But and 8
9 For i For i = 3 3 For i = 3, 4., we get a system of equations wic can be written in a matrix form as ( ) ( ) ( ) ( ) ( ) ( 3 ) ( ) 9 ( 9 ).5 and.5. Following from te eat problem (3), we let Te initial conditions, x,.5 can also be written in a discrete form as follows: i.5.5 (3).5 Substituting te value of and te initial condition into (3), we obtain a tridiagonal system: (33) Hence te approximate solution as a result of solving te tridiagonal system (33), for =,,,... 9, is given as follows: 8
10 xi i x (,5) i Te solution of bot te analytical approac and Cran Nicolson sceme for spatial step, = 5 and = are compared respectively in figure. Also te values of te relative error of bot solutions are plotted in figure. Figure : Comparison of exact solution and approximate solution for = 5 and =. Figure : Relative error of te approximate solution for = 5 and =. Clearly, we see from figure tat, te Cran Nicolson metod as good approximate solution. However, te relative error as a sarp increase at te rigt boundary as sown in figure. We observe tat te error associated wit te numerical sceme diminises as te spatial step size, approaces zero. 83
11 5. Conclusion An implicit Finite Difference approac for solving te parabolic partial differential equation in one dimension as been discussed. We considered te Cran Nicolson sceme due its attractive properties suc as te second order accuracy in bot time and spatial dimensions and stability wic is unconditional. Tese mae it preferable to te explicit metod wic as first order accurate in time dimension and as a restricted stability. We note tat, toug implicit metod are computationally costly due to tat fact tat at eac time step, a wole system of equations is solved, tis is compensated by te accuracy and te unrestricted stability. Te Test case as sown tat te Cran Nicolson metod is a good alternative to solving parabolic DEs, were analytical solution is difficult to obtain or even impracticable to use. To furter improve on te accuracy of te numerical solution to te parabolic DE, we sall explore te avenues of wavelet approaces. References. Arivis, G., Maridais, C., & Nocetto, R. (6). A posteriori error estimates for te Cran Nicolson metod for parabolic equations. Matematics of computation, 5(54), Alan, J. () Advanced Engineering Matematics, st Edition. Harcourt Academic ress, Boston 3. Burden. L. R., Faires, J. D. and Reynolds. C. A. (98). Numerical Analysis. rindle, Weber and Scmidt, page Bruaset, M. A. and Tveito, A. (6). Numerical solution of partial differential equation on parallel computers, Springer-Verlog Berlin Heidelberg. 5. Claes, J. (9). Numerical Solution of artial Differential Equations by te Finite Element Metod. Cambridge University ress, New Yor. 6. Douglas Jr, J., & Kim, S. (). Improved accuracy for locally one-dimensional metods for parabolic equations. Matematical Models and Metods in Applied Sciences, (9), Lapidus, L. and inder, F. G. (999). Numerical solution of artial Differential Equation in Science and Engineering, University of Vermont. 8. LeVeque, J. R. (). Finite Difference Metods for Ordinary and artial Differential Equations: Steady- State and Time-Dependent roblems. SIAM, iladelpia. 9. Kersale E. (3). Analytic Solutions of artial Differential Equations. Scool of Matematics, University of Leeds. ersale.. Kincaid, D. and Ceney, W. (). Numerical Analysis: Matematics of Scientific Computing, 3rd Edition. Broos Cole. Morton, K. W. and Mayer, D. (5). Numerical Solution of artial Differential Equation Second Edition, Cambridge University ress.. Strauss, A. W. (8). artial Differential Equations: An Introduction. Brown University, Jon Wiley & Sons, Florida. 3. Süli, E. (). Finite Element Metods for artial Differential Equations. Lecture Notes. Matematical Institute University of Oxford. 84
Chapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationLecture 21. Numerical differentiation. f ( x+h) f ( x) h h
Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationNUMERICAL DIFFERENTIATION
NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
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 information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
More informationMath 161 (33) - Final exam
Name: Id #: Mat 161 (33) - Final exam Fall Quarter 2015 Wednesday December 9, 2015-10:30am to 12:30am Instructions: Prob. Points Score possible 1 25 2 25 3 25 4 25 TOTAL 75 (BEST 3) Read eac problem carefully.
More informationA Finite Element Primer
A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationarxiv: v1 [math.na] 7 Mar 2019
Local Fourier analysis for mixed finite-element metods for te Stokes equations Yunui He a,, Scott P. MacLaclan a a Department of Matematics and Statistics, Memorial University of Newfoundland, St. Jon
More informationResearch Article Cubic Spline Iterative Method for Poisson s Equation in Cylindrical Polar Coordinates
International Scolarly Researc Network ISRN Matematical Pysics Volume 202, Article ID 2456, pages doi:0.5402/202/2456 Researc Article Cubic Spline Iterative Metod for Poisson s Equation in Cylindrical
More informationInf sup testing of upwind methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 000; 48:745 760 Inf sup testing of upwind metods Klaus-Jurgen Bate 1; ;, Dena Hendriana 1, Franco Brezzi and Giancarlo
More informationNON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE
Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract.
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationSection 3.1: Derivatives of Polynomials and Exponential Functions
Section 3.1: Derivatives of Polynomials and Exponential Functions In previous sections we developed te concept of te derivative and derivative function. Te only issue wit our definition owever is tat it
More informationAMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationInfluence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations
International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 2, pp. 281 302 (2017) ttp://campus.mst.edu/ijde Influence of te Stepsize on Hyers Ulam Stability of First-Order Homogeneous
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More informationNumerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions
Proc. of te 8t WSEAS Int. Conf. on Matematical Metods and Computational Tecniques in Electrical Engineering, Bucarest, October 16-17, 2006 219 Numerical Solution of One Dimensional Nonlinear Longitudinal
More informationQuantum Theory of the Atomic Nucleus
G. Gamow, ZP, 51, 204 1928 Quantum Teory of te tomic Nucleus G. Gamow (Received 1928) It as often been suggested tat non Coulomb attractive forces play a very important role inside atomic nuclei. We can
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationEfficient, unconditionally stable, and optimally accurate FE algorithms for approximate deconvolution models of fluid flow
Efficient, unconditionally stable, and optimally accurate FE algoritms for approximate deconvolution models of fluid flow Leo G. Rebolz Abstract Tis paper addresses an open question of ow to devise numerical
More informationOn the absence of marginal pinching in thin free films
On te absence of marginal pincing in tin free films P. D. Howell and H. A. Stone 6 August 00 Abstract Tis paper concerns te drainage of a tin liquid lamella into a Plateau border. Many models for draining
More informationDigital Filter Structures
Digital Filter Structures Te convolution sum description of an LTI discrete-time system can, in principle, be used to implement te system For an IIR finite-dimensional system tis approac is not practical
More informationUniversity Mathematics 2
University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at
More informationAnalytic Functions. Differentiable Functions of a Complex Variable
Analytic Functions Differentiable Functions of a Complex Variable In tis capter, we sall generalize te ideas for polynomials power series of a complex variable we developed in te previous capter to general
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationlecture 35: Linear Multistep Mehods: Truncation Error
88 lecture 5: Linear Multistep Meods: Truncation Error 5.5 Linear ultistep etods One-step etods construct an approxiate solution x k+ x(t k+ ) using only one previous approxiation, x k. Tis approac enoys
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips
More informationModule 4 Linear Boundary value problems
Module 4 Linear Boundary value problems Lecture Content Hours 1 Finite Difference Metods: Diriclet type boundary 1 condition Finite Difference Metods: Mixed boundary condition 1 3 Sooting Metod 1 4 Sooting
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationRecall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if
Computational Aspects of its. Keeping te simple simple. Recall by elementary functions we mean :Polynomials (including linear and quadratic equations) Eponentials Logaritms Trig Functions Rational Functions
More informationReflection Symmetries of q-bernoulli Polynomials
Journal of Nonlinear Matematical Pysics Volume 1, Supplement 1 005, 41 4 Birtday Issue Reflection Symmetries of q-bernoulli Polynomials Boris A KUPERSHMIDT Te University of Tennessee Space Institute Tullaoma,
More informationarxiv: v3 [math.na] 31 May 2016
Stability analysis of pressure correction scemes for te Navier-Stoes equations wit traction boundary conditions Sangyun Lee a, Abner J. Salgado b a Center for Subsurface Modeling, Institute for Computational
More informationContinuity. Example 1
Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)
More information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
More informationApplication of numerical integration methods to continuously variable transmission dynamics models a
ttps://doi.org/10.1051/ssconf/2018440005 Application of numerical integration metods to continuously variable transmission dynamics models a Orlov Stepan 1 1 Peter te Great St. Petersburg Polytecnic University
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationLecture 10: Carnot theorem
ecture 0: Carnot teorem Feb 7, 005 Equivalence of Kelvin and Clausius formulations ast time we learned tat te Second aw can be formulated in two ways. e Kelvin formulation: No process is possible wose
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationarxiv: v1 [math.na] 20 Jul 2009
STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
More informationSome Review Problems for First Midterm Mathematics 1300, Calculus 1
Some Review Problems for First Midterm Matematics 00, Calculus. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd,
More information1 Lecture 13: The derivative as a function.
1 Lecture 13: Te erivative as a function. 1.1 Outline Definition of te erivative as a function. efinitions of ifferentiability. Power rule, erivative te exponential function Derivative of a sum an a multiple
More informationComputer Derivations of Numerical Differentiation Formulae. Int. J. of Math. Education in Sci. and Tech., V 34, No 2 (March-April 2003), pp
Computer Derivations o Numerical Dierentiation Formulae By Jon H. Matews Department o Matematics Caliornia State University Fullerton USA Int. J. o Mat. Education in Sci. and Tec. V No (Marc-April ) pp.8-87.
More informationError estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs
Interfaces and Free Boundaries 2, 2000 34 359 Error estimates for a semi-implicit fully discrete finite element sceme for te mean curvature flow of graps KLAUS DECKELNICK Scool of Matematical Sciences,
More informationDIRECTLY SOLVING SECOND ORDER LINEAR BOUNDARY VALUE PROBLEMS OF ORDINARY DIFFERENTIAL EQUATIONS. Ra ft Abdelrahim 1, Z. Omar 2
International Journal of Pure and Applied Matematics Volume 6 No. 6, -9 ISSN: - (printed version); ISSN: -95 (on-line version) url: ttp://www.ijpam.eu doi:.7/ijpam.v6i. PAijpam.eu DIRECTLY SOLVING SECOND
More informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
More information