WEIGHTED-NORM FIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PROBLEMS WITH CORNER SINGULARITIES
|
|
- Nickolas Jesse Terry
- 5 years ago
- Views:
Transcription
1 WEIGHTED-NORM FIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PROBLEMS WITH CORNER SINGULARITIES E. LEE, T. A. MANTEUFFEL, AND C. R. WESTPHAL Abstract. A weighted-norm least-squares method is considered for the numerical approximation of solutions that have singularities at the boundary. While many methods suffer from a global loss of accuracy due to boundary singularities, the least-squares method can be particularly sensitive to a loss of regularity. The method we describe here requires only a rough lower bound on the power of the singularity and can be applied to a wide range of elliptic equations. Optimal order discretization accuracy is achieved in weighted H, and functional norms and L accuracy is retained for boundary value problems with a dominant div/curl operator. Our analysis, including interpolation bounds and several Poincaré type inequalities, are carried out in appropriately weighted Sobolev spaces. Numerical results confirm the error bounds predicted in the analysis.. Introduction. In this paper, we develop a method for treating div/curl systems with reduced regularity. While motivated by first-order systems that arise from second-order elliptic boundary value problems, div/curl systems appear in many contexts, for example, in Maxwell s equations and in the vorticity from of Stokes equations. These problems have the fortunate property of a guaranteed smooth solution as long as the data and domain are smooth. However, many problems of interest are posed in non-smooth domains and, as a consequence, lose this property at a finite number of points on the boundary in two dimensions or along curves on the boundary in three dimensions. In the present paper, we study two-dimensional problems that have non-smooth solutions at irregular boundary points, that is, points that are corners of polygonal domains, locations of changing boundary condition type, or both. Similar behavior occurs in problems with discontinuous material coefficients and the methods presented here can easily be extended to that situation. Standard solution techniques that attempt to approximately solve a div/curl system with reduced regularity using H -conforming finite elements will, in general, fail to converge. This phenomenon can be explained by noting that the Sobolev space (H ) is a closed subspace of H(div) H(curl) (see section ), which implies either that (H ) = H(div) H(curl), the case of full regularity, or (H ) is a proper subspace. In this case, the co-dimension is finite and is spanned by so-called singular functions. In the presence of reduced regularity, the solution will, in general, not be in (H ). A standard finite element method using H -conforming elements will converge to the element of (H ) closest to the true solution in the H(div) H(curl) norm. Local mesh refinement will not alter this outcome. If a basis for the singular functions is known, it can be incorporated directly into the finite element space (c.f.[, 6, 4]). In [3, 4], this approach is applied to div/curl systems and shown to restore optimal convergence throughout the domain at a minimal additional cost. For some two-dimensional problems, the singular basis functions are known and can be included in the finite element space. For the other Department of Applied Mathematics, University of Colorado at Boulder; Boulder,CO , USA. {tmanteuf,eunjung}@colorado.edu. Department of Mathematics and Computer Science, Wabash College, PO Box 35, Crawfordsville, IN 47933, USA. westphac@wabash.edu. This work was sponsored by the National Science Foundation under grant number DMS 4873, and the Department of Energy under grant numbers DE FC-ER5479 and DE FG- 3ER5574. This work was sponsored by the National Science Foundation under grant DMS 9875
2 LEE, MANTEUFFEL, AND WESTPHAL problems, or in three dimensions, the exact character of the singular functions is less well understood which makes this approach more difficult to implement. As a different type of remedy for this so-called pollution effect, least-squares methods based on inverse norms can be effective for problems with irregular boundary points, discontinuous coefficients, and data in H. For example, in [5, 8, 3, 9], the functional is posed in terms of H -norms rather than L -norms, resulting in optimal L -approximations to the solution. A more recent approach, called FOSLL, uses an inverse norm induced by the equations, and is shown in [4, ] to be more efficient than the H -norm methods. Graded mesh refinement in weighted Sobolev spaces has been shown to be effective in restoring optimal convergence, in L - and H -norms, in the context of a Galerkin formulation of second-order elliptic problems with reduced regularity (c.f. [], []). However, in that context, the solution is in H. Convergence would occur, although more slowly, without weighting and mesh refinement. In [5], a weighted norm and a sequence of graded meshes is used in an H(div) least-squares functional, arising from a second-order elliptic problem, to restore optimal convergence for the primal variable, in both L - and H -norms, if the flux variable is approximated in a finite element space satisfying the grid decomposition property, for example, Raviart-Thomas elements (c.f.[6]). The flux variable converges in a weighted H(div)-norm. In this paper we examine div/curl systems that lack regularity within the firstorder system least-squares (FOSLS) framework. The basic FOSLS approach is to recast the original system as an appropriate first-order system and apply an L minimization principle over the residual of the equations. If possible, this reformulation is done by minimizing a functional whose quadratic part is equivalent to the product H -norm, indicating that the process is similar to solving a weakly coupled system of Poisson-like equations. This equivalence also guarantees optimal H -accuracy for standard discretizations. For div/curl systems with reduced regularity, as briefly mentioned above, the L based functional fails to be H -equivalent and, as a consequence, standard discretizations suffer from the pollution effect. Here, a weighted-norm leastsquares method is developed that restores optimal convergence using H -conforming finite elements without graded mesh refinement. It replaces the L -norms in the FOSLS functional with weighted L -norms, making the functional norm equivalent to a weighted H -norm. With an appropriate weighting function, this method recovers optimal order accuracy in the weighted L - and H -norms, and retains optimal L convergence even near the singularity. Our method requires only the power of the singularity (not the actual singular solution) to be known a priori and, in practice, can be used with only a rough estimate of the power of the singularity, which can be adaptively determined if unknown (c.f.[4]). The method developed here has some similarity to [5], but considers an H(div) H(curl) functional and considers more general boundary conditions. Most importantly, our analysis admits a more aggressive weighting, resulting in optimal order accuracy in the weighted norms without mesh refinement. In addition, we prove several Poincaré type inequalities in weighted Sobolev spaces under a variety of boundary conditions that, in addition to being necessary for our main result, may be of independent interest. We use Poisson s equation on a domain with a reentrant corner as a model problem and as the formal setting for analysis. The resulting div/curl system is the focus. The analysis in this paper is restricted to two dimensions. However, the approach suggests
3 FOSLS FOR PROBLEMS WITH SINGULARITIES 3 a natural generalization to problems with reduced regularity due to discontinuous coefficients and to problems in three dimensions. This paper is organized as follows. Section contains notation and preliminary discussion. In section 3, the weighted FOSLS functionals are described. Section 4 contains Poincaré inequalities and regularity results in weighted Sobolev spaces. Error bounds for the weighted FOSLS functional are presented in section 5 and section 6 contains computational results.. Singular solutions and preliminaries. For vector function u = (u, u ) t, let the divergence and curl of u be defined in the standard way: u = x u + y u and u = x u y u. Further, define the formal adjoint of the curl operator by ( ) y q q =. x q We use standard notation for Sobolev spaces H k (Ω) d, corresponding inner product (, ) k,ω, and norm k,ω, for k. We drop subscript Ω and superscript d when the domain and dimension are clear by context. Since H (Ω) coincides with L (Ω) we often denote by. Define the subspaces of L (Ω) induced by the divergence and curl of u by H(div) = {u L (Ω) : u < }, H(curl) = {u L (Ω) : u < }. We also make use of the following general inequalities for nonnegative a and b: a + b a + b ( a + b ). (.) Consider the function f(r, θ) = r a in two dimensional polar coordinates. Assume that the origin lies on the boundary of domain Ω w = {(r, θ) : < r < R, < θ < ω < π}, as pictured in figure.. By a direct computation it is clear that f H k (Ω) only for k < a +. ω R Fig... Simple wedge-shaped domain, Ω w. Now, consider Poisson s equation on a domain in R with a corner of interior angle ω. It is well known that, for the case of Dirichlet or Neumann boundary conditions, the solutions of this boundary value problem may include those with radial part of
4 4 LEE, MANTEUFFEL, AND WESTPHAL the form p r π ω in a local polar coordinate system centered at the corner. Thus, for the case of reentrant corners, ω > π, the solution fails to be in H (Ω) and we say that the problem has a singularity (or singular solution). For problems with Dirichlet and Neumann boundary conditions meeting at the corner, solutions may have components of the form p r π (ω). Thus, for mixed boundary conditions, singularities may occur at corners with ω > π/. We now explore this issue in more detail. Define the power of the singularity to be α = π/ω for Dirichlet or Neumann boundary conditions and α = π/(ω) for mixed boundary conditions. The solution to Poisson s equation may be written as p(r, θ) = p (r, θ) + s(r, θ), where p (r, θ) H (Ω) and s(r, θ) H +m (Ω) for m < α. The singular part of the solution has the form s(r, θ) = r α (κ sin(αθ) + κ cos(αθ)), where the values of κ and κ depend on boundary conditions (see [7, 8]). For the FOSLS formulation of this problem, we may similarly decompose unknown u = p as where s(r, θ) has the form u(r, θ) = u (r, θ) + s(r, θ), ( ) s(r, θ) = αr α κ sin(α )θ + κ cos(α )θ. κ cos(α )θ κ sin(α )θ Thus, the unknown u(r, θ) is in H k (Ω) only for k < α. For example, consider Poisson s equation posed on the simple domain in figure.. Let the solution to this boundary value problem in polar coordinates be p = χ(r)r 3 sin(θ/3), where χ(r) is a smooth transition function that is on a platform near the origin and vanishes at the boundaries not adjacent to the origin. Then, p = on Ω and p = r r(r r p) + r θθp ) = (r 3 χ 73 (r) + r 3 χ (r) sin(θ/3), and, thus, it is clear that p L (Ω), but p / H (Ω). We say this problem fails to provide full lifting of the data (from L (Ω) to H (Ω) e.g.). The solution, u = p, is, thus, not in H (Ω). 3. Weighted-norm least squares. As before, let Ω be a domain with a corner of interior angle ω at the origin, and we may, without loss of generality, further assume diam(ω). For f L (Ω), let p satisfy p = f, in Ω, p =, on Γ D, n p =, on Γ N, (3.)
5 FOSLS FOR PROBLEMS WITH SINGULARITIES 5 where n is the outward unit normal to Ω and Ω = Γ D Γ N. When this problem is H regular, the normal FOSLS methodology is to introduce the new unknown, u = p, and rewrite system (3.) as u p =, in Ω, u = f, in Ω, u =, in Ω, u =, on Γ D, n u =, on Γ N, p =, on Γ D, n p =, on Γ N. (3.) Here, is the counter-clockwise unit tangental vector to Ω. Since this system can be posed completely in terms of u, we may decouple the equations in (3.), solve for u first, and then recover p from u. To this end, define the two L -norm functionals, and the spaces, G(u; f) = u + f + u, G (p; u) = u p, V = {v H (Ω) : v = on Γ D, n v = on Γ N }, W = {q H (Ω) : q = on Γ D }. Thus, the two-stage solution process is to minimize G(v; f) over V and then, given the minimizer, u, minimize G over W: () G(u; f) = inf G(v; f), v V () G (p; u) = inf (q; u). q W (3.3) The goal of the FOSLS methodology is, generally, to formulate functionals whose quadratic part is equivalent to the H norm whenever possible. The second stage functional is H equivalent and the solution we seek is always in H. The first stage functional, however, is not always H equivalent. For domains with reentrant corners, there is no H sequence of functions that converges to the solution in the H(div) H(curl) norm. To illustrate, consider the example above where p = χ(r)r 3 sin(θ/3) and u = p. A simple computation reveals that u, u L (Ω), but u / H (Ω). Define the weighted functional by G w (u; f) = w( u + f) + w u, (3.4) where the weight function has the form w = r β for some β >. Define the weighted Sobolev norm, k,β, on Ω in terms of the standard L norm,, by q k,β = ( r β k+ j D j q ) /, (3.5) j k
6 6 LEE, MANTEUFFEL, AND WESTPHAL where D j is the standard distributional derivative of order j. Similarly, define the weighted seminorm by q k,β = ( r β k+ j D j q ) / (3.6) j =k and the associated weighted Sobolev space by Define the div/curl operator, L, and vector f by L = may now write the weighted functional from (3.4) as H k β(ω) = {q : q k,β < }. (3.7) G w (u; f) = Lu f,β. ( ) and f = ( ) f. We The weighted-norm least-squares minimization problem for the first-stage solution is then: find u V such that G w (u; f) = inf v V G w(v; f). The second-stage solution for p remains as described above. We seek values of β that make H (Ω) dense in H(div) H(curl) in the weighted functional norm and result in the most accurate discretizations possible. For the discrete problem, we may choose any finite dimensional subset of H over which to minimize the weighted functional. Let P h denote the space of C piecewise polynomial (or tensor product) elements on triangles (or quadrilaterals) of meshsize h and V h the subspace of P h that satisfies the appropriate boundary conditions on Ω: V h = {v h P h : v h = on Γ D, n v h = on Γ N }. The discrete weighted-norm least-squares minimization problem is, then, to minimize the discrete functional: find u h V h such that G w (u h ; f) = min v h V h G w (v h ; f). (3.8) By unweighting the equations near the singularity, the functional is freed from trying to approximate the solution (which is not in H (Ω)) in the H sense near the singularity. But, away from the singularity, the weighted functional retains the same character as the normal non-weighted functional. We now consider the choice of weight parameter β and its relation to weighted and nonweighted a priori error bounds on the approximated solution. 4. Poincaré Bounds and Regularity Estimates. In this section, we establish several theoretical results in weighted Sobolev spaces and error bounds for the weighted-norm method. Here, we establish several Poincaré bounds in the domain Ω w. We first prove a result for the scalar pure Neumann and pure Dirichlet problems, and then for the scalar mixed boundary condition problem. These results lead to a Poincaré inequality for the vector case.
7 FOSLS FOR PROBLEMS WITH SINGULARITIES 7 Lemma 4. Take Ω = Ω w and let β >, ɛ > and γ = β 3/ ɛ. Further, assume that γ. For functions q Hβ (Ω), with rβ ɛ q, r γ+ q L (Ω), that can be chosen to satisfy r γ q(r, θ) r dr dθ =, (4.) we have the bound: Ω q,β C(ɛ) q,β ɛ, (4.) where C(ɛ) depends on ɛ, β, and Ω and C(ɛ) as ɛ. Proof. For any two points (r, θ) and (r, θ ) in Ω, write q(r, θ) as r q(ˆr, θ) θ q(r, q(r, θ) = q(r, θ ) + dˆr + ˆθ) r ˆr θ ˆθ dˆθ. Multiplying both sides of the equation by r γ+, integrating with respect to r and θ over Ω, and Fubini s theorem yield R γ+ ω q(r, θ) γ + = ω R r r γ+ r q(ˆr, θ) ˆr dˆr dr dθ + ω R θ r γ+ θ q(r, ˆθ) ˆθ dˆθ dr dθ = ω + r ˆr R θ ˆθ r γ+ q(ˆr, θ) ˆr dr dˆr ω R R r γ+ q(r, ˆθ) ˆθ dθ dˆθ dr r r γ+ q(ˆr, θ) ˆr ˆr R ω ω θ ˆθ dr dˆr r γ+ q(r, ˆθ) ˆθ dθ dˆθ dr = ω γ + + R R ω γ+ q(ˆr, θ) ˆr ˆr ˆθ r γ+ q(r, ˆθ) ˆθ dˆr ωrγ+ γ + dˆθ dr ω R r R ω q(ˆr, θ) dˆr ˆr r γ+ q(r, ˆθ) ˆθ dˆθ dr. Note that above q is, by assumption, sufficiently smooth to apply Fubini s theorem. By the triangle inequality we have that R γ+ ω q(r, θ) γ + ω γ + + ω R R ω ˆr γ+ q(ˆr, θ) ωrγ+ ˆr dˆr + γ + r γ+ q(r, ˆθ) ˆθ dˆθ dr. R r q(ˆr, θ) ˆr dˆr
8 8 LEE, MANTEUFFEL, AND WESTPHAL Now, squaring each side and using (a + b + c) 3(a + b + c ) we get, ( q(r, θ) 3 R ( ˆr γ+ q(ˆr, θ) R R (γ+) ˆr dˆr) + 3 q(ˆr, θ) r ˆr dˆr ( (γ + R ω ) + r γ+ q(r, ˆθ) ) R (γ+) ˆθ dˆθ dr. ) (4.3) Multiply each side of (4.3) by r β and integrate with respect to r and θ over Ω. We consider each of the terms on the resulting right-hand side separately. First, ω R R r β ( R ω R R = Rβ+ β ω R Rβ+ β q,γ+. 3 ˆr γ+ q(ˆr, θ) ˆr dˆr) dr dθ r β ˆr γ+4 q(ˆr, θ) ˆr dˆr dr dθ ˆr γ+3 q(ˆr, θ) ˆr ˆr dˆr dθ We now consider the second term in (4.3). Since by the Schwarz inequality, ( R ) q(ˆr, θ) r ˆr dˆr ( R ) = ˆr (ɛ )/ˆr ( ɛ)/ q(ˆr, θ) r ˆr dˆr R R ˆr ɛ dˆr ˆr ɛ q(ˆr, θ) r r ˆr dˆr R = Rɛ ˆr ɛ q(ˆr, θ) ɛ ˆr dˆr, we can apply Fubini s theorem to bound the second term: ω R Rɛ ɛ = Rɛ ɛ = Rɛ ɛβ r ( R r β r ω R R r ω R ˆr ω R = Rɛ ɛβ q,β ɛ Rɛ ɛβ q,β ɛ, q(ˆr, θ) ˆr dˆr) dr dθ r β ˆr ɛ q(ˆr, θ) ˆr dˆr dr dθ r β ˆr ɛ q(ˆr, θ) ˆr dr dˆr dθ ˆr β ɛ q(ˆr, θ) ˆr ˆr dˆr dθ
9 FOSLS FOR PROBLEMS WITH SINGULARITIES 9 The third term can be bounded similarly: ω ( R R ω r β r γ+ q(r, ˆθ) ) ˆθ dˆθ dr dr dθ ( R ) ω ω R r β dr dθ Rω r γ+ q(r, ˆθ) ˆθ = Rβ+ ω ω R r γ+3 q(r, ˆθ) β r ˆθ r dr dˆθ = Rβ+ ω q β,γ+. 3 dr dˆθ Putting the three terms together and substituting γ = β 3/ ɛ we may now write the bound ω R r β q(r, θ) r dr dθ ( 3R ɛ β + 3Rɛ ɛβ + 6(β + ɛ) R ɛ ω ) q β,β ɛ, and the lemma follows by taking the square root of both sides. In what follows, let χ be a smooth function of r where χ = for r < η and χ = for r > η. We take η to be sufficiently small to ensure that supp(χ) Ω. Lemma 4. Take Ω = Ω w and let q be a scalar function in Hβ (Ω), where β >. The following bound holds for χ as defined above: χq,β β (χq),β. (4.4) Proof. (see []): Hardy s inequality for f(t) defined for t > with lim t f(t) = gives f t dt 4 f dt. (4.5) The lemma follows after a change of variables, t = r β, a substitution f(r) = χq(r, θ) for fixed θ, and an integration on both sides with respect to θ. Lemma 4.3 Take Ω = Ω w and let either q Hβ (Ω) with q = on Ω or q (Ω)/R with n q = on Ω. Then H β for β > where C depends only on Ω, and β. Proof. First, if q = on Ω, write q = q(r, θ) as q,β C q,β (4.6) q(r, θ) = θ q(r, ˆθ) ˆθ dˆθ.
10 LEE, MANTEUFFEL, AND WESTPHAL Square both sides and multiply by r β : θ r β q(r, θ) = r β q(r, ˆθ) ˆθ dˆθ ω r β+ q(r, ω ˆθ) r ˆθ Integrate both sides with respect to r and θ over Ω: ω R ω R ω r β q(r, θ) r dr dθ ω r β+ q(r, ˆθ) r ˆθ dˆθ dr dθ ω R ω r β q(r, ˆθ) r ˆθ r dr dˆθ. The lemma follows since the right side is bounded by C q,β. Now, if q Hβ (Ω)/R then it may be chosen to satisfy r γ q r dr dθ = (4.7) Ω for γ chosen as in lemma 4.. By the triangle inequality and lemma 4. we get q,β χ q,β + ( χ)q,β ( ω R ( ) r C (χq),β + r β ( χ) q r dr dθ) η C ( (q),β + q,β ). η Apply lemma 4. with ɛ = to the q,β term on the right side and the lemma follows. For the problem with mixed boundary conditions, consider Ω w partitioned into subdomains Ω = {(r, θ) : r R, θ ω} and Ω = Ω\Ω, as shown in figure 4.. dˆθ. Ω Ω R R R Fig. 4.. Wedge-shaped domain, Ω w, partitioned into subdomains Ω and Ω.
11 FOSLS FOR PROBLEMS WITH SINGULARITIES Lemma 4.4 Consider domain Ω = Ω w, as pictured in figure 4., and let q H β (Ω) vanish on the line segment of Ω corresponding to θ = and r < R. Then, there is a constant, C, dependent only on Ω, β, and R, such that for β >. q,β C q,β (4.8) Proof. For points (r, θ) in Ω we may derive the bound, q,β,ω C q,β,ω, (4.9) completely analogous to the proof of lemma 4.3. Now, consider points (r, θ) in Ω. We may write q = q(r, θ) as q(r, θ) = r r q(ˆr, θ) ˆr θ q( r, dˆr + ˆθ) ˆθ dˆθ, where the point ( r, ) is on the part of Ω where q vanishes. By the Schwarz inequality, the triangle inequality and inequality (.) we have the bound q(r, θ) (R R ) r r q(ˆr, θ) ˆr dˆr + ω θ q( r, ˆθ) ˆθ dˆθ. (4.) We now expand the limits in the integrals, multiply each side by r β, integrate with respect to r over ( R, R), integrate with respect to θ over (, ω), and integrate with respect to r over ( R, R ), by Fubini s theorem: ω R ( R ) r β q(r, θ) dr dθ R (R R ω R R R ) r β R R R R ω R θ q( r, + ω ˆθ) ˆθ dˆθ dr dθ d r. R R q(ˆr, θ) ˆr dˆr dr dθ d r (4.) We use the inequalities, R r R, r ˆr r R, to derive the following simple bounds: ( ˆr R ), r ( R R ) ˆr, ( ) R r r. (4.) R
12 LEE, MANTEUFFEL, AND WESTPHAL By applying the bounds in (4.) and Fubini s theorem, we may now write (4.) as ω R ( R ) r β q(r, θ) r dr dθ R (R ω R ( ) ( ) β ˆr R R ) R ˆr β q(ˆr, θ) R R R ˆr dˆr dθ + ω (R R ω ( ) β R R ) r β q( r, ˆθ) R R ˆθ dˆθ d r β+ (R ω R R ) R β R β ˆr β q(ˆr, θ) ˆr ˆr dˆr dθ R + β ω (R ( ) β R R ω R ) r β q( r, ˆθ) R R r ˆθ r dˆθ d r, which directly implies q,β,ω C q,β,ω, (4.3) where C = β+ (R R )R β R β ((R R ) + ω R ). Combining inequalities (4.9) and (4.3) completes the lemma. We now consider a similar Poincaré inequality for the vector case. Again, consider Ω = Ω w, where Ω is partitioned into Dirichlet and Neumann boundaries, Γ D and Γ N respectively. The following lemma is valid for the pure Dirichlet and Neumann cases and for the mixed boundary condition cases when Γ D includes a part of the boundary adjacent to the origin and ω 3π. Lemma 4.5 Take Ω = Ω w and let u H β (Ω) satisfy u = on Γ D and n u = on Γ N. Assume that for the mixed boundary condition case that ω 3π/. Then there is a constant, C, dependent only on Ω, β and the length of the segments of Γ D and Γ N adjacent to the origin, such that for β >. u,β C u,β (4.4) Proof. First, consider the case when u = on Ω. Denote the part of Ω aligned with θ = as Γ and the part of Ω aligned with θ = ω as Γ. Thus, u = on Γ and x u + y u = on Γ. Since u and x u + y u satisfy the conditions in lemma 4.4, we may use and u,β C u,β x u + y u,β C ( x u + y u ),β.
13 FOSLS FOR PROBLEMS WITH SINGULARITIES 3 Further, take y, since y = corresponds to either ω = π, for which the result holds trivially since the boundary is smooth, or ω = π, which we do not consider. Now, u,β = u,β + u,β = u,β + y x u x u + y u,β ( + x y ) u,β + y x u + y u,β C( u,β + ( x u + y u ),β) C u,β. The case when n u = on Ω is analogous since u = on Γ and n x u +n y u = on Γ. Also, when ω π, 3π, the case for mixed boundary conditions follows similarly using the result of lemma 4.4. The case for mixed boundary conditions when ω = π/, follows from appealing to symmetry in the the pure Dirichlet problem for ω = π. Remark 4.6 Lemma 4.5 can be directly extended to more generally shaped domains. The proof of the scalar Poincaré bounds in lemmas 4., 4., 4.3 and 4.4 are simplified when the domain has the shape of Ω w with only one irregular boundary point. Since we are primarily interested in a local result, proving lemma 4.5 in the simple domain is sufficient for our purposes. Consider the following scalar Poisson problem in Ω w : p = f, in Ω, p =, on Γ D, n p =, on Γ N. (4.5) We refer to system (4.5) as the pure Dirichlet problem when Ω = Γ D ; the pure Neumann problem when Ω = Γ N ; and the mixed boundary condition problem when Γ D includes the part of Ω coinciding with one of either θ = or θ = ω, and Γ N = Ω\Γ D with Γ N. The following regularity results can be found in [3] and []. Lemma 4.7 Assume β < π/ω for the pure Dirichlet problem, < β < π/ω for the pure Neumann problem and β < π/ω for the mixed boundary condition problem. Then, for every f Hβ (Ω), there exists a unique solution to (4.5), p Hβ (Ω) for the pure Dirichlet and mixed boundary condition cases and p H β (Ω)/R for the pure Neumann problem. Moreover, there exists a constant, C, independent of p, such that p,β C f,β. (4.6) Proof. See Chapter of [3] for the Dirichlet and Neumann problems and Chapter of [] for the mixed boundary problem. Define the subspace of functions in Hβ (Ω) satisfying the appropriate boundary conditions by V β = {v H β(ω) : v = on Γ D, n v = on Γ N }.
14 4 LEE, MANTEUFFEL, AND WESTPHAL We now prove a regularity result for functions in V β. Recall that the power of the singularity is defined as α = π/ω for Dirichlet or Neumann boundary conditions and α = π/(ω) for mixed boundary conditions. Lemma 4.8 Consider domain Ω = Ω w. Then there is a positive constant, C, independent of u, such that, for β < α, the following bound holds for all u V β : u,β C Lu,β. (4.7) Proof. From lemma 4.7 we know that any u V β has the decomposition u = φ + ψ, (4.8) where φ, ψ H β (Ω) satisfy { φ = u, in Ω, φ =, on Ω, (4.9) and { ψ = u, in Ω, n ψ =, on Ω. (4.) Then, by applying lemma 4.7 to problems (4.9) and (4.) we have u,β = φ + ψ,β φ,β + ψ,β φ,β + ψ,β C( u,β + u,β ) C Lu,β, which completes the proof. 5. Error Bounds. Let T h = N i= i be a quasi-uniform triangulation of polygonal domain Ω. Let I h represent standard interpolation onto a piecewise polynomial finite element space of degree k. From finite element theory, we have the following interpolation bounds. Lemma 5. Let Ω be a polygonal domain. There exists a constant, C, independent of v, such that, for all v H m (Ω), / v I h v s, Ch m s v m (5.) T h for s m and < m. Here, I h denotes interpolation by a piecewise polynomial of degree k = m. (Note: here the norm s, is the standard H s () norm.) Proof. See [9] or [7]. We now consider a weighted interpolation bound for functions on domains with a polygonal corner at the origin. Define the modified interpolation operator, I h, by { I h I u = h u = n i= u(a i)φ i, if does not intersect the origin, n i= u(a i)φ i, if intersects the origin, where I h is a standard polynomial interpolation operator, φ i are basis functions corresponding to the n + nodal points, a i, and a is the origin, (, ). Thus, the modified interpolation has a value of zero at the origin and resembles I h away from the origin.
15 FOSLS FOR PROBLEMS WITH SINGULARITIES 5 Lemma 5. Let Ω be a polygonal domain. There exists a constant, C, independent of u, such that, for all u Hβ m (Ω) satisfying equation (4.8), u I h u,β, T h / Ch m u m,β, (5.) for < m and β >, where I h is the modified interpolation operator onto piecewise polynomials of degree k = m defined above. Proof. Define K = { (, ) } as the set of elements adjacent to the origin. On T h \K, we have h r min r = x + y r max r min + h with r min = inf{r (x, y) } and r max = sup{r (x, y) } in, and u I h u,β, = u I h u,β, = r β (u I h u) + r (β ) u I h u d rmax β (u I h u) d + rmaxr β min u I h u d Cr β maxh (m ) u m,, + Cr β maxr min hm u m,, = Crmaxh β (m ) ( + r min h ) u m,, Cr β Ch (m ) rmaxr β β min r β D m u d ( Ch (m ) r min + ) β h r β D m u d r min Ch (m ) r β D m u d. maxh (m ) u m,, We now consider the case for which K. Let δ C be a cut-off function defined by δ(r) = {, if r h/3,, if r > h/3, with δ (m) ch m, where δ (m) is the m th derivative of δ. By the triangle inequality, u I h u,β, δu I h (δu),β, + ( δ)u I h (( δ)u),β,. (5.3) By the definition of δ we have I h (( δ)u) = I h (( δ)u) and I h (δu) =. For the second term in (5.3), we apply lemmas 4.4 and 5., the properties in δ, and Fubini s
16 6 LEE, MANTEUFFEL, AND WESTPHAL theorem to obtain ( δ)u I h (( δ)u),β, = ( δ)u I h (( δ)u),β, C (( δ)u I h (( δ)u)),β, Ch β (( δ)u I h (( δ)u)) d Ch (β+m ) D m (( δ)u) d Ch (β+m ) D m j ( δ)d j u d Ch (β+m ) Ch (β+m ) = Ch (m ) j m h 3 h 3 j m j m h 3 h 3 h j m D j u d + j m r(θ) h 3 h β r (β+ j m) D j u d + ( δ)d m u d r(θ) r (β+ j m) D j u d = Ch (m ) u m,β,. h 3 (r/h) β D m u d Using the properties of δ and Fubini s theorem result in a similar bound for the first term in (5.3): δu I h (δu),β, = δu,β, = r β (δu) + r (β ) δu d C r β ( δ u + δ u ) + r (β ) δu d C C C h 3 h 3 h 3 h 3 h r β h u 3 d + C r β u + r (β ) u d h r (β ) u 3 d + C r β u + r (β ) u d r (m ) (r (β m+) u + r (β m) u )d Ch (m ) u m,β,. Thus we have u I h u,β, Ch (m ) T h and the lemma follows. T h u m,β, Ch (m ) u m,β, Lemma 5.3 Assume equation (4.4) holds in Ω. Then, for all u h V h, u h,β Ch η u h,β+η (5.4) for β > and η >. Proof. Using lemma 4.5 and an inverse inequality, we may write u h,β C u h,β+ Ch u h,β+,
17 FOSLS FOR PROBLEMS WITH SINGULARITIES 7 which establishes (5.4) for η =. Repeated application of this inequality thus validates (5.4) for any positive integer. Now consider u h,β = r β u h, r β u h = r β / u h, r β+/ u h u h,β / uh,β+ / Ch u h,β+ /. Taking the square root establishes (5.4) for η = /. Repeating these steps leads to (5.4) for all η = η n = k n / ln for any nonnegative integers k n and l n. For any η >, choose a monotonically decreasing sequence, {η n }, such that η n > η and lim n η n η =. Now, g n = (r β+ηn u h ) is a monotonically increasing function that converges to g = (r β+η u h ) pointwise everywhere. Thus, by the Lebesgue monotone convergence theorem, we have u h,β+η = gdx = lim g n dx = lim u h,β+η n n n and, therefore, u h,β = lim n u h,β lim n Ch ηn u h,β+ηn = (lim n Ch ηn )(lim n u h,β+ηn ) = Ch η u h,β+η, which completes the proof. Define an irregular boundary point of polygonal domain Ω to be a point on Ω where interior angle ω satisfies ω > π when Dirichlet or Neumann boundary conditions are applied on both sides of the point or ω > π/ when one Dirichlet boundary and a Neumann boundary meet at the corner. We now present error bounds for the numerical solution in weighted and unweighted norms. Theorem 5.4 Let Ω be a polygonal domain with one irregular boundary point of interior angle ω and let f L (Ω). Suppose u V satisfy Lu = f. If u h V h is chosen to minimize the weighted functional, G w (u h ; f) = Lu h f,β = inf v h V h Lv h f,β, for β < α, then the approximation error, u u h, satisfies the following bounds: u u h,β Ch α+β u α+β,β, (5.5) G w (u u h ; ) / Ch α+β u α+β,β, (5.6) u u h,β Ch s+β u α+β,β, (5.7) u u h Ch s u α+β,β, (5.8) where s < α, for α + β k + with k the degree of the piecewise polynomial elements in V h.
18 8 LEE, MANTEUFFEL, AND WESTPHAL Proof. By lemmas 4.8 and 5., we have u u h,β C L(u u h ),β C L(u I h u),β C u I h u,β Ch α+β u α+β,β, which establishes both (5.5) and (5.6) since we may write L(u u h ),β = G w (u u h ; ) /. Note that lemmas 4.8 and 5. are satisfied for β < α and α + β. For the weighted L -norm, we write u u h,β = T h u u h,β,. The Cauchy inequality yields, for any ɛ (, ), ( ) ɛ ( u u h,β, = r β u u h d d Ch ɛ ( ( r β u u h ) ) ɛ ɛ d ( r β u u h ) ) ɛ ɛ d = Ch ɛ r β (u u h ). L ɛ () Since H (Ω) is continuously imbedded into L q (Ω) for all q [, ), u u h,β = u u h,β, Ch ɛ T h Thus, by (5.5), we have T h r β (u u h ) L ɛ () Ch ɛ r β (u u h ) L ɛ (Ω) Ch ɛ r β (u u h ) H (Ω) Ch ɛ u u h,β. u u h,β Ch ɛ u u h,β Ch s+β u α+β,β, where any s < α. We now consider the bound on u u h. Let K = { (, ) } and K = T h \K. First, we consider the case β <. If K, then r Ch and r β Ch β. Thus, we have u u h, Ch ( β) u u h,β, Ch ( β) u u h,β,. (5.9) If K, we use the technique above to get u u h, Ch ( ɛ) u u h L ɛ (). Again, since H is continuously imbedded into L q for all q [, ), we have u u h, Ch ɛ u u h L ɛ () K K Ch ɛ u u h L ɛ (K ) Ch ɛ u u h H (K ) = Ch ɛ K r β r β ( u u h + (u u h ) ) dω Ch ( β ɛ) u u h,β,k. (5.)
19 FOSLS FOR PROBLEMS WITH SINGULARITIES 9 Hence by (5.9), (5.), and (5.5) we have u u h = K u u h, + Ch s u α+β,β, where s < α. The proof for β follows analogously. K u u h, h ( β ɛ) u u h,β Remark 5.5 For the optimal finite element convergence of O(h) with respect to the weighted functional and H norms, we select β = α. But theorem 5.4 also requires that β < + α. Thus, when α [/, ), we may use a weighting with β = α and expect optimal rates, but when α (, /), our theory only guarantees at best O(α) convergence using β = + α. Numerical results, however, indicate that values of β larger than the theory allows can be used to recover optimal rates. We explore this in the next section. 6. Computational results. In this section, we present some numerical examples of the weighted-norm procedure to validate the error bounds in the previous section. As a test problem, we minimize the weighted functional on the following L-shaped domain: Ω = (.5,.5) \ [,.5) (.5, ], which yields α = π/ω = /3. Function f is chosen so that the solution of this test problem is u = (χ(r)r 3 sin(θ/3)), where χ(r) = for r < /8, χ(r) = for r > 3/8, and χ(r) is C smooth. Again, note that f L (Ω) but u / H (Ω). Define the following measures of the accuracy of the computed solution, u h : nonweighted functional norm G / = ( u h f + u h ) /, nonweighted L norm of the error ɛ = u u h, nonweighted H seminorm of the error ɛ = u u h, weighted functional norm G / w = G w (u h ; f) /, weighted L norm of the error ɛ w = u u h,β, weighted H seminorm of the error ɛ w = u u h,β. Since α = /3, we choose the optimal weight parameter, β = α = 4/3, for our computations. Table 6. summarizes discretization error and convergence rates for β = 4/3. Asymptotic convergence rates in Ω are found to be approximately O(h) for G / w and ɛ w, O(h ) for ɛ w and O(h 3 ) for ɛ. The approximation does not converge in either the ɛ or G / measures since u / H (Ω). To distinguish between behavior near to and away from the singularity, we consider the error of the solution above on a partitioning of Ω. Define Ω = Ω ( 3 8, 5 8 ) and Ω = Ω\Ω ; see figure 6.. Table 6. summarizes the asymptotic discretization accuracy obtained at the finest mesh size in subdomains Ω and Ω. Away from the singularity we observe optimal accuracy in all measures. As expected, near the singularity, the solution fails to converge in the nonweighted functional and H norms. The nonweighted L error achieves accuracy of approximately O(h 3 ) near the singularity.
20 LEE, MANTEUFFEL, AND WESTPHAL h G / w Ratio Rate ɛ w Ratio Rate e e e e e e e e-.4.3 h ɛ w Ratio Rate ɛ Ratio Rate 8 3.8e- 3.7e- 6.35e e e e e e e e e e e e Table 6. Convergence of discretization error for weighted-norm FOSLS. Ω Ω Fig. 6.. L-shaped domain Ω and subdomains Ω and Ω. G / w G / ɛ w ɛ ɛ w ɛ Ω O(h) O(h) O(h) O(h) O(h ) O(h ) Ω O(h) O() O(h) O() O(h ) O(h 3 ) Ω O(h) O() O(h) O() O(h ) O(h 3 ) Table 6. Accuracy in Ω, Ω, and Ω with β = α. Figure 6. shows the first component of the exact solution, u, and the standard FOSLS approximation u h. Figure 6.3 shows the error of the first component of the approximated solution for the standard FOSLS and the weighted-norm FOSLS methods. We see that the error in the approximation in standard FOSLS is highest near the
21 FOSLS FOR PROBLEMS WITH SINGULARITIES (a) Exact component u. computed by stan- (b) Component u h dard FOSLS. Fig. 6.. Exact solution component u and solution component u h FOSLS on the h = 3 mesh. approximated by standard L error in u L error in u Sfrag replacements..5. PSfrag replacements (a) Standard FOSLS. (b) Weighted-norm FOSLS, β = α. Fig Reduction of the pollution effect by the weighted-norm procedure. Each plot is the error of solution component u h on the h = 3 mesh. singularity, but remains large even away from the corner point. In the weighted-norm FOSLS implementation, the error remains large near the singularity, as we expect, but is now concentrated only near the corner point. The pollution effect is removed by the weighting procedure. There are many boundary value problems not directly covered by the theory presented here that are of interest. For example, Poisson s equation with mixed boundary conditions on the domain used above has a value of α = /3. To recover optimal convergence for this problem, the weighted-norm method requires a value of β larger than theorem 5.4 allows. In other elliptic equations (e.g., Stokes or the linear elasticity equations), the value of α is generally smaller than for Poisson s equation for the same domain and boundary condition type. In each of these cases, a larger β
22 LEE, MANTEUFFEL, AND WESTPHAL value is necessary for optimal convergence. This leads us to consider using larger β than the theory allows. Consider the same example problem as above on uniform mesh sizes of h = /8, /6,..., /5, and values of β ranging from /3 to 3/6. Figure 6.4 plots the convergence rate at the finest level for: the weighted functional norm, G / w ; the weighted L norm, ɛ w; and the L norm, ɛ. While the functional Convergence Rates vs. Β weighted L error weighted functional norm.67 L error Β Fig Convergence rates versus β. The shaded region indicates values of β for which the assumptions of theorem 5.4 are satisfied. norm retains optimal accuracy for large values of β, the solution fails to converge in the weighted and nonweighted L measures for β 3. This indicates that, although the weighted-norm approach seems to be more robust than the theory allows, large values of β should still be used with caution. The method presented here is applicable to a wide range of problems and provides an efficient alternative to more specialized techniques for treating singularities in boundary value problems. Further numerical results for other problems including systems and in three dimensions can be seen in a companion paper. REFERENCES [] T. Apel and B. Heinrich, Mesh refinement and windowing near edges for some elliptic problem, SIAM J. Numer. Anal., 3 (994), pp [] I. Babuska, R. Kellogg, and J. Pitkaranta, Direct and inverse error estimates for finite elements with mesh refinements, Numer. Math., 33 (979), pp [3] M. Berndt, T. Manteuffel, S. McCormick, and G. Starke, Analysis of first-order system least squares (fosls) for elliptic problems with discontinuous coefficients: part i, SIAM J. Numer. Anal., 43 (5), pp [4], Analysis of first-order system least squares (fosls) for elliptic problems with discontinuous coefficients: part ii, SIAM J. Numer. Anal., 43 (5), pp [5] P. Bochev and M. Gunzburger, Analysis of least-squares finite element methods for the stokes equations, Math. Comp., 63 (994), pp [6], On least-squares finite element methods for the poisson equation and their connection to the dirichlet and kelvin principles, SIAM J. Numer. Anal., 43 (5), pp [7] D. Braess, Finite Elements: Theory, Fast Solvers and Applications in Solid Mechanics, Cambridge,. [8] J. Bramble, R. Lazarov, and J. Pasciak, A least-squares approach based on a discrete minus one inner product for first order systems, Math. Comp., 66 (997), pp
23 FOSLS FOR PROBLEMS WITH SINGULARITIES 3 [9] S. Brenner and L. Scott, The Mathematical Theory of Finite Element Methods, Springer- Verlag, 994. [] Z. Cai and S. Kim, A finite element method using singular functions for the poisson equation: Corner singularities, SIAM J. Numer. Anal., 39 (), pp [] Z. Cai, T. Manteuffel, and S. McCormick, First-order system least squares for the stokes equations, with application to linear elasticity, SIAM J. Numer. Anal., 34 (997), pp [], First-order system least squares for velocity-vorticity-pressure form of the stokes equations, with application to linear elasticity, Elect. Trans. Numer. Anal., 3 (997), pp [3] Z. Cai, T. Manteuffel, S. McCormick, and S. Parter, First-order system least squares (FOSLS) for planar linear elasticity: Pure traction problem, SIAM J. Numer. Anal., 35 (998), pp [4] Z. Cai, T. Manteuffel, S. McCormick, and J. Ruge, First-order system LL (FOSLL*): Scalar elliptic partial differential equations, SIAM J. Numer. Anal., 39 (), pp [5] C. Cox and G. Fix, On the accuracy of least squares methods in the presence of corner singularities, Comp. and Maths. with Appls., (984), pp [6] G. Fix, S. Gulati, and G. Wakoff, On the use of singular functions with finite element approximations, J. Comput. Phys., 3 (973), pp [7] P. Grisvard, Elliptic Problems in Nonsmooth Domains, Pitman, 985. [8], Singularities in Boundary Value Problems, Springer-Verlag, 99. [9] S. Kim, T. Manteuffel, and S. McCormick, First-order system least squares (fosls) for spatial linear elasticity: Pure traction, SIAM J. Numer. Anal., 38 (), pp [] V. Kondratiev, Boundary problems for elliptic equations in domains with conical or angular points, Trans. Moscow Math. Soc., 6 (967), pp [] V. Kozlov, V. Maz ya, and J. Rossmann, Spectral Problems Associated with Corner Singularities of Solutions to Elliptic Equations, vol. 85, American Mathematical Society,. [] T. Manteuffel, S. McCormick, J. Ruge, and J. Schmidt, First-order system ll* (fosll*) for general scalar elliptic problems in the plane, SIAM J. Numer. Anal., 43 (5), pp. 98. [3] V. Maz ya, S. Nazarov, and B. Plamenevskij, Asymptotic Theory of Elliptic Boundary Vaule Problems in Singularly Perturbed Domains, Volume I, vol., Birkhäuser,. Operator Theory Advances and Applications. [4] G. Strang and G. Fix, An analysis of the finite element method, Prentice-Hall, Englewood Cliffs, NJ, 973.
A 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 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 informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More 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 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 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 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 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 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 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 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 information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
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 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 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 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 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 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 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 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 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 informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationEfficiency-based h- and hp-refinement strategies for finite element methods
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 7; : 5 [Version: /9/8 v.] Efficiency-based h- and hp-refinement strategies for finite element methods H. De Sterck, Thomas A. Manteuffel,
More informationAn a posteriori error estimator for the weak Galerkin least-squares finite-element method
An a posteriori error estimator for the weak Galerkin least-squares finite-element method James H. Adler a, Xiaozhe Hu a, Lin Mu b, Xiu Ye c a Department of Mathematics, Tufts University, Medford, MA 02155
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 informationFIRST-ORDER SYSTEM LEAST SQUARES FOR THE STOKES EQUATIONS, WITH APPLICATION TO LINEAR ELASTICITY
SIAM J. NUMER. ANAL. c 1997 Society for Industrial Applied Mathematics Vol. 34, No. 5, pp. 1727 1741, October 1997 004 FIRST-ORDER SYSTEM LEAST SQUARES FOR THE STOKES EQUATIONS, WITH APPLICATION TO LINEAR
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 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 informationJournal of Non-Newtonian Fluid Mechanics
J. Non-Newtonian Fluid Mech. 159 (2009) 72 80 Contents lists available at ScienceDirect Journal of Non-Newtonian Fluid Mechanics journal homepage: www.elsevier.com/locate/jnnfm An adaptive mixed least-squares
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 informationFIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR GEOMETRICALLY NONLINEAR ELASTICITY
FIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR GEOMETRICALLY NONLINEAR ELASTICITY T. A. MANTEUFFEL, S. F. MCCORMICK, J. G. SCHMIDT, AND C. R. WESTPHAL Abstract. We present a first-order system least-squares
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 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 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 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 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 informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationFINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS
Numerical Functional Analysis and Optimization, 28(7 8):957 973, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0163-0563 print/1532-2467 online DOI: 10.1080/01630560701493305 FINITE ELEMENT APPROXIMATION
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 informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More 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 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 informationSpectral element agglomerate AMGe
Spectral element agglomerate AMGe T. Chartier 1, R. Falgout 2, V. E. Henson 2, J. E. Jones 4, T. A. Manteuffel 3, S. F. McCormick 3, J. W. Ruge 3, and P. S. Vassilevski 2 1 Department of Mathematics, Davidson
More informationA posteriori error estimation of approximate boundary fluxes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2008; 24:421 434 Published online 24 May 2007 in Wiley InterScience (www.interscience.wiley.com)..1014 A posteriori error estimation
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
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 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 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 informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
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 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 informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
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 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 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 informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
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 informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More information(bu) = f in Ω, (1.1) u = g on Γ I, (1.2)
A DUAL LEAST-SQUARES FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PDES: A NUMERICAL STUDY LUKE OLSON Abstract. In this paper, we develop a least-squares finite element method for linear Partial Differential
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More 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 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 informationComparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes
Comparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes Do Y. Kwak, 1 JunS.Lee 1 Department of Mathematics, KAIST, Taejon 305-701, Korea Department of Mathematics,
More informationA Posteriori Error Estimation Techniques for Finite Element Methods. Zhiqiang Cai Purdue University
A Posteriori Error Estimation Techniques for Finite Element Methods Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 16, 2017 Books Ainsworth & Oden, A posteriori
More informationMath Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space
Math Tune-Up Louisiana State University August, 2008 Lectures on Partial Differential Equations and Hilbert Space 1. A linear partial differential equation of physics We begin by considering the simplest
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
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 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 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 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 informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
More informationA MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS
SIAM J. SCI. COMPUT. Vol. 29, No. 5, pp. 2078 2095 c 2007 Society for Industrial and Applied Mathematics A MULTIGRID METHOD FOR THE PSEUDOSTRESS FORMULATION OF STOKES PROBLEMS ZHIQIANG CAI AND YANQIU WANG
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 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 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 informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
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 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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationINNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD
INNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD Abstract. Finite element methods for the Reissner Mindlin plate theory are discussed. Methods in which both the tranverse displacement and
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 informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationNECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES
NECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES JUHA LEHRBÄCK Abstract. We establish necessary conditions for domains Ω R n which admit the pointwise (p, β)-hardy inequality u(x) Cd Ω(x)
More informationA LEAST-SQUARES APPROACH BASED ON A DISCRETE MINUS ONE INNER PRODUCT FOR FIRST ORDER SYSTEMS
MATHEMATICS OF COMPUTATION Volume 66, Number 219, July 1997, Pages 935 955 S 0025-5718(97)00848-X A LEAST-SQUARES APPROACH BASED ON A DISCRETE MINUS ONE INNER PRODUCT FOR FIRST ORDER SYSTEMS JAMES H BRAMBLE,
More informationEigenvalues of Robin Laplacians on infinite sectors and application to polygons
Eigenvalues of Robin Laplacians on infinite sectors and application to polygons Magda Khalile (joint work with Konstantin Pankrashkin) Université Paris Sud /25 Robin Laplacians on infinite sectors 1 /
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
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 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 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 informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
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 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 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 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 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 informationA first-order system Petrov-Galerkin discretisation for a reaction-diffusion problem on a fitted mesh
A first-order system Petrov-Galerkin discretisation for a reaction-diffusion problem on a fitted mesh James Adler, Department of Mathematics Tufts University Medford, MA 02155 Scott MacLachlan Department
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationTwo new enriched multiscale coarse spaces for the Additive Average Schwarz method
346 Two new enriched multiscale coarse spaces for the Additive Average Schwarz method Leszek Marcinkowski 1 and Talal Rahman 2 1 Introduction We propose additive Schwarz methods with spectrally enriched
More informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
More information