A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
|
|
- Alyson Dennis
- 5 years ago
- Views:
Transcription
1 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 to compute singular solutions of Poisson equations on a polygonal domain with re-entrant angles. Singularities are eliminated and only the regular part of the solution that is in H is computed. The stress intensity factor and the solution can be computed as a post processing step. This method is extended to problems with crack singularities and to a higher-order method for smooth data in [9]. In this paper, we study the Poisson equation with mixed boundary conditions. Examples with various singular points and numerical results are presented. Introduction Assume that Ω R is a polygonal domain. Let Γ D and Γ N be a partition of the boundary of Ω such that Ω = Γ D Γ N and Γ D Γ N =. For simplicity, assume that Γ D is not empty (i.e., meas(γ D ) 0). Let ν denote the outward unit vector normal to the boundary. For a given function f L (Ω), consider the Poisson equation with homogeneous mixed boundary conditions: u = f in Ω, u = 0 on Γ D, (.) u = 0 on Γ ν N, where stands for the Laplacian operator. Solution of (.) has singular behavior near corners even when f is very smooth. Such singular behavior affects the accuracy of the finite element method throughout the whole domain. This work was sponsored in part by KRF C0004. Department of Mathematics, Purdue University, 50 N. University Street, West Lafayette, IN (zcai@math.purdue.edu). Department of Applied Mathematics, Changwon National University, Changwon, , Republic of Korea (sckim@sarim.changwon.ac.kr). Department of Mathematics, Kyungpook National University (skim@knu.ac.kr). Department of Applied Mathematics, Changwon National University (ksr@changwon.ac.kr)
2 Solutions of many elliptic boundary value problems on polygonal domains have a singular function representation: the sum of singular functions and the regular part of the solution. This property has been explored in several ways to design accurate finite element methods in the presence of corner singularities. One approach is the so-called Singular Function Method (SFM) (see, e.g., [3]) that augments singular functions to both the trial and test spaces. Another is the so-called Dual Singular Function Method (DSFM) (see, e.q., [, 4, 3, 6, ]) that augments singular functions to the trial space and the corresponding dual singular functions to the test space. Recently we also use this property in order to calculate accurate finite element approximations to both the solution and the stress intensity factors, that are coefficients of singular functions in the singular function representation of the solution. By using the dual singular functions and a particularly chosen cut-off function, we are able to deduce a well-posed variational problem for the regular part of the solution. Hence, standard finite element approximation of this problem is optimally accurate. The stress intensity factors and the solution can then be computed with optimal accuracy as a post processing step. The Poisson equation with homogeneous Dirichlet boundary conditions was studied in [7]. This note studies mixed boundary conditions. Singular functions for mixed boundary conditions differ from those for Dirichlet boundary conditions in both form and angles of corners. In section, singular function representation of the solution for various boundary conditions is presented and a variational problem for regular part of the solution is derived. In section 3, we introduce a finite element approximation and estimate its error bound. Finally, in section 4, we present two examples with various singular functions and their numerical results. We will use the standard notation and definitions for the Sobolev spaces H t (B) for t 0; the standard associated inner products are denoted by (, ) t,b, and their respective norms and seminorms are denoted by t,b and t,b. The space L (B) is interpreted as H 0 (B), in which case the inner product and norm will be denoted by (, ) B and B, respectively. H D (B) = {u H (B) : u = 0 on Γ D }. Singular Function Representations Let ω,, ω M be internal angles of Ω satisfying π/ < ω < π if boundary condition changes its type, π < ω < π otherwise and denote by v ( =,, M) the corresponding vertices. Let the polar co-ordinates (r, θ ) be chosen at the vertex v so that the internal angle ω is spanned counterclockwise by two half-lines θ = 0 and θ = ω. Below is a list of singular functions at v depending on boundary conditions:
3 D/D If ω > π, there is a singular function of the form s, (r, θ ) = r π ω sin πθ ω ; (.) N/N If ω > π, there is a singular function of the form s, (r, θ ) = r π ω cos πθ ω ; (.) D/N If π < ω 3π, there is a singular function of the form s, (r, θ ) = r π ω sin πθ ω ; (.3) If 3π < ω < π, there are two singular functions of the form s, (r, θ ) = r π ω sin πθ and s ω, 3 (r, θ ) = r 3π ω sin 3πθ ω ; (.4) N/D If π < ω 3π, there is a singular function of the form s, (r, θ ) = r π ω cos πθ ω ; (.5) If 3π < ω < π, there are two singular functions of the form s, (r, θ ) = r π ω cos πθ and s ω, 3 (r, θ ) = r 3π ω cos 3πθ ω. (.6) Here, D/D and N/N mean that type of boundary conditions remains unchanged while D/N and N/D mean that type of boundary conditions changes passing the vertex v. For convenience, we denote index set of singular functions by L. Hence, {} for (.) and (.), L = { } for (.3) and (.5), {, 3 } for (.4) and (.6). 3π + ε ω H (Ω), s, H + π ε ω (Ω), and s, It is easy to see that s, 3 for any ε > 0. Hence, singular functions at the vertex v belong to either H + π ε ω (Ω) H + π ω ε (Ω) 3
4 for D/D and N/N vertex or H + π ω ε (Ω) for D/N and N/D vertex. This indicates that the solution of Poisson equation (.) is in H + π ω ε where ω = max M ˆω and ω ˆω = ω if v is D/D or N/N vertex, if v is D/N or N/D vertex. To deduce an equation for the regular part of the solution, we need to use the so-called dual singular functions that are defined as follows: for l L, s, l (r, θ ) = r ω sin θ and s, l (r, θ ) = r ω are the dual singular functions corresponding to ω s, l (r, θ ) = r sin θ and s, l (r, θ ) = r ω respectively. We will also need cut-off functions. To this end, set ω ω cos ω θ, (.7) cos ω θ (.8) B (t ; t ) = {(r, θ ) : t < r < t and 0 < θ < ω } Ω and B (t ) = B (0; t ). A family of cut-off functions of r, η ρ (r ), is then defined as follows: in B ( ρ { R), ( ) ( ) 3 ( ) } 5 η ρ (r ) = 5 8 4r 3 + 4r ρ R 3 ρ R 4r 3 in B 5 ρ R ( ρ R; ρ R), 0 in Ω \ B (ρ R) where ρ is a parameter in (0, ] and R R is a fixed number so that the η s, l has the same boundary condition as u. We assume that R is small enough so that the intersection of either B (ρ R) and B i (R) or B (R) and B i (ρ i R) for i is empty. It is well known ([, 0, ]) that the solution of problem (.) has the following singular function representation: u = w + M = l L λ, l η ρ (r )s, l (r, θ ) (.9) where w H (Ω) HD (Ω) is the regular part of the solution and λ, l R are the stress intensity factors that can be expressed in terms of u by the following extraction formulas ([5, 4]): λ, l = ( ) fη ρ s, l dx + u (η ρ s, l ) dx. (.0) Ω Ω 4
5 Moreover, the following regularity estimate holds: w + M = l L λ, l C R f. (.) In the remainder of this section, we derive a well-posed problem for w. To this end, assume that ρ in (.9) belongs to (0, ] and denote cut-off functions with bigger supports by η (r ) = η (r ). Choosing η ρ (r ) = η (r ) in (.0) gives Substituting u = w + M λ, l = (u, (η s, l )) + (f, η s, l ). i= λ, l = (w, (η s, l )) + (f, η s, l ) + k L i λ i, k η ρi (r i )s i, k (r i, θ i ) into the above equation yields M i= k L i λ i, k (η ρi s i, k, (η s, l )). (.) When i =, the support of η ρ (r ) for 0 < ρ is B (ρ R) on which η =. Since s, l is harmonic, then for all k L, (η ρi s i, k, (η s, l )) = 0. When i, by the assumption that B i (ρ i R) B (R) = we have that (η ρi s i, k, (η s, l )) = 0, k L i. Hence, we have established the following extraction formulas of λ, l in terms of w: λ, l = (w, (η s, l )) B (R;R) + (f, η s, l ) B (R). (.3) Using (.3) and substituting (.9) into the Poisson equation, we obtain an integrodifferential equation for w: w M = l L (w, (η s, l )) B (R;R) (η ρ s, l ) = f + M (f, η s, l ) B (R) (η ρ s, l ) in Ω. l L = Multiplying the above equation by a test function v HD (Ω) = {v H (Ω) : v = 0 on Γ D }, integrating over Ω, and using integration by parts lead to the following variational problem: finding w H (Ω) HD (Ω) such that a(w, v) = g(v) v H D(Ω), (.4) 5
6 where the bilinear form a(, ) and linear form g( ) are defined by and a(w, v) = a s (w, v) + b(w, v), a s (w, v) = ( w, v), b(w, v) = M = g(v) = (f, v) + l L (w, (η s, l )) B (R;R)( (η ρ s, l ), v) B ( ρ R; ρ R) (.5) M (f, η s, l ) B (R)( (η ρ s, l ), v) B ( ρ R; ρ R) (.6) l L = Note that the second terms in the respective bilinear and linear forms provide a singular correction so that w H (Ω) for f L (Ω). Note also that the bilinear forms a(, ) are not symmetric. Lemma. For any 0 < ρ, we have that with C = with C = 0 7 and that (η s, l ) B (R;R) C R ω (.7) + ω η ρ s, l B (ρ R) C (ρ R) and (η ρ s, l ) B (ρ R) C 3 (ρ R) ω (.8) ω and C 3 = ( Cω +) ω +ω 4(+ω ( ( ) ) + ). Proof. This lemma can be established by an elementary calculation. In a similar fashion as in [7], we can prove the coercivity and continuity of the bilinear form a(, ) and the well posed-ness of problem (.4). Lemma. For 0 < ρ, the bilinear forms a(, ) are continuous and coercive in H D (Ω); i.e. there exist positive constants α, K, and K such that α φ a(φ, φ) + K φ (.9) for all φ H D (Ω) and that a(φ, ψ) K φ ψ (.0) for all φ and ψ in H D (Ω). Theorem. For 0 < ρ, we have that () if f L (Ω), then problem (.4) has a unique solution w H (Ω) H D (Ω). () there exists a positive constant γ such that for any φ H D (Ω). γ φ a(φ, ψ) sup (.) ψ HD (Ω) ψ 6
7 3 Finite Element Approximation This section presents standard finite element approximation on a quasi-uniform grid for w based on the variational problem in (.4). Approximations to the stress intensity factors and the solution of problem (.) can then be calculated according to (.3) and (.9), respectively. Error estimates are established in Theorem 3.. Let T h be a partition of the domain Ω into triangular finite elements; i.e., Ω = K T h K with h = max{diam K : K T h }. Assume that the triangulation T h is regular. Denote continuous piecewise linear finite element space by V h = {φ h C 0 (Ω) : φ h K is linear K T h and φ h = 0 on Γ D } H D(Ω). It is well known that inf ( φ φ h + h φ φ h ) C A h +t φ +t,ω (3.) φ h V h for any φ H D (Ω) H +t (Ω) and 0 t. The finite element approximation to problem (.4) is to find w h V h such that a(w h, v) = g(v) v V h. (3.) Approximations to the λ, l and the solution are calculated as follows: and λ h, l = (w h, (η s, l )) B (R;R) + (f, η s, l ) B (R) (3.3) u h = w h + M = l L λ h, l η ρ (r )s, l (r, θ ). (3.4) In order to establish the error bound in the L norm, we consider the following adoint problem of (.4) with a simplified linear form: find z HD (Ω) such that a(v, z) = (w w h, v) v H D(Ω). (3.5) The next lemma establishes the well-posedness of problem (3.5) and provides the regularity estimate for z. Lemma 3. For 0 < ρ, problem (3.5) has a unique solution z in HD (Ω). Moreover, there is a singular function representation z = w z + M = l L λ z, lη ρ s,l (3.6) where w z H (Ω) HD (Ω) and λz, l R satisfy the regularity estimate w z + M = l L λ z, l C R w w h. (3.7) 7
8 Proof. Similar to Theorem., the adoint problem in (3.5) has a unique solution in H D (Ω) and that there exists a positive constant γ such that γ ψ a(φ, ψ) sup φ HD (Ω) φ ψ H D(Ω). Let z be the solution of (3.5), by the Cauchy-Schwarz inequality we then have that z γ a(φ, z) sup = (w w h, φ) sup φ HD (Ω) φ γ φ HD (Ω) φ γ w w h. (3.8) It is easy to check that the solution, z HD (Ω), of problem (3.5) satisfies z = M = l L ( z, (η ρ s, l )) (η s, l ) (w w h ) in Ω. (3.9) Since the right-hand side of the above equation is at least in L (Ω), so is z. Therefore, z has the singular function representation where w z H (Ω) H D and z = w z + w z + M = M = l L λ z, l η ρ s, l, l L λ z, l C R z. Now, the regularity bound in (3.7) follows from the triangle and Cauchy-Schwarz inequalities, (3.8), and Lemma. that w z + M λ z, l C R z = l L C R ( M ( (η ρ s, l ), z) B(ρ R) (η s, l ) B(R;R) + w w h ) = l L C R ( M + ) w w h. ρ ω γ R C C 3 = l L This proves the inequality in (3.7) with M C C C 3 R = C R ( γ R ρ ω + ) = l L and, hence, the lemma. Now we are ready to establish error bounds for the finite element approximations. 8
9 Theorem 3. (i) For 0 < ρ, there exists a positive constant h 0 such that for all h h 0 (3.) has a unique solution w h in V h. Moreover, let w H (Ω) H D (Ω) be the solution of (.4), then we have the following error estimates: w w h C 4 h f and w w h C 5 h + π ω f. (3.0) (ii) Let λ, l and λ h, l be defined in (.3) and (3.3), respectively. Then λ, l λ h, l C ω R w wh C 6 R ω h + π ω f. (3.) (iii) Let u be the solution of (.) and u h be its approximation defined in (3.4), then we have the following error estimates: u u h C 7 h f and u u h C 8 h + π ω f. (3.) Proof. (i) We first establish error bounds in (3.0) for any solution to problem (3.) that may exist. Then, for f 0, the uniqueness of the solution to problem (.4) and the error bound in (3.0) imply that w h 0. Hence, (3.) has a unique solution w h in V h since it is a finite dimensional problem with the same number of unknowns and equations. To establish error bounds, note first the orthogonality property a(w w h, v) = 0 v V h. (3.3) By choosing v = w w h in (3.5) and using the orthogonality property in (3.3) and the continuity bound in (.0), we have that w w h = a(w w h, z) = a(w w h, z I h z) K w w h z I h z, (3.4) where I h z V h is the nodal interpolant of z. From the triangle inequality, approximation property (3.), the fact that (see []) and Lemma 3., one has η ρ s, l I h (η ρ s, l ) Ch ω, z I h z w z I h w z + M Ch w z + M = Ch = l L l L λ z, l η ρ s, l I h (η ρ s, l ) ω λ z, l C D h π ω w w h. Substituting this into (3.4) and dividing w w h on both sides give w w h K C D h π ω w wh. (3.5) 9
10 Now, it follows Lemma., orthogonality property (3.3), and inequality (3.5) that for any v V h α w w h a(w w h, w w h ) + K w w h = a(w w h, w v) + K w w h K w w h w v + K (K C D h π ω ) w w h, which, together with approximation property (3.), implies the validity of the first error bound in (3.0) with C 4 = α K C A C R for all h h 0. Here, α h 0 = ( K (K C D ) ) ω π. The second error bound in (3.0) is then a direct consequence of (3.5) with C 5 = C 4 K C D. (ii) Note from (.3) and (3.3) that λ, l λ h, l = (w w h, (η s, l )) B (R; R). Hence, (3.) follows from the Cauchy-Schwarz inequality, Theorem 3.(i), and Lemma. that λ, l λ h, l w w h (η s, l ) B (R; R) C 6 R ω h + π ω f with C 6 = C 5C. (iii) It follows from (.9) and (3.4) that u u h = (w w h ) + M = l L (λ, l λ h, l)η ρ s, l. By using the triangle inequality, Lemma., (3.0), and (3.), we have that u u h w w h + M λ, l λ h, l η ρ s, l,b(ρ R) = l L C 4 h f + M ω C 6 ρ (C ρ + C 3 R )h + π ω f. = l L Therefore, the first inequality of (3.) is valid with C 7 = C 4 + C 6 M ω ρ = l L (C ρ + C 3 R )h π ω. In a similar fashion, by Lemma., (3.0), and (3.), we may prove the M ω validity of the second inequality of (3.) with C 8 = C 5 + C 6 C. This completes the proof of the theorem. 0 ρ + = l L
11 4 Numerical Results In this section, numerical results for Poisson equations with mixed boundary conditions are presented. One example is defined on the unit square and the other on a domain with re-entrance corner. u = 0 ν u = 0 (0, 0) (, 0) (a) Example (, ) (, ) u = 0 ν (0, 0) u = 0 (, ) (b) Example Figure : Computational domains and boundary conditions Example. Consider the Poisson equation in (.) with mixed boundary conditions on the unit square Ω = {(x, y) R : 0 < x <, 0 < y < }. Homogeneous Neumann boundary is Γ N = {(x, 0) R : 0 < x < /} and homogeneous Dirichlet boundary is Γ D = Ω \ Γ N (see Fig.(a)). This problem has a geometric singularity at boundary point (/, 0), where the boundary conditions change from Dirichlet to Neumann with an internal angle ω = π. More specifically, the corresponding singular function has the form s = r θ sin( ). Let η be the cut-off function defined in section with R = /4 and choose the righthand side function in (.) to be f = sin(πx) [ π y (y ) + (3y ) ] (η s). Then the exact solution of the underlying problem is u = w + η ρ s, where η ρ is the cut-off function with R = /4 and 0 < ρ and w = sin(πx)y (y ) + (η η ρ ) s is the regular part of the solution. Numerical results are presented in Table,, and 3, respectively, that confirm theoretical estimates.
12 ρ =.00 ρ = 0.50 Mesh Size L -norm Rate L -norm Rate 8.333e e e e e e e e e e Table. The discrete L -norm errors and the convergence rates for w. ρ =.00 ρ = 0.50 H -norm Rate H -norm Rate e e e e e e e e e e Table. The discrete H seminorm and the convergence rates for w. ρ =.00 ρ = 0.50 λ λ h Rate λ λ h Rate e e e e e e e e e e Table 3. The absolute value errors and the convergence rates for λ. Example. In this example, we consider a polygonal domain with a re-entrant corner (see Fig.(b)): Ω = {(x, y) R : < x <, < y < } \ T,
13 where T = {(x, y) R : 0 < x <, x < y < 0}. Homogeneous Neumann boundary is Γ N = {(x, 0) R : 0 < x <, y = 0} and homogeneous Dirichlet boundary is Γ D = Ω \ Γ N. At the origin, the internal angle is ω = 7π and the boundary conditions 4 change from Neumann to Dirichlet. Hence, there are two singular functions at the origin: s = r 7 cos θ 7 and s 3 = r 6 7 cos 6θ 7. Let η ρ be the cut-off function with R = /4 and 0 < ρ and choose the right-hand side function in (.) to be where f = f (η (s + s 3 )), (y )y (x )(6y ) if y 0, f = (y + )y (x + )(6xy + x + y + 6y) if y < 0, The exact solution of the underlying problem is then where the regular part w is given by w = (η η ρ )(s u = w + η ρ s + η ρ s 3 + s 3 ) + (x )(y )y if y 0, (x + )(y + )(x + y)y if y 0. Note that the function w is in H (Ω), but not in H 3 (Ω). Numerical results for the discretization accuracy of the finite element approximation to w are given in Tables 4 and 5. Results for the stress intensity factors are contained in Tables 6 and 7. ρ =.00 ρ = 0.50 Mesh Size L -norm Rate L -norm Rate e e e e e e e e Table 4. The discrete L -norm errors and the convergence rates for w. 3
14 ρ =.00 ρ = 0.50 H -norm Rate H -norm Rate e e e e e e e e Table 5. The discrete H seminorm and the convergence rates for w. ρ =.00 ρ = 0.50 λ λ,h Rate λ λ,h Rate 6.635e e e e e e e e Table 6. The absolute value errors and the convergence rates for λ. ρ =.00 ρ = 0.50 λ 3 λ 3,h Rate λ 3 λ 3,h Rate.50e e e e e e e e Table 7. The absolute value errors and the convergence rates for λ 3. References [] I. Babuska and A. Miller, The post-processing approach in the finite element method - part : The calculation of stress intensity factors, Int. J. Numer. Methods Engrg., 0 (984), -9. [] I. Babuska and M. Suri, The h-p version of the finite element method with quasiuniform meshes, RAIRO Modél. Math. Anal. Numér., (987),
15 [3] M. Bourlard, M. Dauge, M.-S. Lubuma, and S. Nicaise, Coefficients of the singularities for elliptic boundary value problems on domains with conical points III. Finite element methods on polygonal domains, SIAM Numer. Anal., 9 (99), [4] H. Blum and M. Dobrowolski, On finite element methods for elliptic equations on domains with corners, Computing, 8 (98), [5] H. Blum, Numerical treatment of corner and crack singularities, Dissertation, Bonn, 98. [6] S. Brenner, Multigrid methods for the computation of singular solutions and stress intensity factor I: Corner singularities, Math. Comp., 68 (6), (999), [7] Z. Cai and S. Kim, A finite element method using singular functions for the Poisson equation: Corner singularities, SIAM J. Numer. Anal., 39:(00), [8] Z. Cai, S. Kim and B. Shin, Solution methods for the poisson equation: Corner singularities, SIAM J. Sci. Comp., 39:(00), [9] Z. Cai, S. Kim and G. Woo,A finite element method using singular functions for the Poisson equation: Crack singularities, Numer. Lin. Alg. Appli. 9:(00), [0] M. Dauge, Elliptic Boundary Value Problems on Corner Domains, Lecture Notes in Mathematics 34, Springer-Verlag, Berlin-Heidelberg, 988. [] M. Daoua, Equations Intégrales pour un Probleme Singulier dans le Plan, These de Troisieme Cycle, Universite Pierre et Marie Curie, Paris, 977. [] M. Dobrowolski, Numerical Approximation of Elliptic Interface and Corner Problems, Habilitation-schrift, Bonn, 98. [3] G. J. Fix, S. Gulati, and G. I. Wakoff, On the use of singular functions with finite elements approximations, J. Comput. Phy., 3 (973), [4] P. Grisvard, Elliptic Problems in Nonsmooth Domains, Pitman, Boston, MA,
A 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 informationFINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS
J. KSIAM Vol.22, No.2, 101 113, 2018 ttp://dx.doi.org/10.12941/jksiam.2018.22.101 FINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS DEOK-KYU JANG AND JAE-HONG PYO DEPARTMENT
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 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 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 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 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 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 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 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 informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA 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 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 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 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 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 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 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 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 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 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 Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More 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 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 informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
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 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 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 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 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 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 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 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 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 informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More 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 informationA gradient recovery method based on an oblique projection and boundary modification
ANZIAM J. 58 (CTAC2016) pp.c34 C45, 2017 C34 A gradient recovery method based on an oblique projection and boundary modification M. Ilyas 1 B. P. Lamichhane 2 M. H. Meylan 3 (Received 24 January 2017;
More informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More 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 informationarxiv: v1 [math.na] 17 Feb 2016
Adapted numerical methods for the numerical solution of the Poisson equation with L 2 boundary data in non-convex domains Thomas Apel Serge Nicaise Johannes Pfefferer arxiv:62.5397v [math.na] 7 Feb 26
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 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 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 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 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 CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS 1. INTRODUCTION
Trends in Mathematics Information Center for Mathematical Sciences Volume 4, Number 2, December 21, Pages 39 44 A CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS
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 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 informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More 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 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 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 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 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 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 informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
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 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 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 informationOn angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic
On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationThe mortar element method for quasilinear elliptic boundary value problems
The mortar element method for quasilinear elliptic boundary value problems Leszek Marcinkowski 1 Abstract We consider a discretization of quasilinear elliptic boundary value problems by the mortar version
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
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 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 informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More 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 informationNumerische Mathematik
Numer. Math. (1996) 75: 59 77 Numerische Mathematik c Springer-Verlag 1996 Electronic Edition A preconditioner for the h-p version of the finite element method in two dimensions Benqi Guo 1, and Weiming
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
More informationBOUNDARY LAYERS FOR TRANSMISSION PROBLEMS WITH SINGULARITIES
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 14, pp. 1 16. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) BOUNDARY
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 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 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 informationThe maximum angle condition is not necessary for convergence of the finite element method
1 2 3 4 The maximum angle condition is not necessary for convergence of the finite element method Antti Hannukainen 1, Sergey Korotov 2, Michal Křížek 3 October 19, 2010 5 6 7 8 9 10 11 12 13 14 15 16
More informationNodal O(h 4 )-superconvergence of piecewise trilinear FE approximations
Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen
More informationWEIGHTED-NORM FIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PROBLEMS WITH CORNER SINGULARITIES
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
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 informationTHE BEST L 2 NORM ERROR ESTIMATE OF LOWER ORDER FINITE ELEMENT METHODS FOR THE FOURTH ORDER PROBLEM *
Journal of Computational Mathematics Vol.30, No.5, 2012, 449 460. http://www.global-sci.org/jcm doi:10.4208/jcm.1203-m3855 THE BEST L 2 NORM ERROR ESTIMATE OF LOWER ORDER FINITE ELEMENT METHODS FOR THE
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 informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationMEYERS TYPE ESTIMATES FOR APPROXIMATE SOLUTIONS OF NONLINEAR ELLIPTIC EQUATIONS AND THEIR APPLICATIONS. Yalchin Efendiev.
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX MEYERS TYPE ESTIMATES FOR APPROXIMATE SOLUTIONS OF NONLINEAR ELLIPTIC EUATIONS AND THEIR APPLICATIONS
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 informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationPriority Program 1253
Deutsche Forschungsgemeinschaft Priority Program 1253 Optimization with Partial Differential Equations Thomas Apel and Dieter Sirch L 2 -Error Estimates for the Dirichlet and Neumann Problem on Anisotropic
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 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 informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More 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 informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
More 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 informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
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 informationError estimates for a finite volume element method for parabolic equations in convex polygonal domains
Error estimates for a finite volume element method for parabolic equations in convex polygonal domains P Chatzipantelidis, R D Lazarov Department of Mathematics, Texas A&M University, College Station,
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 informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More information