Finite difference method for heat equation
|
|
- Owen Stanley
- 6 years ago
- Views:
Transcription
1 Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore January 25, / 46
2 Heat equation: Initial value problem Partial differential equation, µ > 0 Exact solution u t = µu xx, (x, t) R R + u(x, 0) = f(x), x R 1 + u(x, t) = e y2 /4µt f(x y)dy =: (E(t)f)(x) 4πµt Solution bounded in maximum norm u(t) C = E(t)f C f C = sup f(x) x R 2 / 46
3 Finite difference mesh in space and time: partition space R by a uniform mesh of size h and time R + by a uniform mesh of size t x = h, =..., 2, 1, 0, +1, +2,... t n = n t, n = 0, 1, 2,... t n t Let the numerical approximation be x x U n u(x, t n ) 3 / 46
4 Finite difference in space First derivative: D + x U n = U n +1 U n h Second derivative:, Dx U n = U n U 1 n, D 0 h xu n = U +1 n U 1 n 2h Finite difference in time D + x D x U n = U n 1 2U n + U n +1 h 2 D t + U n = U n+1 U n t 4 / 46
5 Finite difference scheme: forward time and central space (FTCS) i.e., D + t U n = µd + x D x U n, n = 0, 1, 2,... U n+1 U n t Update equation: Solve for U n+1 = µ U n 1 2U n + U n +1 h 2 U n+1 = (E h U n ) = λu n 1 + (1 2λ)U n + λu n +1 where Explicit scheme λ = µ t h 2 U n+1 = H(U n ) 5 / 46
6 Maximum principle If λ 1/2, then U n+1 is given by a convex linear combination of U 1 n, U n, U +1 n. Hence min(u n 1, U n, U n +1) U n+1 min Uk n U n+1 k max(u n 1, U n, U n +1) max Uk n, k Hence solution at any time level n is bounded by initial data Maximum norm min Uk 0 U n max Uk 0,, n k k U = sup U Maximum stability: If λ 1/2, then U n+1 = E h U n U n... U 0 and hence U n U 0 = f, n 6 / 46
7 Maximum principle We have E h V V Eh n V V λ 1/2 is necessary: Consider initial condition f = f(x ) = ( 1) ɛ, 0 < ɛ 1 so that f = ɛ After one time step U 1 = [λ( 1) 1 + (1 2λ)( 1) + λ( 1) +1 ]ɛ = (1 4λ)( 1) ɛ and after n time steps U n = (1 4λ) n ( 1) ɛ If λ > 1/2, then U n = 1 4λ n ɛ as n 7 / 46
8 Maximum principle Remark: The condition λ 1 2 depends on the space step h as t h2 2µ implies that the allowed time step t = t = O ( h 2) The scheme is said to be conditionally stable. 8 / 46
9 Truncation error and consistency If u(x, t) is exact solution then with u n = u(x, t n ), the local truncation error τ n := D t + u n µd x + Dx u n Using Taylor s formula, for some x (x 1, x +1 ) and t n I n = (t n, t n+1 ) τ n = 1 2 tu tt(x, t n ) 1 12 µh2 u xxxx ( x, t n ) Using the PDE, we get u tt = µ 2 u xxxx and from the exact solution, we have u xxxx C f xxxx C. The truncation error is bounded as τ n C t max t I n u tt (t) C + Ch 2 u xxxx (t n ) C Ch 2 max t I n u xxxx (t) C, since λ 1/2 Ch 2 f xxxx C 9 / 46
10 Error estimate Let U n and u n be the numerical and exact solutions, and let λ 1/2. Then there is a constant C such that Proof: Set e n = U n u n. Then U n u n Ct n h 2 f xxxx C, t n 0 D + t e n µd + x D x e n = τ n and hence Iterating over time e n+1 = (E h e n ) t τ n n 1 e n = (Eh n e 0 ) t n 1 = t l=0 l=0 (E n 1 l (E n 1 l h τ l ) h τ l ) since e 0 = U 0 u 0 = 0 10 / 46
11 Using the stability estimate and truncation error estimate n 1 e n t E n 1 l h τ l l=0 n 1 t τ l, l=0 C(n t)h 2 f xxxx C E n 1 l h τ l τ l Ch 2 f xxxx C = Ct n h 2 f xxxx C Remark: The FTCS scheme has first order accuracy in time and second order accuracy in space τ n = O ( t + h 2). Since t = O ( h 2) under the stability condition λ 1/2, the scheme has second order accuracy with respect to the mesh width h. If h is halved, then the error reduces by 1/4. But number of time steps increases by a factor of / 46
12 Symbol of E h FTCS scheme General finite difference scheme U n+1 = λu n 1 + (1 2λ)U n + λu n +1 U n+1 = p a p (λ)u n p How does scheme modify one Fourier mode: U n = gn e iξ g = amplitude g n+1 (ξ) = Ẽh(ξ)g n (ξ) ξ is called the wave number. Symbol or characteristic polynomial of E h Ẽ h (ξ) = p a p e ipξ, ξ R Useful in stability analysis: Fourier stability 12 / 46
13 Necessary condition for maximum stability A necessary condition for stability of the operator E h with respect to the discrete maximum norm is that Ẽh(ξ) 1, ξ R Proof: Assume that E h is stable in maximum norm and that Ẽh(ξ 0 ) > 1 for some ξ 0 R. Then with initial condition f = e iξ0 ɛ, the numerical solution after one time step is U 1 = p ( a p e i( p)ξ0 ɛ = a p e )e ipξ0 iξ0 ɛ = Ẽh(ξ 0 )f p By iterating in time, we get U n = [Ẽh(ξ 0 )] n f and U n = Ẽh(ξ 0 ) n ɛ as n Since f = ɛ, this contradicts the stability of the scheme and hence the theorem is proved. 13 / 46
14 FTCS scheme Symbol U n+1 = λu n 1 + (1 2λ)U n + λu n +1 Ẽ h (ξ) = λe iξ + (1 2λ)e 0 + λe iξ = 1 2λ + 2λ cos(ξ) Since cos(ξ) [ 1, +1], the maximum value is attained for cos(ξ) = 1 max Ẽh(ξ) = 1 4λ 1 = λ 1 ξ 2 which is the condition previously derived for maximum stability. 14 / 46
15 Discrete Fourier transform For any grid function V = {V } R define discrete l 2 norm V 2,h = h Space of bounded grid functions Discrete Fourier transform =+ = V 2 1/2 l 2,h = {V R : V 2,h < } V (ξ) = h + = V e iξ Inverse transform V = 1 +π V (ξ)e iξ dξ 2πh π 15 / 46
16 Discrete Fourier transform Parseval relation V 2 2,h = 1 2πh +π π V (ξ) 2 dξ = 1 2π +π/h π/h V (ξh) 2 dξ Definition: Stability in l 2,h The discrete solution operator E h is said to be stable in the norm 2,h if E h V 2,h V 2,h, V l 2,h 16 / 46
17 Von Neumann stability condition A necessary and sufficient condition for E h to be stable in l 2,h is given by Ẽh(ξ) 1, ξ R Proof: For any V l 2,h, Ê h V (ξ) = h (E h V ) e iξ = h = a p V p e iξ p = a p e ipξ h V p e i( p)ξ p = Ẽh(ξ) V (ξ) and iterating in time Ê n h V (ξ) = [Ẽh(ξ)] n V (ξ) 17 / 46
18 Using Parseval relation, stability of E h in l 2,h is equivalent to +π Ê n h V (ξ) 2 dξ = +π Ẽh(ξ) 2n V (ξ) 2 dξ +π π π π V (ξ) 2 dξ which holds if and only if Ẽh(ξ) n 1, n 0, ξ R Remark: If the problem is in the time interval (0, T ), then a less restrictive notion of stability is given by the condition E n h V C T V, 0 n T/ t If the finite difference scheme satisfies E h V (1 + K t) V, K 0 then this implies Eh n V (1 + K t) n V e Kn t V e KT V 18 / 46
19 Fourier transform in space Define numerical scheme at all x R, not ust at the mesh points U n+1 (x) = (E h U n )(x) = p a p (λ)u n (x x p ) Advantage: All the U n lie in same function space, L 2 (R) or C(R), independently of h. Define usual L 2 norm v = + v(x) 2 dx 1/2 Fourier transform + v(ξ) = v(x)e ixξ dx Then Parseval relation v 2 = 1 2π v 2 Ê h v(ξ) = p a p ( v( ph)(ξ) = p a p e iphξ ) v(ξ) = Ẽh(hξ) v(ξ) 19 / 46
20 For the numerical scheme, we thus find U n = 1 2π [ Ẽ h (hξ)] n Û 0 1 Stability in L 2 norm holds if 2π sup ξ R sup Ẽh(hξ) 1 ξ R Ẽh(hξ) n Û 0 = sup Ẽh(hξ) n U 0 ξ R Accuracy of order r The finite difference scheme E h is accurate of order r if Ẽ h (ξ) = e λξ2 + O ( ξ r+2), ξ 0 Example: FTCS scheme Ẽ h (ξ) = 1 2λ + 2λ cos(ξ) = 1 λξ λξ4 + O ( ξ 6) ( 1 = e λξ λ 1 ) 2 λ2 ξ 4 + O ( ξ 6) FTCS is accurate of order 2; For the choice λ = 1 6 it is accurate of order / 46
21 Convergence in L 2 Assume that the scheme E h with λ = t h = constant and is accurate of order r 2 and stable in L 2. Then U n u n Ct n h r f r+2, t n 0 Proof: Since Ẽh(ξ) is bounded in R, we have from accuracy property Ẽh(ξ) e λξ2 C ξ r+2, ξ R By stability and a n b n = (a b)(a n 1 + a n 2 b b n ) it follows that n 1 [Ẽh(ξ)] n e nλξ2 (Ẽh(ξ) e λξ 2 ) [Ẽh(ξ)] n 1 e λξ2 Cn ξ r+2 =0 Take Fourier transform of heat equation in x variable dû dt (ξ, t) = ξ2 û(ξ, t), t > 0, û(ξ, 0) = f(ξ) 21 / 46
22 so that û(ξ, t) = e ξ2 t f(ξ) Fourier transform of error in numerical solution is From Parseval relation But U n u n 2 = 1 2π and using the fact that (Û n û n )(ξ) = ([Ẽh(hξ)] n e n tξ2 ) f(ξ) Û n û n 2 = 1 2π ([Ẽh(hξ)] n e n tξ2 ) f(ξ) [Ẽh(hξ)] n e n tξ2 Cnh r+2 ξ r+2 2 df = iξ f(ξ), dx d r f dx = ( iξ)r f(ξ), λ = t h 2 = const 22 / 46
23 we get U n u n 1 Cnh r+2 ξ r+2 f(ξ) 2π C(n t)h r f (r+2), h 2 = O ( t) = Ct n h r f r+2 23 / 46
24 Multi-level scheme Higher order accuracy in time D 0 t U n = µd + x D x U n or U n+1 U n 1 = µ U n 1 2U n + U n +1 h 2 2 t Formally second order in t and x, but unstable Dufort-Frankel scheme: stable for any λ U n+1 U n 1 2 t ( U 1 n 2 = µ U n+1 +U n 1 2 h 2 ) + U n +1 Note: You need to get U 1 by some other scheme. Second order accurate (need t/h 0, ok if t/h 2 =constant) τh n = O ( t 2) + O ( h 2) ( ) + t2 t 4 h 2 u tt(x, t n ) + O h 2 24 / 46
25 Initial-boundary value problem u t = µu xx, (x, t) (0, 1) R + u(x, 0) = f(x) u(0, t) = 0 u(1, t) = 0 25 / 46
26 IBVP: Forward Euler scheme Consider a uniform grid of M + 1 points, indexed from 0 to M. D + t U n = µd + x D x U n, = 1,..., M 1 or or U n+1 U n t = U n 1 2U n + U n +1 h 2 U n+1 = λu n 1 + (1 2λ)U n + λu n +1, = 1,..., M 1 U n+1 0 = U n+1 M = 0 Maximum norm U,h = max U 0 M 26 / 46
27 IBVP: Forward Euler scheme Maximum stability The forward Euler scheme for the IBVP is stable in maximum norm iff λ 1/2. λ 1/2 is necessary: Consider initial condition U 0 = f = ( 1) sin(πh), = 0, 1,..., M so that f,h = 1. Then U n = [1 2λ 2λ cos(πh)] n f, = 0, 1,..., M If λ > 1/2, then for sufficiently small h 1 2λ 2λ cos(πh) γ > 1 which implies U n,h γ n f,h = γ n, h 0 27 / 46
28 IBVP: Forward Euler scheme Local truncation error: (for any fixed λ) τ n Ch 2 max t I n u xxxx (t) C, I n = (t n, t n+1 ) Error estimate If λ 1/2, then the FTCS scheme for IBVP has error estimate given by U n u n,h Ct n h 2 max t t n u xxxx (t) C 28 / 46
29 Numerical example: FTCS (forward Euler) PDE u t = µu xx x (0, 1) Initial condition u(x, 0) = f(x) = sin(πx) Exact solution u(x, t) = e µπ2t sin(πx) Try with λ = 0.5, 0.6 Matlab code: heat 1d fe.m 29 / 46
30 Implicit scheme for IBVP: Backward Euler Explicit scheme (Forward Euler/FTCS) D + t U n = µd + x D x U n Implicit scheme (Backward Euler/BTCS) D + t U n = µd + x D x U n+1 or U n+1 U n t = µ U n+1 1 n+1 2U + U n+1 +1 h 2 λu n (1 + 2λ)U n+1 λu n+1 +1 = U n, = 1,..., M 1 U n+1 0 = U n+1 M = 0 In matrix notation where U = [U 1, U 2,..., U M 1 ] and BU n+1 = U n 30 / 46
31 where Implicit Ū n+1 and scheme Ū n are now thought of as vectors with M 1components corresponding to the for interior IBVP: mesh-points Backward and B is the Euler diagonally dominant, symmetric, tridiagonal matrix 1+2λ λ λ 1+2λ λ... B = λ 1+2λ λ λ 1+2λ Clearly the system (9.22) may easily be solved for Ū n+1. Maximum stability The backward Euler scheme is stable in maximum norm for any value of λ. The scheme is unconditionally stable. Proof: At time level n + 1, let the maximum be achieved at grid point = 0. Then U n+1 0 = 1 [ λu n+1 n λ λu + U n ] / 46
32 Implicit scheme for IBVP: Backward Euler and U n+1,h = U n λ 2λ 1 + 2λ [ λ U n+1 n λ U + U n ] U n+1,h λ U n,h which implies that U n+1,h U n,h. Remark: If we write BU n+1 = U as U n+1 = B 1 U n =: E U n then maximum norm stability implies that E U,h U,h = E,h 1 32 / 46
33 Implicit scheme for IBVP: Backward Euler Error estimate The backward Euler scheme has the error estimate U n u n,h Ct n ( t + h 2 ) max t t n u xxxx (t) C Proof: Define the error e n = U n u n. U n+1 U n t u n t u n+1 = µ U n+1 1 The local truncation error can be bounded as n+1 2U + U n+1 +1 h 2 = µ un+1 1 2un+1 + u n+1 +1 h 2 + τ n τ n,h C( t + h 2 ) max t I n 1 u(t) C 4, I n 1 = (t n 1, t n ) 33 / 46
34 Implicit scheme for IBVP: Backward Euler Then the error satisfies e n+1 e n t = µ en+1 1 2en+1 + e n+1 +1 h 2 τ n In terms of matrix B and discrete evolution operator E Iterating over time Be n+1 = e n t τ n = e n+1 = E e n te τ n n 1 e n = E e n 0 t E n l τ l, n 1 l=0 n 1 = t E n l τ l since e 0 = 0 l=0 34 / 46
35 Implicit scheme for IBVP: Backward Euler Therefore n 1 e n,h t E n l τ l,h l=0 n 1 t τ l,h, l=0 (n t) max 0 l n 1 τ l,h t n C( t + h 2 ) max t t n u(t) C 4 from stability of E Remark: The Backward Euler scheme is first order accurate in time and second order accurate in space. 35 / 46
36 Implicit scheme for IBVP: Crank-Nicholson D + t U n = µd + x D x ( ) U n + U n+1 2 λ 2 U n (1 + λ)u n+1 λ 2 U n+1 +1 = λ 2 U n 1 + (1 λ)u n + λ 2 U n +1 U n+1 0 = U n+1 M = 0 In matrix notation where BU n+1 = AU n 36 / 46
37 BŪ n+1 = AŪ n, Implicit scheme for IBVP: Crank-Nicholson both A and B are symmetric tridiagonal matrices, with9.2 B diagoant: The Mixed Initial-Boundary Value 1+λ 1 2 λ 0 and λ λ 1+λ 1 2 λ λ B = , λ 1 λ 1 2 λ.... A = λ 1+λ 1 2 λ λ 1+λ 2 λ 1 λ 1 2 λ λ 1 λ With obvious notation we also have U n+1 = B 1 AU n =: E U n B kh U n+1 = A kh U n, or Maximum stability U n+1 = B 1 kh A khu n = E k U n, The Crank-Nicholson scheme where, is stable similarly into maximum the above, norm for λ 1 Proof: Similar to backward Euler scheme B 1 kh V,h V,h. Remark: Backward Euler scheme The same wasapproach unconditionally to stability stable as forinthe maximum backward Euler norm. for λ 1, since the coefficients on the right are then non-neg (1 + λ) U n+1,h λ U n+1,h + U n,h, 37 / 46
38 Implicit scheme for IBVP: Crank-Nicholson Remark: If λ > 1, then U n+1,h (2λ 1) U n,h But since 2λ 1 > 1, this does not yield stability in maximum norm. Error estimate For λ 1, the Crank-Nicholson scheme has the error estimate U n u n,h = O ( t 2 + h 2) Remark: The condition λ 1 is too restrictive. This can be relaxed if we use a weaker norm for stability. 38 / 46
39 IBVP: Fourier stability Vector V = (V 0, V 1,..., V M ) R M+1 Inner product and corresponding norm Space of vectors (V, W ) h := h M V W =0 V 2,h = (V, V ) h l 0 2,h = {V R M+1 : V 2,h <, V 0 = V M = 0} Orthonormal basis vectors for l 0 2,h : φ p R M+1, p = 1,..., M 1 φ p, = 2 sin(πph), = 0, 1,..., M (φ p, φ q ) h = δ pq φ p are eigenfunctions of finite difference operator D + x D x D + x D x φ p, = 2 h 2 [1 cos(πph)]φ p,, = 1, 2,..., M 1 39 / 46
40 IBVP: Fourier stability Let f l2,h 0 be initial data. Then f = M 1 p=1 ˆf p φ p, ˆfp = (f, φ p ) h Parseval relation f 2,h = ˆf 2 Forward Euler method (FTCS): Consider one time step U 1 = f + µ td + x D x f = M 1 p=1 ˆf p [1 2λ(1 cos(πph))]φ p,, = 1, 2,..., M 1 and U 1 0 = U 1 M = 0 40 / 46
41 IBVP: Fourier stability or more generally U n = M 1 p=1 ˆf p [Ẽh(πph)] n φ p,, = 0, 1,..., M where Ẽh(ξ) is the symbol of the operator E h Ẽ h (ξ) = 1 2λ[1 cos(ξ)] By orthogonality of {φ p } and Parseval s relation U n 2,h = ( M 1 p=1 ) 1/2 ˆf p 2 [Ẽh(πph)] 2n max Ẽh(πph) n f p 2,h Check that Ẽh(ξ) 1 for ξ [0, π] iff λ 1/2. In that case, we have stability U n 2,h f 2,h, iff λ 1/2 41 / 46
42 IBVP: Fourier stability This condition is identical to the condition for maximum norm stability of FTCS scheme. Backward Euler method Ẽ h (ξ) = λ[1 cos(ξ)] and 0 Ẽh(ξ) 1 for all ξ [0, π] and any λ > 0. Hence we have unconditional stability. Crank-Nicholson method Ẽ h (ξ) = 1 λ[1 cos(ξ)] 1 + λ[1 cos(ξ)] and Ẽh(ξ) 1 for ξ [0, π] and for any λ > 0. Thus the CN scheme is unconditionally stable in l 0 2,h. 42 / 46
43 IBVP: Fourier stability Convergence of CN scheme in l 0 2,h Let U n and u n be numerical and exact solutions. Then for any λ > 0 U n u n 2,h Ct n (h 2 + t 2 ) max t t n u(t) C 6 Proof: The local truncation error ( ) u n τ n = D t + u n µd x + Dx + u n+1 2 = [ D + t u n u t (x, t n+ 1 2 ) ] µd + x D x [ ] µ D x + Dx u n+ 1 2 u xx (x, t n+ 1 ) 2 and using Taylor formula we get [ u n + u n+1 2 τ n 2,h C( t 2 + h 2 ) max t I n u(t) C 6 u n+ 1 2 ] 43 / 46
44 IBVP: Fourier stability The error e = U n u n satisfies the difference equation e n+1 = E e n t B 1 τ n or n 1 e n = t E n 1 l B 1 τ l l=0 The result follows from the stability of E and the boundedness of B 1, i.e., E V 2,h V 2,h and B 1 V 2,h V 2,h 44 / 46
45 θ-scheme Symbol of the scheme D + t U n = µd + x D x [(1 θ)u n + θu n+1 ] θ = 0 Forward Euler Explicit O ( t + h 2) θ = 1 Backward Euler Implicit O ( t + h 2) θ = 1/2 Crank-Nicholson Implicit O ( t 2 + h 2) Ẽ h (ξ) = 1 2(1 θ)λ(1 cos ξ) 1 + 2θλ(1 cos ξ) For 0 θ 1, we have Ẽh(ξ) 1 for all ξ. Stability requires that 1 4(1 θ)λ min Ẽ h (ξ) = 1 = (1 2θ)λ 1 ξ 1 + 4θλ 2 θ < 1/2: stable in l 0 2,h only if λ < 1 2(1 2θ) θ 1/2: unconditionally stable in l 0 2,h 45 / 46
46 Numerical example Initial and boundary condition same as in (29) Matlab codes: Backward euler scheme: heat 1d be.m Crank-Nicholson scheme: heat 1d cn.m 46 / 46
Numerical 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 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 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 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 informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I 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 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 informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
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 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 informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationImplicit Scheme for the Heat Equation
Implicit Scheme for the Heat Equation Implicit scheme for the one-dimensional heat equation Once again we consider the one-dimensional heat equation where we seek a u(x, t) satisfying u t = νu xx + f(x,
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 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 informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
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 informationIntroduction to PDEs and Numerical Methods: Exam 1
Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently
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 informationFinite difference approximation
Finite difference approximation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen January
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationSmoothing Effects for Linear Partial Differential Equations
Smoothing Effects for Linear Partial Differential Equations Derek L. Smith SIAM Seminar - Winter 2015 University of California, Santa Barbara January 21, 2015 Table of Contents Preliminaries Smoothing
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 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 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 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 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 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 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 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 informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
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 informationMATH3383. Quantum Mechanics. Appendix D: Hermite Equation; Orthogonal Polynomials
MATH3383. Quantum Mechanics. Appendix D: Hermite Equation; Orthogonal Polynomials. Hermite Equation In the study of the eigenvalue problem of the Hamiltonian for the quantum harmonic oscillator we have
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationu(0) = u 0, u(1) = u 1. To prove what we want we introduce a new function, where c = sup x [0,1] a(x) and ɛ 0:
6. Maximum Principles Goal: gives properties of a solution of a PDE without solving it. For the two-point boundary problem we shall show that the extreme values of the solution are attained on the boundary.
More informationChapter 5: Bases in Hilbert Spaces
Orthogonal bases, general theory The Fourier basis in L 2 (T) Applications of Fourier series Chapter 5: Bases in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University
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 information4 Riesz Kernels. Since the functions k i (ξ) = ξ i. are bounded functions it is clear that R
4 Riesz Kernels. A natural generalization of the Hilbert transform to higher dimension is mutiplication of the Fourier Transform by homogeneous functions of degree 0, the simplest ones being R i f(ξ) =
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 informationPARTIAL DIFFERENTIAL EQUATIONS. Lecturer: D.M.A. Stuart MT 2007
PARTIAL DIFFERENTIAL EQUATIONS Lecturer: D.M.A. Stuart MT 2007 In addition to the sets of lecture notes written by previous lecturers ([1, 2]) the books [4, 7] are very good for the PDE topics in the course.
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 informationAn Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations.
An Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations by Tong Chen A thesis submitted in conformity with the requirements
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
More informationLecture 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 informationLecture 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 informationMATH 220 solution to homework 5
MATH 220 solution to homework 5 Problem. (i Define E(t = k(t + p(t = then E (t = 2 = 2 = 2 u t u tt + u x u xt dx u 2 t + u 2 x dx, u t u xx + u x u xt dx x [u tu x ] dx. Because f and g are compactly
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationChapter 4. MAC Scheme. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes. u = 0
Chapter 4. MAC Scheme 4.1. MAC Scheme and the Staggered Grid. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes equation in the velocity-pressure formulation: (4.1.1)
More informationDiffusion on the half-line. The Dirichlet problem
Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection
More informationIntroduction to PDEs and Numerical Methods Tutorial 4. Finite difference methods stability, concsistency, convergence
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 4. Finite difference methods stability, concsistency, convergence Dr. Noemi Friedman,
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 informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationNumerical PDE. DYKwak. September 6, 2007
Numerical PDE DYKwak September 6, 2007 ii Contents 1 Finite Difference Method 7 1.1 2nd order linear p.d.e. in two variables................ 7 1.2 Finite difference operator....................... 8 1.2.1
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 informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationTMA4212 Numerical solution of partial differential equations with finite difference methods
TMA4 Numerical solution of partial differential equations with finite difference methods Brynjulf Owren January 3, 07 Translated and amended by E. Celledoni Preface The preparation of this note started
More informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
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 informationConvergence of Summation-by-Parts Finite Difference Methods for the Wave Equation
J Sci Comput (2) 71:219 245 DOI 1./s1915-16-297-3 Convergence of Summation-by-Parts Finite Difference Methods for the Wave Equation Siyang Wang 1 Gunilla Kreiss 1 Received: 15 January 216 / Revised: 16
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
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 informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
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 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 informationFOURIER TRANSFORMS. 1. Fourier series 1.1. The trigonometric system. The sequence of functions
FOURIER TRANSFORMS. Fourier series.. The trigonometric system. The sequence of functions, cos x, sin x,..., cos nx, sin nx,... is called the trigonometric system. These functions have period π. The trigonometric
More informationDomain decomposition schemes with high-order accuracy and unconditional stability
Domain decomposition schemes with high-order accuracy and unconditional stability Wenrui Hao Shaohong Zhu March 7, 0 Abstract Parallel finite difference schemes with high-order accuracy and unconditional
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationCharacterization of Hardy spaces for certain Dunkl operators and multidimensional Bessel operators
Characterization of for certain Dunkl operators and multidimensional Bessel operators Instytut Matematyczny, Uniwersytet Wrocławski, Poland Probability & Analysis Będlewo, 4-8.05.2015. Hardy space H 1
More informationJuan Vicente Gutiérrez Santacreu Rafael Rodríguez Galván. Departamento de Matemática Aplicada I Universidad de Sevilla
Doc-Course: Partial Differential Equations: Analysis, Numerics and Control Research Unit 3: Numerical Methods for PDEs Part I: Finite Element Method: Elliptic and Parabolic Equations Juan Vicente Gutiérrez
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 informationPartial Differential Equations (PDEs) and the Finite Difference Method (FDM). An introduction
Page of 8 Partial Differential Equations (PDEs) and the Finite Difference Method (FDM). An introduction FILE:Chap 3 Partial Differential Equations-V6. Original: May 7, 05 Revised: Dec 9, 06, Feb 0, 07,
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationSelected HW Solutions
Selected HW Solutions HW1 1 & See web page notes Derivative Approximations. For example: df f i+1 f i 1 = dx h i 1 f i + hf i + h h f i + h3 6 f i + f i + h 6 f i + 3 a realmax 17 1.7014 10 38 b realmin
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationChapter 7: Bounded Operators in Hilbert Spaces
Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84
More informationNUMERICAL ANALYSIS STUDY GUIDE. Topic # Rootfinding
NUMERICAL ANALYSIS STUDY GUIDE BRYON ARAGAM These are notes compiled during the Summer of 2010 while studying for the UCLA Numerical Analysis qualifying exam. Not all topics listed below are covered, which
More informationNumerical methods for complex systems. S. V. Gurevich
Numerical methods for complex systems S. V. Gurevich December 16, 2008 2 Part I Partial Differential Equations 3 Chapter 1 Intorduction 1.1 Definition, Notation and Classification A differential equation
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
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 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 informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
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 informationConvergence and Error Bound Analysis for the Space-Time CESE Method
Convergence and Error Bound Analysis for the Space-Time CESE Method Daoqi Yang, 1 Shengtao Yu, Jennifer Zhao 3 1 Department of Mathematics Wayne State University Detroit, MI 480 Department of Mechanics
More informationMATH 126 FINAL EXAM. Name:
MATH 126 FINAL EXAM Name: Exam policies: Closed book, closed notes, no external resources, individual work. Please write your name on the exam and on each page you detach. Unless stated otherwise, you
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationarxiv: v1 [physics.comp-ph] 22 Feb 2013
Numerical Methods and Causality in Physics Muhammad Adeel Ajaib 1 University of Delaware, Newark, DE 19716, USA arxiv:1302.5601v1 [physics.comp-ph] 22 Feb 2013 Abstract We discuss physical implications
More informationMATH-UA 263 Partial Differential Equations Recitation Summary
MATH-UA 263 Partial Differential Equations Recitation Summary Yuanxun (Bill) Bao Office Hour: Wednesday 2-4pm, WWH 1003 Email: yxb201@nyu.edu 1 February 2, 2018 Topics: verifying solution to a PDE, dispersion
More informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationIntroduction to Mathematical Modelling: Ordinary Differential Equations and Heat Equations
Introduction to Mathematical Modelling: Ordinary Differential Equations and Heat Equations Knut Andreas Lie SINTEF ICT, Dept Applied Mathematics INF2340 / Spring 2005 p Overview 1 Ordinary differential
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationA second-order asymptotic-preserving and positive-preserving discontinuous. Galerkin scheme for the Kerr-Debye model. Abstract
A second-order asymptotic-preserving and positive-preserving discontinuous Galerkin scheme for the Kerr-Debye model Juntao Huang 1 and Chi-Wang Shu Abstract In this paper, we develop a second-order asymptotic-preserving
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationNotes on Numerical Analysis
Notes on Numerical Analysis Alejandro Cantarero This set of notes covers topics that most commonly show up on the Numerical Analysis qualifying exam in the Mathematics department at UCLA. Each section
More information30 crete maximum principle, which all imply the bound-preserving property. But most
3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More information