A MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS
|
|
- Lesley Lawson
- 5 years ago
- Views:
Transcription
1 SIAM J. SCI. COMPUT. Vol. 29, No. 5, pp c 2007 Society for Industrial and Applied Mathematics A MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS ZHIQIANG CAI AND YANQIU WANG Abstract. The purpose of this paper is to develop and analyze a multigrid solver for the finite element discretization of the pseudostress system associated with the differential operator A γgrad div over 2 2 matrix-valued functions. This system is derived from the pseudostressvelocity formulation [Z. Cai, B. Lee, and P. Wang, SIAM J. Numer. Anal., 42 (2004), pp ] of two-dimensional Stokes problems through the penalty method or natural time discretization for the respective steady- or unsteady-state problems. Here γ>0is a constant associated with either the penalty parameter or the time-step size, and A is a singular map. In this paper, we develop a multigrid solver for the discrete problem using both the Raviart Thomas (RT) and the Brezzi Douglas Marini 1+γ (BDM) finite elements. We show that the multigrid convergence rate is O( 2 1+γ 2 ), where m is +m the number of smoothings. This convergence rate is independent of the mesh size and the number of levels used in multigrid. Moreover, numerical results presented in this paper show that the multigrid convergence rate for the BDM elements does not depend on the parameter γ. Key words. multigrid, mixed finite element, pseudostress-velocity formulation, Stokes problem AMS subject classifications. 65N30, 65N55 DOI / Introduction. The stress-velocity-pressure formulation [7, 15] for incompressible Newtonian flows comes from the original physical equations. Recently, this formulation has been gaining consistent attention because of the arising interest in non-newtonian flows [6, 21]. For some non-newtonian flows, the stress cannot be eliminated and, hence, a formulation containing the stress as a fundamental unknown is unavoidable. In addition, accurate direct calculation of the stress is important for computing the traction on a fluid-solid interface. Although the stress can be obtained in the velocity-pressure method by differentiating the velocity, this degrades the order of the approximation. In this paper we are interested in mixed finite element methods for the Stokes problem. Notice that one disadvantage of using the stress-velocity-pressure formulation is that the symmetry constraint of the stress tensor poses extra difficulty in the discretization. For example, stable conforming mixed elements for symmetric matrices, e.g., [3, 5, 20], employ many degrees of freedom and, consequently, lead to large linear systems even for coarse meshes. There are different ways to avoid this difficulty. One way is to use the idea of pseudostress [11] and a pseudostress-velocity formulation. The pseudostress is nonsymmetric and, hence, any stable pairs for the Darcy flows can be adapted to the pseudostress-velocity formulation. This is documented in [13]. Furthermore, the velocity can be eliminated through the penalty method for the steady-state problem or the natural time discretization for the unsteady-state problem. So one need only solve numerically the pseudostress system associated with a differential operator A γ grad div over 2 2 matrix-valued functions. Here, γ>0 Received by the editors May 30, 2006; accepted for publication (in revised form) March 27, 2007; published electronically September 28, This work was supported in part by the National Science Foundation under grant DMS Department of Mathematics, Purdue University, 150 N. University Street, West Lafayette, IN (zcai@math.purdue.edu, yqwang@math.purdue.edu). 2078
2 MULTIGRID FOR STOKES PROBLEMS 2079 is a constant depending on either the penalty parameter or the time-step size, and A is a projection operator whose kernel intersects the kernel of the divergence operator. The purpose of this paper is to construct and analyze a multigrid solver for this pseudostress system. It must be pointed out that this problem is more delicate to handle than a normal H(div) problem because of the complicated kernel structure. We also mention other methods of circumventing the difficulty posed by the symmetric stress. One is to use a stabilization technique [7]. Another is to use the leastsquares method, which does not require stable pairs of finite elements. Least-squares methods with fast multigrid solvers for the stress-velocity-pressure formulation and its variants, stress-velocity and pseudostress-velocity formulations, were studied in [11]. Multigrid methods are efficient for many problems. A vast amount of research has been done in this area, e.g., [8, 18, 22, 25]. However, the classical techniques for the multigrid method do not work for the H(div)-type problem since standard smoothers cannot damp out the highly oscillating part of divergence-free components. Hence, special smoothers have to be developed in order to achieve an optimal multigrid convergence rate [1, 2, 19, 24]. An additional difficulty in the multigrid analysis for the pseudostress system is the handling of the kernel of A. To overcome the above difficulties, we employ a Helmholtz-like decomposition of the discrete space, which also utilizes the operator A, so that the analysis can be done on separate subspaces. In this paper we consider multigrid solvers for systems discretized using both the Raviart Thomas (RT) and the Brezzi Douglas Marini (BDM) finite elements. We show that the multigrid convergence rate is bounded above by ρ = C(1+γ 2 ) C(1+γ 2 )+2m < 1, where m is the number of smoothings and C is a constant independent of the mesh size or the number of levels used in the multigrid. When γ 1, which is always true for the penalty method of the stationary problems as the penalty parameter goes to zero, the multigrid convergence for the pseudostress system can be considered independent of γ and, hence, the penalty parameter. When γ 1, which corresponds to the small time-step size in the time discretization of unsteady-state problems, ρ is close to one. Nevertheless, numerical results presented in this paper show that the multigrid convergence rate for the BDM elements does not depend on the parameter γ even if γ is small. We give the outline of the paper as follows. In sections 2 and 3, we describe the pseudostress system of the Stokes problem and its finite element discretization. The key Helmholtz-type decomposition of the finite element space is given in section 4. We construct a multigrid solver for the discrete H(div)-type problem and give a convergence rate estimate in section 5. Some supporting numerical results are presented in section 6. Finally, in the appendix, we prove an essential inequality needed by the analysis in section Notation. Let Ω be a convex polygon. Denote by R 2 the set of twodimensional vector-valued functions and by M 2 the set of 2 2 matrix-valued functions. Throughout the paper, we adopt the convention that a Greek character denotes a 2 2 matrix and a bold roman character in lower case denotes a vector. Letting τ =(τ ij ) 2 2 M 2 and v =(v 1,v 2 ) t R 2, define div v = v1 x + v2 y, divτ = v = ( v1 x v 2 x v 1 y v 2 y ), curl v = ( τ11 x τ 21 x ( v 1 y v2 y + τ12 y + τ22 y v 1 x v 2 x ), ).
3 2080 ZHIQIANG CAI AND YANQIU WANG Define the inner product between vectors and 2 2 matrices by u v = u 1 v 1 + u 2 v 2 and σ : τ = σ ij τ ij, 1 i,j 2 respectively. We use the usual notation L 2 (Ω) to denote the set of square integrable functions on Ω and H s (Ω), where s>0 is a real number, to denote the Sobolev space on Ω with index s. Let C0 (Ω) be the space of infinitely differentiable functions on Ω with compact support. Denote by H0(Ω) s the closure of C0 (Ω) under the norm of H s (Ω) and define H s (Ω)=(H0(Ω)) s, the dual space of H0(Ω). s It is natural to extend the above definitions to vector- and matrix-valued functions using product spaces. For example, H s (Ω, R 2 ) and H s (Ω, M 2 ) denote the Sobolev spaces over the set of vector- and 2 2 matrix-valued functions, respectively. Other notations, such as L 2 (Ω, R 2 ), L 2 (Ω, M 2 ), are defined in the same fashion. For simplicity, denote s,ω to be the Sobolev norm with index s over scalar-, vector-, or matrix-valued functions, depending on the type of the function. Similarly, denote by s,ω the Sobolev seminorm and by (, ) the L 2 inner product. Define the space H(div, Ω, M 2 )={σ L 2 (Ω, M 2 ) div σ L 2 (Ω, R 2 )} with the norm σ 2 H(div,Ω,M 2) =(σ, σ)+(divσ, divσ). 2. Pseudostress system. Consider the two-dimensional steady-state Stokes problem Δu + p = f in Ω, (2.1) div u = 0 in Ω, u Ω = 0, where u is the velocity and p is the pressure which satisfies the compatibility condition pdx =0. Ω Let A : M 2 M 2 be a fourth order tensor defined by A τ = τ 1 (tr τ ) I, 2 where tr τ = τ 11 + τ 22 is the trace of τ and I is the 2 2 identity matrix. We immediately notice that A is a projection onto the trace-free subspace of M 2 and Ker(A) ={v I v is a scalar function}. Define the pseudostress by σ = p I + u. Then the Stokes problem (2.1) can be rewritten in the pseudostress-velocity formulation: divσ = f in Ω, (2.2) Aσ u = 0 in Ω, u Ω = 0. The incompressible constraint div u = 0 is enforced through div u =tr( u) =0. Notice that tr σ = 2p also needs to satisfy the compatibility condition tr σ dx =0. Ω
4 MULTIGRID FOR STOKES PROBLEMS 2081 From the well-known regularity estimate of problem (2.1) on convex polygons and the definition of the pseudostress, we have the following regularity results. Let f L 2 (Ω); then the solutions to problem (2.1) and (2.2) satisfy u H 2 (Ω) H 1 0(Ω), p H 1 (Ω)/R, σ H 1 (Ω, M 2 ), and (2.3) u 2,Ω c f 0,Ω, p H 1 (Ω)/R c f 0,Ω, σ 1,Ω c f 0,Ω, where c is a positive constant independent of f. Remark 1. Notice that the symmetric physical stress σ phy = pi+( u+( u) t ) and the pressure p can be expressed algebraically in terms of the pseudostress σ: σ phy = σ +(Aσ) t, p = 1 tr σ. 2 Therefore, they can be computed in a postprocessing procedure in the same accuracy as the approximation of σ. Denote V = L 2 (Ω, R 2 ) and { } Σ = H(div, Ω, M 2 )/span{i} = τ H(div, Ω, M 2 ) tr τ dx =0. Ω Then we have the following variational form for (2.2). Given f V, find σ Σ and u V such that (2.4) { (Aσ, τ )+(divτ, u) = 0 for all τ Σ, (divσ, v) =(f, v) for all v V. Notice that the velocity boundary condition becomes the nature boundary condition in the mixed formulation. According to [10], the existence and uniqueness of problem (divτ,v) (2.4) follows from the well-known inf-sup condition sup τ H 1 (Ω,M 2) τ 1,Ω κ v 0,Ω for all v V and the following lemma [12]. Lemma 2.1. For all τ Σ, we have τ 2 0,Ω C A ( A 1/2 τ 2 0,Ω + divτ 2 1,Ω), where C A 1 is independent of τ. Using the penalty method [10] on problem (2.4) gives a new system (2.5) { (Aσγ, τ )+(divτ, u γ ) = 0 for all τ Σ, (divσ γ, v) 1 γ (u γ, v) =(f, v) for all v V, where γ is a positive constant which is usually very large. One advantage of using the penalty method is that the mixed system (2.5) can be easily decoupled. Define a bilinear form over H(div, Ω, M 2 )by Λ(σ, τ )=(Aσ, τ )+γ(divσ, divτ ) for all σ, τ Σ. By eliminating u γ, system (2.5) can be reduced to the following H(div) problem: (2.6) Λ(σ γ, τ )=γ(f, divτ ) for all τ Σ. By Lemma 2.1, it is easy to see that Λ(, ) is coercive over space Σ. Hence problem (2.6) is well-posed. Then u γ can be easily calculated through u γ = γ(divσ γ f).
5 2082 ZHIQIANG CAI AND YANQIU WANG Finally, we mention the case of time-dependent Stokes problems. It is not hard to see that the implicit time discretization usually leads to the following pseudostressvelocity form: divσ 1 γ u = f in Ω, (2.7) Aσ u = 0 in Ω, u Ω = 0, where γ 1 is related to the time-step size. Again, problem (2.7) can easily be decoupled and will lead to the pseudostress problem (2.6). The only difference is that γ is small in this case. 3. Finite element approximation. In the remainder of this paper, we will focus on solving problem (2.6) using the finite element method and a multigrid solver. To discretize the problem, we need a good finite element approximation to the space Σ. There are several well-known conforming elements for the H(div) space over the vector field R 2. We will consider the lowest order RT element [23] and the lowest order BDM element [9]. It should be easy to extend the results to higher order elements. Let T h be a quasi-uniform triangulation of Ω with characteristic mesh size h. On each triangle T T h, denote P n (T ) to be the set of all polynomials with degrees less than or equal to n. Define {( ) ( ) ( )} 1 0 x RT T = span,, 0 1 y and BDM T =(P 1 (T )) 2. The degrees of freedom for the RT element are the zeroth order moments of the normal components on each edge of T. The degrees of freedom for the BDM element are the zeroth order and the first order moments of the normal components on each edge of T. We use two copies of the RT/BDM element to approximate the space H(div, Ω, M 2 ). Let Define the finite element spaces Σ RT T = {τ M 2 (τ i1,τ i2 ) RT T for i =1, 2}, Σ BDM T = {τ M 2 (τ i1,τ i2 ) BDM T for i =1, 2}. Σ RT h = {τ Σ τ T Σ RT T } and Σ BDM h = {τ Σ τ T Σ BDM T }. In later analysis, we use more common properties of the RT and the BDM elements than their differences. Hence we will use Σ h to denote either Σ RT h or Σ BDM h. It is well known that the necessary continuity requirement for Σ h to be in Σ is that any function in the discrete space should have continuous normal components across each edge of the mesh. Notice that Σ h also inherits the constraint tr σ dx = 0 from Σ. Ω Consider the discrete problem of (2.6): find σ h Σ h such that (3.1) Λ(σ h, τ h )=γ (f, div τ h ) for all τ h Σ h.
6 MULTIGRID FOR STOKES PROBLEMS 2083 Let σ and σ h be solutions to problems (2.4) and (3.1), respectively. The following error bound is established in [13] for the RT element with γ = O(h 1 ): σ σ h 0,Ω + div (σ σ h ) 0,Ω h ( σ 1,Ω + f 1,Ω ). Here and thereafter, we use to denote less than or equal to with a factor c independent of the mesh size h or other parameters appearing in the inequality. The above error bound for the BDM element can be obtained in a similar fashion. Next we give the discretization of problem (2.4) and its error estimate, which will play an important role in later analysis. Define the space V h L 2 (Ω, R 2 ) as follows: V h = {(v 1,v 2 ) t v i T P 0 (T ) for i =1, 2}. Since both the RT element and the BDM element satisfy the discrete inf-sup condition [10], we have divσ h = V h. Consider the discrete problem for system (2.4): given f L 2 (Ω, R 2 ), find σ h Σ h and u h V h such that (3.2) { (Aσh, τ h )+(divτ h, u h ) = 0 for all τ h Σ h, (divσ h, v h )=(f, v h ) for all v h V h. By Lemma 2.1 and the discrete inf-sup condition of the RT/BDM element space, Problem (3.2) is well-posed and admits a unique solution. Next, we will give the mixed finite element error estimate for problem (3.2). To do this, we need to show that Σ h, under the constraint tr σ dx =0,isagood Ω approximation for Σ. Denote by Π h the natural interpolation associated with the degrees of freedom onto Σ h + span{i}. Let P Vh be the L 2 projection onto V h. Define Π h : Σ Σ h by (3.3) Π h τ = Π h τ Ω (tr Π h τ ) dx 2 Ω I for all τ Σ, where Ω = Ω dx. Clearly, Π h is a linear operator and Π h τ Σ h. By the properties of the natural interpolation Π h [10] and the definition of Π h, it is easy to prove [13] the following lemma. Lemma 3.1. The commutative property divπ h = P Vh div holds. Furthermore, for τ Σ H 1 (Ω, M 2 ), we have (3.4) (3.5) τ Π h τ 0,Ω h τ 1,Ω, divτ div(π h τ ) 0,Ω divτ 0,Ω. The following lemma gives error estimates for the mixed finite element approximation (3.2). Lemma 3.2. Let f L 2 (Ω, R 2 ). Define (σ, u) to be the solution of problem (2.4) and (σ h, u h ) to be the solution of problem (3.2). We have σ σ h 0,Ω h σ 1,Ω, u u h 0,Ω h u 1,Ω + h σ 1,Ω. Proof. By using the mixed finite element theory and the properties of operator Π h, it is standard to prove that [10, 5] A(σ σ h ) 0,Ω h σ 1,Ω, u u h 0,Ω h u 1,Ω + h σ 1,Ω.
7 2084 ZHIQIANG CAI AND YANQIU WANG Next, notice that div(σ σ h )=f P Vh f. By Lemma 2.1 and the approximation property of the L 2 projection, we have σ σ h 0,Ω A(σ σ h ) 0,Ω + div(σ σ h ) 1,Ω h σ 1,Ω + sup v H 1 0 (Ω,R2 ) (f P Vh f, v) v 1,Ω (f, v P Vh v) = h σ 1,Ω + sup v H 1 0 (Ω,R2 ) v 1,Ω h σ 1,Ω + h f 0,Ω = h σ 1,Ω + h divσ 0,Ω h σ 1,Ω. This completes the proof of the lemma. 4. A decomposition in Σ h. In this section we introduce a Helmholtz-type decomposition of the discrete space Σ h. This decomposition and its properties borrowed many ideas from the decomposition presented in the work of Arnold, Falk, and Winther on H(div ) problems [1, 2]. The difference is that our decomposition needs to use the projection operator A, which makes the proof more complicated. For the convenience of readers, we give all details of the proof even if some parts are similar to the proof given in [1, 2]. Consider the mixed finite element problem (3.2). Let (σ h, u h ) be the solution to (3.2) with the right-hand side function f L 2 (Ω, R 2 ). Denote div 1 h f = σ h. The operator div 1 h maps L 2 (Ω, R 2 )toσ h and can be considered as a pseudoinverse of the div operator. It is not hard to see that for all τ h Σ h such that divτ h = 0, Define (Adiv 1 h f, τ h)= (divτ h, u h )=0. (4.1) (4.2) Σ 0 h = {τ h Σ h such that divτ h = 0}, Σ 1 h = {τ h div 1 h V h}. Then Σ 0 h and Σ 1 h are orthogonal under Λ(, ). It is easy to see that Σ h = Σ 0 h Σ 1 h. Define the space Ũh H 1 (Ω, R 2 ) as follows: for the RT element, let Ũ h = {(q 1,q 2 ) t H 1 (Ω, R 2 ) q 1,q 2 are P 1 polynomials on each T T h } and for the BDM element, let Ũ h = {(q 1,q 2 ) t H 1 (Ω, R 2 ) q 1,q 2 are P 2 polynomials on each T T h }. Define {( ( ( )} 1 0 y U h = 0) Ũh/span,,. 1) x It is well known that every τ h Σ 0 h has a potential q h Ũh such that τ h = curl q h [10]. The constraint Ω tr τ h dx = 0 further implies that q h can be chosen from U h. In other words, Σ 0 h = curl (U h ). Hence we have the following exact sequence: 0 U h curl Σ h div V h 0.
8 MULTIGRID FOR STOKES PROBLEMS 2085 For each τ h Σ h, we have the decomposition σ h = curl q h + div 1 h v h, where q h U h, v h = divσ h V h and the decomposition is orthogonal under both (A, ) and (div, div ). This decomposition is usually referred to as the discrete Helmholtz decomposition. Next, we study some properties related to the discrete Helmholtz decomposition. Define an operator grad h : V h Σ h as follows: (grad h v h, τ h )= (divτ h, v h ) for all τ h Σ h. The operator grad h is actually the discrete dual operator of div. Notice that for an arbitrary σ h Σ h, Aσ h is not necessarily in Σ h since it may not satisfy the H(div, Ω, M 2 ) continuity requirement on normal components. Therefore we need to define discrete operators A h : Σ h Σ h and Λ h : Σ h Σ h by (A h σ h, τ h )=(Aσ h, τ h ) for all τ h Σ h, (Λ h σ h, τ h )=Λ(σ h, τ h ) for all τ h Σ h. It is equivalent to say that A h = P Σh A, where P Σh is the L 2 projection onto the finite dimensional space Σ h and Λ h σ h = A h σ h γgrad h divσ h. By the property of the projection operator P Σh,wehave (4.3) A h τ h 0,Ω Aτ h 0,Ω for all τ h Σ h. Let T H be a quasi-uniform triangulation of Ω with characteristic mesh size H and let T h be a refinement of T H with characteristic mesh size h. We follow the notation previously introduced with the understanding that spaces and operators with subscript H are defined over the mesh T H. For example, Σ H is defined on T H in the same way that Σ h is defined on T h. Besides, the RT/BDM finite element spaces are known to be nested; that is, Σ H Σ h. Lemma 4.1. Let v h V h. Define σ h = div 1 h v h and σ H = div 1 H v h. Then σ h σ H 0,Ω H v h 0,Ω, div(σ h σ H ) 0,Ω H grad h v h 0,Ω. Proof. Let σ Σ be the solution to { (Aσ, τ )+(divτ, u) = 0 for all τ Σ, (divσ, w) =(v h, w) for all w L 2 (Ω, R 2 ). By the regularity property (2.3), Lemma 3.2, and the fact that h<h,wehave σ h σ H 0,Ω σ σ h 0,Ω + σ σ H 0,Ω H σ 1,Ω H v h 0,Ω. Notice that div(σ h σ H )=v h P VH v h, where P VH V H. Define σ Σ by is the L 2 projection onto { (A σ, τ )+(divτ, u) = 0 for all τ Σ, (div σ, w) =(v h P VH v h, w) for all w L 2 (Ω, R 2 ).
9 2086 ZHIQIANG CAI AND YANQIU WANG Recall that Π h and Π H satisfy divπ h = P Vh div and divπ H = P VH div. By using Lemma 3.1, the Schwarz inequality, and the regularity result (2.3), we have v h P VH v h 2 0,Ω =(div σ, v h P VH v h ) =(div(π h σ Π H σ), v h )= ((Π h σ Π H σ), grad h v h ) Π h σ Π H σ 0,Ω grad h v h 0,Ω H σ 1,Ω grad h v h 0,Ω H v h P VH v h 0,Ω grad h v h 0,Ω. This completes the proof of the lemma. Lemma 4.2. Assume σ h Σ h satisfying Λ(σ h, τ H )=0for all τ H Σ H. Let σ h = curl q h + div 1 h v h, where q h U h and v h = divσ h V h, be the Λ-orthogonal decomposition of σ h. Denote σ 1 h = div 1 h v h. Then (4.4) (4.5) (4.6) (Aσ 1 h, σ 1 h) γ 1 H 2 Λ(σ h, σ h ), (1 + γh 2 ) σ 1 h 2 0,Ω (1 + γ 2 )Λ(σ h, σ h ), q h 2 0,Ω H 2 (Acurl q h, curl q h ). Proof. The proof of inequality (4.6) is a little long and hence will be put in the appendix. Here we give the proof of inequalities (4.4) and (4.5). Define σ h Σ h by (4.7) Λ( σ h, τ h )=(Aσ 1 h, τ h ) for all τ h Σ h. It is clear that Λ( σ h, curl p h )=(Aσ 1 h, curl p h ) = 0 for all p h U h. This implies σ h Σ 1 h. In other words, if we denote div σ h = ṽ h, then σ h = div 1 h ṽh. Define σ H = div 1 H ṽh. By the Schwarz inequality, (4.8) (Aσ 1 h,σ 1 h)=λ( σ h, σ 1 h)=λ( σ h, σ h ) = Λ( σ h σ H, σ h ) Λ( σ h σ H, σ h σ H ) 1/2 Λ(σ h, σ h ) 1/2. We also notice that (4.7) implies that Λ h σ h = A h σ 1 h. Furthermore, by setting τ h = σ h in (4.7) and using the Schwarz inequality, we can easily see that Λ( σ h, σ h ) (Aσ 1 h, σ1 h ). Hence, by the fact that A is a projection, Lemma 4.1, the triangular inequality, and inequality (4.3), we have Λ( σ h σ H, σ h σ H ) σ h σ H 2 0,Ω + γ div( σ h σ H ) 2 0,Ω H 2 γ 1 (γ div σ h 2 0,Ω + γ 2 grad h div σ h 2 0,Ω) (4.9) H 2 γ 1 (γ div σ h 2 0,Ω + A h σ h 2 0,Ω + (A h γgrad h div) σ h 2 0,Ω) H 2 γ 1 (Λ( σ h, σ h )+ Λ h σ h 2 0,Ω) H 2 γ 1 ((Aσ 1 h, σ 1 h)+ A h σ 1 h 2 0,Ω) H 2 γ 1 (Aσ 1 h, σ 1 h). Combining (4.8) and (4.9) gives (Aσ 1 h, σ1 h ) Hγ 1/2 (Aσ 1 h, σ1 h )1/2 Λ(σ h, σ h ) 1/2.Thus (4.4) holds.
10 MULTIGRID FOR STOKES PROBLEMS 2087 Next, we prove (4.5). Notice that Λ(σ h, τ H ) = 0 for all τ H Σ H implies (4.10) (divσ 1 h, divτ H )= γ 1 (Aσ 1 h, τ H ) γ 1 (Acurl q h, τ H ). By definition, divσ 1 h 1,Ω = (divσ 1 h sup, v). v H 1 0 (Ω,R2 ) v 1,Ω For v H 1 0(Ω, R 2 ), let v H be the L 2 projection of v onto V H. Then the following estimates are well known: v v H 0,Ω H v 1,Ω, v H 0,Ω v 0,Ω v 1,Ω. Define τ H = div 1 H v H and τ h = div 1 h v H. We have (Acurl q h, τ h ) = 0. Besides, by Lemma 3.2 and the regularity result (2.3), it is not hard to see that τ H 0,Ω v H 0,Ω. Combining the above and using Lemma 4.1 and (4.10), we have (divσ 1 h sup, v) (divσ 1 h sup, v v H) (divσ 1 h + sup, v H) v H 1 0 (Ω,R2 ) v 1,Ω v H 1 0 (Ω,R2 ) v 1,Ω v H 1 0 (Ω,R2 ) v 1,Ω H divσ 1 h 0,Ω + sup v H 1 0 (Ω,R2 ) (divσ 1 h, divτ H) v 1,Ω = H divσ 1 γ 1 (Aσ 1 h h 0,Ω + sup, τ H) γ 1 (Acurl q h, τ H τ h ) v H 1 0 (Ω,R2 ) v 1,Ω H divσ 1 h 0,Ω + γ 1 (Aσ 1 h, σ 1 h) 1/2 + γ 1 H(Acurl q h, curl q h ) 1/2. Using the above result, Lemma 2.1, inequality (4.4), the inequality 2γ 1 1+γ 2, and the fact that H<O(1), we have and σ 1 h 2 0,Ω (Aσ 1 h, σ 1 h)+ divσ 1 h 2 1,Ω (1 + γ 2 )Λ(σ 1 h, σ 1 h)+ divσ 1 h 2 0,Ω + γ 2 (Acurl q h, curl q h ) (1 + γ 2 )Λ(σ h, σ h ) γh 2 σ 1 h 2 0,Ω (1 + γ 2 )γh 2 (Aσ 1 h, σ 1 h)+γ divσ 1 h 2 0,Ω + γ 1 (Acurl q h, curl q h ) (1 + γ 2 )Λ(σ h, σ h ). Adding up the above inequalities gives the estimate (4.5). This completes the proof of the lemma. 5. A multigrid solver. For simplicity, we rewrite problem (3.1) in the following form: find σ h Σ h such that (5.1) Λ(σ h, τ h )=(F, τ h ) for all τ h Σ h, where F Σ h is defined by (F, τ h )=γ(f, divτ h ).
11 2088 ZHIQIANG CAI AND YANQIU WANG In this section, we construct and analyze a multigrid solver for the discrete problem (5.1). Let T k, k =1,...,K, be a series of nested meshes under uniform refinements and assume that T 1 has characteristic mesh size O(1). Denote the characteristic mesh size of T k to be h k. Define Σ k = Σ hk, the RT/BDM finite element space based on T k. Notice that Σ k 1 Σ k. Similarly, we define spaces U k and V k. Denote div 1 k = div 1 h k. Denote Q k 1 : Σ k Σ k 1 and P k 1 : Σ k Σ k 1 to be the L 2 projection and the Λ projection, respectively. Let R k : Σ k Σ k be a symmetric positive definite operator, which will be defined later. Let m be the smoothing times on each level. Finally, denote A k to be A hk. We define a V-cycle multigrid iteration operator Mg k (, ) : Σ k Σ k Σ k inductively as follows. Algorithm 1. Set Mg 1 (τ, F )=Λ 1 1 F. Assuming that Mg k 1(, ) has been defined, define Mg k (τ, F ) as follows: set τ 0 = τ and define (1) τ l = τ l 1 + R k (F Λ k τ l 1 ) for l =1,...,m; (2) σ m = τ m + Mg k 1 (0, Q k 1 (F Λ k τ m )); (3) σ l = σ l 1 + R k (F Λ k σ l 1 ) for l = m +1,...,2m. Set Mg k (τ, F )=σ 2m. The V-cycle multigrid iteration for solving problem (5.1) over Σ k is defined below. Given an initial guess σ 0 Σ k, we defined a sequence approximating the solution σ h by σ i+1 = Mg k (σ i, F ), i =0, 1,... The error propagation operator ε : Σ k Σ k for the above iteration is [8] σ h σ i+1 = ε(σ h σ i )=Mg k (σ h σ i, 0). Similarly to the proof of Theorem 7.4 in [8] and Theorem 5.1 in [1], one can show that the multigrid solver with R k as the smoother and m smoothing times on each level will have the convergence rate (5.2) Λ(ετ, τ ) C p C p +2m Λ(τ, τ ) for all τ Σ k if the following two assumptions are true: (M.1) The spectrum of I R k Λ k is in [0, 1). That is, 0 Λ((I R k Λ k )τ, τ ) < Λ(τ, τ ) for all τ Σ k. (M.2) There exists a positive constant C p such that Λ((I P k 1 )τ, τ ) C p R k (Λ k τ, Λ k τ ) for all τ Σ k. In the remainder of this section, we will construct a smoother R k which satisfies assumptions (M.1) and (M.2). Let V k be the set of all vertices in T k, including the boundary vertices. For each x V k, define Ω k,x to be the interior of the union of all triangles in T k which have x as a vertex. Let Σ k,x be the RT/BDM space defined on T k such that each function in Σ k,x vanishes outside Ω k,x. Notice that Σ k,x (Σ k + span{i}), but Σ k,x Σ k in general, since Ω tr τ dx may not vanish for τ Σ k,x. Define operator Λ k,x : Σ k,x Σ k,x by (Λ k,x σ, τ )=Λ(σ, τ ) for all σ, τ Σ k,x.
12 MULTIGRID FOR STOKES PROBLEMS 2089 Notice that Λ k,x is symmetric and (Λ k,x σ, σ) = 0 if and only if σ = ci, where c is a constant. Since ci / Σ k,x and Σ k,x is a finite dimensional space, it is easy to see that Λ k,x is symmetric positive definite over Σ k,x, although the spectrum of Λ k,x may not have a uniform lower bound independent of all k and x. Let Q k,x : (Σ k + span{i}) Σ k,x and P k,x : (Σ k + span{i}) Σ k,x be the L 2 projection and the Λ projection, respectively. Here P k,x is well defined since the operator Λ k,x is symmetric positive definite. Define operator R k : Σ k (Σ k + span{i}) using an additive Schwarz scheme: (5.3) Rk = 1 P k,x Λ 1 k = 1 Λ 1 k,x 3 3 Q k,x. x V k x V k Define R k : Σ k Σ k by (5.4) R k = Q Σk Rk, where Q Σk : (Σ k + span{i}) Σ k is the L 2 projection. It is clear that R k is symmetric positive definite over Σ k under the L 2 inner product. Remark 2. For simplicity, in this paper we will analyze only the smoother R k as defined in (5.4). The analysis can easily be extended to the Hiptmair smoother [19], which is usually more time-efficient in terms of a single sweep of smoothing process, and the same multigrid convergence rate estimate holds. Next, we show that R k satisfies assumptions (M.1) and (M.2). The proof of assumption (M.1) follows directly from the Schwarz inequality and the fact that each triangle in T k is covered by exactly three subdomains in {Ω k,x } x Vk. In the following, we concentrate on the proof of assumption (M.2). First, we state the following result. Lemma 5.1. Assume the following assumption holds. (R.1) For all τ k Σ k, there exist a decomposition (I P k 1 )τ k = τ k,x, where τ k,x Σ k,x, x V k and a positive constant C 0 such that Λ(τ k,x, τ k,x ) C 0 Λ((I P k 1 )τ k, τ k ). x V k Then (M.2) holds with C p =3C 0. Proof. Assume (R.1) holds. By using the Schwarz inequality, we have Λ((I P k 1 )τ k, τ k )= Λ(τ k,x, τ k )= Λ(τ k,x, P k,x τ k ) x V k x V k ( x V k Λ(τ k,x, τ k,x ) ) 1/2 Λ 1/2 ( x V k P k,x τ k, τ k 3C 1/2 0 Λ 1/2 ((I P k 1 )τ k, τ k )Λ 1/2 ( R k Λ k τ k, τ k ) = 3C 1/2 0 Λ 1/2 ((I P k 1 )τ k, τ k )Λ 1/2 (R k Λ k τ k, τ k ). The last step comes from the fact that span{i} is in the kernel of Λ k. This completes the proof of the lemma. )
13 2090 ZHIQIANG CAI AND YANQIU WANG Finally, we need to prove that condition (R.1) holds. Lemma 5.2. Condition (R.1) holds with C 0 = c(1 + γ 2 ), where c is a positive constant independent of k or γ. Proof. Let σ k = (I P k 1 )τ k, where τ k Σ k. We decompose σ k into σ k = curl q k + div 1 k v k, where q k U k and v k = divσ k. Denote σ 1 k = div 1 k v k. Clearly, Λ(curl q k, curl q k )+Λ(σ 1 k, σ 1 k)=λ(σ k, σ k ). Therefore we can consider the decomposition of curl q k and σ 1 k separately. Next we need a partition of unity. Let θ x be a continuous piecewise linear function defined on T k such that the support of each θ x is in Ω k,x and θ x W j, (Ω) h j k, j =0, 1, where W j, (Ω) is the Sobolev space with index (j, ) [14]. This kind of partition of unity exists and has been widely used in the analysis of overlapping Schwarz domain decomposition methods. Clearly we have the decomposition curl q k = x V k curl (Π Uk (θ x q k )), where Π Uk is the natural interpolation onto U k associated with the degrees of freedom. By the standard scaling argument, it is not hard to see that curl (Π Uk (θ x q k )) 0,Ω curl (θ x q k ) 0,Ω. Combining the above, Lemma 4.2, and the fact that h k = h k 1 /2, we have Λ(curl (Π Uk (θ x q k )), curl (Π Uk (θ x q k ))) x V k (5.5) = (Acurl (Π Uk (θ x q k )), curl (Π Uk (θ x q k ))) x V k curl (Π Uk (θ x q k )) 2 0,Ω curl (θ x q k ) 0,Ω x V k x V k curl q k 2 0,Ω + h 2 k q k 2 0,Ω (Acurl q k, curl q k ) Λ(σ k, σ k ). Next, for σ 1 k, we consider the decomposition σ1 k = x V k Π Σk (θ x σ 1 k ), where Π Σ k is the interpolation operator as defined in (3.3) on the mesh T k. This decomposition is well defined since Π Σk is a linear operator. By the standard scaling argument, we have Π Σk (θ x σ 1 k) 0,Ω θ x σ 1 k 0,Ω. Therefore, by Lemma 4.2 and the fact that 2h k = h k 1, Λ(Π Σk (θ x σ 1 k), Π Σk (θ x σ 1 k)) x V k (5.6) ( Π Σk (θ x σ 1 k) 2 0,Ω + γ div(π Σk (θ x σ 1 k)) 2 0,Ω) x V k σ 1 k 2 0,Ω + γ( divσ 1 k 2 0,Ω + h 2 k σ1 k 2 0,Ω) (1 + γ 2 )Λ(σ k, σ k ).
14 Finally, we define the decomposition σ k = MULTIGRID FOR STOKES PROBLEMS 2091 x V k (curl (Π Uk (θ x q k ))+Π Σk (θ x σ 1 k)). The stability of this decomposition comes from (5.5) and (5.6). This completes the proof of the lemma. By inequality (5.2) and Lemmas 5.1 and 5.2, the convergence rate of the multigrid 1+γ solver is O( 2 1+γ 2 +m ). 6. Numerical results. We report some numerical results for the multigrid solver for the H(div)-type problem (5.1). Let Ω be the unit square (0, 1) (0, 1). The multigrid solver with a smoother as defined in section 5 is used. The coarsest mesh is as shown in Figure 6.1. We solve problem (5.1) on the 2nd 5th multigrid levels. The right-hand side F Σ h is selected randomly Fig The coarsest mesh. Both the lowest order RT element and the lowest order BDM element are used in the discretization. The experiments are done for γ =10 4, 10 2, 1, 10 2, We tested the iteration numbers needed for reaching a relative residual The results are listed in Tables 6.1 and 6.2. Table 6.1 Iteration numbers using the RT element. Level γ =10 4 γ =10 2 γ =1 γ =10 2 γ = Table 6.2 Iteration numbers using the BDM element. Level γ =10 4 γ =10 2 γ =1 γ =10 2 γ = From the tables we can see that when γ is large, the iteration numbers seem to be independent of γ and the number of levels, which agrees with the theory. When γ is
15 2092 ZHIQIANG CAI AND YANQIU WANG small, the theory suggests that the convergence rate deteriorates. From the numerical results we can see that this is true for the RT element. However, for the BDM element, the convergence rate still seems to be independent of γ. Appendix. Proof of (4.6). As stated before, we assume that Ω is a convex polygon. Let σ h and q h be defined as in Lemma 4.2. Then q h is in H 1 (Ω, R 2 )/R 2 and for all p H ŨH, wehave (A.1) (Acurl q h, curl p H )=Λ(curl q h, curl p H )=Λ(σ h, curl p H )=0. Denote n and t to be the unit outer normal vector and the unit tangential vector along Ω. Define the curl operator for all τ H 1 (Ω, M 2 )by ( τ 11 y curl τ = + ) τ12 x. τ21 y + τ22 x First we notice that the following lemma holds. Lemma A.1. Assume that there exists σ H 1 (Ω, M 2 ) satisfying q h = curl σ and that (S1) Aσ = σ, or equivalently, σ is trace-free; (S2) (σt) Ω = 0; (S3) σ 1,Ω q h 0,Ω. Then (4.6) holds. Proof. Define σ 1 by { (Aσ 1, τ )+(divτ, u) = 0 for all τ Σ, (divσ 1, v) =(divσ, v) for all v V. The regularity result (2.3) implies σ 1 1,Ω divσ 0,Ω. Since σ σ 1 H 1 (Ω, M 2 ) and div(σ σ 1 )=0, there must exist a potential function p H 2 (Ω, R 2 ) such that curl p = σ σ 1 and p 2,Ω σ 1,Ω + σ 1 1,Ω σ 1,Ω q h 0,Ω. Let p H be the natural interpolation of p onto ŨH. Then by the assumptions on σ, (A.1), and integration by parts, q h 2 0,Ω = (curl σ, q h )= (σ, curl q h )= (Aσ, curl q h ) = (Acurl p, curl q h )= (Acurl (p p H ), curl q h ) H p 2,Ω Acurl q h 0,Ω H q h 0,Ω Acurl q h 0,Ω. This completes the proof of (4.6). In the remainder of this section, we will construct a σ that satisfies the assumptions in Lemma A.1. For any scalar function f, define curl f =( f y, f x )t. Then there exists an L 2 - orthogonal Helmholtz decomposition (see [16]) q h = curl φ + ψ, where φ H0 1 (Ω) and ψ H 1 (Ω)/R satisfy φ 1,Ω q h 0,Ω and ψ 1,Ω q h 0,Ω. Then ( ) ( ) φ ψ φ (ψ + C) q h = curl = curl + curl u ψ φ (ψ + C) φ
16 MULTIGRID FOR STOKES PROBLEMS 2093 for all C R and u H 2 (Ω, R 2 ). Define ( ) φ (ψ + C) (A.2) σ = + u. (ψ + C) φ Clearly σ H 1 (Ω, M 2 ), and we need only find C and u such that σ satisfies (S1), (S2), and (S3) in Lemma A.1, which can be translated to (A.3) div u = 2φ, u t Ω =(ψ + C)n Ω, C q h 0,Ω, u 2,Ω q h 0,Ω. Since q h H 1 (Ω, R 2 )/R 2 and φ H0 1 (Ω), we have ψn ds = ψdx = (q h curl φ) dx Ω Ω Ω (A.4) = curl φdx = φt ds = 0. Ω Recall that Ω is a convex polygon. Denote z i, i =1,...,N, to be the corners of Ω. To close the polygon, we set z 0 = z N or z N+1 = z 1. Denote Γ i to be the edge of Ω connecting z i and z i+1. Define n i =(n 1,i,n 2,i ) t and t i =(n 2,i, n 1,i ) t to be the unit outer normal vector and the unit tangential vector on Γ i, respectively. According to the trace theorem on polygons (see Theorem in [17]), we have (A.5) ε 0 Ω s 1 ψ(z i + st i )+ψ(z i st i 1 ) 2 ds for all ε < min i=1,...,n z i z i 1. Let s 0, s be points on Ω and s s 0 dt be the line integral from s 0 to s along Ω in counterclockwise direction. Define (A.6) and C = Ω 2φdx Ω ( s s 0 ψn dt) n ds Ω ( s s 0 n dt) n ds s (A.7) g(s) = (ψ + C)n dt. s 0 By (A.4) and Green s formula, we immediately see that g is continuous on Ω. Furthermore, we clearly have g Γi H 3/2 (Γ i, R 2 ) and C φ 0,Ω + ψ 0, Ω φ 0,Ω + ψ 1,Ω q h 0,Ω, g 3/2,Γi ψ 1/2, Ω + C φ 0,Ω + ψ 1,Ω q h 0,Ω. By (A.5), (A.6), and (A.7), it is not hard to see that Ω 2φdx = g n ds and Ω ε s 1 g t (z i + st i ) n i 1 g 2 t (z i st i 1 ) n i ds = 0 ε 0 < s 1 (ψ + C)(z i + st i ) (ψ + C)(z i st i 1 ) 2 (n i 1 n i ) 2 ds
17 2094 ZHIQIANG CAI AND YANQIU WANG for all ε<min i=1,...,n z i z i 1. Recall that φ H 1 0 (Ω). Combining all the above and according to Theorem 7.1 in [4], there exists u H 2 (Ω, R 2 ) such that div u = 2φ, u Ω = g, and u 2,Ω φ 1,Ω + i=1,...,n g 3/2,Γi φ 1,Ω + ψ 1,Ω q h 0,Ω. Clearly u defined above and C defined in (A.6) satisfy (A.3), and hence σ defined as in (A.2) satisfies (S1), (S2), and (S3). This completes the proof for (4.6). REFERENCES [1] D. N. Arnold, R. S. Falk, and R. Winther, Preconditioning in H(div) and applications, Math. Comp., 66 (1997), pp [2] D. N. Arnold, R. S. Falk, and R. Winther, Multigrid in H(div) and H(curl), Numer. Math., 85 (2000), pp [3] D. N. Arnold, Jr., J. Douglas, and C. P. Gupta, A family of higher order mixed finite element methods for plane elasticity, Numer. Math., 45 (1984), pp [4] D. N. Arnold, L. R. Scott, and M. Vogelius, Regular inversion of the divergence operator with Dirichlet boundary conditions on a polygon, Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 15 (1988), pp [5] D. N. Arnold and R. Winther, Mixed finite elements for elasticity, Numer. Math., 92 (2002), pp [6] F. P. T. Baaijens, Mixed finite element methods for viscoelastic flow analysis: A review, J. Non-Newtonian Fluid Mech., 79 (1998), pp [7] M. Behr, L. P. France, and T. E. Tezduyar, Stabilized finite element methods for the velocity-pressure-stress formulation of incompressible flows, Comput. Methods Appl. Mech. Engrg., 104 (1993), pp [8] J. H. Bramble, Multigrid Methods, Pitman Res. Notes Math. Ser. 294, Longman Scientific & Technical, Harlow, UK, [9] F. Brezzi, J. Douglas, and L. D. Marini, Two families of mixed finite elements for second order elliptic problems, Numer. Math., 47 (1985), pp [10] F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods, Springer-Verlag, New York, [11] Z. Cai, B. Lee, and P. Wang, Least-squares methods for incompressible Newtonian fluid flow: Linear stationary problems, SIAM J. Numer. Anal., 42 (2004), pp [12] Z. Cai and G. Starke, First-order system least squares for the stress-displacement formulation: Linear elasticity, SIAM J. Numer. Anal., 41 (2003), pp [13] Z. Cai, C. Tong, P. Vassilevski, and C. Wang, Mixed Finite Element Methods for Incompressible Flow: Stationary Stokes Equations, manuscript, [14] P. G. Ciarlet, The Finite Element Method for Elliptic Problems, North Holland, Amsterdam, [15] M. I. Gerritsma and T. N. Phillips, Compatible spectral approximations for the velocitypressure-stress formulation of the Stokes problem, SIAM J. Sci. Comput., 20 (1999), pp [16] V. Girault and P. A. Raviart, Finite Element Methods for Navier Stokes Equations, Springer Ser. Comput. Math. 5, Springer-Verlag, New York, [17] P. Grisvard, Singularities in Boundary Value Problems, Research in Applied Mathematics 22, Springer-Verlag, Berlin, [18] W. Hackbusch, Multigrid Methods and Applications, Springer-Verlag, Berlin, [19] R. Hiptmair, Multigrid method for H(div) in three dimensions, Electron. Trans. Numer. Anal., 6 (1997), pp [20] C. Johnson and B. Mercier, Some equilibrium finite element methods for two-dimensional elasticity problems, Numer. Math., 30 (1978), pp [21] R. Keunings, Advances in the computer modeling of the flow of polymeric liquids, Comp. Fluid Dyn. J., 9 (2001), pp [22] J. Mandel, S. McCormick, and R. Bank, Variational multigrid theory, in Multigrid Methods, S. F. McCormick, ed., SIAM, Philadelphia, 1987, pp
18 MULTIGRID FOR STOKES PROBLEMS 2095 [23] P. A. Raviart and J. M. Thomas, A mixed finite element method for 2nd order elliptic problems, in Mathematical Aspects of Finite Element Methods, Lecture Notes in Math. 606, I. Galligani and E. Magenes, eds., Springer-Verlag, 1977, pp [24] P. S. Vassilevski and J. Wang, Multilevel iterative methods for mixed finite element discretizations of elliptic problems, Numer. Math., 63 (1992), pp [25] J. Xu, Iterative methods by space decomposition and subspace correction, SIAM Rev., 34 (1992), pp
Lecture 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 informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
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 informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
MATHEMATICS OF COMPUTATION Volume 75, Number 256, October 2006, Pages 1659 1674 S 0025-57180601872-2 Article electronically published on June 26, 2006 ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED
More informationNONCONFORMING MIXED ELEMENTS FOR ELASTICITY
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 295 307 c World Scientific Publishing Company NONCONFORMING MIXED ELEMENTS FOR ELASTICITY DOUGLAS N. ARNOLD Institute for Mathematics
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 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 informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
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 informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationA u + b u + cu = f in Ω, (1.1)
A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL Abstract. This paper presents analysis of a weighted-norm least squares finite element method for elliptic
More informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
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 informationc 2008 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 46, No. 3, pp. 640 65 c 2008 Society for Industrial and Applied Mathematics A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
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 informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
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 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 informationFIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PLANAR LINEAR ELASTICITY: PURE TRACTION PROBLEM
SIAM J. NUMER. ANAL. c 1998 Society for Industrial Applied Mathematics Vol. 35, No. 1, pp. 320 335, February 1998 016 FIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PLANAR LINEAR ELASTICITY: PURE TRACTION
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationarxiv: v2 [math.na] 23 Apr 2016
Improved ZZ A Posteriori Error Estimators for Diffusion Problems: Conforming Linear Elements arxiv:508.009v2 [math.na] 23 Apr 206 Zhiqiang Cai Cuiyu He Shun Zhang May 2, 208 Abstract. In [8], we introduced
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
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 informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationOverlapping Schwarz Preconditioners for Spectral. Problem in H(curl)
Overlapping Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Technical Report TR2002-83 November 22, 2002 Department of Computer Science Courant Institute of Mathematical
More informationLeast-squares Finite Element Approximations for the Reissner Mindlin Plate
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl., 6, 479 496 999 Least-squares Finite Element Approximations for the Reissner Mindlin Plate Zhiqiang Cai, Xiu Ye 2 and Huilong Zhang
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 informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationarxiv: v3 [math.na] 8 Sep 2015
A Recovery-Based A Posteriori Error Estimator for H(curl) Interface Problems arxiv:504.00898v3 [math.na] 8 Sep 205 Zhiqiang Cai Shuhao Cao Abstract This paper introduces a new recovery-based a posteriori
More informationNumerical Analysis and Computation of Fluid/Fluid and Fluid/Structure Interaction Problems Jason S. Howell
J.S.Howell Research Statement 1 Numerical Analysis and Computation of Fluid/Fluid and Fluid/Structure Interaction Problems Jason S. Howell Fig. 1. Plot of the magnitude of the velocity (upper left) and
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
More informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More informationarxiv: v1 [math.na] 19 Dec 2017
Arnold-Winther Mixed Finite Elements for Stokes Eigenvalue Problems Joscha Gedicke Arbaz Khan arxiv:72.0686v [math.na] 9 Dec 207 Abstract This paper is devoted to study the Arnold-Winther mixed finite
More informationA Robust Approach to Minimizing H(div)-Dominated Functionals in an H 1 -Conforming Finite Element Space
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2000; 00:1 23 [Version: 2002/09/18 v1.02] A Robust Approach to Minimizing H(div)-Dominated Functionals in an H 1 -Conforming Finite
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 informationA posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationMODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 138 145 MODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT JIŘÍ MARYŠKA, OTTO SEVERÝN AND MARTIN VOHRALÍK Abstract. A stochastic discrete
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationA NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES EQUATIONS
MATHEMATICS OF COMPUTATION Volume 74, Number 250, Pages 543 554 S 0025-5718(04)01711-9 Article electronically published on August 31, 2004 A NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationA DECOMPOSITION RESULT FOR BIHARMONIC PROBLEMS AND THE HELLAN-HERRMANN-JOHNSON METHOD
Electronic ransactions on Numerical Analysis. Volume 45, pp. 257 282, 2016. Copyright c 2016,. ISSN 1068 9613. ENA A DECOMPOSIION RESUL FOR BIHARMONIC PROBLEMS AND HE HELLAN-HERRMANN-JOHNSON MEHOD WOLFGANG
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 informationSTABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS
STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS RAYTCHO LAZAROV AND XIU YE Abstract. In this paper, we derive two stabilized discontinuous finite element formulations, symmetric
More informationError estimates for the Raviart-Thomas interpolation under the maximum angle condition
Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation
More informationMIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS
MIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS DOUGLAS N. ARNOLD, RICHARD S. FALK, AND JAY GOPALAKRISHNAN Abstract. We consider the finite element solution
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 informationFind (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ.
A POSTERIORI ESTIMATES FOR THE STOKES EIGENVALUE PROBLEM CARLO LOVADINA, MIKKO LYLY, AND ROLF STENBERG Abstract. We consider the Stokes eigenvalue problem. For the eigenvalues we derive both upper and
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
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 informationA SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY
A SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY J. GOPALAKRISHNAN AND J. GUZMÁN Abstract. We presented a family of finite elements that use a polynomial space augmented
More informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationA priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces
A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces A. Bespalov S. Nicaise Abstract The Galerkin boundary element discretisations of the
More informationMultigrid Methods for Maxwell s Equations
Multigrid Methods for Maxwell s Equations Jintao Cui Institute for Mathematics and Its Applications University of Minnesota Outline Nonconforming Finite Element Methods for a Two Dimensional Curl-Curl
More informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
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 informationA Finite Element Method Using Singular Functions: Interface Problems
A Finite Element Method Using Singular Functions: Interface Problems Seokchan Kim Zhiqiang Cai Jae-Hong Pyo Sooryoun Kong Abstract The solution of the interface problem is only in H 1+α (Ω) with α > 0
More informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
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 informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationQUADRILATERAL H(DIV) FINITE ELEMENTS
QUADRILATERAL H(DIV) FINITE ELEMENTS DOUGLAS N. ARNOLD, DANIELE BOFFI, AND RICHARD S. FALK Abstract. We consider the approximation properties of quadrilateral finite element spaces of vector fields defined
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More informationMATH 676. Finite element methods in scientific computing
MATH 676 Finite element methods in scientific computing Wolfgang Bangerth, Texas A&M University Lecture 33.25: Which element to use Part 2: Saddle point problems Consider the stationary Stokes equations:
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 informationA mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N)
wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad
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 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 informationA WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS
A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This article introduces and analyzes a weak Galerkin mixed finite element method
More informationEXACT DE RHAM SEQUENCES OF SPACES DEFINED ON MACRO-ELEMENTS IN TWO AND THREE SPATIAL DIMENSIONS
EXACT DE RHAM SEQUENCES OF SPACES DEFINED ON MACRO-ELEMENTS IN TWO AND THREE SPATIAL DIMENSIONS JOSEPH E. PASCIAK AND PANAYOT S. VASSILEVSKI Abstract. This paper proposes new finite element spaces that
More informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
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 information1. Fast Solvers and Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl)
DDM Preprint Editors: editor1, editor2, editor3, editor4 c DDM.org 1. Fast Solvers and Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Bernhard Hientzsch 1 1. Introduction.
More informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationDISCRETE EXTENSION OPERATORS FOR MIXED FINITE ELEMENT SPACES ON LOCALLY REFINED MESHES
DISCRETE EXTENSION OPERATORS FOR MIXED FINITE ELEMENT SPACES ON LOCALLY REFINED MESHES MAR AINSWORTH, JOHNNY GUZMÁN, AND FRANCISCO JAVIER SAYAS Abstract. The existence of uniformly bounded discrete extension
More informationDOUGLAS N. ARNOLD, RICHARD S. FALK, AND R. WINTHER. Dedicated to Professor Ivo Babuška on the occasion of his seventieth birthday.
MATHEMATICS OF COMPUTATION Volume 66, Number 219, July 1997, Pages 957 984 S 0025-5718(97)00826-0 PRECONDITIONING IN H (div) AND APPLICATIONS DOUGLAS N ARNOLD, RICHARD S FALK, AND R WINTHER Dedicated to
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 informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions Juan Gabriel Calvo Preprint No. 25-2017 PRAHA 2017 A VIRTUAL
More informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
More informationANDREA TOSELLI. Abstract. Two-level overlapping Schwarz methods are considered for nite element problems
OVERLAPPING SCHWARZ METHODS FOR MAXWELL'S EQUATIONS IN THREE DIMENSIONS ANDREA TOSELLI Abstract. Two-level overlapping Schwarz methods are considered for nite element problems of 3D Maxwell's equations.
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 informationA UNIFORMLY ACCURATE FINITE ELEMENT METHOD FOR THE REISSNER MINDLIN PLATE*
SIAM J. NUMER. ANAL. c 989 Society for Industrial and Applied Mathematics Vol. 26, No. 6, pp. 5, December 989 000 A UNIFORMLY ACCURATE FINITE ELEMENT METHOD FOR THE REISSNER MINDLIN PLATE* DOUGLAS N. ARNOLD
More informationA NOTE ON CONSTANT-FREE A POSTERIORI ERROR ESTIMATES
A NOTE ON CONSTANT-FREE A POSTERIORI ERROR ESTIMATES R. VERFÜRTH Abstract. In this note we look at constant-free a posteriori error estimates from a different perspective. We show that they can be interpreted
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationA Second Elasticity Element Using the Matrix Bubble
Portland State University PDXScholar Mathematics and Statistics Faculty Publications and Presentations Fariborz Maseeh Department of Mathematics and Statistics 2011 A Second Elasticity Element Using the
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 information