Introduction to Domain Decomposition Methods in the numerical approximation of PDEs. Luca Gerardo-Giorda
|
|
- Antonia Palmer
- 5 years ago
- Views:
Transcription
1 Introduction to Domain Decomposition Methods in the numerical approximation of PDEs Luca Gerardo-Giorda
2 Plan 0 Motivation Non-overlapping domain decomposition methods Multi-domain formulation Variational formulation of the multi-domain problem Iterative methods Non-overlapping methods at the finite dimensional level Finite dimensional formulation of the single domain problem Finite dimensional multi-domain formulation Iterative substructuring methods at the finite dimensional level 3 Overlapping Schwarz Methods Classical Schwarz methods Convergence of the Classical Schwarz method 4 Optimal and Optimized Schwarz Methods Optimal Schwarz Method Optimized Schwarz Methods 5 Applications of Optimized Schwarz Methods Maxwell equations Fluid-Structure Interaction problems in haemodynamics Heterogeneous models coupling in electrocardiology
3 3. Overlapping Schwarz Methods
4 When everything started (869) Question: Existence of harmonic function on non-lipschitz domains? Dirichlet principle. The solution of Laplace equation u = 0 on a bounded domain Ω with Dirichlet boundary conditions u = g on Ω is the infimum of the Dirichlet integral Ω v dx over all the functions v satisfying the boundary condition v = g on Ω. K. H. A. Schwarz (843-9) In 869 Schwarz introduces a method to prove that the infimum is actually attained.
5 Classical Alternating Schwarz Method The method is based on decomposing a general domain Ω into two overlapping domains with Lipschitz boundary Ω = Ω Ω.
6 Classical Alternating Schwarz Method The method is based on decomposing a general domain Ω into two overlapping domains with Lipschitz boundary Ω = Ω Ω. u = 0 in Ω u = g on Ω Ω u = u0 on Γ solve on the disk
7 Classical Alternating Schwarz Method The method is based on decomposing a general domain Ω into two overlapping domains with Lipschitz boundary Ω = Ω Ω. u = 0 in Ω u = g on Ω Ω u = u on Γ solve on the rectangle
8 Classical Alternating Schwarz Method The method is based on decomposing a general domain Ω into two overlapping domains with Lipschitz boundary Ω = Ω Ω. u = 0 in Ω u = g on Ω Ω u = u on Γ solve on the disk
9 Classical Alternating Schwarz Method The method is based on decomposing a general domain Ω into two overlapping domains with Lipschitz boundary Ω = Ω Ω. u = 0 in Ω u = g on Ω Ω u = u on Γ solve on the rectangle
10 Classical Alternating Schwarz Method The method is based on decomposing a general domain Ω into two overlapping domains with Lipschitz boundary Ω = Ω Ω. u n = 0 in Ω u n = 0 in Ω u n = g on Ω Ω u n = g on Ω Ω u n = un on Γ solve on the disk u n = un on Γ solve on the rectangle In 869 Schwarz proves convergence of the iterative procedure using the maximum principle.
11 Classical Schwarz Method Model problem u + η u = f on Ω = R lim x u(x) < + Domain decomposition: Ω =(, L) R Ω =(0, ) R Classical Schwarz Algorithm u n + η un = f in Ω u n + η un = f in Ω lim x u n (x) < lim x u n (x) < u n (L, y) = un (L, y) y R u n (0, y) = un (0, y) y R
12 Classical Schwarz Method: Convergence Fourier transform along the interface û(x, k) = e iky u(x, y) dy Error Equation ê n = û û n ê n = û û n (η + k xx )ê n = 0 x < L, k R (η + k xx )ê n = 0 x > 0, k R ê n (L, k) =ên (L, k) k R e n (0, k) =en (0, k) k R The solution is given by k ê n (x, k) =A (k) e +η (x L) + B (k) e k +η (x L) k ê n (x, k) =A (k) e +η x + B (k) e k +η x
13 Classical Schwarz Method: Convergence From boundedness assumption: ê n (x, k) =ên (L, k)e k +η (x L) ê n (x, k) =ên (0, k)e k +η x By induction: where ê n (0, k) = ρ n cs ê 0 (0, k) ê n (L, k) = ρ n cs ê 0 (L, k) ρ cs = ρ cs (k, L, η) =e k +η L < k R For the Classical Schwarz Alternating Method overlapping is mandatory for convergence
14 Classical Schwarz Method ! (k) k Fast convergence for high frequencies, slow convergence for low frequencies
15 Classical Schwarz Method
16 Classical Schwarz Method
17 Classical Schwarz Method
18 Classical Schwarz Method
19 Classical Schwarz Method
20 Classical Schwarz Method
21 Classical Schwarz Method
22 Classical Schwarz Method
23 Classical Schwarz Method
24 Classical Schwarz Method
25 Classical Schwarz Method
26 Nonoverlapping Schwarz Method Pierre-Louis Lions (990) On the Schwarz Alternating Method III: a variant for nonoverlapping subdomains Robin interface conditions n u n + α u n = n u n + α u n n u n + α u n = n u n + α u n ensure convergence also without overlap First of all, it is possible to replace the constants in the Robin conditions by two proportional functions on the interface, or even by local or nonlocal operators.
27 4. Optimal and Optimized Schwarz Methods
28 Modified Schwarz Method u n + η un = f in Ω ( x + S ) u n (L, y) =( x + S ) u n (L, y) y R u n + η un = f in Ω ( x + S ) u n (0, y) =( x + S ) u n (0, y) y R A Fourier Transform provides the error equation (η + k xx )ê n (x, k) =0 x < L, k R ( x + σ (k)) ê n (L, k) =( x + σ (k)) ê n (L, k) k R (η + k xx )ê n (x, k) =0 x < L, k R ( x + σ (k)) ê n (0, k) =( x + σ (k)) ê n (0, k) k R
29 Modified Schwarz Method The solutions are given by ê n (x, k) =ên (L, k) σ (k) k +η e σ (k)+ k +η ê n (x, k) =ên (0, k) σ (k)+ σ (k) k +η (x L) k +η k +η e k +η x By induction: ê n (0, k) = ρ n os ê 0 (0, k) ê n (L, k) = ρ n os ê 0 (L, k) where ρ os = ρ os (k, L, η, σ, σ )= σ (k) σ (k)+ k +η k +η σ (k)+ σ (k) k +η k +η e k +η L ρ os (k, L, η, σ, σ )= σ (k) σ (k)+ k +η k +η σ (k)+ σ (k) k +η k +η ρ cs(k, L, η)
30 Optimal Schwarz Method ρ os (k, L, η, σ, σ )= σ (k) σ (k)+ k +η k +η σ (k)+ σ (k) k +η k +η e k +η L σ (k) = k + η σ (k) = k + η = ρ os (k, L, η, σ, σ ) 0 Convergence in two iterations, independently from initial guess overlap size L (allows nonoverlapping decompositions) problem coefficient η Optimal but not viable k + η is a nonlocal operator û n x = k + η û n û n x = k + η û n k + η Steklov-Poincaré operator
31 Optimized Schwarz Method We need viable approximation of the exact operators: σ (k) = k + η σ (k) = k + η Low frequency approximations (based on Taylor expansion) σ (k) = η + η k + O(k 4 ) σ (k) = η η k + O(k 4 ) Zero-th order (Robin condition) ρ T 0 (k, L, η) = k + η η k + η + η ρ cs (k, L, η) Order (differential operator on the interface) k + η η ρ T (k, L, η) = η k k + η + η + η k ρ cs (k, L, η)
32 Optimized Schwarz Method (L = 0) We need viable approximation of the exact operators: σ (k) = k + η σ (k) = k + η Low frequency approximations (based on Taylor expansion) σ (k) = η + η k + O(k 4 ) σ (k) = η η k + O(k 4 ) Zero-th order (Robin condition) ρ T 0 (k, L, η) = k + η η k + η + η Order (differential operator on the interface) k + η η ρ T (k, L, η) = η k k + η + η + η k
33 Optimized Schwarz Method 0.9 Classical Schwarz Taylor order 0 Taylor order ! (k) k
34 Optimized Schwarz Method (L = 0) ! (k) Classical Schwarz Taylor order 0 Taylor order k
35 Optimized Schwarz Method Uniform approximation σ app (k) =p + q k σ app (k) = p q k k ρ = + η p q k k + η p q k ρ cs (k, L, η) k + η + p + q k k + η + p + q k Optimization as solution of a min-max problem Zero-th order conditions (q = q = 0) min p,p 0 max k min k kmax k + η p k + η p e k +η L k + η + p k + η + p Second order conditions min p j,q j 0 max k min k kmax k + η p q k k + η + p + q k k + η p q k k + η + p + q k e k +η L
36 Optimized Schwarz Method (p = p, q = q = 0) Theorem (M. J. Gander) For L > 0 and k max =, the solution p of the min-max problem is given by the unique root of the equation ρ Opt0 (k min, L, η, p )=ρ Opt0 ( k(l, η, p ), L, η, p ) k(l, η, p) = L(p + L(p η)) L For L = 0 and k max finite, the optimal parameter p is given by p = /4 k min + η kmax + η
37 Optimized Schwarz Method (L = 0) ! (k) Classical Schwarz Taylor order 0 Taylor order Optimized order k Equioscillation: ρ Opt0 (k min, L, η, p )=ρ Opt0 (k max, L, η, p )
38 Optimized Schwarz Method (p = p, q = q ) From: M.J. Gander, Optimized Schwarz Methods, SIAM J. Num. Anal. (006)
39 Optimized Schwarz Method
40 Optimized Schwarz Method
41 Optimized Schwarz Method
42 Optimized Schwarz Method
43 Optimized Schwarz Method
44 Algebraic formulation of the Optimized Schwarz Method The Dirichlet model problem we are considering is η u ν u = f in Ω u Ω = 0 The sequential version of the OSM reads as follows Given λ 0 = nt ν u0 + α u 0 on Γ Solve for p 0. η u p+ ν u p+ = f in Ω n T ν up+ + α u p+ = λ p on Γ. λ p+ = n T ν up+ + α u p+ 3. η u p+ ν u p+ = f in Ω n T ν up+ + α u p+ = λ p+ on Γ 4. λ p+ = n T ν up+ + α u p+ Remark The normal derivatives in the interface conditions emerge naturally in the variational formulation The traces in the interface conditions have to be included by adding a term in the variational formulation
45 Algebraic formulation of the Optimized Schwarz Method Let us consider a non-overlapping decomposition. The bilinear forms within Ω j (j =, ) a j (u, v) = η Ω j uv M j Mass matrix + Ω j ν u v S j Stiffness matrix (n T ν u) v + Ω j Γ (n T ν u) v yield subdomain matrices B j = η M j + S j that already include the normal derivative on the interface Γ.
46 Algebraic formulation of the Optimized Schwarz Method Let us consider a non-overlapping decomposition. The bilinear forms within Ω j (j =, ) a j (u, v) = η Ω j uv M j Mass matrix + Ω j ν u v S j Stiffness matrix yield subdomain matrices B j = η M j + S j that already include the normal derivative on the interface Γ. To include the trace component we add a mass matrix defined on the interfaces nodes. [M Γ ] ij = Γ φ (Γ) j Γ φ(γ) i Γ
47 Algebraic formulation of the Optimized Schwarz Method We also need the restriction matrices, from the nodes in Ω j to the interface Γ. Let N j and N Γ be the number of internal and interface degrees of freedom for Ω j (j =, ), respectively. After a proper reordering of the unknowns, they are the Boolean operators R jγ =[0 Id ]= N Γ N j... N Γ N Γ Their transposes are the zero-extension matrices from the interface Γ to Ω j. For u = uj u Γ : R jγ R T jγ u Γ = u Γ R T jγ R jγ uj u Γ = 0 uγ We can now write the subdomain matrices associated to Robin interface conditions, as A = B + α R T Γ M Γ R Γ A = B + α R T Γ M Γ R Γ
48 Algebraic formulation of the Optimized Schwarz Method We are now in the position of stating the fully discrete formulation of the Optimized Schwarz Method in sequential iterative form. Given λ 0 Solve for p 0 A u p+ = f + R T Γ λp λ p+ = R Γ (f B u p+ )+α M Γ R Γ u p+ 3 A u p+ = f + R T Γ λp+ 4 λ p+ = R Γ (f B u p+ )+α M Γ R Γ u p+
49 Optimized Schwarz Method as a preconditioner The parallel formulation of the Optimized Schwarz Method reads: A u p+ = f + R T Γ λp A u p+ = f + R T Γ λp λ p+ = R Γ (f B u p+ )+α M Γ R Γ u p+ λ p+ = R Γ (f B u p+ )+α M Γ R Γ u p+ Recalling the definition of A and A : A = B + α R T Γ M Γ R Γ A = B + α R T Γ M Γ R Γ, from the linear systems above we have u p+ = A f + A RT Γ λp B u p+ = f + R T Γ λp α R T Γ M Γ R Γ u p+ u p+ = A f + A RT Γ λp B u p+ = f + R T Γ λp α R T Γ M Γ R Γ u p+ We can now replace those quantities in the definitions of the interface variables.
50 Optimized Schwarz Method as a preconditioner u p+ = A f + A RT Γ λp B u p+ = f + R T Γ λp α R T Γ M Γ R Γ u p+ λ p+ = R Γ (f B u p+ )+α M Γ R Γ u p+ B u p+ = f + R T Γ λp α R T Γ M Γ R Γ u p+ = R Γ (f f R T Γ λp + α R T Γ M Γ R Γ u p+ )+α M Γ R Γ u p+ = R Γ R T Γ Id Γ p λ + α R Γ R T Γ Id Γ p+ M Γ R Γu + α M Γ R Γ u p+ = λ p +(α + α ) M Γ R Γ u p+ u p+ = A (f + R T Γ λp ) = λ p +(α + α ) M Γ R Γ A (f + R T Γ λp ) = Id +(α + α ) M Γ R Γ A RT Γ λ p +(α + α ) M Γ R Γ A f
51 Optimized Schwarz Method as a preconditioner u p+ = A f + A RT Γ λp B u p+ = f + R T Γ λp α R T Γ M Γ R Γ u p+ λ p+ = R Γ (f B u p+ )+α M Γ R Γ u p+ B u p+ = f + R T Γ λp α R T Γ M Γ R Γ u p+ = R Γ (f f R T Γ λp + α R T Γ M Γ R Γ u p+ )+α M Γ R Γ u p+ = R Γ R T Γ Id Γ p λ + α R Γ R T Γ Id Γ p+ M Γ R Γu + α M Γ R Γ u p+ = λ p +(α + α ) M Γ R Γ u p+ u p+ = A (f + R T Γ λp ) = λ p +(α + α ) M Γ R Γ A (f + R T Γ λp ) = Id +(α + α ) M Γ R Γ A RT Γ λ p +(α + α ) M Γ R Γ A f
52 Optimized Schwarz Method as a preconditioner Solving for λ p+ and λ p+ amounts to solve: λ p+ = λ p+ = Id +(α + α ) M Γ R Γ A RT Γ λ p +(α + α ) M Γ R Γ A f Id +(α + α ) M Γ R Γ A RT Γ λ p +(α + α ) M Γ R Γ A f Such algorithm is actually a fixed point iteration procedure to solve I I (α + α ) M Γ Σ I (α + α ) M Γ Σ I λ λ =(α + α ) R,Γ A f R,Γ A f where Σ = R,Γ A RT,Γ Σ = R,Γ A RT,Γ.
53 Classical Schwarz Method as a preconditioner Despite not being convergent without overlap in the iterative form, the Classical Schwarz can be used as a preconditioner for a Krylov acceleration procedure to solve the interface equation. In this case, the matrices associated to the Dirichlet Aii A A j = iγ problems on the interface of Ω j (j =, ) are given by 0 Id The parallel version of the Classical Schwarz Method reads as follows.. A u p+ = f + R T,Γ λp A u p+ = f + R T,Γ λp. λ p+ = R,Γ u p+ λ p+ = R,Γ u p+ Usiing. we solve for λ p+ and λ p+ : λ p+ = R,Γ A f + A RT,Γ λp λ p+ = R,Γ A f + A RT,Γ λp Fixed point iteration procedure to solve I R,Γ A R T,Γ R,Γ A R T,Γ I λ λ = R,Γ A f R,Γ A f
54 Numerical comparison 0 Classical Schwarz Classical Schwarz with Krylov Multigrid 0 Optimized Schwarz Optimized Schwarz with Krylov
55 Optimized Schwarz Method The optimized interface condition is given by n u n + Sopt un = n un + S opt un n u n + Sopt un = n un + S opt un where S opt j is a viable approximation of the Steklov-Poincaré operator in Ω j. The choice of the optimal operator is PDE-dependent
56 Optimized Schwarz Method The optimized interface condition is given by n u n + Sopt un = n un + S opt un n u n + Sopt un = n un + S opt un where S opt j is a viable approximation of the Steklov-Poincaré operator in Ω j. PDE Optimized interface operator u + ηu = f () n u +[p (τ )(τ )] u ν u + b u = f () n u + n b p ν u u + ω u = f (3) n u + [(p + iq) (α + iβ)(τ )(τ )] u curl curl u ω u = f (4) (curl u) n + p+iq curl Γ curl Γ ω I (n u n) () : M.J. Gander, SIAMJ.NUM. ANAL. (006) (3) : M.J. Gander, F. Magoulès, F. Nataf, SIAMJ.SCI. COMPUT. (00) () : O. Dubois, DDMIN SCIENCE AND ENGINEERING, (007) (4) : A. Alonso-Rodrìguez, L.GG, SIAMJ.SCI. COMPUT. (006)
57 How to choose parameters in general
58 General settings of an Optimized Schwarz Method When the operator is elliptic, the algorithm is based on a combination of Fluxes and Traces PDE Optimized interface operator u + ηu = f n u +[p (τ )(τ )] u ν u + b u = f n u + n b p ν u u + ω u = f n u + [(p + iq) (α + iβ)(τ )(τ )] u curl curl u ω u = f (curl u) n + p+iq curl Γ curl Γ ω I (n u n) When the operator is hyperbolic, the algorithm is based on a combination of incoming characteristics and outcoming characteristics
59 5. Applications of Optimized Schwarz Methods
60 Optimized Schwarz Method for the full Maxwell system Time-Harmonic Maxwell Equations: iωεe + curl H σe = J iωµh + curl E = 0 Combination of incoming and outcoming characteristic variables. Solve for n =,,... Dirichlet coupling conditions Hyperbolic problems Characteristic variables Maxwell Equation Impedance conditions B nj (E, H) :=n j ε µ E + n j (H n j ). iωεe,n + curl H,n σe,n = J in Ω iωµh,n + curl E,n = 0 in Ω (Bn + S Bn )(E,n, H,n ) = (Bn + S Bn )(E,n, H,n ) on Γ iωεe,n + curl H,n σe,n = J in Ω iωµh + curl E = 0 in Ω (Bn + S Bn )(E,n, H,n ) = (Bn + S Bn )(E,n, H,n ) on Γ, S, S : Tangential operators Possibly pseudodifferential Based on Transparent b.c.
61 Transparent Boundary Conditions Fourier Analysis: Ω = R 3 Ω =(, L) R, Ω =(0, ) R J, compactly supported in Ω, lim r r H, n E, = 0 ω = ω µε, k = k y + k z (tangential frequency) Theorem (V. Dolean, M.J. Gander, L.GG, SIAM J. Sci. Comp. 009) Let the operators S l,l=,, have the Fourier symbol k F(S l )= ( y k z k y k z k ω + i ω) k y k z kz ky. Then, the solution of Maxwell s equations in Ω, with the boundary conditions (Bn + S Bn )(E (0, y, z), H (0, y, z)) = 0 (y, z) R (Bn + S Bn )(E (L, y, z), H (L, y, z)) = 0 (y, z) R coincides with the restriction on Ω, of the solution of Maxwell s equations in R 3.
62 Optimized Schwarz Method for Maxwell equations Approximate the nonlocal part of the symbols. k F(S l ):=γ y k z k y k z l k y k z kz ky γ l C(k z, k y ) # γ γ ρ # (k) 0 P(k) iω εµ P(k)+iω εµ Λ(k, L) s iω εµ k s+iω εµ P(k) s P(k)+s Λ(k, L) 3 k ω εµ+iω εµs P(k) iω εµ P(k)+iω εµ ρ ( k ) s 4 iω εµ k s +iω εµ s iω εµ k s +iω εµ P(k) s P(k)+s P(k) s P(k)+s Λ(k, L) 5 k ω εµ+iω εµs k ω εµ+iω εµs P(k) iω εµ P(k)+iω εµ ρ 4 ( k ) P(k) = k ω εµ Λ(k, L) =e P(k)L
63 Optimized Schwarz Method for Maxwell equations: asymptotics with overlap, L = h without overlap, L = 0 Case ρ parameters ρ parameters k + ω h none none C 6 ω h 3 p = C ω h 3 3 C CC 4 4 ω h p = ω C h 3 (k + ω ) 4 5 C 0 ω h h 5 (k+ ω ) 0 h 5 3 p = (k + ω ) h 3 p = p = C 5 ω 7 5 h C 5 ω 6 5 h p = (k + ω ) h 5 p = (k + ω ) h 5,, (k + ω ) 4 C 8 ω C 4 C(k h p = + ω ) C h h (k + ω ) C h 4 3 p = C 8 ω C 4 h 4 3 p = C 8 ω C 4 3 h 4 p = (k + ω ) h 4 4, 3 8 C 4 3 p = (k + ω ) 8 C 4 3 h 4,
64 Numerical results Asymptotics!" # -,./0&*&%*+3'(4'+'5*(&$'*3'(6)00%+4)64*(%&$/7&%/'!$)(/*+%),' ),'! ),'# ),'8 ),'9 ),': ;<$!! = ;<$!!>8 = ;<$!!>8 = ;<$!!>: = ;<$!!>: =!" # -,./0&*&%*+3'(4'+'5*(&$'+*+!*3'(6)00%+4)64*(%&$/7&%/'!$)(/*+%),' ),'# ),'8 ),'9 ),': ;<$!!=# > ;<$!!=# > ;<$!!=9 > ;<$!!=9 %&'()&%*+, %&'()&%*+,!"!!"!# $!"!# $ 3d tests: Ω =(0, ) 3, ε = µ =, σ = 0 Ω =(0, β) (0, ) Ω =(α, ) (0, ) (0 < α β < ) Uniform mesh, FV scheme (INRIA, Sophia Antipolis) it S (it GS ) with overlap, L = h without overlap, L = 0 h /8 /6 /30 /8 /6 /30 Case 9(3) 9(7) 46() (93) (40) (0) Case 4() 9(4) 3(6) 48(9) 69(36) 98(48) Case 3 6() 8(4) (6) 65(35) 80(4) 66(55) Case 4 5(3) 9(5) (7) 38(8) 60(33) 04(39) Case 5 6(3) 8(4) (6) 70(36) 80(4) 76(55) it S : iterative version it GS : accelerated (GMRES) version
65 Chicken heating in a Whirlpool TC4 Microwave Ω =[0, 0.3] [0, 0.36] [0, 0.0]m Magnetron size m, metallic b.c. ε C = farad σ C = 3 0 m siemens m farad m ε A = σ A = 0 siemens m µ C = µ A = 4π 0 ω = π.45 GHz 7 henry m subdomains h = 0.005m Electrical Field Magnetic Field
66 Fluid Structure Interaction Problem Fluid-structure problem. Find the fluid velocity u, pressure p and the structure displacement η such that D A u ρ f + ρ f ((u w) )u T f = f f in Ω t (0, T ), Dt f u = 0 in Ω t (0, T ), f η ρ s t T s = f s in Ω 0 s (0, T ), () u = η on Σ t (0, T ), t T f n T s n = 0 on Σ t (0, T ), Geometry problem. Given η Σ t, find a map A : Ω 0 f Ω t f A t (x 0 )=x 0 + Ext( η Σ 0 ), such that Ω t f = A t (Ω 0 f ).
67 Time discretization and Fluid Structure interface treatment Given Ω f, u, w, u n, η n and η n, find u n+, p n+, and η n+ such that ρ f δ t u n+ + ρ f ((u w ) )u n+ T n+ f = f n+ f in Ω f, u n+ = 0 in Ω f, ρ s δ tt η n+ n+ n+ T s = f s in Ω 0 s, u n+ = δ t η n+ on Σ, T n+ f n = T n+ s n on Σ, with suitable boundary conditions for the artificial sections Γ t j, Γ j,s. Here: δ t w n+ = w n+ w n t δ tt w n+ = w n+ w n + w n t
68 Partitioned procedures Dirichlet-Neumann partitioned algorithms solve Fluid with the Structure velocity as Dirichlet b.c. solve Structure with the Fluid normal stress as Neumann b.c. Slow convergence and strong relaxation needed in the presence of large added mass effect. Robin-Robin partitioned algorithms At each iteration m: find u m, p m and η m such that Fluid Structure ρ f t u m + ρ f ((u w ) )u m T f,m = f f + ρ f t un in Ω f, u m = 0 in Ω f, η m η n α f u m T f,m n = α f T s,m n, on Σ, t ρ s t η m T s,m = f s + ρ s t ( ηn η n ) in Ω 0 s, α s t η m + T s,mn = α s t ηn + α s u m + T f,m n, on Σ The convergence is driven by the choice of α f and α s
69 . Potential Flow - Generalized String (P/GS) ρ f δ t u + p = 0 in Ω f = {x < 0} R, u = 0 in Ω f, u n = δ t η on Σ = {x = 0} R, p = ρ s H s δ tt η + βh s η GH s yy η, on Σ Robin Interface conditions: linear combinations with (α f, ), (α s, ) αf α f t δ t u x p = t ρ sh s t βh s η + GH s yy η + F (ux n, ηn, η n ), ρs H s t + βh s + α s η GH s yy η = α s t δ t u x + p + F (u n, η n, η n ), t here: u x = u n
70 . Potential Flow - Generalized String (P/GS) Since on Σ x p = ρ f δ t u x we can handle the interface condition in terms of the sole pressure p By applying the divergence operator and the Fourier Transform in the y direction, we get xx p m+ + k p m+ = 0 in (, 0) p m+ + α f t ρ f x p m+ = ρs H s t ρs H s t + βh s + GH s k + α s t + βh s + GH s k α f η m, for x = 0 t η m = p m α s t ρ f x p m for x = 0. We have p m+ (x, k) =A m+ (k) e kx where A m+ (k) is determined by the interface conditions x p m+ = A m+ (k) ke kx
71 . Potential Flow - Generalized String (P/GS) Then ρf + α f tk A m+ (k) = ρs H s t ρ f ρs H s t + βh s + GH s k + α s t + βh s + GH s k α f η m, for x = 0, t η m = ρ f α s tk ρ f A m (k), for x = 0. We thus have A m+ (k) = ρ(k) A m (k) The reduction factor of the iterative algorithm is independent of the iteration and is given by ρ P/GS α (k) = f F GS (k) α s + F GS (k) α s F P (k) α f + F P (k) Here: F GS (k) = ρ sh s t The algorithm converges provided ρs H s 0 α f < + βh s t + GH s tkmax t + βh s t + GH s tk F P (k) = ρ f t k 0 α s < ρ f α f + α s > 0 t k max
72 . Potential Flow - Generalized String (P/GS): parameter optimization We seek for α f and α s such that ρ P/GS (k) is minimal, in A {α(k) =α M + α K k, α M, α K 0} corresponding to a generalized Robin boundary condition: α M u α K Σ u + T f n = f (η) We have: ρ P/GS (k) = α f F GS (k) α s + F GS (k) α s F P (k) α f + F P (k) = 0 for α opt f Thus = F GS (k) = ρ sh s t + βh s t + GH s t k α GS f,m = ρ sh s t + βh s t = α heur f α GS f,k = GH s t
73 . Potential Flow - Generalized String (P/GS): parameter optimization α opt s = F P (k) = ρ f t k = integral boundary condition. Find α s through a minimization problem α opt s (k) = argmin max sup α s (k) A k [0,k max ] α f A α s F P (k) α f + F P (k) = argmin α s (k) A max k [0,k max ] α s t k ρ f If α s,k = 0 then ρ P/GS is unbounded for k = α P s = α s,m We have max k [0,k max ] α s t k ρ f = if α s,m ρ f t k max > otherwise = α P s = ρ f t k max
74 . Stokes - Generalized String (S/GS) ρ f δ t u u + p = 0 in Ω f = {x < 0} R, u = 0 in Ω f, u n = δ t η on Σ = {x = 0} R, p = ρ s H s δ tt η + βh s η GH s yy η, on Σ Reduction factor: ρ S/GS (k) = α f F GS (k) α s + F GS (k) α s F S (k) α f + F S (k) where F S (k) = µγ f (k)(k + γ f (k)) k γ f (k) = ρf µ t + k α opt = F f GS (k) as before = α GS f = ρ sh s t + βh s t + GH s t k α opt s from minimization α S s = µ t ρf t k + µ t (k ) k + ρ f + µ t (k ).
75 3. Stokes - Linear Incompressible Elasticity (S/LIE) ρ f δ t u u + p = 0 in Ω f = {x < 0} R, u = 0 in Ω f, ρ s δ tt η λη + βη + χ = 0 in Ω s =(0, H s ) R, η = 0 in Ω s Reduction factor: ρ S/LIE (k) = α f F LIE (k) α s + F LIE (k) α s F S (k) α f + F S (k) where: ε(k) = kγ s H k +γs s + O(Hs ), F LIE (k) = ε(k)+ t δ(k) = (ε(k)+) γ sh s + O(H δ(k) λγ s γ s H s 6 s ) γ s (k) = ρs λ t + β λ + k. For thin structures there is no significant improvement with respect to the GS model.
76 Numerical results I: rectangular domain Ω f 0 = 6 cm, Ωs 0 = 6 H s cm, ρ f = g/cm 3 ρ s =. g/cm 3 µ = poise E = dyne/cm 3 H s = 0. cm t = 0 3 s ν = 0.3, R = 0.5, β = E ( ν )R P(t) =0 4 dyne/cm for t s Average number of iterations and (in brackets) relative CPU time α f = α s = 0 0 α P s α S s Ref (.00) 5.37 (0.9) 4.6 (0.79) 4.50 (0.77) ρ s = (.00) 5.6 (0.94) 4.87 (0.8) 4.87 (0.8) ρ s = (.00) 5.75 (0.96) 5.00 (0.83) 4.87 (0.8) H s = (.00) 5.87 (0.87) 5.00 (0.74) 5.00 (0.74) H s = (.00) 7.87 (0.73) 7.6 (0.7) 5.75 (0.53) t = (.00) 5.69 (0.8) 4.94 (0.7) 4.87 (0.70) t = (.00) 5.56 (0.80) 4.8 (0.69) 4.4 (0.64) E = (.00) 5.6 (0.86) 4.87 (0.75) 4.87 (0.75) E = (.00) 5.50 (0.73) 4.37 (0.58) 5.00 (0.67) ν = (.00) 4.87 (0.87) 4.00 (0.7) 4.00 (0.7) α heur f α GS f α GS f α GS f Matlab - Fluid: P -bubble/p - Structure: P
77 Numerical results II: simplified carotid Average number of iterations and (in brackets) relative CPU time α f = α s = 0 0 α P s α S s Ref (.00) 9.9 (0.76) 7.00 (0.54) 6.75 (0.5) H s = (.00) 4.08 (0.6) 3.50 (0.58).67 (0.54) α heur f α GS f α GS f α GS f Matlab - Fluid: P -bubble/p - Structure: P
78 The Bidomain model for cardiac electrophysiology Complex problem: coexistence of intra- and extra-cellular media Ohm s Law charge conservation div (D i u i )=χ(c m t u + I ion (u, w)) in Ω (0, T ) div (D e u e )= χ(c m t u + I ion (u, w)) in Ω (0, T ) u = u i u e in Ω (0, T ) t w = R(u, w) in Ω (0, T ) u i intracellular potential u e extracellular potential w recovery variable n T D i,e u i,e = 0 in Ω (0, T ) u(x, 0) =u 0 (x), w(x, 0) =w 0 (x), in Ω Unique Up to a constant = Closing condition Ω u e dx = 0 MAIN ISSUE: Computational complexity
79 Simplifying Bidomain: Monodomain Model Bidomain in terms of u and u e χc m u t λdi + λ u [D i u +(D i + D e ) u e ]= I app. λdi D e + λ u e + I ion (u, w, c) =I app Monodomain λd i = D e D M (x) = λ + λ D i (x) χc m u t (D M (x) u) +I ion (u, w) =I app in Ω (0, T ) t w = R(u, w) in Ω (0, T ) n T D i u = 0 in Ω (0, T ) u(x, 0) =u 0 (x), w(x, 0) =w 0 (x), in Ω Pros Easier to solve Cons Less reliable... Model adaptivity: solve Bidomain only where needed = heterogeneous coupling
80 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Bidomain/Bidomain coupling n T B D i u i,b + αu i,b = n T B D i u i,m + αu i,m n T B D e u e,b + βu e,b = n T B D e u e,m + βu e,m Non-symmetric formulation n T B λd i + λ ( u B + u e,b ) n T B = n T B D e + λ u e,b λd i + λ ( u M + u e,m ) n T B + λα + λ u B + λα β + λ u e,b D e + λ u e,m + λα + λ u M + λα β + λ u e,m n T B D i ( u B + u e,b )+n T B D e u e,b + αu B +(α + β)u e,b = n T B D i ( u M + u e,m )+n T B D e u e,m + αu M +(α + β)u e,m
81 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Bidomain/Bidomain coupling n T B D i u i,b + αu i,b = n T B D i u i,m + αu i,m n T B D e u e,b + βu e,b = n T B D e u e,m + βu e,m Dimensional mismatch Non-symmetric formulation n T B λd i + λ ( u B + u e,b ) n T B = n T B D e + λ u e,b λd i + λ ( u M + u e,m ) n T B + λα + λ u B + λα β + λ u e,b D e + λ u e,m + λα + λ u M + λα β + λ u e,m n T B D i ( u B + u e,b )+n T B D e u e,b + αu B +(α + β)u e,b = n T B D i ( u M + u e,m )+n T B D e u e,m + αu M +(α + β)u e,m L.GG, M. Perego, A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology, MAN, DOI: 0.05/man/00057, Published online: 0 Aug 00
82 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Bidomain/Bidomain coupling n T B D i u i,b + αu i,b = n T B D i u i,m + αu i,m n T B D e u e,b + βu e,b = n T B D e u e,m + βu e,m Dimensional mismatch β = λα Non-symmetric formulation n T B λd i + λ ( u B + u e,b ) n T B = n T B D e + λ u e,b λd i + λ ( u M + u e,m ) n T B + λα + λ u B + λα β + λ u e,b D e + λ u e,m + λα + λ u M + λα β + λ u e,m n T B D i ( u B + u e,b )+n T B D e u e,b + αu B +(α + β)u e,b = n T B D i ( u M + u e,m )+n T B D e u e,m + αu M +(α + β)u e,m L.GG, M. Perego, A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology, MAN, DOI: 0.05/man/00057, Published online: 0 Aug 00
83 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Bidomain/Bidomain coupling n T B D i u i,b + αu i,b = n T B D i u i,m + αu i,m n T B D e u e,b + βu e,b = n T B D e u e,m + βu e,m Dimensional mismatch β = λα Non-symmetric formulation n T B λd i + λ ( u B + u e,b ) n T B = n T B D e + λ u e,b + λα + λ u B λd i + λ ( u M + u e,m ) n T B D e + λ u e,m + λα + λ u M n T B D i ( u B + u e,b )+n T B D e u e,b + αu B +( + λ)αu e,b = n T B D i ( u M + u e,m )+n T B D e u e,m + αu M +( + λ)αu e,m L.GG, M. Perego, A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology, MAN, DOI: 0.05/man/00057, Published online: 0 Aug 00
84 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Bidomain/Bidomain coupling n T B D i u i,b + αu i,b = n T B D i u i,m + αu i,m n T B D e u e,b + βu e,b = n T B D e u e,m + βu e,m Dimensional mismatch β = λα 3 λd i = D e Non-symmetric formulation n T B λd i + λ ( u B + u e,b ) n T B = n T B D e + λ u e,b + λα + λ u B λd i + λ ( u M + u e,m ) n T B D e + λ u e,m + λα + λ u M n T B D i ( u B + u e,b )+n T B D e u e,b + αu B +( + λ)αu e,b = n T B D i ( u M + u e,m )+n T B D e u e,m + αu M +( + λ)αu e,m L.GG, M. Perego, A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology, MAN, DOI: 0.05/man/00057, Published online: 0 Aug 00
85 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Bidomain/Bidomain coupling n T B D i u i,b + αu i,b = n T B D i u i,m + αu i,m n T B D e u e,b + βu e,b = n T B D e u e,m + βu e,m Dimensional mismatch β = λα 3 λd i = D e Non-symmetric formulation n T B λd i + λ ( u B + u e,b ) n T B D e + λ u e,b + λα + λ u B = n T B λd i + λ u M + λα + λ u M n T B D i ( u B + u e,b )+n T B D e u e,b + αu B +( + λ)αu e,b = n T B D i ( u M + u e,m )+n T B D e u e,m + αu M +( + λ)αu e,m L.GG, M. Perego, A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology, MAN, DOI: 0.05/man/00057, Published online: 0 Aug 00
86 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Bidomain/Bidomain coupling n T B D i u i,b + αu i,b = n T B D i u i,m + αu i,m n T B D e u e,b + βu e,b = n T B D e u e,m + βu e,m Dimensional mismatch β = λα 3 λd i = D e 4 Approx of u e,m : Non-symmetric formulation n T B λd i + λ ( u B + u e,b ) n T B D e + λ u e,b + λα + λ u B = n T B λd i + λ u M + λα + λ u M n T B D i ( u B + u e,b )+n T B D e u e,b + αu B +( + λ)αu e,b = n T B D i ( u M + u e,m )+n T B D e u e,m + αu M +( + λ)αu e,m L.GG, M. Perego, A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology, MAN, DOI: 0.05/man/00057, Published online: 0 Aug 00
87 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Bidomain/Bidomain coupling n T B D i u i,b + αu i,b = n T B D i u i,m + αu i,m n T B D e u e,b + βu e,b = n T B D e u e,m + βu e,m Dimensional mismatch β = λα 3 λd i = D e 4 Approx of u e,m : I app = 0, n T D e u e,m = 0, n T D i ( u M + u e,m )=0: u e,m = u M + K + λ Non-symmetric formulation n T B λd i + λ ( u B + u e,b ) n T B = n T B λd i + λ u M D e + λ u e,b + λα + λ u B + λα + λ u M n T B D i ( u B + u e,b )+n T B D e u e,b + αu B +( + λ)αu e,b = n T B D i ( u M + u e,m )+n T B D e u e,m + αu M +( + λ)αu e,m L.GG, M. Perego, A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology, MAN, DOI: 0.05/man/00057, Published online: 0 Aug 00
88 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Bidomain/Bidomain coupling n T B D i u i,b + αu i,b = n T B D i u i,m + αu i,m n T B D e u e,b + βu e,b = n T B D e u e,m + βu e,m Dimensional mismatch β = λα 3 λd i = D e 4 Approx of u e,m : I app = 0, n T D e u e,m = 0, n T D i ( u M + u e,m )=0: u e,m = u M + K + λ K = u rest + λ u e,m = 0 u M = u rest Non-symmetric formulation n T B λd i + λ ( u B + u e,b ) n T B = n T B λd i + λ u M D e + λ u e,b + λα + λ u B + λα + λ u M n T B D i ( u B + u e,b )+n T B D e u e,b + αu B +( + λ)αu e,b = n T B D i ( u M + u e,m )+n T B D e u e,m + αu M +( + λ)αu e,m L.GG, M. Perego, A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology, MAN, DOI: 0.05/man/00057, Published online: 0 Aug 00
89 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Bidomain/Bidomain coupling n T B D i u i,b + αu i,b = n T B D i u i,m + αu i,m n T B D e u e,b + βu e,b = n T B D e u e,m + βu e,m Dimensional mismatch β = λα 3 λd i = D e 4 Approx of u e,m : I app = 0, n T D e u e,m = 0, n T D i ( u M + u e,m )=0: u e,m = (u M u rest ) + λ Non-symmetric formulation n T B λd i + λ ( u B + u e,b ) n T B = n T B λd i + λ u M D e + λ u e,b + λα + λ u B + λα + λ u M n T B D i ( u B + u e,b )+n T B D e u e,b + αu B +( + λ)αu e,b = αu rest L.GG, M. Perego, A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology, MAN, DOI: 0.05/man/00057, Published online: 0 Aug 00
90 Bidomain-Monodomain coupling: Optimized Schwarz Coupling Generalized Robin Coupling conditions on Γ n T B λd i + λ ( u B + u e ) n T B D e + λ u e + λα + λ u B = n T B λd i + λ u M + λα + λ u M n T B D i ( u B + u e )+n T B D e u e + αu B +( + λ)αu e = αu rest Convergence For α > 0, the OSA converges for every initial guess. Proof. Fourier analysis: the interface error reduces as e p+ Γ with ρ(α) <. ρ(α) e p Γ! k h =.e 3! k h =.e L.GG, M. Perego, A. Veneziani, Optimized Schwarz coupling of Bidomain and Monodomain models in electrocardiology, MAN, DOI: 0.05/man/00057, Published online: 0 Aug 00
91 Bidomain-Monodomain coupling: Convergence Convergence (L.GG, M. Perego, A. Veneziani, MAN, 0) For α > 0, the Optimized Schwarz Algorithm converges for every initial guess, and for every Fourier mode k its reduction factor is independent of the iteration and is given by α σi l µm (k) α [ Ψ(k, α)] σ ρ(k, α) = α + σi l µm (k) i l η Ψ(k, α) (k ) η + (k ) α [ Ψ(k, α)] + σi l η Ψ(k, α) (k, ) η + (k ) where µ M (k) = χc m t σ l i + σ l e + σt i σt e σ t i + σ t e σ l i + σ l e σ l i σl e k, Ψ(k, α) = (σt i + σe t )k (σi l + σe l ) η (k) (σi t + σe)k t (σi l + σe) l η + (k) α + σl i η + (k) α + σi l η (k), and η ± (k) = χcm σi l t ± σ l i + σi t k + σe l χcm t χcm t + σi t k σe l + σ t ek χcm t + σ t ek + 4 σ l i σl e χcm. t
92 Bidomain-Monodomain coupling: Parameter Optimization Exact solution: ρ(α exact ) 0 α exact (k) = σi l µm (k) =σi l χc m t σ l i + σ l e + σt i σt e σ t i + σ t e σ l i + σ l e σ l i σl e k unpractical α exact (k) : α exact (k) = σ l i η + (k) Ψ(k, α exact (k)) η (k) Ψ(k, α exact (k)) η + (k) Ideal solution (still difficult to compute) min α R + max k [0,k max ] ρ(k, α) Viable choice: low frequency approximation α opt = α exact (0) α = σ l i χc m t σ l i +σl e σ l i σl e
Comparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping
Comparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping M. El Bouajaji 1, V. Dolean 2, M. J. Gander 3 and S. Lanteri 1 1 Introduction Transmission conditions
More informationFrom the invention of the Schwarz method to the Best Current Methods for Oscillatory Problems: Part 2
and From the invention of the Schwarz method to the Best Current Methods for Oscillatory Problems: Part 2 martin.gander@unige.ch University of Geneva Woudschoten, October 2014 Are There Two Equations?
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 informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationTHE OPTIMISED SCHWARZ METHOD AND THE TWO-LAGRANGE MULTIPLIER METHOD FOR HETEROGENEOUS PROBLEMS. by Neil Greer
THE OPTIMISED SCHWARZ METHOD AND THE TWO-LAGRANGE MULTIPLIER METHOD FOR HETEROGENEOUS PROBLEMS by Neil Greer Submitted for the degree of Doctor of Philosophy Department of Mathematics School of Mathematical
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 informationDG discretization of optimized Schwarz methods for Maxwell s equations
DG discretization of optimized Schwarz methods for Maxwell s equations Mohamed El Bouajaji, Victorita Dolean,3, Martin J. Gander 3, Stéphane Lanteri, and Ronan Perrussel 4 Introduction In the last decades,
More informationOptimized Schwarz Methods for Maxwell equations
Optimized Schwarz Methods for Maxwell equations Victorita Dolean, Martin Gander, Luca Gerardo-Giorda To cite this version: Victorita Dolean, Martin Gander, Luca Gerardo-Giorda. Optimized Schwarz Methods
More informationDirichlet-Neumann and Neumann-Neumann Methods
Dirichlet-Neumann and Neumann-Neumann Methods Felix Kwok Hong Kong Baptist University Introductory Domain Decomposition Short Course DD25, Memorial University of Newfoundland July 22, 2018 Outline Methods
More informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationOPTIMIZED SCHWARZ METHODS
SIAM J. NUMER. ANAL. Vol. 44, No., pp. 699 31 c 006 Society for Industrial and Applied Mathematics OPTIMIZED SCHWARZ METHODS MARTIN J. GANDER Abstract. Optimized Schwarz methods are a new class of Schwarz
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationc 2012 Society for Industrial and Applied Mathematics
SIAM J SCI COMPUT Vol 34, No 4, pp A48 A7 c Society for Industrial and Applied Mathematics OPTIMIZED SCHWARZ METHODS FOR THE TIME-HARMONIC MAXWELL EQUATIONS WITH DAMPING M EL BOUAJAJI, V DOLEAN, M J GANDER,
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationOptimized Schwarz methods for the time-harmonic Maxwell equations with damping
Optimized Schwarz methods for the time-harmonic Maxwell equations with damping Mohamed El Bouajaji, Victorita Dolean, Martin Gander, Stephane Lanteri To cite this version: Mohamed El Bouajaji, Victorita
More informationCoupling Stokes and Darcy equations: modeling and numerical methods
Coupling Stokes and Darcy equations: modeling and numerical methods Marco Discacciati NM2PorousMedia 24 Dubrovnik, September 29, 24 Acknowledgment: European Union Seventh Framework Programme (FP7/27-23),
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More informationEfficient domain decomposition methods for the time-harmonic Maxwell equations
Efficient domain decomposition methods for the time-harmonic Maxwell equations Marcella Bonazzoli 1, Victorita Dolean 2, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier 4 1 Inria Saclay (Defi
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
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 informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationTwo-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations
Two-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations Marcella Bonazzoli 2, Victorita Dolean 1,4, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier
More informationA Parallel Schwarz Method for Multiple Scattering Problems
A Parallel Schwarz Method for Multiple Scattering Problems Daisuke Koyama The University of Electro-Communications, Chofu, Japan, koyama@imuecacjp 1 Introduction Multiple scattering of waves is one of
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 informationToward black-box adaptive domain decomposition methods
Toward black-box adaptive domain decomposition methods Frédéric Nataf Laboratory J.L. Lions (LJLL), CNRS, Alpines Inria and Univ. Paris VI joint work with Victorita Dolean (Univ. Nice Sophia-Antipolis)
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
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 informationASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM C. Canuto 1, L. F. Pavarino 2, and A. B. Pieri 3 1 Introduction Discontinuous Galerkin (DG) methods for partial differential
More information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationA Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems
Contemporary Mathematics Volume 8, 998 B 0-88-0988--03030- A Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems Daoqi Yang. Introduction
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationSolutions of Selected Problems
1 Solutions of Selected Problems October 16, 2015 Chapter I 1.9 Consider the potential equation in the disk := {(x, y) R 2 ; x 2 +y 2 < 1}, with the boundary condition u(x) = g(x) r for x on the derivative
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationAdditive Schwarz method for scattering problems using the PML method at interfaces
Additive Schwarz method for scattering problems using the PML method at interfaces Achim Schädle 1 and Lin Zschiedrich 2 1 Zuse Institute, Takustr. 7, 14195 Berlin, Germany schaedle@zib.de 2 Zuse Institute,
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationA simple FEM solver and its data parallelism
A simple FEM solver and its data parallelism Gundolf Haase Institute for Mathematics and Scientific Computing University of Graz, Austria Chile, Jan. 2015 Partial differential equation Considered Problem
More informationFirst-order overdetermined systems. for elliptic problems. John Strain Mathematics Department UC Berkeley July 2012
First-order overdetermined systems for elliptic problems John Strain Mathematics Department UC Berkeley July 2012 1 OVERVIEW Convert elliptic problems to first-order overdetermined form Control error via
More informationScalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems
Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems Pierre Jolivet, F. Hecht, F. Nataf, C. Prud homme Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann INRIA Rocquencourt
More informationADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013
ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
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 informationHeriot-Watt University
Heriot-Watt University Heriot-Watt University Research Gateway The optimised Schwarz method and the two-lagrange multiplier method for heterogeneous problems in general domains with two general subdomains
More informationDeflated Krylov Iterations in Domain Decomposition Methods
305 Deflated Krylov Iterations in Domain Decomposition Methods Y.L.Gurieva 1, V.P.Ilin 1,2, and D.V.Perevozkin 1 1 Introduction The goal of this research is an investigation of some advanced versions of
More informationPALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations
2013 SIAM Conference On Computational Science and Engineering Boston, 27 th February 2013 PALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations U. Villa,
More information18. Balancing Neumann-Neumann for (In)Compressible Linear Elasticity and (Generalized) Stokes Parallel Implementation
Fourteenth nternational Conference on Domain Decomposition Methods Editors: smael Herrera, David E Keyes, Olof B Widlund, Robert Yates c 23 DDMorg 18 Balancing Neumann-Neumann for (n)compressible Linear
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationOptimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms
Optimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms Marcus Sarkis Worcester Polytechnic Inst., Mass. and IMPA, Rio de Janeiro and Daniel
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 informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationAn Introduction to Domain Decomposition Methods: algorithms, theory and parallel implementation
An Introduction to Domain Decomposition Methods: algorithms, theory and parallel implementation Victorita Dolean, Pierre Jolivet, Frédéric ataf To cite this version: Victorita Dolean, Pierre Jolivet, Frédéric
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationParallel Galerkin Domain Decomposition Procedures for Parabolic Equation on General Domain
Parallel Galerkin Domain Decomposition Procedures for Parabolic Equation on General Domain Keying Ma, 1 Tongjun Sun, 1 Danping Yang 1 School of Mathematics, Shandong University, Jinan 50100, People s Republic
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationOVERLAPPING SCHWARZ ALGORITHMS FOR ALMOST INCOMPRESSIBLE LINEAR ELASTICITY TR
OVERLAPPING SCHWARZ ALGORITHMS FOR ALMOST INCOMPRESSIBLE LINEAR ELASTICITY MINGCHAO CAI, LUCA F. PAVARINO, AND OLOF B. WIDLUND TR2014-969 Abstract. Low order finite element discretizations of the linear
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationBoundary layers for Navier-Stokes equations with slip boundary conditions
Boundary layers for Navier-Stokes equations with slip boundary conditions Matthew Paddick Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie (Paris 6) Séminaire EDP, Université Paris-Est
More informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
More informationSome examples in domain decomposition
Some examples in domain decomposition Frédéric Nataf nataf@ann.jussieu.fr, www.ann.jussieu.fr/ nataf Laboratoire J.L. Lions, CNRS UMR7598. Université Pierre et Marie Curie, France Joint work with V. Dolean
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationSpace-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries
Space-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries Laurence Halpern, Caroline Japhet, and Jérémie Szeftel 3 LAGA, Université Paris XIII, Villetaneuse
More informationThe All-floating BETI Method: Numerical Results
The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationIterative Methods for Stokes/Darcy Coupling
Iterative Methods or Stokes/Darcy Coupling Marco Discacciati Ecole Polytechnique Fédérale de Lausanne Institut d Analyse et Calcul Scientiique Chair o Modelling and Scientiic Computing marco.discacciati@epl.ch
More informationUne méthode parallèle hybride à deux niveaux interfacée dans un logiciel d éléments finis
Une méthode parallèle hybride à deux niveaux interfacée dans un logiciel d éléments finis Frédéric Nataf Laboratory J.L. Lions (LJLL), CNRS, Alpines Inria and Univ. Paris VI Victorita Dolean (Univ. Nice
More informationNumerical Simulations of the Fluid-Structure Coupling in Physiological Vessels Mini-Workshop I
Numerical Simulations of the Fluid-Structure Coupling in Physiological Vessels INRIA Rocquencourt, projet REO miguel.fernandez@inria.fr CRM Spring School, Montréal May 16, 2005 Motivations It is generally
More informationAN OPTIMAL BLOCK ITERATIVE METHOD AND PRECONDITIONER FOR BANDED MATRICES WITH APPLICATIONS TO PDES ON IRREGULAR DOMAINS
AN OPTIMAL BLOCK ITERATIVE METHOD AND PRECONDITIONER FOR BANDED MATRICES WITH APPLICATIONS TO PDES ON IRREGULAR DOMAINS MARTIN J. GANDER, SÉBASTIEN LOISEL, AND DANIEL B. SZYLD Abstract. Classical Schwarz
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationCoupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations
Coupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations P. Bochev, J. Cheung, M. D Elia, A. Frishknecht, K. Kim, M. Parks, M. Perego CM4 summer school, Stanford,
More informationOn domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption
On domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Emerson Roxar,
More informationMéthodes de raffinement espace-temps
Méthodes de raffinement espace-temps Caroline Japhet Université Paris 13 & University of Maryland Laurence Halpern Université Paris 13 Jérémie Szeftel Ecole Normale Supérieure, Paris Pascal Omnes CEA,
More informationFinite Element Methods for Maxwell Equations
CHAPTER 8 Finite Element Methods for Maxwell Equations The Maxwell equations comprise four first-order partial differential equations linking the fundamental electromagnetic quantities, the electric field
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationA Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation
A Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation Y. Boubendir, X. Antoine, C. Geuzaine June 28, 2011 Abstract This paper presents a new non-overlapping domain decomposition
More informationA Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation
A Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation Y. Boubendir, X. Antoine, C. Geuzaine May 8, 2012 Abstract This paper presents a new non-overlapping domain decomposition
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More information1 Fundamentals of laser energy absorption
1 Fundamentals of laser energy absorption 1.1 Classical electromagnetic-theory concepts 1.1.1 Electric and magnetic properties of materials Electric and magnetic fields can exert forces directly on atoms
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationSchur Complement Technique for Advection-Diffusion Equation using Matching Structured Finite Volumes
Advances in Dynamical Systems and Applications ISSN 973-5321, Volume 8, Number 1, pp. 51 62 (213) http://campus.mst.edu/adsa Schur Complement Technique for Advection-Diffusion Equation using Matching Structured
More informationCode: 101MAT4 101MT4B. Today s topics Finite-difference method in 2D. Poisson equation Wave equation
Code: MAT MTB Today s topics Finite-difference method in D Poisson equation Wave equation Finite-difference method for elliptic PDEs in D Recall that u is the short version of u x + u y Dirichlet BVP:
More informationAxioms of Adaptivity (AoA) in Lecture 2 (sufficient for optimal convergence rates)
Axioms of Adaptivity (AoA) in Lecture 2 (sufficient for optimal convergence rates) Carsten Carstensen Humboldt-Universität zu Berlin 2018 International Graduate Summer School on Frontiers of Applied and
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 11, pp. 1-24, 2000. Copyright 2000,. ISSN 1068-9613. ETNA NEUMANN NEUMANN METHODS FOR VECTOR FIELD PROBLEMS ANDREA TOSELLI Abstract. In this paper,
More informationAN INTRODUCTION TO DOMAIN DECOMPOSITION METHODS. Gérard MEURANT CEA
Marrakech Jan 2003 AN INTRODUCTION TO DOMAIN DECOMPOSITION METHODS Gérard MEURANT CEA Introduction Domain decomposition is a divide and conquer technique Natural framework to introduce parallelism in the
More informationA Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term
A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term Peter Sternberg In collaboration with Dmitry Golovaty (Akron) and Raghav Venkatraman (Indiana) Department of Mathematics
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationShifted Laplace and related preconditioning for the Helmholtz equation
Shifted Laplace and related preconditioning for the Helmholtz equation Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Schlumberger Gould Research), Martin Gander (Geneva) Douglas
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More information