PARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM
|
|
- Annabelle Thompson
- 5 years ago
- Views:
Transcription
1 PARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM Daoud S Daoud Dept. of Mathematics Eastern Mediterranean University Famagusta-North Cyprus via Mersin 10-Turkey. e mail: daoud.daoud@emu.edu.tr 1/37 M5: Space-Time Parallel Methods for PDEs
2 OUTLINE 1. The Control Model Problem 2. The Parareal Algorithm 3. The Parreal-Inverse Algorithm 4. The Stability Analysis of Parreal- Control Algorithm 5. Discussion and Conclusions 2/37 M5: Space-Time Parallel Methods for PDEs
3 THE MODEL PROBLEM u t = u xx + p(t)u + φ(x, t), 0 x 1, 0 tu(x, 0) = f(x), 0 x 1 u(0, t) = g 0 (t), 0 < x T u(1, t) = g 1 (t), 0 < x T subject to the over specification condition at a point x 0 in the spatial domain (1) u(x 0, t) = E(t), (2) where f, g 0, g 1, E and φ are known functions while the functions u and p are unknown, and 1 < p(t) < 0 for t [0, T ]. 3/37 M5: Space-Time Parallel Methods for PDEs
4 CONTROL PROBLEM RELATED REFERENCES 1. J.R. Cannon, Yanping Lin, S.Wang, Determination of source parameter in parabolic equations. Meccanica, vol. 27:85-94, J.R. Cannon, Yanping Lin, S.Wang, Numerical procedures for the determination of an unknown coefficient in semi -linear parabolic equations. Inverse problems, vol. 10: , Yanping Lin, An inverse problem for a class of quasilinear parabolic equation, SIAM J. Math. Anal. Vol. 22: , /37 M5: Space-Time Parallel Methods for PDEs
5 The Engineering Applications of The Control Problem 1- Heat Conduction process 2- Thermoealsiticity 3-Chemical Diffusion and Control Theory 5/37 M5: Space-Time Parallel Methods for PDEs
6 THE CONTROL PROBLEM SOLUTION ALGORITHM Algorithm 1: The Solution of the Control Problem Algorithm Step 1: Consider the backward Euler scheme for the solution of (1) given by u n i t where p (t n+1 ) is the predicted value of p(t) at time T n+1 u n+1 i = h 2 (u n+1 i 1 2un+1 i + u n+1 i+1 ) + p (t n+1 )u n+1 i + φ(x 0, t n+1 ) (3) Step 2: Update the value of the function p (t n+1 ) using (2) u(x 0, t) = E(t), such that p(t n+1 ) l+1 = E (t n+1 ) h 2 (u n+1(l) x 0 h for l = 0, 1,... 2un+1(l) x 0 E(t n+1 ) + u n+1(l) x 0 +h ) φ(x 0, t n+1 ) The value of the pair (u n+1, p n+1(l+1) ) will be corrected using Step 1 and Step 2, recursively. Step 3: if u n+1(l+1) u n+1(l) < ε and p n+1(l+1) p n+1(l) < ε, then stop otherwise proceed with step 1 and step 2, ε is a pre specified accuracy. (4) 6/37 M5: Space-Time Parallel Methods for PDEs
7 RELATED REFERENCES TO THE PARAEAL ALGORITHM 1- J.-L. Lions, Y. Maday, and G. Turinici, Resolution d EDP par un scheme en temps parareal, C.R. Acad. Sci. Paris Ser. I math., 332: , G. Bal and Y. Maday, A parareal time discretization for non linear PDEs with application to the pricing of an american put, Proceeding of the Workshop on Domain Decomposition, Zurich, Switzerland, LNCSE series, Springer Verlag, G. Bal, Parallelization in time of(stochastic) ordinary differential equations, Math. Meth. Anal.Num (to appaer) C. Farhat and M. Chandesris, time decomposed parallel time-integrators;theory and feasibility studies for fluid, structure,and fluid-structire applications, Int.J.Numer.Meth. Engng. 58: , Y. Maday and G. Turininc, Parallel in time algorithms for quantum control: Parareal time discretization scheme, Int. J. Quant. Chem. 93: , G.A. Staff, Convergence and stability of the parareal algorithm: A numerical and theoretical investigation, preprint Numerics no. 2/2003, Norwegian University of Science and Technology, Norway. 7- G.A. Staff and E.M. Ronquist, Stability of the parareal algorithm, 7/37 M5: Space-Time Parallel Methods for PDEs
8 PARAREAL ALGORITHM Algorithm 2: Parareal Algorithm Step 1: For k = 0, and over Ω [t n, t n + n t] solve for u n+1 0 using coarse propagator u n+1 0 = G t (u n 0) n = 0, 1,..., N 1, whereu 0 0 = f(x). Step 2: Over each time interval [t n, t n + sδt], solve for v n+1 0, using fine propagator v n+1 0 = F sδt (v n 0 ) = i=s i=1 F iδt (v n 0 ) where v n 0 = u n 0 Step 3: Evaluate the jumps S n 0 = F sδt (u n 0) G t (u n 0), then u n+1 1 = G t (u n 0) + S n 0 8/37 M5: Space-Time Parallel Methods for PDEs
9 Step 4: Correction step; For k 1, and n = 1,..., N 1 u n+1 k = G t (u n k) + F sδt (u n k 1) G t (u n k 1) Step 5: If all jumps are not sufficiently small, the process can be iterated as follows u n+1 k = G t (u n k) + F sδt (u n k 1) G t (u n k 1) u 0 k = u 0 for all k 9/37 M5: Space-Time Parallel Methods for PDEs
10 The Solution of The Control Problem by Parareal Algorithm Algorithm 3: Parareal- Control Problem Algorithm Step 1: Over the domain Ω [t n, t n+1 ] and for k = 1, consider the coarse propagator i.e. u n+1 1 u n 1 t = (u xx ) n p(t n )u n+1 n = 0, 1,... N 1 The solution u n+1 1 denoted by G t (u n 1) p(t n+1 ) correction: Consider the correction of p(t) by the following relation p(t n+1 ) = E (t n+1 ) (u xx ) 1 (x0,t n+1 ) φ(x 0, t n+1 ) E(t n+1 ) Step 2: For k + 1 > 1 and over the domain Ω [t n, t n+1 ] consider the coarse propagator i.e. u n+1 k+1 un k+1 t = (u xx ) n+1 k+1 + p(t n)u n+1 k+1 n = 0, 1,... N 1 10/37 M5: Space-Time Parallel Methods for PDEs
11 p(t n+1 ) correction: Consider the correction of p(t) by the following relation p(t n+1 ) = E (t n+1 ) (u xx ) k+1 (x0,t n+1 ) f(x 0, t n+1 ) E(t n+1 ) The solution u n+1 k+1 is denoted by G t(u n k+1 ) Step 3: consider the fine propagator solution over Ω [t n, t n+l ], l = 1,..., s 1 Solve for u n+l k The solution u n+s k = u n+1 then for for l = 1,... s 1 where s = t δt The solution u n+1 k+1 u n+l 1 k δt = (u xx ) n+l 1 k k is denoted by F sδt (u n k ) + p(t n+l 1 )u n+l 1 k p(t n+l ) = E (t n+l ) u xx,k (x0,t n+l ) φ(x 0, t n+l ), E(t n+l ) is given by u n+1 k+1 = G t(u n k+1) + F sδt (u n k) G t (u n k) (5) 11/37 M5: Space-Time Parallel Methods for PDEs
12 Stability of the Parareal-Control Algorithm Let u(x, t) be the solution of such that u t = u xx + p(t)u(t), (6) u(0, t) = 0, u(1, t) = 0 and u(x, 0) = u 0. (7) The spatial derivative operator is approximated by the second order central difference approximation given by (u xx ) (xi,t) h 2 [u(x i+1, t) 2u(x i, t) + u(x i 1, t)] + O(h 2 ). (8) 12/37 M5: Space-Time Parallel Methods for PDEs
13 Consider the Fourier transform of the model problem, and over the Fourier domain the problem corresponding to (6) is given by û t = Q(ξ, t)û(ξ, t) (9) where Q(ξ, t) = q(ξ) + p(t), (p(t) = p(t)) such that Q(u) = q(ξ) = 2h 2 sin 2 (ξ/2), and h is the space mesh spacing. Then Q(ξ)û(ξ), Q(ξ, t) = q(ξ) + p(t) = 2h 2 sin 2 ( ξ ) + p(t). (10) 2 13/37 M5: Space-Time Parallel Methods for PDEs
14 For the numerical solution of the model problem we defined the coarse time interval t such that t n+1 = t n + t, and the fine solution time step size is give by δt, such that t δt = s, t = T N Over the time interval [t n, t n + t], and corresponding to (9), and (10) the amplification factor (Stability Function) of the considered schemes over the Fourier domain are given by, û n+1 (ξ, t n+1 ) = (1 Q(ξ, t n ) t) 1 û n (ξ, t n ) ( 1 = 1 + (2h 2 sin 2 ( ξ 2 n)) t) ) p(t ûn (ξ, t n ) (11) which is unconditional stable for p(t) < 0,. 14/37 M5: Space-Time Parallel Methods for PDEs
15 The amplification factor of the explicit scheme, over the time interval [t n, t n + sδt] is given by û n+1 (ξ, t n ) = s i=1 (1 + Q(ξ, t n+i 1)δt)û n (ξ, t n ) û n (ξ, t n ) = (12) s i=1 (1 + ( 2h 2 sin 2 ( ξ 2 ) + p(t n+i 1))δt)û n (ξ, t n ) and its conditional stable scheme according to the explicit scheme condition and for any p(t). 15/37 M5: Space-Time Parallel Methods for PDEs
16 û n k 1 û k n 1 û k 1 n 1 2 û k n 2 û k 1 n 2 û k 2 n 2 3 û k n 3 û k 1 n 3 û k 2 n 3 û k 3 n 3 4 û k n 4 û k 1 n 4 û k 2 n 4 û k 3 n 4 û k 4 n 4 The Parareal Iterative Scheme 16/37 M5: Space-Time Parallel Methods for PDEs
17 û n+1 k+1 = G t(û n k+1) + F sδt (û n k) G t (û n k) (13) Subtitling (11) and (12) in (13) then the solution û n+1 k+1 parareal algorithm, algorithm, is given by of the model problem by the û n+1 k+1 = (1 Q(ξ, t n) t) 1 û n k+1 + s (1 + Q(ξ, t n+i 1 )δt)û n k (1 Q(ξ, t n ) t) 1 û n k i=1 (14) At this point we will present the stability studies for the following cases case 1- δt = t (s = 1) for sequential algorithm case 2- δt = t s (s > 1) for Parareal algorithm 17/37 M5: Space-Time Parallel Methods for PDEs
18 Case I, δt = t For the first case and over the time interval [t n, t n + t] with time step δt = t, we obtain the following relation û n+1 k+1 = (1 Q(ξ, t n) t) 1 û n k+1 + (1 + Q(ξ, t n) t)û n k (1 Q(ξ, t n) t) 1 û n k Let (15) F (ξ, t n ) = (1 + Q(ξ, t n ) t) 1 (1 Q(ξ, t n ) t) 1 (16) Then (15) is re written as follows û n+1 k+1 = (1 + Q(ξ, t n) t)f (ξ, t n ) }{{} ûn k+1 + (1 + Q(ξ, t n ) t)(1 F (ξ, t n )) }{{} ûn k (17) where û 0 k+1 = û0, 18/37 M5: Space-Time Parallel Methods for PDEs
19 For the stability of the scheme given by (17) we will consider the approach of analysis by G. Staff Parareal algorithm-a survey of present work, 2003 Then the iterated solution û n+1 k+1 ) û n+1 k+1 = k i=0 ( n i is given by ((1 + Q(ξ, t n ) t)(1 F (ξ, t n ))) i ((1 + Q(ξ, t n ) t)f (ξ, t n )) n i û 0 and then the amplification factor (Stability Function) of the parreal algorithm (15), is given by T (ξ, t n ) = k i=0 ( n i ) ((1 + Q(ξ, t n ) t)(1 F (ξ, t n ))) i (1 + Q(ξ, t n ) t)f (ξ, t n )) n i (18) 19/37 M5: Space-Time Parallel Methods for PDEs
20 The algorithm is stable whenever Consider, sup 1 n N sup 1 k N T (ξ, t n ) 1 (19) T (ξ, t n ) < k i=0 ( n i ) ((1 + Q(ξ, t n ) t)(1 F (ξ, t n ))) i ((1 + Q(ξ, t n ) t)f (ξ, t n )) n i < n ( (1 + Q(ξ, t n ) t)(1 F (ξ, t n )) + (1 + Q(ξ, t n ) t)f (ξ, t n ) ) n i=0 T (ξ, t n ) < ( ((1 + Q(ξ, t n ) t)(1 F (ξ, t n ))) + ((1 + Q(ξ, t n ) t)f (ξ, t n )) ) n 20/37 M5: Space-Time Parallel Methods for PDEs
21 and hence for stability conditions is to show that the stability function, ρ(ξ, t n ), satisfying ρ(ξ, t n ) = (1 + Q(ξ, t n ) t)(1 F (ξ, t n )) (1 + Q(ξ, t n ) t)f (ξ, t n ) < 1 (20) and to investigate the necessary conditions for r and p(t) as well. Consider ρ(ξ, t n ) = (1 Q(ξ, t n ) t) 1 ( 1 + (Q(ξ, t n ) t) 2) < 1 Since Q(ξ, t n ) = 2h 2 sin 2 (ξ) + p(t n ), then 1+ ( 2h 2 sin 2 (ξ/2) + p(t n )) t ) 2 = ( 1 + 4r 2 sin 4 (ξ/2) 4r sin 2 (ξ/2)p(t n ) t + p 2 (t n ) t Therefore (1 + 4r 2 sin 4 (ξ/2) 4r sin 2 (ξ/2)p(t n ) t + p 2 (t n ) t 2) < 1 + 2r sin 2 (ξ/2) ( 2r sin 2 (ξ/2) 2p(t n ) t ) p(t n ) t for 1 < p(t n ) < 0 and for all n = 0, 1,..., N, r = t/h 2 21/37 M5: Space-Time Parallel Methods for PDEs
22 Then for r < 1/4 concluding that 1 + 2r sin 2 (ξ/2) ( 2r sin 2 (ξ/2) 2p(t n ) t ) p(t n ) t < 1 + 2r sin 2 (ξ/2) (2r 2p(t n ) t) p(t n ) t < 1 + 2r sin 2 (ξ/2) (1/2 2p(t n ) t) p(t n ) t < 1 + 2r sin 2 (ξ/2) (1/2 + 2 t) p(t n ) t < since r = t < 1/4, then h r sin 2 (ξ/2) (1/2 + 1/2) p(t n ) t = (1 Q(ξ, t n ) t) Therefore, for r = t h 2 < 1/4 and 1 < p(t) < 0 concluding that ρ(ξ, t n ) = (1 + Q(ξ, t n ) t)(1 F (ξ, t n )) (1 + Q(ξ, t n ) t)f (ξ, t n ) < 1 and then T (ξ, t n ) < 1 for all n 22/37 M5: Space-Time Parallel Methods for PDEs
23 The above analysis and discussions present the proof of the main result given by the following theorem. Theorem 1: For the given control model problem (1) solved by the parareal algorithm (Algorithm 3) given by, u n+1 k+1 = G t(û n k+1) + F δt (û n k) G t (û n k) (21) where G t and F δt are the coarse and fine solvers, respectively, and for s = 1 i.e. t = δt. Then ρ(ξ, t n ) < 1 for all δt = t < 1 h 2 h 2 4 and 1 < p(t) < 0, where ρ(ξ, t n ) is the stability function induced by (15), given by (20). 23/37 M5: Space-Time Parallel Methods for PDEs
24 Case 2: s = t/δt > 1 û n+1 k+1 = (1 Q(ξ, t n) t) 1 û n k+1 + s (1 + Q(ξ, t n+i 1 )δt)û n k (1 Q(ξ, t n ) t) 1 û n k In this case the stability function of the Parareal Scheme is given by i=1 ρ(ξ, t n ) = 2(1 Q(ξ, t n ) t) 1 (22) s (1 + Q(ξ, t n )δt (23) To investigate the necessary condition to satisfy ρ(ξ, t n ) < 1, Consider i=1 24/37 M5: Space-Time Parallel Methods for PDEs
25 For the first term (1 Q(ξ, t n ) t) 1 }{{} [ 2 (1 Q(ξ, t n ) t) i=1(1 + Q(ξ, t n+i 1 )δt) }{{} (1 Q(ξ, t n ) t) 1 < 1 for any r c = t/h 2 and for 1 < p(t) < 0, which in fact represent the stability of the coarse amplification factor ( Stability function). For the second term [ ] s 2 (1 Q(ξ, t n ) t) (1 + Q(ξ, t n+i 1 )δt) (25) i=1 ] (24) 25/37 M5: Space-Time Parallel Methods for PDEs
26 perform the multiplication in (25) concluding that s (1 + Q(ξ, t n+i 1 )δt) 1 + δt i=1 s Q(ξ, t n+i 1 ) + O(δt 2 ) (26) where Q(ξ, t n+i 1 ) = 2h 2 sin 2 ( ξ 2 ) + p(t n+i 1), and substitute (26) in (25) leading to { 2 (1 Q(ξ, tn ) t) [ 1 + δt s i=1 Q(ξ, t n+i 1) + O(δt 2 ) ]} = Q(ξ, t n ) t δt(1 Q(ξ, t n ) t) s i=1 Q(ξ, t n+i 1) 1 + Q(ξ, t n ) t s i=1 Q(ξ, t n+i 1)δt + O(δt 2 ) where Q(ξ, t n+i 1 ) = 2h 2 sin 2 (ξ/2) + p(t n+i 1 ), therefore i=1 26/37 M5: Space-Time Parallel Methods for PDEs
27 1 + Q(ξ, t n ) t δt s ( 2h 2 sin 2 (ξ/2) + p(t n+i 1 )) = i=1 1 + ( 2h 2 sin 2 (ξ/2) + p(t n ) ) t δt where r c = t Therefore 1 2r c sin 2 (ξ/2) + tp(t n ) + h 2, r f = δt h 2, and t δt = s. 1 2r c sin 2 (ξ/2) + tp(t n ) + }{{} <0 s ( 2h 2 sin 2 (ξ/2) + p(t n+i 1 ) ) i=1 s ( 2rf sin 2 (ξ/2) δtp(t n+i 1 ) ) i=1 s ( 2rf sin 2 (ξ/2) + δt ) < 1 + t i=1 27/37 M5: Space-Time Parallel Methods for PDEs
28 then [ ρ(ξ, t n ) = (1 Q(ξ, t n ) t) 1 2 (1 Q(ξ, t }{{} n ) t) ] + Q(ξ, t n+i 1 )δt) }{{ i=1(1 } for 1 < p(t n ) < 0 for all t n [0, T ]. < 1 28/37 M5: Space-Time Parallel Methods for PDEs
29 For the case when s = t/δt > 1 we will have the following theorem Theorem 2: For the given control model problem (1) solved by the parareal algorithm (Algorithm 3) given by, u n+1 k+1 = G t(û n k+1) + F sδt (û n k) G t (û n k) (27) where G t and F sδt are the coarse and fine solvers, respectively, and for s = t/δt > 1. Then ρ(ξ, t n ) < 1 for all r c = t/h 2 and r f = δt/h 2 satisfying the coarse and fine propagator stability constrained, respectively, and 1 < p(t) < 0, where ρ(ξ, t n ) is the stability function induced by (15), given by (23). 29/37 M5: Space-Time Parallel Methods for PDEs
30 Model Problem u t = u xx + p(t)u + φ(x, t) defined over the domain x Ω = (0, 1) and t [0, 1]. With exact solution given by u(x, t) = e t2 (cos(πx) + sin(πx)), p(t) is given by p(t) = 1 + t 2 and φ(x, t) = e t2 (π t 2 2t)(cos(πx) + sin(πx)). Subject to u(0, t) = e t2 u(1, t) = e t2 u(x, 0) = cos(πx) + sin(πx) 30/37 M5: Space-Time Parallel Methods for PDEs
31 For comparison purpose the above model problem has been solved using 1- Control Problem Algorithm (Algorithm 1) 2- Parareal-Control Problem Algorithm (Algorithm 2) 31/37 M5: Space-Time Parallel Methods for PDEs
32 Comments and Discussions 1. Due to the iterative process by the parareal algorithm the solution of the control problem does not require updating of the value of p(t), over each single coarse time step, as required by algorithm Due to inducing of the fine time step, in the parareal algorithm, the numerical results shows how accurate result we will achieve for u(x, t) and p(t) by algorithm 2 over the use of algorithm 1 (for u(x, t) see tables 1-4, and for p(t) see tables 5-8). 3. There is flexibility in this algorithm that we may use different spatial mesh spacing in the coarse solvers, due to its unconditional stable criteria. 4. From the sketch of the stability region by the polar graphs, as shown in figs 1-2, It shows how the stability region, amplification factor, effected for different values of (δt/h 2 ), and p(t) values as well. 5. The considered specific Case 1 is designated for the sequential type of algorithm to solve the Control problem, and by using Case 2 it shows that by Parareal algorithm the constrained on the ratio δt/h 2 is not required. 32/37 M5: Space-Time Parallel Methods for PDEs
33 Coarse r c = 0.2, Fine r f = 0.2 Coarse r c = 0.2, Fine r f = r f =0.2 r c = r f =0.2 r c =0.2 p(t)= /37 M5: Space-Time Parallel Methods for PDEs
34 Coarse r c = 0.25, Fine r f = 0.25 Coarse r c = 0.4, Fine r f = r f =0.25 r c = r f =0.4 r c = Coarse r c = 0.5, Fine r f = r f =0.5 r c = /37 M5: Space-Time Parallel Methods for PDEs
35 r c = 1.25 r c = 1 r c = 0.8 r c = 0.5 r c = 0.25 r c = 0.2 u(x, t) e e e e e e-4 p(t) e e-2 Table 1: The value of u(x, t) and p(t) using 5 corrections only for algorithm 1 k r c = 1.25 r c = 1 r c = 0.5 r c = 0.25 r f = 0.25 r f = 0.25 r f = 0.25 r f = e e e e e e e e e e e e e e e e e e e e-4 Table 2: The error for different s ratios, for r f = 0.25 k r c = 1.25 r c = 1 r c = 0.5 r c = r f = r f = r f = r f = e e e e e e e e e e e e e e e e e e e e-5 Table 3: The error for different s ratios, for r f = k r c = 1 r c = 0.8 r c = 0.4 r c = 0.2 r f = 0.2 r f = 0.2 r f = 0.2 r f = e e e e e e e e e e e e e e e e e e e e-5 Table 4: The error for different s ratios, for r f = /37 M5: Space-Time Parallel Methods for PDEs
36 r c = 1.25 r c = 1 r c = 0.8 r c = 0.5 r c = 0.25 r c = 0.2 u(x, t) e e e e e e-4 p(t) e e-2 Table 5: The value of u(x, t) and p(t) using 5 corrections only for algorithm 1 k r c = 1.25 r c = 1 r c = 0.5 r c = 0.25 r f = 0.25 r f = 0.25 r f = 0.25 r f = e e e e e e e-2 Table 6: The error of p(t) for different s ratios, for r f = 0.25 k r c = 1 r c = 0.8 r c = 0.4 r c = 0.2 r f = 0.2 r f = 0.2 r f = 0.2 r f = e e e e e e e-2 Table 7: The error of p(t) for different s ratios, for r f = 0.2 k r c = 1.25 r c = 1 r c = 0.5 r c = r f = r f = r f = r f = e e e e e e e e e-2 Table 8: The error p(t)for different s ratios, for r f = /37 M5: Space-Time Parallel Methods for PDEs
37 k r c = 1 r c = 1 r c = 1 r c = 0.5 r c = 0.5 r f = r f = 0.2 r f = 0.25 r f = 0.25 r f = e e e e e e e e e e e e e e e e e e e e e e e e e-4 Table 9: The solution for different ratios of Coarse and Fine solvers 37/37 M5: Space-Time Parallel Methods for PDEs
Stability of the Parareal Algorithm
Stability of the Parareal Algorithm Gunnar Andreas Staff and Einar M. Rønquist Norwegian University of Science and Technology Department of Mathematical Sciences Summary. We discuss the stability of the
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More informationRATIO BASED PARALLEL TIME INTEGRATION (RAPTI) FOR INITIAL VALUE PROBLEMS
RATIO BASED PARALLEL TIME INTEGRATION (RAPTI) FOR INITIAL VALUE PROBLEMS Nabil R. Nassif Noha Makhoul Karam Yeran Soukiassian Mathematics Department, American University of Beirut IRISA, Campus Beaulieu,
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More informationParareal Methods. Scott Field. December 2, 2009
December 2, 2009 Outline Matlab example : Advantages, disadvantages, and survey of usage in literature Overview Overview Notation Algorithm Model equation example Proposed by Lions, Maday, Turinici in
More informationA parallel algorithm for the heat equation with derivative boundary conditions
International Mathematical Forum, 2, 2007, no. 12, 565-574 A parallel algorithm for the heat equation with derivative boundary conditions M. Akram University College of Information Technology University
More informationAn Adjoint Approach for Stabilizing the Parareal Method
An Adjoint Approach for Stabilizing the Parareal Method Feng Chen a, Jan S Hesthaven b, Yvon Maday c, Allan Nielsen b a Department of Mathematics, Baruch College, New York, NY 10010, USA b Mathematics
More informationnorges teknisk-naturvitenskapelige universitet Parallelization in time for thermo-viscoplastic problems in extrusion of aluminium
norges teknisk-naturvitenskapelige universitet Parallelization in time for thermo-viscoplastic problems in extrusion of aluminium by E. Celledoni, T. Kvamsdal preprint numerics no. 7/2007 norwegian university
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 29, No. 2, pp. 556 578 c 27 Society for Industrial and Applied Mathematics ANALYSIS OF THE PARAREAL TIME-PARALLEL TIME-INTEGRATION METHOD MARTIN J. GANDER AND STEFAN VANDEWALLE
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationPARALLELIZATION IN TIME OF (STOCHASTIC) ORDINARY DIFFERENTIAL EQUATIONS
PARALLELIZATION IN TIME OF (STOCHASTIC) ORDINARY DIFFERENTIAL EQUATIONS GUILLAUME BAL Abstract. This paper analyzes some properties of the parareal algorithm, which can be used to parallelize the time
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationParallel in Time Algorithms for Multiscale Dynamical Systems using Interpolation and Neural Networks
Parallel in Time Algorithms for Multiscale Dynamical Systems using Interpolation and Neural Networks Gopal Yalla Björn Engquist University of Texas at Austin Institute for Computational Engineering and
More informationMATH 251 Final Examination August 10, 2011 FORM A. Name: Student Number: Section:
MATH 251 Final Examination August 10, 2011 FORM A Name: Student Number: Section: This exam has 10 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all work
More informationHigher-Order Difference and Higher-Order Splitting Methods for 2D Parabolic Problems with Mixed Derivatives
International Mathematical Forum, 2, 2007, no. 67, 3339-3350 Higher-Order Difference and Higher-Order Splitting Methods for 2D Parabolic Problems with Mixed Derivatives Jürgen Geiser Department of Mathematics
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
More informationPartial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators
PDE-SEP-HEAT-1 Partial Differential Equations Separation of Variables 1 Partial Differential Equations and Operators et C = C(R 2 ) be the collection of infinitely differentiable functions from the plane
More informationMATH 819 FALL We considered solutions of this equation on the domain Ū, where
MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,
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 informationAnalytical solution for determination the control parameter in the inverse parabolic equation using HAM
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017, pp. 1072 1087 Applications and Applied Mathematics: An International Journal (AAM Analytical solution
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 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 informationSpacetime Computing. A Dynamical Systems Approach to Turbulence. Bruce M. Boghosian. Department of Mathematics, Tufts University
Spacetime Computing A Dynamical Systems Approach to Turbulence Bruce M. Boghosian, Tufts University Presented at 3rd NA-HPC Roadmap Workshop, Royal Society, London January 26, 2009 Acknowledgements: Hui
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationA numerical method for the wave equation subject to a non-local conservation condition
A numerical method for the wave equation subject to a non-local conservation condition Whye-Teong Ang Division of Engineering Mechanics School of Mechanical and Production Engineering Nanyang Technological
More informationMath 251 December 14, 2005 Answer Key to Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Name Section Math 51 December 14, 5 Answer Key to Final Exam There are 1 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning
More informationNumerical methods for a fractional diffusion/anti-diffusion equation
Numerical methods for a fractional diffusion/anti-diffusion equation Afaf Bouharguane Institut de Mathématiques de Bordeaux (IMB), Université Bordeaux 1, France Berlin, November 2012 Afaf Bouharguane Numerical
More informationConvergence of a first order scheme for a non local eikonal equation
Convergence of a first order scheme for a non local eikonal equation O. Alvarez, E. Carlini, R. Monneau, E. Rouy March 22, 2005 Abstract We prove the convergence of a first order finite difference scheme
More informationOverlapping Schwarz Waveform Relaxation for Parabolic Problems
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03038-7 Overlapping Schwarz Waveform Relaxation for Parabolic Problems Martin J. Gander 1. Introduction We analyze a new domain decomposition algorithm
More informationON NUMERICAL SOLUTIONS OF PARTIAL DIFFERENTIAL EQUATIONS BY THE DECOMPOSITION METHOD. Mustafa Inc
153 Kragujevac J. Math. 26 (2004) 153 164. ON NUMERICAL SOLUTIONS OF PARTIAL DIFFERENTIAL EQUATIONS BY THE DECOMPOSITION METHOD Mustafa Inc Department of Mathematics, Firat University, 23119 Elazig Turkiye
More informationBang bang control of elliptic and parabolic PDEs
1/26 Bang bang control of elliptic and parabolic PDEs Michael Hinze (joint work with Nikolaus von Daniels & Klaus Deckelnick) Jackson, July 24, 2018 2/26 Model problem (P) min J(y, u) = 1 u U ad,y Y ad
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationThe Heat Equation John K. Hunter February 15, The heat equation on a circle
The Heat Equation John K. Hunter February 15, 007 The heat equation on a circle We consider the diffusion of heat in an insulated circular ring. We let t [0, ) denote time and x T a spatial coordinate
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
More informationNumerical schemes of resolution of stochastic optimal control HJB equation
Numerical schemes of resolution of stochastic optimal control HJB equation Elisabeth Ottenwaelter Journée des doctorants 7 mars 2007 Équipe COMMANDS Elisabeth Ottenwaelter (INRIA et CMAP) Numerical schemes
More informationA NUMERICAL ITERATIVE SCHEME FOR COMPUTING FINITE ORDER RANK-ONE CONVEX ENVELOPES
A NUMERICAL ITERATIVE SCHEME FOR COMPUTING FINITE ORDER RANK-ONE CONVEX ENVELOPES XIN WANG, ZHIPING LI LMAM & SCHOOL OF MATHEMATICAL SCIENCES, PEKING UNIVERSITY, BEIJING 87, P.R.CHINA Abstract. It is known
More informationMath 251 December 14, 2005 Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Math 251 December 14, 2005 Final Exam Name Section There are 10 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning of each question
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationSpectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions
Applied Mathematical Sciences, Vol. 1, 2007, no. 5, 211-218 Spectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions M. Javidi a and A. Golbabai b a Department
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationA Nonlinear Inverse Parabolic Problem
Applied Mathematical Sciences, Vol. 1, 2007, no. 24, 1181-1186 A Nonlinear Inverse Parabolic Problem N. Azizi and R. Pourgholi Department of Mathematics Islamic Azad University of Karaj Branch, Iran azizin41@yahoo.com,
More informationFinite Element Decompositions for Stable Time Integration of Flow Equations
MAX PLANCK INSTITUT August 13, 2015 Finite Element Decompositions for Stable Time Integration of Flow Equations Jan Heiland, Robert Altmann (TU Berlin) ICIAM 2015 Beijing DYNAMIK KOMPLEXER TECHNISCHER
More informationGreedy control. Martin Lazar University of Dubrovnik. Opatija, th Najman Conference. Joint work with E: Zuazua, UAM, Madrid
Greedy control Martin Lazar University of Dubrovnik Opatija, 2015 4th Najman Conference Joint work with E: Zuazua, UAM, Madrid Outline Parametric dependent systems Reduced basis methods Greedy control
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 informationNational Taiwan University
National Taiwan University Meshless Methods for Scientific Computing (Advisor: C.S. Chen, D.L. oung) Final Project Department: Mechanical Engineering Student: Kai-Nung Cheng SID: D9956 Date: Jan. 8 The
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationNumerical methods in molecular dynamics and multiscale problems
Numerical methods in molecular dynamics and multiscale problems Two examples T. Lelièvre CERMICS - Ecole des Ponts ParisTech & MicMac project-team - INRIA Horizon Maths December 2012 Introduction The aim
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationIterative Domain Decomposition Methods for Singularly Perturbed Nonlinear Convection-Diffusion Equations
Iterative Domain Decomposition Methods for Singularly Perturbed Nonlinear Convection-Diffusion Equations P.A. Farrell 1, P.W. Hemker 2, G.I. Shishkin 3 and L.P. Shishkina 3 1 Department of Computer Science,
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
More informationApplication of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms
Future Generation Computer Systems () 65 7 Application of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms R. Naidoo a,b, S. Baboolal b, a Department
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 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 informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationNinth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673
Ninth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673 A NEW NUMERICAL APPROACH FOR THE SOLUTION OF CONTAMINANT TRANSPORT EQUATION Mohamed El-Gamel Department of Mathematical
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationScientific Data Computing: Lecture 3
Scientific Data Computing: Lecture 3 Benson Muite benson.muite@ut.ee 23 April 2018 Outline Monday 10-12, Liivi 2-207 Monday 12-14, Liivi 2-205 Topics Introduction, statistical methods and their applications
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationParallel-in-Time for Parabolic Optimal Control Problems Using PFASST
Zuse Institute Berlin Takustr. 7 14195 Berlin Germany SEBASTIAN GÖTSCHEL, MICHAEL L. MINION Parallel-in-Time for Parabolic Optimal Control Problems Using PFASST ZIB Report 17-51 (May 2017) Zuse Institute
More informationGreedy Control. Enrique Zuazua 1 2
Greedy Control Enrique Zuazua 1 2 DeustoTech - Bilbao, Basque Country, Spain Universidad Autónoma de Madrid, Spain Visiting Fellow of LJLL-UPMC, Paris enrique.zuazua@deusto.es http://enzuazua.net X ENAMA,
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 informationNumerical methods for Mean Field Games: additional material
Numerical methods for Mean Field Games: additional material Y. Achdou (LJLL, Université Paris-Diderot) July, 2017 Luminy Convergence results Outline 1 Convergence results 2 Variational MFGs 3 A numerical
More informationSplitting Enables Overcoming the Curse of Dimensionality
Splitting Enables Overcoming the Curse of Dimensionality Jerome Darbon Stanley Osher In this short note we briefly outline a new and remarkably fast algorithm for solving a a large class of high dimensional
More informationNumerical analysis of the one-dimensional Wave equation subject to a boundary integral specification
Numerical analysis of the one-dimensional Wave equation subject to a boundary integral specification B. Soltanalizadeh a, Reza Abazari b, S. Abbasbandy c, A. Alsaedi d a Department of Mathematics, University
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 informationASYMPTOTIC THEORY FOR WEAKLY NON-LINEAR WAVE EQUATIONS IN SEMI-INFINITE DOMAINS
Electronic Journal of Differential Equations, Vol. 004(004), No. 07, pp. 8. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ASYMPTOTIC
More informationNumerical Solutions of the Burgers System in Two Dimensions under Varied Initial and Boundary Conditions
Applied Mathematical Sciences, Vol. 6, 22, no. 3, 563-565 Numerical Solutions of the Burgers System in Two Dimensions under Varied Initial and Boundary Conditions M. C. Kweyu, W. A. Manyonge 2, A. Koross
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
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 informationOptimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients
Optimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients Victorita Dolean, Martin Gander, Erwin Veneros To cite this version: Victorita Dolean, Martin Gander, Erwin Veneros. Optimized
More informationOptimized Schwarz methods with overlap for the Helmholtz equation
Optimized Schwarz methods with overlap for the Helmholtz equation Martin J. Gander and Hui Zhang Introduction For the Helmholtz equation, simple absorbing conditions of the form n iω were proposed as transmission
More informationNew Multigrid Solver Advances in TOPS
New Multigrid Solver Advances in TOPS R D Falgout 1, J Brannick 2, M Brezina 2, T Manteuffel 2 and S McCormick 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, P.O.
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationA Time-Parallel Implicit Method for Accelerating the Solution of Nonlinear Structural Dynamics Problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2008; 0:1 25 [Version: 2000/01/19 v2.0] A Time-Parallel Implicit Method for Accelerating the Solution of Nonlinear
More informationA Nonlinear PDE in Mathematical Finance
A Nonlinear PDE in Mathematical Finance Sergio Polidoro Dipartimento di Matematica, Università di Bologna, Piazza di Porta S. Donato 5, 40127 Bologna (Italy) polidoro@dm.unibo.it Summary. We study a non
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationCompact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems
Compact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems Daniel Ritter 1, Ulrich Rüde 1, Björn Gmeiner 1, Rochus Schmid 2 Copper Mountain, March 18th, 2013
More informationJournal of Computational and Applied Mathematics. Multigrid method for solving convection-diffusion problems with dominant convection
Journal of Computational and Applied Mathematics 226 (2009) 77 83 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationA 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 informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationA FINITE DIFFERENCE DOMAIN DECOMPOSITION ALGORITHM FOR NUMERICAL SOLUTION OF THE HEAT EQUATION
mathematics of computation volume 57, number 195 july 1991, pages 63-71 A FINITE DIFFERENCE DOMAIN DECOMPOSITION ALGORITHM FOR NUMERICAL SOLUTION OF THE HEAT EQUATION CLINT N. DAWSON, QIANG DU, AND TODD
More informationA CCD-ADI method for unsteady convection-diffusion equations
A CCD-ADI method for unsteady convection-diffusion equations Hai-Wei Sun, Leonard Z. Li Department of Mathematics, University of Macau, Macao Abstract With a combined compact difference scheme for the
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationParallel-in-Time Optimization with the General-Purpose XBraid Package
Parallel-in-Time Optimization with the General-Purpose XBraid Package Stefanie Günther, Jacob Schroder, Robert Falgout, Nicolas Gauger 7th Workshop on Parallel-in-Time methods Wednesday, May 3 rd, 2018
More informationAn improved collocation method based on deviation of the error for solving BBMB equation
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 6, No. 2, 2018, pp. 238-247 An improved collocation method based on deviation of the error for solving BBMB equation Reza
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 information