Finite Difference Methods Assignments
|
|
- Madison Dixon
- 5 years ago
- Views:
Transcription
1 Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling
2 Assignment 1: A one-dimensional eat equation Subtask 1 a) Use te forward Euler sceme to solve te one-dimensional eat problem u t = λu xx 0 <x<1, t > 0 u(x, 0) = f(x) (1) u(0,t) = u(1,t)=0 were Compute and plot te values u n f(x) =sinπx +0.3sin2πx 0.1sin3πx λ =0.1, =0.1, k =0.05 for n =0, 10, 20, 30, 40, 50. b) Te same problem as in a) but wit k =0.1. Make comments on te results. c) Solve te same problem as in a) but wit λ =0.05. Coose k as large as possible but keep λk/ 2 1/2. Plot te results and compare tem wit a). How does reduction in te value of te parameter λ affect te results? Make experiments wit different values of λ. d) Solve te same problem as in a) wit te backward Euler metod. Tis gives a system of equations wit N 1 unknowns at every time step. Te resulting system is a tridiagonal system and can be solved wit MATLAB routines. Begin by considering te coefficient matrix. Verify by experiments te unconditional stability of te sceme. e) Solve te same problem as in d) but wit te Crank Nicolson metod. Subtask 2 Our simple model problem (1) can be generalized in different ways. a) Cange te boundary values to u(0,t)=g 1 (t) and u(1,t)=g 2 (t) were g 1 and g 2 are cosen in a proper way and compute te solution wit te forward Euler sceme. 1 Institutionen för teknisk databeandling
3 b) Let λ be a function of x and/ or t. Try wit λ = e x, λ = e t, λ = e x+t and oters values in combination wit backward Euler. Are tere any computational complications in comparison wit te case wen λ constant? c) Let λ =2/(2+u 2 (x, t)) i.e. let λ depend on te temperature u and compute tesolutionwitteforwardeulermetod. d) Try wit λ = 1. Wat appens? e) Try wit u t = λu xx + au x + bu + f and replace u x wit (u n +1 un 1 )/2 in forward Euler and similarly for backward Euler and Crank Nicolson. Assignment 2: A two-dimensional eat equation Nowweaveteequationu t = λ 1 u xx + λ 2 u yy in te region 0 <x,y<1, t>0. Te constants λ 1 and λ 2 are bot positive. Suppose we ave te initial temperature u(x, y, 0) = f(x, y), e.g. = sin πx sin πy and tat u(x, y, t) = 0 on te boundary. a) Define forward Euler, implement it and compute te solution for some combination of x, y, and k. Try to find a stability condition for λ 1 = λ 2. b) Implement te backward Euler metod. How many unknowns do we get on eac time level? c) Try te ADI-metod (ADI=Alternating Direction Implicit) given by /2 i i u n i k/2 /2 i k/2 /2 i+1, = λ 1 /2 i+1, = λ 1 2/2 i 2 x 2/2 i 2 x + /2 i 1, + /2 i 1, + λ 2 u n i,+1 2u n i + u n i, 1 2 y i,+1 + λ 2un+1 i + i, y Here we ave used an extra time level wit time index n +1/2. Instead of solving one very large linear system we solve a number of smaller systems. 2 Institutionen för teknisk databeandling
4 Assignment 3: A yperbolic problem Subtask 1 We study te PDE-problem u t = u x 0 <x<1, 0 <t u(x, 0) = f(x) =sin2πx u(x, t) =u(x +1,t) wic as te solution u(x, t) =f(x + t) a) Use te leap-frog sceme u n 1 2k = un +1 un 1 2 and try te following combinations of and k 1) =0.1, k = ) =0.1, k = ) =0.05, k = ) =0.05, k = ) =0.1, k =0.1 Use te exact solution for t = k. Compare te computed solution wit te correct solution. Try to find ow te quality depends on te value of λ = k/. In cases of instability, on wic time level is it first observed? For wic x? Make observations of te propagation of perturbations. Subtask 2 a) We cange te boundary condition to u(1, t) = 0. Te PDE-problem is still well posed and as te solution { f(x + t) for x + t 1, u(x, t) = 0 oterwise. But te leap-frog sceme needs a boundary condition also on te left. Coose simple extrapolation 0 =2u n 1 u n 1 2. Wic effects does tis numerical boundary condition ave? 3 Institutionen för teknisk databeandling
5 b) To reduce te perturbations we cange te difference sceme to = u n 1 +2k un +1 u n 1 2 δ(u n un un 1 4u n un 1 2 ) Tis is leap-frog wit a dissipative term and it can be applied in x 2,x 3,..., x N 2. In x 1 and x N 1 te simple leap-frog sceme is used. Use =0.05 and k =0.04 and coose different values of δ between 0 and 0.2 and try to find an optimal value of δ. Assignment 4: A gas flow problem One dimensional gas flow can be described by te yperbolic system ρ m e t + m ρu 2 + p (e + p)u x = (2) Here ρ is te density of te gas, m te momentum, and e te total energy. Te velocity is u = m/ρ. Under te assumption tat we ave an ideal gas, te pressure is given by p =0.4(e 1 2 m2 /ρ). Te system can be written in te form U t + F (U) x =0, were te vectors U and F (U) can be identified in (2). Te initial state is U(0,x)= for x<0 and U(0,x)= for x>0 Subtask 1 Implement and solve te problem in te interval 2 <x<2uptot =0.6. Use te leap-frog wit a dissipative term and approximately 400 grid points 4 Institutionen för teknisk databeandling
6 in space. Try different values of δ to find a solution wit as few oscillations as possible. To obtain stability t sould be cosen suc tat t 0.2 x. In Figure 1 te solution, at a time t>0(nott =0.6), is sown. Note tat suc a good resolution cannot be obtained wit te metods used in tis assignment. 1.4 Density 4 Pressure 1.6 Velocity Figure 1: Te solution of te gas flow example (2) Subtask 2 Explain ow te upper limit of t is found. Te limit as someting to do wit te spectral radius of te Jacobian, F/ U, wic is around 4.7. Assignment 5: An elliptic problem Subtask 1 A simple elliptic PDE-problem is { uxx + u yy = f(x, y) 0 <x,y<1 u(x, y) = g(x, y) on te boundary Tis is te Diriclet problem for Poisson s equation. Te step lengts are x = y = 1/N and te most natural difference approximation is te five-point formula u i+1, 2u i, + u i 1, x 2 + u i,+1 2u i, + u i, 1 y 2 = f(x i,y ) 5 Institutionen för teknisk databeandling
7 Wen tis is applied we get a linear system wic is sparse, i.e. most of te elements of te coefficient matrix are zeros. a) First, find te linear system for N = 4 and ten do te same for a general value of N. How is te coefficient matrix canged if x y? b) Compute te solution and plot it for N =10, 20, 30. Make own coices of f(x, y) andg(x, y). Assignment 6: Air Quality Modeling and te Advection Diffusion Equation I Introduction An important area in modern Environmental Engineering is te study of various air pollutants. Te concentrations of tese pollutants are described by air quality models wic are often formulated as partial differential equations. Wit te use of models te ope is to predict ow peak concentrations will cange in response to predefined canges in te source of pollution. Consider now te concentration u(x, y, z, t) of a gaseous compound. If we ave knowledge of te concentration at time t =0,wewouldliketobeabletopredict future concentrations. Let, for instance, u be te concentration of a noxious gas formed at an industrial plant. At t = 0 a concentrated cloud of te pollutant is released into te air surrounding te plant. Of great importance is te knowledge of ow concentrated te gas will be wen it reaces te nearby residential areas. To model te air quality (i.e. equation, te concentration of u) we use te continuity t = J (3) were is te nabla operator, =(,, ). Equation (3) states tat te x y z time development of u is related to te flux J (amount per area and time) of te gas. Te total flux is made up of two terms. First we ave te dispersive effect of diffusion, given by Fick s first law, J diffusion = D u, were D denotes te diffusion constant. In addition to diffusion te gas is transported by te wind troug a process called advection. Tis leads to a flux, J advection = v u, werev 6 Institutionen för teknisk databeandling
8 is te wind vector. If we now combine te expressions for te flux wit equation (3) we get te advection diffusion equation, t were is te Laplace operator, u = 2 u x 2 = D u (v u) (4) + 2 u y u z 2. In some special cases equation (4) can be solved analytically. However, in most cases, we ave to rely on numerical metods. Task 1 To start off gently we can simplify equation (4) in te following way. First we leave out te diffusion term (i.e. D = 0). If in addition to tis we only let te wind blow in te x-direction (v =(v, 0, 0)) we get te following equation: wic can be expanded into t = (v u) x t = v x u v x Finally, if v is te same for all x, tenequation(5)turnsinto t = v x To solve equation (6) numerically te following explicit finite difference sceme 1 can be used: u n k = v un un 1 Use te sceme above and te initial condition { 5 x < 1 u 0 (x) = 0 oterwise to solve (6) for te area defined by { x 15 0 <t 4 1 Standard notation is used in all te scemes. Tis means tat and n are indices for space and time, respectively, and k are step lengts for space and time, respectively. (5) (6) 7 Institutionen för teknisk databeandling
9 Let v = 1 and coose =0.1 andk =0.05 for te step lengts. Plot te concentration profile at t = 4 togeter (in te same figure as sown in figure (2 B)) wit te profiles you obtain in Task 2 and 3. Explain te concentration profile at t = 4. Is tis wat you would expect wit only advection wit a constant wind velocity? Te exact solution of 6 is u(x, t) = u 0 (x vt) (sown in figure(2 A). Does your numerical solution differ from it? Note tat te problem above (equation (6) togeter wit te initial condition) is not a well-posed problem. To make it a well-posed problem we ave to add boundary conditions. Here, and in te oter tasks we will use natural constraints assuming tat te solution is constant outside a predefined interval. If we solve te problem using a large enoug region, like x 15, tis won t lead to any complications. Figure 2: A Te exact solution to equation (6). B Numerical solutions to Task 1, 2 and 3. Task 2 If we now complicate te situation given in Task 1 somewat by letting te wind (still ust blowing in te x-direction) be a function of x but still ignoring diffusion, ten our concentration u can be described by equation (5), given above. An explicit sceme for tis problem is: u n k = v u n u n 1 u n v v 1 (7) a) Deduce te order of accuracy for te sceme (7). 8 Institutionen för teknisk databeandling
10 b) Use te same step lengts, initial conditionandregionasintask1tosolve equation (5) wit te sceme (7). Let te wind velocity be defined as: { x v(x) = 2 /(1 + x 2 ) x>0 0 oterwise Plot te concentration profile at t = 4 togeter (in te same figure) wit te profiles of Task 1 and 3. Task 3 If we enter diffusion into equation (6), we obtain te following equation t = D 2 u x v 2 x wic can be solved using te sceme, (8) u n k = D un +1 2u n + u n 1 2 v un u n 1 (9) a) Solve equation (8), using te same initial condition, region and te same wind as in Task 1. Use =0.1, k =0.004 and D =1. Plot te concentration profile at t = 4 togeter (in te same figure) wit te profiles you obtained in Task 2 and 3. Explain te differences between te tree profiles. b) Te stability conditions for te finite difference scemes can be examined using, for example, te Fourier metod. Te stability condition for te sceme (9) is k (2D + v) 1 2 Te coice of step lengts in Task 3a terefore gives a stable difference sceme. Now, solve te same problem as in Task 3a but use te step lengts ( =0.1, k =1/210) and ( =0.1, k =1/209). Wit tese coices of step-lengts it is not possible to find te solution at exactly t = 4, owever, you can still study te penomenon of unstable difference scemes. Explain te results. 9 Institutionen för teknisk databeandling
11 Assignment 7: Air Quality Modeling and te Advection Diffusion Equation II Task 1 In assignment 6 we ave only looked at cases were te concentration depends on one space variable and time, u = u(x, t). Let s now study wat appens if te wind vector as two components, v =(v 1,v 2, 0). If we set D =0tegeneral advection diffusion equation (4) turns into: t = v 1 x v 2 y u v 1 x u v 2 y (10) Equation (10) can be solved using a generalization of te sceme (7): (i, ) u n (i, ) k = v 1 (i, ) un (i, ) u n (i 1,) v 2 (i, ) un (i, ) u n (i, 1) u n (i, ) v 1(i, ) v 1 (i 1,) u n (i, ) v 2(i, ) v 2 (i, 1) (11) Use tis sceme to solve (10) wit te following specifications. Solve for x 30, y 30, t 15, step lengts =0.5, k =0.1 and te initial condition: { 50(1 + cos πr u 0 (x, y) = ) r oterwise were { r 2 =(x x 0 ) 2 +(y y 0 ) 2 (x 0,y 0 )=(0, 20) Te wind is defined by { R = x2 + y 2 v =( (y +0.5x)/R, (x 0.5y)/R) Make plots of te solution at t =0,t =5,t = 10, and t = 15. Make a velocity plot of te wind, i.e. plot arrows indicating te wind direction (use quiver in matlab). Try to explain your results. Cange te step size to k = Wat appens? 10 Institutionen för teknisk databeandling
12 Task 2 Add te dissipation term wit D = 1.0. Approximate te Laplace operator wit second order centered differences. Re-do te calculations wit = 0.5 and k =0.01. Compare wit task 1 and explain bot te pysical and te numerical beavior. Task 3 Suggest a better numerical metod for tis problem and motivate your coice, (possibly by sowing tat it works better). 11 Institutionen för teknisk databeandling
Finite 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 1. Rationalize the denominator and fully simplify the radical expression 3 3. Solution: = 1 = 3 3 = 2
MTH - Spring 04 Exam Review (Solutions) Exam : February 5t 6:00-7:0 Tis exam review contains questions similar to tose you sould expect to see on Exam. Te questions included in tis review, owever, are
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 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 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 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 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 Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
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 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 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 informationLab 6 Derivatives and Mutant Bacteria
Lab 6 Derivatives and Mutant Bacteria Date: September 27, 20 Assignment Due Date: October 4, 20 Goal: In tis lab you will furter explore te concept of a derivative using R. You will use your knowledge
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
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 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 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 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 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 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 informationQuantization of electrical conductance
1 Introduction Quantization of electrical conductance Te resistance of a wire in te classical Drude model of metal conduction is given by RR = ρρρρ AA, were ρρ, AA and ll are te conductivity of te material,
More informationPreface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser
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 informationNotes on wavefunctions II: momentum wavefunctions
Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles
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 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 informationSome Applications of Fractional Step Runge-Kutta Methods
Some Applications of Fractional Step Runge-Kutta Metods JORGE, J.C., PORTERO, L. Dpto. de Matematica e Informatica Universidad Publica de Navarra Campus Arrosadia, s/n 3006 Pamplona Navarra SPAIN Abstract:
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 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 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 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 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 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 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 informationChapter Seven The Quantum Mechanical Simple Harmonic Oscillator
Capter Seven Te Quantum Mecanical Simple Harmonic Oscillator Introduction Te potential energy function for a classical, simple armonic oscillator is given by ZÐBÑ œ 5B were 5 is te spring constant. Suc
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 informationPrecalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!
Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.
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 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 informationMAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points
MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim
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 informationThird part: finite difference schemes and numerical dispersion
Tird part: finite difference scemes and numerical dispersion BCAM - Basque Center for Applied Matematics Bilbao, Basque Country, Spain BCAM and UPV/EHU courses 2011-2012: Advanced aspects in applied matematics
More information1. State whether the function is an exponential growth or exponential decay, and describe its end behaviour using limits.
Questions 1. State weter te function is an exponential growt or exponential decay, and describe its end beaviour using its. (a) f(x) = 3 2x (b) f(x) = 0.5 x (c) f(x) = e (d) f(x) = ( ) x 1 4 2. Matc te
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 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 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 informationLesson 6: The Derivative
Lesson 6: Te Derivative Def. A difference quotient for a function as te form f(x + ) f(x) (x + ) x f(x + x) f(x) (x + x) x f(a + ) f(a) (a + ) a Notice tat a difference quotient always as te form of cange
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 information1. Which one of the following expressions is not equal to all the others? 1 C. 1 D. 25x. 2. Simplify this expression as much as possible.
004 Algebra Pretest answers and scoring Part A. Multiple coice questions. Directions: Circle te letter ( A, B, C, D, or E ) net to te correct answer. points eac, no partial credit. Wic one of te following
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 information0.1 Differentiation Rules
0.1 Differentiation Rules From our previous work we ve seen tat it can be quite a task to calculate te erivative of an arbitrary function. Just working wit a secon-orer polynomial tings get pretty complicate
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 informationRobotic manipulation project
Robotic manipulation project Bin Nguyen December 5, 2006 Abstract Tis is te draft report for Robotic Manipulation s class project. Te cosen project aims to understand and implement Kevin Egan s non-convex
More informationEntropy and the numerical integration of conservation laws
Pysics Procedia Pysics Procedia 00 2011) 1 28 Entropy and te numerical integration of conservation laws Gabriella Puppo Dipartimento di Matematica, Politecnico di Torino Italy) Matteo Semplice Dipartimento
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 informationMathematics 123.3: Solutions to Lab Assignment #5
Matematics 3.3: Solutions to Lab Assignment #5 Find te derivative of te given function using te definition of derivative. State te domain of te function and te domain of its derivative..: f(x) 6 x Solution:
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 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 informationFunctions of the Complex Variable z
Capter 2 Functions of te Complex Variable z Introduction We wis to examine te notion of a function of z were z is a complex variable. To be sure, a complex variable can be viewed as noting but a pair of
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 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 informationMATH 1020 Answer Key TEST 2 VERSION B Fall Printed Name: Section #: Instructor:
Printed Name: Section #: Instructor: Please do not ask questions during tis exam. If you consider a question to be ambiguous, state your assumptions in te margin and do te best you can to provide te correct
More informationMath 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0
3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,
More informationParabolic PDEs: time approximation Implicit Euler
Part IX, Capter 53 Parabolic PDEs: time approximation We are concerned in tis capter wit bot te time and te space approximation of te model problem (52.4). We adopt te metod of line introduced in 52.2.
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
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 information= h. Geometrically this quantity represents the slope of the secant line connecting the points
Section 3.7: Rates of Cange in te Natural and Social Sciences Recall: Average rate of cange: y y y ) ) ), ere Geometrically tis quantity represents te slope of te secant line connecting te points, f (
More informationWe name Functions f (x) or g(x) etc.
Section 2 1B: Function Notation Bot of te equations y 2x +1 and y 3x 2 are functions. It is common to ave two or more functions in terms of x in te same problem. If I ask you wat is te value for y if x
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
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 informationIEOR 165 Lecture 10 Distribution Estimation
IEOR 165 Lecture 10 Distribution Estimation 1 Motivating Problem Consider a situation were we ave iid data x i from some unknown distribution. One problem of interest is estimating te distribution tat
More information