SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
|
|
- Tamsyn Riley
- 5 years ago
- Views:
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
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 informationANALYSIS 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 informationAN 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 informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More 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 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 informationNUMERICAL 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 informationGraded 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 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 informationNUMERICAL 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 informationA 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 informationNUMERICAL 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 informationA 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 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 informationMixed 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 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 informationTWO 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 informationOn 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 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 informationDGFEM. 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 informationCONVERGENCE 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 informationFourier 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 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 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 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 informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More 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 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 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 informationSuperconvergence and H(div) Projection for Discontinuous Galerkin Methods
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 00:1 6 [Version: 2000/07/27 v1.0] Superconvergence and H(div) Projection for Discontinuous Galerkin Methods Peter Bastian
More informationOptimal 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 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 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 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 informationAdaptive 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 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 informationInteraction 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 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 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 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 informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationA 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 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 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 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 informationDiscrete 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 informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationWEAK GALERKIN FINITE ELEMENT 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 informationA 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 informationON 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 informationOn 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 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 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 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 informationGALERKIN 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 informationA 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 informationFinite Element Clifford Algebra: A New Toolkit for Evolution Problems
Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/
More informationA 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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationEquations 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 informationAN 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 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 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 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 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 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 informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
More informationDiscontinuous 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 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 informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationPiecewise 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 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 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 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 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 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 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 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 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 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 informationApplied 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 informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationA 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 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 informationON 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 informationA 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 informationFourier 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 informationNodal 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 informationPartially 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 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 informationResearch 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 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 informationConvergence 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 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 informationOn 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 informationDISCONTINUOUS 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 informationFinite 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 informationElectronic 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 informationA 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 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 information