SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS

Size: px
Start display at page:

Download "SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS"

Transcription

1 Proceedings of ALGORITMY 2009 pp SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar nonstationary convection-diffusion equation with a nonlinear convection and a linear diffusion. We carry out the space semi-discretization with the aid of the nonsymmetric interior penalty Galerkin (NIPG) method and the time discretization by the time discontinuous Galerkin method linearized by extrapolation from previous time interval. The resulting scheme is unconditionally stable, has a high order of accuracy with respect to space and time coordinates and requires only solutions of linear algebraic problems at each time step. We derive a priori error estimate in the L 2 -norm. Key words. a priori error estimate, time discontinuous Galerkin method 1. Introduction. We numerically solve a nonstationary nonlinear convectiondiffusion equation, which represents a model problem for the system of the compressible Navier-Stokes equations. The class of discontinuous Galerkin (DG) methods seems to be one of the most promising candidates to construct high order accurate schemes for solving convection-diffusion problems. For a survey about DG methods, see [1] or [3]. An analysis of DG methods was presented in many papers, see, e.g. [6], [8], [11], [10]. In [8] we carried out the space semi-discretization of the scalar convection-diffusion equation with the aid of the discontinuous Galerkin finite element method and derived a priori error estimates. Within this contribution, we deal with the time discretization of the resulting system of ordinary differential equations. In contrary to [9], where we used the so-called backward difference formulae (BDF) approach, here we employ the time discontinuous Galerkin method, see [15]. Since this scheme is implicit and we would like to avoid solving strongly nonlinear problem at each time level, we recomend linearization in a similar way as it was done before for BDF in [9]. We present a formulation of the arbitrary order linearized time discontinuous Galerkin scheme and derive a priori error estimate for the second order scheme. 2. Continuous problem. Let Ω R d (d = 2 or 3) be a bounded polyhedral domain and T > 0. We set Q T = Ω (0, T ). By Ω and Ω we denote the closure and boundary of Ω, respectively. Let us consider the following initial-boundary value problem: Find u : Q T R such that (1) (2) (3) u t + f(u) = u + g in Q T, u Ω (0,T ) = u D, u(x, 0) = u 0 (x), x Ω. In (1) (3), f = (f1,..., f d ), f s C 2 (R), f s (0) = 0, s = 1,..., d represents convective terms, g C([0, T ]; L 2 (Ω)) represents volume sources. The Dirich- This work is partly supported by the Grant No. 316/2006/B-MAT/MFF of the Grant Agency of the Charles University Prague and also by the project LC06052 financed by MSMT (Necas Center for Mathematical Modeling). Charles University Prague, Faculty of Mathematics and Physics, Sokolovská 83, Prague, Czech Republic, (vlasakmila@gmail.com) 1

2 2 M. VLASAK let boundary condition is given over Ω by u D, which is the trace of some u C([0, T ]; H 1 (Ω)) L (Q T ) and u 0 L 2 (Ω) is an initial condition. We use the standard notation for Lebesgue, Sobolev and Bochner function spaces (see, e. g. [12]). In order to introduce the concept of a weak solution, we define (u, w) = uw dx, u, w L 2 (Ω), Ω a(u, w) = u w dx, u, w H 1 (Ω), Ω b(u, w) = f(u) w dx, u H 1 (Ω) L (Ω), w L 2 (Ω), Ω Definition 2.1. We say that a function u is a weak solution of (1) (3) if the following conditions are satisfied (4) a) u u L 2 (0, T ; H0 1 (Ω)), u L (Q T ), d b) (u(t), w) + b(u(t), w) + a(u(t), w) = (g(t), w) dt for all w H0 1 (Ω) in the sense of distributions on (0, T ), c) u(0) = u 0 in Ω. By u(t) we denote the function on Ω such that u(t)(x) = u(x, t), x Ω. With the aid of techniques from [13] and [14], it is possible to prove that there exists a unique weak solution. We shall assume that the weak solution u is sufficiently regular, namely, (5) u W 1, (0, T ; H p+1 (Ω)) W 2, (0, T ; H 1 (Ω)), an integer p 1 will denote a given degree of polynomial approximations. Such a solution satisfies problem (1) (3) pointwise. 3. Space semi-discretization. We discretize problem (4) in space with the aid of the discontinuous Galerkin finite element method with nonsymmetric treatment of stabilization terms and interior and boundary penalties. This approach is called the NIPG variant of the DGFE method, see [1]. We derived the space discretization of (1) (3) by the NIPG variant of DGFE method in [8] hence here we present only the final expressions Triangulation. Let T h (h > 0) be a partition of the domain Ω into a finite number of closed d-dimensional mutually disjoint simplices K i.e., Ω = K T h K. By K we denote the boundary of element K T h and set h K = diam(k), h = max K Th h K. We set Γ the faces of T h (Γ = K T h K). By ρ K we denote the radius of the largest d-dimensional ball inscribed into K and by K we denote the d-dimensional Lebesgue measure of K. Furthermore, we use the following notation: n = (n 1,..., n d ) a normal vector to Γ which is well defined almost everywhere (on Ω we use outer normal, inside of Ω we use one (arbitrary but fixed) direction at every point of Γ) Broken Sobolev spaces. We define the so-called broken Sobolev space in the following way (6) H s (Ω, T h ) = {w; w K H s (K) K T h }

3 ALGORITMY 2009 TIME DG 3 and define there the norm (7) ( w H s (Ω,T h ) = K T h w 2 H s (K) ) 1/2 and the seminorm (8) ( w Hs (Ω,T h ) = K T h w 2 H s (K) ) 1/2. For w H 1 (Ω, T h ), we introduce the following notation on Γ \ Ω: (9) w R (x) = lim w(x + δn) δ 0+ wl(x) = lim w(x + δn) δ 0 w = 1 2 (w R + w L ), [w] = w L w R and on Ω we put (10) w L (x) = lim w(x + δn) δ 0 w = wl, [w] = w L 3.3. Space discretization. For u, w H 2 (Ω, T h ) we set A h (u, w) = ( ) (11) u w dx u n[w] w n[u] ds (12) (13) K T h + b h (u, w) = Γ Γ\ Ω K T h K σ[u] [w] ds H(u L, u R, n) [w] ds + K l h (w) (t) = (g(t), w) Γ Ω H(u L, u D, n) w L ds f(u) w dx, u, w H 1 (Ω, T h ), u L (Ω) Ω ( w n u D (t) σu D (t) w) ds. The penalty parameter function σ in (11) and (13) along the face e Γ is defined by σ e = 1/(h K + h K), e = K K. The function H(,, ) in the face integrals in (12) is called the numerical flux, well-known from the finite volume method and it approximates the terms f(u) n. Now we define the space of discontinuous piecewise polynomial functions (14) S h = S p, 1 (Ω, T h ) = {w; w K P p (K) K T h }, where P p (K) denotes the space of all polynomials on K of degree p, where the integer p 1 is a given degree of approximation.

4 4 M. VLASAK We find that the exact solution of (4) with property (5) satisfies the identity ( ) u (15) t (t), w h + A h (u(t), w h ) + b h (u(t), w h ) = l h (w h ) (t) for all w h S h and all t (0, T ). The (semi)-discrete problem (15) represents a system of ordinary differential equations (ODEs) which is solved by a suitable solver in the next section. 4. Time discretization. Since problem (15) is stiff, it is necessary to solve it with a method having a large stability domain. In [9] we employ the well known backward difference formulae (BDF). Within this contribution, we present a new approach based on the time discontinuous Galerkin method. Since this method is implicit, we introduce linearization from previous time interval in such a way that linear part of the problem is treated implicitly and the nonlinear one explicitly. Therefore, the resulting semi-implicit scheme leads to a highly stable method which requires only a solution of the linear algebraic problem at each time step. We consider a partition 0 = t 0 < t 1 <... < t r = T consisting of time intervals = (t m 1, t m ], m = 1,..., r of the length = τ m and τ = max m=1,...,r τ m. We set the space time fully discrete space (16) q S h,τ {w L 2 (Q T ) : w Im = t s z s, z s S h }. s=0 For simplicity we set (17) (18) w± m = lim w(t m + δ), δ 0± {w} m = w+ m w m. To avoid the nonlinearity in our problem we define ŵ Im as prolongation of w from previous time interval 1, it means there exists z P q ( 1 ) such that (19) z(t) = w Im 1 (t), t 1, ŵ Im (t) = z(t), t. Now we are able to define fully discrete solution U S h,τ. Definition 4.1. We define the approximate solution of problem from Definition (2.1) a) c) as functions U S h,τ satisfying the conditions a) (U, w) + A h (U, w) + b h (Û, w)dt + ({U} m 1, w+ m 1 ) = l h (w)dt (20) w S h,τ, m = 2,..., r, b)(u 0, w) = (u 0, w) w S h. We should remark on some strange character of this scheme. This scheme is one step method as it is usual for time discontinuous Galerkin, but on the other hand we need to use the information from the whole previous time interval for construction of the prolongation, so this method is not self started.

5 ALGORITMY 2009 TIME DG 5 5. Error estimates. Our goal is to analyse the error estimates of the approximate solution U obtained by the method (20) for the piecewise linear approximation in time (q = 1). In the next we set (21) w 2 = A h (w, w) w H 2 (Ω, T h ) Parabolic projection. For the next purpose we need projection π : L 2 (Q T ) S h,τ well known for the time discontinuous Galerkin. This projection satisfies (22) (πv v, t j w)dt = 0 w S h, j = 0,..., q 1, πv m = Πv m, where Π be the standard L 2 (Ω) S h orthogonal projection. About the projection π it is possible to prove following estimates Lemma 5.1. Let w W q+1, (, S h ), then (23) w(t) πw(t) L 2 (Ω) Cτ q+1 m, t, w(t) πw(t) Cτ q+1 m, t, where C depends on time derivatives of w. Proof. The lemma can be prooved by standard scaling argument. Since we have simplicial mesh we obtain standard estimates for projection Π under sufficient regularity of function w (24) w Πw L 2 (Ω) Ch p+1 w Πw Ch p, we can estimate projection error of arbitrary function w W q+1, (, H 1 (Ω)) W 1, (, H p+1 (Ω)) by (25) w(t) πw(t) L 2 (Ω) C(h p+1 + τ q+1 m ), t, w(t) πw(t) C(h p + τ q+1 m ), t Extrapolation. We set ŵ as extrapolation from previous time interval. For q = 1 we get (26) ŵ(t) = (1 + t t m 1 )w m 1 t t m 1 w+ m 2 t. As standard result from interpolation theory with additional assumption (27) τ m C and under sufficient regularity w W 2, (0, T, H 1 (Ω)) we get (28) w ŵ L2 (Ω) Cτ 2, w ŵ Cτ 2, where constant C depends on time derivatives of w. This can be prooved by standard scaling argument.

6 6 M. VLASAK 5.3. Properties of the forms A h and b h. Here we want to summarize the properties of the forms A h and b h. Lemma 5.2. Let u W 1, (0, T, H p+1 (Ω)). Then it holds that A h (u(t) Πu(t), w) Ch p w w S h, A h (v, w) C v w v, w S h. A h (u(t) πu(t), w) C(h p + τ q+1 ) w w S h Proof. The proof can be found in [8] in Lemma 9. Lemma 5.3. Let u W 1, (0, T, H p+1 (Ω)). Let the numerical fluxes H be Lipschitz continuous, conservative and consistent. Then it holds that b h (v, w) b h ( v, w) C( v v L2 (Ω) + v v ) w v, v L 2 (Ω), w S h, b h (u, w) b h (Πu, w) Ch p+1 w w S h, b h (v, w) b h ( v, w) C v v L2 (Ω) w v, v, w S h, b h (u, w) b h (û, w) Cτ 2 w w S h. Proof. The proof can be found in [8] in Lemma 5 and Lemma Main result. In the sequel we use the notation U u = U πu + πu u = ξ + η. Theorem 5.4. Let u be the exact solution of problem (4) satisfying (5). Let the mesh be regular (h K /ρ K C), assumption (27) be satisfied and the numerical fluxes H be Lipschitz continuous, conservative and consistent. Let U be the approximate solution defined by (20). Then (29) max U s u(t s ) 2 C(h 2p + τ 4 + e 1 2 s=0,...,r L 2 (Ω) + e0 2 L 2 (Ω) )et C Proof. Since U s u(t s ) = ξ s +η s, in virtue of (25), it is only sufficient to estimate ξ s. Let us integrate (4) over and substract this equation from (20). Then we have ( ) (ξ, w) + A h (ξ, w)dt + ({ξ} m 1, w+ m 1 ) = (η, w)dt + ({η} m 1, w+ m 1 ) I m (30) A h (η, w)dt + b h (u, w) b h (Û, w)dt w S h,τ. Now let us have a look at all the parts of (30). (31) (η, w)dt + ({η} m 1, w+ m 1 ) = (η m, w m ) (η+ m 1, w+ m 1 ) (η, v )dt + ({η} m 1, w m 1 I m = (η m, w m ) (η m 1, w+ m 1 ) (η, w )dt Using properties (22) of π we obtain (32) (η m, w m ) (η m 1, w+ m 1 ) (η, w )dt = 0. + )

7 ALGORITMY 2009 TIME DG 7 Estimating linear form A h using Lemma 5.2 we obtain A h (η, w)dt 1 (33) w 2 dt + τ m C(h 2p + τ 4 ). 4 Now we should estimate nonlinear form b h. (34) b h (u, v) b h (Û, w) b h(u, w) b h (û, w) + b h (û, w) b h (Πû, w) + b h (Πû, w) b h (Û, w) We estimate the individual terms using Lemma 5.3. (35) b h (u, w) b h (û, w) Cτ 2 w b h (û, w) b h (Πû, w) Ch p+1 w b h (Πû, w) b h (Û, w) C Πû Û L 2 (Ω) w With the aid of the definition of extrapolation of û and Û we get (36) Πû Û L 2 (Ω) ( = 1 + t t ) m 1 ( 1 + t t ) n 1 τ n 1 (Πu m 1 U m 1 ξ m 1 L2 (Ω) + t t n 1 ) t t m 1 τ n 1 (Πu m 2 U m 2 + ) L 2 (Ω) Πu m 2 U m 2 + L2 (Ω). We estimate last norm by (24) and (25) (37) Πu m 2 U+ m 2 L2 (Ω) Πu m 2 u m 2 L2 (Ω) + u m 2 πu m 2 + L2 (Ω) + πu m 2 + U+ m 2 L2 (Ω) Ch p+1 + C(h p+1 + τ 2 ) + ξ+ m 2 L 2 (Ω) C(h p+1 + τ 2 ) + {ξ} m 2 L2 (Ω) + ξ m 2 L2 (Ω). When we use assumption (27), we could complete estimate for form b h b h (u, w) b h (Û, w)dt 1 (38) w 2 dt + τ m C(h 2p+2 + τ 4 ) 4 +τ m C( ξ m 1 2 L 2 (Ω) + {ξ} m 2 2 L 2 (Ω) + ξm 2 2 L 2 (Ω) ). When we apply w = 2ξ to the left-hand side of (30), we obtain (39) 2(ξ, ξ) + 2 ξ 2 dt + 2({ξ} m 1, ξ+ m 1 ) = ξ m 2 L 2 (Ω) ξm L 2 (Ω) + 2 ξm L 2 (Ω) 2(ξm 1, ξ+ m 1 ) + 2 ξ 2 dt = ξ m 2 L 2 (Ω) ξm 1 2 L 2 (Ω) + {ξ} m 1 2 L 2 (Ω) + 2 ξ 2 dt. If we apply all these estimates with w = 2ξ together we gain (40) ξ m 2 L 2 (Ω) ξm 1 2 L 2 (Ω) + {ξ} m 1 2 L 2 (Ω) τ m q(h, τ) + τ m C( ξ m 1 2 L 2 (Ω) + {ξ} m 2 2 L 2 (Ω) + ξm 2 2 L 2 (Ω) ),

8 8 M. VLASAK where q(h, τ) = O(h 2p + τ 4 ). Summing over m = 2,..., k r and using (27), we gain (41) ξ k 2 L 2 (Ω) ξ1 2 L 2 (Ω) + {ξ} k 1 2 L 2 (Ω) k 1 ( ) T q(h, τ) + C τ s+1 ξ s 2 L 2 (Ω) + {ξ} s 1 2 L 2 (Ω) + ξs 1 s=1 2 L 2 (Ω) k 1 T q(h, τ) + τ 2 C ξ 0 2 L 2 (Ω) + C (τ s+1 + τ s+2 ) ξ s 2 L 2 (Ω) + τ s+1 {ξ} s 1 2 L 2 (Ω) s=1 k 1 T q(h, τ) + τ 2 C ξ 0 2 L 2 (Ω) + C τ s+1 ξ s 2 L 2 (Ω) + τ s+1 {ξ} s 1 2 L 2 (Ω). Applying Gronwall s lemma we obtain ( ) k ξ k 2 L 2 (Ω) + {ξ} k 1 2 L 2 (Ω) C ξ 0 2 L 2 (Ω) + ξ1 2 L 2 (Ω) + T q (42) (1 + τ s C) which proves our theorem. s=1 s=2 (C ξ 0 2 L 2 (Ω) + ξ1 2 L 2 (Ω) + T q ) e T C, REFERENCES [1] 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(5): , [2] P. G. Ciarlet. The Finite Elements Method for Elliptic Problems. North-Holland, Amsterdam, New York, Oxford, [3] B. Cockburn, G. E. Karniadakis, and C.-W. Shu (eds.). Discontinuous Galerkin methods. In Lecture Notes in Computational Science and Engineering 11. Springer, Berlin, [4] V. Dolejší, M. Feistauer. Error estimates of the discontinuous Galerkin method for nonlinear nonstationary convection-diffusion problems. Numer. Func. Anal. Optimiz. 26(25-26): , [5] V. Dolejší, M. Feistauer, J. Hozman: Analysis of Semi-Implicit DGFEM for Nonlinear Convection-Diffusion Problems, Comput. Methods Appl. Mech. Engrg. (submitted). [6] V. Dolejší, M. Feistauer, and V. Kučera, V. Sobotíková. An optimal L (L 2 )-error estimate of the discontinuous Galerkin method for a nonlinear nonstationary convection-diffusion problem. IMA J. Numer. Anal., (submitted). [7] V. Dolejší, M. Feistauer, and C. Schwab. A finite volume discontinuous Galerkin scheme for nonlinear convection diffusion problems. Calcolo, 39:1 40, [8] V. Dolejší, M. Feistauer, and V. Sobotíková. Analysis of the Discontinuous Galerkin Method for Nonlinear Convection Diffusion Problems, Comput. Methods Appl. Mech. Eng. 194: , [9] V. Dolejší and M. Vlasák. Analysis of a BDF DGFE Scheme for Nonlinear Convection- Diffusion Problems, Numer. Math. 110: , [10] P. Houston, J. Robson, and E. Süli. Discontinuous hp-finite element methods for advectiondiffusion problems. SIAM J. Numer. Anal., 39(6): , [11] P. Houston, C. Schwab, and E. Süli. Discontinuous Galerkin finite element approximation of quasilinear elliptic boundary value problems I: The scalar case. IMA J. Numer. Anal., 25: , [12] A. Kufner, O. John, and S. Fučík. Function Spaces. Academia, Prague, [13] P. L. Lions. Mathematical Topics in Fluid Mechanics. Oxford Science Publications, [14] K. Rektorys. The Method of Discretization in Time and Partial Differential Equations. Reidel, Dodrecht, [15] V. Thomée. Galerkin finite element methods for parabolic problems. 2nd revised and expanded ed.. Spriger, Berlin, 2006.

RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS

RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS Proceedings of ALGORITMY 2016 pp. 113 124 RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS VÍT DOLEJŠÍ AND FILIP ROSKOVEC Abstract.

More information

ANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION

ANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION Proceedings of EQUADIFF 2017 pp. 127 136 ANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION MILOSLAV FEISTAUER, ONDŘEJ BARTOŠ, FILIP ROSKOVEC, AND ANNA-MARGARETE

More information

AN OPTIMAL UNIFORM A PRIORI ERROR ESTIMATE FOR AN UNSTEADY SINGULARLY PERTURBED PROBLEM

AN OPTIMAL UNIFORM A PRIORI ERROR ESTIMATE FOR AN UNSTEADY SINGULARLY PERTURBED PROBLEM INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 11, Number 1, Pages 24 33 c 2014 Institute for Scientific Computing and Information AN OPTIMAL UNIFORM A PRIORI ERROR ESTIMATE FOR AN UNSTEADY

More information

Different Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method

Different 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 information

Local discontinuous Galerkin methods for elliptic problems

Local 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 information

ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS

ELLIPTIC 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 information

NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL

NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper

More information

Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons

Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées

More information

1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying

1. 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 information

NUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM

NUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM Proceedings of ALGORITMY 212 pp. 44 415 NUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM NAVEED AHMED AND GUNAR MATTHIES Abstract. In this paper, we combine

More information

A space-time Trefftz method for the second order wave equation

A space-time Trefftz method for the second order wave equation A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil

More information

NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE

NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING

More information

A space-time Trefftz method for the second order wave equation

A space-time Trefftz method for the second order wave equation A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of

More information

hp-version Discontinuous Galerkin Finite Element Methods for Semilinear Parabolic Problems

hp-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 information

Mixed Discontinuous Galerkin Methods for Darcy Flow

Mixed Discontinuous Galerkin Methods for Darcy Flow Journal of Scientific Computing, Volumes 22 and 23, June 2005 ( 2005) DOI: 10.1007/s10915-004-4150-8 Mixed Discontinuous Galerkin Methods for Darcy Flow F. Brezzi, 1,2 T. J. R. Hughes, 3 L. D. Marini,

More information

Discontinuous Galerkin Methods: Theory, Computation and Applications

Discontinuous 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 information

TWO DIMENSIONAL SIMULATION OF FLUID-STRUCTURE INTERACTION USING DGFEM

TWO DIMENSIONAL SIMULATION OF FLUID-STRUCTURE INTERACTION USING DGFEM Proceedings of ALGORITMY 212 pp. 85 94 TWO DIMENSIONAL SIMULATION OF FLUID-STRUCTURE INTERACTION USING DGFEM JAROSLAVA HASNEDLOVÁ-PROKOPOVÁ,, MILOSLAV FEISTAUER, JAROMÍR HORÁČEK, ADAM KOSÍK, AND VÁCLAV

More information

On discontinuity capturing methods for convection diffusion equations

On discontinuity capturing methods for convection diffusion equations On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,

More information

On 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 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 information

DGFEM. Theory and Applications to CFD

DGFEM. Theory and Applications to CFD DGFEM. Theory and Applications to CFD Miloslav Feistauer Charles University Prague in cooperation with V. Dolejší, V. Kučera, V. Sobotíková and J. Prokopov á Presented at Humboldt-Universität, July 2009

More information

CONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS

CONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS J. Korean Math. Soc. 34 (1997), No. 3, pp. 515 531 CONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS S. K. CHUNG, A.K.PANI AND M. G. PARK ABSTRACT. In this paper,

More information

Fourier Type Error Analysis of the Direct Discontinuous Galerkin Method and Its Variations for Diffusion Equations

Fourier Type Error Analysis of the Direct Discontinuous Galerkin Method and Its Variations for Diffusion Equations J Sci Comput (0) 5:68 655 DOI 0.007/s095-0-9564-5 Fourier Type Error Analysis of the Direct Discontinuous Galerkin Method and Its Variations for Diffusion Equations Mengping Zhang Jue Yan Received: 8 September

More information

A 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 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 information

Superconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients

Superconvergence 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 information

A 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 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 information

An Equal-order DG Method for the Incompressible Navier-Stokes Equations

An Equal-order DG Method for the Incompressible Navier-Stokes Equations An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce

More information

A-posteriori error analysis of hp-version discontinuous Galerkin finite element methods for second-order quasilinear elliptic problems

A-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 information

STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS

STABILIZED 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 information

A local-structure-preserving local discontinuous Galerkin method for the Laplace equation

A 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 information

Superconvergence and H(div) Projection for Discontinuous Galerkin Methods

Superconvergence and H(div) Projection for Discontinuous Galerkin Methods INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 00:1 6 [Version: 2000/07/27 v1.0] Superconvergence and H(div) Projection for Discontinuous Galerkin Methods Peter Bastian

More information

Optimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method

Optimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method Report no. 03/06 Optimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method Emmanuil H. Georgoulis and Endre Süli Oxford University Computing Laboratory Numerical

More information

arxiv: v1 [math.na] 29 Feb 2016

arxiv: 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 information

On an Approximation Result for Piecewise Polynomial Functions. O. Karakashian

On 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 information

arxiv: v1 [math.na] 15 Nov 2017

arxiv: 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 information

Adaptive backward difference formula discontinuous Galerkin finite element method for the solution of conservation laws

Adaptive backward difference formula discontinuous Galerkin finite element method for the solution of conservation laws Adaptive backward difference formula discontinuous Galerkin finite element method for the solution of conservation laws V. Dolejší and P. Kůs Charles University Prague, Faculty of Mathematics and Physics,

More information

PARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS

PARTITION 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 information

Interaction of Incompressible Fluid and Moving Bodies

Interaction of Incompressible Fluid and Moving Bodies WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics

More information

arxiv: v1 [math.na] 27 Jan 2016

arxiv: 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 information

Yongdeok Kim and Seki Kim

Yongdeok 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 information

AMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50

AMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50 A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in

More information

Discontinuous Galerkin Methods

Discontinuous 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 information

A SIMPLE TUTORIAL ON DISCONTINUOUS GALERKIN METHODS. Jennifer Proft CERMICS, ENPC. J. Proft CERMICS, ENPC

A 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 information

BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM

BUBBLE 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 information

BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS

BUBBLE 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 information

1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods

1 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 information

Discrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem

Discrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant

More information

WEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES

WEAK 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 information

WEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS

WEAK 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 information

A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS

A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS Z. LU Communicated by Gabriela Marinoschi In this paper, we discuss a

More information

ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen

ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department

More information

On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions

On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions Proceedings of the 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August -, 5 pp36-41 On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability

More information

A Two-grid Method for Coupled Free Flow with Porous Media Flow

A 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 information

Hybridized Discontinuous Galerkin Methods

Hybridized 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 information

A Mixed Nonconforming Finite Element for Linear Elasticity

A 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 information

GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS

GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND CHARALAMBOS MAKRIDAKIS Abstract. We consider discontinuous as well as continuous Galerkin methods for the time discretization

More information

A Posteriori Error Estimation for hp-version Time-Stepping Methods for Parabolic Partial Differential Equations

A Posteriori Error Estimation for hp-version Time-Stepping Methods for Parabolic Partial Differential Equations Numerische Mathematik manuscript No. (will be inserted by the editor) A Posteriori Error Estimation for hp-version Time-Stepping Methods for Parabolic Partial Differential Equations Dominik Schötzau 1,

More information

Finite Element Clifford Algebra: A New Toolkit for Evolution Problems

Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/

More information

A Least-Squares Finite Element Approximation for the Compressible Stokes Equations

A 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 information

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,

More information

Equations paraboliques: comportement qualitatif

Equations paraboliques: comportement qualitatif Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Equations paraboliques: comportement qualitatif par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 25/6 1 Contents

More information

AN EQUILIBRATED A POSTERIORI ERROR ESTIMATOR FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD

AN EQUILIBRATED A POSTERIORI ERROR ESTIMATOR FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD AN EQUILIBRATED A POSTERIORI ERROR ESTIMATOR FOR THE INTERIOR PENALTY DISCONTINUOUS GALERIN METHOD D. BRAESS, T. FRAUNHOLZ, AND R. H. W. HOPPE Abstract. Interior Penalty Discontinuous Galerkin (IPDG) methods

More information

ANALYSIS 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 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 information

Discontinuous Galerkin methods for fractional diffusion problems

Discontinuous Galerkin methods for fractional diffusion problems Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation

More information

1. Introduction. We consider the model problem: seeking an unknown function u satisfying

1. 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 information

Numerical Methods for the Navier-Stokes equations

Numerical 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 information

PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this

More information

An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element

An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University

More information

Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations

Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,

More information

Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems

Analysis 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 information

A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations

A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for

More information

Piecewise Smooth Solutions to the Burgers-Hilbert Equation

Piecewise Smooth Solutions to the Burgers-Hilbert Equation Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang

More information

STABILIZED HP-DGFEM FOR INCOMPRESSIBLE FLOW. Math. Models Methods Appl. Sci., Vol. 13, pp , 2003

STABILIZED 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 information

A posteriori error estimates for non conforming approximation of eigenvalue problems

A 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 information

Lecture Note III: Least-Squares Method

Lecture 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 information

ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS

ENERGY 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 information

b i (x) u + c(x)u = f in Ω,

b 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 information

A posteriori error estimates of fully discrete finite element method for Burgers equation in 2-D

A posteriori error estimates of fully discrete finite element method for Burgers equation in 2-D A posteriori error estimates of fully discrete finite element method for Burgers equation in -D Tong Zhang a,b, Jinyun Yuan b, Damazio Pedro b a School of Mathematics & Information Science, Henan Polytechnic

More information

Energy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations

Energy 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 information

A 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 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 information

A Matlab Tutorial for Diffusion-Convection-Reaction Equations using DGFEM

A 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 information

Applied Mathematics Letters. Nonlinear stability of discontinuous Galerkin methods for delay differential equations

Applied Mathematics Letters. Nonlinear stability of discontinuous Galerkin methods for delay differential equations Applied Mathematics Letters 23 21 457 461 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Nonlinear stability of discontinuous Galerkin

More information

ETNA Kent State University

ETNA Kent State University Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.

More information

A DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR TIME DEPENDENT PARTIAL DIFFERENTIAL EQUATIONS WITH HIGHER ORDER DERIVATIVES

A DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR TIME DEPENDENT PARTIAL DIFFERENTIAL EQUATIONS WITH HIGHER ORDER DERIVATIVES MATHEMATCS OF COMPUTATON Volume 77, Number 6, April 008, Pages 699 730 S 005-5718(07)0045-5 Article electronically published on September 6, 007 A DSCONTNUOUS GALERKN FNTE ELEMENT METHOD FOR TME DEPENDENT

More information

ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE

ASYMPTOTICALLY 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 information

ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM

ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM Internat. J. Math. & Math. Sci. Vol. 22, No. 3 (999 587 595 S 6-72 9922587-2 Electronic Publishing House ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR

More information

A stabilized finite element method for the convection dominated diffusion optimal control problem

A stabilized finite element method for the convection dominated diffusion optimal control problem Applicable Analysis An International Journal ISSN: 0003-6811 (Print) 1563-504X (Online) Journal homepage: http://www.tandfonline.com/loi/gapa20 A stabilized finite element method for the convection dominated

More information

Fourier analysis for discontinuous Galerkin and related methods. Abstract

Fourier analysis for discontinuous Galerkin and related methods. Abstract Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the

More information

Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations

Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen

More information

Partially Penalized Immersed Finite Element Methods for Parabolic Interface Problems

Partially Penalized Immersed Finite Element Methods for Parabolic Interface Problems Partially Penalized Immersed Finite Element Methods for Parabolic Interface Problems Tao Lin, Qing Yang and Xu Zhang Abstract We present partially penalized immersed finite element methods for solving

More information

An 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 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 information

Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations

Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen

More information

A posteriori error estimates applied to flow in a channel with corners

A 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 information

Convergence analysis of a finite volume method for the Stokes system using non-conforming arguments

Convergence analysis of a finite volume method for the Stokes system using non-conforming arguments IMA Journal of Numerical Analysis (2005) 25, 523 548 doi:10.1093/imanum/dri007 Advance Access publication on February 7, 2005 Convergence analysis of a finite volume method for the Stokes system using

More information

Abstract. 1. Introduction

Abstract. 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 information

On angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic

On angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen

More information

DISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS

DISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS DISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS CHI-WANG SHU Abstract. In these lectures we give a general survey on discontinuous Galerkin methods for solving

More information

Finite volume method for nonlinear transmission problems

Finite volume method for nonlinear transmission problems Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings

More information

Electronic Transactions on Numerical Analysis Volume 32, 2008

Electronic Transactions on Numerical Analysis Volume 32, 2008 Electronic Transactions on Numerical Analysis Volume 32, 2008 Contents 1 On the role of boundary conditions for CIP stabilization of higher order finite elements. Friedhelm Schieweck. We investigate the

More information

A STABILIZED NONCONFORMING QUADRILATERAL FINITE ELEMENT METHOD FOR THE GENERALIZED STOKES EQUATIONS

A STABILIZED NONCONFORMING QUADRILATERAL FINITE ELEMENT METHOD FOR THE GENERALIZED STOKES EQUATIONS INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number 2, Pages 449 457 c 2012 Institute for Scientific Computing and Information A STABILIZED NONCONFORMING QUADRILATERAL FINITE ELEMENT

More information

Generalized Newtonian Fluid Flow through a Porous Medium

Generalized 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 information