Space-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries
|
|
- Amberlynn Willis
- 6 years ago
- Views:
Transcription
1 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 9343, France, halpern@math.univ-paris3.fr LAGA, Université Paris XIII, Villetaneuse 9343, France; CSCAMM, University of Maryland College Park, MD 74 USA, japhet@cscamm.umd.edu, the first two authors are partially supported by french ANR COMMA and GdR MoMaS. 3 Département de mathématiques et applications, Ecole Normale supérieure, 45 rue d Ulm, 753 Paris cedex 5 France, Jeremie.Szeftel@ens.fr Introduction In many fields of applications it is necessary to couple models with very different spatial and time scales and complex geometries. Amongst them are ocean-atmosphere coupling and far field simulations of underground nuclear waste disposal. For such problems with long time computations, a splitting of the time interval into windows is essential. This allows for robust and fast solvers in each time window, with the possibility of nonconforming space-time grids, general geometries, and ultimately adaptive solvers. Optimized Schwarz Waveform Relaxation OSWR methods were introduced and analyzed for linear advection-reaction-diffusion problems with constant coefficients in [, 3, 9]. All these methods rely on an algorithm that computes independently in each subdomain over the whole time interval, exchanging space-time boundary data through optimized transmission operators. They can apply to different space-time discretization in subdomains, possibly nonconforming and need a very small number of iterations to converge. Numerical evidences of the performance of the method with variable smooth coefficients were given in [9]. An extension to discontinuous coefficients was introduced in [ 4], with asymptotically optimized Robin transmission conditions in some particular cases. In [, 6], semi-discretization in time in one dimension was performed using discontinuous Galerkin, see [8, ]. In [7], we extended the analysis to the two dimensional case. We obtained convergence results and error estimates for rectangular or stripsubdomains. For the space discretization, we extended numerically the nonconforming approach in [5] to advection-diffusion problems and optimized order transmission
2 76 L. Halpern et al. conditions, to allow for non-matching grids in time and space on the boundary. The space-time projections between subdomains were computed with an optimal projection algorithm without any additional grid, as in [ 5]. In [7], two dimensional simulations with continuous coefficients were presented. We present here new results in two directions: we extend the proof of convergence of the OSWR algorithm to nonoverlapping subdomains with curved interfaces. We also present simulations for two subdomains, with piecewise smooth coefficients and a curved interface, for which no error estimates are available. We finally present an application to the porous media equation. We consider the advection-diffusion-reaction equation, t u + bu ν u+cu = f in R N,T, with initial condition u, and N =. The advection, diffusion and reaction coefficients b, ν and c, are piecewise smooth, we suppose ν ν > a.e.. The Continuous OSWR Algorithm We consider a decomposition into nonoverlapping subdomains Ω i,i {,..., I}, organized as depicted in Fig.. The interfaces between the subdomains are supposed to be flat at infinity. For any i {,..., I}, Ω i is the boundary of Ω i, n i the unit exterior normal vector to Ω i, N i is the set of indices of the neighbors of Ω i.for j N i, Γ i,j is the common interface. Ω i Ω i Fig.. Decomposition in subdomains. Left: Robin transmission conditions, right: second order transmission conditions. Following [,, 3, 4], we introduce the boundary operators S i,j acting on functions defined on Γ i,j : S i,j ϕ = p i,j ϕ + q i,j t ϕ + Γi,j r i,j ϕ s i,j Γi,j ϕ, with respectively Γ and Γ the gradient and divergence operators on Γ. p i,j,q i,j, r i,j, s i,j are real parameters. q i,j =, will be referred to as a Robin operator. We introduce the coupled problems t u i + b i u i ν i u i +c i u i = f in Ω i,t νi ni b i n i ui + S i,j u i = νj ni b j n i uj + S i,j u j on Γ i,j,t,j N i.
3 Space-Time Nonconforming OSWR for Heterogeneous Problems 77 As coefficients ν and b are possibly discontinuous on the interface, we note, for s Γ i,j, ν i s = lim ε νs εn i. The same notation holds for b. Under regularity assumptions, solving is equivalent to solving for i {,..., I} with u i the restriction of u to Ω i. We now introduce an algorithm to solve. An initial guess g i,j is given in L,T Γ i,j for i {,..., I},j N i. We solve iteratively t u k i + b i u k i ν i u k i +c iu k i = f in Ω i,t, νi ni b i n i u k i + S i,j u k i = νj ni b j n i u k j + S i,j u k j on Γ i,j,t,j N i. 3 with the convention ν i ni b i n i u i + S i,j u i = g i,j, j N i. Theorem. Assume b i W, Ω i N, ν i W, Ω i, p i,j W, Γ i,j with p i,j > a.e.. Ifq i,j =,orifq i,j = q > with r i,j W, Γ i,j N, r i,j = r j,i on Γ i,j, s i,j W, Γ i,j, s i,j >, s i,j = s j,i on Γ i,j, the algorithm 3 converges in each subdomain to the solution of problem. Proof. We first need some results in differential geometry. For every j N i, the normal vector n i can be extended in a neighbourhood of Γ i,j as a smooth function ñn i with length one. Let ψ i,j C Ω i, such that ψ i,j in a neighbourhood of Γ i,j, ψ i,j in a neighbourhood of Γ i,k for k N i,k j and j N i ψ i,j > on Ω i. Let ñn i be defined on a neighbourhood of the support of ψ i,j. We can extend the tangential gradient and divergence operators in the support of ψ i,j by: Γi,j ϕ := ϕ ñn i ϕñn i, Γi,j ϕ := ϕ ϕ ñn i ñn i. It is easy to see that Γi,j ϕ Γi,j = Γi,j ϕ, Γi,j ϕ Γi,j = Γi,j ϕ and for ϕ and χ with support in suppψ i,j,wehave Γi,j ϕ χdx= Ω i ϕ Γi,j χdx. Ω i 4 Now we prove Theorem. The key point is to obtain energy estimates for the homogeneous problem, i.e. for f = u =. We sketch the proof in the most difficult case q i,j = q>. For the geometry, we consider the case depicted in the right part of Fig.. In that case Ω i has at most two neighbours with We set ϕ i = ϕ L Ω i, ϕ i = ν i ϕ L Ω, ϕ i i, = ϕ L Ω i, ϕ i,, = ϕ W, Ω i and β i = pi,j+p j N i ψ i,j β i,j with β i,j = j,i.. We multiply the first equation of 3byβ i uk i, integrate on Ω i,t then integrate by parts in space,
4 78 L. Halpern et al. t β iu k i t i + β i u k i τ, i dτ β i b i β i u k i dx dτ Ω i + c i + Ω i b i βi u k i dx dτ ν i β i u k i dx dτ Ω i βi,j ν i ni u k i b i n i u k i Γ i,j uk i dσ dτ =. 5. We multiply the first equation of 3by t u k i, integrate on Ω i,t and integrate by parts in space, t u k i i dτ + uk i t i + c i u k i + b i u k i t u k i dx dτ Ω i ν i ni u k i t u k i dσ dτ =. 6 Γ i,j 3. We multiply the first equation of 3by Γi,j ψi,j r i,j u k i integrate on Ω i,t integrate by parts in space to obtain t u k i Γi,j ψi,jr i,j u k i dx dτ+ b i u k i Γi,j ψi,jr i,j u k i dx dτ Ω i Ω i t + c i u k i Γi,j ψi,jr i,j u k i dx dτ ν i ni u k i Γi,j r i,j u k i dσ dτ Ω i Γ i,j t 4 Γi,j u k i i dτ C ν i u k i i + β iu k i i dτ We multiply the first equation of 3 by Γi,j ψi,j s i,j Γi,j u k i integrate on Ω i,t, integrate by parts in space using 4. Using that ν i u k i Γi,j ψi,j s i,j Γi,j u k i dx dτ Ω i t Γi,j,u k i i dτ C ν i u k i i dτ, we obtain ψ i,j si,j Γi,j u k i t i + Γi,j u k i i dτ t + ψi,j s i,j c i Γi,j u k i dx dτ + ν i ni u k i Γ i,j s i,j Γi,j u k i dσ dτ Ω i Γ i,j t b i u k i Γi,j ψi,j s i,j Γi,j u k i dx dτ + C ν i u k i i dτ. 8 Ω i We add 6, 7 and 8, multiply the result by q, and add it to 5. We use ab a ε + ε b in the integral terms in the right-hand side, simplify with the left-hand side, and obtain
5 Space-Time Nonconforming OSWR for Heterogeneous Problems 79 β i u k i t i + q uk i t i +q ψ i,j si,j Γi,j u k i t i + β i u k i τ, i dτ + q t u k i i dτ + q 8 Γi,j u k i i dτ t q ν i ni u k i t u k i + Γ i,j r i,j u k i Γ i,j s i,j Γi,j u k i dσ dτ Γ i,j t βi,j ν i ni u k i b i n i u k i Γ i,j uk i dσ dτ q b i i,, + c i i, u k i t i + C β i u k i i dτ + q ν i u k i i dτ. 9 Recalling that s i,j = s j,i on Γ i,j and r i,j = r j,i on Γ i,j, we use now the identity: νi ni u k i b i n i u k i + S i,j u k i νi ni u k i b i n i u k i S j,i u k i =4 βi,j ν i ni u k i b i n i u k i uk i + qν i ni u k i tu k i + Γ i,j r i,j u k i 4 Γi,j s i,j Γi,j u k i +qp i,j p j,i b i n i t u k i + Γi,j r i,j u k i Γi,j s i,j Γi,j u k i u k i +p i,j + p j,i p i,j p j,i b i n i u k i. Replacing into 9, we obtain β i u k i t i + q u k i t i +q ψ i,j si,j Γi,j u k i t i + β i u k i τ, i dτ + q t u k i i dτ + νi ni u k i 4 b i n i u k i S j,iu k i dσ dτ Γ i,j + q 8 Γi,j u k i i dτ νi ni u k i b i n i u k i + S i,j u k i dσ dτ 4 Γ i,j t + p i,j +p j,i p i,j +p j,i +b i n i u k i dσ dτ+ q Γ i,j b i i,, + c i i, u k i t i + q p i,j +p j,i +b i n i t u k i + Γi,j r i,j u k i Γi,j s i,j Γi,j u k i u k i dσ dτ Γ i,j + C β i u k i i dτ + q ν i u k i i dτ. In order to estimate the fourth term in the right-hand side of, we observe that p i,j +p j,i +b i n i u k i tu k i dσ dτ = p i,j +p j,i +b i n i u k i Γ i,j t dσ. Γ i,j By the trace theorem in the right-hand side, we write: p i,j + p j,i +b i n i u k i t u k i dσ dτ C u k i t i ν i u k i t i, Γ i,j
6 8 L. Halpern et al. and t u k i t i = t u k i u k i t u k i i Ω i u k i i, we obtain q t p i,j + p j,i +b i n i u k i t u k i dσ dτ Γ i,j q t u k i i dτ + q 8 4 uk i t i +C β i u k i i dτ. 3 Moreover, integrating by parts and using the trace theorem, we have: q Γi,j s i,j Γi,j u k i i,j + p j,i +b i n i u k i dσ dτ Γ i,j q 6 Γi,j u k i i dτ + C Γi,j u k i i dτ + β i u k i i dτ. 4 Using, we estimate the third term in the right-hand side of by: q b i i,, + c i i, u k i t i q t u k i i dτ + C β i u k i i dτ. 5 8 Replacing 4, 3 and 5 in, then using the transmission conditions, we have: β i u k i t i + q uk i t i +q ψ i,j si,j Γi,j u k i t i + β i u k i τ, i dτ + q t u k i 4 i dτ + q 6 Γi,j u k i i dτ + νi ni u k i b i n i u k i S j,i u k i dσ dτ 4 Γ i,j νj ni u k j b j n i u k j + S i,j u k j dσ dτ 4 Γ i,j + C β i u k i i dτ + q ν i u k i i dτ. We now sum up over the interfaces j N i, then over the subdomains for i I, and over the iterations for k K, the boundary terms cancel out, and with αt = 4 i {,...,I} j N νj i Γ i,j ni u i b j n i u j + S i,juj dσ dτ, we obtain for any t,t, β i u k i t i + q ν i u k i t i + ν β i u k i i dτ k {,...,K} i {,...,I} αt+c β i u k i i dτ + q ν i u k i i dτ. k {,...,K} i {,...,I}
7 Space-Time Nonconforming OSWR for Heterogeneous Problems 8 We conclude with Gronwall s lemma that the sequence converges in L,T; H Ω i. 3 Numerical Results We recall the discrete time nonconforming Schwarz waveform relaxation method developed in [7]. Let T i be the time partition in subdomain Ω i, with N i +intervals In, i and time step kn i. We define interpolation operators I i and projection operators P i in each subdomain as in [7], and we solve t I i Ui k+ bu i k ν i Ui k+c i Ui k = P i f in Ω i,t, νi ni b i n i U k i + S i,j Ui k = P i ν j ni b j n i U k j + S i,j Uj k on Γi,j,T, with S i,j U = p i,j U + q i,j t I i U+ Γi,j r i,j U s i,j Γi,j U, and S i,j U = p i,j U + q i,j t I j U+ Γi,j r i,j U s i,j Γi,j U. The coefficients p i,j and q i,j are defined through an optimization procedure, see [], restricted to values such that the subdomain problems are well-posed. The time semi-discrete analysis was performed in [7] in the case b =. For the space discretization, we use the nonconforming approach in [ 5] extended to problem and order transmission conditions, to allow non-matching grids in time and space on the boundary. We have implemented the algorithm with P finite elements in space in each subdomain. Time windows are used in order to reduce the number of iterations of the algorithm. To reduce the number of parameters and following [ ], we choose r i,j = Π Γi,j b j with Π Γi,j the tangential trace on Γ i,j, and s i,j = ν j even though the present analysis does not cover this case. The optimized parameters are constant along the interface. They correspond to a mean value of the parameters obtained by a numerical optimization of the convergence factor. We first give an example of a multidomain solution with one time window. The physical domain is Ω =,,, the final time is T =. The initial value is u =.5e x.55 +y.7 and the right-hand side is f =. The domain Ω is split into two subdomains Ω =,.5, and Ω =.5,,. The reaction factor c is zero, the advection and diffusion coefficients are bb =,, ν =. y, and b =.,, ν =.sinxy. The mesh size over the interface and time step in Ω are h =/3 and k =/8, while in Ω, h = /4 and k = /94. On Fig., we observe, at final time T =, a very good behavior of the multidomain solution after 5 iterations. The relative error with the one domain solution is of the same order as the error of the scheme. We analyze now the precision in time. The space mesh is conforming and the converged solution is such that the residual is smaller than 8. We compute a variational reference solution on a time grid with 4,96 time steps. The nonconforming solutions are interpolated on the previous grid to compute the error. We start with a time grid with 8 time steps for the left domain and 94 time steps for the right domain. Thereafter the time steps are divided by several times. Figure 3 left shows
8 8 L. Halpern et al. y x Fig.. Nonconforming DG-OSWR solution after 5 iterations. conforming finner grid conforming coarser grid nonconforming domain nonconforming domain k Order Robin Error 3 4 Error Number of refinements 5 5 Iterations Fig. 3. Error between variational and DG-OSWR solutions versus the refinement in time left, and versus the iterations right. the norms of the error in L I; L Ω i versus the number of refinements, for both subdomains. First we observe the order in time for the nonconforming case. This fits the theoretical estimates, even though we have theoretical results only for Robin transmission conditions. Moreover, the error obtained in the nonconforming case, in the subdomain where the grid is finer, is nearly the same as the error obtained in the conforming finer case.
9 Space-Time Nonconforming OSWR for Heterogeneous Problems 83 The computations are done using Order transmissions. Indeed, the error between the multidomain and variational solutions decreases much faster with the Order transmissions conditions than with the Robin transmissions conditions as we can see on Fig. 3 right, in the conforming case. We now consider advection-diffusion equations with discontinuous porosity: ω t u + bu ν u =. The physical domain is Ω =,,, the final time is T =.5. The initial value is u =.5e x.7 +y.5. Domain Ω is split into two subdomains Ω,.5 and Ω,.5 with sinπs 8, s, <s<a parametrization of the interface, as in Fig. 4. The advection and diffusion coeffi y x Fig. 4. Domain decomposition with Ω left and Ω right. cients are b = sin π y cosπx, 3cos π y sinπx, ν =.3, ω =., and b = b, ν =., ω =. We consider first a conforming grid in space. The mesh size over the interface is h =/4 and time step in Ω is k =/8, while in Ω, k =/94. On Fig. 5, we observe, at final time T =.5, that the approximate solution computed using ten time windows and 3 iterations in each time window is close to the variational solution computed in one time window on the conforming finer space-time grid as shown on the error. We now consider nonconforming grids in space as shown on Fig. 4. The mesh size over the interface and time step in Ω are h =/4 and k =/8, while in Ω, h =/8 and k =/94. On Fig. 6, we observe, at final time T =.5, that the approximate solution computed using 5 iterations in one time window is close to the variational solution computed on the conforming finer space-time grid. On Fig. 7 we observe the precision versus the mesh size and time step. The converged solution is such that the residual is smaller than 8. A variational reference solution is
10 84 L. Halpern et al. computed on a time grid with,48 time steps and 384 mesh grid. The space-time nonconforming solutions are interpolated on the previous grid to compute the error. We start with a time grid with 3 time steps and 4 mesh size for the left domain and time steps and mesh size for the right domain and divide by the time step and mesh size several times. Figure 7 shows the norms of the error in L I; L Ωi versus the time steps, for both subdomains. We observe the order for the nonconforming space-time case, even though we have theoretical results only for the time semi-discretized case in [7]. 4 x y x Fig. 5. Error between variational and DG-OSWR solutions, at final time, after time windows and 3 iterations per window. 4 Conclusions We have analyzed the continuous algorithm for variable discontinuous coefficients and general decompositions. We have shown numerically that the method preserves the order of the one domain scheme in the case of discontinuous variable coefficients, nonconforming grids in space and time and a curved interface. An analysis of the influence of the decomposition in time windows is in progress.
11 Space-Time Nonconforming OSWR for Heterogeneous Problems x Fig. 6. DG-OSWR solution at final time, after 5 iterations. Error versus the time step Domain Domain Slope Error y Time step Fig. 7. Error curves versus the time step.
12 86 L. Halpern et al. References. D. Bennequin, M.J. Gander, and L. Halpern. A homographic best approximation problem with application to optimized Schwarz waveform relaxation. Math. Comput., 78:85 3, 9.. E. Blayo, L. Halpern, and C. Japhet. Optimized Schwarz waveform relaxation algorithms with nonconforming time discretization for coupling convection-diffusion problems with discontinuous coefficients. In O.B. Widlund and D.E. Keyes, editors, Decomposition Methods in Science and Engineering XVI, volume 55 of Lecture Notes in Computational Science and Engineering, pp Springer Berlin, Heidelberg, New York, M.J. Gander and L. Halpern. Optimized Schwarz waveform relaxation methods for advection reaction diffusion problems. SIAM J. Numer. Anal., 45: , M.J. Gander, L. Halpern, and M. Kern. Schwarz waveform relaxation method for advection diffusion reaction problems with discontinuous coefficients and non-matching grids. In O.B. Widlund and D.E. Keyes, editors, Decomposition Methods in Science and Engineering XVI, volume 55 of Lecture Notes in Computational Science and Engineering, pp Springer Berlin, Heidelberg, New York, M.J. Gander, C. Japhet, Y. Maday, and F. Nataf. A new cement to Glue nonconforming grids with Robin interface conditions : The finite element case. In R. Kornhuber, R.H.W. Hoppe, J. Périaux, O. Pironneau, O.B. Widlund, and J. Xu, editors, Domain Decomposition Methods in Science and Engineering, volume 4 of Lecture Notes in Computational Science and Engineering, pp Springer Berlin, Heidelberg, New York, L. Halpern and C. Japhet. Discontinuous Galerkin and nonconforming in time optimized Schwarz waveform relaxation for heterogeneous problems. In U. Langer, M. Discacciati, D.E. Keyes, O.B. Widlund, and W. Zulehner, editors, Decomposition Methods in Science and Engineering XVII, volume 6 of Lecture Notes in Computational Science and Engineering, pp. 9. Springer Berlin, Heidelberg, New York, L. Halpern, C. Japhet, and J. Szeftel. Discontinuous Galerkin and nonconforming in time optimized Schwarz waveform relaxation. In Proceedings of the Eighteenth International Conference on Domain Decomposition Methods, 9. in electronic form. These proceedings in printed form. 8. C. Johnson, K. Eriksson, and V. Thomée. Time discretization of parabolic problems by the discontinuous Galerkin method. RAIRO Modél. Math. Anal. Numér., 9, V. Martin. An optimized Schwarz waveform relaxation method for the unsteady convection diffusion equation in two dimensions. Appl. Numer. Math., 5:4 48, 5.. V. Thomée. Galerkin Finite Element Methods for Parabolic Problems. Springer Berlin, Heidelberg, New York, 997.
Mé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 informationNONCONFORMING IN TIME DOMAIN DECOMPOSITION METHOD FOR POROUS MEDIA APPLICATIONS
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira Eds Lisbon, Portugal,14-17 June 2010 NONCONFORMING IN TIME DOMAIN DECOMPOSITION METHOD FOR POROUS
More informationSpace-time domain decomposition for a nonlinear parabolic equation with discontinuous capillary pressure
Space-time domain decomposition for a nonlinear parabolic equation with discontinuous capillary pressure Elyes Ahmed, Caroline Japhet, Michel Kern INRIA Paris Maison de la Simulation ENPC University Paris
More informationOptimized Schwarz Waveform Relaxation: Roots, Blossoms and Fruits
Optimized Schwarz Waveform Relaxation: Roots, Blossoms and Fruits Laurence Halpern 1 LAGA, Université Paris XIII, 99 Avenue J-B Clément, 93430 Villetaneuse, France, halpern@math.univ-paris13.fr 1 Introduction:
More informationA posteriori error estimates for space-time domain decomposition method for two-phase flow problem
A posteriori error estimates for space-time domain decomposition method for two-phase flow problem Sarah Ali Hassan, Elyes Ahmed, Caroline Japhet, Michel Kern, Martin Vohralík INRIA Paris & ENPC (project-team
More informationAdvection Diffusion Problems with Pure Advection Approximation in Subregions
Advection Diffusion Problems with Pure Advection Approximation in Subregions M. J. Gander, L. Halpern 2, C. Japhet 2, and V. Martin 3 Université de Genève, 2-4 rue du Lièvre, CP 64, CH-2 Genève. Switzerland.
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 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 informationTIME SPACE DOMAIN DECOMPOSITION AND REACTIVE TRANSPORT IN POROUS MEDIA
XVIII International Conference on Water Resources CMWR 2 J. Carrera (Ed) c CIMNE, Barcelona, 2 TIME SPACE DOMAIN DECOMPOSITION AND REACTIVE TRANSPORT IN POROUS MEDIA Florian Haeberlein and Anthony Michel
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 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 informationSpace-Time Domain Decomposition Methods for Transport Problems in Porous Media
1 / 49 Space-Time Domain Decomposition Methods for Transport Problems in Porous Media Thi-Thao-Phuong Hoang, Elyes Ahmed, Jérôme Jaffré, Caroline Japhet, Michel Kern, Jean Roberts INRIA Paris-Rocquencourt
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 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 informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
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 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 informationTWO-SCALE SIMULATION OF MAXWELL S EQUATIONS
ESAIM: PROCEEDINGS, February 27, Vol.6, 2-223 Eric Cancès & Jean-Frédéric Gerbeau, Editors DOI:.5/proc:27 TWO-SCALE SIMULATION OF MAWELL S EQUATIONS Hyam Abboud,Sébastien Jund 2,Stéphanie Salmon 2, Eric
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 informationFinite volume method for nonlinear transmission problems
Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings
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 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 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 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 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 informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
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 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 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 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 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 informationOptimized Schwarz methods in spherical geometry with an overset grid system
Optimized Schwarz methods in spherical geometry with an overset grid system J. Côté 1, M. J. Gander 2, L. Laayouni 3, and A. Qaddouri 4 1 Recherche en prévision numérique, Meteorological Service of Canada,
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationDETERMINATION OF THE BLOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION
DETERMINATION OF THE LOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION y FRANK MERLE and HATEM ZAAG Abstract. In this paper, we find the optimal blow-up rate for the semilinear wave equation with a power nonlinearity.
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
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 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 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 informationThe mortar element method for quasilinear elliptic boundary value problems
The mortar element method for quasilinear elliptic boundary value problems Leszek Marcinkowski 1 Abstract We consider a discretization of quasilinear elliptic boundary value problems by the mortar version
More informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
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 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 informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationInternal Stabilizability of Some Diffusive Models
Journal of Mathematical Analysis and Applications 265, 91 12 (22) doi:1.16/jmaa.21.7694, available online at http://www.idealibrary.com on Internal Stabilizability of Some Diffusive Models Bedr Eddine
More informationTwo new enriched multiscale coarse spaces for the Additive Average Schwarz method
346 Two new enriched multiscale coarse spaces for the Additive Average Schwarz method Leszek Marcinkowski 1 and Talal Rahman 2 1 Introduction We propose additive Schwarz methods with spectrally enriched
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow
ANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow Michel Kern Inria Paris Rocquencourt Maison de la Simulation C2S@Exa Days, Inria Paris Centre, Novembre 2016 M. Kern
More informationarxiv: v1 [math.na] 6 Nov 2017
Efficient boundary corrected Strang splitting Lukas Einkemmer Martina Moccaldi Alexander Ostermann arxiv:1711.02193v1 [math.na] 6 Nov 2017 Version of November 6, 2017 Abstract Strang splitting is a well
More informationFEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
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 informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationComparison 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 informationarxiv: v1 [math.na] 16 Dec 2015
ANALYSIS OF A NEW HARMONICALLY ENRICHED MULTISCALE COARSE SPACE FOR DOMAIN DECOMPOSITION METHODS MARTIN J. GANDER, ATLE LONELAND, AND TALAL RAHMAN arxiv:52.05285v [math.na] 6 Dec 205 Abstract. We propose
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationEfficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems
Efficient smoothers for all-at-once multigrid methods for Poisson and Stoes control problems Stefan Taacs stefan.taacs@numa.uni-linz.ac.at, WWW home page: http://www.numa.uni-linz.ac.at/~stefant/j3362/
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 informationDOMAIN DECOMPOSITION APPROACHES FOR MESH GENERATION VIA THE EQUIDISTRIBUTION PRINCIPLE
DOMAIN DECOMPOSITION APPROACHES FOR MESH GENERATION VIA THE EQUIDISTRIBUTION PRINCIPLE MARTIN J. GANDER AND RONALD D. HAYNES Abstract. Moving mesh methods based on the equidistribution principle are powerful
More informationA non-standard Finite Element Method based on boundary integral operators
A non-standard Finite Element Method based on boundary integral operators Clemens Hofreither Ulrich Langer Clemens Pechstein June 30, 2010 supported by Outline 1 Method description Motivation Variational
More informationEquilibrium analysis for a mass-conserving model in presence of cavitation
Equilibrium analysis for a mass-conserving model in presence of cavitation Ionel S. Ciuperca CNRS-UMR 58 Université Lyon, MAPLY,, 696 Villeurbanne Cedex, France. e-mail: ciuperca@maply.univ-lyon.fr Mohammed
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 informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationarxiv: v1 [math.na] 14 May 2016
NONLINEAR PRECONDITIONING: HOW TO USE A NONLINEAR SCHWARZ METHOD TO PRECONDITION NEWTON S METHOD V. DOLEAN, M.J. GANDER, W. KHERIJI, F. KWOK, AND R. MASSON arxiv:65.449v [math.na] 4 May 6 Abstract. For
More informationPostprint.
http://www.diva-portal.org Postprint This is the accepted version of a chapter published in Domain Decomposition Methods in Science and Engineering XXI. Citation for the original published chapter: Gander,
More informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
More informationAn additive average Schwarz method for the plate bending problem
J. Numer. Math., Vol. 10, No. 2, pp. 109 125 (2002) c VSP 2002 Prepared using jnm.sty [Version: 02.02.2002 v1.2] An additive average Schwarz method for the plate bending problem X. Feng and T. Rahman Abstract
More informationOptimal Schwarz Waveform Relaxation for the One Dimensional Wave Equation. GANDER, Martin Jakob, HALPERN, Laurence, NATAF, F.
Article Optimal Schwarz Waveform Relaxation for the One Dimensional Wave Equation GANDER, Martin Jakob, HALPERN, Laurence, NATAF, F. Reference GANDER, Martin Jakob, HALPERN, Laurence, NATAF, F. Optimal
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 informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationOn the rate of convergence of Schwarz waveform relaxation methods for the time-dependent Schrödinger equation
On the rate of convergence of Schwarz waveform relaxation methods for the time-dependent Schrödinger equation X. Antoine a,b, E. Lorin d,c a Institut Elie Cartan de Lorraine, Université de Lorraine, F-54506
More informationAdditive Average Schwarz Method for a Crouzeix-Raviart Finite Volume Element Discretization of Elliptic Problems
Additive Average Schwarz Method for a Crouzeix-Raviart Finite Volume Element Discretization of Elliptic Problems Atle Loneland 1, Leszek Marcinkowski 2, and Talal Rahman 3 1 Introduction In this paper
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
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 informationA FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem
A FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem Leszek Marcinkowski 1 and Nina Dokeva 2 1 Department of Mathematics, Warsaw University, Banacha 2, 02 097 Warszawa, Poland,
More information/00 $ $.25 per page
Contemporary Mathematics Volume 00, 0000 Domain Decomposition For Linear And Nonlinear Elliptic Problems Via Function Or Space Decomposition UE-CHENG TAI Abstract. In this article, we use a function decomposition
More informationCut Finite Elements for Convection in Fractured Domains
Cut Finite Elements for Convection in Fractured Domains Erik Burman, Peter Hansbo, Mats G. Larson and Karl Larsson arxiv:1801.06103v2 [math.na] 30 May 2018 Abstract We develop a cut finite element method
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 informationParameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition
Parameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition P. A. Farrell 1,A.F.Hegarty 2, J. J. H. Miller 3, E. O Riordan 4,andG.I. Shishkin 5 1
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 informationStochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows
Stochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows Ivan otov Department of Mathematics, University of Pittsburgh
More informationThe Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 29, 1999 Abstract This paper deals with the construction of wavelet approximatio
The Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 9, 1999 Abstract This paper deals with the construction of wavelet approximation spaces, in the framework of the Mortar method. We
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationOn the performance of the Algebraic Optimized Schwarz Methods with applications
On the performance of the Algebraic Optimized Schwarz Methods with applications Report 12-1-18 Department of Mathematics Temple University October 212 This report is available at http://www.math.temple.edu/~szyld
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 informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationTWO-SCALE CONVERGENCE ON PERIODIC SURFACES AND APPLICATIONS
TWO-SCALE CONVERGENCE ON PERIODIC SURFACES AND APPLICATIONS Grégoire ALLAIRE Commissariat à l Energie Atomique DRN/DMT/SERMA, C.E. Saclay 91191 Gif sur Yvette, France Laboratoire d Analyse Numérique, Université
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationSome Geometric and Algebraic Aspects of Domain Decomposition Methods
Some Geometric and Algebraic Aspects of Domain Decomposition Methods D.S.Butyugin 1, Y.L.Gurieva 1, V.P.Ilin 1,2, and D.V.Perevozkin 1 Abstract Some geometric and algebraic aspects of various domain decomposition
More informationOn solving linear systems arising from Shishkin mesh discretizations
On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More informationComparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes
Comparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes Do Y. Kwak, 1 JunS.Lee 1 Department of Mathematics, KAIST, Taejon 305-701, Korea Department of Mathematics,
More information