Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
|
|
- Richard Bennett
- 5 years ago
- Views:
Transcription
1 Opuscula Matematica Vol. 26 No Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating multidimensional parabolic semilinear problems in a very efficient way. Te metod reaces te fourt order in time and it can be combined wit standard spatial discretizations of any order to obtain unconditinally convergent numerical algoritms. Te main teoretical results wic guarantee tis property are explained ere, as well as te metod caracteristics wic guarantee a very strong reduction of computational cost in comparison wit classical discretization metods. Keywords: fractional step metods, non-linear parabolic problems, convergence. Matematics Subject Classification: 65M6, 65M INTRODUCTION In tis paper we deal wit te construction of numerical metods wic efficiently integrate te following non-linear parabolic problems: find y(t) : [t, T ] H, a solution of { y (t) = L(t)y(t) + f(t) + g(t, y(t)), y(t ) = y. Here H is a Hilbert space of functions defined on a certain domain Ω R n, L(t) is a linear second order elliptic differential operator, f(t) is te source term, g(t, y(t)) is a non-linear function wic meets some additional requirements and y y ( x) is te initial condition. Te numerical integration of tis problem can be viewed as a double process of discretization wic adequately combines a time integrator (typically an adapted ODE Solver) wit a discretization of te spatial variables. It is well known tat te use of classical implicit scemes for te numerical integration of tis type of problems poses mainly two problems. On te one and, te convergence of te numerical sceme is usually obtained by imposing excessively restrictive stability requirements wic force (1) 47
2 48 Blanca Bujanda, Juan Carlos Jorge us to use low order metods (like te implicit Euler rule) or complicated fully implicit metods for reacing iger orders of convergence in time. On te oter and, we must use slow iterative metods to resolve te internal stage equations because tey are defined as te solution of large non-linear systems. Also, in tis type of scemes, te convergence is deduced by imposing strong stability requirements (for example, B-stability). In [5] a semiexplicit metod, wic is used to determinate te contribution of te linear term, is combined in an additive way wit an explicit metod used for te contribution of te non-linear term. Tus, linear systems only appear in te computation of te internal stages. In [7], splitting metods are applied in te numerical integration of some non-linear advention-diffusion-reaction problems; ere te advention terms are integrated wit an explicit metod. In order to avoid tese problems, in [3] a new class of linearly implicit metods is developed; tey are built by combining a standard spatial discretization of Finite Differences or Finite Elements type wit a suitable time discretization, wic is deduced from a Fractional Step Runge Kutta (FSRK) metod to deal wit te linear terms and an explicit Runge Kutta metod to define te contribution of te non-linear term. By applying tis new linearly implicit metods (see [1, 2, 3]) we obtain numerical algoritms wic are convergent under a single requirement of te linear absolute stability type. It is well known (see [8]) tat te Fractional Step metods are very efficient in integrating multidimensional linear parabolic problems if we suitably decompose te linear elliptic operator and te source term; te advantages of tese metods manifest temselves also in tese new linearly implicit metods. Tus in te resolution of te internal stages only we must solve linear systems wit very simple matrices and it is not necessary to apply classical iterative metods to solve tem. For example, if te elliptic linear operator L(t) is n d i(t) 2 y x 2 i + n v i(t) y x i + n k i(t), and we realize a spatial discretization by applying classical Finite Differences, te matrices wic appear at eac stage are tridiagonal; consequently, te computational cost of te final metod is of te same order as tat of an explicit metod. 2. NUMERICAL SCHEME By applying firstly a spatial discretization process of type Finite Difference, Finite Elements,... to problem (1) we obtain a family of Initial Value Problems wic depend on te parameter (, ] used to discretize in space as follows: Find Y (t) : [t, T ] V, a solution of { Y (t) = L (t)y (t) + f (t) + g (t, Y (t)), Y (t ) = y. (2) For eac it is common to consider a finite dimensional space V ; suc space is te space of discrete functions on a mes in Finite Differences, te subspace of H of piecewise polynomial functions in a classical Finite Element discretization. We
3 New efficient time integrators for non-linear parabolic problems 49 suppose tat suc spatial semidiscretization process is uniformly convergent of order q, i.e., for sufficiently smoot functions y(t) tere is π y(t) Y (t) C q, t [t, T ]. (3) Here, te operator π denotes te restriction to te mes nodes if Finite Differences used and if we use te standard Finite Elements, ten π are suitable projections in te space V. In order to integrate efficiently tis problem by applying te linearly implicit numerical metods wic are first appeared in [3], we must adequately decompose te spatial discretization of te elliptic operator and te source term in n addends as follows: L (t) = n L i(t) and f (t) = n f i(t); te FSRK metod will act on tese terms. Te contribution of te non-linear term g (t, U (t)) is given by te explicit RK metod. From tis process, we obtain te following sceme: s Y m+1 = Y m +τ were i Y m,i = Y m+τ b ki i ( ) L ki(t m,i )Y m,i + f ki(t m,i ) + τ s ( ) i 1 a kj ij L kj(t m,j )Y m,j + f kj(t m,j ) + τ b n+1 i g (t m,i, Y m,i ), (4) a n+1 ij g (t m,j, Y m,j ). In tis algoritm, we denote wit Y m te numerical approximations to te exact solution y(t m ) at te times t m = t + m τ (τ is te time step), U = U (t ), k i, k j {1,, n}, n is te number of te levels of te metod, t m,i = t + (m + c i ) τ and te intermediate approximations U m,i for i = 1,, s are te internal stages of te metod. By setting some coefficients to zero in te latest formula as follows A k = (a k ij) were a k ij = b k = (b k j ) were b k j = a n+1 ij if k = n + 1 y i > j, a kj ij if k = k j and i j, oterwise, b n+1 j if k = n + 1, b kj j if k = k j, oterwise, we can organize tese coefficients in a table wic can be considered as an extension of te tables introduced by Butcer for standard RK metods, in te following way C e A 1 A 2... A n A n+1 (b 1 ) T (b 2 ) T... (b n ) T (b n+1 ) T were C = diag(c 1,, c s ) and e = (1,, 1).
4 41 Blanca Bujanda, Juan Carlos Jorge We will use (C, A i, b i ) n+1 to refer ro tis metod sortly and (C, Ai, b i ) to denote every standard RK metod involved in (4). To study te convergence of numerical sceme (4) we introduce, as usual, te global error at te time t = t m by E m = π y(t m ) Y m ; tus we say tat (4) is uniformly convergent of order p in time and of order q in space, if E m satisfies E m C( q + τ p ). (5) To get tis bound, we begin by studying te stability of numerical sceme (4); tus we rewrite it in a more compact way by using te following tensorial notation: given M (m ij ) R s s and v (v i ) R s, we denote M (m ij I V ) V s s and v (v i I V ) V s; we group te evaluations of te source terms f i(t) and of te linear operators L i (t) for i = 1,, n and for all m = 1, 2,, Fi m = (f i(t m,1 ),, f i (t m,s )) T V s and ˆL m i = diag (L i(t m,1 ),, L i (t m,s )) ( ) T L(V, V ). We also group te contribution of te stages Ỹ m = Y m,1,, Y m,s ( T V s and te nonlinear term Ĝm (Ỹ m) = g (t m,1, Y m,1 ),, g (t m,s, Y m,s )) V s. In tis way we obtain n n (Ī τ A i ˆL m m i )Ỹ = ē Y m + τ A i Fi m + τ A n+1 Ĝm (Ỹ m ), Y m+1 = Y m + τ n (b i ) T (ˆLm i Ỹ m + F m i ) + τ (b n+1 ) T Ĝ m (Ỹ m ). Wen te operator ( Ī n Aj ˆL j) m is invertible 1), te numerical solution can be written as were Y m+1 = R ( τ ˆL m 1,, τ ˆL m n) Y m + S ( τ ˆL m 1,, τ ˆL m n, τf m 1,, τf m n + T (τ ˆL m 1,, τ ˆL m n, τ Ĝm (Ỹ m )), (6) R(τ ˆL m 1,, τ ˆL m n) = Ī + n (b i ) T τ ˆL n 1ē, m i (Ī A j τ ˆL j) m (7) ) and S ( τ ˆL m 1,, τ ˆL m n, τ F1, m, τ Fn) m = (8) n ( = τ (b i ) T Fi m + ˆL m (Ī n i τ A j ˆL m ) 1 ( n j τ A k Fk m ) ) T (τ ˆL m 1,, τ ˆL m n, τ Ĝm (Ỹ m )) = (9) n = τ (b i ) T (Ī n ˆLm i A j ˆL m 1τA n+1 j) Ĝ m (Ỹ m ) + τ (b n+1 ) T Ĝ m (Ỹ m ). k=1 1) In [4] it is proved tat tis operator is invertible and tat its inverse operator is bounded independently of and τ.
5 New efficient time integrators for non-linear parabolic problems 411 Note tat first term (7) is te transition operator if we consider linear omogeneous problems; te contribution of te source term appears in second term (8) and, finally, te contribution of te nonlinear term g (t, Y (t)) is defined by tird addend (9). To study te stability of te numerical sceme we must bound tese tree addends (see [1]). Te contribution of te linear terms, given by (7), can be bounded by applying te following result (see [4]): Teorem 2.1. Let an A-stable FSRK metod be given suc tat all teir stages are implicit (i.e., n k=1 ak ii, for i = 1,, s) or its first stage only is explicit (i.e., n k=1 ak ii, for i = 2,, s and aki 11 = for all i = 1,, s) and it also satisfies (,,, 1) T A i = (b i ) T and a k1 ss, and let {L i (t)} n be a linear maximal coercive system of operators suc tat: a) for eac t [, T ] te system of operators {L i (t)} n is commutative, b) tere exist n constants M i suc tat L i (t )Y L i (t)y t t C L i (t)y t, t [t, T ]. Ten tere exists a constant γ, independent of τ, suc tat R( τ ˆLm 1,, τ ˆL m n) e γτ (1) olds. In [4], te following bound for second addend (8) is also proved to old: S( τ ˆLm 1,, τ ˆL m n, τ F m 1,, τ F m n) C τ n F m i To finis te study of te stability, we must bound tird term (9); in [1], te following result is proved, rewritten for our purposes: Teorem 2.2. Let (4) be a linearly implicit metod meeting te conditions of Teorem 2.1. If we use suc metod to integrate nonlinear parabolic problem (2) in time, were te nonlinear part g (t, Y ) satisfies g (t, X ) g (t, Y ) L X Y, t [t, T ], ten:. T (τ ˆL m 1,, τ ˆL m n, τ Ĝm ( X m )) T (τ ˆL m 1,, τ ˆL m n, τ Ĝm (Ỹ m )) C τ X m Y m. (11) To study te consistency of te numerical sceme, we define te local error of te time discretization metod as e m = Y (t m ) Ŷ m, were Ŷ m is obtained wit one step of (4) taking as starting value Y m 1 = Y (t m 1 ). We say tat (4) is uniformly consistent of order p if e m C τ p+1, m. (12) In [2] a teorem wic can be rewritten for tis case as follows is proved:
6 412 Blanca Bujanda, Juan Carlos Jorge Teorem 2.3. Let us assume tat problem (2) satisfies te smootness requirements Y (p+1) i 1 (t) C, L (ϱ1) i 1 (t) L(ϱ l 1) i l 1 (t)y (ϱ l) i l (t) C, t [t, T ], l {2,, p}, (i 1,, i l ) {1,, n + 1} l, (ϱ 1,, ϱ l 1 ) {,, p 1} l 1 and ϱ l {r + 1,, p} suc tat 2 l + l ϱ k p + 2, were Y i (t) = L i(t)y (t) + f i (t), i = 1,, n + 1 wit L (n+1) (t) = g Y (t, Y (t)) and f (n+1) (t) = g Y (t, Y (t)) Y (t) + g (t, Y (t)). Let (C, A i, b i ) n+1 be a linearly implicit FSRK metod satisfying te reductions (C) k e k A j (C) k 1 e =, j = 1,, n + 1, k = 1,, r wit r = E [ ] p 1 2 togeter wit te order conditions (b i1 ) T (C) ρ1 e = 1 ρ 1 + 1, Ten (12) is true. k=1 (b i1 ) T (C) ρ1 A i2 (C) ρ2 A i l (C) ρ l e = l 1 (l j + 1) + l {2,, p}, (i 1,, i l ) {1,, n + 1} l, (ρ 1,, ρ l 1 ) {,, p 1} l 1 and ρ l {r + 1,, p 1} suc tat 1 l +, l ρ k k=j l ρ k p. To obtain te convergence of te sceme (bound (5)), we decompose te global error as follows: k=1 E m π y(t m ) Y (t m ) + Y (t m ) Ŷ m + Ŷ m Y m were we ave used te intermediate approximations Y (t m ) and Ŷ m. Te space discretization is convergent of order q, tus by applying (3) we obtain E m C q + e m + Ŷ m Y m. Now we use te fact tat te time discretization is consistent of order p (i.e. it satisfies (12)) E m C q + C τ p+1 + Ŷ m Y m. Finally, we bound te tird addend it by using (6) (9) as follows Ŷ m Y m = R ( τ ˆL m 1,, τ ˆL m n) ( π y(t m 1 ) Y (t m 1 ) ) + + T (τ ˆL m 1,, τ ˆL m n, τ Ĝm (π y(t m 1 ))) T (τ ˆL m 1,, τ ˆL m n, τ Ĝm (Y m 1 )).
7 New efficient time integrators for non-linear parabolic problems 413 Tus, by applying (1) and (11), we obtain E m C q + C τ p+1 + ( e γ τ + C τ ) E m 1. From tis result is easy to conclude tat sceme (4) satisfies (5). 3. A NEW FOURTH ORDER LINEARLY IMPLICIT METHOD Now we present te main ideas of te construction process for a linearly implicit tree-level FSRK metod; suc metod attains fourt order in time wen it is applied to te numerical integration of te semilinear problems of type (1) combined wit spatial discretizations. To te best our knowledge, tis is te only fourt order metod of tis class wic appears in te literature. Te construction process for tis kind of metods is very complicated due to two main drawbacks; te first one is related to te ig number of order conditions tat we must impose, wic makes us cope wit large and complicated non-linear systems. Te second one is related to te additional requirements tat we must impose in order to guarantee te stability of te final metod. Precisely, te order conditions tat we must impose are (see [2]): (r i) A i e = C e for i = 1, 2, 3 (α 1) (b 1 ) T e = 1 (α 2) (b 1 ) T C e = 1 2 (α 3) (b 1 ) T A 1 C e = 1 6 (α 4) (b 1 ) T C C e = 1 3 (α 5) (b 1 ) T C C C e = 1 4 (α 6) (b 1 ) T C A 1 C e = 1 8 (α 7) (b 1 ) T A 1 C C e = 1 12 (α 8) (b 1 ) T A 1 A 1 C e = 1 24 (β 1) (b 2 ) T e = 1 (β 2) (b 2 ) T C e = 1 2 (β 3) (b 2 ) T A 2 C e = 1 6 (β 4) (b 2 ) T C C e = 1 3 (β 5) (b 2 ) T C C C e = 1 4 (β 6) (b 2 ) T C A 2 C e = 1 8 (β 7) (b 2 ) T A 2 C C e = 1 12 (β 8) (b 2 ) T A 2 A 2 C e = 1 24 (γ 1) (b 3 ) T e = 1 (γ 2) (b 3 ) T C e = 1 2 (γ 3) (b 3 ) T A 3 C e = 1 6 (γ 4) (b 3 ) T C C e = 1 3 (γ 5) (b 3 ) T C C C e = 1 4 (γ 6) (b 3 ) T C A 3 C e = 1 8 (γ 7) (b 3 ) T A 3 C C e = 1 12 (γ 8) (b 3 ) T A 3 A 3 C e = 1 24 ( ij) (b i ) T A j C e = 1 6 for {i, j} {1, 2, 3} ( ij) (b i ) T C A j C e = 1 8 for {i, j} {1, 2, 3} ( ij) (b i ) T A j C C e = 1 12 for {i, j} {1, 2, 3} ( ijk ) (b i ) T A 1j A k C e = 1 24 for {i, j, k} {1, 2, 3} One may ceck tat we need at least nine stages to obtain a linearly implicit metod wic fulfils tese order conditions. Restrictions (r 1 ), (r 2 ) and (r 3 ) require te assumption tat te first stage of te metod (C, A 1, b 1 ) is explicit; we must note tat in tis way te computational cost of te final metod will be similar to te cost of an eigt stage metod wit all stages implicit. To obtain a stable metod, it is convenient to impose te restrictions (,,, 1) T A i = (b i ) T for i = 1, 2, 3 (see [3]); in order to get te numerical solution directly from te calculation of te last stage Y m+1 (= Y m,9 ) and, consequently, to reduce ligtly te computational cost of te final metod, we impose te same restriction on te tird level (,,, 1) T A 3 = (b 3 ) T.
8 414 Blanca Bujanda, Juan Carlos Jorge Te computational cost of te final metod can be reduced (wen te coefficients of te linear operator do not depend on time) by assuming tat a 1 ii = a2 jj = a for i = 3, 5, 7, 9 and for j = 2, 4, 6, 8; besides, using tese restrictions we simplify te study of te stability of te FSRK metod. Tus, te table of te deduced metod will ave te following structure: C e A 1 A 2 (b 1 ) T (b 2 ) T = = c 2 c 3 c 4 c 5 c 6 c 7 c 8 c 9 a 1 21 a a 1 31 a a 2 32 a 1 41 a 1 43 a 2 42 a a 1 51 a 1 53 a a 2 52 a 2 54 a 1 61 a 1 63 a 1 65 a 2 62 a 3 64 a a 1 71 a 1 73 a 1 75 a a 2 72 a 3 74 a 3 74 a 1 81 a 1 83 a 1 85 a 1 87 a 2 82 a 3 84 a 3 84 a b 1 1 b 1 3 b 1 5 b 1 7 a b 2 2 b 2 4 b 2 6 b 2 8 b 1 1 b 1 3 b 1 5 b 1 7 a b 2 2 b 2 4 b 2 6 b 2 8 A 3 (b 3 ) T = a 3 21 a 3 31 a 3 32 a 3 41 a 3 42 a 3 43 a 3 51 a 3 52 a 3 53 a 3 54 a 3 61 a 3 62 a 3 63 a 3 64 a 3 65 a 3 71 a 3 72 a 3 73 a 3 74 a 3 75 a 3 76 a 3 81 a 3 82 a 3 83 a 3 84 a 3 85 a 3 86 a 3 87 b 3 1 b 3 2 b 3 3 b 3 4 b 3 5 b 3 6 b 3 7 b 3 8 b 3 1 b 3 2 b 3 3 b 3 4 b 3 5 b 3 6 b 3 7 b 3 8 We begin te construction process for te metod by resolving te following order conditions: (r 1 ), (r 2 ), (α 1 ),, (α 8 ), (β 1 ),, (β 8 ), ( 12 ), ( 21 ), ( 21 ), ( 12 ), ( 112 ), ( 121 ) ( 211 ), ( 221 ), ( 212 ) and ( 122 ); wic involve te coefficients of te FSRK metods (C, A i, b i ) 2 only. Tus we obtain a four order family of FSRK metods depending on te free parameters a, c 4, c 5, c 6, c 7, a 1 65, a 1 83 and a To perform, in an efficient way, te numerical integration of te linear stiff problems, it is necessary to impose additional restrictions of type A-stability. In order to introduce tis concept in te simplest way we apply a two-level FSRK metod to te numerical integration of te following test problem y (t) = (λ 1 + λ 2 ) y(t) wit Re(λ i ), for i = 1, 2;
9 New efficient time integrators for non-linear parabolic problems 415 obtaining te recurrence: y m+1 = ( τ λ i (b i ) T ( I 2 τ λ j A j) 1 e )y m. By substituting in tis expression z i for τ λ i for i = 1, 2 we obtain a rational function of two complex variables; wic we call it te Amplification Function associated wit te FSRK metod 2 ( 2 R(z 1, z 2 ) = 1 + z i (b i ) T I z j A j) 1 e. We say tat an FSRK metod is A-stable iff R(z 1, z 2 ) 1, z i C being Re(z i ), i = 1, 2. As te FSRK metod as nine stages, wit te first one of tem explicit, and attains fourt order, te amplification function of te metod can be written as follows: R(z 1, z 2 ) = R 1 (z 1 )R 2 (z 2 ) + Rest were R i (z i ) = 1 + (1 4 a) z i (1 68 a + 12 a2 )z 2 i + ( a + 6 a2 4 a 3 )z 3 i (1 a z i ) ( a a2 4 a 3 + a 4 )z 4 i (1 a z i ) 4 are te amplification functions associated wit te eac standard RK metod (C, A i, b i ) for i = 1, 2 and were: Rest = e 41 z 4 1 z 2 + e 32 z 3 1 z e 23 z 2 1 z e 14 z 1 z e 42 z 4 1 z e 33 z 3 1 z 3 2 (1 az 1 ) 4 (1 az 2 ) e 24 z 2 1 z e 43 z 4 1 z e 34 z 3 1 z e 44 z 4 1 z 4 2 (1 az 1 ) 4 (1 az 2 ) 4, e 41 = E 4,1 1 24, e 14 = E 1,4 1 24, e 23 = E 2,3 1 12, e 32 = E 3,2 1 12, e 42 = E 4,2 4 a e 32 4 a e , e 24 = E 2,4 4 a e 23 4 a e , e 33 = E 3,3 4 a e 23 4 a e , e 43 = E 4,3 4 ae 42 1 a 2 e 41 4 a e a 2 e 32 1 a 2 e , e 34 = E 3,4 4 ae 24 1 a 2 e 14 4 a e a 2 e 23 1 a 2 e , e 44 = E 4,4 4 a e 43 4 a e 34 1 a 2 e 42 2 a 3 e a 2 e 33 4 a 3 e 32 1 a 2 e 24 4 a 3 e 23 2 a 3 e ,
10 416 Blanca Bujanda, Juan Carlos Jorge were E i1,i 2 = (b k1 ) T A k2 A kj e k (k 1,,k j) {1,2} j n l ( k)=i l, l=1,2 wit n l ( k) te number of te times tat te index l appears in te integer vector k. Given tis FSRK metod of Alternating Directions type attains fourt order, and te standard RK metod (C, A 1, b 1 ) may be reduced to a DIRK metod of five stages wit te first stage explicit, te standard RK metod (C, A 2, b 2 ) may be reduced to a four-stage SDIRK metod, we can see in [6] tat suc metods (C, A i, b i ) for i = 1, 2 are L-stable for a = Taking into account tat te strong stability and, particularly, te L-stability are very interesting properties wen we carry out te numerical integration of te parabolic problems wit irregular or incompatible data, we ave fixed a = Now, if we were able to make null te Rest we would obtain an L-stable FSRK metod. We ave cecked tat it is impossible to cancel te Rest because by substituting te previous values we observe tat e 44 as te constant value Neverteless, we can determine te remaining free parameters in order to minimize to some extent te contribution of te Rest to te amplification function. Due to te complexity, size and nonlinearity of te last expressions, we ave taken a 1 87 = in order to sorten tem; now, we can ceck tat e 14 = e 41, e 23 = e 32, e 24 = e 42 and e 43 = e 34, simplifying te subsequent study. By cancelling e 43 = e 34 and e 33 we can fix te parameters a 1 65 and a By using tis option, te remaining coefficients of te Rest ave te following values: e 14 = , e 23 = and e 24 = Wit te imposed restrictions, ave obtained a family of FSRK metods wic depend on te parameters c 4, c 5, c 6 c 7, attain fourt order and are A-stable. In te last part of te construction of tis metod we solve te remaining conditions. In suc order conditions, some coefficients of te last explicit level (C, A 3, b 3 ) appear. Namely we must solve: (C, A 3, b 3 ): (r 3 ), (γ 1 ),, (γ 8 ), ( 13 ), ( 23 ), ( 31 ), ( 32 ), ( 13 ), ( 23 ), ( 31 ), ( 32 ), ( 13 ), ( 23 ), ( 31 ), ( 32 ), ( 113 ), ( 131 ), ( 311 ), ( 331 ), ( 313 ), ( 133 ), ( 223 ), ( 232 ), ( 322 ), ( 332 ), ( 323 ), ( 233 ), ( 123 ), ( 132 ), ( 213 ), ( 231 ), ( 321 ) and ( 312 ). By resolving tese order conditions we obtain a family of metods wic depend on te parameters: a 3 63, a 3 73, a 3 74, a 3 76, a 3 82, a 3 83, a 3 84, a 3 85, a 3 86, a 3 87, c 4, c 5, c 6 and c 7. Suc coefficients are fixed in order to minimize to some extent te main term of te local error and also to simplify coefficients. Tus, te final metod is given by: A 1 (b 1 ) T = = B A
11 New efficient time integrators for non-linear parabolic problems 417 A 2 (b 2 ) T = B A A 3 (b 3 ) T = C A Ce=(, , , 1 3, , 2 3, 1 3, ,1)T. 4. NUMERICAL EXPERIMENT In tis section we sow a numerical test were we ave integrated te following non-linear convection-diffusion problem wit a non-linear reaction term r(y) = k 1 y + k 2 y + y3 (1+y 2 ) 2 y t = y v 1 y x1 v 2 y x2 r(y) + f, (x 1, x 2, t) Ω [, 5], y(x 1,, t) = y(x 1, 1, t) =, x 1 [, 1] and t [, 5], y(, x 2, t) = y(1, x 2, t) =, x 2 [, 1] and t [, 5], y(x 1, x 2, ) = x 3 1(1 x 1 ) 3 x 3 2(1 x 2 ) 3, (x 1, x 2 ) Ω, wit Ω = (, 1) (, 1), v 1 = 1 + x 1 x 2 e t, v 2 = 1 + x 1 t, k 1 = k 2 = 1+(x1+x2)2 e t 2 and f = 1 3 e t x 1 (1 x 1 )x 2 (1 x 2 ). We ave carried out te time discretization using te metod described in te previous section and we ave discretized in space on a rectangular mes by using a central difference sceme. Te combination of metods of finite differences type of order q in te spatial discretization stage and pt order time integrators of te type described in tis paper provides a totally discrete sceme wose global error is [y(t m )] Y m = O(q + τ p ); ere [v] denotes te restriction of te function v to te mes node.
12 418 Blanca Bujanda, Juan Carlos Jorge In order to apply te integration metod described ere, we ave taken for i = 1, 2 and j = 1,, 9 L i (t m,j )Y m,j = δ xix i Y m,j [v i ] δˆxi Y m,j g (t m,j, Y m,j ) = 3(Y m)2 (Y m)4 (1 + (Y m)2 ) 3 Y m,j (Y m,j ) 3 (1 + (Y m,j ) 2 ) 2 [k i ] Y m,j 3(Y m )2 (Y m )4 2(1 + (Y m )2 ) 3 Y m,j, were δ xix i and δˆxi denote te classical central differences. Note tat, by using tis decomposition in eac stage of (4) we must solve only one linear system wose matrix is tridiagonal; tus te computational complexity of computing te internal stages of tis metod is of te same order as in te cases of an explicit metod. Note tat tese operators L i (t) are not commutative: neverteless, te numerical results obtained are correct. Te numerical maximum global errors ave been estimated as follows E N,τ = max Y N,τ (x 1i, x 2i, t m ) Y x 1i,x 2i,t m were Y N,τ (x 1i, x 2i, t m ) are te numerical solutions obtained at te mes point (x 1i, x 2i ) and te time point t m = m τ, on a rectangular mes wit (N + 1) (N + 1) nodes and wit time step τ, wile Y is te numerical solution calculated at te same mes points and time steps, but by using a spatial mes wit (2N + 1) (2N + 1) nodes, alving te mes size, and wit time step τ 2. We compute te numerical orders of convergence as log 2 E N,τ E 2N,τ/2. Table 1. Numerical Errors (E N,τ ) N = 8 N = 16 N = 32 N = 64 N = 128 N = 256 N = E E E E E E E-6 In Table 1, we sow te numerical errors and, in Table 2, teir corresponding numerical orders of convergence. In order to guarantee tat te contribution to te error of te spatial and temporal part are of te same size, we ave assumed te relation Nτ C =.1 between te time step τ and te mes size 1 N. Table 2. Numerical orders of convergence N = 8 N = 16 N = 32 N = 64 N = 128 N =
13 New efficient time integrators for non-linear parabolic problems 419 REFERENCES [1] B. Bujanda, J.C. Jorge, Stability results for linearly implicit Fractional Step discretizations of non-linear time dependent parabolic problems, to appears in Appl. Num. Mat. [2] B. Bujanda, J.C. Jorge, Order conditions for linearly implicit Fractional Step Runge Kutta metods (Preprint), Universidad Pública de Navarra. [3] B. Bujanda, J.C. Jorge, Efficient linearly implicit metods for nonlinear multidimensional parabolic problems, J. Comput. Appl. Mat. 164/165 (24), [4] B. Bujanda, J.C. Jorge, Stability results for fractional step discretizations of time dependent coefficient evolutionary problems, Appl. Numer. Mat. 38, (21), [5] G.J. Cooper, A. Sayfy, Aditive metods for te numerical solution of ordinary differential equations, Mat. of Comp. 35, (198), [6] E. Hairer, G. Wanner, Solving ordinary differential equations II, Springer-Verlag, [7] W. Hundsdorfer, J.G. Verwer, Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations, Springer, 23. [8] N.N. Yanenko, Te metod of fractional steps, Springer, Blanca Bujanda Universidad Pública de Navarra Dpto. Matemática e Informática Juan Carlos Jorge jcjorge@unavarra.es Universidad Pública de Navarra Dpto. Matemática e Informática Received: October 3, 25.
Some Applications of Fractional Step Runge-Kutta Methods
Some Applications of Fractional Step Runge-Kutta Metods JORGE, J.C., PORTERO, L. Dpto. de Matematica e Informatica Universidad Publica de Navarra Campus Arrosadia, s/n 3006 Pamplona Navarra SPAIN Abstract:
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More informationAn Approximation to the Solution of the Brusselator System by Adomian Decomposition Method and Comparing the Results with Runge-Kutta Method
Int. J. Contemp. Mat. Sciences, Vol. 2, 27, no. 2, 983-989 An Approximation to te Solution of te Brusselator System by Adomian Decomposition Metod and Comparing te Results wit Runge-Kutta Metod J. Biazar
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationNumerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions
Proc. of te 8t WSEAS Int. Conf. on Matematical Metods and Computational Tecniques in Electrical Engineering, Bucarest, October 16-17, 2006 219 Numerical Solution of One Dimensional Nonlinear Longitudinal
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 informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
More informationResearch Article Cubic Spline Iterative Method for Poisson s Equation in Cylindrical Polar Coordinates
International Scolarly Researc Network ISRN Matematical Pysics Volume 202, Article ID 2456, pages doi:0.5402/202/2456 Researc Article Cubic Spline Iterative Metod for Poisson s Equation in Cylindrical
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationN igerian Journal of M athematics and Applications V olume 23, (2014), 1 13
N igerian Journal of M atematics and Applications V olume 23, (24), 3 c N ig. J. M at. Appl. ttp : //www.kwsman.com CONSTRUCTION OF POLYNOMIAL BASIS AND ITS APPLICATION TO ORDINARY DIFFERENTIAL EQUATIONS
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationRecent Progress in the Integration of Poisson Systems via the Mid Point Rule and Runge Kutta Algorithm
Recent Progress in te Integration of Poisson Systems via te Mid Point Rule and Runge Kutta Algoritm Klaus Bucner, Mircea Craioveanu and Mircea Puta Abstract Some recent progress in te integration of Poisson
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationThird order Approximation on Icosahedral Great Circle Grids on the Sphere. J. Steppeler, P. Ripodas DWD Langen 2006
Tird order Approximation on Icosaedral Great Circle Grids on te Spere J. Steppeler, P. Ripodas DWD Langen 2006 Deasirable features of discretisation metods on te spere Great circle divisions of te spere:
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationLecture 2: Symplectic integrators
Geometric Numerical Integration TU Müncen Ernst Hairer January February 010 Lecture : Symplectic integrators Table of contents 1 Basic symplectic integration scemes 1 Symplectic Runge Kutta metods 4 3
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationSome Review Problems for First Midterm Mathematics 1300, Calculus 1
Some Review Problems for First Midterm Matematics 00, Calculus. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd,
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationFlapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using the differential transform method
Meccanica 2006) 41:661 670 DOI 10.1007/s11012-006-9012-z Flapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using te differential transform metod Ozge Ozdemir Ozgumus
More informationOn convexity of polynomial paths and generalized majorizations
On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationParabolic PDEs: time approximation Implicit Euler
Part IX, Capter 53 Parabolic PDEs: time approximation We are concerned in tis capter wit bot te time and te space approximation of te model problem (52.4). We adopt te metod of line introduced in 52.2.
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationWeak second order explicit exponential Runge-Kutta methods for stochastic differential equations
CSSE-44 September 8, 015 Weak second order explicit exponential Runge-Kutta metods for stocastic differential equations Yosio Komori, David Coen and Kevin Burrage, Department of Systems Design and Informatics
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationInfluence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations
International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 2, pp. 281 302 (2017) ttp://campus.mst.edu/ijde Influence of te Stepsize on Hyers Ulam Stability of First-Order Homogeneous
More informationA CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Number 3, September 2003 A CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION GH. MICULA, E. SANTI, AND M. G. CIMORONI Dedicated to
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationMA455 Manifolds Solutions 1 May 2008
MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More 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 information1 + t5 dt with respect to x. du = 2. dg du = f(u). du dx. dg dx = dg. du du. dg du. dx = 4x3. - page 1 -
Eercise. Find te derivative of g( 3 + t5 dt wit respect to. Solution: Te integrand is f(t + t 5. By FTC, f( + 5. Eercise. Find te derivative of e t2 dt wit respect to. Solution: Te integrand is f(t e t2.
More informationOnline Appendix for Lerner Symmetry: A Modern Treatment
Online Appendix or Lerner Symmetry: A Modern Treatment Arnaud Costinot MIT Iván Werning MIT May 2018 Abstract Tis Appendix provides te proos o Teorem 1, Teorem 2, and Proposition 1. 1 Perect Competition
More informationMATH1131/1141 Calculus Test S1 v8a
MATH/ Calculus Test 8 S v8a October, 7 Tese solutions were written by Joann Blanco, typed by Brendan Trin and edited by Mattew Yan and Henderson Ko Please be etical wit tis resource It is for te use of
More informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More informationApplication of numerical integration methods to continuously variable transmission dynamics models a
ttps://doi.org/10.1051/ssconf/2018440005 Application of numerical integration metods to continuously variable transmission dynamics models a Orlov Stepan 1 1 Peter te Great St. Petersburg Polytecnic University
More information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationGeneric maximum nullity of a graph
Generic maximum nullity of a grap Leslie Hogben Bryan Sader Marc 5, 2008 Abstract For a grap G of order n, te maximum nullity of G is defined to be te largest possible nullity over all real symmetric n
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationOSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix
Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationQuasiperiodic phenomena in the Van der Pol - Mathieu equation
Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van
More informationMath 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0
3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,
More informationApproximation of the Viability Kernel
Approximation of te Viability Kernel Patrick Saint-Pierre CEREMADE, Université Paris-Daupine Place du Marécal de Lattre de Tassigny 75775 Paris cedex 16 26 october 1990 Abstract We study recursive inclusions
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationDigital Filter Structures
Digital Filter Structures Te convolution sum description of an LTI discrete-time system can, in principle, be used to implement te system For an IIR finite-dimensional system tis approac is not practical
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationlecture 35: Linear Multistep Mehods: Truncation Error
88 lecture 5: Linear Multistep Meods: Truncation Error 5.5 Linear ultistep etods One-step etods construct an approxiate solution x k+ x(t k+ ) using only one previous approxiation, x k. Tis approac enoys
More informationTwo Step Hybrid Block Method with Two Generalized Off-step Points for Solving Second Ordinary Order Differential Equations Directly
Global Journal of Pure and Applied Matematics. ISSN 0973-768 Volume 2, Number 2 (206), pp. 59-535 Researc India Publications ttp://www.ripublication.com/gjpam.tm Two Step Hybrid Block Metod wit Two Generalized
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationNotes on Multigrid Methods
Notes on Multigrid Metods Qingai Zang April, 17 Motivation of multigrids. Te convergence rates of classical iterative metod depend on te grid spacing, or problem size. In contrast, convergence rates of
More informationMath 312 Lecture Notes Modeling
Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a
More information2.3 More Differentiation Patterns
144 te derivative 2.3 More Differentiation Patterns Polynomials are very useful, but tey are not te only functions we need. Tis section uses te ideas of te two previous sections to develop tecniques for
More informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
More informationEnergy-preserving variant of collocation methods 1
European Society of Computational Metods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Industrial and Applied Matematics (JNAIAM) vol.?, no.?, 21, pp. 1-12 ISSN 179 814 Energy-preserving
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips
More informationEfficient algorithms for for clone items detection
Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire
More information