Chapter 5 FINITE DIFFERENCE METHOD (FDM)
|
|
- Moses Young
- 6 years ago
- Views:
Transcription
1 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 equations by finite difference equations. Tese finite difference approximations are algebraic in form; tey relate te value of te dependent variable at a point in te solution region to te values at some neigboring points. Tus a finite difference solution basically involves tree steps: Dividing te solution region into a grid of nodes. Approximating te given differential equation by finite difference equivalent tat relates te dependent variable at a point in te solution region to its values at te neigboring points. Solving te difference equations subject to te prescribed boundary conditions and/or initial conditions. Te course of action taken in tree steps is dictated by te nature of te problem being solved, te solution region, and te boundary conditions. Te most commonly used grid patterns for two-dimensional problems are sown in Fig. 5.. (a) (b) (c) (d) Fig. 5. Common grid patterns (a) Rectangular grid, (b) skew grid, (c) triangular grid, and (d) circular grid. Capter 5 Page
2 MEE7 Computer Modeling Tecniques in Engineering 5. Finite Element Scemes Before finding te finite difference solutions to specific PDEs, we will look at ow one constructs finite difference approximations from a given differential equation. Tis essentially involves estimating derivatives numerically. Consider a function f(x) sown in Fig.5., we can approximate its derivative, slope or te tangent at P by te slop of te arc PB, given te forward-difference formula, x x x x + x x Fig.5. Estimates for te derivative of f(x) at P using forward, backward, and central differences. f '( x ) f ( x + x) f ( x x ) (5.) or te slop of te arc AP, yielding te backward-difference formula, f ( x) f ( x x) f '( x ) x or te slop of arc AB, resulting in te central-difference formula, f ( x + x) f ( x x) f '( x ) x We can also estimate te second derivative of f(x) at P as (5.) (5.) f ''( x ) f ( x + x) f ( x ) + f ( x x) ( x) (5.) Any approximation of a derivative in terms of values at a discrete set of points is called finite difference approximation. A more general approac is using Taylor s series. According to te well-known expansion, Capter 5 Page
3 MEE7 Computer Modeling Tecniques in Engineering f ( x x) f ( x ) xf '( x ) x f x x f x! ( ) ''( )! ( ) '''( ) L (5.5) and f ( x x) f ( x ) xf '( x ) x f x x f x! ( ) ''( )! ( ) '''( ) + + L (5.6) upon adding tese expansions, f ( x + x) + f ( x x) f ( x ) + ( x) f ''( x ) + O( x) (5.7) were O ( x) is te error introduced by truncating te series. Assuming tat tese terms are negligible, we can obtain f ''( x ) f ( x + x) f ( x) + f ( x x) x wic is Eq. (5.). Subtracting Eq. (5.6) from Eq. (5.5) and neglecting terms of te order ( x) yields f '( x ) f ( x + x) f ( x x) x wic is Eq.(5.). Tis sows tat te leading errors in Eqs. (5.) and (5.) are of te order ( x). Higer order finite difference approximations can be obtained by taking more terms in Taylor series expansion. Example 5.. To apply te difference metod to find te solution of a function (x,t), we divide te solution region in te x-t plane into equal rectangles or meses of sides x and t as in Fig.5.. We let te coordinates (x,t) of a typical grid point or node be xi x, i,,,... tj t, j,,,... and te value if at P be P ( i x, j t) ( i, j) (5.) Wit tis notation, te central difference approximations of te derivatives of at te (i,j)t node are Capter 5 Page
4 MEE7 Computer Modeling Tecniques in Engineering x / i, j ( i+, j) ( i, j), x (5.9a) t / i, j (, i j+ ) (, i j), t (5.9b) xx tt / / i, j i, j ( i+, j) ( i, j) + ( i, j), (5.9c) ( x) (, i j+ ) (, i j) + (, i j ), (5.9d) ( t) Fig.5. FD mes for two independent variables x and t. Example 5. Laplace s equation in two dimensions may ten be expressed as x y (5.) Eq. (5.) is known as te five-point equal arm difference equation using Fig. 5.. Te FD solution procedure of Poisson s or Laplace s equations may ten be summarized as follows: Divide te domain of interest (in wic te potential is to be determined) into suitable fine grid. Instead of a solution for (x,y), wic provides its continuous variation for a given carge distribution ρ(x,y), te FD solution will provide discrete values of at te nodes of te establised grid. Capter 5 Page
5 MEE7 Computer Modeling Tecniques in Engineering Apply te difference equation at eac node of te grid to obtain, for example, N equations in te N unknown node potentials. Solve te resulting system of equations, eiter iteratively or using one of te direct metods. Fig. 5. Geometry of te five-point star used in -D difference equations. 5. Difference Equation at Interface between Two Dielectric Media In many engineering applications, interfaces between two different dielectric media are encountered. For tis, we will derive a special case difference equation tat sould be satisfied at nodes on te interface between two dielectrics. Fig. 5.5 illustrates te geometry of an interface, and te difference equation in tis case may be obtained from Gauss s law for te electric field, Fig. 5.5 Geometry of grid nodes at te interface between medium of ε and medium of ε. Capter 5 Page 5
6 MEE7 Computer Modeling Tecniques in Engineering E ds q s ε (5.) q in eq. (5.) because tere is no free carge enclosed by te surface s. Substituting E, we obtain ε ε ε ds dc n dc s c c (5.) Te surface integration on te left-and side of equation (5.) was replaced by a contour integration, because in Fig. 5.5 we are dealing wit a two-dimensional case, and te solution of is independent of te axial independent variable z. / n denotes te normal derivative of on te contour c. Carrying out detailed integration of / n along c, we obtain ( εdc) ε + ε + ( ε) n c + ε + ε + ( ε) (5.) Rearranging te terms in equation (5.), we obtain ε + ε + ( ε + ε ) + ( ε + ε ) ( ε + ε ) (5.) Te equation (5.) may be used at interfaces between two dielectrics. Its use in engineering problems will be illustrated by following examples. 5. FDM for PDEs 5.. FD of Diffusion PDEs k (5.5) t x were k is a constant. Te equivalent FD approximation is ( i, j+ ) ( i, j) ( i+, j) ( i, j) + ( i, j) k, (5.6) t ( x) were x i x, i,,,..., n, t j t, j,,,... In Eq. (5.6), we ave used te forward difference formula for te derivative wit respect to t and central difference formula for tat wit respect to x. If we let Capter 5 Page 6
7 MEE7 Computer Modeling Tecniques in Engineering r t k( x) (5.7) ten Eq. (5.6) can be written as (, i j+ ) r( i+, j) + ( r) (, i j) + r( i, j) (5.) Tis explicit formula can be used to compute ( xt, + t) explicitly in terms of (x,t). Tus te values of along te first time row (see Fig. 5.), t t, can be calculated in terms of te boundary and initial conditions, ten te values of along te second time row, t t, are calculated in terms of te first time row, and so on. In order to ensure a stable solution or reduce errors, care must be exercised in selecting te value of r in Eq. (5.7) and Eq.(5.) is valid only if te coefficient (-r) in Eq.(5.) is nonnegative or o < r /. If we coose r /, Eq. (5.) becomes (, i j+ ) [ ( i+, j) + ( i, j)] (5.9) 5.. FD of Helmoltz Equation + k (5.) were E z for TM modes or H z for TE modes in waveguide problems, wile k is te wave number given by k ω µε β. To apply te FD metod, we discretize te cross section of te wave guide by a suitable mes. Applying central difference approximation for partial derivatives to Eq. (5.) gives ( i +, j) + ( i, j) + ( i, j + ) + ( i, j ) ( k ) ( i, j) (5.) 5.. FD of Wave Equation u (5.) x t were u is te speed of te wave. An equivalent FD formula is u ( i+, j) ( i, j) + ( i, j) (, i j+ ) (, i j) + (, i j), ( x) ( t) were x i x, i,,,..., n, t j t, j,,,... Tis equation can be written as Capter 5 Page 7
8 MEE7 Computer Modeling Tecniques in Engineering ( i, j + ) ( r) ( i, j) + r[ ( i +, j) + ( i, j)] ( i, j ) (5.) were (i,j+) is an approximation to (x,t) and r is te aspect ratio given by u t r x (5.) Equation (5.) is an explicit formula for te wave equation. For te solution algoritm in Eq.(5.) to be stable, te aspect ratio r < will be cosen. If we coose r, Eq. (5.) becomes (, i j+ ) ( i+, j) + ( i, j) (, i j ) (5.5) Unlike te single-step scemes of Eqs.(5.), te two-step scemes of Eqs. (5.) and (5.5) require tat te values of at times j and j- be known to get at time j+. Tus, we must derive a separate algoritm to start te solution of Eq. (5.) or (5.5); tat is, we must compute (i,) and (i,). To do tis, we utilize te prescribed initial condition. For example, suppose te initial condition on te FDE in Eq. (5.) is t t We use te backward-difference formula ( x, ) ( i, ) ( i, ) t t or (,) i (, i ) (5.6) Substituting Eq. (5.6) into Eq. (5.) and taking j (i.e.,at t ), we can obtain r (,) i ( r) (,) i + [ ( i,) + ( i+,) (5.7) Using te starting formula Eq. (5.7) togeter wit te prescribed boundary and initial conditions, te value of at any grid point (i,j) can be obtained directly from Eq. (5.). Accuracy and Stability of FD Solutions Tere are tree sources of errors tat are nearly unavoidable in numerical solution of pysical problems: () modeling errors, () truncation (or discretization) errors, () round-off errors. Capter 5 Page
9 MEE7 Computer Modeling Tecniques in Engineering Eac of tese error types will affect accuracy and terefore degrade te solution. Fig. 5.6 sows us te error as a function of te mes size. Fig. 5.6 Error as a function of te mes size 5.5 Example of D FDM Consider te rectangular region sown in Fig. 5.7a. Te electric potential is specified on te conducting boundaries. Use te finite difference representation to solve for te potential distribution witin tis region. Fig. 5.7a Rectangular geometry and boundary condition for te electric potential problem. Capter 5 Page 9
10 MEE7 Computer Modeling Tecniques in Engineering Fig.5.7b Geometry of x finite difference mes. Solution Te electric potential everywere in te rectangular region sould satisfy Laplace s equation. Using a numerical solution means we will define in te rectangular region of interest by calculating its values at discrete points, te nodes of a mes. Te step-by-step solution procedure includes te following: () Layout a coarse square mes and identify te nodes at wic te electric potential is t be calculated. Te geometry of a mes is sown in Fig. 5.7b. Te value of (mes size) in tis case is 5cm. () Replace Laplace s equation by its finite difference representation. i+, j + i, j + i, j+ + i, j i, j (5.) i,j are te discrete values of te potential at points (nodes) witin te domain of interest. () Apply te difference equation in step at eac node. At node, or (. 5) (5.9) At node, (. 5) Capter 5 Page
11 MEE7 Computer Modeling Tecniques in Engineering or (5.) + At node, 5 + (. ) or (5.) () Eqs. (5.9) to (5.) are tree equations in te tree unknowns,, and. Tese tree equations may be solved using one of te metods described in calculus courses. Te results are.79, 7., 6.79 (5) Wit te coarse mes we used, we do not expect to get accurate final results. Redoing te problem wit a smaller value of sould improve te accuracy of te solution. Fig.5.7c sows te mes geometry for.5cm, wic is alf te mes size used in te previous calculations. In tis case, owever, we ave twenty-one unknown values of te potential at te various nodes. (6) Once again, applying te difference equation of step at te various nodes results in te following x matrix: Capter 5 Page
12 MEE7 Computer Modeling Tecniques in Engineering Fig.5.7c Geometry of.5-cm mes wit twenty-one potential nodes. Table 5. Comparison between Finite Difference and Analytical Results Potential values Percentage error Potential values Percentage error Analytical solution ( 5 cm ) (.5 cm ) Table 5. compares te results for 5cm, and.5cm. from tis comparison, it is clear tat te.5cm results agree better wit te analytical solution available for tis simple geometry. As expected, te accuracy of te FD results improves wit te reduction in size. Clearly, any furter reduction in results in a larger-size matrix; ence, a compromise sould be made between te desired accuracy and te computational time and effort required. Te solution for te electric potential at te various nodes is given by Capter 5 Page
13 MEE7 Computer Modeling Tecniques in Engineering EXAMPLE In te 6 m rectangular region sown in Figure.(a), te electric potential is zero on te boundaries. Te carge distribution, owever, is uniform and given by p v ε. Solve Poisson s equation to determine te potential distribution in te rectangular region. Solution To determine te potential distribution in te rectangular region, we use Poisson s equation. + x y ρ ε v wit zero potential on te boundaries. Figure. Geometry of te 6 m rectangular region and te m mes. By establising te rectangular grid sown in Figure.(b), we realize tat we ave six nodes and, ence, six unknown potentials for wic to solve. Replacing by its finite difference representation, we obtain (, +, +, +,, ) + i+ j i j i j+ i j i j Capter 5 Page
14 MEE7 Computer Modeling Tecniques in Engineering It sould be noted tat altoug te altoug te mes size was not explicitly used in solving Laplace s equation in te previous example, is included as a part of te matrix formation in solving Poisson s equation. In SI system of units, sould be in meters. By applying te preceding difference equation at te various nodes in Figure.(b), we obtain te following matrix equation: 6 5 Instead of solving te resulting six equations, we may note some symmetry considerations in Figure.(b). It is clear tat 6 5 and tat Taking tese symmetry considerations into account, te number of equations reduces to two, and we obtain te following solution:.56, 5.7 To improve te accuracy of te potential distribution, finer mes suc as te one sown in Figure. is required. Because of te large number of nodes in tis case, symmetry sould be used, and a solution for only one-quarter of te rectangular geometry is desired. Te application of te difference equation at nodes,,, and 5 sould proceed routinely, wereas special care sould be exercised at te boundary nodes, 6, 9,, and, and also at te corner node. For example, applying te difference equation at node 6 yields ) ( a Or ) ( Capter 5 Page
15 MEE7 Computer Modeling Tecniques in Engineering In equation + ( ), symmetry was used to x y complete te five-point star difference equation. Specifically te potential at node a to te rigt of 6 was taken equal to. Similarly at te corner node, we obtain 5 ( ) + b c Because of symmetry, b and 9 c, ence, ( + 9 ) + Fig.. Te finer mes solution and symmetry consideration of example. Te matrix equation for te twelve nodes sown in Figure. is ten Capter 5 Page 5
16 MEE7 Computer Modeling Tecniques in Engineering Te coefficient in te coefficient matrix (to te left) of equation appears wenever symmetry consideration is used at boundary and corner nodes. It sould be noted tat te coefficient matrix in equation is te same for bot Laplace s and Poissons equations. Te constant vector on te rigt-and side of equation, owever, depends on te carge distribution witin and te potential at te boundaries of te region of interest. Furtermore, if instead of a uniform carge distribution we ave a given carge distribution ), ( y x v ρ, te constants vector on te rigt-and side of equation sould reflect te value of ), ( y x v ρ calculated at eac node. Solution of equation gives , Capter 5 Page 6
NUMERICAL 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 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 informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
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 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 informationExercise 19 - OLD EXAM, FDTD
Exercise 19 - OLD EXAM, FDTD A 1D wave propagation may be considered by te coupled differential equations u x + a v t v x + b u t a) 2 points: Derive te decoupled differential equation and give c in terms
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 informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
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 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 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 informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
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 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 informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
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 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 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 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 informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More informationQuantum Numbers and Rules
OpenStax-CNX module: m42614 1 Quantum Numbers and Rules OpenStax College Tis work is produced by OpenStax-CNX and licensed under te Creative Commons Attribution License 3.0 Abstract Dene quantum number.
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 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 informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationGRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES
GRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES Cristoper J. Roy Sandia National Laboratories* P. O. Box 5800, MS 085 Albuquerque, NM 8785-085 AIAA Paper 00-606 Abstract New developments
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 informationCalculus I - Spring 2014
NAME: Calculus I - Spring 04 Midterm Exam I, Marc 5, 04 In all non-multiple coice problems you are required to sow all your work and provide te necessary explanations everywere to get full credit. In all
More informationSection 15.6 Directional Derivatives and the Gradient Vector
Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,
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 informationFall 2014 MAT 375 Numerical Methods. Numerical Differentiation (Chapter 9)
Fall 2014 MAT 375 Numerical Metods (Capter 9) Idea: Definition of te derivative at x Obviuos approximation: f (x) = lim 0 f (x + ) f (x) f (x) f (x + ) f (x) forward-difference formula? ow good is tis
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More information7.1 Using Antiderivatives to find Area
7.1 Using Antiderivatives to find Area Introduction finding te area under te grap of a nonnegative, continuous function f In tis section a formula is obtained for finding te area of te region bounded between
More informationb 1 A = bh h r V = pr
. Te use of a calculator is not permitted.. All variables and expressions used represent real numbers unless oterwise indicated.. Figures provided in tis test are drawn to scale unless oterwise indicated..
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 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 informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
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 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 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 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 informationarxiv: v1 [math.na] 3 Nov 2011
arxiv:.983v [mat.na] 3 Nov 2 A Finite Difference Gost-cell Multigrid approac for Poisson Equation wit mixed Boundary Conditions in Arbitrary Domain Armando Coco, Giovanni Russo November 7, 2 Abstract In
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 informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
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 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 informationCalculus I Practice Exam 1A
Calculus I Practice Exam A Calculus I Practice Exam A Tis practice exam empasizes conceptual connections and understanding to a greater degree tan te exams tat are usually administered in introductory
More information1 Upwind scheme for advection equation with variable. 2 Modified equations: numerical dissipation and dispersion
1 Upwind sceme for advection equation wit variable coefficient Consider te equation u t + a(x)u x Applying te upwind sceme, we ave u n 1 = a (un u n 1), a 0 u n 1 = a (un +1 u n ) a < 0. CFL condition
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 informationCHAPTER (A) When x = 2, y = 6, so f( 2) = 6. (B) When y = 4, x can equal 6, 2, or 4.
SECTION 3-1 101 CHAPTER 3 Section 3-1 1. No. A correspondence between two sets is a function only if eactly one element of te second set corresponds to eac element of te first set. 3. Te domain of a function
More informationHigher Derivatives. Differentiable Functions
Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.
More informationMath 34A Practice Final Solutions Fall 2007
Mat 34A Practice Final Solutions Fall 007 Problem Find te derivatives of te following functions:. f(x) = 3x + e 3x. f(x) = x + x 3. f(x) = (x + a) 4. Is te function 3t 4t t 3 increasing or decreasing wen
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 information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More information(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?
Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
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 informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationThese errors are made from replacing an infinite process by finite one.
Introduction :- Tis course examines problems tat can be solved by metods of approximation, tecniques we call numerical metods. We begin by considering some of te matematical and computational topics tat
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationNumerical Solution to Parabolic PDE Using Implicit Finite Difference Approach
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,
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 informationChapter 2 Limits and Continuity
4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationMath 1241 Calculus Test 1
February 4, 2004 Name Te first nine problems count 6 points eac and te final seven count as marked. Tere are 120 points available on tis test. Multiple coice section. Circle te correct coice(s). You do
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 information1watt=1W=1kg m 2 /s 3
Appendix A Matematics Appendix A.1 Units To measure a pysical quantity, you need a standard. Eac pysical quantity as certain units. A unit is just a standard we use to compare, e.g. a ruler. In tis laboratory
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 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 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 informationName: Answer Key No calculators. Show your work! 1. (21 points) All answers should either be,, a (finite) real number, or DNE ( does not exist ).
Mat - Final Exam August 3 rd, Name: Answer Key No calculators. Sow your work!. points) All answers sould eiter be,, a finite) real number, or DNE does not exist ). a) Use te grap of te function to evaluate
More informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
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 informationSECTION 1.10: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES
(Section.0: Difference Quotients).0. SECTION.0: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES Define average rate of cange (and average velocity) algebraically and grapically. Be able to identify, construct,
More information2.3 Algebraic approach to limits
CHAPTER 2. LIMITS 32 2.3 Algebraic approac to its Now we start to learn ow to find its algebraically. Tis starts wit te simplest possible its, and ten builds tese up to more complicated examples. Fact.
More informationConvexity and Smoothness
Capter 4 Convexity and Smootness 4.1 Strict Convexity, Smootness, and Gateaux Differentiablity Definition 4.1.1. Let X be a Banac space wit a norm denoted by. A map f : X \{0} X \{0}, f f x is called a
More informationf a h f a h h lim lim
Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point
More informationLogarithmic functions
Roberto s Notes on Differential Calculus Capter 5: Derivatives of transcendental functions Section Derivatives of Logaritmic functions Wat ou need to know alread: Definition of derivative and all basic
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 informationUNIVERSITY OF MANITOBA DEPARTMENT OF MATHEMATICS MATH 1510 Applied Calculus I FIRST TERM EXAMINATION - Version A October 12, :30 am
DEPARTMENT OF MATHEMATICS MATH 1510 Applied Calculus I October 12, 2016 8:30 am LAST NAME: FIRST NAME: STUDENT NUMBER: SIGNATURE: (I understand tat ceating is a serious offense DO NOT WRITE IN THIS TABLE
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 informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More informationCHAPTER 2 MODELING OF THREE-TANK SYSTEM
7 CHAPTER MODELING OF THREE-TANK SYSTEM. INTRODUCTION Te interacting tree-tank system is a typical example of a nonlinear MIMO system. Heiming and Lunze (999) ave regarded treetank system as a bencmark
More informationChapter 1D - Rational Expressions
- Capter 1D Capter 1D - Rational Expressions Definition of a Rational Expression A rational expression is te quotient of two polynomials. (Recall: A function px is a polynomial in x of degree n, if tere
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 informationExcursions in Computing Science: Week v Milli-micro-nano-..math Part II
Excursions in Computing Science: Week v Milli-micro-nano-..mat Part II T. H. Merrett McGill University, Montreal, Canada June, 5 I. Prefatory Notes. Cube root of 8. Almost every calculator as a square-root
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim1887@aol.com rev 1 Aug 8, 216 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
More information1.5 Functions and Their Rates of Change
66_cpp-75.qd /6/8 4:8 PM Page 56 56 CHAPTER Introduction to Functions and Graps.5 Functions and Teir Rates of Cange Identif were a function is increasing or decreasing Use interval notation Use and interpret
More informationMath Test No Calculator
Mat Test No Calculator MINUTES, QUESTIONS Turn to Section of your answer seet to answer te questions in tis section. For questions -, solve eac problem, coose te best answer from te coices provided, and
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 informationSFU UBC UNBC Uvic Calculus Challenge Examination June 5, 2008, 12:00 15:00
SFU UBC UNBC Uvic Calculus Callenge Eamination June 5, 008, :00 5:00 Host: SIMON FRASER UNIVERSITY First Name: Last Name: Scool: Student signature INSTRUCTIONS Sow all your work Full marks are given only
More informationA.P. CALCULUS (AB) Outline Chapter 3 (Derivatives)
A.P. CALCULUS (AB) Outline Capter 3 (Derivatives) NAME Date Previously in Capter 2 we determined te slope of a tangent line to a curve at a point as te limit of te slopes of secant lines using tat point
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 informationEXTENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS WITH APPLICATION TO AN AEROACOUSTIC PROBLEM
SIAM J. SCI. COMPUT. Vol. 26, No. 3, pp. 821 843 c 2005 Society for Industrial and Applied Matematics ETENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS
More informationBob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk
Bob Brown Mat 251 Calculus 1 Capter 3, Section 1 Completed 1 Te Tangent Line Problem Te idea of a tangent line first arises in geometry in te context of a circle. But before we jump into a discussion of
More information