Control Errors in CFD!
|
|
- Lewis Newman
- 5 years ago
- Views:
Transcription
1 Control Errors in CFD! Bernard Müller Department of Scientic Computing, Information Tecnolog, Uppsala Universit, Box 337, S Uppsala, Sweden, bernd@tdb.uu.se, URL ttp:// bernd/ Abstract Error control in computational uid dnamics (CFD) as been crucial for reliabilit and ecienc of numerical ow simulations. Te roles of truncation and rounding errors in difference approximations are discussed. Truncation error control is reviewed for ODEs. For dierence approximations of PDEs, discretization error control b Ricardson extrapolation is outlined. Applications to anisotropic grid adaptation in CFD are sown. Alternative approaces of error control in CFD are mentioned. 1 Introduction Computational uid dnamics (CFD) as been maturing using wind tunnel measurements for validation. However, experiments are often not available, in particular in new multidisciplinar applications. Terefore, te assessment of matematical models as become crucial. Tus, numerical errors in te discretization of matematical models ave to be minimized. Error control is essential for te reliabilit of CFD. Moreover, error control is te prerequisite for adaptive metods, wic are decisive for ecienc. Eriksson et al. [] cite te pilosoper Ludwig Wittgenstein ( ): On wat ou cannot compute wit error control, ou must be silent. Since Wittgenstein was an aeronautical engineer before turning to pilosop, e was well aware of te need for reliable computations in aeronautics. Truncation and Rounding Errors Let us rst consider te approximation of te rst derivative df(x) of a given function f(x) b nite dierences. As an example, we take te second-order central dierence approximation D f(x) = f(x + ) f(x ) (1) were is te step size. Inserting te Talor series expansions around x into (1), we obtain te truncation error T (x) := D f(x) df(x) = f (x) +O( 4 ) 6 () sowing te second-order accurac. Te truncation error indicates ow to get iger accurac: -renement Te grid size is reduced. For example, if te grid size is alved, i.e. is used instead of, te truncation error of te second-order central dierence approximation is reduced b a factor of = 1 4, 1 i.e. T (x) T (x). p-renement Te order of accurac p is increased. E.g. instead of te second-order central dierence approximation wit p =, we can use a iger order metod, e.g. p = 4, for wic te truncation error is lower. For a dierence metod of order p, te truncation error will be reduced b a factor of p, wen te grid size is alved. Tus, te accurac is te more increased b grid renement te larger p is. However, no matter wat te order of p is, we cannot see te expected error reduction as. Ten, we even observe an increase of te error. Te reason for te error increase as is caused b rounding errors, wic lead to cancellation errors. Te errors of second-, fourt- and sixt-order central difference approximations of dexp() are sown in 1
2 1 Error of nd, 4t and 6t order central FDM for d exp() / d x 1 nd order FDM 4t order FDM 6t order FDM 1 Inlupp 1,. central single central double error log(absolute error of FDM for d exp() / d x log() Figure 1: Error of second-, fourt- and sixt-order central dierence approximations in double precision. of dexp() Figure : Error of second-order central dierence approximation of dexp() wit single and double precision. Fig. 1. Te nite dierences were computed in double precision. Te error of D exp(), te second-order central dierence approximation of dexp(), is sown in Fig.. Te inuence of single and double precision on te rounding error is clearl visible. An real number u is represented as a oating point number fl(u) on te computer. Te following relation olds: were ɛ ɛ m fl(u) = u(1 + ɛ) { 1 7 for single precision 1 16 for double precision were ɛ m is macine epsilon. Wen computing te second-order central dierence approximation, te dominant rounding error is te cancellation error caused b computing te dierence of almost equal numbers. Neglecting te error caused b te oating point division b, we get on te computer instead of (1) D fl(f(x)) = = (D f(x))[1 + were f = f(x + ) f(x ) and ɛ ± te relative error of fl(f(x ± )). We see tat te relative cancellation error is large, if f(x±) f are large and cannot be reduced sucientl b multiplication b ɛ ±, for wic ɛ ± ɛ m olds. Tus, cancellation errors and rounding errors in general can be better kept under control b using double precision tan b single precision. However, even wit double precision, we sould be aware of rounding errors, wic can pollute numerical results, cf. Fig.s 1 and. In CFD, cancellation errors occur, wen we are dealing wit ows, for wic te canges are minute compared wit te mean ow. An example is te computation of compressible low Mac number ow. Wen solving te compressible Euler and Navier-Stokes equations for Mac numbers M 1, te termodnamic quantities like densit ρ, pressure p and temperature T onl cange ver sligtl compared wit teir stagnation values. Te consequence is tat te computation of e.g. p can lead to large cancellation errors, wic make te results useless. Fig. 3 fl(f(x + )) fl(f(x )) sows te expansion wave in a tube computed f(x + ) ɛ + f(x ) wit te axismmetric Navier-Stokes equations at low Mac numbers [1], [9]. Te ɛ ] f f re-
3 (x,) extrapolated value (x,) Figure 3: Erroneous pressure distribution of an expansion fan in a tube computed wit conventional formulation of te Navier-Stokes equations leading to cancellation errors in single precision at M = O(1 6 ) [1], [9], : computed result, : analtical solution. sults wit single precision are completel erroneous owing to cancellation errors. However, if te Navier-Stokes equations are formulated in terms of te perturbation variables wit respect to stagnation variables like ρ = ρ ρ, te cancellation errors can be minimized [1], [9]. Computations at M = O(1 1 ) can be performed wit single precision using te perturbation formulation, wic is matematicall equivalent to te original Navier-Stokes equations but formulated in a form better suited for computation, cf. [1], [9]. You simpl solve for e.g. densit cange instead of densit wit te same conservation laws. 3 Truncation Error Control for ODEs Let us consider an initial value ordinar differential equation (ODE) { = f(x, ), x > () = (3) were x denotes time and = (x) is te unknown function to be determined. We want to control te step size of a numerical metod to solve ODE (3) based on te local truncation Figure 4: Basic idea of Ricardson extrapolation. error. We sall consider Ricardson extrapolation and Runge-Kutta-Felberg metods. 3.1 Ricardson Extrapolation Te Ricardson extrapolation is a general metod to improve te accurac of numerical results. Te basic idea of extrapolating results wit step sizes and to = can be summarized as follows, cf. Fig. 4: A metod of order p is used, meaning tat we know te teoretical beaviour of te truncation error as. Compute a solution in x k using step size : (x k, ) Compute a solution in x k using step size : (x k, ) Given two values wit dierent error and knowing tat te truncation error beaves as O( p ) it is possible to extrapolate an approximation of te exact value (x k, ). For a metod of order p, te exact value (solution) can be expressed as (x k ) = (x k, ) + c p + O( q ), 3
4 were c is a constant and p and q are integers wit p < q. Solving te problem wit step sizes and ields (x k ) = (x k, ) + c p + O( q ) (x k ) = (x k, ) + c() p + O( q ) Subtracting tese equations, we obtain te leading term of te truncation error c p = (x k, ) (x k, ) p 1 + O( q ). (4) Tus, adding te leading error term to (x k, ), we obtain a better approximation of te exact solution (x k ), because (x k ) = (x k, )+ (x k, ) (x k, ) p +O( q ) 1 Te algoritm for truncation error control based on Ricardson extrapolation becomes: 1. Compute k using a metod of order p wit step sizes and.. Estimate te truncation error b te leading term c p wit Ricardson extrapolation (4). 3. If error c p < tol were tol is a prescribed tolerance, accept k. Compute new step size k+1, were k+1 > k. 4. If error c p > tol, discard k and compute new k less ten te previous one. Compute k. Continue until error c p < tol. If te truncation error is T (x k ) c p, we coose te new step size k+1 suc tat te estimated truncation error is approximatel equal to te desired tolerance tol, i.e. T (x k+1 ) tol. Tus, te new step size is computed from k+1 = ( ) Θ tol 1/p k, T (x k ) were Θ < 1 is a safet factor. 3. Runge-Kutta-Felberg Metods Instead of computing (x k, ) wit two dierent step sizes, one can compute (x k, ) wit two metods aving dierent orders of accurac. Te dierence of te two results can be used to estimate te truncation error. Te computational work is not doubled, if Runge- Kutta-Felberg metods are used. For example, 7 stages are performed and two dierent weigtings of te stages ield two metods of orders 4 and 5, respectivel. Ten, te truncation error can be estimated b te dierence of te two results. Tis error control is e.g. used in te MATLAB ODE-routine ode45. 4 Discretization Error Control for PDEs Consider an initial boundar value problem for a partial dierential equation (PDE) like te Navier-Stokes equations Lu =, (5) were L denotes te dierential operator and u = u(x, t) is te unknown function to be determined. Here, x and t denote te spatial coordinate and time, respectivel. As an example, we consider te scalar linear 1D advection equation wit initial and boundar conditions u t + au x =, x [, 1], t >, (6) u(x, ) = f(x), x [, 1], u(, t) = g(t), t, were a > is a constant advection velocit. We coose a nite dierence metod to discretize te PDE (5) L u =, (7) e.g. te rst-order explicit upwind sceme wit v = u for (6) v n+1 j v n j t + a vn j vn j 1 x =, (8) 4
5 were j = 1,..., N = 1 x, n =,..., t t 1. Insert te exact solution u of te PDE (5) into te FDM (7) to obtain te truncation error T = L u = L u Lu. (9) Inserting u(x j, t n ) for v n j into (8), we obtain T (x j, t n ) = u tt(x j, t n ) t + O( t ) a u xx(x j, t n ) x + O( x )(1) = c x + O( x ) + d t + O( t ) Te task is to estimate te leading error terms, i.e. in our example c x and d t. Following Ferm and Lötstedt [3] and Lötstedt et al. [8], we can use Ricardson extrapolation to approximate te truncation error b te discretization error: 1. Compute v n+1 b solving L v n = on ne grid wit step size = x and time step k = t.. Restrict te ne grid solution v to te coarse grid wit step size, i.e. use restriction operator (N assumed even) I vn = [vn, vn, vn 4,..., vn N ]T. Compute discretization error τ = L I vn (11) on te coarse grid wit time step k, i.e. in te example τ (x j, t n ) = vn+1 j vj n t j =, 4,..., N. + a vn j vn j, x (1) 3. Approximate te dierence of te truncation errors on te ne and coarse grids b te dierence of te discretization errors on te ne and coarse grids (same k = t). Ten, te leading spatial truncation error is approximated b Ricardson extrapolation: c p τ = L I vn I L v n p 1. (13) 5 In our example, we ave p = 1. Since L v n =, te rigt and side of (13) is te dierence of te spatial discretizations on te coarse and ne grids, i.e. for our example c τ (x j, t n ) = a vn j vn j a vn j vn j 1 x x 4. Ceck discretization error τ and rene or coarsen grid, if needed. Recalculate starting from 1. until te tolerances tol are met. Tus: Wile τ tol refine, rene grid, i.e. use. Wile τ tol coarsen, coarsen grid, i.e. use. 5. After completing te grid adaptation in computing v n+1 wit te required level of te spatial discretization error, control temporal truncation error as for ODEs, cf. section 3. 5 Anisotropic Grid Adaptation In multi dimensions, te discretization errors in eac coordinate direction are computed b 1. coarsening in x l -direction, l = 1,, 3 in 3D,. computing te dierence approximations of te x l -derivatives, l = 1,, 3, on coarse grids using te ne grid solution, 3. comparing wit ne grid discretizations of te respective derivatives. If te error tolerance for renement (coarsening) is not met, ceck weter renement (coarsening) in one, two or all tree x l - directions is needed. Ten, an anisotropic grid is obtained. Te procedure can be easil extended to - nite dierence metods on general, i.e. non- Cartesian, structured grids. Te extension to structured nite volume metods is outlined.
6 1.5 u u x. Figure 5: Anisotropic adapted grid wit u velocit contours in te transformed domain, adaptation based on discretization error control for incompressible ow over a at plate at Re = 5 b Ferm and Lötstedt [3]. in detail b Ferm and Lötstedt [3]. Dening coarse grid cells b te union of ne grid cells allows to generalize Ferm and Lötstedt's approac to unstructured nite volume metods. Te data structure is considerabl simplied, if te renement and coarsening is not done locall, but blockwise. Terefore, Ferm and Lötstedt use anisotropic grid adaptation wit a block-structured nite volume metod for te D compressible and incompressible Navier-Stokes equations [3]. Examples b Ferm and Lötstedt [3] are presented ere. Stead incompressible ow over a at plate at Renolds number Re = 5 is computed using te articial compressibilit metod. Wile te initial grid ad 34 cells, te adapted grid wit cells is sown in Fig. 5 togeter wit te u- velocit contours. Te renement tolerance was tol refine =.17. Note tat te - coordinate is stretced in Fig. 5. Te u- velocities and te corresponding errors are compared at dierent stations in Fig.s 6 and 7, respectivel. Compared to te initial grid, te error is considerabl reduced wit te adapted anisotropic grid Figure 6: u velocit proles (smbols ) in transformed -coordinate compared wit te Blasius solution (solid lines), adaptation based on discretization error control for incompressible ow over a at plate at Re = 5 b Ferm and Lötstedt [3]. u Figure 7: Error of u velocit contours in te transformed domain wit initial grid (solid lines) and adapted grid (dotted lines), adaptation based on discretization error control for incompressible ow over a at plate at Re = 5 b Ferm and Lötstedt [3]. 6
7 outflow, < outflow, > 1 TE LE TE Figure 8: Anisotropic adapted grid wit u velocit contours in te transformed domain, adaptation based on discretization error control for laminar subsonic ow over a NACA1 airfoil b Ferm and Lötstedt [3] LE and TE denote leading and trailing edges, respectivel x Fig. 8 sows te adapted grid based on discretization error control for laminar subsonic ow over a NACA1 airfoil at 1.5 angle of attack, Mac number M =.5 and Renolds number Re = 3 b Ferm and Lötstedt [3]. Te contour plots in Fig.s 9 and 1 were obtained b tat approac wit te renement tolerance tol refine =.58, 16 cells in te initial grid and 96 cells in te adapted grid (Fig. 8). Te analsis of te dependence of te solution error u u on te discretization error τ and te results b Ferm and Lötstedt [3] indicate tat te discretization error is a useful error indicator. 6 Conclusions A subtle analsis of certain ows, e.g. boundar laer ow, reveals te relation between te discretization errors and te errors of te ow variables [3]. Te conclusion is tat it often suces to control te discretization errors. Oter a posteriori error estimates are based on Talor expansion [4], [5], [7]. u(x) = u()+x T u()+ 1 xt ( T u())x+... (14) Figure 9: u velocit contours wit anisotropic grid adaptation based on discretization error control for laminar subsonic ow over a NACA1 airfoil b Ferm and Lötstedt [3] x Figure 1: Pressure contours wit anisotropic grid adaptation based on discretization error control for laminar subsonic ow over a NACA1 airfoil b Ferm and Lötstedt [3]. 7
8 Te adjoint metod optimizes te grid for solving a PDE b solving te adjoint problem. Tis approac is used wit te Galerkin nite element metod in [6], [], [11], [1]. Oter errors like oscillations at moving uid interfaces in gas mixtures [9], ux discretization errors due to source terms [9], numerical instabilities due to te 'carbuncle' penomenon [1], etc. are fundamental to te numerical metods used and require numerical analsis to understand te reasons for te errors and devise cures. Te most important recommendation is to be critical to te results produced b a CFD code and to use not onl psical insigt but also all te available knowledge on numerical analsis to get reliable results ecientl. Terefore, control errors in CFD! ACKNOWLEDGEMENTS Te contributions b m colleagues Lars Ferm and Per Lötstedt are gratefull acknowledged. I tank Hideaki Aiso, NAL, Japan, for te invitation to te t ANSS and Yoko Takakura for te ospitalit at Toko Noko Universit. References [1] Becker, R., Rannacer, R.: An Optimal Control Approac to a posteriori Error Estimation in Finite Element Metods. Preprint 1-3, IWR Heidelberg, 1. [] Eriksson, K., Estep, D., Hansbo, P., Jonson, C.: Introduction to Adaptive Metods for Dierential Equations. Acta Numerica (1995), pp [3] Ferm, L., Lötstedt, P.: Anisotropic Grid Adaptation for Navier-Stokes' Equations. Tecnical Report -3, Department of Information Tecnolog, Uppsala Universit, June. [4] Jasak, H., Gosman, A.D.: Automatic Resolution Control for te Finite Volume Metod. Part 1: A-posteriori Error Estimates. Submitted to Numerical Heat Transfer,. [5] Jasak, H., Gosman, A.D.: Automatic Resolution Control for te Finite Volume Metod. Part : Adaptive Mes Renement and Coarsening. Submitted to Numerical Heat Transfer,. [6] Jonson, C.: On Computabilit and Error Control in CFD. International J. for Numerical Metods in Fluids, Vol. (1995), pp [7] Lepage, C.Y., Sueric-Gulick, F., Habasi, W.G.: Anisotropic 3-D Mes Adaptation on Unstructured Hbrid Meses. AIAA Paper -859, 4t ASM&E, Jan., Reno, Nevada. [8] Lötstedt, P., Söderberg, S., Ramage, A., Hemmingsson-Fränden,L.: Implicit Solution of Hperbolic Equations wit Space- Time Adaptivit. BIT, Vol. 4 (), pp [9] Müller, B., Jenn, P.: Improving te Reliabilit of Low Mac Number Flow Computations. CFD J., Vol. 9, No. 1 (1), pp [1] Pandol, M., D'Ambrosio, D.: Numerical Instabilities in Upwind Metods: Analsis and Cures for te Carbuncle Penomenon. J. Comput. Psics, Vol. 166 (1), pp [11] Rannacer, R.: Error Control in Finite Element Computations. In Bulgak, H. and Zenger, C. (eds), Proc. Summer Scool Error Control and Adaptivit in Scientic Computing, Kluwer Academic Publisers (1998), pp [1] Sesterenn, J., Müller, B., Tomann, H.: On te Cancellation Problem in Calculating Compressible Low Mac Number Flows. J. Comput. Psics, Vol. 151 (1999), pp
Numerical 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 informationRunge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t))
Runge-Kutta metods Wit orders of Taylor metods yet witout derivatives of f (t, y(t)) First order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous
More informationChapter 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 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 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 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 informationMA2264 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL. By Dr.T.Kulandaivel Department of Applied Mathematics SVCE
MA64 -NUMERICAL METHODS UNIT V : INITIAL VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS B Dr.T.Kulandaivel Department of Applied Matematics SVCE Numerical ordinar differential equations is te part
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 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 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 informationDifferentiation. Area of study Unit 2 Calculus
Differentiation 8VCE VCEco Area of stud Unit Calculus coverage In tis ca 8A 8B 8C 8D 8E 8F capter Introduction to limits Limits of discontinuous, rational and brid functions Differentiation using first
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More information2.3 Product and Quotient Rules
.3. PRODUCT AND QUOTIENT RULES 75.3 Product and Quotient Rules.3.1 Product rule Suppose tat f and g are two di erentiable functions. Ten ( g (x)) 0 = f 0 (x) g (x) + g 0 (x) See.3.5 on page 77 for a proof.
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 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 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 informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
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 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 informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
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 informationAn Adaptive Model Switching and Discretization Algorithm for Gas Flow on Networks
Procedia Computer Science 1 (21) (212) 1 1 1331 134 Procedia Computer Science www.elsevier.com/locate/procedia International Conference on Computational Science, ICCS 21 An Adaptive Model Switcing and
More informationPolynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions
Concepts: definition of polynomial functions, linear functions tree representations), transformation of y = x to get y = mx + b, quadratic functions axis of symmetry, vertex, x-intercepts), transformations
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 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 informationA Feed-Back Approach to Error Control in Finite Element. Methods: Basic Analysis and Examples. Roland Becker and Rolf Rannacher 1
A Feed-Back Approac to Error Control in Finite Element Metods: Basic Analysis and Examples Roland Becker and Rolf Rannacer 1 Summary. A rened approac to residual-based error control in nite element (FE)
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 informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
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 informationA Numerical Scheme for Particle-Laden Thin Film Flow in Two Dimensions
A Numerical Sceme for Particle-Laden Tin Film Flow in Two Dimensions Mattew R. Mata a,, Andrea L. Bertozzi a a Department of Matematics, University of California Los Angeles, 520 Portola Plaza, Los Angeles,
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 informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
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 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 informationarxiv: v1 [math.na] 11 Apr 2016
Hig order approximation to non-smoot multivariate functions Anat Amir David Levin arxiv:164.281v1 [mat.na] 11 Apr 216 April 12, 216 Abstract Approximations of non-smoot multivariate functions return low-order
More informationA Hybridizable Discontinuous Galerkin Method for the Compressible Euler and Navier-Stokes Equations
A Hybridizable Discontinuous Galerkin Metod for te Compressible Euler and Navier-Stokes Equations J. Peraire and N. C. Nguyen Massacusetts Institute of Tecnology, Cambridge, MA 02139, USA B. Cockburn University
More informationMTH-112 Quiz 1 Name: # :
MTH- Quiz Name: # : Please write our name in te provided space. Simplif our answers. Sow our work.. Determine weter te given relation is a function. Give te domain and range of te relation.. Does te equation
More informationGrid-independent large-eddy simulation of compressible turbulent flows using explicit filtering
Center for Turbulence Researc Proceedings of te Summer Program 2010 203 Grid-independent large-edd simulation of compressible turbulent flows using explicit filtering B D. You, S. T. Bose AND P. Moin Te
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 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 informationImproved Rotated Finite Difference Method for Solving Fractional Elliptic Partial Differential Equations
American Scientific Researc Journal for Engineering, Tecnolog, and Sciences (ASRJETS) ISSN (Print) 33-44, ISSN (Online) 33-44 Global Societ of Scientific Researc and Researcers ttp://asrjetsjournal.org/
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 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 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 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 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 informationMATH1901 Differential Calculus (Advanced)
MATH1901 Dierential Calculus (Advanced) Capter 3: Functions Deinitions : A B A and B are sets assigns to eac element in A eactl one element in B A is te domain o te unction B is te codomain o te unction
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
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 informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
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 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 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 information2.3. Applying Newton s Laws of Motion. Objects in Equilibrium
Appling Newton s Laws of Motion As ou read in Section 2.2, Newton s laws of motion describe ow objects move as a result of different forces. In tis section, ou will appl Newton s laws to objects subjected
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat 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 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 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 informationBlanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating
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 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 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 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 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 information1 The Euler Forward scheme (schéma d'euler explicite)
TP : Finite dierence metod for a European option M2 Modélisation aléatoire - Université Denis-Diderot Cours EDP en Finance et Métodes Numériques. December 2016 1 Te Euler Forward sceme (scéma d'euler explicite)
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 informationWind Turbine Micrositing: Comparison of Finite Difference Method and Computational Fluid Dynamics
IJCSI International Journal of Computer Science Issues, Vol. 9, Issue 1, No 1, January 01 ISSN (Online): 169-081 www.ijcsi.org 7 Wind Turbine Micrositing: Comparison of Finite Difference Metod and Computational
More informationProvable Security Against a Dierential Attack? Aarhus University, DK-8000 Aarhus C.
Provable Security Against a Dierential Attack Kaisa Nyberg and Lars Ramkilde Knudsen Aarus University, DK-8000 Aarus C. Abstract. Te purpose of tis paper is to sow tat tere exist DESlike iterated cipers,
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
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 informationAPPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS. Introduction
Acta Mat. Univ. Comenianae Vol. LXVII, 1(1998), pp. 57 68 57 APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS A. SCHMIDT Abstract. Te pase transition between solid and liquid in an
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 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 informationSolving Second Order Linear Dirichlet and Neumann Boundary Value Problems by Block Method
IAENG International Journal o Applied Matematics 43: IJAM_43 04 Solving Second Order Linear Diriclet and Neumann Boundar Value Problems b Block Metod Zanaria Abdul Majid Mod Mugti Hasni and Norazak Senu
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 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 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 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 informationINTRODUCTION TO CALCULUS LIMITS
Calculus can be divided into two ke areas: INTRODUCTION TO CALCULUS Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and
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 informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More information1 ode.mcd. Find solution to ODE dy/dx=f(x,y). Instructor: Nam Sun Wang
Fin solution to ODE /=f(). Instructor: Nam Sun Wang oe.mc Backgroun. Wen a sstem canges wit time or wit location, a set of ifferential equations tat contains erivative terms "/" escribe suc a namic sstem.
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 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 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 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 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 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 information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationEstimation Approach to Ratio of Two Inventory Population Means in Stratified Random Sampling
American Journal of Operational Researc 05, 5(4): 96-0 DOI: 0.593/j.ajor.050504.03 Estimation Approac to Ratio of Two Inventor Population Means in Stratified Random Sampling Subas Kumar Yadav, S. S. Misra,*,
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 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 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 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 informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
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 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 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 information