32.1 Weak mixed formulation
|
|
- Corey Hudson
- 5 years ago
- Views:
Transcription
1 Part VII, Chapter 32 arcy s equations arcy s equations consist of the following PEs in the domain : 1 σ + p = f, a.e. in, (32.1a) σ = g, a.e. in. (32.1b) Theunknownsaretheprimalvariablepandthedualvariableσ;pisalsocalled the potential and σ the flux. The data are, f, and g; is a second-order tensor if the material is anisotropic and it may depend on x if the material is non-homogeneous. The PEs (32.1), together with appropriate boundary conditions, are often used to model porous media flows, e.g., fluid flows in aquifers and petroleum reservoirs. In this context, σ is the seepage velocity, p the pressure, and the material permeability; (32.1a) is called arcy s law, and (32.1b) expresses mass conservation. In this chapter, we derive well-posed weak formulations for(32.1) with various boundary conditions. Then we study mixed finite element approximations using H(div)-conforming spaces for the dual variable Weak mixed formulation We assume that f L 2 (), g L 2 (), and that takes symmetric values and the eigenvalues of are bounded from below and from above uniformly in by λ and λ, We consider irichlet, Neumann, and mixed irichlet Neumann conditions for (32.1). We will see that contrary to the primal formulation studied in Chapter 21, irichlet conditions on p are enforced weakly, and Neumann conditions on σ are enforced strongly irichlet boundary condition In this section we consider the irichlet condition p = a d on. Let us first proceed formally by assuming that all the functions are smooth enough.
2 422 Chapter 32. arcy s equations Multiplying (32.1a) by a smooth R d -valued test function τ, integrating over, integrating by parts the pressure term, and using the irichlet boundary condition, we infer that (τ ( 1 σ) p τ)dx = τ f dx (τ n)a d ds. Furthermore, multiplying (32.1b) by a smooth scalar-valued test function q and integrating over gives q σdx = qgdx. Since f L 2 () and g L 2 (), the volume integrals make sense if we assume that σ,τ H(div;) = {ς L 2 () ς L 2 ()} and that p,q L 2 ().Togivearigorousmeaningtotheboundaryintegral,letusrecall that the normal component of fields in H(div;) can be defined through the map γ n : H(div;) H 1 2( ), see (B.71). Assuming that a d H 2( ), 1 the boundary integral is understood as γ n (τ),a d, where, denotes the duality pairing between H 1 2( ) and H 2( ). 1 We can now define the linear forms F d (τ) := τ f dx γ n(τ),a d and G d (q) := qgdx and the bilinear forms (recall that the tensor is symmetric) a(ς,τ) := 1 τ ς dx, b(ς,q) := q ς dx, (32.2) leading to the following problem: Find σ V := H(div;) and p Q := L 2 () such that a(σ,τ)+b(τ,p) = F d (τ), τ V, b(σ,q) = G d (q), q Q. (32.3) The above assumptions imply that F d V, G d Q Q, a B(V,V ), and b B(V,Q). Proposition 32.1 (Well-posedness). Assume f L 2 (), g L 2 (), and a d H 1 2( ). Then (32.3) is well-posed. Moreover, the pair (σ,p) satisfies the PEs (32.1) a.e. in, and p H 1 () satisfies the boundary condition γ 0 (p) = a d a.e. on where γ 0 : H 1 () H 1 2( ) is the trace map. Proof. (1) We apply Theorem Set B = : V Q so that ker(b) = {v V v = 0}. We have already seen that the bilinear forms a and b are continuous. Moreover, the inequality a(ς,ς) λ 1 ς 2 L 2 () implies that the bilinear form a is coercive in ker(b) since ς L2 () = ς H(div;) if ς ker(b); as a result, the two conditions in (30.5) hold. Moreover, (30.6) follows from Lemma 32.2 below. Hence, all the required assumptions hold. (2) Let us recover the PEs (32.1) and the boundary condition. Upon testing
3 Part VII. PEs in Mixed Form 423 (32.3)withanarbitraryfunctionτ C0 ()andq = 0,weinferthat 1 σ+ p,τ = τ f dx, where, is the duality product in the distribution sense. Since C0 () is dense in L 2 (), this proves that 1 σ + p = f L 2 (). Since σ,f L 2 (), this implies that p L 2 (), so that p H 1 (). A similar distribution argument shows that σ = g in L 2 (). Since p H 1 (), we can use the integration by parts formula (p τ +τ p)dx = γ n (τ),γ 0 (p). We then infer that γ n (τ),γ 0 (p) a d = 0 for all τ H(div;). The surjectivity of γ n from Theorem B.113 implies that φ,γ 0 (p) a d = 0 for all φ H 1 2( ), which proves that γ 0 (p) = a d in H 2( ). 1 Lemma 32.2 (Surjectivity of divergence). Let be a Lipschitz domain. The operator : H(div;) L 2 () is surjective: there is β > 0 such that q τ dx inf sup β. (32.4) q L 2 () q L2 () τ H(div;) τ H(div;) Proof. Let q L 2 (). Let φ H0() 1 be such that ( φ, ψ) L2 () = (q,ψ) L2 () for all ψ H0(). 1 Setting ς q = φ, we have ς q H(div;), ς q = q, and ς q H(div;) c q L 2 (). Setting β = c 1, (32.4) follows from q τ dx q ςq dx sup τ H(div;) τ H(div;) ς q H(div;) β q L 2 (). Remark 32.3 (Regularity of ς q ). It is known from elliptic regularity that the field ς q from the above proof is in H s () with s ( 1 2,1] (see 21.4). The embedding argument from Theorem B.101 then implies that ς q L t () where 1 t 1 2 s d, i.e., t 2 d d 2s (e.g., t [1, ) for s = 1 and d = 2) Neumann boundary condition We now consider the Neumann boundary condition σ n = a n on. We still look for σ H(div;), we assume a n H 1 2( ) and we interpret the boundary condition as γ n (σ) = a n. Since γ n (σ),1 = σdx, we infer that a n and g must satisfy the compatibility condition a n,1 = gdx. (32.5) It is convenient to lift the Neumann boundary condition, by considering a field σ n H(div;) so that γ n (σ n ) = a n (recall that γ n is surjective onto H 1 2( )),andtomakethechangeofvariableσ 0 = σ σ n sothatσ 0 satisfies a homogeneous Neumann boundary condition. Let us set H 0 (div;) := {ς H(div;) γ n (ς) = 0} = ker(γ n ); this is a Hilbert space equipped with the natural norm. Moreover, only the gradient of the primal variable p now
4 424 Chapter 32. arcy s equations appears in the problem, so that we conventionally enforce pdx = 0 and we set L 2 0() := {q L 2 () qdx = 0}. The weak formulation is as follows: Find σ 0 V := H 0 (div;) and p Q := L 2 0() such that a(σ 0,τ)+b(τ,p) = F n (τ), τ V, (32.6) b(σ 0,q) = G n (q), q Q, with linear forms F n (τ) := τ (f 1 σ n )dx and G n (q) := q(g σ n )dx; the bilinear forms a and b are the same as in (32.3) for irichlet conditions. Proposition 32.4 (Well-posedness). Assume that f L 2 (), g L 2 (), a n H 1 2( ) and that the compatibility condition (32.5) holds. Then (32.6) is well-posed. Moreover, the pair (σ := σ 0 +σ n,p) satisfies the PEs (32.1) a.e. in, and the boundary condition γ n (σ) = a n holds in H 1 2( ). Proof. See Exercise Remark 32.5 (Choice of σ n ). One possibility to define σ n H(div;) is to set σ n = φ where φ H () 1 = {q H 1 () qdx = 0} solves the pure Neumann problem φ rdx = grdx+ a n,γ 0 (r) for all r H () 1 (this problem is well-posed owing to (32.5)). The compatibility condition (32.5) implies that any test function r H 1 () can be considered. Taking first r C0 () yields σ n = φ = g, whence γ n (σ n ) a n,γ 0 (r) = 0 for all r H 1 (). The trace map γ 0 being surjective onto H 1 2( ), we conclude that γ n (σ n ) = a n. Note that this construction of the lifting σ n yields G n = Mixed irichlet Neumann boundary conditions ( ) Let d n be a partition of the boundary with d 0 and n 0. We want to enforce the mixed irichlet Neumann conditions p = a d on d and σ n = a n on n. A rigorous mathematical setting for these conditions entails some subtleties. We have seen in that a Neumann condition on is enforced using the normal trace operator γ n which maps onto the dual space H 1 2( ). Here, the Neumann condition is to be enforced on n only, and this leads us to consider the restriction to n of a linear form in H 1 2( ). Let H ( n ) be composed of the functions θ defined on n whose zero-extension to, say θ, is in H 1 2( ). Then, for all ς H(div;), the restriction γ n (ς) n is defined in H ( n ) by using the duality product γ n (ς) n,θ n := γ n (ς), θ. (32.7) The restriction γ n (ς) d is defined in H ( d ) similarly.
5 Part VII. PEs in Mixed Form 425 We assume that a d H ( d ) and that a n H ( n ). We lift the Neumann condition by proceeding as in Remark We set σ n = φ where φ Hd 1() := {q H1 () γ 0 (q) d = 0} solves the mixed irichlet Neumann problem φ rdx = a n,γ 0 (r) n n for all r Hd 1 (). Observing that σ n = 0 and that for all θ H ( n ) there is r θ Hd 1() such that γ 0 (r θ ) = θ (the zero-extension of θ to ), we infer that γ n (σ n ) n a n,θ n = γ n (σ n ), θ φ r θ = 0, i.e., γ n (σ n ) n = a n. We make the change of variable σ 0 = σ σ n. Let us set H n (div;) := {ς H(div;) γ n (ς) n = 0}; this is a Hilbert space equipped with the natural norm. The weak formulation is as follows: Find σ 0 V := H n (div;) and p Q := L 2 () such that a(σ 0,τ)+b(τ,p) = F dn (τ), τ V, (32.8) b(σ 0,q) = G n (q), q Q, with linear form F dn (τ) := τ (f 1 σ n )dx γ n (τ) d,a d d, and G n, a, b are defined as above. Proposition 32.6 (Well-posedness). Assume f L 2 (), g L 2 (), a d H ( d ), and a n H ( n ). Then (32.8) is well-posed. Moreover, the pair (σ := σ 0 +σ n,p) satisfies the PEs (32.1) a.e. in, the irichlet condition γ 0 (p) d = a d hold a.e. on d, and the Neumann condition γ n (σ) n = a n holds in H ( n). Proof. We just sketch the differences with respect to the above proofs. The surjectivityofthedivergencefollowsbydefiningφ H 1 ()suchthat φ = r, φ d = 0, n φ n = 0 and setting ς r = φ (for mixed conditions, the regularity result holds only for s (0, 1 2 ) (see Jochmann [289]), which still verify regularity implies ς r L t () for all t [1,2 d d d 2s ], note that 2d 2s > 2). To recover the irichlet boundary condition, we first infer that γ n (τ),γ 0 (p) = γ n (τ) d,a d d, τ H n (div;). Let ψ C0 ( d ) and let ψ be the zero-extension of ψ to ; there is τ ψ H n (div;) such that γ n (τ) = ψ. This implies that d ψ(γ 0 (p) a d )ds = 0, which in turn gives γ 0 (p) d = a d since ψ is arbitrary Primal, dual, and dual mixed formulations Recall that p is called the primal variable and σ the dual variable. In this section, we consider alternative formulations where either p or σ is eliminated.
6 426 Chapter 32. arcy s equations We focus on homogeneous irichlet boundary conditions for simplicity. The material extends to other types of (non-homogeneous) boundary conditions. The primal formulation of the PEs (32.1) with the boundary condition p = 0 on is obtained by eliminating σ using that σ = (f p). This leads to the following formulation: { Find p H 1 0 () such that a (p,r) := p rdx = F (r), r H0(), 1 (32.9) withlinearformf (r) := ( r f+rg)dx.thisproblemhasbeenanalyzed in Chapter 21 (notice that F H 1 (), see Remark 21.7). In particular, p is the unique minimizer of the energy functional E (q) := 1 2 a (q,q) F (q) over H0(), 1 see Remark The dual formulation is obtained by eliminating p using divergence-free test functions in arcy s law (observe that p τ dx = 0 if τ is divergencefree since p = 0), and enforcing the mass conservation equation explicitly. This leads to the following formulation: { Find σ H(div;) with σ = g such that a (σ,τ) := σ 1 τ (32.10) dx = F (τ), τ H(div = 0;), with the space H(div = 0;) := {ς H(div;) ς = 0} and the linear form F (τ) := τ f dx. The well-posedness of (32.10) can be established by lifting the divergence constraint using σ g H(div;) such that σ g = g, making the change of variable σ 0 = σ σ g H(div = 0;), and observing that the bilinear form a is coercive in H(div = 0;) equipped with the natural norm. Moreover, defining the complementary energy functional E (ς) := 1 2 (ς f) 1 (ς f)dx and observing that (32.10) amounts to E (σ)(τ) = 0 for all τ H(div = 0;), we infer that σ solving (32.10) is the unique maximizer of E on the affine subspace of H(div;) with divergence equal to g, The dual formulation is seldom used for approximation purposes since it requires to manipulate divergence-free vector fields; an interesting application is presented in We can now relate the primal and dual formulations (32.9) and (32.10) to the mixed formulation (32.3), which in the present context is called the dual mixed formulation. Proposition 32.7 (Equivalence, energy identity). The primal formulation (32.9), the dual formulation (32.10), and the dual mixed formulation (32.3) are equivalent in the sense that the solutions p from (32.9) and (32.3) coincide, the solutions σ from (32.10) and (32.3) coincide, and 1 σ + p = f holds. Moreover, the following energy identity holds: min E (q) = E (p) = E (σ) = max E (ς). (32.11) q H0 1() ς H(div;), ς Proof. See Exercise 32.3.
7 Part VII. PEs in Mixed Form 427 Remark 32.8 (Lagrangian). Proposition implies that the pair (σ, p) solving the dual mixed formulation (32.3) is the unique saddle point of the Lagrangian L(ς,q) = (1 2 ς 1 ς q ς)dx F d (ς) G d (q). Since L(ς,q) = E (ς) 1 2 f f dx + q(g ς)dx, we infer that E (p) = E (σ) = L(σ,p) 1 2 f f dx. Remark 32.9 (Elasticity). The above formalism can be applied to the linear elasticity equations, see Chapter 27. enoting by u the displacement and the stress tensor, arcy s law is replaced by the constitutive equation S: = 1 2 ( u+( u)t ), where S is the (fourth-order) compliance tensor, and the mass conservation equation is replaced by the equilibrium equation = g. Weak mixed formulations are discussed in Gatica [232, 2.4.3] Approximation of the mixed formulation In this section we analyze an H(div)-conforming approximation of the weak formulation (32.3) focusing on irichlet boundary conditions for simplicity iscrete problem and well-posedness Let(T h ) h>0 beashape-regularsequenceofaffinesimplicialmeshes,letp d k (T h) be the H(div)-conforming Raviart Thomas finite element space from of order k 0, and let P b k (T h) be the broken finite element space built using P k polynomials; recall that Pk(T d h ) = {ς h H(div;) ψk(ς d h K ) RT k, K T h }, P b k(t h ) = {q h L 2 () ψ g K (q h K) P k, K T h }, (32.12a) (32.12b) where ψk d is the contravariant Piola transformation and ψg K the pullback by the geometric map T K. Note that ς h K RT k and q h K P k since the mesh is affine. It is possible to work with rectangular meshes using Cartesian Raviart Thomas elements for the dual variable and Q k polynomials for the primal variable. Moreover, other boundary conditions can be treated; for mixed irichlet Neumann conditions, boundary faces are assumed to be located in either d or n, and an H n (div;)-conforming subspace is built by taking fields in P d (T h ) with zero normal component on n. The discrete counterpart of (32.3) is Find σ h V h := Pk d(t h) and p h Q h := Pk b(t h) such that a(σ h,τ h )+b(τ h,p h ) = F(τ h ), τ h V h, (32.13) b(σ h,q h ) = G(q h ), q h Q h. Let ker(b h ) = {ς h V h b(ς h,q h ) = 0, q h Q h } and recall that ker(b) = {ς H(div;) ς = 0}.
8 428 Chapter 32. arcy s equations Lemma (iscrete divergence). ker(b h ) ker(b), and there is a uniform constant β = (β 2 Jh d 2 L(L 2 ;L 2 ) + 1) 1 2 > 0, where Jh d is the L2 - stable commuting projection from 15.3, such that inf sup q h ς h dx β. (32.14) q h Q h ς h V h q h L 2 () ς h H(div;) Proof. Thatker(B h ) ker(b)followsfromlemma13.11which,inparticular, shows that : V h = Pk d(t h) Pk b(t h) = Q h. Let us now prove (32.14). Let q h Q h. Using Lemma 32.2 we know that there is ς qh H(div;) such that ς qh = q h and β ς qh H(div;) q h L2 (). Let us set ςh := J h d(ς q h ). Theorem 15.4 implies that ςh = J h b( ς q h ) = Jh b(q h) = q h (since Jh b leaves Pk b(t h) pointwise invariant). Since ςh L 2 () Jh d L(L 2 ;L 2 ) ς qh L 2 (), we infer that β ςh H(div;) q h L2 (), whence sup ς h V h q h ς h dx ς h H(div;) q h ςh dx ςh H(div;) = q h 2 L 2 () ςh β q h L 2 (). H(div;) Corollary (Well-posedness). (32.13) is well-posed. Proof. We apply Proposition Condition (31.4a) follows from the coercivity of a on ker(b) and ker(b h ) ker(b); condition (31.4b) is just (32.14). Remark (Variant of (32.14)). Adapting the last step of the proof of Lemma using the L 2 -norm for ς h instead of the H(div)-norm, we infer that inf qh Q h sup q h ς h dx ςh V h q h L 2 () ς h β L 2 () L := β J d 2 h 1 L(L 2 ;L 2 ) > β. Remark (Fortin operator). Since : H(div;) L 2 () is surjective, it admits a bounded right inverse which we denote ( ) 1. Then Π h = Jh d ( ) 1 Ih b ( ) is a Fortin operator, see Exercise Error analysis One can derive an error estimate from the general setting of Lemma 31.3, see Exercise However, we can exploit the particular structure of arcy s equations to derive a more specific estimate which distinguishes the approximation of the dual variable and of its divergence, both with constants independent of the inf-sup constant β in (32.14) (this is not the case in Exercise 32.7). Let I b h denote the L2 -orthogonal projector onto P b k (T h), see Theorem (Error estimate). Let (σ,p) and (σ h,p h ) solve (32.3) and (32.13), respectively. Let V h,g := {ς h V h B h ς h = Ih b (g)}. Then, σ σ h L 2 () c 1 inf σ ς h L ς h V 2 (), h,g (σ σ h ) L 2 () = inf σ φ h L 2 (), Φ h Pk b(t h) (32.15a) (32.15b) p p h L 2 () c 3 inf σ ς h L ς h V 2 () +2 inf p q h L h,g q h Q 2 (), (32.15c) h
9 Part VII. PEs in Mixed Form 429 with c 1 = λ λ and c 3 = c1 λ β with β L 2 L := β J d 2 h 1 L(L 2 ;L 2 ). Moreover, if σ H r (), σ H r (), and p H r () with r (0,k +1], then σ σ h L 2 () ch r σ H r (), p p h L2 () ch r ( σ Hr () + p Hr ()). (σ σ h ) L 2 () ch r σ H r (), Proof. (1) We first observe that the following Galerkin orthogonality holds: a(σ σ h,τ h )+b(τ h,p p h ) = 0, τ h V h, (32.16a) b(σ σ h,q h ) = 0, q h Q h. (32.16b) Since B h σ h I b h (g) Q h and (B h σ h I b h (g),q h) L2 () = b(σ h σ,q h ) = 0 for all q h Q h owing to (32.16b), we infer that σ h V h,g. (2) Let ς h V h,g, then σ h ς h ker(b h ) ker(b). Since (32.16a) implies that a(σ σ h,τ h ) = 0 for all τ h ker(b), we infer that λ 1 σ σ h 2 L 2 () a(σ σ h,σ σ h ) = a(σ σ h,σ ς h )+a(σ σ h,ς h σ h ) = a(σ σ h,σ ς h ) λ 1 σ σ h L2 () σ ς h L2 (). Thus, σ σ h L 2 () c 1 σ ς h L 2 (), and (32.15a) holds. Since σ h = B h σ h = I b h (g) and Ib h is an L2 -orthogonal projection, (32.15b) also holds. (3) Let q h Q h. Owing to Remark 32.12, we infer that there is τ h V h such that τ h = p h q h and β L 2 τ h L 2 () p h q h L 2 (). We infer that p h q h L 2 () = b(τ h,p h q h ) = b(τ h,p h p)+b(τ h,p q h ) = a(σ σ h,τ h )+b(τ h,p q h ) λ 1 σ σ h L2 () τ h L2 () + τ h L2 () p q h L2 (), owing to (32.16a). The above properties of τ h combined with (32.15a) and the triangle inequality imply that (32.15c) holds. (4) Finally, if σ H r () with r > 0, then σ L p () with p > 2 owing to the Embedding Theorem B.101. Since σ = g L 2 (), we infer using (10.1b) that σ is in the domain of the Raviart Thomas interpolation operator I d h. Taking ς h = I d h (σ), which is in V h,g since I d h (σ) = Ib h ( σ) = Ib h (g), the convergence rate on σ σ h L 2 () results from Corollary The other rates follow from this one and the approximations properties of Ih b. Remark (Proof variant). If g = 0, we can also take ς h = Jh d(σ) since Jh d(σ) = J h b( σ) = J h b(g) = 0 so that ς h V h,0. This way, there is no need to use the operator Ih d and invoke Theorem B.101. We now use duality techniques to derive an improved-order error estimate on the primal variable. One difference with respect to the primal formulation analyzed in is that we now bound the (discrete) error I b h (p) p h,
10 430 Chapter 32. arcy s equations instead of the full error p p h. As in , we assume that the following smoothing property holds: There is s (0,1] and a uniform constant c sp such that the (adjoint) solution z H 1 0() of the PE ( z) = φ satisfies the a priori bound z H 1+s () c sp φ L 2 (). Sufficient conditions for this smoothing property to hold follow from elliptic regularity, see We also make the (so-called) multiplier assumption on the field stating that the map H s () ξ ξ H s () is bounded, see (21.29). Theorem (Potential supercloseness). Under the above smoothing property and multiplier assumption, the following holds: I b h(p) p h L2 () ch s σ σ h H(div;). (32.17) Proof. Let z H 1+s () H 1 0() be the adjoint solution with data φ = I b h (p) p h and set ξ = z. We observe that I b h(p) p h 2 L 2 () = (p p h,i b h(p) p h ) L2 () = (p p h, ξ) L2 () = (p p h, J d h(ξ)) L 2 () = ( 1 (σ σ h ),J d h(ξ)) L 2 () = ( 1 (σ σ h ),J d h(ξ) ξ) L2 () +( 1 (σ σ h ),ξ) L2 (), where we have used (p I b h (p),q h) L2 () = 0 for all q h P b k (T h) and the definition of the adjoint solution in the first line, and J d h (ξ) = J b h ( ξ) = ξ since ξ P b k (T h) and (32.16a) in the second line. The first term on the right-hand side, say T 1, is bounded as T 1 λ 1 σ σ h L 2 () J d h(ξ) ξ L 2 () c σ σ h L 2 ()h s ξ H s (), and ξ H s () c z H 1+s () owing to the multiplier assumption. For the second term, say T 2, (32.16b) implies that (σ σ h, z) L2 () = ( (σ σ h ),z) L2 () = ( (σ σ h ),z I b h(z)) L2 (). This leads to the bound T 2 c (σ σ h ) L 2 ()h z H 1 (), and we conclude using the smoothing property z H 1+s () c sp I b h (p) p h L 2 (). Remark (L -bounds). L -bounds on the dual and primal errors are derived, e.g., by Gastaldi and Nochetto [231] Hybridization ( ) Hybridization was introduced by Fraejis de Veubeke in 1965 (see [225] for a reprint) as a computationally effective technique that allows one to transform the original (symmetric indefinite) linear system on the primal and dual variables into a symmetric positive definite system on an auxiliary variable that
11 Part VII. PEs in Mixed Form 431 can be interpreted as the Lagrange multiplier associated with the continuity constraint on the normal component of the dual variable. Let us focus on simplicial Raviart Thomas elements so that the discrete spaces are V h = Pk d(t h) and Q h = Pk b(t h). Let Λ h := Pk b(f h ) be spanned by piecewise polynomials of order at most k on each mesh interface. Recall the broken Raviart Thomas space V hy h := P d,b k (T h ) spanned by piecewise RT k polynomials in each mesh cell. The hybridized version of (32.13) is Find σ h V hy h, p h Q h, and λ h Λ h such that a(σ h,τ h )+b(τ h,p h )+c(τ h,λ h ) = F(τ h ), τ h V hy h, (32.18) b(σ h,q h ) = G(q h ), q h Q h, c(σ h,µ h ) = 0, µ h Λ h, with bilinear form c(τ h,µ h ) := F Fh F [[τ h]] n F µ h ds. The third equation in (32.18) implies that σ h has continuous normal component, i.e., σ h V h. Restricting test functions in the first equation to V h then shows that the pair (σ h,p h ) coincides with the unique solution of (32.13). The crucial advantage of (32.18) is that (σ h,p h ) can be eliminated locally. This operation is called static condensation. To state our result, we need some notation. For all µ h Λ h and all K T h, we define µ K = (µ h F ) F FK, where µ h is extended by zero on boundary faces; thus, µ K is a piecewise polynomial of order at most k on the boundary of K. Furthermore, we write ψ = (ψ σ,ψ p ) V k := RT k P k, and we define the local bilinear form c K (ψ,ζ) := ( 1 ψ σ,ζ σ ) L2 (K) (ψ p, ζ σ ) L2 (K) (ζ p, ψ σ ) L2 (K). Let Φ K : L 2 (K) L 2 (K) V k be such that Φ K (f,g) V k solves c K (Φ K (f,g),ψ) = (f,ψ σ ) L 2 (K) (g,ψ p ) L 2 (K) for all ψ V k. Let Θ K : L 2 ( K) V k be such that Θ K (µ) V k solves c K (Θ K (µ),ψ) = (µ,ψ σ n K ) L2 ( K) for all ψ V k. Proposition (Static condensation). The pair (σ h,p h ) solves (32.13) if and only if (σ h K,p h K ) = Φ K (f K,g K )+Θ K (λ K ) where λ h Λ h (extended by zero on boundary faces) is the unique solution of 1 ΘK(λ σ K ),ΘK(µ σ K )) L 2 (K) = l(µ h ), µ h Λ h, (32.19) K T h ( with l(µ h ) := K T h (f,θk σ(µ K)) L2 (K) (g,θ p K (µ K)) L2 (K). The algebraic realization of this problem leads to a symmetric positive definite matrix. Plugging this formula into the third equation of (32.18) leads to a linear problem in the Lagrange multipliers whose algebraic realization is a symmetric positive definite matrix; more details can be found in Boffi et al. [60, ] and Cockburn and Gopalakrishnan [150]. Remark (One unknown per cell). Implementation of lowest-order Raviart Thomas methods with one unknown per mesh cell and relations to other methods (including Finite Volumes and Mimetic Finite ifferences) are discussed in Younes et al. [469], Vohralík and Wohlmuth [454].
12 432 Chapter 32. arcy s equations (i) Show that Φ K and Θ K are defined by solving well-posed (local) irichlet problems in mixed form. Exercises Exercise 32.1 (Neumann condition). Prove Proposition (Hint: for the surjectivity of the divergence, pose a pure Neumann problem.) Exercise 32.2 (Integration by parts). Let Hd 1() and H n(div;) be defined in Prove that ( q ς+q ς)dx = 0 for all q H1 d () and all ς H n (div;). (Hint: observe that γ 0 (q) n H ( n ).) Exercise 32.3 (Primal, dual formulations). Prove Proposition Exercise 32.4 (Primal mixed formulation). Consider the problem of finding p H 1 () such that p = f and γ 0 (p) = g with f L 2 () and g H 1 2( ). erive a mixed formulation of this problem with unknowns (p,λ) H 1 () H 1 2( ) and show that it is well-posed; this method is introduced in Babuška [30]. (Hint: set b(v,µ) = µ,γ 0 (v) and observe that B = γ 0 : H 1 () H 1 2( ).) Recover the PE and the boundary condition. Exercise 32.5 (Fortin operator). Justify Remark (Hint: use arguments similar to those of the proof of Lemma ) Exercise 32.6 (Inf-sup condition). The goal is to prove the inf-sup condition (32.14) using the canonical Raviart Thomas interpolation operator. (i) o this by using elliptic regularity. (Hint: pose a irichlet problem.) (ii) o this by using the surjectivity of : H 1 () L 2 (). Exercise 32.7 (Error estimate). Prove that σ σ h H(div;) c 1 inf ς h V h σ ς h H(div;), p p h L2 () c 3 inf ς h V h σ ς h H(div;) +2 inf q h Q h p q h L2 (), no 1+ with c 1 = (1+ λ λ )(1+ 1 β ) and c 3 = c 1 λ β. Assuming that σ H r (), σ L H r (), and p H r () with r (0,k + 2 1], prove that σ σ h H(div;) ch r ( σ Hr ()+ σ Hr ()) and p p h L2 () ch r ( σ Hr ()+ σ Hr ()+ p H r ()). (Hint: use the commuting projection Jh d.)
13 Part VII. PEs in Mixed Form 433 Solution to exercises Exercise 32.1 (Neumann condition). There are only two differences with the proof of Proposition The first one concerns the divergence operator which now maps from H 0(div;) to L 2 0 () (owing to the divergence theorem). To prove that this operator is surjective, let q L2 0 () and let us solve the pure Neumann problem φ H 1 () such that φ = q and φdx = 0. Then ς q = φ has all the desired properties. The second difference concerns the recovery of the boundary condition which is now a simple consequence of γ n(σ) = γ n(σ 0) + γ n(σ n) = a n. Exercise 32.2 (Integration by parts). Observe that ( q ς+q ς)dx = γ n(ς),γ0(q). 1 Now γ 0(q) n H 2 00 ( n) since its zero-extension to is γ0(q) (since q H1 d ()) which is in H 1 2( ). This implies that γ n(ς),γ 0(q) = γ n(ς) n,γ 0(q) n n, and this last quantity vanishes since ς H n(div;). Exercise 32.3 (Primal, dual formulations). Let p H 1 0 () solve (32.9) and define σ = 1 (f p). Since ( p τ + p τ)dx = 0 for all τ H(div;), the first equation of (32.3) holds; moreover, we infer that σ qdx = gqdx for all q H1 0 (), implying that σ = g so that the second equation of (32.3) holds. Let now (σ,p) H(div;) L 2 () solve (32.3). Then σ = g and taking a divergence-free test function τ in the first equation of (32.3) shows that (32.10) holds. Finally, let σ H(div;) with σ = g solve (32.10); let p L 2 () be defined such that p τ dx = ( 1 σ f) τ dx for all τinh(div;) (note that the righthand side vanishes if τ is divergence-free). Since : H(div;) L 2 () is surjective owing to Lemma 32.2, p is well-defined, and its definition implies that the first equation of (32.3) holds; the second one follows from σ = g. Finally, since p = f 1 σ, the energy identity results from E (p) = 1 p pdx = 1 (ς 1 (ς f) f) = E (σ). 2 2 Exercise 32.4 (Primal mixed formulation). The weak mixed formulation is Find p V := H 1 () and λ Q := H 1 2( ) such that a(p,q) + b(q,λ) = f(q), q V, b(p,µ) = g(µ), µ Q, with a(p,q) = p qdx, b(v,µ) = µ,γ0(v), f(q) = fqdx, and g(µ) = µ,g. All these forms are bounded. Moreover, by reflexivity of H 1 2( ), B = γ 0 so that ker(b) = ker(γ 0) = H 1 0 () owing to Theorem B.106, whence we infer that the bilinear form a is coercive in ker(b). Furthermore, B is surjective still by Theorem B.106, so that well-posedness follows from Theorem Finally, the second equation in the weak mixed formulation implies that µ,γ 0(p) g = 0 for all µ H 1 2( ), so that we recover the boundary condition γ 0(p) = g a.e. on, and the PE follows by taking q arbitrary in H 1 0 () in the first equation. Exercise 32.5 (Fortin operator). Let us verify the two properties stated in Lemma 18.9(i) with W = H(div;), W h = P d k (T h), V = L 2 (), and V h = P b k (T h). The operator Π h maps from H(div;) to P d k (T h) as required. Moreover, we infer that, for all q h P b k (T h) and all v H(div;), q h Π h (v)dx = q h (J d h (( ) 1 (I b h ( v))))dx = q h J i h (Ib h ( v))dx = q h I b h ( v)dx = q h vdx, since J d h = J i h ( ), J i h leaves Pb k (T h) pointwise invariant, and q h is in P b k (T h). Moreover, using the stability of all the operators, one can see that Π h (v) H(div;) c v L 2 (). Exercise 32.6 (Inf-sup condition). (i) Let q h Q h. Let φ H 1 0 () solve φ = q h and set ς qh = φ. Elliptic regularity implies that ς qh H s () with s > 1 2, i.e., ςq h ˇV d (), see (13.14c) with p = 2. This means that ς qh is in the domain of canonical interpolation operator I d h. Moreover, Id K (ς q h K) L 2 (K) c( ς qh K L 2 (K) +h s K ς q h K H s (K) ) for all K T h. Summing over mesh cells and since h K is bounded by the diameter of, we infer that I d h (ςq h ) L 2 () c ςq h H s () c q L 2 (). We can now conclude as in the proof of Lemma 32.10, but this time we take ς h = Id h (ςq h ).
14 434 Chapter 32. arcy s equations (ii) Let again q h Q h. We use the hint to infer that there is ς qh H 1 () such that ς qh = q h and ς qh H 1 () c q L 2 (). Since I d K (ς q h K) L 2 (K) c( ς qh K L 2 (K) + h 1 K ς q h K H 1 (K) ) for all K T h, we can now conclude as in Step (i). Exercise 32.7 (Error estimate). The error bound on the dual variable follows from Lemma 31.3 since a = λ, α = λ, and b = 1. For the primal variable, we proceed as in the proof of Theorem and use the bound on the dual variable. To bound the best-approximation error on the dual variable, we use the commuting projection J d h and Theorem 15.4 to infer that σ J d h (σ) H(div;) σ J d h (σ) L 2 () + σ J b h ( σ) L 2 () c inf σ ς h ς h P k d L 2 () + inf σ φ h (T h ) φ h P k b L 2 (), (T h ) and we conclude invoking Corollary 14.8 (with p = 2 and x = d,b).
Mixed 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 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 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 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 informationTraces and Duality Lemma
Traces and Duality Lemma Recall the duality lemma with H / ( ) := γ 0 (H ()) defined as the trace space of H () endowed with minimal extension norm; i.e., for w H / ( ) L ( ), w H / ( ) = min{ ŵ H () ŵ
More informationON LEAST-SQUARES FINITE ELEMENT METHODS FOR THE POISSON EQUATION AND THEIR CONNECTION TO THE DIRICHLET AND KELVIN PRINCIPLES
ON LEAST-SQUARES FINITE ELEMENT METHODS FOR THE POISSON EQUATION AND THEIR CONNECTION TO THE DIRICHLET AND KELVIN PRINCIPLES PAVEL BOCHEV 1,3 AND MAX GUNZBURGER 2 Abstract. Least-squares finite element
More informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More 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 informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More 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 informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More 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 informationCoupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements
Coupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements V.J. Ervin E.W. Jenkins S. Sun Department of Mathematical Sciences Clemson University, Clemson, SC, 29634-0975, USA. Abstract We study
More informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More 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 informationQUADRILATERAL H(DIV) FINITE ELEMENTS
QUADRILATERAL H(DIV) FINITE ELEMENTS DOUGLAS N. ARNOLD, DANIELE BOFFI, AND RICHARD S. FALK Abstract. We consider the approximation properties of quadrilateral finite element spaces of vector fields defined
More informationA 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 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 informationDiscontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints
J A N U A R Y 0 1 8 P R E P R N T 4 7 4 Discontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints gor Voulis * and Arnold Reusken nstitut für Geometrie und Praktische
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 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 informationON LEAST-SQUARES FINITE ELEMENT METHODS FOR THE POISSON EQUATION AND THEIR CONNECTION TO THE DIRICHLET AND KELVIN PRINCIPLES
SIAM J. NUMER. ANAL. Vol. 43, No. 1, pp. 340 362 c 2005 Society for Industrial and Applied Mathematics ON LEAST-SQUARES FINITE ELEMENT METHODS FOR THE POISSON EQUATION AND THEIR CONNECTION TO THE DIRICHLET
More informationWeak Formulation of Elliptic BVP s
Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed
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 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 informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
More informationNumerische Mathematik
Numer. Math. (2012) 122:61 99 DOI 10.1007/s00211-012-0456-x Numerische Mathematik C 0 elements for generalized indefinite Maxwell equations Huoyuan Duan Ping Lin Roger C. E. Tan Received: 31 July 2010
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationA DECOMPOSITION RESULT FOR BIHARMONIC PROBLEMS AND THE HELLAN-HERRMANN-JOHNSON METHOD
Electronic ransactions on Numerical Analysis. Volume 45, pp. 257 282, 2016. Copyright c 2016,. ISSN 1068 9613. ENA A DECOMPOSIION RESUL FOR BIHARMONIC PROBLEMS AND HE HELLAN-HERRMANN-JOHNSON MEHOD WOLFGANG
More informationApproximation in Banach Spaces by Galerkin Methods
2 Approximation in Banach Spaces by Galerkin Methods In this chapter, we consider an abstract linear problem which serves as a generic model for engineering applications. Our first goal is to specify the
More informationA multipoint flux mixed finite element method on hexahedra
A multipoint flux mixed finite element method on hexahedra Ross Ingram Mary F. Wheeler Ivan Yotov Abstract We develop a mixed finite element method for elliptic problems on hexahedral grids that reduces
More informationPart VIII, Chapter 38. Residual-based stabilization Model problem
Part VIII, Chapter 38 Residual-based stabilization This chapter, as the next one, is concerned with the approximation of riedrichs systemsusingh 1 -conformingfiniteelementsinastandardgalerkin setting.
More informationUniform inf-sup condition for the Brinkman problem in highly heterogeneous media
Uniform inf-sup condition for the Brinkman problem in highly heterogeneous media Raytcho Lazarov & Aziz Takhirov Texas A&M May 3-4, 2016 R. Lazarov & A.T. (Texas A&M) Brinkman May 3-4, 2016 1 / 30 Outline
More informationA FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS. 1. Introduction
A FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS CHRISTOPHER HARDER, DIEGO PAREDES 2, AND FRÉDÉRIC VALENTIN Abstract. We aim at proposing novel
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More 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 informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationTime domain boundary elements for dynamic contact problems
Time domain boundary elements for dynamic contact problems Heiko Gimperlein (joint with F. Meyer 3, C. Özdemir 4, D. Stark, E. P. Stephan 4 ) : Heriot Watt University, Edinburgh, UK 2: Universität Paderborn,
More informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
More information1 Constrained Optimization
1 Constrained Optimization Let B be an M N matrix, a linear operator from the space IR N to IR M with adjoint B T. For (column vectors) x IR N and y IR M we have x B T y = Bx y. This vanishes for all y
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 informationICES REPORT A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams
ICS RPORT 15-17 July 2015 A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams by Omar Al-Hinai, Mary F. Wheeler, Ivan Yotov The Institute for Computational ngineering
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 informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
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 informationControllability of the linear 1D wave equation with inner moving for
Controllability of the linear D wave equation with inner moving forces ARNAUD MÜNCH Université Blaise Pascal - Clermont-Ferrand - France Toulouse, May 7, 4 joint work with CARLOS CASTRO (Madrid) and NICOLAE
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More 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 informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationMIXED FINITE ELEMENT METHODS FOR LINEAR VISCOELASTICITY USING WEAK SYMMETRY
MIXED FINITE ELEMENT METHODS FOR LINEAR VISCOELASTICITY USING WEAK SYMMETRY MARIE E. ROGNES AND RAGNAR WINTHER Abstract. Small deformations of a viscoelastic body are considered through the linear Maxwell
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 SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY
A SECOND ELASTICITY ELEMENT USING THE MATRIX BUBBLE WITH TIGHTENED STRESS SYMMETRY J. GOPALAKRISHNAN AND J. GUZMÁN Abstract. We presented a family of finite elements that use a polynomial space augmented
More 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 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 informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
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 informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationFinite Element Clifford Algebra: A New Toolkit for Evolution Problems
Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/
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 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 informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationCSE386C METHODS OF APPLIED MATHEMATICS Fall 2014, Final Exam, 9:00-noon, Tue, Dec 16, ACES 6.304
CSE386C METHODS OF APPLIED MATHEMATICS Fall 2014, Final Exam, 9:00-noon, Tue, Dec 16, ACES 6.304 1. Let A be a closed operator defined on a dense subspace D(A) of a Hilbert space U. (a) Define the adjoint
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 informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
More 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 information37.1 Basic ideas. (37.1a) Symmetry: A k = (A k ) T for all k {1:d}, a.e. in D, (37.1b) Positivity: µ 0 > 0 s.t. K+K T X 2µ 0 I m a.e. in D. (37.
Part VIII, Chapter 37 Friedrichs systems Symmetric positive systems of first-order linear PEs were introduced in 1958 by Friedrichs[232] to study transonic flows. Friedrichs wanted to handle within a single
More informationMixed Hybrid Finite Element Method: an introduction
Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola
More informationu = f in Ω, u = q on Γ. (1.2)
ERROR ANALYSIS FOR A FINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET BOUNDARY CONTROL PROBLEMS S. MAY, R. RANNACHER, AND B. VEXLER Abstract. We consider the Galerkin finite element approximation of
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
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 informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
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 informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
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 informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
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 information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationOverview of normed linear spaces
20 Chapter 2 Overview of normed linear spaces Starting from this chapter, we begin examining linear spaces with at least one extra structure (topology or geometry). We assume linearity; this is a natural
More informationMimetic Finite Difference methods
Mimetic Finite Difference methods An introduction Andrea Cangiani Università di Roma La Sapienza Seminario di Modellistica Differenziale Numerica 2 dicembre 2008 Andrea Cangiani (IAC CNR) mimetic finite
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More 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 informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More 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 informationMATH 263: PROBLEM SET 1: BUNDLES, SHEAVES AND HODGE THEORY
MATH 263: PROBLEM SET 1: BUNDLES, SHEAVES AND HODGE THEORY 0.1. Vector Bundles and Connection 1-forms. Let E X be a complex vector bundle of rank r over a smooth manifold. Recall the following abstract
More informationNew Helmholtz-Weyl decomposition in L r and its applications to the mathematical fluid mechanics
New Helmholtz-Weyl decomposition in L r and its applications to the mathematical fluid mechanics Hideo Kozono Mathematical Institute Tohoku University Sendai 980-8578 Japan Taku Yanagisawa Department of
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 informationAPPROXIMATION OF THE EIGENVALUE PROBLEM FOR THE TIME HARMONIC MAXWELL SYSTEM BY CONTINUOUS LAGRANGE FINITE ELEMENTS
APPROXIMATION OF THE EIGENVALUE PROBLEM FOR THE TIME HARMONIC MAXWELL SYSTEM BY CONTINUOUS LAGRANGE FINITE ELEMENTS ANDREA BONITO, AND JEAN-LUC GUERMOND,3 Abstract. We propose and analyze an approximation
More informationLectures on. Weak Solutions of Elliptic Boundary Value Problems
Lectures on Weak Solutions of Elliptic Boundary Value Problems S. Kesavan The Institute of Mathematical Sciences, CIT Campus, Taramani, Chennai - 600 113. e mail: kesh@imsc.res.in 2 1 Some abstract variational
More informationOutline. 1 Motivating examples. 2 Hilbert complexes and their discretization. 3 Finite element differential forms. 4 The elasticity complex 1 / 25
Outline 1 Motivating examples Hoge theory, Hilbert complexes, an finite element ifferential forms Douglas N. Arnol 2 Hilbert complexes an their iscretization 3 Finite element ifferential forms 4 The elasticity
More information10. Smooth Varieties. 82 Andreas Gathmann
82 Andreas Gathmann 10. Smooth Varieties Let a be a point on a variety X. In the last chapter we have introduced the tangent cone C a X as a way to study X locally around a (see Construction 9.20). It
More informationA Second Elasticity Element Using the Matrix Bubble
Portland State University PDXScholar Mathematics and Statistics Faculty Publications and Presentations Fariborz Maseeh Department of Mathematics and Statistics 2011 A Second Elasticity Element Using the
More information