1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
|
|
- Ronald Morgan
- 5 years ago
- Views:
Transcription
1 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 for the Stokes equations based on a weak Galerkin (WG) method introduced in [15]. Discrete divergence free bases are constructed explicitly for the lowest order weak Galerkin elements in two and three dimensional spaces. These basis functions can be derived on general meshes of arbitrary shape of polygons and polyhedrons. With the divergence free basis derived, the discrete divergence free WG scheme can eliminate pressure variable from the system and reduces a saddle point problem to a symmetric and positive definite system with many fewer unknowns. Numerical results are presented to demonstrate the robustness and accuracy of this discrete divergence free WG method. Key words. Weak Galerkin, finite element methods, the Stokes equations, divergence free. AMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50 1. Introduction. The Stokes problem seeks unknown functions u and p satisfying (1.1) (1.2) (1.3) A u + p = f in Ω, u = 0 in Ω, u = 0 on Ω, where Ω is a polygonal domain in R d with d = 2, 3 and A is a symmetric and positive definite d d matrix-valued function in Ω. For the nonhomogeneous boundary condition u = g on Ω, one can use the standard procedure by letting u = u 0 + u g. u g is a known function satisfying u g = g on Ω and u 0 is zero at Ω and satisfies (1.1)-(1.2) with different right hand sides. The weak form in the primary velocity-pressure formulation for the Stokes problem (1.1) (1.3) seeks u [H 1 0 (Ω)] d and p L 2 0(Ω) satisfying (1.4) (1.5) (A u, v) ( v, p) = (f, v), v [H 1 0 (Ω)] d ( u, q) = 0, q L 2 0(Ω). In the standard finite element methods for the Stokes and the Navier-Stokes equations, both pressure and velocity are approximated simultaneously. The primitive system is a large saddle point problem. Numerical solvers for such indefinite Computer Science and Mathematics Division Oak Ridge National Laboratory, Oak Ridge, TN, 37831,USA (mul1@ornl.gov). This research was supported in part by National Science Foundation Grant DMS Division of Mathematical Sciences, National Science Foundation, Arlington, VA (jwang@ nsf.gov). The research of Wang was supported by the NSF IR/D program, while working at the Foundation. However, any opinion, finding, and conclusions or recommendations expressed in this material are those of the author and do not necessarily reflect the views of the National Science Foundation. Department of Mathematics, University of Arkansas at Little Rock, Little Rock, AR (xxye@ualr.edu). This research was supported in part by National Science Foundation Grant DMS
2 2 systems are usually less effective and robust than solvers for definite systems. On the other hand, the divergence-free finite element method, discrete or exact, computes numerical solution of velocity by solving a symmetric positive definite system in a divergence-free subspace. It eliminates the pressure from the coupled equations and hence significantly reduces the size of the system. The divergence-free method is particularly attractive in the cases where the velocity is the primary variable of interest, for example, the groundwater flow calculation. The main tasks in the implementation of the divergence-free method are to understand divergence-free subspaces, weakly or exactly, and to construct bases for them. Many finite element methods, continuous [2, 3, 8] and discontinuous [1, 4, 11, 12, 16], have been developed and analyzed for the Stokes and the Navier-Stokes equations. Divergence-free basis for different finite element methods have been constructed [5, 6, 7, 9, 10, 17, 18, 19]. A weak Galerkin finite element method was introduced in [15] for the Stokes equations in the primal velocity-pressure formulation. This method is designed by using discontinuous piecewise polynomials on finite element partitions with arbitrary shape of polygons/polyhedra. Weak Galerkin methods were first introduced in [13, 14] for second order elliptic equations. In general, weak Galerkin finite element formulations for partial differential equations can be derived naturally by replacing usual derivatives by weakly-defined derivatives in the corresponding variational forms, with an option of adding a stabilization term to enforce a weak continuity of the approximating functions. Therefore the weak Galerkin method developed in [15] for the Stokes equations naturally has the form: find u h V h and p h W h satisfying (1.6) (1.7) (A w u h, w v) + s(u h, v) ( w v, p h ) = (f, v), ( w u h, q) = 0 for all the test functions v V h and q W h where V h and W h will be defined later. The stabilizer s(u h, v) in (1.6) is parameter independent. Let D h be a discrete divergence free subspace of V h such that ( w v, q) = 0 for ant q W h. Then the discrete divergence free WG formulation is to find u h D h satisfying (1.8) (A w u h, w v) + s(u h, v) = (f, v), v D h. System (1.8) is symmetric and positive definite with many fewer unknowns. The main purpose of this paper is to construct bases for D h in two and three dimensional spaces. A unique feature of these divergence free basis functions is that they can be obtained on general meshes such as hybrid meshes or meshes with hanging nodes. Numerical examples in two dimensional space are provided to confirm the theory. Although the Stokes equations is considered, the divergence free basis can be use for solving the Navier-Stokes equations. 2. A Weak Galerkin Finite Element Method. In this section, we will review the WG method for the Stokes equations introduced in [15] with k = 1. Let T h be a partition of the domain Ω consisting mix of polygons satisfying a set of conditions specified in [14]. In addition, we assume that all the elements T T h are convex. Denote by F h the set of all edges in 2D or faces in 3D in T h, and let F 0 h = F h\ Ω be the set of all interior edges or faces. We define a weak Galerkin finite element space for the velocity as follows V h = { v = {v 0, v b } : {v 0, v b } T [P 1 (T )] d [P 0 (e)] d, e T, v b = 0 on Ω }.
3 We would like to emphasize that there is only a single value v b defined on each edge in 2D and face in 3D. For the pressure variable, we have the following finite element space W h = { q : q L 2 0(Ω), q T P 0 (T ) }. For a given v V h, a weak gradient and a weak divergence are defined locally on each T T h as follows. Definition 2.1. A weak gradient, denoted by w, is defined as the unique polynomial ( w v) [P 0 (T )] d d for v V h satisfying the following equation, (2.1) ( w v, q) T = (v 0, q) T + v b, q n T, q [P 0 (T )] d d, 3 Definition 2.2. A weak divergence, denoted by ( w ), is defined as the unique polynomial ( w v) P 0 (T ) for v V h that satisfies the following equation (2.2) ( w v, ϕ) T = (v 0, ϕ) T + v b n, ϕ T, ϕ P 0 (T ). Denote by Q 0 the L 2 projection operator from [L 2 (T )] d onto [P 1 (T )] d and denote by Q b the L 2 projection from [L 2 (e)] d onto [P 0 (e)] d. Let (u; p) be the solution of (1.1)-(1.3). Define Q h u = {Q 0 u, Q b u} V h. let Q h be the local L 2 projections onto P 0 (T ). We introduce three bilinear forms as follows s(v, w) = T Q bv 0 v b, Q b w 0 w b T, a(v, w) = h 1 (A w v, w w) T + s(v, w), b(v, q) = ( w v, q) T. Algorithm 1. A numerical approximation for (1.1)-(1.3) can be obtained by seeking u h = {u 0, u b } V h and p h W h such that (2.3) (2.4) a(u h, v) b(v, p h ) = (f, v 0 ), b(u h, q) = 0, q W h. v V h Define (2.5) v 2 = a(v, v). The following optimal error estimates have been derived in [15]. Theorem 2.1. Let (u; p) [H 1 0 (Ω) H 2 (Ω)] d (L 2 0(Ω) H 1 (Ω)) and (u h ; p h ) V h W h be the solution of (1.1)-(1.3) and (2.3)-(2.4), respectively. Then, the following error estimates hold true (2.6) (2.7) Q h u u h + Q h p p h Ch( u 2 + p 1 ), Q 0 u u 0 Ch 2 ( u 2 + p 1 ).
4 4 Define a discrete divergence free subspace D h of V h by (2.8) D h = {v V h ; b(v, q) = 0, q W h }. By taking the test functions from D h, the weak Galerkin formulation (2.3)-(2.4) is equivalent to the following divergence-free weak Galerkin finite element scheme. Algorithm 2. A discrete divergence free WG approximation for (1.1)-(1.3) is to find u h = {u 0, u b } D h such that (2.9) a(u h, v) = (f, v 0 ), v = {v 0, v b } D h. System (2.9) is symmetric and positive definite with many fewer unknowns. It can be solved effectively by many existing solvers. The main task of this paper is to construct basis for D h. In the next two sections, discrete divergence free bases will be constructed explicitly for two and three dimensional spaces. 3. Construction of Discrete Divergence Free Basis for Two Dimensional Space. For a given partition T h, let V 0 h be the set of all interior vertices. Let N F = card(f 0 h ), N V = card(v 0 h ) and N K = card(t h ). It is known based on the Euler formula that for a partition consisting of convex polygons, then (3.1) N F + 1 = N V + N K. For a mesh T h with hanging nodes, the relation in (3.1) is still true if we treat the hanging nodes as vertices. First we need to derive a basis for V h. For each T T h and any v = {v 0, v b } V h, v 0 is a vector function with two components and each component is a linear function. Therefore there are six linearly independent linear functions Φ j+1, Φ j+2,, Φ j+6 in V h such that they are nonzero only at the interior of element T. For each e i F 0 h, v b is a constant vector function. Thus, there are two linearly independent constant functions Ψ i,1 and Ψ i,2 in V h which take nonzero value only on e i. Then it is easy to see that (3.2) V h = span{φ 1,, Φ 6NK, Ψ 1,1, Ψ 1,2,, Ψ NF,1, Ψ NF,2}. For a given function v = {v 0, v b } V h, it is easy to see that v 0 can be spanned by the basis functions Φ i and v b by the basis functions Ψ j,k. Next, we will find the dimension of D h. Since the dimension for pressure space W h is N K 1, it follows from (3.1) that (3.3) dim(d h ) = dim(v h ) dim(w h ) = 6N K + 2N F N K + 1 = 6N K + N F + N V. Lemma 3.1. The basis functions Φ 1,, Φ 6NK of V h in (3.2) are in D h and linearly independent. Proof. Let Φ i = {Φ i,0, Φ i,b }. The definition of Φ i implies Φ i,b = 0. For any q W h, it follows from (2.2), q = 0 and Φ i,b = 0, b(φ i, q) = ( w Φ i, q) T = ( (Φ i,0, q) T + Φ i,b n, q T ) = 0,
5 which proves the lemma since the linear independence of Φ 1,, Φ 6NK is obvious. For any e i F 0 h,let ψ i,1 and ψ i,2 be two basis functions of V h associated with e i. Let n ei and t ei be a normal vector and a tangential vector to e i respectively. Define Υ i = C 1 Ψ i,1 + C 2 Ψ i,2 such that Υ i ei = t ei. Obviously Υ i V h is only nonzero on e i. Lemma 3.2. Functions Υ 1,, Υ NF V h are in D h and linearly independent. Proof. Let Υ i = {Υ i,0, Υ i,b }. For any q W h, it follows from (2.2) and q = 0, b(υ i, q) = ( w Υ i, q) T = ( (Υ i,0, q) T + Υ i,b n, q T ) = Υ i,b n, q T = 0, 5 where we use the fact t ei n = 0. Since Υ i is only nonzero on e i, Υ 1,, Υ NF linearly independent. We completed the proof. are For a given interior vertex P i Vh 0, assume that there are r elements having P i as a vertex which form a hull H Pi as shown in Figure??. Then there are r interior edges e j (j = 1,, r) associated with H Pi. Let n ej be a normal vector on e j such that normal vectors n ej j = 1,, r are counterclockwise around vertex P i as shown in Figure??. For each e j, let Ψ j,1 and Ψ j,2 be the two basis functions of V h which is only nonzero on e j. Define Θ j = C 1 Ψ j,1 + C 2 Ψ j,2 V h such that Θ j ej = n ej. Define Λ i = r j=1 1 e Θ j j. Lemma 3.3. Functions Λ 1,, Λ NV V h are in D h and linearly independent. Proof. Suppose that there exist constants c 1,, c NV such that N V i=1 c iλ i = 0. Let Λ l be associated with a hull H Pl such that there exists e m as one of the interior edges of H Pl and edge e m has a boundary node as one of its end points. Since Λ i = 0 on e m for i l, we have N V 0 = c i Λ i = e m c l Λ l = e m c l Θ m = c l n em, e m i=1 which implies c l = 0. By this way, we can prove that all c i = 0 and Λ 1,, Λ NV are linearly independent. Next, we will show that b(λ i, q) = 0 for all q W h. Let q j W h such that q j = 1 on T j T h and q j = 0 otherwise. So we only need to show that (3.4) b(λ i, q j ) = 0, j = 1,, N K. Let Λ i and q j be associated with hull H Pi and element T j respectively. If T j T h is not in H Pi, we easily have b(λ i, q j ) = 0. If T j H Pi, let e s and e s+1 be its two edges
6 6 in H Pi shown in Figure??, b(λ i, q j ) = ( w Λ i, q j ) T = ( w Λ i, q j ) Tj We proved Λ i D h for i = 1,, N V. = (Λ i,0, q) Tj + Λ i,b n, q T j 1 = e s e s n 1 e s n + e s+1 e s+1 n e s+1 n = 0. Theorem 3.4. Let D h be defined in (2.8). Then for two dimensional space, D h is spanned by the following basis functions, (3.5) D h = Span{Φ 1,, Φ 6NK, Υ 1,, Υ NF, Λ 1,, Λ NV }. N V Proof. The number of the functions in the right hand side of (3.5) is 6N K + N F + which is equal to dim(d h ) due to (3.3). Next, we prove that Φ 1,, Φ 6NK, Υ 1,, Υ NF, Λ 1,, Λ NV are linearly independent. Since Φ i take zero value on all f F h, Φ i will be linearly independent to all Υ l and Λ m. Suppose (3.6) C 1 Υ C NF Υ NF + C NF +1Λ C NF +N V Λ NV = 0. Multiplying (3.6) by Υ i and integrating over e i, we have C i e i = 0, where we use the fact t ei n e = 0. Thus we can obtain C i = 0 for i = 1,, N F. By Lemma 3.3, we can prove C i = 0 for i = N F + 1,, N F + N V. The proof of the lemma is completed. 4. Construction of Discrete Divergence Free Basis for Three Dimensional Space. Let T h be a partition of Ω R 3 consisting polyhedrons without hanging nodes. Recall N F = card(f 0 h ), N V = card(v 0 h ) and N K = card(t h ). Denote by E h all the edges in T h and let E 0 h = E h\ Ω. Let N E = card(e 0 h ). It is known based on the Euler formula that for a partition consisting of convex polyhedrons, then (4.1) N V + N F + 1 = N E + N K. For each T T h and any v = {v 0, v b } V h, v 0 is a vector function with three components and each component is a linear function. Therefore there are twelve linearly independent linear functions Φ j+1, Φ j+2,, Φ j+12 in V h such that they are nonzero only at the interior of element T. For each face f i Fh 0, v b is a constant vector function with three component. Thus there are three linearly independent constant vector functions Ψ i,1, Ψ i,2 and Ψ i,3 in V h which take nonzero value only on the face f i. Then it is easy to see that (4.2) V h = span{φ 1,, Φ 12NK, Ψ 1,1, Ψ 1,2, Ψ 1,3, Ψ NF,1, Ψ NF,2, Ψ NF,3}.
7 7 Since the dimension for pressure space W h is N K 1, (4.1) implies (4.3) dim(d h ) = dim(v h ) dim(w h ) = 12N K + 3N F N K + 1 = 12N K + 2N F + N E N V. Lemma 4.1. The functions Φ 1,, Φ 12NK in (4.2) are in D h and linearly independent. Proof. Let Φ i = {Φ i,0, Φ i,b }. The definition of Φ i implies Φ i,b = 0. For any q W h, it follows from (2.2), q = 0 and Φ i,b = 0 that for any q W h, b(φ i, q) = ( w Φ i, q) T = ( (Φ i,0, q) T + Φ i,b n, q T ) = 0, which finished the proof of the lemma since the linear independence of Φ 1,, Φ 12NK is obvious. For any face f i F 0 h, let n be a unit normal vector of f i and let t 1 and t 2 be two linearly independent unit tangential vectors to the face f i. Define Υ i,1 = C 1,1 Ψ i,1 + C 1,2 Ψ i,2 + C 1,3 Ψ i,3 and Υ i,2 = C 2,1 Ψ i,1 + C 2,2 Ψ i,2 + C 2,3 Ψ i,3 such that Υ i,1 fi = t 1 and Υ i,1 fi = t 2 respectively. Obviously Υ i,1 and Υ i,2 are in V h and only nonzero on f i. Lemma 4.2. Functions Υ 1,1, Υ 1,2, Υ NF,1, Υ NF,2 V h are in D h and linearly independent. Proof. Let Υ i,j = {Υ i,j,0, Υ i,j,b } with j = 1, 2. For any q W h, it follows from (2.2) and q = 0 for j = 1, 2, b(υ i,j, q) = ( w Υ i,j, q) T = ( (Υ i,j,0, q) T + Υ i,j,b n, q T ) = Υ i,j,b n, q T = 0, where we use the fact t j n = 0 with j = 1, 2. Since Υ i,j is only nonzero on f i, Υ 1,1, Υ 1,2, Υ NF,1, Υ NF,2 are linearly independent. We completed the proof. For a given interior edge E i E h, assume there are r elements having E i as one of their edges which form a solid denoted by S Ei shown in Figure 4.1. Then there are r interior faces f j (j = 1,, r) in S Ei. Let n j be a unit normal vector on the face f j such that normal vectors n j j = 1,, r form oriented loop around the interior edge E i shown in Figure 4.1. For each f j, let Ψ j,1, Ψ j,2 and Ψ j,3 be the three basis functions of V h which are only nonzero on f j. Define Θ j = C 1 Ψ j,1 + C 2 Ψ j,2 + C 3 Ψ j,3 V h such that Θ j fj = n j. Define Λ i = r j=1 1 f Θ j j. Lemma 4.3. Functions Λ 1,, Λ NE V h are in D h. Proof. Let Λ i and q j be associated with hull S Ei and element T j respectively. If T j T h is not in S Ei, we easily have b(λ i, q j ) = 0. If T j S Ei, let faces f s and f s+1
8 8 be its two faces in S Ei shown in Figure 4.1, b(λ i, q j ) = ( w Λ i, q j ) T = ( w Λ i, q j ) Tj We proved Λ i D h for i = 1,, N E. = (Λ i,0, q) Tj + Λ i,b n, q T j 1 = f s f s n 1 f s n + f s+1 f s+1 n f s+1 n = 0 Unfortunately, Λ 1,, Λ NE are linearly dependent. Let P i be an interior vertex in T h and G Pi be a hull formed by the elements T T h sharing P i. Let e j Eh 0, j = 1,, t with P i as one of its end point and Λ j, j = 1,, t be the discrete divergence free functions associated with e j. With appropriate choosing n j in defining Λ j, one can prove that t j=1 Λ j = 0. However, if we eliminate one function from {Λ 1,, Λ t } randomly, say Λ 1, we will prove that {Λ 2,, Λ t } are linearly independent in the following lemma. Lemma 4.4. Functions Λ 2,, Λ t are linearly independent. Proof. Let f Fh 0 be an interior face in G P i with e 1 and e 2 as its two edges in G Pi. The definition of Λ j implies that only Λ 1 and Λ 2 are nonzero on f. Suppose that there exist constants c 2,, c t such that t i=2 c iλ i = 0. Then we have 0 = t i=2 f c i Λ i = f c 2 Λ 2 = f c 2 Θ f = c 2 n f, which implies c 2 = 0. By this way, we can prove that all c i = 0 and Λ 2,, Λ t are linearly independent. We start with {Λ 1,, Λ NE } and eliminate one function for each G Pi for i = 1,, N V. With renumbering the functions, we end up with N E N V discrete divergence free functions: {Λ 1,, Λ NE N V }. Lemma 4.5. Functions {Λ 1,, Λ NE N V } are linearly independent. Proof. The proof of the lemma is similar to the proofs of Lemma 3.3 and Lemma 4.4. Theorem 4.6. Let D h be defined in (2.8). Then for three dimensional space, D h is spanned by the following basis functions, (4.4) D h = Span{Φ 1,, Φ 12NK, Υ 1,1, Υ 1,2,, Υ NF,1, Υ NF,2, Λ 1,, Λ NE N V }. Proof. The number of the functions in the right hand side of (4.4) is 12N K +2N F + N E N V which is equal to dim(d h ) due to (4.3). Similar to the proof of Theorem 3.4, we can prove that {Φ 1,, Φ 6NK, Υ 1,1, Υ 1,2,, Υ NF,1, Υ NF,2, Λ 1,, Λ NE N V } are linear independent. 5. Numerical Experiments. In this section, we shall report several results of numerical examples for two dimensional Stokes equations. The divergence-free finite element scheme introduced in Algorithm 2 is used. The main purpose if to numerically validate the accuracy and efficiency of the WG scheme.
9 Let v h D h and q h W h, the error for the WG-FEM solution is measured in three norms defined as follows: v h 2 : = ( ) w v h 2 dt + +h 1 T (v 0 v b ) 2 ds, (A discrete H 1 -norm), v 0 2 : = T T T v 0 2 dx, (Element-based L 2 -norm). 9 Table 5.1 Test Case 1: Numerical error and convergence rates for the Stokes equation with homogeneous boundary conditions on the uniform rectangular meshes. h u h Q h u order u 0 Q 0 u order 1/ e e-01 1/ e e e / e e e / e e e / e e e / e e e Table 5.2 Test Case 1: Numerical error and convergence rates for the Stokes equation with homogeneous boundary conditions on the mixed polygonal meshes. h u h Q h u order u 0 Q 0 u order e e e e e e e e e e e e e e e e e e e e e e e Test case 1. The domain is set as Ω = (0, 1) (0, 1). Let the exact solution u and p as follows, ( ) 10x u = 2 y(x 1) 2 (2y 1)(y 1) 10xy 2 (2x 1)(x 1)(y 1) 2 and p = 10(2x 1)(2y 1). It is easy to check that homogeneous Dirichlet boundary condition is satisfied for this testing. The right hand side function f is given to match the exact solutions. The first test shall be performed on the uniform rectangular meshes and the mixed polygonal meshes. The uniform rectangular meshes are generated by partition the domain Ω into n n sub-rectangles. The mesh size is denoted by h = 1/n. Moreover, the WG divergence free algorithm is also test on the mixed polygonal type meshes. We start with the initial mesh shown as the Figure 5.1, which contains the mixture of triangles and quadrilaterals. The next level of mesh is to refine the previous level of mesh by connecting the mid-point on each edge. The mesh size in this case is also denoted by h. The error profile is reported in Table for the rectangular meshes and mixed polygonal meshes, respectively. Both of the tables show the same convergence rate as the theoretical conclusion, which is O(h) in the H 1 norm and O(h 2 ) in the L 2 norm.
10 Test case 2. The domain is given by Ω = (0, 1) (0, 1). Let the exact solutions u and p as follows, ( ) x(1 x)(1 2y) u =, and p = 2(y x). y(1 y)(1 2x) The Dirichlet boundary condition and the right hand side function is set to match the above exact solutions. It is easy to check that the exact solution u satisfies the non-homogeneous boundary condition. For this testing, the WG divergence free algorithm is perform on the triangular grids. The uniform triangular girds are generated by: (1) partition the domain into n n sub-rectangles; (2) divide each square element into two triangles by the diagonal line with a negative slope. The mesh size is denoted by h = 1/n. For the calculation of the pressure p h, we shall make use of the basis function v V h \D h. This basis function is corresponding to the velocity v b related of the normal direction on each edge. Let v V h \D h, the pressure p h is computed as follows, b(v, p h ) = a(u h, v) (f, v 0 ). Beside testing two norms of the error in velocity, we also measure the L 2 error in pressure. The numerical results in Table 5.3 show an O(h) convergence in the norm for velocity, O(h 2 ) convergence in the L 2 -norm for velocity, and O(h) convergence in the L 2 norm for pressure, which are confirmed by Theorem 2.1. Table 5.3 Test Case 2: Numerical error and convergence rates for the Stokes equation with nonhomogeneous boundary conditions. h u h Q h u u 0 Q 0 u p h p e e e e e e e e e e e e e e e e e e e e e e e e-03 Conv.Rate e e-01 Acknowledgement. We offer our gratitude to professor Eric Lord and professor David Singer for their help on obtaining Equation (4.1).
11 11 REFERENCES [1] B. Cockburn, G. Kanschat, D. Schötzau, and C. Schwab, Local discontinuous Galerkin methods for the Stokes system, SIAM J. Numer. Anal., 40 (2002) [2] M. Crouzeix and P. A. Raviart, Conforming and nonconforming finite element methods for solving the stationary Stokes equations, RAIRO Anal. Numer., 7 (1973) [3] V. Girault and P.A. Raviart, Finite Element Methods for the Navier-Stokes Equations: Theory and Algorithms, Springer-Verlag, Berlin, [4] V. Girault, B. Rivière, and M.F. Wheeler, A discontinuous Gelerkin method with nonconforming domain decomposition for Stokes and Navier-Stokes problems, Math. Comp., 74 (2004) [5] D. Griffiths, Finite element for incompressible flow, Math. Meth. in Appl. Sci., 1 (1979) [6] D. Griffiths, The construction of approximately divergence-free finite element, The Mathematics of Finite Element an Its Applications III, Ed. J.R. Whiteman, Academic Press, [7] D. Griffiths, An approximately divergence-free 9-node velocity element for incompressible flows, Inter. J. Num. Meth. in Fluid, 1 (1981) [8] M. D. Gunzburger, Finite Element Methods for Viscous Incompressible Flows, A Guide to Theory, Practice and Algorithms, Academic, San Diego, [9] K. Gustafson and R. Hartman, Divergence-free basis for finite element schemes in hydrodynamics, SIAM J. Numer. Anal., 20 (1983) [10] K. Gustafson and R. Hartman, Graph theory and fluid dynamics, SIAM J. Alg. Disc. Meth., 6 (1985) [11] O. A. Karakashian and W. N. Jureidini, A nonconforming finite element method for the stationary Navier-Stokes equations, SIAM J. Numerical Analysis, 35 (1998) [12] J. Liu and C. Shu, A high order discontinuous Galerkin method for 2D incompressible flows, J. Comput. Phys., 160 (2000) [13] J. Wang and X. Ye, A weak Galerkin finite element method for second-order elliptic problems, J. Comp. and Appl. Math, 241 (2013) [14] J. Wang and X. Ye, A Weak Galerkin mixed finite element method for second-order elliptic problems, Math. Comp., 83 (2014), [15] J. Wang and X. Ye A Weak Galerkin Finite Element Method for the Stokes Equations, Advances in Computational Mathematics, DOI /s , arxiv: [16] J. Wang and X. Ye, New finite element methods in computational fluid dynamics by H(div) elements, SIAM J. Numer. Anal., 45 (2007) [17] J. Wang, X. Wang and X. Ye Finite Element Methods for the Navier-Stokes equations by H(div) Elements, Journal of Computational Mathematics, 26 (2008) [18] X. Ye and C. Hall, A Discrete divergence free basis for finite element methods, Numerical Algorithms, 16 (1997) [19] X. Ye and C. Hall, The Construction of null basis for a discrete divergence operator, J. Computational and Applied Mathematics, 58 (1995)
12 12 n es+1 e 1 e 2 e s+1 e s P i T j n es n e1 n e2 (a) (b) Fig (a) A 2D hull H Pi ; (b) A hull with hanging node. n 3 n 2 E i f s n 4 f s+1 n 1 T j Fig A 3D Hull S Ei Fig Example 1: Level 1 of mixed polygonal mesh.
AMS 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 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 informationNew class of finite element methods: weak Galerkin methods
New class of finite element methods: weak Galerkin methods Xiu Ye University of Arkansas at Little Rock Second order elliptic equation Consider second order elliptic problem: a u = f, in Ω (1) u = 0, on
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 informationWeak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationA 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 informationSTABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS
STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS RAYTCHO LAZAROV AND XIU YE Abstract. In this paper, we derive two stabilized discontinuous finite element formulations, symmetric
More 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 informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More 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 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 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 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 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 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 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 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 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 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 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 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 informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
More 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 informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationUNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 10, Number 3, Pages 551 570 c 013 Institute for Scientific Computing and Information UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT
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 informationarxiv: v2 [math.na] 1 Aug 2018
SIMPLIFIED WEAK GALERKIN AND NEW FINIE DIFFERENCE SCHEMES FOR HE SOKES EQUAION YUJIE LIU AND JUNPING WANG arxiv:1803.0010v [math.na] 1 Aug 018 Abstract. his article presents a simplified formulation for
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More 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 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 informationMODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 138 145 MODEL OF GROUNDWATER FLOW IN FRACTURED ENVIRONMENT JIŘÍ MARYŠKA, OTTO SEVERÝN AND MARTIN VOHRALÍK Abstract. A stochastic discrete
More 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 informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
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 informationNumerische Mathematik
Numer. Math. (2003) 94: 195 202 Digital Object Identifier (DOI) 10.1007/s002110100308 Numerische Mathematik Some observations on Babuška and Brezzi theories Jinchao Xu, Ludmil Zikatanov Department of Mathematics,
More informationarxiv: v1 [math.na] 14 Apr 2011
A WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND-ORDER ELLIPIC PROBLEMS JUNPING WANG AND XIU YE arxiv:1104.2897v1 [math.na] 14 Apr 2011 Abstract. In this paper, authors shall introduce a finite element method
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 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 informationA mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N)
wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad
More 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 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 informationSuperconvergence and H(div) Projection for Discontinuous Galerkin Methods
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 00:1 6 [Version: 2000/07/27 v1.0] Superconvergence and H(div) Projection for Discontinuous Galerkin Methods Peter Bastian
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 informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationA C 0 -Weak Galerkin Finite Element Method for the Biharmonic Equation
J Sci Comput (2014) 59:473 495 DOI 10.1007/s10915-013-9770-4 A C 0 -Weak Galerkin Finite Element Method for the Biharmonic Equation Lin Mu Junping Wang Xiu Ye Shangyou Zhang Received: 3 December 2012 /
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 informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationarxiv: v1 [math.na] 5 Mar 2013
WEAK GALERKIN FINIE ELEMEN MEHODS FOR HE BIHARMONIC EQUAION ON POLYOPAL MESHES LIN MU, JUNPING WANG, AND XIU YE arxiv:133.97v1 [math.na] 5 Mar 13 Abstract. A new weak Galerkin (WG) finite element method
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 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 informationIt is known that Morley element is not C 0 element and it is divergent for Poisson equation (see [6]). When Morley element is applied to solve problem
Modied Morley Element Method for a ourth Order Elliptic Singular Perturbation Problem Λ Wang Ming LMAM, School of Mathematical Science, Peking University Jinchao u School of Mathematical Science, Peking
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 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 informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationAn Adaptive Mixed Finite Element Method using the Lagrange Multiplier Technique
An Adaptive Mixed Finite Element Method using the Lagrange Multiplier Technique by Michael Gagnon A Project Report Submitted to the Faculty of the WORCESTER POLYTECHNIC INSTITUTE In partial fulfillment
More informationANALYSIS AND CONVERGENCE OF A COVOLUME METHOD FOR THE GENERALIZED STOKES PROBLEM
MATHEMATICS OF COMPUTATION Volume 66, Number 217, January 1997, Pages 85 104 S 0025-5718(97)00792-8 ANALYSIS AND CONVERGENCE OF A COVOLUME METHOD FOR THE GENERALIZED STOKES PROBLEM S. H. CHOU Abstract.
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 informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationHybrid Discontinuous Galerkin methods for incompressible flow problems
Hybrid Discontinuous Galerkin methods incompressible flow problems Christoph Lehrenfeld, Joachim Schöberl MathCCES Computational Mathematics in Engineering Workshop Linz, May 31 - June 1, 2010 Contents
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More 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 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 informationICES REPORT Direct Serendipity Finite Elements on Convex Quadrilaterals
ICES REPORT 17-8 October 017 Direct Serendipity Finite Elements on Convex Quadrilaterals by Todd Arbogast and Zhen Tao The Institute for Computational Engineering and Sciences The University of Texas at
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 informationResearch Article A Stabilized Low Order Finite-Volume Method for the Three-Dimensional Stationary Navier-Stokes Equations
Mathematical Problems in Engineering Volume 212, Article ID 297269, 14 pages doi:1.1155/212/297269 Research Article A Stabilized Low Order Finite-Volume Method for the Three-Dimensional Stationary Navier-Stokes
More informationA Posteriori Error Estimates for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems
DOI 10.1007/s10915-013-9771-3 A Posteriori Error Estimates for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems Long Chen Junping Wang Xiu Ye Received: 29 January 2013 / Revised:
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 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 information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More 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 informationA NUMERICAL APPROXIMATION OF NONFICKIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXX, 1(21, pp. 75 84 Proceedings of Algoritmy 2 75 A NUMERICAL APPROXIMATION OF NONFICIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA R. E. EWING, Y. LIN and J. WANG
More informationIFE for Stokes interface problem
IFE for Stokes interface problem Nabil Chaabane Slimane Adjerid, Tao Lin Virginia Tech SIAM chapter February 4, 24 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 / 2 Problem statement
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 information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
More informationStochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows
Stochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows Ivan otov Department of Mathematics, University of Pittsburgh
More 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 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 informationA Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
More informationNONCONFORMING MIXED ELEMENTS FOR ELASTICITY
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 295 307 c World Scientific Publishing Company NONCONFORMING MIXED ELEMENTS FOR ELASTICITY DOUGLAS N. ARNOLD Institute for Mathematics
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationOn the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2006; 22:741 751 Published online 13 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cnm.846
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 informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationA Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains
A Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains Martin Costabel Abstract Let u be a vector field on a bounded Lipschitz domain in R 3, and let u together with its divergence
More informationA FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem
A FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem Leszek Marcinkowski 1 and Nina Dokeva 2 1 Department of Mathematics, Warsaw University, Banacha 2, 02 097 Warszawa, Poland,
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik im Forschungsverbund Berlin e.v. Preprint ISSN 946 8633 Computational comparison between the Taylor Hood and the conforming Crouzeix Raviart element
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
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 information