An Equal-order DG Method for the Incompressible Navier-Stokes Equations
|
|
- Jeffery Holland
- 5 years ago
- Views:
Transcription
1 An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp , 009 Abstract We introduce and analyze a discontinuous Galerkin method for the incompressible Navier-Stokes equations that is based on finite element spaces of the same polynomial order for the approximation of the velocity and the pressure. Stability of this equal-order approach is ensured by a pressure stabilization term. A simple elementby-element post-processing procedure is used to provide globally divergence-free velocity approximations. For small data, we prove the existence and uniqueness of discrete solutions and carry out an error analysis of the method. A series of numerical results are presented that validate our theoretical findings. 1 Introduction In this paper, we introduce and study, theoretically as well as computationally, a new discontinuous Galerkin (DG) method for the stationary School of Mathematics, University of Minnesota, Vincent Hall, Minneapolis, MN 55455, USA, cockburn@math.umn.edu Department of Mathematics, Texas A&M University, Mailstop 3368, College Station, TX 77843, USA, kanschat@tamu.edu. Supported in part by NSF through award no. DMS and by award no. US-C , made by ing Abdullah University of Science and Technology (AUST). Mathematics Department, University of British Columbia, Mathematics Road, Vancouver V6T 1Z, Canada, schoetzau@math.ubc.ca. Supported in part by the Natural Sciences and Engineering Research Council of Canada (NSERC). 1
2 incompressible Navier-Stokes equations, namely, ν u + (u u) + p = f in Ω, u = 0 in Ω, (1.1) u = 0 on Γ = Ω. Here, ν is the kinematic viscosity, u the velocity field, p the pressure, and f an external body force. For simplicity, we take Ω to be a polygonal domain in R. This paper is the last in a series of papers devoted to the study of DG methods for incompressible flows. Indeed, in [9], a local discontinuous Galerkin (LDG) method for the Stokes system was introduced. The method was then extended to the Oseen problem in [5]; see also the review [6]. The extension to the stationary incompressible Navier-Stokes was carried out later in [7], where an unexpected difficulty, introduced by the presence of the nonlinearity of the equations, was uncovered. It can be stated, roughly speaking, as follows: it is not possible to devise locally conservative DG methods which are also energy-stable, unless the approximate velocity is globally divergence-free. In [7], it was shown how to overcome this difficulty by using a sequence of suitably defined DG solutions of Oseen problems. Let us sketch it. A classic argument to find the solution of problem (1.1) consists in generating a sequence {(u k, p k )} k 0 of solutions of Oseen problems ν u k+1 + (w )u k+1 + p k+1 = f in Ω, u k+1 = 0 in Ω, u k+1 = 0 on Γ, where w = u k. This sequence converges to the solution (u, p) of (1.1) provided the data satisfies a small data condition of the form ν f 0 < ρ Ω, (1.) for a positive constant ρ Ω only depending on Ω; see, for example, [15]. The method proposed in [7] is a discrete version of this idea. It consists in generating a sequence {(u k h, pk h )} k 0 of approximations to the above solutions which then converge to the approximation sought. In this case, (u k+1 h, p k+1 h ) is given by a DG discretization of the above Oseen problem where the flow field w = P (u k h ) is now taken to be a post-processed velocity, obtained from u k h in such a way that it turns out to be globally divergencefree. The post-processing operator P can be computed in a simple, elementby-element fashion. DG methods in which the post-processing P can be
3 taken as the identity were briefly mentioned in [7] and explicitly further explored in [8, 11]; these methods provide exactly divergence-free velocities. In [7], the case of approximate velocities which are piecewise polynomials of degree k and approximate pressures that are piecewise polynomials of degree k 1 was considered. In this paper, we treat the more complicated case of equal-order spaces where the approximate pressures are also piecewise polynomials of degree k. Although the elements are now no longer matched properly with respect to their approximation properties, equalorder approaches are often observed to be more computationally efficient than their mixed-order counterparts; see for example the numerical tests in [9] for Stokes flow problems. For the non-linear Navier-Stokes equations, the main difficulty in the analysis of the equal-order case stems from the fact that, in order to compensate for using a bigger polynomial space for the approximate pressure, the numerical traces for the velocity and pressure cannot be taken as in [7], that is, as û p h = {u h } and p h = {p h }, but have to be modified as follows: û p h = {u h } + γ E [p h n] and p h = {p h }, where the stabilization function γ E is strictly positive. This function then introduces a pressure stabilization term in the variational formulation of the method. As a direct consequence, the above-mentioned post-processing projection P depends not only on u h but also on p h, since û p h depends on u h and p h. If the problem were linear, this would not cause any major difficulty, as was shown in [5, 9] and [16]. However, the presence of the nonlinearity renders the analysis of this modification quite delicate. The main contribution of this paper is to show that the optimal order of convergence can be obtained for the DG norm error of the approximation in the velocity. On the other hand, the theoretical and observed convergence rates for the L -norm errors in the pressure are suboptimal, as is expected for equal-order methods. The paper is organized as follows. In Section, we describe the equalorder DG method. In Section 3, we state and discuss our main results: the existence and uniqueness of DG approximations and a-priori error estimates. Sections 4, 5, 6 and 7 are devoted to the proofs of these results. In Section 8, we display numerical results showing that the orders of convergence predicted by our theoretical estimates are sharp. We end in Section 9 with some concluding remarks. 3
4 Equal-order DG Discretization.1 Notation We consider subdivisions T h of (global) mesh size h of the domain into triangles or rectangles. We assume that the meshes do not have hanging nodes and are shape-regular. We denote by E h the set of all edges of T h, by E I h the set of all interior edges, and by EB h the set of all boundary edges. As usual, we write h E for the length of edge E. Finally, we denote by E() the set of all elemental edges of. Using the space P k () of polynomials of degree up to k on, we define the discontinuous finite element space P k h = { v L (Ω) : v P k (), T h }. For a given polynomial order k 1, we introduce the discrete velocity and pressure spaces by V h = (Ph k ), Q h = { q Ph k : q dx = 0 }. (.1) In order to define derivative operators on these spaces, we introduce average and jump operators on edges. To do so, let + and be two adjacent elements of T h that share an interior edge E = + E I h. Let ϕ be any piecewise smooth function (vector- or scalar-valued), and let us denote by ϕ ± the traces of ϕ on E taken from within the interior of ±. On the edge E, we then define the mean value { } as {ϕ } = 1 (ϕ+ + ϕ ). Furthermore, for piecewise smooth vector-valued function, we define the jumps across E as [v n] = v + n + + v n, [v n] = v + n + + v n, [qn] = q + n + + q n. where n ± denotes the outward unit normal vector on the boundary ±. For a boundary edge E = Γ E B h, we analogously define {ϕ } = ϕ+, [v n] = v + n, [v n] = v + n and [qn] = q + n, where we now denote by ϕ +, v +, q + the traces of ϕ, v, q on, respectively, and by n the unit outward normal on Γ. Ω 4
5 . The Oseen problem We begin the description of the equal-order method by defining the DG discretization of the Oseen problem ν u + (w )u + p = f in Ω, u = 0 in Ω, u = 0 on Γ. Here, w is a generic divergence-free flow field in the space J(Ω) = { w H div (Ω) : w 0 in Ω }, where H div (Ω) is the space of vector functions in L (Ω) whose divergence is in L (Ω). It reads as follows: for w J(Ω), find u V h and p Q h such that a h (u, v) + o h (w; u, v) + b h (v, p) = f v dx v V h, Ω (.) b h (u, q) + c h (p, q) = 0 q Q h. Here, the forms a h, o h and b h correspond to ν times the Laplacian, the advection term and the divergence operator, respectively. The form c h is associated with the pressure jumps and plays a stabilization role. These forms are specified below...1 The discrete Laplacian Since DG formulations for the Laplacian have been intensely researched during the last years, see for example [], we refrain from listing special discretizations and will simply make assumptions on the method used. To state them we need to introduce the standard DG norm on the space V h : u 1,h = T h u L () + E E h E κ 0 h E [u n] ds, where κ 0 is a positive parameter depending on the polynomial degree k, but not on the mesh size h. Assumption.1 With u the solution to the vector-valued equation ν u = f in Ω, u = 0 on Γ, we make the following assumptions on the discretizing DG bilinear form a h : 5
6 Boundedness: there is a positive constant c a independent of ν and the mesh size h such that for any u and v in V h a h (u, v) ν c a u 1,h v 1,h. (.3a) Stability: there is a positive constant α independent of ν and the mesh size h such that for any u V h a h (u, u) ν α u 1,h. (.3b) Consistency: for the solution u above and any v V h there holds a h (u, v) = f v dx. (.3c) Approximation property: if the solution u above belongs to H s+1 (Ω), for some exponent s 1, and P Vh u is its L -projection onto V h, then there is a positive constant C independent of ν and the mesh size h such that for any v V h Ω a h (u P Vh u, v) Ch min{s,k} ν 1 u s+1 ν 1 v 1,h (.3d) We remark, that the interior penalty methods and the LDG methods are among those DG schemes for which these assumptions hold. Note also that, while the energy norm error estimates from [5, 7, 9] were proven for the LDG discretization of the Laplacian, their proofs rely on these assumptions only; see also [, 10]. For the optimal L -norm or negativeorder norm error estimates derived in [5, 9], also adjoint consistency of the bilinear form a h is needed: for the solution u above and any v V h there holds a h (v, u) = f v dx. In the following, we will not make use of any duality argument. Therefore, adjoint consistency is not necessary for the analysis presented in this article... The advection form For the convective term, we take the standard upwind flux introduced in [13], namely, o h (w; u, v) = ( ) u w : v dx (w n)u v ds, T h Ω 6
7 where u is the upwind value of u on with respect to w, namely, { u limɛ 0 u(x ɛ w(x)), x \ Γ, (x) = 0, x Γ, with Γ denoting the inflow part of Γ: Γ = { x Γ : w(x) n(x) < 0 }. As observed in [13], the following result holds. Proposition. If w J(Ω), we have O h (w; u, u) 0 u V h. Let us further recall the Lipschitz continuity of the form o h with respect to the convective term; cf. [7, Proposition 4.]. To that end, we define We now have the following result. V(h) = H 1 0 (Ω) + V h. Proposition.3 Let w 1, w J(Ω), u V(h) and v V h. Then we have o h (w 1 ; u, v) o h (w ; u, v) c o w 1 w 1,h u 1,h v 1,h, with a constant c o > 0 that is independent of ν and the mesh size h...3 The discrete divergence We take b h (u, q) = u q dx {u } [qn] ds T h c h (p, q) = ν 1 γ 0 h E E E I E h E E E h [pn] [qn] ds, where γ 0 is a positive parameter. These forms correspond to a discretization of the incompressibility condition on the exact velocity similar to that used in [7], namely, u h q dx + û p h n q ds = 0. (.4) 7
8 Here, the numerical fluxes associated with the incompressibility constraint, û p h and p h, are defined as follows. If E E I h is an interior edge, we take û p h = {u h } + γ E [p h n], p h = {p h }, (.5a) for a stabilization function which, in accordance with the discussion in [5, Appendix], we take as γ E = ν 1 γ 0 h E. (.5b) Note that if we set γ 0 = 0, we obtain the numerical traces used in [7] for the mixed-order DG method. On boundary edges E E B h, we set û p h = 0, p h = p h. (.5c) After integration by parts, it can be seen that the form b h is bounded. Proposition.4 For functions v V h and q Q h, the form b h is bounded by b h (v, q) c b v 1,h q 0, with a constant c b > 0 that is independent of ν and the mesh size h. Moreover, the form b h satisfies the following stability condition on V h Q h. This condition has been shown in [16] and in a slightly different form in [9]. To state it, we need to introduce the semi norm h on Q h given by p h = γ 0 h E [p] L (E). E E I h Proposition.5 There exist constants β 1 > 0 and β > 0 independent of ν and the mesh size h such that b h (v, q) sup β 1 q 0 β q 0 v V h v h q Q h. 1,h.3 Post-processing Next, we introduce the local post-processing. Let the function (u, p) be piecewise smooth, and let û p be the numerical flux associated with (u, p) as in (.5). We define the post-processed velocity w = P (u, p) as follows. The restriction of w to is the function in BDM k () defined via the following moments on and its edges E: w n ϕ ds = û p n ϕ ds ϕ P k (E), E E(), (.6a) E E w v dx = u v dx v M k (). (.6b) 8
9 Here, when is a triangle, we have BDM k () = ( P k () ), M k () = { ϕ + (b ψ) : (ϕ, ψ) P k 1 () P k () }, with b being the elemental polynomial bubble in P 3 () vanishing on. When is a rectangle, we have BDM k () = ( P k () ) span { (x k+1 y), (x y k+1 ) }, M k () = ( P k () ). Note that these moments are closely related to those of the standard Brezzi- Douglas-Marini (BDM) projection; see [4, Section III.3]. In fact, it can be easily seen that, if (u, p) are in H (Ω) H 1 (Ω), then P (u, p) is simply the BDM projection of u (cf. [4, Section III.3]). Evidently, the post-processed velocity w is well-defined and can be computed in an element-by-element fashion. Moreover, it belongs to the space Ṽ h = { v H div (Ω) : v BDM k (), T h }. The main property of w = P (u h, p h ) is that it is exactly divergencefree provided that (u h, p h ) V h Q h satisfies the weak incompressibility condition in (.). Proposition.6 Let (u h, p h ) V h Q h satisfy the second equation in (.) and let w = P (u h, p h ) Ṽh be the post-processed velocity associated with (u h, p h ). Then w is exactly divergence-free and belongs to J(Ω). Proof For q Q h, we obtain w q dx = ( Ω T h = T h ( ) w q dx + w n q ds ) u h q dx + û p h n q ds, by the definition of w in (.6). Finally, rewriting the weak incompressibility condition in (.) in the flux formulation (.4), we obtain Ω w q dx = 0. We claim that w Q h. As a consequence, we conclude that w 0 in Ω. 9
10 To prove the claim, we note that w P k () for all T h and w dx = w n ds = n ds = 0, Ω Γ û p h Γ which in turn is a consequence of the definition of w in (.6) and the fact that we have û p h = 0 on Γ. This completes the proof. The following stability result will be crucial in the analysis of the method. It is analogous to the stability result of Proposition 4.6 in [7]. Proposition.7 For piecewise polynomial functions u and p, the postprocessing operator P admits the stability estimate P ( (u, p) 1,h c s u 1,h + α 1 ν p ) 1 h, with a constant c s > 0 that is independent of ν and the mesh size h. In the estimate of Proposition.7, the coercivity constant α of the bilinear form a h is introduced for convenience only; this allows us to avoid certain maximum terms in the proof of existence and uniqueness of discrete solutions. The proof of this result will be presented in Section 4..4 The Navier-Stokes equations We are finally ready to present the equal-order DG method for the incompressible Navier-Stokes equation (1.1). It is nothing else than the DG discretization of the Oseen problem (.) where the convective velocity field w h is the post-processed velocity w h = P (u h, p h ): find u h V h and p h Q h such that a h (u h, v) + o h (w h ; u h, v) + b h (v, p h ) = f v dx v V h, Ω (.9) b h (u h, q) + c h (p h, q) = 0 q Q h. As in [7], we compute the solution of this nonlinear system by the Picard iteration mentioned in Section 1. We rewrite it as follows. For w Ṽh J(Ω), let ( T u (w), T p (w) ) V h Q h be the solution of a h (T u (w), v) + o h (w; T u (w), v) + b h (v, T p (w)) = f v dx v V h, b h (T u (w), q) + c h (T p (w), q) = 0 q Q h. Furthermore, define T as the post-processing applied to the Oseen solution, namely Ω T : Ṽh J(Ω) Ṽh J(Ω), T (w) = P (Tu (w), T p (w)). 10
11 Obviously, the field w h in (.9) must be a fixed point of T. To find it, we set w 0 h = 0, and generate the sequence {wk h } k 0 by using the Picard iteration 3 Main Results w k+1 h = T (w k h). (.10) 3.1 Existence and uniqueness of approximate solutions Our first main result is concerned with the convergence of the Picard iteration and with the existence and uniqueness of the approximate solutions provided by the equal-order DG method. To state it, we need to introduce the following discrete small data condition: ν f 0 < ρ 1 for ρ 1 = α c o c p c s, (3.1) where c p > 0 is the constant in the Poincaré inequality v 1,h c 1 p v 0 v V(h). (3.) For a proof of this inequality, see, for example, [1, Lemma.1] or [3]. Theorem 3.1 Consider the equal-order DG method given by (.9) for which the assumptions (.3a) and (.3b) on the bilinear form a h hold. Under the small data condition (3.1), the Picard iteration in (.10) converges to a unique fixed point w h Ṽh J(Ω), and the method has a unique solution w h Ṽh J(Ω) and (u h, p h ) V h Q h. Moreover, the solution satisfies the stability bound ( uh 1,h + α 1 ν p h h) 1 c p αν f 0 and w h 1,h c pc s αν f 0. The proof of Theorem 3.1 is presented in Section Inf-sup stability Next, we present the inf-sup condition that we use to derive our error estimates. To that end, we introduce the global form A h (w)(u, p; v, q) = a h (u, v) + o h (w; u, v) + b h (v, p) b h (u, q) + c h (p, q), and the norm (u, p) = ( ν u 1,h + ν 1 p 0 + ν 1 p ) 1 h. 11
12 As a consequence of the generalized inf-sup condition for b h in Proposition.5 and the introduction of the pressure stabilization form c h, the following stability result holds for the form A h. Proposition 3. Assume the small data condition (3.1), and let w J(Ω) be bounded by w 1,h c sc p αν f 0. Then there is a constant c A > 0 independent of ν and the mesh size h such that inf (u,p) V h Q h A h (w)(u, p; v, q) sup (v,q) V h Q h (u, p) (v, q) c 1 A. The proof of Proposition 3. is given in Section Error estimates Our second main result are a-priori error estimates for the equal-order DG method. To derive them, we need a small data condition which is slightly different than that in (3.1), namely: ν f 0 < ρ for ρ = ( c A c o c p c s max{1, α 1 } ) 1. (3.3) In the following, we assume for simplicity that the strictest of the small data conditions in (1.), (3.1) and (3.3) holds, namely: ν f 0 < min{ρ Ω, ρ 1, ρ }. (3.4) This also ensures the existence and uniqueness of Navier-Stokes solutions and their DG approximations. We then have the following error estimates. Theorem 3.3 Assume that the small data condition (3.4) holds. Let the solution (u, p) of the incompressible Navier-Stokes equations (1.1) belong to (u, p) H s+1 (Ω) H s (Ω), s 1. Let (u h, p h ) V h Q h be the solution of the equal-order DG method given by (.9) for which the assumptions (.3) on the bilinear form a h hold, and let w h = P (u h, p h ) be the post-processed exactly divergence-free velocity in Ṽh J(Ω). Then we have the error estimates ν 1 u wh 1,h + (u u h, p p h ) Ch min{s,k} ( ν 1 u s+1 + ν 1 p s ), with a constant C > 0 that is independent of ν and the mesh size h. 1
13 The proof of Theorem 3.3 is given in Section 7. Remark 3.4 For smooth solutions, the result in Theorem 3.3 predicts convergence rates of the order O(h k ) for the error measured in the norm. This rate is optimal with respect to the DG norm error in the velocity. With respect to the L -norm error in the pressure, this result is suboptimal by one order since polynomials of degree k are used to approximate the pressure variable. The numerical results in Section 8 indicate that these suboptimal rates in the pressure are in fact observed in practice. For the DG norm error in the post-processed velocity, rates of the order O(h k ) are predicted, which is verified in our numerical experiments as well. In this sense, the error estimate of Theorem 3.3 is sharp. 4 Proof of Proposition.7 For each element T h, we introduce the space { V() = u BDM k () : u v dx = 0 and define u = E E() E u n ds. } v M k (), Obviously, is a norm on V(). This follows from the fact that all the interior BDM moments are zero for functions in V(); see [4]. Lemma 4.1 There is a constant C > 0 independent of the mesh size h such that h 1 u L () + u H 1 () 1 Ch u for any u V(). Proof On the reference element, the equivalence of norms on finite dimensional spaces shows that there is C > 0 such that ũ L ( e ) + ũ H 1 ( e ) C ũ e for all ũ V( ). Let us then consider an arbitrary triangle and u V(). We define ũ( x) = 1 det(df ) DF u(f ( x)), 13
14 where F : is the affine elemental mapping and DF its Jacobian. This mapping has the property that ũ V( ); see [4, Section III.1]. Furthermore, scaling arguments show that u 0, C ũ 0, e, u 1, Ch 1 ũ 1, e, with a constant C > 0 only depending on the shape-regularity of the mesh; cf. [4, Section III.1]. Next E E() be an edge of. The corresponding edge on the reference element is denoted by Ẽ. Using the properties of the mapping u ũ and the density argument from [15, Equation (3.4.)], we conclude that ũ n e L ( E) e = sup ee ũ n e ϕ d s eϕ P k ( E) e ϕ L ( E) e Ch 1 sup ϕ P k (E) E u n ϕ ds ϕ L (E) Ch 1 u n L (E). Here, we have also used that, if ϕ( x) = ϕ(f ( x)), then ϕ L (E) Ch 1 ϕ L ( E) e. The estimates above allow us the conclude that h 1 u L () + u H 1 () Ch 1 This completes the proof. Ch 1 ũ e Ch 1 u. ( ũ L ( e ) + ũ H 1 ( e ) For λ L (E h ), we now let U λ be defined elementwise as the unique polynomial in BDM k () satisfying: U λ n ϕ ds = λ n ϕ ds ϕ P k (E), E E(), E E U λ v dx = 0 v M k (). Clearly, we have that U λ Ṽh. Proposition 4. There is a constant C > 0 independent of the mesh size h such that U λ 1,h C E λ L (E). E E h h 1 ) 14
15 Proof For any element, we have U λ V(). Hence, Lemma 4.1 and the shape-regularity of the mesh yield U λ H 1 () h 1 U λ. T h T h Let now E E() be an edge of. The inverse estimate U λ L (E) Ch 1 U λ L (), the shape-regularity of the mesh and the result in Lemma 4.1 for the L - norm show that h 1 E U λ L (E) Ch 1 U λ, It then follows readily that U λ 1,h U λ. T h h 1 Testing the defining moments of U λ with ϕ = U λ n and applying the Cauchy-Schwarz inequality, we obtain that, on each edge E, U λ n ds λ n ds. E Taking into account the shape-regularity of the mesh, the assertion follows from the above estimates. We are now ready to complete the proof of Proposition.7. Let (u, p) be a piecewise polynomial function and let w = P (u, p) be the post-processed field obtained from (u, p). Taking into account the form of the flux (.5) associated with the incompressibility, we have where and E E E U 1 n ϕ ds = U 1 v dx = U n ϕ ds = E E E w = U 1 + U, {u } n ϕ ds ϕ P k (E), E E(), u v dx v M k (), γ E [p] n ϕ ds ϕ P k (E), E E() E I h, U n ϕ ds = 0 ϕ P k (E), E E() E B h, U v dx = 0 v M k (), 15
16 The function U 1 is the post-processed velocity used in the mixed-order approach in [7]. Hence, from [7, Proposition 4.6], there is a constant C > 0 independent of the mesh size h such that U 1 1,h C u 1,h. Moreover, Proposition 4. implies that there is constant C > 0 independent of the mesh size h such that U 1,h C h 1 E γ E [p] L (E). E E I h Since γ E = ν 1 γ 0 h E, see (.5b), we have w 1,h C ( u 1,h + ν p h) cs ( u 1,h + α 1 ν p h). This shows the result of Proposition.7. 5 Proof of the Theorem 3.1 To prove the existence and uniqueness of the approximate solutions of the equal-order DG method, we proceed as in [7] and mimic the proof of existence and uniqueness of the exact solution to problem (1.1); see also [15, Theorem ]. Thus, we show that, under the small data condition (3.1), the Picard iteration (.10) is a contraction on a suitably defined bounded set. We begin by obtaining key estimates on the solution of the Oseen problem (.). Lemma 5.1 If w J(Ω), then the Oseen problem (.) admits a unique solution (u h, p h ) V h Q h with the stability estimate ( uh 1,h + α 1 ν p h h) 1 c p αν f 0. Proof The existence and uniqueness of the solution (u h, p h ) V h Q h of the Oseen problem (.) can be obtained as in [5]; see also [16] for the Stokes equations. To obtain the estimates, we proceed as follows. Taking v = u h in the first equation of the formulation (.), q = p h in the second, and adding the resulting equations, we get a h (u h, u h ) + o h (w; u h, u h )+c h (p h, p h ) = f u h dx. Ω 16
17 Using the ellipticity assumption (.3b) on the form a h and the stability of the form o h (w; u h, u h ) guaranteed by Proposition., we get that να u h 1,h + ν 1 p h h Ω f u h dx, and, after applying the Cauchy-Schwarz inequality, that να u h 1,h + ν 1 p h h f 0 u h 0 c p f 0 u h 1,h, by Poincaré s inequality (3.). The estimate now follows. We are now ready to state and prove the main result of this subsection. Lemma 5. Under the small data condition (3.1), the mapping T defining the Picard iteration (.10) is a contraction on the set { h = w Ṽh J(Ω) : w 1,h c } sc p αν f 0. Proof First, we note that by the stability estimate of Proposition.7, we have that, for any w h, T (w) 1,h = P (T u (w), T p (w)) 1,h c s ( Tu (w) 1,h + α 1 ν T p (w) h c s c p αν f 0, by Lemma 5.1. Hence, T maps h into itself. Now, let w 1 and w be in h, and set (u 1, p 1 ) = (T u (w 1 ), T p (w 1 )) and (u, p ) = (T u (w ), T p (w )). It can be readily seen that the difference (u, p) = (u 1 u, p 1 p ) solves a h (u, v) + o h (w 1 ; u, v) + b h (v, p) b h (u, q) + c h (p, q) = L h (v), for any (v, q) V h Q h. Here, L h (v) is the functional L h (v) = o h (w 1 ; u, v) o h (w ; u, v). Thus, if we take v = u and q = p, we obtain a h (u, u) + o h (w 1 ; u, u) + c h (p, p) = L h (u), and hence, proceeding as in the proof of the previous result, that να u 1,h + ν 1 p h L h(u). ) 1 17
18 Let us now estimate the functional L h (u). In view of the Lipschitz continuity of the Oseen form o h in Proposition.3, we have L h (u) c o w 1 w 1,h u 1,h u 1,h, and by the estimate of solutions of the discrete Oseen problem (.), Lemma 5.1, applied to u, we obtain This implies that and hence that L h (u) c oc p αν f 0 w 1 w 1,h u 1,h. να u 1,h + ν 1 p h c oc p αν f 0 w 1 w 1,h u 1,h, ( ) u 1,h + α 1 ν p 1 h c oc p α ν f 0 w 1 w 1,h. Finally, by the stability property of the post-processing P of Proposition.7, we obtain T (w 1 ) T (w ) 1,h c s ( u 1,h + α 1 ν p h c oc p c s α ν f 0 w 1 w 1,h, By the small data condition (3.1), we conclude that T is a contraction on h. This completes the proof. Banach s fixed point theorem now implies that the Picard iteration in (.10) converges to a unique fixed point w h in h : w h = T (w h ). Thus, the unique solution (u h, p h ) V h Q h of the Oseen problem (.) associated with the convective field w h is the unique solution of (.9). This completes the proof of Theorem Proof of Proposition 3. To prove Proposition 3., fix (u, p) V h Q h. By the coercivity of a h in (.3b), the non-negativity of o h in Proposition., and the definition of h, we see that A h (w)(u, p; u, p) = a h (u, u) + o h (w; u, u) + c h (p, p) να u 1,h + ν 1 p h. (6.1) ) 1 18
19 Further, due to the generalized inf-sup condition for b h in Proposition.5, there is a velocity v V h such that b h ( v, p) β 1 p 0 β p 0 p h, v 1,h p 0. From the boundedness of a h and b h, (.3a) and Proposition.4, the continuity property of o h in Proposition.3, and the bound for v, we obtain A h (w)(u, p; v, 0) = a h (u, v) + o h (w; u, v) + b h ( v, p) β 1 p 0 β p 0 p h a h (u, v) o h (w; u, v) β 1 p 0 β p 0 p h (νc a + c o w 1,h ) u 1,h p 0. The small data condition (3.1) and the bound on w in Proposition 3. imply that c o w 1,h c oc s c p αν f 0 αν. Therefore, A h (w)(u, p; v, 0) β 1 p 0 β p 0 p h c 1 ν u 1,h p 0, (6.) with c 1 = c a +α. Next, we set (v, q) = (u, p)+δν 1 ( v, 0) with δ > 0 still to be chosen. From the bounds in (6.1), (6.), and the geometric-arithmetic inequality, we readily conclude that A h (w)(u, p; v, q) =A h (w)(u, p; u, p) + δν 1 A h (w)(u, p; v, 0) ( α δ ε ) ( 1c 1 ν u 1,h + 1 δ β ) ε ν 1 p h ( + δ β 1 c 1 β ) ν 1 p ε 1 ε 0, with parameters ε 1 > 0, ε > 0 still at our disposal. Choosing ε 1 and ε sufficiently large shows that A h (w)(u, p; v, q) (α δc ) ν u 1,h + δβ 1 3 ν 1 p 0 + (1 δc 3 ) ν 1 p h, with constants c > 0 and c 3 > 0 that are independent of ν and the mesh size h. Next, by choosing δ sufficiently small, we deduce the existence of a constant c 4 > 0 independent of ν and the mesh size h such that A h (w)(u, p; v, q) c 4 (u, p). (6.3) On the other hand, due the bound for v, we have (v, q) (u, p) + δν 1 ( v, 0) (u, p) + δν 1 v 1,h (u, p) + δν 1 p 0 c 5 (u, p). The bounds in (6.3) and (6.4) imply the result of Proposition 3.. (6.4) 19
20 7 Proof of Theorem Preliminaries To prove Theorem 3.3, let (u, p) denote the solution of the Navier-Stokes equations. We recall from [7, Equation (1.5)] that the velocity field u satisfies the stability bound u 1,h c p ν f 0. (7.1) Further, let (u h, p h ) V h Q h be the DG approximation from (.9) and w h = P (u h, p h ) Ṽh J(Ω) the post-processed velocity field. We note that, under the regularity assumptions of Theorem 3.3, we have the error equation A h (u)(u, p; v, q) A h (w h )(u h, p h ; v, q) = 0 (7.) for all (v, q) V h Q h. Here, we have used the consistency of all the DG forms involved; cf. assumption (.3c) on the bilinear form a h. 7. Proof of the error estimates Let us now bound the errors e u = u u h, e p = p p h and e w = u w h. To that end, we proceed in several steps. Step 1 (decomposition of the errors): As usual, we decompose the errors into two parts: e u = η u + ξ u = (u Π Vh u) + (Π Vh u u h ), e p = η p + ξ p = (p Π Qh p) + (Π Qh p p h ), where Π Vh and Π Qh denote the L -projections onto V h and Q h, respectively. Using the triangle inequality and the approximation properties of the L -projections, we obtain that (e u, e p ) (η u, η p ) + (ξ u, ξ p ) C E s (u, p) + (ξ u, ξ p ), (7.3) where we set, for notational convenience, E s (u, p) = h min{s,k} ( ν 1 u s+1 + ν 1 p s ). Hence, to bound the error in u and p, we need to further estimate the term (ξ u, ξ p ) on the right-hand side of (7.3). Step (a preliminary bound for (ξ u, ξ p )): We claim that (ξ u, ξ p ) ( c s max{1, α 1 } ) 1 ν 1 ew 1,h + C E s (u, p), (7.4) 0
21 with a constant C > 0 that is independent of ν and the mesh size h. Proof To prove this bound, we start from the inf-sup condition in Proposition 3.: there is a test function (v, q) V h Q h such that (v, q) = 1 and c 1 A (ξ u, ξ p ) A h (w h )(ξ u, ξ p ; v, q). Note that, by Lemma 5. and the small data condition (3.4), Proposition 3. is applicable. We then rewrite the right-hand side of the above inequality by using the error equation (7.) and obtain that A h (w h )(ξ u, ξ p ; v, q) =o h (w h ; ξ u, v) + o h (w h ; u h, v) o h (u; u, v) a h (η u, v) b h (v, η p ) + b h (η u, q) c h (η p, q). We continue to express the difference between two of the above trilinear forms as follows: o h (w h ; u h, v) o h (u; u, v) = o h (w h ; ξ u, v) o h (w h ; η u, v) Combining these identities yields where + o h (w h ; u, v) o h (u; u, v) c 1 A (ξ u, ξ p ) T 1 + T + T 3, (7.5) T 1 = a h (η u, v) + b h (v, η p ) + b h (η u, q) + c h (η p, q), T = o h (w h ; η u, v), T 3 = o h (w h ; u, v) o h (u; u, v). Let us now bound these terms. For T 1, we use the approximation property (.3d) for a h, the analogous approximation properties for b h and c h derived in [9, Corollary 3.8] and the fact that (v, q) = 1. We obtain that T 1 C E s (u, p) (v, q) C E s (u, p), with a constant C > 0 independent of ν and the mesh size h. To estimate T, we use the continuity bound of o h in Proposition.3, the stability result for P in Proposition.7 and the a-priori bound for (u h, p h ) in Theorem 3.1: ( T c o c s uh 1,h + α 1 ν p h ) 1 h η u 1,h v 1,h c oc s c p αν f 0 ν 1 ηu 1,h ν 1 v 1,h. 1
22 Hence, the small data condition (3.4), the fact that (v, q) = 1 and the approximation properties of the L -projection Π Vh yield T C E s (u, p), with a constant C > 0 independent of ν and the mesh size h. To estimate T 3, we proceed similarly and apply the continuity of o h in Proposition.3, the stability bound for u in (7.1) and the fact that (v, q) = 1. This results in T 3 c o e w 1,h u 1,h v 1,h c o c p ν f 0 ν 1 ew 1,h. From (7.5) and the bounds for the terms T 1, T and T 3, we now deduce that (ξ u, ξ p ) c A c o c p ν f 0 ν 1 ew 1,h + C E s (u, p), for a constant C > 0 independent of ν and the mesh size h. Using the small data condition (3.4), we see that c A c o c p ν f 0 ( c s max{1, α 1 } ) 1, which concludes the proof of (7.4). Step 3 (a preliminary bound for e w ): Next, we show that there holds ν 1 ew 1,h c s max{1, α 1 } (ξu, ξ p ) + C E s (u, p), (7.6) with a constant C > 0 independent of ν and the mesh size h. Proof Take an arbitrary function ṽ Ṽh. By the triangle inequality, e w 1,h u ṽ 1,h + ṽ w h 1,h. The second term on the right-hand side of this inequality can be estimated by employing the property that P (ṽ, 0) = ṽ, the stability bound for P in Proposition.7 and the fact that p h = 0 for p H 1 (Ω). This yields ṽ w h 1,h = P (ṽ u h, p h ) 1,h c s ( ṽ uh 1,h + α 1 ν p h h ) 1 c s ν ( 1 max{1, α 1 } ν ṽ uh 1,h + ν 1 p p h h = c s ν 1 max{1, α 1 } (ṽ uh, e p ). ) 1
23 Since (ṽ u h, e p ) = (ṽ u, 0) + (e u, e p ), a repeated application of the triangle inequality results in ṽ w h 1,h c s ν 1 max{1, α 1 } ( (u ṽ, 0) + (eu, e p ) ) c s ν 1 max{1, α 1 } (eu, e p ) + C u ṽ 1,h c s ν 1 max{1, α 1 } (ξu, ξ p ) + C (η u, η p ) + C u ṽ 1,h, with a constant C > 0 independent of ν and the mesh size h. Therefore, we conclude that ν 1 ew 1,h c s max{1, α 1 } (ξu, ξ p ) + C (η u, η p ) + C inf ν 1 u ṽ 1,h. ev V e h The bound in (7.6) now follows from the approximation properties of the L -projections and the ones of the space Ṽh. Step 4 (conclusion of the result in Theorem 3.3): We are now ready to prove our error estimates. Indeed, by combining the results (7.4) and (7.6) of Steps and 3, respectively, we obtain (ξ u, ξ p ) ( c s max{1, α 1 } ) 1 ν 1 ew 1,h + C E s (u, p) 1 (ξ u, ξ p ) + C E s (u, p). Absorbing the first term on the right-hand side of this inequality in the left-hand side shows that (ξ u, ξ p ) C E s (u, p), (7.7) with a constant C > 0 independent of ν and the mesh size h. The error estimate for (e u, e p ) now follows immediately from (7.3) and (7.7). To prove the estimate for the error e w = u w h in the post-processed velocity, we now only need to combine the bounds in (7.6) and (7.7): e w 1,h C (ξ u, ξ p ) + C E s (u, p) C E s (u, p), with a constant C > 0 independent of ν and the mesh size h. This completes the proof of Theorem Numerical Experiments In this section, we carry out numerical experiments that confirm our theoretical findings. The test problem we use was introduced by ovasznay [1]. 3
24 On the square Ω = [ 1, 3 ] [0, ], this solution is where λ = u 1 (x, y) = 1 e λx cos(πy), u (x, y) = λ π eλx sin(πy), p(x, y) = 1 eλx + p, 8π 3 ν 1 + ν + 16π, p = e λx dx. 1 We impose this function as Dirichlet boundary condition on the boundary of Ω and choose the viscosity ν = 1/10. In our implementation, we use the spaces V h and Q h in (.1) for different values of k on sequences of uniform meshes made of squares of size h = l for some integer l. As a DG discretization for the Laplacian, we choose the symmetric interior penalty form a h with a penalty parameter νk(k + 1)/h. It satisfies the properties listed in Assumption.1 with κ 0 = k(k + 1). The pressure stabilization parameter γ E in (.5b) is chosen as γ E = ν 1 h/10 independent of k. We use standard DG terms to incorporate the inhomogeneous Dirichlet boundary conditions. In Table 1, we report the history of convergence of the method. There, we denote by e u and e w the errors e u = u u h and e w = u w h, respectively. We see the expected convergence of order O(h k ) for the energy norm of the error, as well as for the standard DG norm 1,h of the errors in u h and w h, in agreement with Theorem 3.3. Additionally, we see that the post-processing improves the error. The norm of the pressure converges with the order O(h k ) for k = 1 and k = 3; a slightly higher order is observed for k =. In Table, we report the maximum norm of the (piecewise) divergence of the approximate velocities. Note that u h does not converge faster than e u 1,h, and that the post-processed divergence is of the order of the computational accuracy. 9 Concluding Remarks Let us end by pointing out that what has been done here for the two dimensional case can also be carried out for the three-dimensional one. Tetrahedra and cubic prisms could be easily handled. 4
25 h (e u, e p ) order e u 1,h e w 1,h e p 0 order k = 1 (linear) e e e-0 5.7e e e-1 6.7e-1.87e e e e e e e e e k = (quadratic) e e e- 1.45e e e-.3e-.07e e e e e e e e e-6.6 k = 3 (cubic) e e e-3 1.9e e e e-4.3e e e e-5.89e e e e e Table 1: History of convergence of the errors for the velocities u h, w h and the pressure p h for the ovasznay flow (ν = 1/10). Polynomials of degree k and meshes made of squares of size h are used. sup u h sup w h h k = 1 k = k = 3 k = 1 k = k = e e e e-07.80e e e-1.16e- 1.0e-3 4.0e e e e e-3 1.7e-4.36e e e e- 1.16e e-5.43e e e-10 Table : History of convergence of the divergence of each of the two approximate velocities inside the elements for the ovasznay flow (ν = 1/10). Polynomials of degree k and meshes made of squares of size h are used. 5
26 References [1] D. N. Arnold, An interior penalty finite element method with discontinuous elements, SIAM J. Numer. Anal. 19 (198), [] D. N. Arnold, F. Brezzi, B. Cockburn, and L. D. Marini, Unified analysis of discontinuous Galerkin methods for elliptic problems, SIAM J. Numer. Anal. 39 (00), [3] S. Brenner, Poincaré-Friedrichs inequalities for piecewise H 1 functions, SIAM J. Numer. Anal. 41 (003), [4] F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods, Springer Verlag, [5] B. Cockburn, G. anschat, and D. Schötzau, Local discontinuous Galerkin methods for the Oseen equations, Math. Comp. 73 (004), [6], The local discontinuous Galerkin methods for linear incompressible flow: A review, Computer and Fluids 34 (005), [7], A locally conservative LDG method for the incompressible Navier-Stokes equations, Math. Comp. 74 (005), [8], A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations, J. Sci. Comput. 31 (007), [9] B. Cockburn, G. anschat, D. Schötzau, and C. Schwab, Local discontinuous Galerkin methods for the Stokes system, SIAM J. Numer. Anal. 40 (00), [10] J. Gopalakrishnan and G. anschat, Application of unified DG analysis to preconditioning DG methods, Proceedings of the Second M.I.T. Conference on Computational Fluid and Solid Mechanics, Cambridge, MA, USA (.J. Bathe, ed.), Elsevier, 003, pp [11] G. anschat and D. Schötzau, Energy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations, Int. J. Numer. Meth. Fluids 57 (008), [1] L.I.G. ovasznay, Laminar flow behind a two-dimensional grid, Proc. Camb. Philos. Soc. 44 (1948),
27 [13] P. Lesaint and P. A. Raviart, On a finite element method for solving the neutron transport equation, Mathematical Aspects of Finite Elements in Partial Differential Equations (C. de Boor, ed.), Academic Press, 1974, pp [14] I. Perugia and D. Schötzau, An hp-analysis of the local discontinuous Galerkin method for diffusion problems, J. Sci. Comput. 17 (00), [15] A. Quarteroni and A. Valli, Numerical Approximation of Partial Differential Equations, Springer Verlag, New York, [16] D. Schötzau, C. Schwab, and A. Toselli, Stabilized hp-dgfem for incompressible flow, Math. Models Methods Appl. Sci. 13 (003),
An Equal-Order DG Method for the Incompressible Navier-Stokes Equations
J Sci Comput 009) 40: 188 10 DOI 10.1007/s10915-008-961-1 An qual-order DG Method for the Incompressible Navier-Stokes quations Bernardo Cockburn Guido anschat Dominik Schötzau Received: 4 September 008
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 informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationKey words. Incompressible magnetohydrodynamics, mixed finite element methods, discontinuous Galerkin methods
A MIXED DG METHOD FOR LINEARIZED INCOMPRESSIBLE MAGNETOHYDRODYNAMICS PAUL HOUSTON, DOMINIK SCHÖTZAU, AND XIAOXI WEI Journal of Scientific Computing, vol. 40, pp. 8 34, 009 Abstract. We introduce and analyze
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 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 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 informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More 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 informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationA 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 informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationA 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 PROJECTION-BASED ERROR ANALYSIS OF HDG METHODS BERNARDO COCKBURN, JAYADEEP GOPALAKRISHNAN, AND FRANCISCO-JAVIER SAYAS
A PROJECTION-BASED ERROR ANALYSIS OF HDG METHODS BERNARDO COCBURN, JAYADEEP GOPALARISHNAN, AND FRANCISCO-JAVIER SAYAS Abstract. We introduce a new technique for the error analysis of hybridizable discontinuous
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationA Hybrid Discontinuous Galerkin Method 2 for Darcy-Stokes Problems 3 UNCORRECTED PROOF
1 A Hybrid Discontinuous Galerkin Method 2 for Darcy-Stokes Problems 3 Herbert Egger 1 and Christian Waluga 2 4 1 Center of Mathematics, Technische Universität München, Boltzmannstraße 3, 85748 5 Garching
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 informationIMA Preprint Series # 2146
AN ANALYSIS OF THE MINIMAL DISSIPATION LOCAL DISCONTINUOUS GALERKIN METHOD FOR CONVECTION-DIFFUSION PROBLEMS By Bernardo Cockburn and Bo Dong IMA Preprint Series # 2146 ( November 2006 ) INSTITUTE FOR
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 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 informationANALYSIS OF AN INTERFACE STABILISED FINITE ELEMENT METHOD: THE ADVECTION-DIFFUSION-REACTION EQUATION
ANALYSIS OF AN INTERFACE STABILISED FINITE ELEMENT METHOD: THE ADVECTION-DIFFUSION-REACTION EQUATION GARTH N. WELLS Abstract. Analysis of an interface stabilised finite element method for the scalar advectiondiffusion-reaction
More informationDiscontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
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 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 informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationHybridized DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More 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 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 informationA posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More 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 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 information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
More informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationA connection between grad-div stabilized FE solutions and pointwise divergence-free FE solutions on general meshes
A connection between grad-div stabilized FE solutions and pointwise divergence-free FE solutions on general meshes Sarah E. Malick Sponsor: Leo G. Rebholz Abstract We prove, for Stokes, Oseen, and Boussinesq
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More 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 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 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 informationA SIMPLE TUTORIAL ON DISCONTINUOUS GALERKIN METHODS. Jennifer Proft CERMICS, ENPC. J. Proft CERMICS, ENPC
A SIMPLE TUTORIAL ON DISCONTINUOUS GALERKIN METHODS Jennifer Proft Outline Definitions A simple example Issues Historical development elliptic equations hyperbolic equations Discontinuous vs. continuous
More informationA MIXED FINITE ELEMENT METHOD WITH EXACTLY DIVERGENCE-FREE VELOCITIES FOR INCOMPRESSIBLE MAGNETOHYDRODYNAMICS
A MIXED FINITE ELEMENT METHOD WITH EXACTLY DIVERGENCE-FREE VELOCITIES FOR INCOMPRESSIBLE MAGNETOHYDRODYNAMICS CHEN GREIF, DAN LI, DOMINIK SCHÖTZAU, AND XIAOXI WEI Abstract. We introduce and analyze a mixed
More informationSTABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS
STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS RAYTCHO LAZAROV AND XIU YE Abstract. In this paper, we derive two stabilized discontinuous finite element formulations, symmetric
More 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 informationGeneralized Newtonian Fluid Flow through a Porous Medium
Generalized Newtonian Fluid Flow through a Porous Medium V.J. Ervin Hyesuk Lee A.J. Salgado April 24, 204 Abstract We present a model for generalized Newtonian fluid flow through a porous medium. In the
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS ERIK BURMAN AND BENJAMIN STAMM Abstract. In this paper we give an analysis of a bubble stabilized discontinuous Galerkin
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
More informationA posteriori error estimates for a Maxwell type problem
Russ. J. Numer. Anal. Math. Modelling, Vol. 24, No. 5, pp. 395 408 (2009) DOI 0.55/ RJNAMM.2009.025 c de Gruyter 2009 A posteriori error estimates for a Maxwell type problem I. ANJAM, O. MALI, A. MUZALEVSKY,
More informationarxiv: v1 [math.na] 15 Nov 2017
Noname manuscript No. (will be inserted by the editor) An HDG method with orthogonal projections in facet integrals Issei Oikawa arxiv:1711.05396v1 [math.na] 15 Nov 2017 Received: date / Accepted: date
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationAn hp-adaptive Mixed Discontinuous Galerkin FEM for Nearly Incompressible Linear Elasticity
An hp-adaptive Mixed Discontinuous Galerkin FEM for Nearly Incompressible Linear Elasticity Paul Houston 1 Department of Mathematics, University of Leicester, Leicester LE1 7RH, UK (email: Paul.Houston@mcs.le.ac.uk)
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationICES REPORT Analysis of the DPG Method for the Poisson Equation
ICES REPORT 10-37 September 2010 Analysis of the DPG Method for the Poisson Equation by L. Demkowicz and J. Gopalakrishnan The Institute for Computational Engineering and Sciences The University of Texas
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 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 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 informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
More informationHybrid Discontinuous Galerkin methods for incompressible flow problems
Hybrid Discontinuous Galerkin methods incompressible flow problems Christoph Lehrenfeld, Joachim Schöberl MathCCES Computational Mathematics in Engineering Workshop Linz, May 31 - June 1, 2010 Contents
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationA-posteriori error analysis of hp-version discontinuous Galerkin finite element methods for second-order quasilinear elliptic problems
Report no. NA-0/1 A-posteriori error analysis of hp-version discontinuous Galerkin finite element methods for second-order quasilinear elliptic problems Paul Houston 1, Endre Süli, and Thomas P. Wihler
More informationA mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N)
wwwijmercom Vol2, Issue1, Jan-Feb 2012 pp-464-472 ISSN: 2249-6645 A mixed finite element approximation of the Stokes equations with the boundary condition of type (D+N) Jaouad El-Mekkaoui 1, Abdeslam Elakkad
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationNew DPG techniques for designing numerical schemes
New DPG techniques for designing numerical schemes Jay Gopalakrishnan University of Florida Collaborator: Leszek Demkowicz October 2009 Massachusetts Institute of Technology, Boston Thanks: NSF Jay Gopalakrishnan
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationhp-version Discontinuous Galerkin Finite Element Methods for Semilinear Parabolic Problems
Report no. 3/11 hp-version Discontinuous Galerkin Finite Element Methods for Semilinear Parabolic Problems Andris Lasis Endre Süli We consider the hp version interior penalty discontinuous Galerkin finite
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More 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 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 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 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 informationSTABILIZED HP-DGFEM FOR INCOMPRESSIBLE FLOW. Math. Models Methods Appl. Sci., Vol. 13, pp , 2003
Mathematical Models and Methods in Applied Sciences c World Scientiic Publishing Company STABILIZED HP-DGFEM FOR INCOMPRESSIBLE FLOW D. SCHÖTZAU1 C. SCHWAB 2 A. TOSELLI 2 1 Department o Mathematics, University
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 informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More 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 informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationA Matlab Tutorial for Diffusion-Convection-Reaction Equations using DGFEM
MIDDLE EAST TECHNICAL UNIVERSITY Ankara, Turkey INSTITUTE OF APPLIED MATHEMATICS http://iam.metu.edu.tr arxiv:1502.02941v1 [math.na] 10 Feb 2015 A Matlab Tutorial for Diffusion-Convection-Reaction Equations
More informationA DISCONTINUOUS GALERKIN METHOD FOR ELLIPTIC INTERFACE PROBLEMS WITH APPLICATION TO ELECTROPORATION
A DISCONTINUOUS GALERIN METHOD FOR ELLIPTIC INTERFACE PROBLEMS WITH APPLICATION TO ELECTROPORATION GRÉGORY GUYOMARC H AND CHANG-OC LEE Abstract. We present a discontinuous Galerkin (DG) method to solve
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More information