REGULARITY AND MULTIGRID ANALYSIS FOR LAPLACE-TYPE AXISYMMETRIC EQUATIONS. Hengguang Li. IMA Preprint Series #2417.
|
|
- Oliver Williams
- 5 years ago
- Views:
Transcription
1 REGULARITY AND MULTIGRID ANALYSIS FOR LAPLACE-TYPE AXISYMMETRIC EQUATIONS By Hengguang Li IMA Preprint Series #417 February 013) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY OF MINNESOTA 400 Lind Hall 07 Church Street S.E. Minneapolis, Minnesota Phone: Fax: URL:
2 REGULARITY AND MULTIGRID ANALYSIS FOR LAPLACE-TYPE AXISYMMETRIC EQUATIONS HENGGUANG LI Abstract. Consider axisymmetric equations associated with Laplace-type operators. We establish full regularity estimates in high-order Kondrat vetype spaces for possible singular solutions due to the non-smoothness of the domain and to the singular coefficients in the operator. Then, we show suitable graded meshes can be used in high-order finite element methods to achieve the optimal convergence rate, even when the solution is singular. Using these results, we further propose multigrid V-cycle algorithms solving the system from linear finite element discretizations on optimal graded meshes. We prove the multigrid algorithm is a contraction, with the contraction number independent of the mesh level. Numerical tests are provided to verify the theorem. 1. Introduction Let Ω R + := {r, z), r > 0} be a bounded polygonal domain in the rz-plane and the integer β = 0 or 1. Let Γ 0 be the interior part of Ω {r = 0} and Γ 1 := Ω\Γ 0. See Figure 1 for example. We consider equations with parameter β, 1) L β u β := r + r 1 r + z βr β )u β = f β in Ω, with the β-dependent boundary conditions ) 3) u 0 = 0 on Γ 1, r u 0 = 0 on Γ 0 ; u 1 = 0 on Ω. Equation 1) plays an important role in the study of various equations in 3D axisymmetric domains. For instance, denote by R 3 Ω := Ω Γ 0 ) [0, π) the domain obtained by the rotation of Ω meridian domain) about the z-axis Figure 1). Recall the Cartesian coordinates x, y, z) and the cylindrical coordinates r, θ, z) on Ω. For any ṽ L Ω), its Fourier coefficients with respect to θ are defined by 4) v r, z) = F ṽ) := 1 π π π Let w be the solution of the 3D Poisson equation 5) ṽr, θ, z)e iθ dθ, Z. w = f in Ω, w = 0 on Ω. Then, the first three Fourier coefficients w, 1, of w are given by equation 1) as follows: w 0 = u 0, given f 0 = F 0 f) β = 0); w 1 = u 1, given f 1 = F 1 f) β = 1); and w 1 = u 1, given f 1 = F 1 f) β = 1). In particular, when the 3D data are axisymmetric, u 0 is the meridian trace of w and L 1 u 1 = f 1 is regarded as Date: February 1, 013. H. Li was supported in part by the NSF grant DMS
3 H. LI Figure 1. The 3D axisymmetric domain Ω left); the meridian polygonal domain Ω right). the azimuthal Stoes equation [1, ]. We also mention that L 1 coincides with the principal part of the azimuthal Maxwell operator [6, 15]. This dimension reduction from 3D to D has the potential for substantial computational savings in numerical approximations of equations on 3D axisymmetric domains. However, the coordinate transformation introduces new differential operators with singular coefficients as in 1) and new function spaces with singular or degenerate weights. The development of robust numerical methods for these equations calls for rigorous numerical analysis, which has recently drawn much attention from the scientific community. For example, [6] provided a comprehensive introduction on spectral methods for various axisymmetric equations. Under the assumption that the solution has sufficient regularity, finite element analysis for axisymmetric problems was discussed in [1, 5, 1, 1, 6]. In the case that the solution possesses singularities, the singular expansion of the solution for equation 1) was studied in [17, 7]. As the most relevant results, we mention the papers [3] and [15], both for equation 1) with β = 0. In [3], a second-order regularity estimate was established in a class of weighted Kondrat ev-type) spaces for singular solutions. Consequently, the author proposed new linear finite element methods approximating singular solutions in the optimal rate. In [15], the multigrid V-cycle algorithm was studied for finite element spaces augmented with non-polynomial functions. The result was derived only for solutions with required regularity. This paper contains our systematic study for equation 1) on its regularity, finite element approximation, and multigrid analysis. Motivated by [18, 4, ], we first introduce high-order Kondrat ev-type weighted spaces Definition.) to handle possible singular solutions from the non-smoothness of the domain and from the singular coefficients. Using appropriate isomorphic mappings in weighted spaces, we prove the full-regularity estimates Theorem 3.5) for singular solutions of equation 1). Based on this regularity result and properties of interpolation operators in weighted spaces, we give an explicit construction of nested graded meshes, such that the associated high-order finite element methods achieve the optimal convergence rate for singular solutions Theorem 3.16). Then, we analyze the system from the linear finite element discretization on optimal graded meshes. With the growth rate of the condition numbers for nested subspaces Lemma 4.5), we provide smoothing properties of the Richardson smoother and approximation properties of the numerical solution in various weighted spaces. This leads to our convergence result Theorem 4.1) on the multigrid V-cycle algorithm for singular solutions of equation 1).
4 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 3 For equation 1), our high-order regularity result and optimal high-order finite element methods approximating singular solutions are new. These results generalize the results in [3] low-order regularity and linear finite element results only for the case β = 0). The uniform convergence of the proposed multigrid method on graded meshes is also new for the axisymmetric equation 1), which recovers the classical multigrid result for elliptic equations with full-regularity solutions [7, 11]. We mention that high-order Fourier coefficients ) of w in 5) are determined by equations similar to equation 1) with β = 1 [6]. Therefore, our analysis shall provide building blocs for approximating general D Fourier coefficients of the 3D problem 5). The multigrid method studied here shall be a crucial ingredient in developing efficient multigrid solves for more complex axisymmetric problems. The rest of the paper is organized as follows. In Section, we introduce weighted Sobolev spaces for our problem. We also summarize estimates for functions in these spaces. In Section 3, we establish high-order regularity results in Kondrat evtype weighted spaces for the possible singular solution of equation 1). Then, we give high-order finite element methods approximating the singular solutions in the optimal rate. In Section 4, we first describe the multigrid V-cycle algorithm for equation 1). Then, we show the proposed multigrid algorithm converges uniformly, independent of the mesh level, even when the solution is singuar. In Section 5, we include multiple numerical tests for the multigrid method to verify the theory. Acnowledgements. The author would lie to than Jay Gopalarishnan for discussions and suggestions on this research.. Preliminaries and notation We introduce appropriate weighted spaces to study equation 1). Useful estimates regarding these spaces will be collected in the second part of this section..1. Function spaces. We first recall a class of weighted Sobolev spaces from [6]. Definition.1. Type I weighted spaces). For an integer m 0, define L 1Ω) := {v, v rdrdz < }, H1 m Ω) := {v, c α v L 1Ω), α m}, Ω where the muti-index α = α 1, α ) is a pair of nonnegative integers, α = α 1 + α, and c α = r α1 z α. The norms and the semi-norms for any v H1 m Ω) are v H1 m Ω) := c α v) rdrdz, v H1 m Ω) := c α v) rdrdz. α m Ω Furthermore, we define two spaces H m + Ω) and H m Ω). For H m + Ω), if m is not even, α =m H m + Ω) := {v H m 1 Ω), i 1 r v {r=0} = 0, 1 i < m }, 6) v Hm + Ω) = v H m 1 Ω); if m is even, besides the condition in 6), we require Ω m 1 r v) r 1 drdz < for any v H+ m Ω), and the corresponding norm is v H m + Ω) = v H1 mω) + r m 1 v) r 1 drdz ) 1/. Ω Ω
5 4 H. LI For H m Ω), if m is not odd, 7) H m Ω) := {v H1 m Ω), r i v {r=0} = 0, 0 i < m 1 v Hm Ω) = v H m 1 Ω); if m is odd, besides the condition in 7), we require Ω m 1 r v) r 1 drdz <, for any v H m Ω), and the corresponding norm is v H m Ω) = v H1 mω) + r m 1 v) r 1 drdz ) 1/. Thus, we denote different subspaces: Ω } H 1 +,0Ω) := H 1 +Ω) {v Γ1 = 0}, H 1,0Ω) := H 1 Ω) {v Ω = 0}. We now introduce another type of weighted spaces for our analysis on singular solutions of equation 1). Definition.. Type II weighted spaces). Let Q i be the ith vertex of Ω and define the vertex set Q := {Q i } I i=1. Denote by L the smallest distance from a vertex to any disjoint edge of Ω. Let Bx, r 0 ) be the ball centered at x with radius r 0. Let ϑ C Ω\Q) be a function, such that ϑ = x Q i in V i := Ω BQ i, L/) and ϑ L/ in Ω\ I i=1 V i. Note that V i and V j are disjoint if i j. Thus, we define for µ R and for any open set G Ω, with the semi-norm and norm K m µ,1g) := {v, ϑ µ+ α α c v L 1G), α m}. v Kµ,1 m G) := ϑ m µ c α v L G), 1 α =m v K m µ,1 G) := m In addition, for 0 l m, we define the subspaces of K m µ,1ω): l=0 v K l µ,1 G). 8) Kµ,+Ω) m := {v Kµ,1Ω), m Ω i 1 r ϑ µ+l v) ) r 1 drdz <, 1 i l }, 9) Kµ, Ω) m := {v Kµ,1Ω), m Ω i r ϑ µ+l v) ) r 1 drdz <, 0 i l 1 }. The corresponding norms are v K m µ,+ Ω) = v Kµ,1 m Ω) + l m 1 i l Ω v K m µ, Ω) = v Kµ,1 m Ω) + l m 0 i l 1 i 1 r ϑ µ+l v) ) r 1 drdz) 1/, Ω i r ϑ µ+l v) ) r 1 drdz) 1/. Recall the Fourier coefficients v of ṽ in 4). We then have the following isomorphism Theorem II.3.1) [6]. Proposition.3. For m 0, ṽ v ) Z defines an isomorphic mapping 10) H m Ω) Π Z H m ) Ω), where the space H) m Ω) is defined as follows. For > m 1, w H m ) Ω) = w H1 mω) + m r m w L Ω). 1
6 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 5 If m 1, for 0 j 1, H m ) Ω) = Hm + Ω) { j rw Γ0 H m ) Ω) = Hm Ω) { j rw Γ0 = 0} and w H m ) Ω) = w H m + Ω) = 0} and w H m ) Ω) = w H m Ω) even ); odd ). Remar.4. For ṽ H m Ω), H) m Ω) are natural spaces to study the traces of its Fourier coefficients on Ω. We shall show that the Kondrat ev-type spaces Kµ,1Ω) m are the appropriate spaces to study the trace when ṽ only has reduced regularity. Note that for m 0, H m 0) Ω) = Hm + Ω), H m ±1) Ω) = Hm Ω). Thus, we shall consider equation 1) in H m + Ω) if β = 0 and in H m Ω) if β = 1. Inheriting from the 3D Poisson equation 5), we have the boundary condition u β Γ1 = 0. The restrictions on u β in type I weighted spaces lead to different boundary conditions on Γ 0 := Ω\Γ 1 in ) and 3). Since the Fourier coefficients may be complex even though the original function is real, we shall denote the complex extensions of all the spaces above without change of notation. The Fourier transform 4) leads to the variational formulation for equation 1) [6]: Find u 0 H+,0Ω) 1 for any v 0 H+,0Ω) 1 resp. u 1 H,0Ω) 1 for v 1 H,0Ω)), 1 such that 11) a β u β, v β ) := r u β r v β + z u β z v β + β r u βv β )rdrdz = f β, v β, Ω where f 0 H 1 +,0Ω) resp. f 1 H 1,0Ω) ), the dual space of H 1 +,0Ω) resp. H 1,0Ω)) with the pivot space L 1Ω). For f β L 1Ω), the right hand side is the L 1-inner product. Then, we have the well-posedness result. Corollary.5. The variational formulation 11) defines a unique solution u 0 H 1 +,0Ω) resp. u 1 H 1,0Ω)) and there is a constant C > 0, independent of u β and f β, such that 1) u 0 H 1 + Ω) C f 0 H 1 +,0 Ω) and u 1 H 1 Ω) C f 1 H 1,0 Ω). Proof. For β = 0, the well-posedness follows from the isomorphism 10) and the 3D Poincaré inequality; for β = 1, it follows from the fact that a β, ) is coercive and continuous in H,0Ω). 1 The estimates in 1) are immediate consequences of the well-posedness... Weighted estimates. We now recall several estimates from []. These estimates give relations between various spaces and will be used to formulate our high-order regularity results in the next section. We distinguish the vertices on and off the z-axis as follows. A vertex on the z-axis will be denoted by Q z ; A vertex away from the z-axis will be denoted by Q r. Recall the neighborhood V := BQ, L/) Ω of the vertex Q. For a vertex Q z, let Γ V = Γ 0 V and we denote by Ṽ := V Γ V) [0, π) Ω the neighborhood of Q z in Ω. For an integer l 1, V/l := B Q, L/l) ) Ω and Ṽ/l := V/l Γ V/l) [0, π). Throughout the paper, by H, we mean the dual space of H. As in Definition.1, we also use the multi-index α = α 1, α, α 3 ) for a 3D domain, such that α = α 1 + α + α 3 and α = x α1 y α z α3. The generic constant C > 0 in our analysis below may be different at different occurrences. It will depend on the computational
7 6 H. LI domain, but not on the functions involved in the estimates or the mesh level in the finite element algorithms. The following lemma, derived as Lemmas.8 and.9 in [], contains useful weighted estimates in usual Sobolev spaces in the 3D neighborhood Ṽ of a vertex Q z. We refer to [] for the detailed proof. Lemma.6. Let Ṽ Ω be the neighborhood of a vertex Q z and ρ be the distance to Q z. Suppose α m ρ a+ α α v L <, for m 0 and a R. Then, for Ṽ) any 0 l m, C ρ a+l v H l Ṽ) ρ a+ α α v L Ṽ) C ρ a+l v H l Ṽ). α m l m Recall the multi-index α = α 1, α ) and the notation α c from Definition.1. Then, the following lemma concerns the connection between the two types of weighted spaces in the D neighborhood V of a vertex Q z in the rz-plane Lemmas.10 and.11 in []). Lemma.7. Let V Ω be the neighborhood of a vertex Q z and ρ be the distance to Q z. Suppose v K m a,1v). Then, for 0 l m and a R, C ρ a+l v H l 1 V) v K m a,1 V) C l m ρ a+l v H l 1 V). The following two lemmas concern the local behavior of functions in weighted spaces near the vertices. Lemma.8. Recall ϑ in Definition.. In the neighhood V of a vertex Q r away from the z-axis, for any v H 1 +,0Ω) resp. v H 1,0Ω)), we have ϑ 1 v L 1 V) C v H 1 + V), resp. ϑ 1 v L 1 V) C v H 1 V)). Proof. On V, both H 1 + and H 1 resp. L 1) are equivalent to the usual Sobolev space H 1 resp. L ), since r is bounded away from 0. Therefore, it suffices to show for any v H 1 V) {v Γ1 = 0}, 13) ϑ 1 v L V) C v H 1 V). However, the estimate in 13) is well nown based on a local Poincaré inequality. See [0, 19, 5, 4]. Lemma.9. Recall ϑ in Definition.. Let v H 1 0 Ω) and let Ṽ = V Γ V) [0, π) be the 3D neighborhood of a vertex Q z on the z-axis. Then, ϑ 1 v L Ṽ) C v H 1 Ṽ). Proof. Ṽ can be characterized in the spherical coordinates ρ, θ, φ) centered at Q z : Ṽ = {ρ, ω), 0 < ρ < L/, ω ω Q z}, where ω Q z S is the polygonal domain on the unit sphere S. Then, for any v H 1 Ṽ) {v Ω = 0}, v = vx + vy + vz = vρ + v φ ρ + vθ ρ sin φ,
8 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 7 and v ds C vφ + v θ ω Q z ω Q z sin ) sin φdθdφ, φ which is just the Poincaré inequality on ω Q z and ds = sin φdθdφ is the volume element on ω Q z. Thus, we obtain Ṽ v ρ dxdydz = L/ 0 ω Q z which completes the proof. v dsdρ C L/ 0 ω Q z v ρ + v φ ρ + vθ ) ρ ρ sin dsdρ φ = C v dxdydz, 3. Regularity and finite element analysis The solution u β for equation 1) is well defined in weighted spaces Corollary.5). When the 3D domain Ω is smooth, using the isomorphism in 10), we will have full regularity estimates in type I weighted spaces. These estimates will not hold when Ω is a general polygonal domain, since the regularity of u β is determined by the smoothness of Ω and by the operator L β. In this section, we show type II weighted spaces are suitable to study singular solutions. We derive high-order regularity estimates in these spaces for singular solutions and propose finite element algorithms approximating singular solutions in the optimal rate Regularity analysis. Before formulating our main regularity result Theorem 3.5), we first need several local regality estimates. All these results are written for different values of the parameter β, namely, 0 or 1. Lemma 3.1. Near a vertex Q r away from the z-axis, there exists η β > 0, such that for any 0 a β < η β, if f β K m a β 1,1V), m 0, the solution u β of equation 11) satisfies 14) u 0 K m+ a 0 +1,1 V/) C f 0 K m a0 1,1 V) + f 0 H 1 +,0 Ω) ) 15) u 1 K m+ a 1 +1,1 V/) C f 1 K m a1 1,1 V) + f 1 H 1,0 Ω) ). Proof. We only show 14), since the proof of 15) follows similarly. We apply a localization argument. Let ζ be a smooth cutoff function, such that ζ = 1 on V/ and ζ = 0 outside V. Thus, ζu 0 has the Dirichlet boundary condition on V. Then, 16) where r + r 1 r + z)ζu 0 = ζf 0 + g 0 in V, g 0 = u 0 r + z)ζ + r ζ r u 0 + z ζ z u 0 ) + r 1 u 0 r ζ. Note that the support of g 0 is within V\V/ and K m µ,1 = H m on V\V/. Therefore, using the interior regularity estimate in H m [14] and 1), we have g 0 K m a0 1,1 V) = g 0 K m a0 1,1 V\V/) C g 0 H m V\V/) C u 0 H V\V/) C f 0 H m V\V/4) + f 0 H 1 +,0 Ω) ) C f 0 K m a0 1,1 V\V/4) + f 0 H 1 +,0 Ω) ) By Corollary.5, equation 16) has a unique solution ζu 0 H 1 +V). Since r is bounded away from 0 on V, ζu 0 H 1 V) and the space K m µ,1v) is the same as the Ṽ
9 8 H. LI Kondrat ev space K m µ V) Definition 1.1 in [3]). In addition, the regularity of ζu 0 is determined by the principle part r z of the operator in 16). Then, based on the D regularity estimates in the Kondrat ev spaces for the Laplace operator [18, 4, 16], for the solution ζu 0 H 1 V), there exists η 0 > 0, such that for any 0 a 0 < η 0, ζu 0 K m+ a 0 +1,1 V) C ζf 0 + g 0 K m a0 1,1 V) C ζf 0 K m a0 1,1 V) + g 0 K m a0 1,1 V) ) C f 0 K m a0 1,1 V) + f 0 H 1 +,0 Ω) ). The lemma is thus proved due to the definition of the function ζ. Remar 3.. Near the vertex Q r, since r is away from 0, the principle part of the operator L β coincides with the Laplace operator r z. Therefore, the regularity index η β is the same for β = 0 and β = 1. Let ω be the interior angle at Q r. Then, 17) η 0 = η 1 = π/ω, which is the first eigenvalue of the operator pencil associated with the Laplace operator [16, 18, 19, 4]. Define ũ β x, y, z) := 1 π u β r, z)e iβθ and f β x, y, z) := 1 π f β r, z)e iβθ. Then, we see that, by the change of variables, ũ β satisfies the 3D Poisson equation, 18) ũ β = f β in Ω ũ β = 0 on Ω, β = 0, 1. In addition, its Fourier coefficients F ũ β ) = { uβ if = β, 0 if β. We now give a regularity estimate for ũ β near a vertex on the z-axis. Lemma 3.3. For a vertex Q z on the z-axis, there is η β > 0 for all m 0, such that for any 0 a β < η β, the solution ũ β of equation 18) satisfies ϑ aβ 1+ α α ũ β L Ṽ/))1/ α m+ C α m ϑ a β+1+ α α fβ L Ṽ))1/ + f β H 1 Ω)). Proof. We use a localization augment similar to the one in Lemma 3.1. Let ζ be a smooth cutoff function, such that ζ = 1 on Ṽ/ and ζ = 0 outside Ṽ. Then, 19) where ζũ β = ζ f β + g β, g β = ũ β ζ + x ũ β x ζ + y ũ β y ζ + z ũ β z ζ. Since g β vanishes near Q z, using the well-posedness of the Poisson problem 18), the usual interior regularity estimate [14], and the expression of g β above, we have ϑ aβ+1+ α α g β ) 1/ L C ũβ Ṽ) H Ṽ\Ṽ/)) 0) α m C f β H m 1 Ṽ\Ṽ/4) + f β H 1 Ω) ) C ϑ aβ+1+ α α fβ L Ṽ))1/ + f β H 1 Ω)). α m
10 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 9 Based on the 3D regularity estimates in the Kondrat ev spaces for equation 19) Theorem 1. in []) and 0), there exists η β > 0, such that for any 0 a β < η β, the solution ζũ β satisfies ϑ aβ 1+ α α ζũ β ) L Ṽ))1/ α m+ C α m C α m C α m ϑ a β+1+ α α ζ f β + g β ) L Ṽ))1/ ϑ a β+1+ α α ζ f β ) L Ṽ))1/ + α m ϑ a β+1+ α α fβ L Ṽ))1/ + f β H 1 Ω)). The lemma is thus proved due to the definition of ζ. ϑ a β+1+ α α g β L Ṽ))1/) Remar 3.4. Different from the situation in Lemma 3.1 and Rema 3., near the vertex Q z, the index η β is determined by the operator pencil of the 3D Laplace operator in a bounded region on the unit sphere [16, 18]. To be more precise, in the neighborhood Ṽ of Qz, let ω be the projection of Ṽ on the unit sphere. We write equation 18) in the polar coordinations, ũ β = ρ ρ ρ ) +ρ ρ + )ũ β, where is the Laplace-Beltrami operator on ω, = cot φ φ + φ + sin φ θ. For any vr, z) L 1V), define ṽ β := 1 π ve iβθ. The smallest eigenvalue λ β of with respect to ṽ β is given by λ 0 ṽ 0 = ṽ 0 = cot φ φ + φ )ṽ 0, λ 1 ṽ 1 = ṽ 1 = cot φ φ + φ sin φ)ṽ 1, which is determined by the associated Legendre s differential equation [13, 7]. Then, the explicit formula for the index η β near Q z is 1) η β = λ β + 1/4. In contrast to the case near Q r in 17), η 0 and η 1 may be different near Q z, since λ β are eigenvalues associated with different differential operators. Combining the local estimates in the lemmas above, we derive the global regularity estimate for equation 1). Theorem 3.5. Let u 0 H 1 +,0Ω) and u 1 H 1,0Ω) be the solutions of equation 1) with β = 0 and 1, respectively. For m 0, there exist Υ 0 > 0 and Υ 1 > 0, such that for any 0 a 0 < Υ 0 and 0 a 1 < Υ 1, we have ) 3) u 0 K m+ a 0 +1,1 Ω) C f 0 K m a0 1,+ Ω), ϑ a1 r 1 u 1 L 1 Ω) + u 1 K m+ a 1 +1,1 Ω) C f 1 K m a1 1, Ω), given that f 0 K m a 1,+Ω) and f 1 K m a 1, Ω). Proof. Let η i β be the upper bound of the parameter a β for the ith vertex Q i in Lemmas 3.1 and 3.3. Define Υ 0 = min i η i 0) and Υ 1 = min i η i 1). We only show 3), since the proof for ) follows similarly.
11 10 H. LI Recall that the weighted space K m a,1 resp., K m a,+ and K m a, ) is equivalent to the weighted space H m 1 resp., H m + and H m ) in a subdomain Ω sub Ω that is away from the vertex set. Based on the isomorphism in 10), the well-posedness and the usual interior regularity estimate for the 3D Poisson equation 18), we have 4) ϑ a1 r 1 u 1 L 1 Ω sub ) + u 1 K m+ a 1 +1,1 Ω sub) C f 1 Hm Ω ) + f 1 H 1 Ω) C f 1 K m a1 1, Ω ) + f 1 K 0 a1 1, Ω) ), where Ω sub Ω Ω such that Ω is away from the vertex set, and Ω = Ω Ω {r = 0}) ) [0, π) is from the rotation of Ω about the z-axis. Let V be the neighborhood of a vertex Q r away from the z-axis. By Lemma 3.1, the fact that r is bounded away from 0 on V, Lemma.8, and a 1 0, we have ϑ a1 r 1 u 1 L 1 V/) + u 1 K m+ a 1 +1,1 V/) C f 1 K m a1 1,1 V) + f 1 H 1,0 Ω) ) 5) = C f 1 K m a1 1,1 V) + sup v H 1,0 Ω),v 0 f 1, v) L 1 / v H 1,0 Ω)) C f 1 K m a1 1, V) + f 1 K 0 a1 1, Ω) ). Now, let V be the small neighborhood of a vertex Q z on the z-axis. By Lemma.7, we first have for any 0 l m, 6) ϑ 1 a1+l f 1 H l 1 V) C f 1 K m a1 1,1 V). Then, for f 1 Ka m 1 1, V), 6) and the condition in 9) i r ϑ 1 a1+l f 1 ) ) r 1 drdz <, 0 i l 1)/ Ω lead to ϑ 1 a1+l f 1 H V). l Then, by Lemma.6, the isomorphism in 10), and the definitions of the weighted spaces in 7) and 9), we have ϑ 1 a1+ α α f1 L Ṽ) C ϑ 1 a1+l f1 H l Ṽ) 7) α m C 0 l m 0 l m ϑ 1 a1+l f 1 H l V) C f 1 K m a 1 1, V). Then, by Lemma.7, the isomorphism in 10), Lemma.6, Lemma 3.3, Lemma.9, 7), and a 1 0, we have 8) ϑ a1 r 1 u 1 L 1 V/) + u 1 K m+ a 1 +1,1 V/) C ϑ a1 r 1 u 1 L 1 V/) + ϑ a1 1+l u 1 H1 l V/))1/) C l m+ C α m C α m l m+ ϑ a1 1+l ũ 1 H l Ṽ/))1/ ϑ a1+1+ α α f1 L Ṽ))1/ + sup v H 1 0 Ω),v 0 ϑ a1+1+ α α f1 L Ṽ))1/ + ϑ a1+1 f1 L Ω) C f 1 K m a1 1, V) + f 1 K 0 a1 1, Ω) ). Then, the proof is completed by combining 4), 5), and 8). f ) 1, v) L Ω)/ v H1 Ω) )
12 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 11 Remar 3.6. We used the global parameter Υ β in Theorem 3.5 in order to simplify the exposition. The proof in fact shows that Theorem 3.5 holds as long as for the ith vertex Q i, we choose the parameter 0 a β < ηβ i, where ηi β is determined as in Remars 3. and 3.4. This will result in the analogue of Theorem 3.5 in Kondrat ev-type spaces with vector subindices, in which the function ϑ is allowed to have different exponents near different vertices to capture the local behavior of the solution. See [4] for such a presentation for elliptic equations. 3.. Finite element analysis. Let T n = {T i } be a triangulation of Ω with shaperegular triangles T i. Denote by V β n the continuous Lagrange finite element space of order m with the corresponding Dirichlet boundary conditions described in ) and 3), respectively. Then, we have [1] 9) V 0 n H 1 +,0Ω) and V 1 n H 1,0Ω). In turn, the finite element solution u β n V β n for equation 11) is defined by 30) a β u β n, vn) β = f β vnrdrdz, β vn β V β n. Ω We denote a point in the rz-plane by x = r, z). Let Π : C 0 Ω) V β be the usual nodal interpolation operator, such that Πvx i ) = vx i ), where x i is the ith node in the triangulation. Note that we use the same notation for the interpolation operators into different spaces V β n, β = 0 and 1. In the text below, by a triangle T, we mean the closed set including both the interior and boundary of the triangle. We first recall two approximation results from Lemmas 6.1, 6. and 6.3 in [6], which are based on embedding theorems in weighted spaces. Lemma 3.7. Let x i be a node in the triangulation, which is on the z-axis. Let K i be the union of triangles that contain x i. Then, for v H 1 K i ), we have 31) for v H 1 K i ) H 1 Ω), 3) v Πv H 1 1 K i) Ch m K i v H 1 K i) ; r 1 v Πv) L 1 K i) Ch m K i v H where h Ki denotes the diameter of K i and m 1. 1 K i), Lemma 3.8. Let T T n be a triangle that does not intersect the z-axis. Suppose min x T rx) Ch T. Then, for v H1 T ), we have 33) v Πv H 1 1 T ) Ch m T v H 1 T ), 34) r 1 v Πv) L 1 T ) Ch m T v H 1 T ), where h T denotes the diameter of T and m 1. Remar 3.9. By Lemma 3.7 and Lemma 3.8, one can derive the convergence result for the finite element solution in 30) on quasi-uniform meshes u 0 u 0 n H 1 + Ω) Ch m u 0 H + Ω) and u 1 u 1 n H 1 Ω) Ch m u 1 H Ω), given that u β is sufficiently regular and h is the mesh size. As discussed above, the regularity of the solution, however, depends on the differential operator L β and the domain Ω. When the solution u β is singular e.g., u β / H1 Ω)), the linear finite element method on the quasi-uniform mesh only gives a suboptimal convergence rate [6].
13 1 H. LI Based on our regularity estimates in Theorem 3.5 and Remar 3.6, we use the following process to generate a sequence of graded triangulations approximating the singular solutions of equation 1). Definition The κ-refinement). Let κ 0, 1/] be the grading parameter and T be a triangulation of Ω such that no two vertices of Ω belong to the same triangle of T. Then the κ-refinement of T, denoted by κt ), is obtained by dividing each edge AB of T into two parts as follows. If neither A nor B is in the vertex set Q, then we divide AB into two equal parts. Otherwise, if A Q, we divide AB into AC and CB such that AC = κ AB. This will divide each triangle of T into four triangles and leads to a finer mesh. Suppose the initial mesh T 0 satisfies the above conditions. Then the nth level mesh with κ-refinement is obtained recursively by T n = κt n 1 ), n = 1,,.... This graded process was proposed in [3] for the linear finite element method approximating the axisymmetric Poisson equation. We here generalize it to highorder finite element methods for the new equation 1). We need the following notation to carry out the analysis on graded meshes. Let n be the number of κ-refinements of the domain. Thus, the final triangulation is T n. Let T i,j T j, j n, be the union of triangles in T j that contain the vertex Q i Q. Note that T i,j T i,l for j l and i T i,j occupies the neighborhood of the vertex set Q in the triangulation T j. Recall the regularity estimate for the solution and the parameter Υ β in Theorem 3.5. We choose the grading parameter in the κ-refinement, for different values of β, 35) κ β = min1/, m/a β ), for any 0 < a β < Υ β, where m 1 is the degree of piecewise polynomials in the finite element space V β n. Then, based on analysis on T n \ i T i,0, on i T i,j 1 \ i T i,j, and on T i,n, our error estimates are summarized in the following lemmas. Lemma For κ β defined in 35), let U T n be the union of triangles that intersect T n \ i T i,0. Then, u 0 Πu 0 n H 1 + T n\ it i,0) C nm u 0 H 1 U) u 1 Πu 1 n H 1 T n\ it i,0) C nm u 1 H 1 U). Proof. Assume U is away from the vertices of the domain this is true when n > ). Then, based on Definition 3.10, the mesh size on U is O n ). Summing up the estimates in 31), 3), 33), and 34) completes the proof. For the estimates on T i,0, the union of initial triangles containing the vertex Q i, we consider the new coordinate system that is a simple translation of the old rz-coordinate system, now with Q i at the origin. Then, for a subset G T i,0 and 0 < λ < 1, we define the dilations of G and of a function as follows, for r, z) G, 36) G λ := G/λ = {r λ, z λ ) = λ 1 r, λ 1 z)}, v λ r λ, z λ ) := vr, z). Then, we have the following scaling estimates from Lemma 4.5, []. Lemma 3.1. Suppose G λ is in the neighborhood V of Q i. If Q i is on the z-axis, v λ K m a,1 G λ ) = λ a 3/ v K m a,1 G), r 1 λ v λ L 1 G λ ) = λ 1/ r 1 v L 1 G);
14 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 13 if Q i is not on the z-axis, Cλ a 1 v K m a,1 G) v λ K m a,1 G λ ) Cλ a 1 v K m a,1 G). Now, we are ready to give estimates on the region T i,j 1 \T i,j. Lemma Let U T n be the union of triangles that intersect G := T i,j 1 \T i,j. Suppose G V. Let h be the mesh size on U and ξ = sup x G ϑx). Let a 0. Then, for v K a+1,1 U), 37) and for v H 1 U) K a+1,1 U), 38) v Πv H 1 + G) Cξ a h/ξ) m v K a+1,1 U), v Πv H 1 G) Cξ a h/ξ) m v K a+1,1 U). Proof. Recall the new coordinate system with Q i as the origin. Recall the dilations in 36) and the parameter L in Definition.. Note that Πv) λ = Πv λ ). Then, we choose λ = ξ/l, such that G λ V. If Q i is on the z-axis, by Lemma 3.1, the definitions of the weighted spaces, 31), and 3), we have r 1 v Πv) L 1 G) + v Πv H 1 1 G) C r 1 v Πv) L 1 G) + v Πv K 1 1,1 G) = λ 1/ r 1 λ v λ Πv λ )) L 1 G λ ) + v λ Πv λ ) K 1 1,1 G λ )) Cλ 1/ r 1 λ v λ Πv λ )) L 1 G λ ) + v λ Πv λ ) H 1 1 G λ )) Cλ 1/ h/λ) m v λ H 1 U λ ) Cλ1/ h/λ) m v λ K 1,1 U λ) Ch/ξ) m v K 1,1 U) Cξa h/ξ) m v K U). a+1,1 If Q i is not on the z-axis, the proof is similar. With the corresponding estimate in Lemma 3.1, the definitions of the weighted spaces, 33), and 34), we have r 1 v Πv) L 1 G) + v Πv H 1 1 G) C r 1 v Πv) L 1 G) + v Πv K 1 1,1 G) C v Πv K 1 1,1 G) C v λ Πv λ ) K 1 1,1 G λ ) C v λ Πv λ ) H 1 1 G λ ) Ch/λ) m v λ H 1 U λ ) Ch/λ) m v λ K 1,1 U λ) Ch/ξ)m v K 1,1 U) Cξa h/ξ) m v K U). a+1,1 Thus, the estimate in 38) is proved. The estimate in 37) follows similarly using Lemma 3.1, the definitions of the weighted spaces, 31), and 33). Lemma For κ β defined in 35), let U T n be the union of triangles that intersect G := T i,j 1 \T i,j. Then, for 0 < a β < Υ β, u 0 Πu 0 H 1 + G) C nm u 0 K U), a 0 +1,1 u 1 Πu 1 H 1 G) C nm u 1 K U). a 1 +1,1 Proof. Definition 3.10 shows that the mesh size on U is Oκ j 1 β j 1 n ). Using the notation of Lemma 3.13, we have ξ = Oκ j 1 β ) on T i,j 1 \T i,j. Therefore, using
15 14 H. LI Lemma 3.13, we have u 1 Πu 1 ) H 1 G) Cκ j 1)a1 1 j 1 n)m u 1 K a 1 +1,1 U) C j 1)m j 1 n)m u 1 K a 1 +1,1 U) = C nm u 1 K a 1 +1,1 U). Then, we have proved the second estimate in this lemma. The first estimate can be proved similarly by Lemma 3.13 and the observation on the mesh size for the regions G and U. The following lemma gives the error bounds on the last patch T i,n of triangles that have the vertex Q i as the common node. Lemma For κ β defined in 35), let U T n be the union of triangles that intersect T i,n. Suppose T i,0 V near Q i. Then, u 0 Πu 0 H 1 + T i,n) C nm u 0 K a 0 +1,1 U) u 1 Πu 1 ) H 1 T i,n) C nm u 1 K a 1 +1,1 U) + ϑ a1 r 1 u 1 L 1 U)). Proof. Definition 3.10 shows that the mesh on U has the size Oκ n β ). Let ζ be a smooth function, such that ζ = 0 in the small neighborhood of the vertex Q i and ζ = 1 on a region containing all the other nodes in T i,n. Let v = u 1 ζu 1. Note that Πu 1 = Πζu 1 ) and T i,0 = T i,n /λ, where λ = κ n β. If Q i is on the z-axis, by Lemma 3.1, we have and r 1 ζu 1 L 1 T i,n) = λ 1/ r 1 λ ζ λu 1 ) λ L 1 T i,0) Cλ 1/ r 1 λ u 1) λ L 1 T i,0) = C r 1 u 1 L 1 T i,n), ζu 1 K m 1,1 T i,n) = λ 1/ ζ λ u 1 ) λ K m 1,1 T i,0) Cλ 1/ u 1 ) λ K m 1,1 T i,0) = C u 1 K m 1,1 T i,n). Therefore, by Lemma.9, 10), Lemma 3.7, Lemma 3.1, and the estimates above, r 1 u 1 Πu 1 ) L 1 T i,n) + u 1 Πu 1 H 1 1 T i,n) r 1 v L 1 T i,n) + r 1 ζu 1 Πu 1 ) L 1 T i,n) + v H 1 1 T i,n) + ζu 1 Πu 1 K 1 1,1 T i,n) 39) C r 1 u 1 L 1 T i,n) + u 1 K 1 1,1 T i,n) +λ 1/ r 1 λ ζ λu 1 ) λ Πu 1 ) λ ) L 1 T i,0) + ζ λ u 1 ) λ Πu 1 ) λ K 1 1,1 T i,0)) C r 1 u 1 L 1 T i,n) + u 1 K 1 1,1 T i,n) + λ 1/ ζ λ u 1 ) λ H 1 U λ ) ) C r 1 u 1 L 1 T i,n) + u 1 K 1 1,1 T i,n) + λ 1/ u 1 ) λ K 1,1 U λ) ) = C r 1 u 1 L 1 T i,n) + u 1 K 1,1 U)). If Q i is not on the z-axis, r is bounded away from 0, by Lemma 3.1, r 1 ζu 1 L 1 T i,n) = r 1 ϑ)ϑ 1 ζu 1 L 1 T i,n) C ζu 1 K 0 1,1 T i,n) C ζ λ u 1 ) λ K 0 1,1 T i,0) C u 1 ) λ K 0 1,1 T i,0) C u 1 K 0 1,1 T i,n); ζu 1 K m 1,1 T i,n) C ζ λ u 1 ) λ K m 1,1 T i,0) C u 1 ) λ K m 1,1 T i,0) C u 1 K m 1,1 T i,n).
16 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 15 Thus, by Lemma.8, Lemma 3.8, Lemma 3.1, and the estimates above, r 1 u 1 Πu 1 ) L 1 T i,n) + u 1 Πu 1 H 1 1 T i,n) r 1 v L 1 T i,n) + r 1 ζu 1 Πu 1 ) L 1 T i,n) + v H 1 1 T i,n) + ζu 1 Πu 1 H 1 1 T i,n) 40) C u 1 K 1 1,1 T i,n) + ζu 1 Πu 1 K 1 1,1 T i,n)) C u 1 K 1 1,1 T i,n) + ζ λ u 1 ) λ Πu 1 ) λ K 1 1,1 T i,0)) C u 1 K 1 1,1 T i,n) + ζ λ u 1 ) λ K 1,1 Ti,0)) C u 1 K 1,1 Ti,n). Note that max x Ti,n ϑx) = Oκ n β ). Combining 39) and 40), we have r 1 u 1 Πu 1 ) L 1 T i,n) + u 1 Πu 1 H 1 1 T i,n) C r 1 u 1 L 1 T i,n) + u 1 K 1,1 U)) Cκ na1 1 ϑ a1 r 1 u 1 L 1 T i,n) + u 1 K a 1 +1,1 U)) C nm ϑ a1 r 1 u 1 L 1 T i,n) + u 1 K a 1 +1,1 U)). This completes the proof of the second estimate of the lemma. The first estimate follows from a similar process using Lemmas 3.1,.8,.9, 3.7, and 3.8. The local estimates above lead to the main result on the convergence rate of our finite element method. Theorem Let u β be the solution of equation 1) and u β n V β n be the finite element solution defined in 30). Recall κ β from 35) for 0 < a β < Υ β as defined in Theorem 3.5. Then, for f 0 Ka m 1 Ω) and f 0 1,+ 1 Ka m 1 Ω), 1 1, u 0 u 0 n H 1 + Ω) CdimV 0 n) m/ f 0 K m 1 a 0 1,+ Ω), u 1 u 1 n H 1 Ω) CdimV 1 n) m/ f 1 K m 1 a 1 1, Ω). where dimv β n) = O4 n ) is the dimension of V β n and m 1. Proof. The proof follows by summing up the estimates in Lemmas 3.11, 3.14, and 3.15, together with the regularity results in Theorem 3.5. Remar The graded meshes from the κ-refinements are nested and consist of shape-regular triangles. Various numerical tests for the linear finite element methods m = 1) approximating singular solutions in the case β = 0 can be found in [3], which convincingly verify Theorem As mentioned in Remar 3.6, instead of a uniform α β for all vertices, we can choose different regularity indices a β near different vertices Q i. This will result in different grading parameters κ β close to the vertices. Nevertheless, as long as 0 < a β < ηβ i for each Q i, we shall have the optimal convergence rate indicated in Theorem This flexibility can help to improve the shape regularity of the triangulation. 4. The multigrid algorithm From now on, we concentrate on the multigrid analysis for the linear finite element approximation of equation 1) on graded meshes given in the last section.
17 16 H. LI 4.1. The multigrid V-cycle algorithm. We denote the β-dependent linear finite element spaces V β 0 Vβ 1 Vβ K that are defined on graded meshes T in the last section, 0 K, where the parameter κ β is chosen as in 35) for m = 1. For simplicity, we assume that there is only one singular vertex Q for the solution. Namely, the solution is in H1 or ηβ i > 1) except in the neighborhood of Q. Therefore, the graded mesh κ β < 0.5) is implemented only for triangles touching Q, while the usual quasiuniform κ = 0.5) decomposition is performed for other triangles. If we modify the weight function in the definition of the space Kµ,1Ω), m letting ϑ be the distance to Q in its neighborhood V and be 1 in the neighborhoods of other vertices, then Theorem 3.16 in this case reads 41) 4) u 0 u 0 H 1 + Ω) CdimV 0 ) 1/ f 0 K 0 a0 1,+ Ω), u 1 u 1 H 1 Ω) CdimV 1 ) 1/ f 1 K 0 a1 1, Ω). In this section, we use this modified function ϑ and also define := r, z ) t. Let T T be the union of all the triangles in T touching Q. Define the jth layer of T, 43) L j := T j 1 \T j, 0 j, where T 1 := T 0. Then, we define the piecewise-constant function on T 44) ω Lj := κ β ) j, where κ β = 1/a β < 0.5 is the grading parameter for Q. Then, we define the following mesh-dependent weighted inner product 45) v, w ) := Ω ω v w rdrdz, v, w V β, and the norm induced by the inner product, v := ωv rdrdz, v V β. Ω Let I 1 : Vβ 1 Vβ be the coarse-to-fine operator, which is the natural injection. The fine-to-coarse operator I 1 : V β Vβ 1 is the transpose with respect to the inner product in 45), I 1 v, w 1 ) 1 := v, I 1w 1 ), v V β, w 1 V β 1. Let A β : Vβ Vβ be the operator associated with equation 30), 46) A β v, w ) := a β v, w ), v, w V β. In addition, let f β Vβ be the function, such that f β, v ) = f β, v ) L 1 Ω), v V β. Thus, we are ready to formulate the abstract multigrid V-cycle algorithm for equation 1) with the Richardson smoother.
18 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 17 Algorithm 4.1. The V-cycle algorithm). The th level V-cycle algorithm produces MG, f β, z 0, l) as an approximate solution for A β uβ = f β with initial guess z 0 V β, where l denotes the number of pre-smoothing and postsmoothing steps. For = 0, we define MG0, f β 0, z 0, l) = A β 0 ) 1 f β 0. For 1, the approximate solution MG, f β, z 0, l) is computed recursively in three steps: Pre-smoothing. For 1 j l, compute z j by z j = z j 1 + γ f β Aβ z j 1), where the constant γ is the damping factor that we will specify later. Coarse grid correction. Let g 1 = I 1 f β Aβ z l) be the restriction of the residue on the 1st level. Define q V β 1 by Then, we compute z l+1 by q := MG 1, g 1, 0, l). z l+1 = z l + I 1q. Post-smoothing. For l + j l + 1, compute z j by z j = z j 1 + γ f β Aβ z j 1). We then define MG, f β, z 0, l) = z l+1. Remar 4.. Let E : V β Vβ be the error-propagation operator for the th level V-cycle algorithm defined above. Namely, E u β z 0) = u β MG, f β, z 0, l). The relation below follows from a straightforward calculation: 47) E = R l Id I 1P 1 + I 1E 1 P 1 )R l, where Id is the identity operator on V β, R : V β Vβ measuring the effect of the smoothing step is defined by 48) R := Id γ A β, and P 1 : V β Vβ 1 is the transpose of I 1 i.e., for any v V β, 49) with respect to the bilinear forms, a β P 1 v, w 1 ) = a β v, I 1w 1 ), w V β 1. Note that P 1 is in fact the projection on V β 1 with respect to a β, ). Then, for any v 1, w 1 V β 1, we have a β P 1 I 1v 1, w 1 ) = a β I 1v 1, I 1w 1 ) = a β v 1, w 1 ). Then, for any v V β 1, w V β, 50) a β I 1 v 1, Id I 1P 1 ) )w = a β I 1v 1, w ) a β I 1v 1, I 1P 1 w ) = a β v 1, P 1 w ) a v 1, P 1 I 1P 1 w ) ) = 0.
19 18 H. LI Figure. Reference triangles ˆT 1 and ˆT corresponding to different triangles touching the z-axis. For the sae of economy of notation, we use the same notation for the operators Id, I 1, I 1, E, R, and P 1, even though V 0 and V1 are different subspaces of H1 1 Ω) see 9)). Let {x i } be the set of nodes in the triangulation T and d i, be the diameter of the support of the basis function associated with the node x i. Throughout the text, by A B, we mean there are generic constants C 1, C > 0, such that C 1 B A C B. For a triangle T, r max T ) := max x T rx) and r min T ) := min x T rx). For a triangle T T away from the z-axis, we denote by ˆT its standard reference triangle of diameter one. If T touches the z-axis, its reference triangle ˆT will be either ˆT 1 or ˆT Figure ), to preserve the intersection set between T and the z-axis. For v V β and x T, we define ˆv ˆx) := v x), where ˆx ˆT is the image of x after the affine mapping. We now investigate necessary properties of the inner products and of the operators in Algorithm 4.1. Lemma 4.3. Let T i T be a triangle and h i, be its diameter. Then, vrdrdz h i, v x l ) max dl,, rx l ) ), v V β, T i x l T i where x l, 1 l 3, is the vertex of T i and d l, is defined as above. Proof. Case I: T i {r = 0} =. Thus, r max T i )/r min T i ) < M with a uniform constant M > 0 for any T i. We map T i to the reference triangle ˆT. Based on the norm equivalence on finite dimensional spaces, we have vrdrdz r max T i ) vdrdz = r max T i ) T i ˆv dˆrdẑ T i T i ˆT r max T i )h i, ˆv ˆx l ) h i, v x l ) max dl,, rx l ) ), x l T i ˆx l ˆT where we used the fact that r max T i ) max d l,, rx l ) ) for any triangle away from the z-axis. Case II: T i {r = 0} =. Let {λ m } m=1,,3 be the barycentric coordinates associated with the three vertices of T i. 1. T i {r = 0} is a line segment.) We map T i to the the first reference ˆT 1. Without loss of generality, we assume the first node of T i is away from the z-axis.
20 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 19 Let r 1 be the distance from its first node to the z-axis. Note that r 1 h i,. Based on the norm equivalence on finite dimensional spaces, we have vrdrdz = r 1 vλ 1 drdz = r 1 T i ˆv ˆλ 1 dˆrdẑ T i T i ˆT 1 h 3 i, ˆv ˆx l ) h i, v x l ) max dl,, rx l ) ). ˆx l ˆT 1 x l T i. T i {r = 0} is a point.) We map T i to the second reference ˆT. Suppose that its third vertex is on the z-axis. Denote by r 1 and r the distances form the first and second vertices of T i to the z-axis, respectively. Then, for any x T i, rx) = r 1 λ 1 x) + r λ x) with r 1 r h i,. Using the norm equivalence on finite element spaces, we have vrdrdz = vr 1 λ 1 + r λ )drdz h i, T i ˆv ˆλ 1 + ˆλ )dˆrdẑ T i T i ˆT h 3 i, ˆv ˆx l ) h i, v x l ) max dl,, rx l ) ). ˆx l ˆT x l T i Hence, the lemma is proved for any T i T. Lemma 4.4. The weighted inner product in 45) satisfies v, v ) dimv β ) 1 i v x i ) max d i,, rx i ) ), v V β. Proof. Recall dimv β ) 4. Thus, for any v V β, we have dimv β ) 1 dimv β ) i=1 dimv β ) i=1 v x i ) max d i,, rx i ) ) κ) Ii) Ii) ) κ Ii) )v x i ) max d i,, rx i ) ), where the function I : N N {0} is such that Ii) = j if the node x i is in the jth layer L j of T. By Definition 3.10, for a node x i L j, the diameter of a triangle T l touching x i satisfies h l, j κ j. Therefore, using Lemma 4.3, we have dimv β ) 1 dimv β ) i=1 v x i ) max d i,, rx i ) ) κ) Ii) Ii) ) κ Ii) )v x i ) max di,, rx i ) ) T l T x i T l ω Tl h l, v x i ) max di,, rx i ) ) T l T x i T l ω Tl vrdrdz = ωv rdrdz = v, v ), T l T T l Ω which completes the proof.
21 0 H. LI Lemma 4.5. The spectral radus of A β defined in 46) satisfies 0 < ρa β ) CdimVβ ). Proof. According to the definition, A β is symmetric positive definite SPD) with respect to the inner product, ). Therefore, it has a discrete spectrum and all its eigenvalues are real and positive. Thus, we only need to show for any V β v 0, 51) A β v, v ) v, v ) CdimV β ). Recall that near the singular vertex Q, the grading parameter κ β = 1/a β < 0.5. If a triangle T i is in the jth layer L j of T, 0 j, its diameter h i, κ j j, and so is the diameter d l, of the support of the basis function associated with the vertex x l of T i. Case I: T i {r = 0} =. Recall that r max T i )/r min T i ) M is bounded with M > 0 for any T i. Then, i, r maxt i ) v v )rdrdz r max T i ) v v drdz Ch T i T i r max T i ) 5) vx l ) vx l )rx l ), x l T i x l T i T i v drdz where we used the usual H 1 -L inverse inequality in the finite element space and the norm equivalence on finite dimensional spaces. Note r max T i ) 1 Ch 1 i,. By the norm equivalence on finite dimensional spaces, r 1 vdrdz r max T i ) 1 vdrdz = r max T i ) 1 T i ˆv dˆrdẑ T i T i ˆT Ch i, ˆv ˆx l ) C v 53) x l ) max dl,, rx l ) ). x l T i ˆx l ˆT Case II: T i {r = 0} =. Let {λ m } m=1,,3 be the barycentric coordinates associated with the three vertices of T i. 1. T i {r = 0} is a line segment.) We map T i to ˆT 1. Suppose the first node of T i is away from the z-axis. Let r 1 be the distance from the first node to the z-axis. Then, r 1 h i,. Based on the norm equivalence on finite dimensional spaces, 54) ˆv ˆv )ˆλ 1 dˆrdẑ + ˆv ˆλ 1 dˆrdẑ C ˆv dˆrdẑ, ˆT 1 ˆT 1 ˆT 1 since each side above defines the square of a norm for ˆv. Then, v v )rdrdz = r 1 v v )λ 1 drdz T i T i r 1 T i h i, ˆv ˆv )ˆλ 1 dˆrdẑ Cr 1 ˆv dˆrdẑ ˆT 1 ˆT 1 r 1 ˆv ˆx l ) v 55) x l ) max dl,, rx l ) ), ˆx l ˆT 1 x l T i where we used the scaling argument, the norm equivalence on finite dimensional spaces, and 54).
22 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 1 For v V 1, v Ti is a linear function that vanishes on the z-axis. Therefore, v Ti = ν r for a constant ν. Then, we have r 1 vdrdz = νrdrdz = r 1 T i ν ˆλ 1 dˆrdẑ h 3 i,ν T i T i ˆT 1 h i, vx l ) v 56) x l ) max dl,, rx l ) ), x l T i x l T i where we used the norm equivalence on finite dimensional spaces.. T i {r = 0} is a point.) We map T i to ˆT. Suppose that its first two vertices are not on the z-axis. Denote by r 1 and r the distances from the first vertex and the second vertex of T i to the z-axis, respectively. For any point x T i, rx) = r 1 λ 1 x) + r λ x) and r 1 r h i,. Using the norm equivalence on finite dimensional spaces, we have 57) ˆv ˆv )ˆλ 1 + ˆλ )dˆrdẑ + ˆT ˆv ˆλ 1 + ˆλ )dˆrdẑ C ˆT ˆv dˆrdẑ, ˆT since each side of 57) defines the square of a norm for ˆv. Similarly, we have v v )rdrdz h i, T i h i, ˆv ˆv )ˆλ 1 + ˆλ )dˆrdẑ T i ˆT Ch i, ˆv dˆrdẑ h i, ˆv ˆx l ) v 58) x l ) max dl,, rx l ) ). ˆT ˆx l ˆT x l T i For v V 1, we can write the linear function v Ti = l ν lφ l, where for l = 1,, ν l is a constant and φ l is the linear basis function associated with the lth vertex of T i. It was shown that T i r 1 φ l drdz < in [1]. Therefore, T i r 1 v drdz is bounded. We then have r 1 vdrdz = r 1 λ 1 + r λ ) 1 vdrdz h 1 i, T i ˆλ 1 + ˆλ ) 1ˆv dˆrdẑ T i T i ˆT h i, ˆv ˆx l ) = h i, vx l ) v 59) x l ) max dl,, rx l ) ), ˆx l ˆT x l T i x l T i where the ey step is to use the norm equivalence between ˆT ˆλ 1 +ˆλ ) 1ˆv dˆrdẑ) 1/ and ˆx l ˆT ˆv ˆx l) ) 1/. This can be done because the first term is finite for all linear functions that vanish on the z-axis. Hence, based on 5), 55), and 58), we have for any T i T 60) v v )rdrdz C v x l ) max dl,, rx l ) ), v V β. T i x l T i 61) Summing up 53), 56), and 59), we obtain that for any v V 1 and T i T, r 1 vdrdz C v x l ) max dl,, rx l ) ). T i x l T i
23 H. LI Consequently, using 60), 61), and Lemma 4.4, for any v V β, we have A β v, v ) = a β v, v ) = v v + βr v)rdrdz T T i T i C v x l ) max dl,, rx l ) )) T i T x l T i v x i ) max di,, rx i ) ) dimv )v, v ). 1 i dimv β ) Remar 4.6. Based on Lemma 4.5, we can choose the damping factor 6) γ = O dimv β ) 1) for the multigrid V-cycle, such that the spectral radius 63) 0 < ργ A β ) The multigrid analysis. We now concentrate on the convergence analysis of the multigrid V-cycle defined in Algorithm 4.1 for the axisymmetric equation 1). We first have the smoothing property for the Richardson smoother. Lemma 4.7. Recall the operator R from 48). Then, a β Id R )R l ) 1 v, v l a β Id R l ) )v, v, v V β. Proof. Since R = Id γ A β, by 63), for any 0 j i, we have a β Id R )Rv i ), v = γ RA i β v, A β v ) γ R j Aβ v, A β v ) = a β Id R )R j v, v ). Then, the lemma is proved by the following telescopic cancellation l)a β Id R )R l v, v ) l 1 j=0 a β Id R )R j v, v ) a β Id R l )v, v ). Lemma 4.8. The error-propagation operator E in 47) is symmetric positive semi-definite with respect to a β, ) for 0. Proof. We prove it by induction. For = 0, the lemma holds. Assume the statement holds for 1, 1. We first show the symmetry of E. By 47), 48), 49), and the assumption for E 1, we have a β E v, w ) = a β R l Id P 1 + E 1 P 1 )R l v, w ) = a β R l v, Id P 1 )Rw l ) + aβ E 1 P 1 Rv l, Rw l ) = a β R l v, Id P 1 )Rw l ) + aβ R l v, E 1 P 1 Rw l ) = a β v, R l Id P 1 + E 1 P 1 )R l w ) = aβ v, E w ).
24 REGULARITY AND MG METHODS FOR AXISYMMETRIC EQUATIONS 3 We now show that E is positive and semi-definite. By 47), 49), and the assumption on E 1, a β E v, v ) = a β Rv l, Rv l ) a β Id E 1 )P 1 Rv l, P 1 Rv l ) = a β Id P 1 )Rv l, Id P 1 )Rv l ) +a β E 1 P 1 Rv l, P 1 Rv l ) 0. In the next two lemmas, we give approximation results needed in the analysis. Lemma 4.9. Let E be H 1 +,0Ω) for β = 0 and be H 1,0Ω) for β = 1. For any v V β and 1, we have Id I 1P 1 )v CdimV β ) 1/ Id I 1P 1 )v E. Proof. We now give the proof for the case β = 1. Consider the auxiliary variational problem: Find ζ H,0Ω), 1 such that 64) a 1 ζ, ξ) = ωid I 1P 1 )v ξrdrdz, ξ H,0Ω), 1 Ω where a 1 ζ, ξ) = L 1 ζ, ξ) L 1 Ω) is the bilinear form in 11) and ω is the function defined in 44). Let ζ 1 V β 1 be the finite element solution of 64) on the 1st level mesh. Recall the layer L j of the triangulation T in 43). Then, since the regularity index 0 < a 1 < 1, by Theorem 3.5, 4), 44), and 64), we have 65) ζ I 1ζ 1 H 1 Ω) = ζ ζ 1 H 1 Ω) CdimV1 1) 1 L 1 ζ K 0 a 1 1, Ω) dimv 1 ) 1 L 1 ζ K 0 a 1 1, Ω) = dimv1 ) 1 CdimV 1 ) 1 j=0 ϑ 1 a1 L 1 ζ L 1 Lj) j=0 κ 1 a1)j 1 L 1 ζ L 1 Lj) = CdimV1 ) 1 = CdimV 1 ) 1 Id I 1P 1 )v. Then, by 50), 65), and Lemma 4.4, we have Id I 1P 1 )v = a 1 ζ, Id I 1P 1 ) )v = a 1 ζ I 1 ζ 1, Id I 1P 1 ) )v ζ I 1ζ 1 H 1 Ω) Id I 1P 1 )v H 1 Ω) CdimV 1 ) 1/ Id I 1P 1 Thus, we obtain for β = 1, Id I 1P 1 The proof for β = 0 follows similarly. κ j 1 j L 1 ζ L 1 Lj) j=0 )v Id I 1P 1 )v H 1 Ω). )v CdimV 1 ) 1/ Id I 1P 1 )v H 1 Ω). Lemma Recall the space E from Lemma 4.9. For any v V β Id I 1P 1 )v E CdimV β ) 1/ A β v. and 1,
25 4 H. LI Proof. Using the Cauchy-Schwarz inequality, 50), and Lemma 4.9, we have Id I 1P 1 )v E = a β Id I 1P 1 )v, Id I 1P 1 ) )v = a β Id I 1P 1 ) )v, v = Id I 1P 1 )v, A β v ) which completes the proof. Id I 1P 1 )v A β v CdimV β ) 1/ Id I 1P 1 )v E A β v, Then, we have the estimate for the error-propagation operator E in 47). Theorem Recall the number of smoothing steps l in Algorithm 4.1. Then, 66) a β E v, v ) C l + C a βv.v ), v V β, where C > 0 is independent of. Proof. Recall the space E defined in Lemma 4.9. By Lemma 4.10, 6), and Lemma 4.7, we first have a β Id P 1 )Rv l, Id P 1 )Rv l ) = Id P 1 )Rv l ) E Id R )Rv l, Rv l ) CdimV β ) 1 a β A β Rl v, Rv l ) = CdimV β ) 1 γ 1 a β 67) Ca β Id R )R l v, v ) C Id R l )v, v ). l a β We now give the proof by induction. For = 0, 66) holds since E 0 = 0. Assume 66) holds for 1, 1. Let c = Cl + C) 1 where C is the constant in 67). Note that c = 1 c)cl 1. Then, by 47), 49), the assumption on E 1, and 67), a β Id P 1 a β E v, v ) = a β Id P 1 )Rv l, Id P 1 )Rv l ) )R l v, Id P 1 +a β E 1 P 1 Rv l, P 1 Rv l ) )Rv l ) + caβ P 1 R l v, P 1 R l v ) = 1 c)a β Id P 1 )Rv l, I P 1 )Rv l ) + caβ R l v, R l v ) 1 c)cl 1 a β Id R l ) )v, v + caβ R l v, R l v ) which completes the proof. = ca β Id R l )v, v ) + caβ R l v, R l v ) = ca β v, v ), Hence, the multigrid V-cycle algorithm is a contraction with contraction number strictly less than 1, independent of the mesh level. Theorem 4.1. With the choice of the damping factor γ in 63), the multigrid V-cycle scheme defined in Algorithm 4.1 solving the axisymmetric equation 30) converges uniformly, E v E C l + C v E, v V β, where C is independent of the mesh level and E is the space in Lemma 4.9. Proof. By Lemma 4.8, let 0 α 1 α n be the eigenvalues of E and ν 1,, ν n be the corresponding orthonormal eigenvectors with respect to a β, ). Then, for any v V β, we write v = 1 i n ω i ν i. Therefore, a β E v, v ) =
REGULARITY AND MULTIGRID ANALYSIS FOR LAPLACE-TYPE AXISYMMETRIC EQUATIONS
MATHEMATICS OF COMPUTATION Volume 84, Number 93, May 015, Pages 1113 1144 S 005-5718014)0879- Article electronically published on September 8, 014 REGULARITY AND MULTIGRID ANALYSIS FOR LAPLACE-TYPE AXISYMMETRIC
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 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 informationANALYSIS OF A MODIFIED SCHRÖDINGER OPERATOR IN 2D: REGULARITY, INDEX, AND FEM
ANALYSIS OF A MODIFIED SCHRÖDINGER OPERATOR IN 2D: REGULARITY, INDEX, AND FEM HENGGUANG LI AND VICTOR NISTOR Version: 2.; Revised: April 2 by V.; Run: January 4, 2008 Abstract. Let r = (x 2 +x2 2 )/2 be
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 informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
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 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 informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationVARIATIONAL AND NON-VARIATIONAL MULTIGRID ALGORITHMS FOR THE LAPLACE-BELTRAMI OPERATOR.
VARIATIONAL AND NON-VARIATIONAL MULTIGRID ALGORITHMS FOR THE LAPLACE-BELTRAMI OPERATOR. ANDREA BONITO AND JOSEPH E. PASCIAK Abstract. We design and analyze variational and non-variational multigrid algorithms
More informationMULTIGRID METHODS FOR MAXWELL S EQUATIONS
MULTIGRID METHODS FOR MAXWELL S EQUATIONS A Dissertation Submitted to the Graduate Faculty of the Louisiana State University and Agricultural and Mechanical College in partial fulfillment of the requirements
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationweak Galerkin, finite element methods, interior estimates, second-order elliptic
INERIOR ENERGY ERROR ESIMAES FOR HE WEAK GALERKIN FINIE ELEMEN MEHOD HENGGUANG LI, LIN MU, AND XIU YE Abstract Consider the Poisson equation in a polytopal domain Ω R d (d = 2, 3) as the model problem
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationThomas Apel 1, Ariel L. Lombardi 2 and Max Winkler 1
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANISOTROPIC MESH REFINEMENT IN POLYHEDRAL DOMAINS: ERROR ESTIMATES WITH DATA IN
More informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
More informationINTRODUCTION TO MULTIGRID METHODS
INTRODUCTION TO MULTIGRID METHODS LONG CHEN 1. ALGEBRAIC EQUATION OF TWO POINT BOUNDARY VALUE PROBLEM We consider the discretization of Poisson equation in one dimension: (1) u = f, x (0, 1) u(0) = u(1)
More informationANALYSIS OF THE FINITE ELEMENT METHOD FOR TRANSMISSION/MIXED BOUNDARY VALUE PROBLEMS ON GENERAL POLYGONAL DOMAINS
ANALYSIS OF THE FINITE ELEMENT METHOD FOR TRANSMISSION/MIXED BOUNDARY VALUE PROBLEMS ON GENERAL POLYGONAL DOMAINS HENGGUANG LI, ANNA MAZZUCATO, AND VICTOR NISTOR Abstract. We study theoretical and practical
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
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 informationThe continuity method
The continuity method The method of continuity is used in conjunction with a priori estimates to prove the existence of suitably regular solutions to elliptic partial differential equations. One crucial
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 informationSubdiffusion in a nonconvex polygon
Subdiffusion in a nonconvex polygon Kim Ngan Le and William McLean The University of New South Wales Bishnu Lamichhane University of Newcastle Monash Workshop on Numerical PDEs, February 2016 Outline Time-fractional
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 informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationFinite Elements. Colin Cotter. January 18, Colin Cotter FEM
Finite Elements January 18, 2019 The finite element Given a triangulation T of a domain Ω, finite element spaces are defined according to 1. the form the functions take (usually polynomial) when restricted
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationFinite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem
Journal manuscript No. (will be inserted by the editor Finite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem Constantin Christof Christof Haubner Received:
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationIMA Preprint Series # 2114
IMPROVING THE RTE OF CONVERGENCE OF HIGH ORDER FINITE ELEMENTS ON POLYHEDR II: MESH REFINEMENTS ND INTERPOLTION By Constantin Bacuta Victor Nistor and Ludmil T. Zikatanov IM Preprint Series # 2114 ( pril
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationNEUMANN AND MIXED PROBLEMS ON CURVILINEAR POLYHEDRA
NEUMANN AND MIXED PROBLEMS ON CURVILINEAR POLYHEDRA Monique Dauge We prove regularity results in L p Sobolev spaces. On one hand, we state some abstract results by L p functional techniques : exponentially
More informationMultigrid finite element methods on semi-structured triangular grids
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, -5 septiembre 009 (pp. 8) Multigrid finite element methods on semi-structured triangular grids F.J.
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationR T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T
2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationQUASI-OPTIMAL RATES OF CONVERGENCE FOR THE GENERALIZED FINITE ELEMENT METHOD IN POLYGONAL DOMAINS. Anna L. Mazzucato, Victor Nistor, and Qingqin Qu
QUASI-OPTIMAL RATES OF CONVERGENCE FOR THE GENERALIZED FINITE ELEMENT METHOD IN POLYGONAL DOMAINS By Anna L. Mazzucato, Victor Nistor, and Qingqin Qu IMA Preprint Series #2408 (September 2012) INSTITUTE
More informationCONSTRAINED PERCOLATION ON Z 2
CONSTRAINED PERCOLATION ON Z 2 ZHONGYANG LI Abstract. We study a constrained percolation process on Z 2, and prove the almost sure nonexistence of infinite clusters and contours for a large class of probability
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
More informationFinite Element Methods for Maxwell Equations
CHAPTER 8 Finite Element Methods for Maxwell Equations The Maxwell equations comprise four first-order partial differential equations linking the fundamental electromagnetic quantities, the electric field
More 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] 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 informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 41-69, 2010. Copyright 2010,. ISSN 1068-9613. ETNA ANALYSIS OF THE FINITE ELEMENT METHOD FOR TRANSMISSION/MIXED BOUNDARY VALUE PROBLEMS ON
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationAnalysis of a high order trace finite element method for PDEs on level set surfaces
N O V E M B E R 2 0 1 6 P R E P R I N T 4 5 7 Analysis of a high order trace finite element method for PDEs on level set surfaces Jörg Grande *, Christoph Lehrenfeld and Arnold Reusken * Institut für Geometrie
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 informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationControllability of the linear 1D wave equation with inner moving for
Controllability of the linear D wave equation with inner moving forces ARNAUD MÜNCH Université Blaise Pascal - Clermont-Ferrand - France Toulouse, May 7, 4 joint work with CARLOS CASTRO (Madrid) and NICOLAE
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
More informationMATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5
MATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5.. The Arzela-Ascoli Theorem.. The Riemann mapping theorem Let X be a metric space, and let F be a family of continuous complex-valued functions on X. We have
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationarxiv: v2 [math.na] 17 Jun 2010
Numerische Mathematik manuscript No. (will be inserted by the editor) Local Multilevel Preconditioners for Elliptic Equations with Jump Coefficients on Bisection Grids Long Chen 1, Michael Holst 2, Jinchao
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationOversampling for the partition of unity parallel finite element algorithm
International conference of Computational Mathematics and its application Urumqi, China, July 25-27, 203 Oversampling for the partition of unity parallel finite element algorithm HAIBIAO ZHENG School of
More informationChapter One. The Calderón-Zygmund Theory I: Ellipticity
Chapter One The Calderón-Zygmund Theory I: Ellipticity Our story begins with a classical situation: convolution with homogeneous, Calderón- Zygmund ( kernels on R n. Let S n 1 R n denote the unit sphere
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More information(x, y) = d(x, y) = x y.
1 Euclidean geometry 1.1 Euclidean space Our story begins with a geometry which will be familiar to all readers, namely the geometry of Euclidean space. In this first chapter we study the Euclidean distance
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationA Mixed Method for Axisymmetric Div-Curl Systems
Portland State University PDXScholar Mathematics and Statistics Faculty Publications and Presentations Fariborz Maseeh Department of Mathematics and Statistics 2008 A Mixed Method for Axisymmetric Div-Curl
More informationHigh order Galerkin appoximations for parametric second order elliptic partial differential equations
Research Collection Report High order Galerkin appoximations for parametric second order elliptic partial differential equations Author(s): Nistor, Victor; Schwab, Christoph Publication Date: 2012 Permanent
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf-Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2018/19 Part 4: Iterative Methods PD
More informationIMPROVING THE RATE OF CONVERGENCE OF HIGH ORDER FINITE ELEMENTS ON POLYHEDRA II: MESH REFINEMENTS AND INTERPOLATION
IMPROVING THE RTE OF CONVERGENCE OF HIGH ORDER FINITE ELEMENTS ON POLYHEDR II: MESH REFINEMENTS ND INTERPOLTION CONSTNTIN BCUT, VICTOR NISTOR, ND LUDMIL T. ZIKTNOV Version: 3.4; Revised: 07/18/07; Run:
More informationNotes on Complex Analysis
Michael Papadimitrakis Notes on Complex Analysis Department of Mathematics University of Crete Contents The complex plane.. The complex plane...................................2 Argument and polar representation.........................
More informationWe simply compute: for v = x i e i, bilinearity of B implies that Q B (v) = B(v, v) is given by xi x j B(e i, e j ) =
Math 395. Quadratic spaces over R 1. Algebraic preliminaries Let V be a vector space over a field F. Recall that a quadratic form on V is a map Q : V F such that Q(cv) = c 2 Q(v) for all v V and c F, and
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationMultigrid Method ZHONG-CI SHI. Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China. Joint work: Xuejun Xu
Multigrid Method ZHONG-CI SHI Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China Joint work: Xuejun Xu Outline Introduction Standard Cascadic Multigrid method(cmg) Economical
More informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
More informationDiscrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem
The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant
More informationarxiv: v2 [math.ag] 24 Jun 2015
TRIANGULATIONS OF MONOTONE FAMILIES I: TWO-DIMENSIONAL FAMILIES arxiv:1402.0460v2 [math.ag] 24 Jun 2015 SAUGATA BASU, ANDREI GABRIELOV, AND NICOLAI VOROBJOV Abstract. Let K R n be a compact definable set
More informationW if p = 0; ; W ) if p 1. p times
Alternating and symmetric multilinear functions. Suppose and W are normed vector spaces. For each integer p we set {0} if p < 0; W if p = 0; ( ; W = L( }... {{... } ; W if p 1. p times We say µ p ( ; W
More informationNumerische Mathematik
Numer. Math. (2012) 122:61 99 DOI 10.1007/s00211-012-0456-x Numerische Mathematik C 0 elements for generalized indefinite Maxwell equations Huoyuan Duan Ping Lin Roger C. E. Tan Received: 31 July 2010
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationAnisotropic mesh refinement in polyhedral domains: error estimates with data in L 2 (Ω)
Anisotropic mesh refinement in polyhedral domains: error estimates with data in L 2 (Ω) Thomas Apel Ariel L. Lombardi Max Winkler February 6, 2014 arxiv:1303.2960v1 [math.na] 12 Mar 2013 Abstract. The
More informationAnalysis in weighted spaces : preliminary version
Analysis in weighted spaces : preliminary version Frank Pacard To cite this version: Frank Pacard. Analysis in weighted spaces : preliminary version. 3rd cycle. Téhéran (Iran, 2006, pp.75.
More informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
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 informationEQUIVALENCE OF TOPOLOGIES AND BOREL FIELDS FOR COUNTABLY-HILBERT SPACES
EQUIVALENCE OF TOPOLOGIES AND BOREL FIELDS FOR COUNTABLY-HILBERT SPACES JEREMY J. BECNEL Abstract. We examine the main topologies wea, strong, and inductive placed on the dual of a countably-normed space
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationIMA Preprint Series # 2143
A BASIC INEQUALITY FOR THE STOKES OPERATOR RELATED TO THE NAVIER BOUNDARY CONDITION By Luan Thach Hoang IMA Preprint Series # 2143 ( November 2006 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY
More informationPoint estimates for Green s matrix to boundary value problems for second order elliptic systems in a polyhedral cone
Maz ya, V. G., Roßmann, J.: Estimates for Green s matrix 1 ZAMM Z. angew. Math. Mech. 00 2004 0, 1 30 Maz ya, V. G.; Roßmann, J. Point estimates for Green s matrix to boundary value problems for second
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
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 informationSome notes on Coxeter groups
Some notes on Coxeter groups Brooks Roberts November 28, 2017 CONTENTS 1 Contents 1 Sources 2 2 Reflections 3 3 The orthogonal group 7 4 Finite subgroups in two dimensions 9 5 Finite subgroups in three
More informationA Perron-type theorem on the principal eigenvalue of nonsymmetric elliptic operators
A Perron-type theorem on the principal eigenvalue of nonsymmetric elliptic operators Lei Ni And I cherish more than anything else the Analogies, my most trustworthy masters. They know all the secrets of
More informationAnalysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient
Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient Yifeng Xu 1 Jun Zou 2 Abstract Based on a new a posteriori error estimator, an adaptive finite element method is proposed
More informationOutline of the course
School of Mathematical Sciences PURE MTH 3022 Geometry of Surfaces III, Semester 2, 20 Outline of the course Contents. Review 2. Differentiation in R n. 3 2.. Functions of class C k 4 2.2. Mean Value Theorem
More informationVISCOSITY SOLUTIONS. We follow Han and Lin, Elliptic Partial Differential Equations, 5.
VISCOSITY SOLUTIONS PETER HINTZ We follow Han and Lin, Elliptic Partial Differential Equations, 5. 1. Motivation Throughout, we will assume that Ω R n is a bounded and connected domain and that a ij C(Ω)
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 informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More information