A SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS
|
|
- Bertina Warren
- 6 years ago
- Views:
Transcription
1 A SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS CONSTANTIN BACUTA AND KLAJDI QIRKO Abstract. We investigate new PDE discretization approaces for solving variational formulations wit different types of trial and test spaces. Te general mixed formulation we consider assumes a stability LBB condition and a data compatibility condition at te continuous level. We expand on te Bramble-Pasciak s least square formulation for solving suc problems by providing new ways to coose approximation spaces and new iterative processes to solve te discrete formulations. Our proposed metod as te advantage tat a discrete inf sup condition is automatically satisfied by natural coices of test spaces (first) and corresponding trial spaces (second). In addition, for te proposed iterative solver, a nodal basis for te trial space is not required. Applications of te new approac include discretization of first order systems of PDEs, suc as div curl systems and time-amonic Maxwell equations. 1. Introduction Over te last two decades, tere ave been many advances in applying finite element least squares metods to approximate first order systems of PDEs, [6, 7, 8, 14, 15, 16, 17, 18, 26, 27, 28]. However, wen compared to te more establised field of finite element metods for elliptic problems, a unified teoretical framework for least squares approximation of solutions of first order systems of PDEs is missing. Our proposed framework provides powerful preconditioning tecniques, as efficient error estimators, is suitable to multilevel tecniques, and leads to robust and easy to implement solvers. We combine known teory and discretization tecniques for approximating elliptic problems and for symmetric saddle point problems (see [1, 5, 11, 12, 13, 23, 29, 31, 32, 36]) to obtain a unified framework for discretizing variational formulations wit different types of test and trial spaces. In particular, te framework can be applied to least square approximation for a large class of first order systems of PDEs. For te applications we consider, te solution spaces are L 2 type spaces, and te data can reside in weak negative norm spaces. We require tat te test spaces be H 1 type spaces wit suitable boundary conditions, and 2000 Matematics Subject Classification. 74S05, 74B05, 65N22, 65N55. Key words and prases. least squares, saddle point systems, mixed metods, multilevel metods, Uzawa type algoritms, conjugate gradient, cascadic algoritm, dual DPG. Te work was supported by NSF, DMS
2 2 CONSTANTIN BACUTA AND KLAJDI QIRKO te discrete test spaces be conforming finite element spaces built using te action of te continuous differential operator associated wit a given problem. Among te advantages of te metod are te following: te discretization leads to saddle point variational formulation wit automatic discrete inf sup condition, and assembly of stiffness matrices for te trial spaces is avoided. Te general abstract problem tat we plan to discretize using a Saddle Point Least Squares Metod (SPLS) is: Find p Q suc tat (1.1) b(v, p) = f, v, for all v V or B p = f, were V and Q are infinite dimensional Hilbert spaces and b(, ) is a continuous bilinear form on V Q, tat satisfies a standard inf sup condition, and f V belongs to te range of B. In te special case wen te operator B associated wit te form b(, ) is injective, our metod can be viewed as a conforming Petrov-Galerkin metod. From te point of view of coosing te discrete spaces, our metod can be caracterized as te dual of te Discontinuous Petrov-Galerkin (DPG) metod introduced by Demkowicz and Gopalakrisnan in [20, 21]. Wile bot metods ave strong connections wit te least squares and minimum residual tecniques, te proposed discretization process stands apart from te DPG approac because of te opposite order and different ways in wic te trial and test spaces are cosen. We propose te following main steps of our saddle point least square discretization metod: Step 1 Reduce te general problem (1.1) to a saddle point least square formulation, using te natural inner product a 0 (, ) on V V as te (1, 1) form of te saddle point system (see problem (2.12)). Step 2 Coose a standard conforming approximation space V for te variational space V. Step 3 Construct a discrete trial space M Q using te operator B associated wit te form b(, ). For example, take M := C 1 BV, or M := Q C 1 BV, were C 1 is te Riesz representation operator for te space Q and Q is a projection from Q to a subspace M. Te pair (V, M ) will automatically satisfy a discrete inf sup condition. Step 4 Write te discrete version of SPLS formulation on (V, M ), see (3.3), and replace a 0 (, ) by an equivalent form a prec (, ) on V V. Step 5 Solve te new discrete SPLS problem using an Uzawa type iterative process tat requires only te action of te preconditioner associated wit a prec (, ) and te action of C 1 B or Q C 1 B. Te paper is organized as follows. In Section 2, we introduce notation and review basic abstract results needed to describe te metod. We also include ere te first step of reduction of a mixed problem to a Saddle Point Problem (SPP). In Section 3, we present te discretization part of te metod (Step 2 and Step 3) and discuss te coice of discrete spaces
3 SADDLE POINT LEAST SQUARES 3 and teir approximability. Uzawa type Iterative solvers witout trial space bases are presented in Section 4. Te special coice of spaces wit discrete inf sup condition and teir approximability properties are presented in Section 5. In Section 6 we present an example of SPLS discretization for a div curl system. Te Appendix (Section 8) contains some important functional analysis results needed for te proofs of te paper. 2. Notation and Background In tis section, we start wit a review of te notation of te classical SPP teory and introduce te spaces, te operators and te norms for te general abstract case. We let V and Q be two Hilbert spaces wit inner products a 0 (, ) and (, ) respectively, wit te corresponding induced norms V = = a 0 (, ) 1/2 and Q = = (, ) 1/2. Te dual pairings on V V and Q Q are denoted by,. Here, V and Q denote te duals of V and Q, respectively. Wit te inner products a 0 (, ) and (, ), we associate te operators A : V V and C : Q Q defined by and Au, v = a 0 (u, v) for all u, v V Cp, q = (p, q) for all p, q Q. Te operators A 1 : V V and C 1 : Q Q are called te Rieszcanonical isometries and satisfy te following propperties (2.1) a 0 (A 1 u, v) = u, v, A 1 u V = u V, u V, v V, (2.2) (C 1 p, q) = p, q, C 1 p = p Q, p Q, q Q. Next, we suppose tat b(, ) is a continuous bilinear form on V Q, satisfying te inf-sup condition. More precisely, we assume tat (2.3) inf p Q sup v V and (2.4) sup p Q sup v V b(v, p) p v = m > 0 b(v, p) p v = M <. Here, and trougout tis paper, te inf and sup are taken over nonzero vectors. Wit te form b, we associate te linear operators B : V Q and B : Q V defined by Bv, q = b(v, q) = B q, v for all v V, q Q. Let V 0 be te kernel of B or C 1 B, i.e., V 0 = Ker(B) = {v V Bv = 0} = {v V C 1 Bv = 0}. Due to (2.4), V 0 is a closed subspace of V.
4 4 CONSTANTIN BACUTA AND KLAJDI QIRKO We notice ere tat C 1 B : V Q and A 1 B : Q V are dual to eac oter in te Hilbert sense, and consequently, te Scur complement S 0 := C 1 BA 1 B is a symmetric operator. In addition, S 0 is a positive definite operator wit σ 0 (S 0 ) [m 2, M 2 ], see Lemma 8.1 or [2] Te general problem. Assume tat b : V Q R is a bilinear form satisfying (2.3) and (2.4), and let f V be given. Many variational formulations of PDEs, including first order systems, can be written in te mixed form: Find p Q suc tat (2.5) b(v, p) = f, v, for all v V. Te operator form of (2.5) is to solve te following equation for p (2.6) B p = f. Te existence and te uniqueness of (2.5) was first studied by Aziz and Babuška in [1] and is known as te Babuška s Lemma. Lemma 2.1. (Babuška) Let b : V Q R be a bilinear form satisfying (2.3) and (2.4), and let f V. Te problem: Find p Q suc tat (2.7) b(v, p) = f, v, for all v V as a unique solution if and only if (2.8) f, v = 0, for all v V 0. If (2.8) olds, and p is te solution of (2.7), ten (2.9) m p p S0 := (S 0 p, p) 1/2 = F V M p. A proof of te Lemma can be found in [1, 2]. To justify our saddle point least squares terminology, we write (2.6) in te equivalent form (2.10) A 1 B p = A 1 f. Using tat te Hilbert transpose of A 1 B : V Q is te operator C 1 B : V Q, te normal equation for solving (2.10) in te least square sense is: (2.11) C 1 BA 1 B p = C 1 BA 1 f. Since S 0 := C 1 BA 1 B is a symmetric positive definite operator, te problem (2.11) as a unique solution p. Ten, we consider te well posed saddle point problem: Find (u, p) (V, Q) suc tat (2.12) a 0 (u, v) + b(v, p) = f, v for all v V, b(u, q) = 0 for all q Q. One can immediately ceck tat p is te unique solution of (2.11) if and only if (u, p) is te unique solution of (2.12), were u = A 1 (B p F ).
5 SADDLE POINT LEAST SQUARES 5 Proposition 2.2. In te presence of te continuous inf sup condition (2.3) and te compatibility condition (2.8), we ave tat p is te unique solution of (2.5) if and only if (u = 0, p) is te unique solution of (2.12). Remark 2.3. Even in te absence of te compatibility condition (2.8), but assuming te continuous inf sup condition (2.3), we ave tat (2.12) as unique solution regardless of te fact tat (2.5) is wellposed or not. Te above arguments justify te name of te system (2.12) as te SPLS variational formulation of (1.1) A note on approximating te solution by Uzawa type algoritms. Given a relaxation parameter α (0, 2/M 2 ), te Uzawa algoritm for approximating te solution (u = 0, p) of (2.12) can be described as follows. Start wit any p 0 Q. For j = 1, 2,... Compute u j V, q j, p j Q by a 0 (u j, v) = f, v b(v, p j 1 ), v V q j = C 1 Bu j, p j = p j 1 + α q j Te convergence of te algoritm is discussed in many publications, see e.g., [5, 13, 22, 23, 30]. It is easy to ceck tat (2.13) p p j = (I αc 1 BA 1 B )(p p j 1 ) = (I αs 0 )(p p j 1 ). By using te spectral bounds for S 0, we obtain a convergence rate γ := I αs 0 = max{ 1 αm 2, 1 αm 2 } < 1, provided tat α ( ) 2 0, M. 2 It is easy to ceck tat for j 1, we also ave (2.14) u j+1 = (I αa 1 B C 1 B)u j. Due to te compatibility assumption (2.8), eac residual representation u j in te Uzawa algoritm, belongs to V0. Using te spectral properties of S 1 := A 1 B C 1 B : V0 V0 (see Lemma (8.1) (v)) and (2.14), we ave tat te series u j is absolutely convergent. Consequently, if we j=1 start te algoritm wit p 0 = 0, ten p k = αc 1 B( k u j ) and p k p = C 1 B w wit w = α j=1 u j V0. Tus, te Uzawa algoritm not only approximates te solution of (2.7), but also finds a special representation of te solution: p = C 1 B w wit w V0. Similar remarks old for Uzawa Gradient (UG) and Uzawa Conjugate Gradient (UCG) algoritms. 3. SPLS discretization In ligt of Propsition 2.2, we are led to te corresponding SPLS discretization of (1.1). We let V V and M Q be finite dimensional j=1
6 6 CONSTANTIN BACUTA AND KLAJDI QIRKO approximation spaces and consider te restrictions of te forms a 0 (, ) and b(, ) to te discrete spaces V and M. Assume tat te following discrete inf sup condition olds for te pair (V, M ). (3.1) inf p M b(v, p ) sup v V p v = m > 0. We define te corresponding discrete operators A, C, B, and B. For example, A is te discrete version of A, and is defined by A u, v = a 0 (u, v ), for all u V, v V. We define V,0 to be te kernel of B V,0 := {v V b(v, q ) = 0, for all q M }, V,1 := V,0 = A 1 B (M ), and S,0 := C 1 B A 1 B. Remark 3.1. If V,0 V 0, ten te compatibility condition (2.8) implies a discrete compatibility condition. Consequently, under te discrete stability assumption (3.1), te problem: Find p M suc tat as unique solution. b(v, p ) = f, v, for all v V, In general, V,0 is not contained in V 0. Tus, te compatibility condition (2.8) does not imply a discrete compatibility condition, and te problem: Find p M suc tat (3.2) b(v, p ) = f, v, for all v V, or B p = f, migt not be well posed. Neverteless, under te assumption (3.1), te following saddle point discrete variational problem : Find (u, p ) V M suc tat (3.3) a 0 (u, v ) + b(v, p ) = f, v for all v V, b(u, q ) = 0 for all q M, always as a unique solution. Using te corresponding discrete operators, te problem (3.3) is equivalent to: Find (u, p ) V M suc tat u (3.4) = A 1 (f B T p ) C 1 B A 1 B p = C 1 B A 1 f. Since te operator S,0 = C 1 B A 1 B is a symmetric positive definite on M, it is invertible, and te second equation of (3.4) as te unique solution p. Tus, one way to solve te system (3.3), would be to solve te second equation of (3.4) for p, and ten find u from te first equation of (3.4). Since C 1 B and A 1 B are dual to eac oter in te Hilbert sense and te problem (3.2) is also equivalent wit A 1 B p = A 1 f, we ave tat te solution of te second equation of (3.4) is in fact te least squares solution of (3.2). Te component u of te solution of (3.3) or (3.4) becomes te
7 SADDLE POINT LEAST SQUARES 7 representation on V of te residual associated wit te least square solution of (3.2). Anoter reason for wic te saddle point least squares discretrization of (1.1) is te variational formulation (3.3) is tat (3.3) is te natural discretization of te SPLS formulation (2.12) of (1.1). Using tat (3.3) is te discrete variational formulation of (2.12), and based on te classical error analysis for SPP teory, we can find standard estimates for 0 u and p p, see [9, 11, 32, 36]. If we assume discrete stability, ten te second component p of te solution of (3.3) provides a good approximation to te least square solution p of (1.1) even in te absence of te compatibility condition (2.8). In te case wen te compatibility condition (2.8) olds, we are dealing wit a special saddle point problem, and a sarp error estimate for p p can be proved using te Xu-Zikatanov argument in [36]. Teorem 3.2. Let b : V Q R satisfy (2.3) and (2.4), and assume tat f V is given and satisfies (2.8). Assume tat p is te solution of (1.1) and V V, M Q are cosen suc tat te discrete inf sup condition (3.1) olds. Ten, if (u, p ) is te solution of (3.3), te following error estimate olds: (3.5) 1 M u p p M m inf p q. q M Te proof of te estimate can be found in te Appendix. Te rigt part of te estimate (3.5) is an improvement of te similar estimate presented in [10] tat provides a bound tat is linear in m 2. Tis improvement is significant if m depends on. Remark 3.3. All te considerations made so far in tis section make sense if te form a 0 (, ), as an inner product on V, is replaced by anoter inner product wic gives rise to an independent of equivalent norm on V. Certainly, te definition of A, S,0, and m, will cange accordingly wit te new inner product, but te error estimate (3.5) remains valid wit different estimating constants tat factor in te norm equivalence constants. Tis observation leads to an efficient preconditioning approac for SPLS discretization. More precisely, if te Uzawa algoritm is involved to solve (3.3), ten te action of A 1 can be replaced by te action of any equivalent preconditioner. From te first equation of (3.3) and (3.5) we get tat te solution (u, p ) satisfies: u = f B p V M 2 inf p q, q M were f is te functional v f, v on V. Tis gives an estimate for te residual associated wit te discrete equation (3.2). Remark 3.4. Te Bramble-Pasciak least squares metod in [10] for discretizing (2.7), focuses on solving only te variational formulation of te m
8 8 CONSTANTIN BACUTA AND KLAJDI QIRKO second equation of (3.4), i.e., Find p M suc tat (3.6) (C 1 B A 1 B p, q ) = f, A 1 B q, for all q M, wic can be reformulated as: Find p M suc tat (3.7) b(a 1 B q, p ) = f, A 1 B q, for all q M, Te approac in [10] to numerically solve (3.7) is based on using families of stable pairs {(V, M )}, bases for bot V and M, and on replacing te action of A 1 on V by an approximate inverse. Our approac is based on iteratively solving te coupled saddle point system (3.3) wose p component of te solution is te least squares solution of (3.2) and on using bases and matrix assembly only for te test space V. Remark 3.5. Te DPG metod, [20, 21], for discretizing (2.7), also reduces to a saddle point formulation tat is similar to (3.3). Te details are presented in Teorem 2.4 of [20] and a proof can be found in [24]. For te DPG approac, te trial spaces Q or M are sougt as a product between an interior component and and an interface compomonent at bot continuous and discrete levels. Once a finite dimensional subspace V r of V (wit certain properties) is cosen, te discrete test space, using our notation, is V = A 1 r B M, were A 1 r is te Riesz representation operator for te Hilbert space (V r, a 0 (, )). In our approac, we coose te discrete trial space V V first, a finite dimensional subspace M Q, and ten define M := Q C 1 BV, were C 1 is te Riesz representation operator for te space Q and Q is a projection from Q to M, (see Section 5 for details). Given te connection between approximating te solution of te mixed variational problem (1.1) and solving te discrete problem (3.3), one can just rely on te available tecniques and solvers for symmetric saddle point problems of type (2.12) on stable families of spaces {(V, M )} in order to solve (1.1). Neverteless, if it is difficult to find families of stable pairs or tey exist but are difficult to implement, ten we present next alternative ways of approximating te solution of (1.1). 4. Iterative solvers witout basis for te trial space Te computational callenge we face wen solving (3.3) on (V, M ) is tat one migt not be able to find stable pairs (V, M ). Even in te case tat stable pairs are available, a global linear system corresponding to (3.3) migt be difficult to assemble and solve. It is possible to solve (3.3) witout even aving explicit bases for M by using te Uzawa (U), Uzawa Gradient (UG), or Uzawa Conjugate Gradient (UCG) algoritms. Following [3], we ave tat te standard U and UG algoritms can be rewritten suc tat tey differ only by te way te parameter α is cosen. For te Uzawa ( ) algoritm, 2 we ave to coose a fixed number α = α 0 in te interval 0,. For te M 2 UG algoritm, te parameter α is cosen to impose te ortogonality of
9 SADDLE POINT LEAST SQUARES 9 consecutive residuals associated wit te second equation in (3.3). Te first step for Uzawa is identical wit te first step of UG. We combine te two algoritms in: Algoritm 4.1. (U-UG) Algoritms Step 1: Set u 0 = 0 V, p 0 M, compute u 1 V, q 1 M by a 0 (u 1, v) = f, v b(v, p 0 ), for all v V (q 1, q) = b(u 1, q), for all q M. Step 2 : For j = 1, 2,..., compute j, α j, p j, u j+1, q j+1 by (U UG1) a 0 ( j, v) = b(v, q j ), v V (Uα) α j =α 0 for te Uzawa algoritm or (UGα) α j = (q j, q j ) for te UG algoritm b( j, q j ) (U UG2) p j =p j 1 + α j q j (U UG3) u j+1 =u j + α j j (U UG4) (q j+1, q) = b(u j+1, q), for all q M. Here, f is te restriction of f to V. To obtain te UCG algoritm, te UG algoritm is modified as in [9, 35] as follows: First, we define d 1 := q 1 in Step 1, and ten modify Step 2 by replacing b(, q j ) wit b(, d j ), were {d j } is a sequence of conjugate directions: Algoritm 4.2. (UCG) Algoritm Step 1: Set u 0 = 0 V, p 0 M. Compute u 1 V, q 1, d 1 M by a 0 (u 1, v) = f, v b(v, p 0 ), v V (q 1, q) = b(u 1, q), for all q M, d 1 := q 1. Step 2 For j = 1, 2,..., compute j, α j, p j, u j+1, q j+1, β j, d j+1 by (UCG1) a 0 ( j, v) = b(v, d j ), v V (UCGα) α j = (q j, q j ) b( j, q j ) (UCG2) p j =p j 1 + α j d j (UCG3) u j+1 =u j + α j j (UCG4) (q j+1, q) = b(u j+1, q), for all q M (UCGβ) β j = (q j+1, q j+1 ) (q j, q j ) (UCG6) d j+1 =q j+1 + β j d j
10 10 CONSTANTIN BACUTA AND KLAJDI QIRKO Eac one of te described algoritms converges to te discrete solution (u, p ) of (3.3). In addition, te following sarp error estimation result was proved by one of te autors in a sligtly more general context in [3]. Teorem 4.3. If (u, p ) is te discrete solution of (3.3), and (u j+1, p j ) is te j t iteration for U, UG, or UCG, ten (u j+1, p j ) (u, p ) and (4.1) 1 M 2 q j+1 p j p 1 m 2 q j+1, m M 2 q j+1 u j+1 u M m 2 q j+1. Besides convergence of te iteration processes, te result entitles q j+1 as a computable, robust, efficient, and uniform-modulo m estimator for all tree algoritms. Under stability presence, we can use te Teorem 3.2 and te estimates (4.1) to build adaptive or multilevel algoritms for SPLS discretization. A cascadic algoritm for solving symmetric SPPs was introduced by one of te autors in [3]. One major advantage of solving te system (3.3) as te SPLS discretization of (1.1), using one of tese tree algoritms is tat (4.1) remains valid, if te rigt correction of constant factors is made, wen a 0 (, ) is replaced by a uniform equivalent form a prec (, ). Similar to te continuous case in Section 2.2, an important observation ere is tat if te starting initial guess is p 0 = 0, ten te p j - iterates remain in te space C 1 B V and tat {p j } approximates and represent te solution p in te form C 1 B w, wit w V,0. In addition, as presented in te next section, for certain spaces M a basis is not needed for solving te variational problems associated wit q j, and at eac step of te U-type iterative processes, only te action of A 1 or a preconditioner requires an inversion process. 5. Construction of special discrete spaces Next, we will present a metod to address te lack of stability of te approximation spaces. Let V be a subset of V, and let M be a finite dimensional subspace of Q tat as good approximability properties. Typical examples of spaces M are te spaces of piecewise polynomials. We equip M wit an inner product tat could differ from te restriction of te Q inner product on M, but induces an equivalent norm (independent of ). For convenience, we denote te inner product on M by (, ). If Q : Q M is te ortogonal projection onto M, we simply define te space M by M := Q C 1 BV. We consider te restriction of te form a 0 (, ) to V V and te restriction of b(, ) to V M, and define te discrete operators A, C, B, and B
11 SADDLE POINT LEAST SQUARES 11 for te pair (V, M ). For any q M, v V, we ave b(v, q ) = B v, q = (C 1 B v, q ). On te oter and, since q M Q, and v V V, we ave b(v, q ) = (C 1 Bv, q ) = ( Q C 1 Bv, q ). From te above identities, we get (5.1) C 1 B v = Q C 1 Bv for all v V. Tis implies tat C 1 B is onto M and, using tat V and M are finite dimensional spaces, te discrete inf-sup condition (3.1) olds for some m > 0 tat migt depend on. Tus, te problem (3.3) as unique solution (u, p ), and any of te tree Uzawa type algoritms presented in Section 4 can be applied to approximate it. In addition, in ligt of 5.1, we ave tat te residual q j for eac Uzawa type algoritm satisfies q j = C 1 B u j = Q C 1 Bu j. Consequently, te computation of q j involves te computation Bu j, te action of C 1 at te continuous level (often just a multiplication operator) to find C 1 Bu j, and a projection onto te space M - tat is usually a standard finite element space. Tus, for eac of te Uzawa type algoritms, a basis for M is not needed for solving te variational problems associated wit q j Approximability of te space M. Due to Teorem 3.2, in order to expect small discretization error p p, besides stability, one needs to investigate te minimization problem inf p q in te special case q M wen M is a proper subspace of M. Te continuous inf sup condition (2.3) tat we assume, guaranties tat A 1 B is injective and as closed range, see Lemma (8.1) (iii). Tus, te dual operator C 1 B is onto Q, and for any p Q, tere exists u V suc tat C 1 Bu = p. For p = C 1 Bu Q, and any q = Q C 1 Bv M, we ave tat (5.2) p q = C 1 Bu Q C 1 Bv C 1 Bu C 1 Bv + Q C 1 Bv C 1 Bv M u v + Q C 1 Bv C 1 Bv. In order to get good SPLS approximation properties for te solution p of (1.1), it would be enoug to ask for some regularity of a solution u of C 1 Bu = p, for approximation properties of V, and for approximation properties of te projection Q : C 1 BV M. Note tat te argument remains valid regardless of te discrete stability presence.
12 12 CONSTANTIN BACUTA AND KLAJDI QIRKO 5.2. Te special projection case. If Q C 1 Bv = 0 implies C 1 Bv = 0, for any v V, (i.e., Q is injective on C 1 BV ), ten V,0 V 0 and, see Remark 3, te variational formulation (3.2) is well posed, as unique solution p M, and using Proposition 2.2 for te discrete pair (V, M ), we ave tat (u = 0, p ) is te solution of (3.3). Tis will always be te case if we coose M = C 1 BV, and Q to be te identity operator on C 1 BV. In tis special case we do ave M = M = C 1 BV and, if p is te solution of (1.1) and p is te solution of (3.2), from (1.1) and (3.2) we obtain tat 0 = b(v, p p ) = (C 1 Bv, p p ), for all v V. Tus, we simply ave tat p is te ortogonal projection of p onto M, and consequently, using te representation p = C 1 Bw for some w V, we ave (5.3) p p = inf p q = inf C 1 Bw C 1 Bv q M v V M inf w v. v V Compared wit (3.5), te discretization error estimate (5.3) as te advantage tat is independent of te stability constant m and, in addition, it reduces to approximability of functions in V by discrete functions in V. If one of te Uzawa type iterative metods of Section 4 is applied to find te solution (0, p ) of (3.3) ten, according to Teorem (4.3), te iteration error satifies p j p 1 m 2 q j+1. If te discretization error order is available, say p p = O( α ), and an estimate about m is also available, te iteration error can matc te discretization error by imposing te stopping criteria (5.4) q j+1 c 0 m 2 α, were c 0 is a constant independent of. Consequently, using one of te Uzawa algoritms, we can approximate te solution p by an iterate p j C 1 BV up to optimal discretization error order, regardless of (discrete B-) stability presence. From te numerical experiments we performed so far for concrete SPLS discretization problems, projecting on smoot spaces M could lead to better approximation of te continuous solution p, and projecting on coarser spaces MH could lead to stability. It is wort noticing ere, tat if one cooses M suc tat Q q c q, for all q C 1 BV, ten stability of ( V, C 1 ) ( BV implies stability for te pair V, Q C 1 BV ). A typical coice for te space M used to define te projection Q is te space of continuous piecewise polynomials of degree m wit respect
13 SADDLE POINT LEAST SQUARES 13 to a mes T. On M we can coose te standard inner product or we can consider te lumping inner product introduced in [4]. To define te lumping inner product, assume tat {φ 1, φ 2,..., φ m } is a nodal basis for M. We can define (, ) l by (5.5) (φ i, φ j ) l := δ j i (1, φ i), i, j = 1, 2,..., m, and extend it to M M by m m (5.6) α i φ i, β j φ j := i=1 j=1 l m α i β i (1, φ i ). i=1 Using te identity (5.1), see [4], te computation of C 1 m m C 1 B v = i=1 (C 1 Bv, φ i ) φ i = (1, φ i ) i=1 B v becomes b(v, φ i ) (1, φ i ) φ i. We note ere tat by using te lumping inner product, one avoids mass matrix inversion at eac iterative step. All of te above provides a more detailed picture of te five step SPLS discretization summarized at te end of introduction. A concrete example of ow to apply te metod is presented next. 6. SPLS discretization of a div curl system Here, we apply te SPLS discretization metod for a model div curl problem on a polyedral domain Ω R 3. For a given data, we are looking to find te vector function L 2 (Ω) suc tat (6.1) = j in Ω (µ) = g in Ω (µ) n = σ on Γ := Ω, were µ is a given parametric scalar L 2 function tat is strictly positive on Ω. Te variational formulation we adopt for (6.1) is similar to te approac of [10]. We multiply te first equation in (6.1) by w H 1 0 (Ω), multiply te second equation by ϕ H 1 (Ω)/R and after we integrate by parts, we obtain (6.2) (, w) = (j, w) for all w H 1 0 (Ω) (, µ ϕ) = G, ϕ := ( g, ϕ) + (σ, ϕ) Γ for all ϕ H 1 (Ω)/R, were (, ) denotes te standard L 2 type inner product. If we define V := H 1 0 (Ω) H1 (Ω)/R, Q := L 2 (Ω), and te form b(v, p) on (V, Q), by b((w, ϕ), ) := (, w + ϕ), for all (w, ϕ) V, Q, te variational formulation for (6.1) becomes: Find Q suc tat (6.3) b((w, ϕ), ) = F, (w, ϕ) := (j, w) + G, ϕ, for all (w, ϕ) V. Te coice of inner products on te spaces V and Q is essential in building robust solvers wit respect to te function µ. By coosing te weigted
14 14 CONSTANTIN BACUTA AND KLAJDI QIRKO (, ) µ inner product on Q = L 2 (Ω) and te weigted inner product on V := H 1 0 (Ω) H1 (Ω)/R induced by te norm (w, ϕ) 2 V := a µ 1(w, w) + a µ (ϕ, ϕ) := µ 1 w 2 + µ ϕ 2, Ω Ω it can be proved tat, in tis case, te form b(, ) satisfies a continuous inf sup condition wit a constant independent of µ, see [10]. Te C 1 B operator tat appears in te SPLS discretization is (6.4) C 1 B (w, ϕ) = µ 1 curl w + ϕ. By using a Helmoltz decomposition, it is easy to ceck tat C 1 B is onto Q = L 2 (Ω). If te data (j, g, σ) is suc tat te compatibility condition (2.8) is satisfied, ten (6.2) is a well-posed problem and our SPLS discretization can be applied. For SPLS discretization, we can cose V V to be te standard vector space of continuous piecewise functions of degree m wit te appropriate boundary conditions for eac component of V. For M we consider two coices: Case 1) M := C 1 BV, and Case 2) M := Q (C 1 BV ), were Q is te ortogonal projection onto M - te continuous piecewise polynomials of degree m. On M we can consider te standard inner product or te lumping inner product defined in (5.6). In bot cases, we automatically ave a positive discrete inf sup constant m. Sufficient conditions tat can establis stability for tis problem remain to be investigated for various coices of discrete spaces. For te 2D version of (6.1), it is easy to ceck tat we do ave (free) stability for bot types of discretization. Tis is simply due to te fact tat te (curl, curl) form is equivalent wit te Diriclet form (, ) on te space H0 1(Ω). In general, establising stability for te family {(V, M = C 1 BV )} migt not be an easy problem, even if we focus on V to be a standard conforming space of piecewise polynomial of degree k. For example, if we consider te simple case of te operator B being te div : H 1 0 (Ω) L2 (Ω) and C 1 to be just te identity operator, ten te problem of discrete B-stability is equivalent wit establising te classical Stokes stability for (V, div V ). From te works of Scott and Vogelius on divergence stability, [33, 34], for te case of general quasi-uniform meses, te divergence stability remains a 30 year old open problem. For many problems of interest, te B-stability can be establised based on paricular properties of te operator B. Preliminary numerical results for div curl systems and te Maxwell equations, suggest te presence of B-stability for particular coices of test spaces and projection operators. A more detailed investigation of B stability for te 3D div curl or Maxwell problem remains to be investigated in te near future Numerical results. We performed numerical experiments for approximating (6.1) wit Uzawa Conjugate Gradient (UCG) and V te conforming P 1 elements. For M we consider bot cases described above in tis
15 SADDLE POINT LEAST SQUARES 15 section. In Case 2) we project on M te vector space of continuous P 1 functions, wit te componentwise inner product given by (5.6). Te domain Ω was cosen te unit cube, and te function µ was cosen to be te restriction to Ω of { 1, if x < 1 µ = 2 µ 0, if x 1 2, for various values of te constant µ 0. First we consider te case wen µ 0 = 1 (constant µ) and data cosen suc tat te exact solution is = (x(1 x)yz, xy(1 y)z, xyz(1 z)). For Case 1), using a power metod for te discrete Scur complement S,0 we estimated m and imposed te stopping criterion given by (5.4) wit c 0 = 0.4 and α = 1. For Case 2) we just imposed a fixed number of iterations on eac level. Te numerical results are presented below in Table 1 and Table 2. Level µ = 1 P 1 C 1 BP 1 comp Order # of it. m Table 1: UCG for SPLS-P 1 -discretization Case 1) Level µ = 1 P 1 Q C 1 B(P 1 ) comp Order # of it Table 2: UCG for SPLS-P 1 -discretization Case 2) Next, we consider te case of discontinuous µ (µ 0 = 50) and data cosen suc tat te exact solution is = (x(1 x)( 1 2 x)yz, x(1 2 x)y(1 y)z, x(1 x)yz(1 z)). 2 Te numerical results are presented below in Table 3 and Table 4. For Case 1), wit µ aving 2D ceckerboard jump discontinuities (1 on wite squares, µ 0 on te black squares and constant 1 in te z- direction), we obtained similar results. Te estimates for te discrete inf sup
16 16 CONSTANTIN BACUTA AND KLAJDI QIRKO Level µ 0 = 50 P 1 C 1 BP 1 comp Order # of it. m Table 3: UCG for SPLS-P 1 -discretization Case 1) Level µ 0 = 50 P 1 Q C 1 B(P 1 ) comp Order # of it Table 4: UCG for SPLS-P 1 -discretization Case 2) constant m are presented in Table 5. Wile we can notice instability wit respect to (or Level) we can also notice numerical stability wit respect to µ 0. Level µ Table 5: Estimates for m for different values, µ 0 Remark 6.1. We note tat for Case 1), te discretization error satisfies (5.3) and for te conforming coice P 1 for V we ave = O(). Ten, since comp + comp, and te iteration error comp matces te discretization error (see Teorem 4.3) we also expect O() for comp. Tis is numerically observed in Table 1 and Table 3. For Case 2), te discretization error satisfies an estimate of type (5.2). In tis case, our estimate for te discretization error is not optimal, because te functions C 1 Bv for v V are piece wise constants, and te most we can expect from Q C 1 Bv C 1 Bv is O( 1/2 ). Neverteless, te computation we performed in tis case, (see Table 2 and Table 4), sow
17 SADDLE POINT LEAST SQUARES 17 tat te comp is at least O(). Te supper convergence beavior in tis case remains to be teoretically and numerically investigated. 7. Conclusions We presented a saddle point least squares metod to discretize variational formulations wit different types of trial and test spaces. To te autors knowledge, te proposed SPLS approac is different from te DPG formulations as presented in [20, 21, 19], were a trial space is cosen firstly, and a (close to optimal) test space tat provides stability of te pairs is cosen secondly. In te first order systems case, te essential differences between te proposed SPLS discretization and te classical least squares finite element metod approac, as described in [6, 7, 8, 16, 17] or te FOSL & FOSLL approaces in [14, 15, 18, 26, 27, 28], are tat for SPLS discretization te test spaces and te trial spaces are of a different nature. Tis makes SPLS discretization more general, but requires special attention in finding stable spaces for discretization or in finding efficient iterative solvers. Tanks: Te autors would like to tank to Jay Gopalakrsnan for valuable discussions tat relate DPG and te SPLS and to Francisco Sayas for valuable software and programming support for obtaining te numerical solvers. 8. Appendix 8.1. Functional analysis results. For a bounded linear operator T : X Y between two Hilbert spaces X and Y, we denote by T t te Hilbert transpose of T. If X = Y, we say tat T is symmetric if T = T t. For a bounded linear operator T : X X, we denote te spectrum of te operator T by σ 0 (T ). Te following lemma provides important properties of norms and operators to be used in tis paper. A proof of it can be found in [2]. Lemma 8.1. Let A, C, B, and B be te operators associated wit te spaces V, Q and te connecting form b(, ). Assume tat (2.3) and (2.4) are satisfied. i) Te operators C 1 B : V Q and A 1 B : Q V are symmetric to eac oter, i.e., (8.1) (C 1 Bv, q) = a 0 (v, A 1 B q), v V, q Q, consequently, (C 1 B) t = A 1 B and (A 1 B ) t = C 1 B. ii) Te Scur complement on Q is te operator S 0 := C 1 BA 1 B : Q Q. Te operator S 0 is symmetric and positive definite on Q, satisfying (8.2) (S 0 p, p) = sup v V b(v, p) 2 v 2.
18 18 CONSTANTIN BACUTA AND KLAJDI QIRKO Consequently, m 2, M 2 σ 0 (S 0 ) and (8.3) σ 0 (S 0 ) [m 2, M 2 ]. iii) Te following estimate olds (8.4) p S0 := (S 0 p, p) 1/2 = A 1 B p V m p for all p Q. Consequently, A 1 B : Q V as closed range, V 1 := A 1 B (Q) is a closed subspace of V and A 1 B : Q V 1 is an isomorpism. iv) Te Scur complement on V is defined as te operator S := A 1 B C 1 B : V V. Te operator S is symmetric and non-negative definite on V, wit Ker(S) = V 0, S(V) = V 1, and satisfies (8.5) a 0 (Su, v) = (C 1 Bu, C 1 Bv), u, v V. v) Te Scur complement on V 1 = V 0 is te restriction of S to V 1, i.e., S 1 := A 1 B C 1 B : V 1 V 1. Te operator S 1 is symmetric and positive definite on V 1, satisfying (8.6) σ 0 (S 1 ) = σ 0 (S 0 ) [m 2, M 2 ] Proof of Teorem 3.2. Proof. Let p be te solution of (1.1) and assume tat(u, p ) is te solution of (3.3). First, we notice tat te operator T : Q Q defined by T p = p is linear and idempotent. To prove tat T is idempotent we consider te problem: Find (ũ, p ) V M suc tat (8.7) a 0 (ũ, v ) + b(v, p ) = b(v, p ) for all v V, b(ũ, q ) = 0 for all q M. Due to te assumption (3.1), we get tat (8.7) as unique solution and by noticing tat (ũ, p ) = (0, p ) solves te problem, we can conclude tat T p = p and consequently, T 2 = T. According to Kato [25], and Xu and Zikatanov [36], we ave Next, for any q M we ave (8.8) I T L(Q,Q) = T L(Q,Q). p p = (I T )p = (I T )(p q ) I T p q To estimate T we use (3.1): = T p q. T p 1 b(v, T p) sup = 1 b(v, p ) sup. m v V v m v V,1 v Solving for b(v, p ) from te first equation of (3.3), and using te fact tat f, v = b(v, p) for all v V, we furter get
19 SADDLE POINT LEAST SQUARES 19 (8.9) T p 1 m M m p. b(v, p) a 0 (u, v ) sup = 1 b(v, p) sup v V,1 v m v V,1 v Tus, T M m, and from (8.8) we obtain te rigt side of (3.5). To prove te left side of (3.5), from (1.1) and te first equation of (3.3) we get Ten, b(v, p p ) = a 0 (u, v ), for all v V. a 0 (u, v ) b(v, p p ) u = sup = sup M p p, v V v v V v wic concludes te validity of te estimate (3.5). References [1] A. Aziz and I. Babuška. Survey lectures on matematical foundations of te finite element metod. Te Matematical Foundations of te Finite Element Metod wit Applications to Partial Differential Equations, A. Aziz, editor, [2] C. Bacuta. Scur complements on Hilbert spaces and saddle point systems. J. Comput. Appl. Mat., 225(2): , [3] C Bacuta. Cascadic multilevel algoritms for symmetric saddle point systems. Comput. Mat. Appl., 67(10): , [4] C. Bacuta, B. McCracken, and L. Su. Residual reduction algoritms for nonsymmetric saddle point problems. J. Comput. Appl. Mat., 235(6): , [5] M. Benzi, G. Golub, and J. Liesen. Numerical solutions of saddle point problems. Acta Numerica, pages 1 137, [6] P. Bocev, Z. Cai, T. A. Manteuffel, and S. F. McCormick. Analysis of velocity-flux first-order system least-squares principles for te Navier-Stokes equations. I. SIAM J. Numer. Anal., 35(3): , [7] P. Bocev and M.D. Gunzburger. Least-squares finite element metods. In International Congress of Matematicians. Vol. III, pages Eur. Mat. Soc., Züric, [8] P. Bocev and M.D. Gunzburger. Least-squares finite element metods, volume 166 of Applied Matematical Sciences. Springer, New York, [9] D. Braess. Finite Elements. Teory, Fast Solvers, and Applications in Solid Mecanics. Cambridge University Press, Cambridge, [10] J. H. Bramble and J. Pasciak. A new approximation tecnique for div-curl systems. Mat. Comp., 73: , [11] S. Brenner and L.R. Scott. Te Matematical Teory of Finite Element Metods. Springer-Verlag, New York, [12] F. Brezzi. On te existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers. Rev. Française Automat. Informat. Recerce Opérationnelle Sér. Rouge, 8(R-2): , [13] F. Brezzi and M. Fortin. Mixed and Hybrid Finite Element Metods. Springer-Verlag, New York, [14] Z. Cai, C.-O. Lee, T. A. Manteuffel, and S. F. McCormick. First-order system least squares for te Stokes and linear elasticity equations: furter results. Iterative metods for solving systems of algebraic equations.
20 20 CONSTANTIN BACUTA AND KLAJDI QIRKO [15] Z. Cai, C.-O. Lee, T. A. Manteuffel, and S. F. McCormick. First-order system least squares for linear elasticity: numerical results. SIAM J. Sci. Comput., 21(5): , Iterative metods for solving systems of algebraic equations (Copper Mountain, CO, 1998). [16] Z. Cai, T. A. Manteuffel, and S. F. McCormick. First-order system least squares for te Stokes equations, wit application to linear elasticity. SIAM J. Numer. Anal., 34(5): , [17] Z. Cai, T. A. Manteuffel, S. F. McCormick, and S. V. Parter. First-order system least squares (FOSLS) for planar linear elasticity: pure traction problem. SIAM J. Numer. Anal., 35(1): (electronic), [18] Z. Cai, T. A. Manteuffel, S. F. McCormick, and J. Ruge. First-order system LL (FOSLL ): scalar elliptic partial differential equations. SIAM J. Numer. Anal., 39(4): , [19] W. Damen, C. Huang, C. Scwab, and G. Welper. Adaptive Petrov-Galerkin metods for first order transport equations. SIAM J. Numer. Anal., 50(5): , [20] L. Demkowicz and J. Gopalakrisnan. A class of discontinuous Petrov-Galerkin metods. Part I: te transport equation. Comput. Metods Appl. Mec. Engrg., 199(23-24): , [21] L. Demkowicz and J. Gopalakrisnan. A primal DPG metod witout a first-order reformulation. Comput. Mat. Appl., 66(6): , [22] M. Fortin and R. Glowinski. Augmented Lagrangian Metods Applications to te numerical solutions of boundary value problems, volume 15. Nort-Holland, [23] V. Girault and P.A. Raviart. Finite Element Metods for Navier-Stokes Equations, volume 15. Springer-Verlag, Berlin, [24] Jay Gopalakrisnan. Five lectures on dpg metods. arxiv , [25] T. Kato. Estimation of iterated matrices, wit application to te Von Neumann condition. Numer. Mat., 2:22 29, [26] B. Lee, T. A. Manteuffel, S. F. McCormick, and J. Ruge. First-order system leastsquares for te Helmoltz equation. SIAM J. Sci. Comput., 21(5): , Iterative metods for solving systems of algebraic equations (Copper Mountain, CO, 1998). [27] E. Lee and T.A. Manteuffel. FOSLL* metod for te eddy current problem wit tree-dimensional edge singularities. SIAM J. Numer. Anal., 45(2): , [28] T. A. Manteuffel, S. F. McCormick, J. Ruge, and J. G. Scmidt. First-order system. [29] P Monk. Finite element metods for Maxwell s equations. Numerical Matematics and Scientific Computation. Oxford University Press, New York, [30] R.H. Nocetto and J. Pyo. Optimal relaxation parameter for te Uzawa metod. Numer. Mat, 98: , [31] A. Quarteroni and A. Valli. Domain Decomposition Metods for Partial Differential Equations. Clarendon Press, Oxford, [32] F.J. Sayas. Infimum-supremum. Bol. Soc. Esp. Mat. Apl. S ema, (41):19 40, [33] L. R. Scott and M. Vogelius. Conforming finite element metods for incompressible and nearly incompressible continua. In Large-scale computations in fluid mecanics, Part 2 (La Jolla, Calif., 1983), volume 22 of Lectures in Appl. Mat., pages Amer. Mat. Soc., Providence, RI, [34] L. R. Scott and M. Vogelius. Norm estimates for a maximal rigt inverse of te divergence operator in spaces of piecewise polynomials. RAIRO Modél. Mat. Anal. Numér., 19(1): , [35] R. Verfürt. A combined conjugate gradient-multigrid algoritm for te numerical solution of te Stokes problem. IMA J. Numer. Anal., 4(4): , [36] J. Xu and L. Zikatanov. Some observations on Babuška and Brezzi teories. Numer. Mat., 94(1): , 2003.
21 SADDLE POINT LEAST SQUARES 21 University of Delaware, Department of Matematics, 501 Ewing Hall address: University of Delaware, Department of Matematics, 501 Ewing Hall address:
A SADDLE POINT LEAST SQUARES METHOD FOR SYSTEMS OF LINEAR PDES. Klajdi Qirko
A SADDLE POINT LEAST SQUARES METHOD FOR SYSTEMS OF LINEAR PDES by Klajdi Qirko A dissertation submitted to the Faculty of the University of Delaware in partial fulfillment of the requirements for the degree
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 informationSchur Complements on Hilbert Spaces and Saddle Point Systems
Schur Complements on Hilbert Spaces and Saddle Point Systems Constantin Bacuta Mathematical Sciences, University of Delaware, 5 Ewing Hall 976 Abstract For any continuous bilinear form defined on a pair
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationarxiv: v1 [math.na] 31 Oct 2014
arxiv:1411.0634v1 [math.na] 31 Oct 014 A NOTE ON STABILITY AND OPTIMAL APPROXIMATION ESTIMATES FOR SYMMETRIC SADDLE POINT SYSTEMS CONSTANTIN BACUTA Abstract. We establish sharp well-posedness and approximation
More informationarxiv: v1 [math.na] 20 Jul 2009
STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed
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 informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
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 informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationHYBRIDIZED GLOBALLY DIVERGENCE-FREE LDG METHODS. PART I: THE STOKES PROBLEM
MATHEMATICS OF COMPUTATION Volume 75, Number 254, Pages 533 563 S 0025-5718(05)01804-1 Article electronically publised on December 16, 2005 HYBRIDIZED GLOBALLY DIVERGENCE-FREE LDG METHODS. PART I: THE
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 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 informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
More informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationarxiv: v1 [math.na] 6 Dec 2010
MULTILEVEL PRECONDITIONERS FOR DISCONTINUOUS GALERKIN APPROXIMATIONS OF ELLIPTIC PROBLEMS WITH JUMP COEFFICIENTS BLANCA AYUSO DE DIOS, MICHAEL HOLST, YUNRONG ZHU, AND LUDMIL ZIKATANOV arxiv:1012.1287v1
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
More informationERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS
ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS YONGYONG CAI, AND JIE SHEN Abstract. We carry out in tis paper a rigorous error analysis
More informationA UNIFORM INF SUP CONDITION WITH APPLICATIONS TO PRECONDITIONING
A UNIFORM INF SUP CONDIION WIH APPLICAIONS O PRECONDIIONING KEN ANDRE MARDAL, JOACHIM SCHÖBERL, AND RAGNAR WINHER Abstract. A uniform inf sup condition related to a parameter dependent Stokes problem is
More informationCOMPUTATIONAL COMPARISON BETWEEN THE TAYLOR HOOD AND THE CONFORMING CROUZEIX RAVIART ELEMENT
Proceedings of ALGORITMY 5 pp. 369 379 COMPUTATIONAL COMPARISON BETWEEN E TAYLOR HOOD AND E CONFORMING OUZEIX RAVIART ELEMENT ROLF KRAHL AND EBERHARD BÄNSCH Abstract. Tis paper is concerned wit te computational
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More informationDownloaded 11/15/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 55, No. 6, pp. 2787 2810 c 2017 Society for Industrial and Applied Matematics EDGE ELEMENT METHOD FOR OPTIMAL CONTROL OF STATIONARY MAXWELL SYSTEM WITH GAUSS LAW IRWIN YOUSEPT
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationNumerical Analysis of the Double Porosity Consolidation Model
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real 21-25 septiembre 2009 (pp. 1 8) Numerical Analysis of te Double Porosity Consolidation Model N. Boal
More informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
More informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationUniform estimate of the constant in the strengthened CBS inequality for anisotropic non-conforming FEM systems
Uniform estimate of te constant in te strengtened CBS inequality for anisotropic non-conforming FEM systems R. Blaeta S. Margenov M. Neytceva Version of November 0, 00 Abstract Preconditioners based on
More informationFigure 1. Figure. Anticipating te use of two-level domain decomposition preconditioners, we construct a triangulation of in te following way. Let be d
Lower Bounds for Two-Level Additive Scwarz Preconditioners wit Small Overlap * Susanne C. Brenner Department of Matematics University of Sout Carolina Columbia, SC 908 Summary. Lower bounds for te condition
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
More informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationSTABILITY OF DISCRETE STOKES OPERATORS IN FRACTIONAL SOBOLEV SPACES
STABILITY OF DISCRETE STOKES OPERATORS IN FRACTIONAL SOBOLEV SPACES JEAN-LUC GUERMOND,, JOSEPH E PASCIAK Abstract Using a general approximation setting aving te generic properties of finite-elements, we
More informationarxiv: v1 [math.na] 7 Mar 2019
Local Fourier analysis for mixed finite-element metods for te Stokes equations Yunui He a,, Scott P. MacLaclan a a Department of Matematics and Statistics, Memorial University of Newfoundland, St. Jon
More informationADAPTIVE MULTILEVEL INEXACT SQP METHODS FOR PDE CONSTRAINED OPTIMIZATION
ADAPTIVE MULTILEVEL INEXACT SQP METHODS FOR PDE CONSTRAINED OPTIMIZATION J CARSTEN ZIEMS AND STEFAN ULBRICH Abstract We present a class of inexact adaptive multilevel trust-region SQP-metods for te efficient
More informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
More informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
More informationarxiv: v3 [math.na] 31 May 2016
Stability analysis of pressure correction scemes for te Navier-Stoes equations wit traction boundary conditions Sangyun Lee a, Abner J. Salgado b a Center for Subsurface Modeling, Institute for Computational
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
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 informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationA SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR
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 informationDual-Primal Isogeometric Tearing and Interconnecting Solvers for large-scale systems of multipatch continuous Galerkin IgA equations
www.oeaw.ac.at Dual-Primal Isogeometric Tearing and Interconnecting Solvers for large-scale systems of multipatc continuous Galerkin IgA equations C. Hofer, U. Langer RICAM-Report 2015-44 www.ricam.oeaw.ac.at
More informationChemnitz Scientific Computing Preprints
G. Of G. J. Rodin O. Steinbac M. Taus Coupling Metods for Interior Penalty Discontinuous Galerkin Finite Element Metods and Boundary Element Metods CSC/11-02 Cemnitz Scientific Computing Preprints Impressum:
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationStokes Equation. Chapter Mathematical Formulations
Capter 5 Stokes Equation 5.1 Matematical Formulations We introduce some notations: Let v = (v 1,v 2 ) T or v = (v 1,v 2,v 3 ) T. Wen n = 2 we define curlv = v 2 x v 1 y, ) curlη =. ( η y, η x In fact,
More informationBALANCING DOMAIN DECOMPOSITION FOR PROBLEMS WITH LARGE JUMPS IN COEFFICIENTS
MATHEMATICS OF COMPUTATION Volume 65, Number 216 October 1996, Pages 1387 1401 BALANCING DOMAIN DECOMPOSITION FOR PROBLEMS WITH LARGE JUMPS IN COEFFICIENTS JAN MANDEL AND MARIAN BREZINA Abstract. Te Balancing
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationLinearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization
Linearized Primal-Dual Metods for Linear Inverse Problems wit Total Variation Regularization and Finite Element Discretization WENYI TIAN XIAOMING YUAN September 2, 26 Abstract. Linear inverse problems
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationarxiv: v1 [math.na] 11 May 2018
Nitsce s metod for unilateral contact problems arxiv:1805.04283v1 [mat.na] 11 May 2018 Tom Gustafsson, Rolf Stenberg and Jua Videman May 14, 2018 Abstract We derive optimal a priori and a posteriori error
More informationAPPROXIMATION BY QUADRILATERAL FINITE ELEMENTS
MATHEMATICS OF COMPUTATION Volume 71, Number 239, Pages 909 922 S 0025-5718(02)01439-4 Article electronically publised on Marc 22, 2002 APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS DOUGLAS N. ARNOLD,
More informationSmoothed projections in finite element exterior calculus
Smooted projections in finite element exterior calculus Ragnar Winter CMA, University of Oslo Norway based on joint work wit: Douglas N. Arnold, Minnesota, Ricard S. Falk, Rutgers, and Snorre H. Cristiansen,
More informationMIXED AND STABILIZED FINITE ELEMENT METHODS FOR THE OBSTACLE PROBLEM
MIXED AND STABILIZED FINITE ELEMENT METHODS FOR THE OBSTACLE PROBLEM TOM GUSTAFSSON, ROLF STENBERG, AND JUHA VIDEMAN Abstract. We discretize te Lagrange multiplier formulation of te obstacle problem by
More informationFINITE ELEMENT EXTERIOR CALCULUS: FROM HODGE THEORY TO NUMERICAL STABILITY
BULLETIN (New Series) OF THE AMERICAN MATHEMATICAL SOCIETY Volume 47, Number 2, April 2010, Pages 281 354 S 0273-0979(10)01278-4 Article electronically publised on January 25, 2010 FINITE ELEMENT EXTERIOR
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationSOME IMPLEMENTATIONS OF PROJECTION METHODS FOR NAVIER STOKES EQUATIONS
SOME IMPLEMENTATIONS OF PROJECTION METHODS FOR NAVIER STOKES EQUATIONS JEAN-LUC GUERMOND Abstract. Tis paper is concerned wit te implementation of spatially discrete versions of Corin Temam s projection
More informationA New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions
Numerisce Matematik manuscript No. will be inserted by te editor A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions Wang Ming 1, Zong-ci Si 2, Jincao Xu1,3 1 LMAM, Scool of Matematical
More informationMixed Finite Element Methods for Incompressible Flow: Stationary Stokes Equations
Mixed Finite Element Metods for Incompressible Flow: Stationary Stoes Equations Ziqiang Cai, Carles Tong, 2 Panayot S. Vassilevsi, 2 Cunbo Wang Department of Matematics, Purdue University, West Lafayette,
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
More informationON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD. January 2000 Revised April 2000
ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD JAN MANDEL AND RADEK TEZAUR January 2000 Revised April 2000 Abstract In te Dual-Primal FETI metod, introduced by Farat et al [5], te domain is
More informationKey words. Sixth order problem, higher order partial differential equations, biharmonic problem, mixed finite elements, error estimates.
A MIXED FINITE ELEMENT METHOD FOR A SIXTH ORDER ELLIPTIC PROBLEM JÉRÔME DRONIOU, MUHAMMAD ILYAS, BISHNU P. LAMICHHANE, AND GLEN E. WHEELER Abstract. We consider a saddle point formulation for a sixt order
More informationOn convergence of the immersed boundary method for elliptic interface problems
On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many
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 informationFINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES
FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES MICHAEL HOLST AND CHRIS TIEE ABSTRACT. Over te last ten years, te Finite Element Exterior Calculus (FEEC) as
More informationA Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems
J Sci Comput (07 7:95 8 DOI 0.007/s095-06-096-4 A Weak Galerkin Metod wit an Over-Relaxed Stabilization for Low Regularity Elliptic Problems Lunji Song, Kaifang Liu San Zao Received: April 06 / Revised:
More informationEXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD
MATHEMATICS OF COMPUTATION Volume 68, Number 228, Pages 1379 1396 S 0025-5718(99)01093-5 Article electronically publised on February 24, 1999 EXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD
More informationAnalysis of time-dependent Navier-Stokes flow coupled with Darcy
Analysis of time-dependent Navier-Stokes flow coupled wit Darcy flow Ayçıl Çeşmelioğlu and Béatrice Rivière Abstract Tis paper formulates and analyzes a weak solution to te coupling of time-dependent Navier-Stokes
More informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationGrad-div stabilization for the evolutionary Oseen problem with inf-sup stable finite elements
Noname manuscript No. will be inserted by te editor Grad-div stabilization for te evolutionary Oseen problem wit inf-sup stable finite elements Javier de Frutos Bosco García-Arcilla Volker Jon Julia Novo
More informationAdvancements In Finite Element Methods For Newtonian And Non-Newtonian Flows
Clemson University TigerPrints All Dissertations Dissertations 8-3 Advancements In Finite Element Metods For Newtonian And Non-Newtonian Flows Keit Galvin Clemson University, kjgalvi@clemson.edu Follow
More informationH(div) conforming and DG methods for incompressible Euler s equations
H(div) conforming and DG metods for incompressible Euler s equations Jonny Guzmán Filánder A. Sequeira Ci-Wang Su Abstract H(div) conforming and discontinuous Galerkin (DG) metods are designed for incompressible
More informationCONSTRUCTIVELY WELL-POSED APPROXIMATION METHODS WITH UNITY INF SUP AND CONTINUITY CONSTANTS FOR PARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICS OF COMPUTATION Volume 82, Number 284, October 203, Pages 923 952 S 0025-578(203)02697-X Article electronically publised on April 23, 203 CONSTRUCTIVELY WELL-POSED APPROXIMATION METHODS WITH
More informationSome Applications of Fractional Step Runge-Kutta Methods
Some Applications of Fractional Step Runge-Kutta Metods JORGE, J.C., PORTERO, L. Dpto. de Matematica e Informatica Universidad Publica de Navarra Campus Arrosadia, s/n 3006 Pamplona Navarra SPAIN Abstract:
More informationA Study on Using Hierarchical Basis Error Estimates in Anisotropic Mesh Adaptation for the Finite Element Method
A Study on Using Hierarcical Basis Error Estimates in Anisotropic Mes Adaptation for te Finite Element Metod Lennard Kamenski arxiv:06.603v3 [mat.na] 20 Jan 202 Abstract: A common approac for generating
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 informationEffect of Numerical Integration on Meshless Methods
Effect of Numerical Integration on Mesless Metods Ivo Babuška Uday Banerjee Jon E. Osborn Qingui Zang Abstract In tis paper, we present te effect of numerical integration on mesless metods wit sape functions
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More informationarxiv: v1 [math.na] 20 Nov 2018
An HDG Metod for Tangential Boundary Control of Stokes Equations I: Hig Regularity Wei Gong Weiwei Hu Mariano Mateos Jon R. Singler Yangwen Zang arxiv:1811.08522v1 [mat.na] 20 Nov 2018 November 22, 2018
More informationSMAI-JCM SMAI Journal of Computational Mathematics
SMAI-JCM SMAI Journal of Computational Matematics Compatible Maxwell solvers wit particles II: conforming and non-conforming 2D scemes wit a strong Faraday law Martin Campos Pinto & Eric Sonnendrücker
More informationarxiv: v1 [math.na] 19 Mar 2018
A primal discontinuous Galerkin metod wit static condensation on very general meses arxiv:1803.06846v1 [mat.na] 19 Mar 018 Alexei Lozinski Laboratoire de Matématiques de Besançon, CNRS UMR 663, Univ. Bourgogne
More information