Lecture Notes on Numerical Schemes for Flow and Transport Problems
|
|
- Prosper Harrison
- 5 years ago
- Views:
Transcription
1 Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute of Technology ITB Summer School 07: Mathematical Modelling on Industrial and Environmental Problems Institut Teknologi Bandung, Bandung, Indonesia -5 July 07
2 Contents Introduction Finite Difference Method for Diffusion Equation 3. The Euler explicit method Stability Backward time center space Theta-method and the Crank-Nicolson method Finite Difference Method for the Convection Equation 7 3. Upwind method Finite Volume Method for Conservative Problems 0 4. Examples on the approximation for f(u) The Riemann Problem Finite volume method for Burgers equation Introduction In this course you will learn numerical methods, mainly finite difference method for the linear partial differential equations, and a little bit on finite volume method for the nonlinear equations. We will also discuss ways to handle boundary conditions: Dirichlet as well as Neumann types. We hope that this will give you a strong background in dealing with various problems that needs numerical solutions. In principal, finite difference is a method to approximate a differential with a ratio of differences. Application of finite difference approximation to differential equation will give us a difference equation. Basic finite difference formulas are as follows: Forward difference: x f x f(x + x) f(x ) x Backward difference: Center difference: x f x f(x ) f(x x) x x f x f(x + x) f(x x) x + O( x), () + O( x), () + O( x ). (3) Task: Apply Taylor expansion of f(x) around x, and derive all finite difference formulas above. There are several important properties of difference equations for partial differential equations that must be considered before numerical computations are performed. For a well-posed differential equation with an analytical solution u(x, t), and the corresponding difference equation with numerical solution u(x, t n ). Relation between differential equation and the corresponding difference equation is the concept of consistency. Relation between the analytical solution u(x, t) and the approximate solution u(x, t n ) is the concept of convergence. These connections are illustrated in the diagram in Figure. The modified differential equation (MDE) can be obtained if we express each terms in the difference equation, in a Taylor series at a base point. Effectively, this changes the difference equation back into a partial differential equation. If the MDE tends to the original partial differential equation as grid spacing tends to zero, the difference equation is called consistent. Terms appearing in the MDE that do not appear in the original partial differential equation are called as
3 Figure : A connection diagram between the differential equation and its difference equation counterpart. the local discretization error or truncation error terms. Order of the local discretization error is the order of the difference scheme. Here are explanation of some basic terminology of a numerical scheme: The order of a difference equation is the order of the local discretization error. A difference equation is consistent with the partial differential equations if the difference between them vanishes as the size of x, t tend to zero. A difference equation is stable if it produces a bounded solution. A finite difference method is convergent if the solution of the difference equation approaches the exact solution of the partial differential equation. Lax equivalence theorem: It states that for a consistent finite difference method for a wellposed linear initial value problem, the method is convergent if and only if it is stable. Proving convergence of a numerical solution is often difficult, it is much easier showing stability and consistency of a difference scheme. Therefore in practice, Lax equivalence theorem is often used. In summary, the concepts of consistency, order, stability and convergence must be considered when solving a partial differential equation by finite difference methods. Consistency and order can be determined from the modified differential equation. Stability can be investigated by applying the von Neumann method. Convergence can be ensured through Lax equivalence theorem by checking consistency and stability. Finite Difference Method for Diffusion Equation In this section, we will discuss common methods of finite difference for the diffusion equation. Consider a diffusion equation with thermal diffusivity κ u t = κu xx, for 0 < x < L, t > 0, (4) (Boundary condition) u(0, t) = 0, u(l, t) = 0, for t > 0, (5) (Initial condition) u(x, 0) = f(x), for 0 < x < L. (6) In this section, we will discuss various finite difference scheme, starting from the simplest FTCS which is explicit and of order one, the BTCS implicit method, and finally the order two, the Crank-Nicolson scheme and theta method. 3
4 . The Euler explicit method On the computational domain Ω = { (x, t) R 0 < x < L, t > 0 }, we introduce a grid, with x denotes the spatial grid size and t the time step. The grid points (x, t n ) are defined as x = ( ) x, =,,, N x +, x = L/N x, (7) t n = (n ) t, n =,, (8) We discretize (4) by replacing the time derivative t u(x, t n ) by the forward difference approx- Figure : Stencil of the Euler explicit (or FTCS) method for the diffusion equation. imation, and the space derivative u(x x, t n ) by the center difference approximation. Writing the numerical approximation of u(x, t) at grid points (x, t n ) as u n, we obtain the difference equation of (4) as t (un+ u n ) = κ x (un + u n + u n ). (9) The scheme can be simplified as u n+ = u n + S(u n + u n + u n ), with S = κ t x. (0).. Stability To analyze the stability of the FTCS scheme, consider a numerical solution of the form u n = e ikx iωt n, () with k is the wave number and ω the wave frequency. Different representation is u n = λ n e ikx, with λ = e iω t. () Substituting () into the scheme (0) will give us λ = + S(e ik x + e ik x ) = 4S sin (k x/) Thus at each time step, the Fourier mode e ikx is amplified by the factor λ. Finite solution u n for all n will be obtained if λ, which is satisfied if S = κ t x. (3) Task: prove that (3) is the stability condition of the Euler explicit scheme (0). We stress here that condition (3) is the necessary condition for the stability of (9). 4
5 Exercise:. As mention before, Local discretization error is the residual d n found as the difference between the differential equation (4) evaluated at grid point (x, t n ) and the scheme (0). Show that d n = κ ( κ t 6 x) u 4x n, which is of the order O( t, x ). Then, discuss the consistency of (0) as the discrete equation of (4).. Write a numerical code for the explicit Euler equation (0). Solve the diffusion equation (4) on 0 < x < with κ =, using boundary conditions (5) and initial condition f(x) =. For computation, use x = 0., t = 0.005, and plot the numerical solution u(x, t = ). Compare the numerical result with the analytical solution (as obtained form the separation variables method). 3. The FTCS method has an interesting feature when S = 6. In this special case the O( t) and O( x ) of the local discretization error d n cancels exactly, and the method becomes O( t, t x, x 4 ). Check the improved accuracy from your numerical results, by comparing the error of numerical solution u(x, ), at several positions, i.e. x = L/4, L/, 3L/4. Aside: Condition (5) is of Dirichlet type. This condition is easily satisfied by taking u n = 0, and u n N x+ = 0, for all n. And we simply implement (0) for =, 3,, N x. Another type of boundary condition that often met is the Neumann boundary condition, which is as follows u x (0, t) = 0, u x (L, t) = 0, for t > 0. (4) Implementation of this boundary conditions needs extra attention. Apply the order-two center difference to approximate (4), from which we can get the relation u( x, t n ) = u( x, t n ). Using this, we can implement (0) for =, for all n. Next, the right boundary condition u x (L, t) = 0 can be solved in analogous way. Note that, in this Neumann boundary conditions case, values u n and un N x+ should be computed, at all time step tn. % (Algorithm for diffusion equation u t = κu xx ) % Initialition % Specify the computational domain, spatial 0 x L and time 0 t T. % Choose the grid size x and t. % Compute N x + and N t, the total number of grid points in x and t. % Define an initial condition u(x, 0). u( : N x +, ) = 0; u(0.5 N x 4 : 0.5 N x +4, ) = ; % This is a square wave initial condition, centered in the middle. % Specify the value of C Courant number. C = κ t/ x ; % Stability condition: Courant number C /. % Computing u(x, t n+ ) using the difference equation (??). for n = : N t u(, n + ) = 0 for = : N x u(, n + ) = u(, n) + S (u(, n) u(, n) + u(, n)); end u(n x +, n + ) = 0 end 5
6 % Plotting the result xplot=0: x:l for n=:nt figure() plot(xplot, u(:,n+)) pause(0.) drawnow end. Backward time center space Implementation of the backward approximation for time derivative t u(x, t n+ ) and center difference approximation for x u(x, t n+ ) will give us the following implicit scheme t (un+ u n ) = κ x (un+ + un+ + u n+ ). (5) The scheme can be simplified as Su n+ + + ( + S)un+ Cu n+ ) = un, with S = κ t x. (6) Figure 3: Stencil of the BTCS implicit method for the diffusion equation. Exercise: The numerical properties of the implicit scheme (6), and its implementation for solving (4-6).. Show that the local discretization error is d n = κ ( κ t x) u 4x n, which is of order O( t, x ).. Show that the scheme (6) is unconditionally stable. 3. In order to compute u n, for every time step, we need to solve a system of equation Aun+ = u n, with ( + S) S 0 0 S ( + S) S 0 A m m = S ( + S) The vector u n has u n as its element. What is the size of the square matrix A, is it N x, N x or N x +? (Remember that if we deal with Dirichlet boundary condition (5), we already have u n = 0 and un N x+ = 0 for all n.) 4. Write a numerical code for solving (4-6) using the implicit scheme (6). Try several combination of x and t to make sure that the BTCS is an unconditionally stable scheme. 6
7 .3 Theta-method and the Crank-Nicolson method In this section we will move on to the second order scheme. This second order scheme is a convex combination of the explicit scheme and implicit scheme. For a parameter θ, with 0 θ, a combination of θ implicit and ( θ) explicit will result in the following scheme t (un+ u n ) = κ θ x (un + u n + u n ) + κ θ x (un+ + un+ + u n+ ). (7) The scheme (7) with θ = is called the Crank-Nicolson scheme. This is a second order scheme with accuracy O( t, x ). Moreover, it is an unconditionally stable scheme. Figure 4: Stencil of the theta-method. Exercise: Numerical properties of the θ scheme (7). Find the local discretization error d n, and show that (7) is an order two scheme with accuracy O( t, x ).. Show that the scheme (7) is unconditionally stable. 3 Finite Difference Method for the Convection Equation In this section we will discuss numerical methods for the equation u t + du x = 0, subect to the initial condition u(x, 0) = f(x), with x R. We recall that the exact solution of this initial value problem is given by u(x, t) = f(x dt). This solution represents a signal propagates undisturbed in shape with velocity d. The signal propagates to the right if d > 0, and to the left if d < 0. We will use this analytical solution as a benchmark for the numerical schemes. Note that these schemes are easily generalized for more complex problems (but analysis are not as easy as this linear case). 3. Upwind method For discussions in this sub section, we address the following problem u t + du x = 0, for 0 < x < L, t > 0, (8) (Initial condition) u(x, 0) = f(x), for 0 < x < L. (9) Here, we restrict to the case of d > 0 and discuss the simplest finite difference scheme, which is the upwind method. On the computational domain Ω = { (x, t) R 0 < x < L, t > 0 }, we introduce a grid, with x denotes the spatial grid size and t the time step. The grid points (x, t n ) are defined as x = ( ) x =,,, N x +, x = L/N x, (0) t n = (n ) t n =,, () 7
8 For the case d > 0 in which signal propagates to the right, upwind method is Forward Time Backward Space (FTBS). Applying FTBS to the convection equation (8) will yield u n+ or explicitly u n t + d un un x = 0. () u n+ = ( C)u n + Cu n, (3) with C d t x a Courant number. From the implementation of forward time and backward space approximation, we can guess already that this upwind method has an accuracy O( t, x). To support this argument, you should show that the local discretization error of the upwind approximation is indeed of the order O( t, x). Please try! Figure 5: Stencil for the upwind method. From the initial condition (9), values of u, for =,,, N x + is known. The FTBS scheme (3) can be used to compute u n, for =,, N x +, for all t n, n =,,, but not u n, because (3) will need un 0 which is not defined. Here, we need to take a numerical boundary condition: u n, for all n. For the first try, one can ust take un = 0, for all n. Later on you can experiment with various left boundary u n = g(tn ), with g(t) plays role as a signal enters from the left boundary. By doing this, we can compute u n, for =,, N x +, for all n using the FTBS scheme (3). % (Algorithm for convection equation u t + du x = 0) % Initialition % Specify the computational domain, spatial 0 x L and time 0 t T. % Choose the grid size x and t. % Compute N x + and N t, the total number of grid points in x and t. % Define an initial condition u(x, 0). u( : N x +, ) = 0; u( : 4, ) = ; % This is a square wave initial condition. % Specify the value of C Courant number. C = d t/ x; % Stability condition: Courant number C. % Computing u(x, t n+ ) using the difference equation (3). for n = : N t u(, n + ) = 0;% Zero left boundary for = : N x + u(, n + ) = ( C)u(, n) + Cu(, n); end end 8
9 % Plotting the result figure() surf(u ) % A surface plot in Matlab command, compare with surf(u) drawnow Exercise. Numerical aspect of the upwind FTBS scheme. Derive the stability condition of the FTBS method. Show that the local discretization error is d(d t x)u xx n, which is of order one.. Write a numerical code for the FTBS scheme (3). For the first simulation take a combination of t and x such that the Courant number C = d t x is exctly one, since d > 0. Using this set up, you should get a square wave propagates to the right. Make sure your wave propagates with the correct velocity d. Also check that the wave influx that enters from the left side of the domain matches with your chosen numerical left boundary. 3. Perform another simulation using different t in such a way that C =., explain the result. Since the Courant number is larger than one, your numerical solution should be unstable, in which the numerical value u n increases as time progresses. 4. Conduct another computation using different t in such a way that the Courant number C is strictly less than one. Your simulation should demonstrate that very soon discontinuity of the square wave disappear, due to diffusion effect. This is a direct consequence of the error diffusive term d(d t x)u xx n. This damping effect will only be observed when d t x 0. For a stable scheme, this damping effect will be observed only when C = d t x strictly less than one. Please check this statement using your numerical code. Exercise. Another idea is to implement the forward time and center space (FTCS) approximation. This scheme seems preferable for two reasons: it may accommodate both cases d > 0 and d < 0, and it has second order accuracy in space: FTCS scheme is of order O( t, x ). Implement stability analysis and argue that this FTCS scheme is always unstable, which means this scheme is useless.. Upwind scheme for arbitrary d. t (un+ u n ) = d+ x (un u n ) d x (un + u n ), (4) with d + max(d, 0) dan d min(d, 0). Implement this scheme using the initial condition: a hump with peak located in the middle of the domain. Take a positif value of d and simulate the propagation of the initial hump to the right. Then, without changing the scheme, take the negative value of d, and show that the initial hump is now propagates to the left. Note: higher stable scheme are Lax O( t, x ), Lax-Wendroff O( t, x ), and leapfrog O( t, x ). Interested reader may consult reference [3, ] for further study. 9
10 4 Finite Volume Method for Conservative Problems In this section we will discuss another variant of numerical method, which is called the finite volume method. Finite volume methods are closely related to finite difference methods, and a finite volume method can often be interpreted directly as a finite difference approximation to the differential equation. However, finite volume methods are derived on the basis of the integral form of the conservation law, a starting point that turns out to have many advantages. Whereas in finite difference methods all derivatives appear in the equation, each is approximated by a ratio of differences. On the computational domain [0, L] which is divided into N x cells of homogeneous length x, a staggered partition points are x / = 0, x 3/ = x,, x +/ = x,, x Nx+/ = N x x = L (5) We assume that u(x, t n ) denotes the averaged value of u in cell C = [x, x + ]. Within time interval t, this value changes according to the net flux from the left and right boundaries of cell C according to ( u(x, t n+ ) x = u(x, t n ) x t f(x ), t n ) f(x +, t n ). (6) Equation (6) is exactly the conservative principle in discrete form of the quantity u on cell C. What we need to do next is ust a good approximation scheme for flux f(x, t n ). For further discussion, we use the following notations u n = u(x, t n ), and f n = f(x, t n ), for all x and t n in the computational domain. 4. Examples on the approximation for f(u). If the flux f(u) = κu x, equation (8) reduces to diffusion equation u t = κu xx. For the approximation, we take an average f n + (f n + + f n ), the discrete equation (6) reduces to u n+ u n = κ t x (f n + f n + f n ), (7) which is exactly the finite difference approximation (0).. If the flux f(u) = du, equation (8) reduces to convection equation u t + du x = 0. Here we take the upwind approximation: if d > 0, f + f = du, and when d > 0, f + f + = du +. Implementing this on (6) will give us the finite difference approximation (4). Prove this! 4. The Riemann Problem We start the discussion by addressing the basic conservation equation of variable u with flux function f(u) as follows u t + f(u) x = 0, for x R, t > 0, (8) { ul, if x < 0, u r, if x > 0. (9) u(x, 0) = 0
11 Writing b(u) = f (u), equation (8) can be written as u t + b(u) u = 0, for x R, t > 0. (30) x Equation (30) is written in the form of convection equation with b(u) represents the signal velocity. Some observation, if u(x, t) is a solution, then so is u(ax, at), for arbitrary a > 0, implying that the solution is a similarity solution of the form u(x, t) = û(x/t). Solutions of (8,9) are propagated along the characteristic x/t = b(u), with b(u) = f (u) represents the velocity. For the initial condition (9), we can distinguish two cases, b(u l ) > b(u r ) and b(u l ) < b(u r ), to be discussed below. Case b(u l ) > b(u r ). The characteristic emanating from x < 0 have a smaller slope than the slope of the characteristic coming from x > 0. As a consequence, there is an area in which the characteristic intersect, which would lead to multi valued solutions. Here, we have a discontinuous solution. We can easily verify that the following is a solution of the Riemann problem (8,9), and it is given by u(x, t) = { ul, if x/t < s, u r, if x/t > s, (3) where s is defined by the Rangkine Hugoniot condition s = f(u r) f(u l ) u r u l (3) Solution of (3) is known as the shock wave, and s as the shock speed. This solution contains discontinuity, and it is considered as weak solution. A typical shock wave solution and the corresponding characteristic are shown in Figure 6. The shock speed formula (3) is determined such that (3) is a solution of the corresponding integral form of (8), see [] for further study. Figure 6: Shock wave and the corresponding characteristic. Case b(u l ) < b(u r ). In this case the characteristic emanating from x < 0 have a larger slope than those emanating from x > 0, so we have separating characteristic, see Figure 7 (right). In this case we do not expect a discontinuous solution. Solution of (30) is given by u l, if x/t < b(u l ) u(x, t) = w(x/t), if b(u l ) < x/t < b(u r ) (33) u r, if x/t > b(u r ).
12 Figure 7: Rarefaction wave and the corresponding characteristic. where w(η) should satisfy b(w(η)) = η. This solution is called rarefaction wave, and its typical profile is depicted in Figure (7) (left). Despite the fact that the initial condition (9) is discontinuous, the rarefaction wave (33) is continuous. Please check! Exercise:. Burgers equation is the (8) with the flux f(u) = u, written explicitly as follows u t + ( ) x u = 0, subect to the piecewise constant initial condition { α, if x < 0, u(x, 0) = β, if x > 0. In this case b(u) = f (u) = u, therefore two type of solutions are distinguish by the relation between b(α) = α and b(β) = β. If α > β, the solution is a shock wave, formulate this shock wave, and determine the shock speed. If α < β, the solution is a rarefaction wave, formulate the solution.. The kinematic LWR model for traffic flow is equation (8) with the chosen flux. The most common flux function used is the Greenshield flux, in which after normalization f(u) = u( u). Similarly, there are two type solutions of the kinematic wave equation u t + (u( u)) = 0, x subect to the initial condition { α, if x < 0, u(x, 0) = β, if x > 0. In this case b(u) = f (u) = u, therefore two type of solutions are distinguish by the relation between b(α) = α and b(β) = β. If α < β, the solution is a shock wave, formulate this shock wave, and determine the shock speed. If α > β, the solution is a rarefaction wave, formulate the solution.
13 4.3 Finite volume method for Burgers equation Consider the Burgers equation which is the conservative equation with f(u) = u, written explicitly as u t + ( ) x u = 0, x R, t > t n. (34) Basically, the most important step in implementing the finite volume method (7) is on the approximation of f n. The approximation, which is based on analytical solution, holds for + each cell V = [x /, x +/ ] in the computational domain, see Figure 8. Therefore, we first discuss the analytical solution of (34) with the following initial condition { u u(x, t n n ) =, x < x +/ u n +, x > x +/. (35) This analytical formula will be used to find the suitable app Figure 8: Sketch of the spatial domain in the finite volume method, and illustration of the piecewise constant numerical solution at time level t n. If u n > un +, the solution is a shock wave propagating with speed s n = (un + u n +). Let us denote η = x x + t t n. Then, the shock wave solution is given by u R (η; u n, u n +) = { u n, if η < s n u n +, if η > sn. (36) 3
14 If u n < un +, the solution is a rarefaction wave as given by u n, if η < un u R (η; u n, u n +) = η, if u n < η < un + (37) u n +, if η > un +. By taking the numerical flux f x + = u R (η; un, un + ), we obtained the following (Show this!) If u n > un +, then f n + = { (un ), if s n > 0 (un + ), if s n < 0 (38) If u n < un + f n + = (un ), if u n > 0 0, if u n < 0 < un + (un + ), if u n + < 0 (39) We resume here that the finite volume method is (7) with flux given by (38,39), this scheme is also known as Godunov method. Exercise: Figure 9: Shock wave solutions at time t = calculated using Godunov method (left) and Upwind method (right), in comparison with the analytical solution.. Solve the Burgers equation subect to the initial condition {, if x < 0, u(x, 0) =, if x > 0. Analytical solution will be a shock wave, with discontinuity that propagates with speed s =.5. In this case the scheme reduces to u n+ = u n t ( x (un ) ) (un ) (40) Implement this Godunov method, and write a numerical code to simulate the propagation of a shock wave solution. Use x =.5 0, x =.5 0 and plot numerical solution together with the analytical solution at time t =. 4
15 . Alternatively, Burgers equation can be written as u t + uu x = 0. (4) Implementation of the upwind scheme, and since u > 0 we obtain the following scheme u n+ = u n t ( x un u n u n ). (4) Write a numerical code to simulate the propagation of a shock wave solution. Use x =.5 0, x =.5 0 and plot numerical solution together with the analytical solution at time t =. It is shown in (9) that the upwind method cannot capture the shock speed correctly. Clearly the Godunov method (40) is much better compare to the upwind scheme (4). The Godunov scheme can capture the speed of the shock wave. This is one of characteristic features of a conservative scheme. Moreover, the scheme (4) is considered as the non-conservative scheme. References [] Matthei, R.M.M., Rienstra, S.W., ten Thie Boonkkamp, J.H.M., Partial Differential Equations, Modelling, Analysis, Computations, SIAM Monographs, Philadelphia, 005 [] Strauss, W., Partial Differential Equations, an Introduction, John Wiley & Sons, 99. [3] Hoffmann, J., Numerical Methods for Engineers and Scientists, Mc.Graw-Hill, 99. 5
Lecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationLast time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler:
Lecture 7 18.086 Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler: U j,n+1 t U j,n = U j+1,n 2U j,n + U j 1,n x 2 Expected accuracy: O(Δt) in time,
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationNumerical Methods for Hyperbolic Conservation Laws Lecture 4
Numerical Methods for Hyperbolic Conservation Laws Lecture 4 Wen Shen Department of Mathematics, Penn State University Email: wxs7@psu.edu Oxford, Spring, 018 Lecture Notes online: http://personal.psu.edu/wxs7/notesnumcons/
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationFinite Difference Methods (FDMs) 2
Finite Difference Methods (FDMs) 2 Time- dependent PDEs A partial differential equation of the form (15.1) where A, B, and C are constants, is called quasilinear. There are three types of quasilinear equations:
More informationApproximations of diffusions. Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 3b Approximations of diffusions Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.1 04/10/2018 1 Learning objectives Become aware of the existence of stability conditions for the
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationLecture 4.5 Schemes for Parabolic Type Equations
Lecture 4.5 Schemes for Parabolic Type Equations 1 Difference Schemes for Parabolic Equations One-dimensional problems: Consider the unsteady diffusion problem (parabolic in nature) in a thin wire governed
More informationBurgers equation - a first look at fluid mechanics and non-linear partial differential equations
Burgers equation - a first look at fluid mechanics and non-linear partial differential equations In this assignment you will solve Burgers equation, which is useo model for example gas dynamics anraffic
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More information0.3.4 Burgers Equation and Nonlinear Wave
16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave
More informationFinite Differences: Consistency, Stability and Convergence
Finite Differences: Consistency, Stability and Convergence Varun Shankar March, 06 Introduction Now that we have tackled our first space-time PDE, we will take a quick detour from presenting new FD methods,
More information12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis
ATH 337, by T. Lakoba, University of Vermont 113 12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis 12.1 Formulation of the IBVP and the minimax property of its
More informationLecture 17: Initial value problems
Lecture 17: Initial value problems Let s start with initial value problems, and consider numerical solution to the simplest PDE we can think of u/ t + c u/ x = 0 (with u a scalar) for which the solution
More information2.2. Methods for Obtaining FD Expressions. There are several methods, and we will look at a few:
.. Methods for Obtaining FD Expressions There are several methods, and we will look at a few: ) Taylor series expansion the most common, but purely mathematical. ) Polynomial fitting or interpolation the
More informationMethod of Lines. Received April 20, 2009; accepted July 9, 2009
Method of Lines Samir Hamdi, William E. Schiesser and Graham W. Griffiths * Ecole Polytechnique, France; Lehigh University, USA; City University,UK. Received April 20, 2009; accepted July 9, 2009 The method
More informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationq t = F q x. (1) is a flux of q due to diffusion. Although very complex parameterizations for F q
! Revised Tuesday, December 8, 015! 1 Chapter 7: Diffusion Copyright 015, David A. Randall 7.1! Introduction Diffusion is a macroscopic statistical description of microscopic advection. Here microscopic
More informationNumerical Algorithms for Visual Computing II 2010/11 Example Solutions for Assignment 6
Numerical Algorithms for Visual Computing II 00/ Example Solutions for Assignment 6 Problem (Matrix Stability Infusion). The matrix A of the arising matrix notation U n+ = AU n takes the following form,
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationMath Partial Differential Equations 1
Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationChapter 10 Exercises
Chapter 10 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rl/fdmbook Exercise 10.1 (One-sided and
More informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationDissipation and Dispersion
Consider the problem with periodic boundary conditions Dissipation and Dispersion u t = au x 0 < x < 1, t > 0 u 0 = sin 40 πx u(0, t) = u(1, t) t > 0 If a > 0 then the wave is moving to the left and if
More informationNumerical Integration of Linear and Nonlinear Wave Equations
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Dissertations, Theses, and Student Research Papers in Mathematics Mathematics, Department of 12-2004 Numerical Integration
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationAvailable online at Eng. Math. Lett. 2014, 2014:17 ISSN: WAVE ATTENUATION OVER A SUBMERGED POROUS MEDIA I.
Available online at http://scik.org Eng. Math. Lett. 04, 04:7 ISSN: 049-9337 WAVE ATTENUATION OVER A SUBMERGED POROUS MEDIA I. MAGDALENA Industrial and Financial Mathematics Research Group, Faculty of
More informationMulti-Factor Finite Differences
February 17, 2017 Aims and outline Finite differences for more than one direction The θ-method, explicit, implicit, Crank-Nicolson Iterative solution of discretised equations Alternating directions implicit
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More information3.3. Phase and Amplitude Errors of 1-D Advection Equation
3.3. Phase and Amplitude Errors of 1-D Advection Equation Reading: Duran section 2.4.2. Tannehill et al section 4.1.2. The following example F.D. solutions of a 1D advection equation show errors in both
More informationBasics of Discretization Methods
Basics of Discretization Methods In the finite difference approach, the continuous problem domain is discretized, so that the dependent variables are considered to exist only at discrete points. Derivatives
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationModule 3: BASICS OF CFD. Part A: Finite Difference Methods
Module 3: BASICS OF CFD Part A: Finite Difference Methods THE CFD APPROACH Assembling the governing equations Identifying flow domain and boundary conditions Geometrical discretization of flow domain Discretization
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationAn Optimal Dimension of Submerged Parallel Bars as a Wave Reflector
BULLETIN of the Malaysian Mathematical Sciences Society http://math.usm.my/bulletin Bull. Malays. Math. Sci. Soc. (2) 32(1) (2009), 55 62 An Optimal Dimension of Submerged Parallel Bars as a Wave Reflector
More informationdifference operators.
2D1263 : Scientific Computing Lecture 8 (2) 2D1263 : Scientific Computing Lecture 8 (1) Difference operators Let xi m i=0 be equidistant points in [a, b], x i = a + (b a)i/m, and define the forward difference
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Modified Equation of a Difference Scheme What is a Modified Equation of a Difference
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 informationNumerical methods Revised March 2001
Revised March 00 By R. W. Riddaway (revised by M. Hortal) Table of contents. Some introductory ideas. Introduction. Classification of PDE's.3 Existence and uniqueness.4 Discretization.5 Convergence, consistency
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationFinite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More informationFundamentals Physics
Fundamentals Physics And Differential Equations 1 Dynamics Dynamics of a material point Ideal case, but often sufficient Dynamics of a solid Including rotation, torques 2 Position, Velocity, Acceleration
More informationNumerical Schemes Applied to the Burgers and Buckley-Leverett Equations
University of Reading Numerical Schemes Applied to the Burgers and Buckley-Leverett Equations by September 4 Department of Mathematics Submitted to the Department of Mathematics, University of Reading,
More informationHomework for Math , Fall 2016
Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant
More informationTraffic Flow Problems
Traffic Flow Problems Nicodemus Banagaaya Supervisor : Dr. J.H.M. ten Thije Boonkkamp October 15, 2009 Outline Introduction Mathematical model derivation Godunov Scheme for the Greenberg Traffic model.
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More information2 Multidimensional Hyperbolic Problemss where A(u) =f u (u) B(u) =g u (u): (7.1.1c) Multidimensional nite dierence schemes can be simple extensions of
Chapter 7 Multidimensional Hyperbolic Problems 7.1 Split and Unsplit Dierence Methods Our study of multidimensional parabolic problems in Chapter 5 has laid most of the groundwork for our present task
More informationAdvection in two dimensions
Lecture 0 Advection in two dimensions 6. Stability of multiple terms (in multiple dimensions) When we analyzed the stability of time-stepping methods we tended to consider either a single damping term
More informationChapter 5. Numerical Methods: Finite Differences
Chapter 5 Numerical Methods: Finite Differences As you know, the differential equations that can be solved by an explicit analytic formula are few and far between. Consequently, the development of accurate
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )
More informationAMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
More informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More information2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
3 Quantitative Properties of Finite Difference Schemes 31 Consistency, Convergence and Stability of FD schemes Reading: Tannehill et al Sections 333 and 334 Three important properties of FD schemes: Consistency
More informationComputational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. #12 Fundamentals of Discretization: Finite Volume Method
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationChapter 3 Burgers Equation
Chapter 3 Burgers Equation One of the major challenges in the field of comple systems is a thorough understanding of the phenomenon of turbulence. Direct numerical simulations (DNS) have substantially
More informationChapter 10. Numerical Methods
Chapter 0 Numerical Methods As you know, most differential equations are far too complicated to be solved by an explicit analytic formula Thus, the development of accurate numerical approximation schemes
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
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 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 informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More information