A CRITICAL COMPARISON OF TWO DISCONTINUOUS GALERKIN METHODS FOR THE NAVIER-STOKES EQUATIONS USING SOLENOIDAL APPROXIMATIONS
|
|
- Phyllis Ray
- 5 years ago
- Views:
Transcription
1 7th Workshop on Numerical Methods in Applied Science and Engineering (NMASE 8) Vall de Núria, 9 a de enero de 28 c LaCàN, A CRITICAL COMPARISON OF TWO DISCONTINUOUS GALERKIN METHODS FOR THE NAVIER-STOKES EQUATIONS USING SOLENOIDAL APPROXIMATIONS A. Montlaur, J. Peraire 2,A.Huerta 3 and S. Fernández-Méndez 3 Laboratori de Càlcul Numèric (LaCàN), web: Universitat Politècnica de Catalunya, E-834 Barcelona, SPAIN : Escola Politècnica Superior de Castelldefels 3: Departament de Matemàtica Aplicada III {adeline.de.montlaur, antonio.huerta, sonia.fernandez}@upc.edu 2: Aerospace Computational Design Laboratory, web: Department of Aeronautics and Astronautics, Massachusetts Institute of Technology, Cambridge, MA 239, USA peraire@mit.edu Key words: Compact Discontinuous Galerkin, solenoidal, Navier-Stokes, incompressible flow, Interior Penalty Method, fractional-step method. Abstract. This paper compares two methods to solve incompressible problems, in particular the Navier-Stokes equations, using a discontinuous polynomial interpolation that is exactly divergence-free in each element. The first method is an Interior Penalty Method, whereas the second method follows the Compact Discontinuous Galerkin [] approach for the diffusive part of the problem. In both cases the Navier-Stokes equations are then solved using a fractional-step method, using an implicit method for the diffusion part and a semiimplicit method for the convection. Numerical examples compare the efficiency and the accuracy of the two proposed methods. INTRODUCTION Since its introduction the Discontinuous Galerkin (DG) method has become popular for its high-order accuracy. DG methods are actually particulary well suited for the construction of high-order accurate space discretization of problems on general unstructured grids, by the mean of high-order polynomial approximations within elements. Another advantage of the DG discretization is its compactness since the coupling is strictly restricted to the elements sharing a face. Recently, the DG method has been successfully used to treat convection-diffusion problems such as the solution of the Navier-Stokes equations. For example Cockburn et al. developed an accurate and stable scheme, the Local Discontinuous Galerkin (LDG) method for the Stokes, Oseen and Navier-Stokes equations [2, 3]. One drawback of this scheme is
2 the loss of the compactness of the discretization form, due to the fact that the degrees of freedom in one element are connected not only to those in the neighboring elements but also to those in the elements neighboring the immediate neighbors. In [] a Compact Discontinuous Galerkin (CDG) method is applied to elliptic problems. The CDG method is very similar to the LDG method but it eliminates the coupled degrees of freedom between the non-neighboring elements. Bassi and Rebay also proposed methods to solve the Navier-Stokes equations using high order DG methods [4] but especially apply it to compressible problems. To solve incompressible problems, an interesting possibility is to use explicit divergence-free bases. In [5] a DG method with solenoidal approximation has been applied to the Stokes equations and an Interior Penalty Method (IPM) employed to construct the weak form. The aim of this paper is to formulate the Compact Discontinuous Galerkin method and the Interior Penalty Method for the incompressible Navier-Stokes equations using solenoidal approximations as seen in [5] and then to present a critical comparison between them. 2 TWO FORMULATIONS OF THE NAVIER-STOKES PROBLEM Let R nsd be an open, bounded domain with piecewise linear boundary andn sd the number of spatial dimensions. Suppose that is partitioned in n el disjoint subdomains K, which for example correspond to different materials, with also piecewise linear boundaries K that define an internal interphase Γ. The strong form for the homogeneous steady Navier-Stokes problem can be written as: 2ν s u + p + (u u) =f in (a) u = in (b) u = u D on Γ D (c) pn +2ν(n s )u = t on Γ N (d) where f L 2 () is a source term, u the flux velocity, p its pressure, ν the kinematic viscosity and s = 2 ( + T ). 2. The Compact Discontinuous Galerkin method In the CDG approach [], the problem is first transformed into a first order in derivative problem, introducing the gradient of the velocity as an additional variable. This additional variable is then eliminated at an elemental level, using the so-called lifting operators. Introducing the velocity gradient σ =2ν s u, the Navier-Stokes equations can be 2
3 rewritten as the following system of first order in derivative equations σ =2ν s u in, (2a) σ + p + (u u) =f in, (2b) u = in, (2c) u = u D on Γ D, (2d) n σ pn = t on Γ N. (2e) 2.. DG formulation of the Navier-Stokes equations Multiplying the equations (2) by smooth test functions τ, v and q respectively and integrating by parts over an arbitrary subset K with outward normal unit vector n the weak problem is σ : τ dx = 2ν u τdx +2ν u nds, (3a) K K K [ σ : s ] v p v dx [σ : (v n) pv n] ds K K (u u) vdx + u n u vds = f vdx (3b) K K K u q + u nqds = (3c) Note that the above equations are well defined for functions (σ, u,p)and(τ, v,q)in Σ V Q where the spaces K Σ = {τ [L 2 ()] n2 sd τ K [H (K)] n2 sd K T h }, (4a) V = {v [L 2 ()] nsd v K [H (K)] nsd K T h }, (4b) Q = {q L 2 () qdx =,q K H (K) K T h }, (4c) are associated with the triangulation T h = {K} of. In addition the finite element subspaces Σ h Σ, V h V and Q h Q are introduced as Σ h = {τ [L 2 ()] n2 sd τ K [P k (K)] n2 sd K T h }, (5a) V h = {v [L 2 ()] nsd v K [P k (K)] nsd K T h }, (5b) Q h = {q L 2 () qdx =,q K P k (K) K T h }, (5c) where P k (K) is the space of polynomial functions of degree at most k onk. To write expressions applicable over the whole domain, the jump and average operators need to be defined. Two adjacent elements K + and K of the triangulation T h are 3 K
4 considered and e = K + K is their common edge. Assume that n ± denote the unit vectors normal to K ± at any point on the edge e. The average operator is defined as { κ + + +κ on Γ, { } = on. and the jump operator as { p + n + + p n on Γ p n = p n on { n + v + + n v on Γ n v = n v on { n+ v + + n v on Γ n v = n v on { n+ σ + + n σ on Γ n σ = n σ on for scalars, or for vectors, for second order tensors. Using these definitions, the DG formulation of the weak form (3) is: find σ h Σ h, u h V h, p h Q h such that for all K T h σ h : τ = 2ν u h h τ +2ν û σ h τ n ds, (6a) Γ [σ h : s hv p h h v]dx [ˆσ h : v n ˆp h v n ]ds Γ (u h u h ) h vdx + ˆf(u h, u h ) n vds = f vdx (6b) Γ u h h q + û p h nq ds = (6c) Γ where û σ h, ˆσ h,ˆp h, ˆf(u h, u h )andû p h are numerical fluxes to be defined. Note that if we use the integration by parts formula v τdx = τ s vdx ( n v :{τ }+{v} n τ )ds Γ v τ nds (7) valid for all τ [H (T h )] n2 sd and u [H (T h )] nsd, equation (6a) can be written as σ h : τ dx =2ν τ s hu h dx 2ν ( n u h :{τ } {û σ h u h } n τ )ds Γ +2ν (û σ h u h ) τ nds (8) 4
5 Let us now describe the different numerical fluxes of equation (6). The convective numerical flux. The convective flux ˆf(u h, u h ) is defined in a similar way as Lomtev and Karniadakis do in [6]. A standard Roe-splitting flux is used ˆf(u, u) n = [ (F x (u e ), F y (u e )) + ( F x (u i ), F y (u i ) )] n 2 2 R Λ L(ue u i ) where u i and u e are respectively the velocities in the interior and in the exterior of the element under consideration, ( ) ( ) u 2 uv where F x (u) = and F y (u) =, uv v 2 and A, jacobian matrix of F, is written in terms of the left and right eigenvectors: A = RΛL, where Λ contains the corresponding eigenvalues in the diagonal. The diffusive numerical fluxes. Following the definition of the diffusive fluxes by Cockburn et al. for their LDG method [2], if a face e lies inside the domain ˆσ h = {σ h } C n u h + C 2 n σ h, û σ h = {u h } C 2 n u h on Γ, and if e lies on the boundary, ˆσ h = σ h C (u h u D ) n, û σ h = u D on Γ D. C is a positive constant and C 2 is a vector that has to be determined for each interior edge of the domain according to C 2 = 2 (SK K + n+ + S K+ K n ) where S K K {, } is a switch defined for each element edge. The switches always satisfy + that S K K + + SK+ K =. There is several possible choices of the switches, see for example [2]. One possibility is the natural switch, which takes into account the element numbering to set S K± K.Another alternative is to use a consistent switch that satisfies < S K± K < n sd +. (9) e K Calculating C 2 depending on the numbering of the nodes for each element K is the option chosen here and is a consistent switch. The numerical fluxes related to the incompressibility constraint. If a face e lies inside the domain û p h = {u h}, ˆp h = {p h } on Γ, 5
6 whereas on the boundary û p h = u D, ˆp h = p h on Γ D. Finally in order to account for the natural condition on Γ N, we must choose ˆσ h and ˆp h such that n ˆσ h ˆp h n = t on Γ N and set û p h = u and ûσ h = u on Γ N 2..2 CDG formulation of the steady Navier-Stokes equations Using the expression of the fluxes previously defined, equation (8) can now be written σ h : τ dx =2ν τ s hu h dx 2ν ( n u h :{τ } + C 2 n u n τ )ds Γ +2ν (u D u h ) τ nds () Γ D To obtain an expression for σ h as a function of u h, lifting operators are introduced following the method introduced by Arnold et al. [7]. For all e Γ,r e :[L 2 (e)] n2 sd Σ h,l e :[L 2 (e)] nsd Σ h and for each e Γ D,rD e :[L 2(e)] nsd Σ h,aredefinedas r e (v) : τ dx = v :{τ }ds, τ Σ h, e l e (v) : τ dx = v n τ ds, τ Σ h, () e rd(v) e : τ dx = v n τ ds, τ Σ h. e Note that here the lifting operators are defined edgewise, which is where CDG differs from LDG, see []. The lifting operators using by the LDG methods are defined as: r :[L 2 (Γ)] n2 sd Σh,l :[L 2 (Γ)] nsd Σ h and for each e Γ D,r D :[L 2 (Γ D )] nsd Σ h,are defined as r(v) : τ dx = v :{τ }ds, τ Σ h, Γ l(v) : τ dx = v n τ ds, τ Σ h, (2) Γ r D (v) : τ dx = v n τ ds, τ Σ h. Γ D 6
7 From these definitions the correspondence between the CDG and LDG lifting operators is clear: for all v [L 2 (Γ)] n2 sd r(v) = e Γ r e (v), l(v) = e Γ l e (v), r D (v) = e Γ r e D(v). (3) Using this notation, the equation (6a) defining σ h in terms of u h can be rewritten as σ h =2ν s hu h σ e (4) where σ e =2ν(r e ( n u h )+l e (C 2 n u h ) rd(u e D )+rd(u e h )). (5) Following the rationale detailed in [] for an elliptic problem, eliminating the additional variable σ h, the CDG scheme for the steady Navier-Stokes equations is: find u h V h and p h V h such that A CDG h (u h, v)+o h (u h, u h, v)+b h (v,p h )=Fh CDG (v), v V h (6) where A CDG h (u, v) = 2ν s u : s vdx 2ν( n u :{ s v} + { s u} : n v )ds Γ Γ D 2ν(C 2 n u : n s v + C 2 n s u : n v )ds Γ + 2ν(r e ( n u )+l e (C 2 n u )+rd(u)) e :(r e ( n v )+l e (C 2 n v )+rd(v))dx e e Γ + C n u n v ds Γ Γ D O h (u, u; v) = (u u) h vdx + ˆf(u, u) n vds Γ Γ D B h (v,p)= p h vdx + {p} v n ds Γ Γ D Fh CDG (v) = f vdx 2νn u D : s vds Γ D 2ν(n u D ) :(r e ( n v )+l e (C 2 n v ))dx e Γ 2νn v : rd(u e D )dx + C (n u D ) (n v)ds + v tds e Γ Γ D Γ N Using the following identity for any v [H ( )] nsd and any q L 2 () q vdx + {q} n v ds = v qdx qn {v}, Γ 7 Γ
8 the incompressibility equation (6c) can be rewritten as B h (u,q)=g h (q), (7) where G h (q) = qu D nds. (8) Γ D Note that the compactness of the CDG scheme compared to the LDG comes from the product of the lifting operators (in A CDG h (u, v)). Indeed in the CDG scheme the lifting operators are associated to individual edges and there are no connectivities between nonneighbor elements, whereas in the LDG scheme, the product of the lifting operator implies connectivities between non-neighbor elements. 2.2 The Interior Penalty Method Following the approach presented in [5] for the Stokes equations, an Interior Penalty Method formulation is proposed for the Navier-Stokes equations. The convective numerical flux and the numerical flux related to the incompressibility constraint are the same as the ones used in the CDG approach. Concerning the diffusive numerical fluxes, a penalty parameter C is introduced to ensure the coercivity of the bilinear form A IPM h (u, v) = 2ν s u : s vdx 2ν( n u :{ s v} + { s u} : n v )ds Γ Γ D + C n u n v ds Γ Γ D Using the same discrete spaces as the ones previously defined, the IPM scheme for the steady Navier-Stokes equations is: find u h V h and p h P h such that A IPM h (u h, v)+o h (u h, u h, v)+b h (v,p h )=Fh IPM (v), v V h (9) where O h (u, u, v) andb h (v,p) are defined as for the CDG method and Fh IPM (v) = f vdx 2νn u D : s vds+ C (n u D ) (n v)ds+ v tds Γ D Γ D Γ N See that this formulation is compact in the sense that the degrees of freedom of one element are only connected to those of the immediate neighbors. 3 SOLENOIDAL INTERPOLATION SPACES Following approaches such as the ones from Crouzeix and Raviart [8] or Griffiths [9], who were the first ones to introduce solenoidal spaces, and the approach proposed in 8
9 [, ], which uses for the interpolation of the velocity piecewise polynomials (exactly) divergence-free, a similar rationale is pursued following the one presented in [5]. Knowing that any function in [H (K)] nsd, can be written as the sum of a solenoidal part and an irrotational one, the functional space for the velocity can be split in the direct sum: [H ()] nsd = S I where S := { } v [H ()] nsd v K = for i =,...,n el, I { } v [H ()] nsd v K = for i =,...,n el. Note also that the solution of the original problem () or (6) belongs to S. The finite counterparts of S and I are S h = { } v [H ()] nsd v K [P k (K)] nsd, v i = for i =,...,n el, I h { } v [H ()] nsd v K [P k (K)] nsd, v i = for i =,...,n el, such that I h I. Note that the following relations and inclusions are verified: V h = S h I h, V h [H ()] nsd, Q h L 2 () and S h S. Under these circumstances, problem (6) can be split in two uncoupled problems. The first one solves for the divergence-free velocities and hybrid pressures. Following [2] the space of hybrid pressures (pressures along the edges in 2D or faces in 3D) is simply: P := { p p :Γ Γ D R and p = n v for some v S}. (2) Note that the finite dimensional subspace associated with the hybrid pressures, P h P, can be defined directly from (2) restricting velocities to S h. The approximations obtained in this case are u h S h and p h P h solution of { Ah (u h, v)+bh S (v, p h )+O h (u h, u h ; v) =F h (v) v S h, (2a) Bh S (u h, q) =G h (q) q P h, where Bh S (v, p) = p h v n ds Γ Γ D and A h (u h, v) =A CDG h (u h, v) ora h (u h, v) =A IPM h (u h, v), and F h (v) =Fh CDG (v) or F h (v) =Fh IPM (v), depending on which method is used. The second problem, which requires the solution of the previous one, evaluates the interior pressures: find p h Q h B I h(v,p h )=F h (v) A h (u h, v) B S h (v, p h ) O h (u h, u h ; v) v I h, (2b) where Bh(v,p)= I p h vdx. (22) It is important to note that equation (2b) can be solved element by element and pressure is the only unknown. 9
10 4 SOLUTION OF THE NAVIER-STOKES EQUATIONS BY A FRAC- TIONAL STEP METHOD To solve the non-linear Navier-Stokes equations an iterative method in time is used. Starting with the weak form obtained for the steady case (2a), the weak form associated to the unsteady problem is now: given u h S h find u h S h and p h P h solution of ( uh ) t, v + A h (u h, v)+bh S (v, p h )+O h (u h, u h ; v) =F h (v) v S h, (23) Bh S (u h, q) =G h (q) q P h Once again the expressions of A h (u h, v) andf h (v) depend on the choice of the method, CDG or IPM. Let us denote by u n h and pn h the velocity and pressure at time tn,wheret n = t n + t. Following the method developed in [3], the fields u n+ h and p n+ h are calculated from u n h but here an implicit Euler method with a pressure increment is used to derive u h for the diffusive part (Stokes problem). Indeed, it has been observed, by Goda in [4], that adding an old value of the pressure gradient in the first sub-step and then accordingly correcting the velocity in the second sub-step increases the accuracy of the scheme and improves the convergence properties. The idea was then made popular by Van Kan [5] who proposed a second-order incremental pressure-correction scheme. Using the Backward Difference Formula of second order to approximate the time derivative, and a semi-implicit scheme for the nonlinear convective term of the Navier-Stokes equation, the incremental pressurecorrection reads as follows. The first step is: from u n h, un h and p n h,findu h [H ()] nsd such that ( 3u h 4u n h + un h, v 2 t ) + A h (u h, v)+o h (2u n h u n h, u n h; v) =Fh(v) Bh S (v, p n h) v S h, (24a) then the second step is: from u h and p n h,findun+ h ( 3u n+ h S h and p n h P h such that ), v + Bh S (v, p n+ h p n h)= v S h, 3u h 2 t Bh S (u n+ h, q) =G h (q) q P h. (24b) Denoting by u and p the vectors of nodal values of the velocity u and hybrid pressure p, the following matrix form is obtained: the first sub-step is u =(3M +2 t(a + O(2u n u n )) (M(4u n u n )+2 t(f C T p n )) and the second substep p n+ = p n (CM C T ) Cu G, t u n+ = u 2 3 tm C T (p n+ p n ). (25a) (25b)
11 where M is the mass matrix, A the diffusion matrix, O the convection matrix, C T the edge pressure matrix and F and G vectors taking into account the force term and the boundary conditions. After the convergence of the iterative process, the interior pressure p is retrieved by a direct method solving equation (2b). 5 NUMERICAL EXAMPLES 5. Driven cavity example A standard benchmark test for the Navier-Stokes equations is first considered. A plane flow of an isothermal fluid in a lid-driven cavity is modelled in a 2D square domain = ], [ ], [, with zero body force and one moving wall. A continuous velocity (u =(, ) T for. x.9 u =(x, ) T for x. and.9 x ) is imposed on the exterior upper boundary {y =}, and a zero velocity u =(, ) T is enforced on the three other sides. Calculations are made for Re =,, 4. Figure shows the results, which fit to the expected solution. As the Reynolds number is increased, the position of the main vortex moves towards the center of the cavity, both CDG method and IPM giving very similar results. In order to further compare them, the velocity profiles at the vertical centerline are shown in figure 2 for different values for Re. First, it can be noticed that as the Reynolds number increases, the boundary layers are more obvious and the variations in the velocity are sharper. Second, results for CDG method and IPM are almost identical. To compare more precisely the results of both methods an analytical example is taken in the next section. 5.2 General results on a Stokes example A Stokes example with analytical solution is now considered to compare the two proposed methodologies and in particular to study the properties of the diffusion matrix and the error of both methods. The Stokes equations are solved in a 2D square domain =], [ ], [ with Dirichlet boundary conditions on three edges, and a Neumann boundary condition on the fourth edge {y =}. A body force 2( 2y)x ( +2y)x 3 + 2( 4y +6y 2 4y 3 +)x 2 +( (y 3y 2 +2y 3 ))x + 4y +2y 2 8y 3 f = 8( 6y +6y 2 )x 3 + 2( +6y 6y 2 )x 2 +(4 + 48(y 2 y 3 ) + 24(y 4 y))x 2y 2 +24y 3 2y 4 is imposed in order to have the polynomial exact solution ( ) x 2 ( x) 2 (2y 6y 2 +4y 3 ) u =, y 2 ( y) 2 (2x 6x 2 +4x 3 ) p = x( x).
12 (a) Re = (b) Re = (c) Re = 4 Figura : Velocity streamlines for the CDG method (left) and the IPM (right), for different Reynolds numbers Re =,, 4, k = 2, 45 elements, C = (CDG), C =(IPM). 2
13 .5 CDG IPM.5 CDG IPM (a) Re = (b) Re =.5 CDG IPM (c) Re = 4 Figura 2: Velocity profiles at the vertical centerline for the CDG and IPM methods, for different values of Re 3
14 5 7 Condition number 4 3 CDG IPM Condition number 6 CDG IPM 2-5 C C 5 (a) k =2 (b) k =4 Figura 3: Condition number of the diffusion matrix function of the stabilization parameter C for the CDG and IPM methods, for a quadratic approximation of the velocity (left) and a fourth order of approximation of the velocity (right) First the influence of the stabilization parameter C on the conditioning of the diffusion matrix is studied. Figure 3 shows the values of the condition number of the diffusion matrix for two different orders of interpolation and for the CDG and IPM methods. Note that unlike the Interior Penalty Method, the CDG methods (and it would also be the case for the LDG method) allow to choose a value of C as small as wanted for the interior edges, see []. It gives thus more freedom to chose this parameter and the value of the condition number is rather constant for values of C less or equal than. On the other hand, the value of C for the IPM has to be big enough to ensure the coercivity of the bilinear form, for example C for a quadratic velocity and C 4 for a fourth order interpolation of the velocity. Then the IPM gives slightly better results than the CDG method. Then, increasing the value of C increases dramatically the condition number, equally for both methods. Since there is a larger range of valid values for the stabilization parameter C for the CDG method let us study how its choice affects the accuracy of the solution of the Stokes equation. Figure 4 shows the results for the velocity and pressure errors for three different values of C =,,. It can be seen that though a value of or for C gives very similar results, the value of C = give worse results, especially for a coarse mesh in the velocity case and for any size of the mesh in the pressure case. Note that the accuracy is only weakly dependent on the value of C since optimal convergence rate are obtained in any case. These results lead us to choose a small value for C for the CDG method C = for example, which will be optimal from both the condition numbers and the error estimates. Note that if C = is chosen, a non-zero C has to be taken on the boundaries in order to treat properly the boundary conditions. Concerning the IPM, given that any value 4
15 -2 velocity error L C = C = C = -2 h - pressure error L 2 2. C = C = C = -5-2 h - Figura 4: Comparison of the errors obtained with the CDG method for different values of C,fora quadratic approximation of the velocity (left) and a linear interpolation of the pressure (right). -2 velocity error L CDG IPM -2 h - pressure error L CDG IPM -2 h - Figura 5: Comparison of the errors obtained with the CDG and the IPM methods, for a quadratic approximation of the velocity (left) and a linear interpolation of the pressure (right). of C bigger than the minimum value necessary to ensure the coercivity of the bilinear form will deteriorate the conditioning of the diffusion matrix, it is preferable to chose the smallest value possible of C. Both CDG and IPM are now compared from the accuracy point of view, for a quadratic approximation of the velocity and a linear interpolation of the pressure and choosing C = for CDG and C = for IPM. As it can be seen in figure 5 the results are very similar. Depending on the mesh size, one or the other method gives a better result but in any case both methods reach optimal convergence rates for velocity and pressure errors. 6 CONCLUSIONS A critical comparison of the Compact Discontinuous Galerkin method and the Interior Penalty Method is presented here. The two methods are derived using different rationales. The CDG method first transforms the Navier-Stokes problem into a first order in derivative problem and then eliminate the additional variable using the so-called lifting 5
16 operators. The IPM rationale is more straight-forward and deals directly with the second order in derivative problem. But both problems present the major advantage -relative to the LDG method for example- that they are compact formulations: the degrees of freedom of one element are only connected to those of the immediate neighbors. The CDG presents the other advantage that the stabilization parameter can take any value, whereas the penalty parameter in the IPM has to be taken big enough to ensure the coercivity of the bilinear form. Apart from this difference in the choice of the characteristic parameters both methods present very similar results, from the conditioning and the accuracy points of view, both methods reach optimal convergence rates for velocity and pressure errors. REFERENCIAS [] J. Peraire and P.O. Persson, A Compact Discontinuous Galerkin (CDG) method for elliptic problems, In review (26) [2] B. Cockburn, G. Kanschat, I. Perugia and D. Schotzau, Superconvergence of the Local Discontinuous Galerkin method for elliptic problems on cartesian grids, SIAM J. Numer. Anal., 39,, (2) [3] B. Cockburn, G. Kanschat and D. Schotzau, The local discontinuous Galerkin method for linearized incompressible fluid flow: a review, Comp. Fluids, 34, 49-56,(25) [4] F. Bassi and S. Rebay, A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier-Stokes equations,j. Comput. Phys., 3, , (997) [5] A. Montlaur, S. Fernandez-Mendez and A. Huerta, A discontinuous Galerkin method with divergence-free interpolation for the incompressible Stokes equations, Int. J. Numer. Meth. Fluids, inpress [6] I. Lomtev and G. Karniadakis, A Discontinuous Galerkin method for the Navier- Stokes equations, Int. J. Numer. Methods Fluids, 29, , (999) [7] 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, , (22) [8] M. Crouzeix and P.-A. Raviart, Conforming and nonconforming finite element methods for solving the stationary Stokes equations. I, Rev. Française Automat. Informat. Recherche Opérationnelle Sér. Rouge, vol. 7, no. R-3, pp ,
17 [9] D. F. Griffiths, An approximately divergence-free 9-node velocity element (with variations) for incompressible flows,int. J. Numer. Methods Fluids, vol., no. 4, pp , 98. [] G.A. Baker, W.N. Jureidini and O.A. Karakashian, Piecewise solenoidal vector fields and the Stokes problem, SIAM J. Numer. Anal., 27, 6, , (99) [] O.A. Karakashian and W.N. Jureidini, A nonconforming finite element method for the stationary Navier-Stokes equations, SIAM J. Numer. Anal., 35,, 93-2, (998) [2] B. Cockburn and J. Gopalakrishnan, Incompressible finite elements via hybridization. Part I the Stokes system in two space dimensions, SIAMJ.Numer.Anal., 43, 4, (25) [3] J. Donea, S. Giuliani, H. Laval and L. Quartapelle, Finite element solution of the unsteady Navier-Stokes equations by a fractional step method, Comput. Meth. Appl. Mech. Eng.,3, 53-73, (98) [4] K. Goda, A multistep technique with implicit difference schemes for calculating twoor three-dimensional cavity flows, J. Comput. Phys., 3, 76-95, (979) [5] J. van Kan, A second-order accurate pressure-correction scheme for viscous incompressible flow, SIAM J. Sci. Stat. Comput.,7,3, 87-89, (986) Con el apoyo de Universitat Politècnica de Catalunya y E.T.S. d Enginyers de Camins, Canals i Ports de Barcelona 7
Discontinuous Galerkin methods for the Stokes equations using divergence-free approximations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 27; : 9 [Version: 22/9/8 v.] Discontinuous Galerkin methods for the Stokes equations using divergence-free approximations
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 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 informationTHE COMPACT DISCONTINUOUS GALERKIN (CDG) METHOD FOR ELLIPTIC PROBLEMS
SIAM J. SCI. COMPUT. Vol. 30, No. 4, pp. 1806 1824 c 2008 Society for Industrial and Applied Mathematics THE COMPACT DISCONTINUOUS GALERKIN (CDG) METHOD FOR ELLIPTIC PROBLEMS J. PERAIRE AND P.-O. PERSSON
More informationPseudo-divergence-free element free Galerkin method for incompressible fluid flow
International Workshop on MeshFree Methods 23 Pseudo-divergence-free element free Galerkin method for incompressible fluid flow Y. Vidal () and A. Huerta (2) Abstract: Incompressible modelling in finite
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 information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More 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 informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
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 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 informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
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 informationHybrid Discontinuous Galerkin methods for incompressible flow problems
Hybrid Discontinuous Galerkin methods incompressible flow problems Christoph Lehrenfeld, Joachim Schöberl MathCCES Computational Mathematics in Engineering Workshop Linz, May 31 - June 1, 2010 Contents
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
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 informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More 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 informationInterior penalty tensor-product preconditioners for high-order discontinuous Galerkin discretizations
Interior penalty tensor-product preconditioners for high-order discontinuous Galerkin discretizations Will Pazner Brown University, 8 George St., Providence, RI, 9, U.S.A. Per-Olof Persson University of
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 DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
More 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 informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
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 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 informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
More 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 informationextended Hybridizable Discontinuous Galerkin (X-HDG) for Void and Bimaterial Problems
extended Hybridizable Discontinuous Galerkin (X-HDG) for Void and Bimaterial Problems Ceren Gürkan, Esther Sala-Lardies, Martin Kronbichler, and Sonia Fernández-Méndez Abstract A strategy for the Hybridizable
More informationTwo-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5
Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Jingwei Zhu May 14, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationRANS Solutions Using High Order Discontinuous Galerkin Methods
RANS Solutions Using High Order Discontinuous Galerkin Methods Ngoc Cuong Nguyen, Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 2139, U.S.A. We present a practical
More informationHIGH-ORDER DISCONTINUOUS GALERKIN METHODS FOR INCOMPRESSIBLE FLOWS
HIGH-ORDER DISCONTINUOUS GALERKIN METHODS FOR INCOMPRESSIBLE FLOWS Adeline de Montlaur Doctoral Thesis Advisors: Antonio Huerta and Sonia Fernández-Méndez Barcelona, September 2009 Escola Politècnica Superior
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
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 informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationMultigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations. Todd A.
Multigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations by Todd A. Oliver B.S., Massachusetts Institute of Technology (22) Submitted to the
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationSub-Cell Shock Capturing for Discontinuous Galerkin Methods
Sub-Cell Shock Capturing for Discontinuous Galerkin Methods Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 39, U.S.A. A shock capturing strategy for higher order
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More 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 informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More informationOn Dual-Weighted Residual Error Estimates for p-dependent Discretizations
On Dual-Weighted Residual Error Estimates for p-dependent Discretizations Aerospace Computational Design Laboratory Report TR-11-1 Masayuki Yano and David L. Darmofal September 21, 2011 Abstract This report
More informationA SIMPLE TUTORIAL ON DISCONTINUOUS GALERKIN METHODS. Jennifer Proft CERMICS, ENPC. J. Proft CERMICS, ENPC
A SIMPLE TUTORIAL ON DISCONTINUOUS GALERKIN METHODS Jennifer Proft Outline Definitions A simple example Issues Historical development elliptic equations hyperbolic equations Discontinuous vs. continuous
More informationA 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 informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
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 informationUnified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems.
Unified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems. Stefan Girke WWU Münster Institut für Numerische und Angewandte Mathematik 10th
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationA super convergent hybridisable discontinuous Galerkin method for linear elasticity
A super convergent hybridisable discontinuous Galerkin method for linear elasticity arxiv:82.223v [math.na] 6 Feb 28 Ruben Sevilla Zienkiewicz Centre for Computational Engineering, College of Engineering,
More informationHybridizable Discontinuous Galerkin with degree adaptivity for the incompressible Navier-Stokes equations
Hybridizable Discontinuous Galerkin with degree adaptivity for the incompressible Navier-Stokes equations Giorgio Giorgiani, Sonia Fernández-Méndez, Antonio Huerta To cite this version: Giorgio Giorgiani,
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 informationEnergy Norm A Posteriori Error Estimation for Mixed Discontinuous Galerkin Approximations of the Stokes Problem
idgenössische Technische Hochschule Zürich cole polytechnique fédérale de Zurich Politecnico federale di Zurigo Swiss Federal Institute of Technology Zurich nergy Norm A Posteriori rror stimation for Mixed
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationSecond-order finite-element projection method for 3D flows
Second-order finite-element projection method for 3D flows J.-L. Guermond I and L. Quartapelle 2 1 LIMSI-CNRS, BP 133, 91403 Orsay, France, (guermond@limsi.fr) 2 Dipartimento di Fisica del Politecnico
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 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 informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationLocking in the incompressible limit: pseudo-divergence-free element free Galerkin
Locking in the incompressible limit: pseudodivergencefree element free Galerkin Yolanda Vidal Pierre Villon and Antonio Huerta Departament de Matemàtica Aplicada III, ES de Ingenieros de Caminos, Canales
More informationMultigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations. Todd A.
Multigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations by Todd A. Oliver B.S., Massachusetts Institute of Technology (22) Submitted to the
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
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 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 informationAn arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem
An arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem Yifan Wang University of Houston, Department of Mathematics
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationNEWTON-GMRES PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE NAVIER-STOKES EQUATIONS
NEWTON-GMRES PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE NAVIER-STOKES EQUATIONS P.-O. PERSSON AND J. PERAIRE Abstract. We study preconditioners for the iterative solution of the
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
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 informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
More informationDiscontinuous Galerkin Solution of the Navier-Stokes Equations on Deformable Domains
Discontinuous Galerkin Solution of the Navier-Stokes Equations on Deformable Domains P.-O. Persson a J. Bonet b and J. Peraire c, a Department of Mathematics, University of California, Berkeley, Berkeley,
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
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,a Paul F. Fischer b C. Ross Ethier a a University of Toronto, Department of Mechanical
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 informationChapter 1. Introduction and Background. 1.1 Introduction
Chapter 1 Introduction and Background 1.1 Introduction Over the past several years the numerical approximation of partial differential equations (PDEs) has made important progress because of the rapid
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationA Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization. Laslo Tibor Diosady
A Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization by Laslo Tibor Diosady B.A.Sc., University of Toronto (2005) Submitted to
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 informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationMultigrid Algorithms for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations
Multigrid Algorithms for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations Khosro Shahbazi,a Dimitri J. Mavriplis b Nicholas K. Burgess b a Division of Applied
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 informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
HONOM 2011 in Trento DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology DLR Braunschweig 11. April 2011 1 / 35 Research
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
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 informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
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 informationNumAn2014 Conference Proceedings
OpenAccess Proceedings of the 6th International Conference on Numerical Analysis, pp 198-03 Contents lists available at AMCL s Digital Library. NumAn014 Conference Proceedings Digital Library Triton :
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
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 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 informationA STUDY OF MULTIGRID SMOOTHERS USED IN COMPRESSIBLE CFD BASED ON THE CONVECTION DIFFUSION EQUATION
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationAn Equal-Order DG Method for the Incompressible Navier-Stokes Equations
J Sci Comput 009) 40: 188 10 DOI 10.1007/s10915-008-961-1 An qual-order DG Method for the Incompressible Navier-Stokes quations Bernardo Cockburn Guido anschat Dominik Schötzau Received: 4 September 008
More information