A New Hyperbolic Navier-Stokes Formulation and Efficient Reconstructed Discontinuous Galerkin Discretizations
|
|
- Rodney Gordon
- 5 years ago
- Views:
Transcription
1 AIAA AVIATION Forum June Atlanta Georgia 28 Fluid Dynamics Conference.254/ A New Hyperbolic Navier-Stokes Formulation and Efficient Reconstructed Discontinuous Galerkin Discretizations Lingquan Li * Jialin Lou Hong Luo North Carolina State University Raleigh NC USA Hiroaki Nishikawa National Institute of Aerospace Hampton VA USA Downloaded by Hiroaki Nishikawa on July DOI:.254/ In this paper we propose a new hyperbolic Navier-Stokes formulation and develop efficient highorder reconstructed discontinuous Galerkin schemes for solving the compressible Navier-Stokes equations. The new formulation uses the gradients of the primitive variables as additional variables to form a hyperbolic viscous system. This formulation enables simple and economical solution reconstruction and paves the way to the development of a family of truely efficient high-order hyperbolic Navier-Stokes schemes: more efficient than discontinuous Galerkin schemes for comparable orders of accuracy despite the increase in the number of variables/equations due to the hyperbolic formulation at the partial-differential-equation level. With 2 degrees of freedom per cell the method achieves on irregular grids third-order error convergence in the solution and gradients and fourthorder in the inviscid terms. The eigen-structure of the hyperbolic viscous system is analyzed and a suitable dissipation matrix is derived. A class of reconstructed discontinuous Galerkin discretizations are constructed for the hyperbolic Navier-Stokes sytem and solved by an implicit solver based on the first-order backward Euler pseudo-time integration scheme with exact or approximate Jacobians. The methods are implemented in a three-dimensional unstructured code verified by the method of manufactured solutions on hexahedral regular irregular and highly distorted prismatic grids and tested for a Couette flow problem. I. Introduction Today Computational Fluid Dynamics (CFD is applied to a variety of applications in science and engineering. Second-order algorithms dominate practical three-dimensional CFD codes but there is a growing need for higherorder accuracy especially in high-fidelity time-dependent turbulent-flow simulations over complex geometries. Beside many challenges e.g. high-order grid generation and adaptation high-performance computing robust and efficient solver developments etc. there exits a classic challenge about how to deal with the viscous terms in Navier-Stokes (NS equations. Although high-order methods have been well developed for hyperbolic systems such as the inviscid terms they have been less studied for the viscous terms. In particular the widely-used and studied finite-volume (FV and Discontinuous Garlerkin (DG methods become less certain and advantageous for elliptic or parabolic problems requiring extra efforts needed to deal with gradients in the viscous fluxes. Commonly used high-order viscous methods are the BR2 [] and DDG [2] methods. A reconstructed discontinuous Galerkin (rdg method using an inter-cell reconstruction is also developed [3]. All of these methods however are substantially more demanding in computational efforts than the classical continuous finite element methods that are more suited for solving elliptic problems. In this work we explore an alternative to these high-order viscous discretization methods: the hyperbolic Navier-Stokes (HNS method in combination with the rdg method to drive the cost of high-order schemes down to the level at which applications to practical problems would be within reach with currently available resources. Ultimately the developed method is expected to provide an efficient and accurate high-order unstructured-grid solver for scalable high-order CFD computations on fully adaptive grids with ever-increasing parallelism. The HNS method is a first-order hyperbolic system reformulation approach which is originally developed for diffusion equations [4] and the advection-diffusion equations [5]. In the HNS method diffusion/viscous terms are reformulated as a hyperbolic system and discretized by methods for hyperbolic systems e.g. upwind methods. Thus methods developed for hyperbolic systems including high-order methods are directly applicable to diffusion and *Graduate Research Assistant Department of Mechanical and Aerospace Engineering AIAA Student Member. lli32@ncsu.edu Graduate Research Assistant Department of Mechanical and Aerospace Engineering AIAA Student Member. Professor Department of Mechanical and Aerospace Engineering Associate Fellow AIAA. Associate Research Fellow Associate Fellow AIAA. of 28 Copyright 28 by authors. Published by the Inc. with permission.
2 Downloaded by Hiroaki Nishikawa on July DOI:.254/ viscous terms. Moreover the method has been demonstrated to bring various other advantages: high-order and noisefree gradients on severely distorted grids upgraded accuracy in the advection/inviscid schemes and convergence acceleration by the elimination of numerical stiffness due to second-order derivatives. These advantages have been demonstrated for various applications including the two-dimensional compressible/incompressible NS equations [6 7 8] three-dimensional compressible NS equations [ ] with a proper handling high-reynolds-number boundary-layer flows [ 5] unsteady problems [6 7] magnetohydrodynamic models [8 9] third-order partialdifferential equations [2] and diffusion with tensor coefficients [2] and with discontinuous coefficients [22]. In terms of the discretization method the hyperbolic method has been applied to the residual-distribution methods [4 5] the finite-volume methods [ ] the high-order active-flux method [6] high-order compact methods [26] and high-order rdg methods [ ]. This paper is a sequel to the high-order hyperbolic rdg method development and presents the extension to the compressible NS equations. For diffusion equations a hyperbolic system is formed typically with the solution gradients introduced as additional variables which are called the gradient variables. In three dimensions it requires three additional gradient variables per solution variable meaning three additional equations are added to each equation in a target system. This increase in variables and equations presents a challenge in applying the hyperbolic method to DG methods. For example if the three-dimensional compressible NS system is formulated as a hyperbolic system with 2 equations instead of 5 [2] even a second-order DG scheme method would result in 8 discrete equations. To minimize the number of discrete equations we have explored an efficient reconstruction of recycling the gradient variables to form high-order polynomials of the solution variables for model equations [ ]. The recycling strategy dramatically reduces the number of discrete equations and leads to efficient hyperbolic DG schemes. For the NS equations it can bring the number of discrete equations down to 2 from 8 without losing second-order accuracy and even with an upgrade of the inviscid accuracy to third-order. For example the HNS formulation in Ref.[2] involves the density gradient as gradient variables. If DG(P method is used the density gradient (ρ x ρ y ρ z and their gradients i.e. second-order derivatives (ρ xx ρ xy ρ xz ρ yx ρ yy ρ yz ρ zx ρ zy ρ zz would be introduced as discrete unknowns. Then these quantities can be used to build a quadratic polynomial for the density and thus only the cell-average needs to be stored for the density. This particular construction is denoted by DG(PP2+DG(P and leads to a third-order inviscid and second-order viscous scheme. However this efficient construction is not straightforward for other conservative or primitive variables because other gradient variables are the velocity gradients scaled by the viscosity and the heat fluxes. The derivatives of the velocity and temperature which are needed to build the polynomials for the velocity and temperature can be quite complicated if obtained from the gradient variables as they will involve derivatives of the viscosity. To make the recycling trivial for all variables we here propose a new HNS formulation HNS2G with the gradients of the primitive variables as additional unknowns. A similar hyperbolic formulation has been employed for nonlinear or variable-coefficient diffusion problems [ ]; this paper presents for the first time its extension to the compressible NS equations. To further reduce the cost of the hyperbolic DG schemes we focus on the rdg method which has been widely used for compressible inviscid and viscous flow problems as well as for turbulent problems [ ]. The rdg method is a general framework that includes the FV and DG methods as special cases: define an underlying high-order polynomial of on each element reconstruct a higher-order polynomial and solve a Riemann problem that arises from the discontinuous representation of solution on each element interface. An rdg scheme with a polynomial of degree n and a reconstructed polynomial of degree n + is denoted by rdg(pnpn+. As discussed in the previous papers [ ] high-order schemes that combine the hyperbolic approach and the rdg method require less number of discrete unknowns and achieve higher-order accuracy than conventional DG schemes. The reconstruction step is a key to accuracy and stability of the rdg schemes and a recently developed high-order reconstruction technique called the variational reconstruction (VR provides excellent performance of the rdg schemes as demonstrated for compressible inviscid flow [4 42]. This VR reconstruction method is used in the present paper for computing high-order derivatives of the gradient variables in the new HNS system; derivatives of the primitive variables are provided by the gradient variables and their derivatives that are already available. As a result the rdg schemes will have the same number of degrees of freedom as a FV scheme applied to the HNS formulation or conventional P DG schemes for the NS equations i.e. 2 and this is true for arbitrarily high-order accuracy. In this paper we present a new HNS formulation and demonstrate that it enables an efficient construction of highorder rdg schemes for solving the compressible NS equations with 2 degrees of freedom for arbitrarily high-order accuracy. The paper focuses on schemes up to fourth-order accuracy in the inviscid terms which eliminate both dispersion and dissipation errors from a second-order scheme and are expected to provide significant improvements for high-fidelity turbulent-flow simulations. An implicit solver is used to solve the residual equations based on an approximate or exact Jacobian with a pseudo time term which is equivalent to the first-order backward Euler pseudotime integration scheme. The pseudo time step is taken to be infinity in many cases and thus it can be considered 2 of 28
3 simply as an implicit nonlinear solver. The HNS rdg solver is developed for the 3D NS equations and implemented in a 3D code but here we focus on 2D problems and perform accuracy verification studies and solve a simple viscousflow problem to demonstrate that the HNS rdg schemes provide accurate solutions for a realistic viscous boundary conditions. The remainder of this paper is organized as follows. Section 2 describes the target NS equations Section 3 presents the new HNS formulation and the construction of dissipation matrices are given Section 4 discusses the discretization methods Section 5 describes the implicit solver Section 6 presents numerical results and Section 7 closes the paper with remarks. II. We consider the compressible NS equations: Compressible Navier-Stokes Equations Downloaded by Hiroaki Nishikawa on July DOI:.254/ U t + F k x k = G k x k (2. ρ ρu ρv ρw ρu ρu 2 + p ρuv ρuw U = ρv F x = ρuv F y = ρv 2 + p F z = ρvw (2.2 ρw ρuw ρvw ρw 2 + p ρe u(ρe + p v(ρe + p w(ρe + p τ xx τ xy G x = τ yx G y = τ yy (2.3 τ zx τ zy uτ xx + vτ xy + wτ xz q x uτ yx + vτ yy + wτ yz q y τ xz G z = τ yz (2.4 τ zz uτ zx + vτ zy + wτ zz q z where v = (u v w is the velocity vector t is the physical time ρ is the density p is the pressure e is the specific total energy and the viscous stress tensor τ and the heat flux q are given under Stokes hypothesis by τ xx τ xy τ xz τ = τ yx τ yy τ yz = 2 3 µ(div vi + µ ( grad v + (grad v T q = τ zx τ zy τ zz q x q y q z = µ grad T (2.5 P r(γ where I is the identity matrix T is the temperature γ is the ratio of specific heats P r is the Prandtl number µ is the viscosity which is taken as a constant in this paper and the superscript T denotes the transpose. In this paper we focus on the spatial discretization and therefore consider only the spatial part in the rest of the paper: F k x k = G k x k. (2.6 For unsteady problems the physical time derivative can be discretized in time and then treated as a source term as demonstrated in Ref.[]. For many practical problems implicit-time stepping schemes would be required and then the steady solver developed in this paper will serve as a nonlinear solver for a system of unsteady residual equations that needs to be solved at every time step. Unsteady extensions will be discussed in a subsequent paper. 3 of 28
4 III. New Hyperbolic Navier-Stokes System: HNS2G The HNS formulation used in Ref.[2] uses the velocity gradients scaled by the viscosity the heat fluxes and the density gradients as additional variables thus resulting in 2 equations and thus termed HNS2. As we will illustrate later the scaled velocity gradients and the heat fluxes introduce complications in the efficient construction of hyperbolic rdg schemes. To simply the construction we propose to form a hyperbolic viscous system with the gradients of the primitive variables as additional variables. 3.. HNS2G The new HNS formulation is given in the form of a preconditioned conservative system: where U and F k are redefined here as U P τ + F k = S (3. x k Downloaded by Hiroaki Nishikawa on July DOI:.254/ ( P = diag U = [ρ ρv ρe g u g v g w h r] T (3.2 ρv T ν ρ r T ρv v + pi µ v τ τ τ v T (ρe + p (µ v τ τ τv T µ h γ(γ h F = ui vi wi T I ρi (3.3 τ τ τ = 3 4 (g + gt 2 tr(gi µ v = 4 3 µ µ h = γ µ (3.4 P r S = [ g u g v g w h r] T (3.5. (3.6 T v T v T v T v T v T v T v T v T v T h T h T h T r T r T r Note that the system is equivalent to the target steady Navier-Stokes equations (2.6 in the pseudo steady state or as soon as the pseudo time term is dropped with an infinite pseud-otime step and the additional variables are equivalent to the gradients of the primitive variables (thus called the gradient variables: r = ρ = (ρ x ρ y ρ z (3.7 g ux g uy g uz g = [gu T gv T gw] T T = [( u T ( v T ( w T ] T = g vx g vy g vz (3.8 g wx g wy g wz h = T = (h x h y h z. (3.9 To emphasize that this system uses the gradients of the primitive variables as additional variables we term the new system HNS2G. The preconditioned conservative system simplifies the construction of a conservative scheme and a numerical flux [8] and the viscous system is now hyperbolic in the pseudo time as we will discuss below. The relaxation times T r T v and T h are defined as T r = L2 r ν ρ T v = ρl2 r µ v T h = ρl2 r µ h (3. where the length scale L r is defined as in Refs.[ 5] and ν ρ is an artificial viscosity associated with the artificial hyperbolic mass diffusion added to the continuity equation [7 2]. 4 of 28
5 The flux projected along an arbitrary direction n = (n x n y n z T is split into the inviscid part viscous part and artificial hyperbolic mass diffusion part: F n = F x n x + F y n y + F z n z = F i n + F v n + F a n (3. Downloaded by Hiroaki Nishikawa on July DOI:.254/ where where τ n = ρu n ν ρu n v T ρ r n + pn µ v τ τ τ n u n (ρe + p µ h µ v τ nv γ(γ h n F i n = F v n = un F a n = vn (3.2 wn T n ρn τ nx τ ny τ nz u n = v n = un x + vn y + wn z h n = h n = h x n x + h y n y + h z n z (3.3 = τ n = τ xx n x + τ xy n y + τ xz n z τ yx n x + τ yy n y + τ yz n z τ zx n x + τ zy n y + τ zz n z τ nv = τ n v = τ nx u + τ ny v + τ nz w. (3.4 The inviscid part is hyperbolic so it can be estimated by using classic Riemann solvers. The HLLC flux [43] is used in the present work. In the HNS2G system the viscous part is also hyperbolic as shown in the next section Dissipation Matrix for Viscous Flux An important difference between HNS2 and HNS2G is that the viscosity does not appear in the flux for HNS2 but does appear in the case of HNS2G. Therefore the flux Jacobian which is required to construct a dissipation term involves derivatives of the viscosity. This is the complication we get in return for a simplfied high-order scheme construction as discussed later. However we can avoid the complication by considering the flux Jacobian with a frozen viscosity [27]: where ( A v n = P Fv n U = A A 2 A 3 A 4 A 5 A 2 A = µ v ρ τ nv (3.5 (3.6 µ v ρ τ nx µ v ρ τ ny µ v ρ τ nz µ v n x 3 4 µ vn y 3 4 µ vn z A 2 = 2 µ vn y 3 4 µ vn x (3.7 2 µ vn z 3 4 µ vn x µ v ( un x + 2 vn y + 2 wn z 3 4 µ v (un y + vn x 3 4 µ v (un z + wn x 5 of 28
6 Downloaded by Hiroaki Nishikawa on July DOI:.254/ µ vn y 2 µ vn x A 3 = 3 4 µ vn x µ v n y 3 4 µ vn z (3.8 2 µ vn z 3 4 µ vn y 3 ( 4 µ v (un y + vn x µ v 2 un x vn y + 2 wn z 3 4 µ v (vn z + wn y 3 4 µ vn z 2 µ vn x A 4 = 3 4 µ vn z 2 µ vn y ( µ vn x 3 4 µ vn y µ v n z 3 4 µ v (un z + wn x 3 ( 4 µ v (vn z + wn y µ v 2 un x + 2 vn y wn z A 5 = (3.2 µ hn x µ hn y µ hn z γ(γ γ(γ γ(γ 6 of 28
7 Downloaded by Hiroaki Nishikawa on July DOI:.254/ un x n x un y n y un z n z vn x n x vn y n y vn z n z wn x n x A 2 = wn y n y (3.2 wn z n z γ(γ (v 2 e n x γ(γ un x γ(γ vn x γ(γ wn x γ(γ n x ρt h ρt h ρt h ρt h ρt h γ(γ (v 2 e n y γ(γ un y γ(γ vn y γ(γ wn y γ(γ n y ρt h ρt h ρt h ρt h ρt h γ(γ (v 2 e n z γ(γ un z γ(γ vn z γ(γ wn z γ(γ n z ρt h ρt h ρt h ρt h ρt h where the viscosity is frozen in computing the derivative Fv n U. The eigenvalues of the matrix can be obtained as where λ = a nv λ 2 = a nv λ 3 = a mv λ 4 = a mv (3.22 λ 5 = a lv λ 6 = a lv λ 7 = a h λ 8 = a h λ = (3.23 a nv = νv 3νv νh a mv = a lv = a h =. (3.24 T v 4T v T h corresponding to the normal viscous shear viscous and heating waves respectively with a set of linearly independent left and right eigenvectors l and r. Therefore the system is hyperbolic. Then the upwind flux can be constructed 7 of 28
8 with the dissipation term acting on two states U R and U L Downloaded by Hiroaki Nishikawa on July DOI:.254/ A v n U = 8 λ k (l T k Ur k k= ρ[(a nv a mv u n n + a mv v] ρ(a nv a mv u n u n + ρa mv (v v + ρa ( h T γ(γ + µ2 v τnn τ nn ρa h P r n + + τ n τ n τ nn τ nn P r m + [( a nv 4 = 3 a mv τ nn n x + 43 ] a mv τ nx n [( a nv 4 3 a mv τ nn n y + 43 ] a mv τ ny n [( a nv 4 3 a mv τ nn n z + 43 ] a mv τ nz n [ a h h n + γ(γ µ ( v τnn u n µ h P r n + + τ ] n v τ nn u n n P r m + where U = U R U L P r n = a nv a h (3.25 P r m = a mv a h. (3.26 This dissipation term is similar to the one for HNS2 [2] but has different scalings for the viscous stresses and heat fluxes. Finally this dissipation term will be multiplied by P from the left to cancel the effect of the preconditioning [8]. To ensure strong coupling among the system in the discrete level an artificial hyperbolic dissipation is added to the upwind scheme for the viscous flux [7]. Consider the following hyperbolic diffusion system: U P τ + Fd x x + Fd y y + Fd z z = Sd ( of 28
9 where F d x = µ v g ux µ v g vx µ v g wx u v w F d y = µ v g uy µ v g vy µ v g wy u v w F d z = µ v g uz µ v g vz µ v g wz u v w S = g ux g uy g uz g vx g vy g vz g wx g wy g wz. (3.28 The flux projected along an arbitrary direction n is given by F d n = µ v g un µ v g vn µ v g wn un vn wn (3.29 where g un = g ux n x + g uy n y + g uz n z g vn = g vx n x + g vy n y + g vz n z g wn = g wx n x + g wy n y + g wz n z. (3.3 9 of 28 Downloaded by Hiroaki Nishikawa on July DOI:.254/
10 Downloaded by Hiroaki Nishikawa on July DOI:.254/ The flux Jacobian with a frozen viscosity is given by µ v n µ v n µ v n un x nx un y ny un z nz vn x nx A av n = P Fd vn n U = y ny vn z nz. (3.3 wn x nx wn y ny wn z nz The eigenvalues are λ...3 = a nv λ = a nv λ =. (3.32 The resultant absolute Jacobian is u v w n 2 x n x n y n x n z n x n y n 2 y n y n z n x n z n y n z n 2 z n 2 x n x n y n x n z A av n = a nv n x n y n 2 y n y n z n x n z n y n z n 2 z. n 2 x n x n y n x n z n x n y n 2 y n y n z n x n z n y n z n 2 z (3.33 of 28
11 Therefore the artificial dissipation term is given by u ρ + (ρu v ρ + (ρv w ρ + (ρw A av n U = a nv g un n = a nv g vn n g wn n Finally the upwind viscous flux across an interface is given by ρ u ρ v ρ w g un n g vn n g wn n. (3.34 F v n(u L U R = 2 (Fv n(u L + F v n(u R 2 P ( A v n + A av n (U R U L. (3.35 Downloaded by Hiroaki Nishikawa on July DOI:.254/ Dissipation Matrix for Artificial Mass Diffusion Flux The remaining component is the artificial mass diffusion term which is shared by HNS2 and HNS2G and required to obtain high-order accurate density gradients [7 2]. The corresponding flux is given by and the flux Jacobian F a n = ( ν ρ ρn x ρn y ρn z T (3.36 A a n = P Fa n U = ν ρ n x ν ρ n y ν ρ n z O T ρ T ρ T ρ nx ny nz The dissipation term is given as described in Ref.[7 2] by where A a n U = a ρ Therefore the artificial mass diffusion flux is given by a ρ = ρ r n n. (3.37 (3.38 νρ T ρ. (3.39 F a n(u L U R = 2 (Fa n(u L + F a n(u R 2 P A a n (U R U L. (3.4 IV. Discretization In this section we present efficient rdg schemes for the HNS2G system. First we discuss the underlying DG discretization which can be very expensive for HNS2G. Then we illustrate the recycling procedure that dramatically reduces the number of discrete unknwons/equations by taking advantage of the gradient variables. Finally we introduce the rdg method to bring the discretization discretization further down and develop very efficient HNS rdg schemes. 4.. Discontinous Galerkin We begin by formulating the discontinuous Galerkin discretization with a test function Ψ over an element Ω e with a boundary Γ e : multiply the HNS2G system by Ψ and integrate by parts over the element P d dτ ΨU h dω = Ω e Ω e [ ] Ψ F k (U h + ΨS(U h dω ΨF k (U h n k dγ Ψ V p h x. (4. k Γ e of 28
12 where P can be assumed to be constant over the element without affecting accuracy and thus has been taken out of the integral and U h is the solution we seek for in the broken Soblev space V n h : U h V p h V p h = {v h [L 2 (Ω] k : v h Ωi [V k p ] Ω i Ω} (4.2 where k is the dimension of the unknown vector and V p is the space of all polynomials of degree p. The solution can be expressed as a polynomial of the form: U h (x y z τ = C(x y zv(τ (4.3 where C is a basis matrix and V is a vector of unknown polynomial coefficients. For example in the P DG method we have U h = (ρ ρv ρe g x g y g z h x h y h z r x r y r z T (4.4 Downloaded by Hiroaki Nishikawa on July DOI:.254/ V = ( ρ ρv ρe g x g y g z h x h y h z r x r y r z T (4.5 C = I C sub C sub C sub C sub = (4.6 where g x = (g ux g vx g wx g y = (g uy g vy g wy g z = (g uz g vz g wz the overbar indicates the cell-averaged value =.5(x max x min =.5(y max y min and =.5(z max z min. x max x min y max y min z max and z min are the maximum and minimum coordinates in the element Ω e in x- y- and z- directions. Note that the variables (g u g v g w have been rearranged in the vector U h for convenience; compare Equation (4.4 and Equation (3.2. Choose the basis matrix for the test function: Ψ = C T (4.7 then P d [ ] C C T T CVdΩ = F k (U h + C T S(U h dω C T F k (U h n k dγ. (4.8 dτ Ω e Ω e x k Γ e which can be written into the semi-discrete form: where P M = P Ω C T CdΩ R(U h = e P M dv dτ = R(U h (4.9 Ω e [ ] C T F k (U h + C T S(U h dω C T F k (U h n k dγ (4. x k Γ e The interface flux is evaluated at each quadrature point for a given pair of states U L and U R computed by the solution polynomials from two adjacent elements by the HLLC flux [43] for the inviscid part and the sum of the upwind viscous flux (3.35 and the artificial mass diffusion flux (3.4 for the viscous part. This discretization results in a large number of degrees of freedom (or equivalently discrete equations per cell with derivatives added to the vector V: 8 for P 2 for P 2 etc. for the HNS2G system. As shown in the previous papers [ ] it is possible to dramatically reduce the number of degrees of freedom by taking advantage of the gradient variables and a reconstruction technique as described in the next sections Efficient Hyperbolic Discontinuous Galerkin: Hyperbolic DG The first step to reduce the required degrees of freedom is to use or recycle the gradient variables to represent high-order derivatives of the primal variables. This is a simple procedure for model equations as demonstrated in Ref.[ ]. However it is somewhat complicated for the HNS system because the primal variables are the conservative variables (ρ ρu ρv ρw ρe and the gradient variables correspond to the gradients of the primitive 2 of 28
13 variables (ρ u v w T. Thus the derivatives of the conservative variables cannot be directly expressed by these gradient variables; it can be expressed indirectly but the mass matrix will then depend on the solution. Such a construction may be explored in future but in this paper we consider a simplified procedure where we keep the P conservative variables in U h construct high-order polynomials in the primitive variables and use them to evaluate the flux and source term: [ ] C T R(W h = F k (W h + C T S(W h dω C T F k (W h n k dγ (4. x k Γ e Ω e where W h is a set of higher-order polynomials of the primitive variables: W h (x y z τ = C w (x y zv w (τ (4.2 Downloaded by Hiroaki Nishikawa on July DOI:.254/ with W = (ρ r u g u v g v w g w T h and C w is a basis matrix for the primitive variables. In this way higherorder polynomials can be easily constructed for the primal variables (ρ u v w T by directly using the gradient variables (r g u g v g w h. For example for the density we can construct the following: where ρ r x r y r z B = h = B B 2 B 3 B 4 ρ r x r y r z B 2 = x x c B 3 = y y c B 4 = z z c. (4.3 The solution has just been upgraded to linear for the density without introducing additional unknowns. Similarly for the x-velocity we have u B B 2 B 3 B 4 ρu/ρ g ux g uy = g ux g uy (4.4 g uz h g uz and for the temperature T h x h y h z B B 2 B 3 B 4 = h (γ ρ γ [ ] ρe 2ρ (ρu2 + ρv 2 + ρw 2 h x h y h z (4.5 Note that the cell-averaged values of the primitive variables are algebraically computed from the conservative variables that we keep in U h. It follows from these examples that the basis matrix C w is a block diagonal matrix: C 4 4 C 4 4 C w = C 4 4 C 4 4 (4.6 C 4 4 where C 4 4 = B B 2 B 3 B 4 (4.7 for V w defined by V w = (ρ u v w T T (4.8 3 of 28
14 where with ρ = (ρ r x r y r z u = (u g ux g uy g uz v = (v g vx g vy g vz (4.9 w = (w g wx g wy g wz T = (T h x h y h z (4.2 u = ρu/ρ v = ρv/ρ w = ρw/ρ T = (γ ρ γ [ ρe ] 2ρ (ρu2 + ρv 2 + ρw 2. (4.2 This construction is called HNS(PP+P. Note that HNS(PP+P has 2 unknowns and achieve secondorder accuracy in the inviscid terms which is comparable to the P DG scheme applied to the original Navier-Stokes equations. In effect the above efficient construction brings the number of unknowns from 8 to 2 for the P level of accuracy. However as pointed out in Ref.[28] it is first-order accurate in the viscous limit. Further improvement comes from the rdg method. Downloaded by Hiroaki Nishikawa on July DOI:.254/ Efficient Hyperbolic Reconstructed Discontinuous Galerkin: Hyperbolic rdg In the rdg method higher-order accuracy is achieved by upgrading the order of the polynomial used to evaluate the flux and source terms in the residual through higher-order derivative reconstructions. The efficient construction in the previous section is in fact a variant of the rdg method but no reconstruction techniques were needed because the derivatives are directly available through the gradient variables. To further improve accuracy we now consider reconstructing higher-order derivatives by a reconstruction technique and construct higher-order reconstructed polynomials Wh R: Wh R (x y z τ = C R w(x y zv R w(τ (4.22 where C R w is a basis matrix and V R w the unknown coefficient vector containing now reconstructed derivatives. For example if the second derivatives ρ c xx ρ c yy ρ c zz ρ c xy ρ c xz ρ c yz are reconstructed from the gradient variables r x r y and r z at the centroid as indicated by the superscript then we can construct a higher-order polynomial for the density and the density gradient variables as follows: ρ r x r y r z R h B B 2 B 3 B 4 B 5 B 6 B 7 B 8 B 9 B B 2 B 3 B 4 = B 3 B 2 B 4 B 4 B 2 B 3 ρ r x r y r z ρ c xx 2 ρ c yy 2. (4.23 ρ c zz 2 ρ c xy ρ c xz ρ c yz where B 5 = 2 B2 2 Ω e Ω e 2 B2 2dΩ B 6 = 2 B2 3 Ω e Ω e 2 B2 3dΩ B 7 = 2 B2 4 Ω e Ω e 2 B2 4dΩ B 8 = B 2 B 3 B 2 B 3 dω B 9 = B 2 B 4 B 2 B 4 dω B = B 3 B 4 B 3 B 4 dω. Ω e Ω e Ω e Ω e Ω e Ω e The density is now represented by a quadratic polynomial and the corresponding gradient variables by linear polynomials. Applying a similar reconstruction to other primitive variables we arrive at the following basis matrix: C 4 C 4 Wh R (x y z = C 4 C 4 VR (4.24 C 4 4 of 28
15 where C 4 = B B 2 B 3 B 4 B 5 B 6 B 7 B 8 B 9 B B 2 B 3 B 4 B 3 B 2 B 2 B 4 B 4 B 3 (4.25 for V R w defined by V R w = (ρ 2 u 2 v 2 w 2 T 2 T (4.26 where ρ 2 = (ρ r x r y r z ρ c xx 2 ρ c yy 2 ρ c zz 2 ρ c xy ρ c xz ρ c yz (4.27 Downloaded by Hiroaki Nishikawa on July DOI:.254/ and similarly for other variables. It is emphasized that the reconstructed polynomials are used only in the flux and source terms not in the weak formulation: P M dv dτ = R(WR h (4.28 and thus the mass matrix M is still based on the underlying polynomial basis matrix C defined for the conservative variables. For the derivative reconstruction we employ the variational reconstruction for computing the derivatives [4 42]. The variational reconstruction for computing the gradient of a function is denoted by PP(VR and it is applied to the gradient variables here to obtain their gradients that are equivalent to the second derivatives of the density. The resulting scheme is called HNS(PP2+PP(VR which now outperforms the conventional second-order P DG: it has the same 2 degrees of freedom and achieves third-order accuracy in the inviscid terms second-order in the viscous terms and second-order in the gradients. An even more efficient scheme is HNS(PP3+PP2(VR with ρ r x r y r z R h B B 2 B 3 B 4 B 5 B 2 B 6 B 7 B B 3 8 B B 4 9 B B 4 = B + B 2 B5 c B 2 + B 3 B6 c B 3 + B 4 B7 c B 4 + B 2 B8 c + B 3 B5 c B 5 + B 2 B9 c + B 4 B5 c B 6 + B 3 B8 c + B 2 B6 c B 7 + B 2 B c + B 3 B9 c + B 4 B8 c B 8 + B 4 B9 c + B 2 B7 c B 9 + B 3 B c + B 4 B6 c B 2 + B 4 B c + B 3 B7 c B 3 B 4 B 2 B 2 B 3 B 5 B 6 B 7 B 8 B 5 B 9 B 5 B 6 B 8 B B 9 B 8 B 7 B 9 B B 6 B 7 B T ρ r x r y r z ρ c xx 2 ρ c yy 2 ρ c zz 2 ρ c xy ρ c xz ρ c yz ρ c xxx 3. (4.29 ρ c yyy 3 ρ c zzz 3 ρ c xxy 2 ρ c xxz 2 ρ c yyx 2 ρ c xyz ρ c zzx 2 ρ c yyz 2 ρ c zzy 2 5 of 28
16 Scheme DoF/Cell Inviscid Viscous Gradient DG(P 2 O(h 2 O(h 2 O(h DG(P2 5 O(h 3 O(h 3 O(h 2 HNS(PP+P 2 O(h 2 O(h O(h HNS(PP2+PP(VR 2 O(h 3 O(h 2 O(h 2 HNS(PP3+PP2(VR 2 O(h 4 O(h 3 O(h 3 Table. Comparison of the HNS schemes and conventional DG schemes for the 3D NS equations. DoF/Cell is the degrees of freedom within each cell. Downloaded by Hiroaki Nishikawa on July DOI:.254/ where B = 6 B3 2 Ω e B 4 = 2 B2 2B 3 Ω e B 8 = 2 B2 4B 2 Ω e Ω e 6 B3 2dΩ B 2 = 6 B3 3 Ω e Ω e 2 B2 2B 3 dω B 5 = 2 B2 2B 4 Ω e Ω e 6 B3 3dΩ B 3 = 6 B3 4 Ω e B 7 = B 2 B 3 B 4 Ω e Ω e B 2 B 3 B 4 dω Ω e 2 B2 4B 2 dω B 9 = 2 B2 3B 4 Ω e Ω e 2 B2 2B 4 dω B 6 = 2 B2 3B 2 Ω e Ω e 2 B2 3B 4 dω B 2 = 2 B2 4B 3 Ω e Ω e 6 B3 4dΩ Ω e 2 B2 3B 2 dω Ω e 2 B2 4B 3 dω (4.3 and the superscript c indicates that the basis function is evaluated by the centroid coordinates. In this case we reconstruct not only the gradient of the gradient variables but also their second derivatives which correspond to third-order derivatives of the density so that a cubic polynomial can be reconstructed for the density. Here a higher-order version of the variational reconstruction PP2(VR is used to compute the gradients and second-derivatives of the gradient variables. The resulting scheme achieves with the same 2 degrees of freedom fourth-order accuracy in the inviscid terms third-order accuracy in the viscous terms and the solution gradients. Compare it with a conventional P DG scheme which requires 2 degrees of freedom for second- and first-order accurate solutions and gradients respectively; or a conventional P 2 DG which requires 5 degrees of freedom for third- and second-order accurate solutions and gradients respectively. Table shows the comparison between the HNS-rDG schemes and conventional DG schemes applied to the original NS system. HNS(PP+P is comparable to the conventional P DG scheme but has a lower order of accuracy in the viscous terms. This is improved by HNS(PP2+PP(VR which increases the order of accuracy by one order for all while keeping the same 2 degrees of freedom: one-order higher accuracy for the inviscid and gradients and matching the second-order accuracy in the viscous terms with the conventional P DG scheme. A further improvement is achieved by HNS(PP3+PP2(VR: it gives fourth-order accuracy in the inviscid terms and third-order accuracy in the viscous terms and gradients with the same 2 degrees of freedom which is more efficient than the conventional P 2 DG scheme. Thee HNS-rDG schemes may be considered as high-order FV schemes developed in the rdg framework. It can be systematically extended to higher-order in more than one way. For example HNS(PP4+P2P3(VR may be developed by a linear (instead of quadratic VR reconstruction method or HNS(PP4+PP3(VR by a cubic VR reconstruction method. In the former the basis matrix used in the weak formulation may no longer be diagonal and the scheme will be a mixed FV-DG scheme. Furthermore if the primitive variables are used as a set of working variables in the weak formulation HNS-rDG schemes will be again mixed-fv-dg schemes. Such variants should be explored in future and compared with the schemes developed in this paper. V. Implicit Solver We discretize the pseudo time term by the first-order backward Euler scheme: P M V Vm τ = R(ŴR h (5. 6 of 28
17 where τ is a pseudo time step and m is the previous pseudo time level and solve the resulting global system of residual equations for V by the following implicit method: V k+ = V k + V (5.2 where k is the iteration counter the correction V is given as a solution to the linearized system: ( P M τ and R is the residual without the reconstruction: R V = V R(ŴR h (5.3 R = R(W h. (5.4 Downloaded by Hiroaki Nishikawa on July DOI:.254/ The Jacobians for the inviscid HLLC and viscous upwind fluxes are computed exactly. However the residual Jacobian ignores the VR reconstructed part and thus approximate for rdg schemes; the solver reduces to Newton s method only when the pseudo time step is taken as infinity and no reconstruction is performed. The linear system is relaxed by the symmetric Gauss-Seidel relaxation scheme. In this study only one symmetric Gauss-Seidel relaxation is performed per iteration. VI. Numerical Results 6.. 2D Manufactured solution We verify the HNS rdg schemes implemented in a 3D code by methods of manufactured solutions for 2D problems on hexahedral regular irregular and highly distorted prismatic mesh. The second levels of each set of meshes are shown in Figure. The following functions are made the exact solutions by introducing source terms into the Navier-Stokes equations: ρ = cr + crs sin(crx x + cry y u = cu + cus sin(cux x + cuy y v = cv + cvs sin(cvx x + cvy y (6. w = p = cp + cps sin(cpx x + cpy y A constant viscosity µ = ρ inf v inf L/Re is used here for a verification purpose. The subscript denotes the free stream values. The free stream state is taken as ρ inf = cr and v inf = sqrt(cu 2 + cv 2 here. cr crs crx cry cu cus cux cuy cv cvs cvx cvy cp cps cpx and cpy are all taken as constants to generate a smooth solution over the domain. The Reynolds number is taken as and 5. The artificial mass diffusion coefficient ν ρ is taken as h 4 where h = /nelem for 2D problem and L r is taken as /2π. The log of residual and L2 error norms are plotted in Figure 2 and Figure 3. From these two figures we can see that although the HNS(PP+P method will come to a steady solution when all the residuals drop to at least four orders of magnitude HNS(PP3+PP2(VR will come to a steady state after all the residuals drop at least to seven orders of magnitude. From the comparison we can see that the solution becomes steady after all the residuals drop to the order of magnitude of the log of error. The mesh convergence is shown in Figure 4 and Figure 5. For viscous limit when Re = all the primitive variables and auxiliary variables converge to first second and third order of accuracy for HNS(PP+P HNS(PP2+PP(VR and HNS(PP3+PP2(VR respectively. However when it comes to inviscid limit for Re = 5 the gradient variables still converge to first second and third but an over convergence occurs for primitive variables as they converge to second- third- and forth-order of accuracy for these three methods respectively. Table 2 shows the number of iterations versus /h for Re = on hexahedral mesh. The number of iterations is taken at five six seven orders of magnitude residual reductions for HNS(PP+P HNS(PP2+PP(VR and HNS(PP3+PP2(VR respectively. From the table we can see that the iteration increases linearly with /h for all the three schemes demonstrating the convergence acceleration feature of the hyperbolic method; it increases quadratically for conventional viscous schemes in such a low-reynolds-number problem Couette flow Couette flow is a laminar viscous flow between two parallel plates. The bottom plate with a fixed temperature (T is stationary while the top plate with a fixed temperature (T is moving at a constant speed of U in the positive x- 7 of 28
18 /h PP+P PP2+PP(VR PP3+PP2(VR Table 2. Iteration numbers versus /h Re=. direction. The distance between two plates is H. Under some simplifications with the costant µ and low speed to ensure the nearly incompressible effectsthe analytical solution is given [44] as u = U z w = (6.2 H Downloaded by Hiroaki Nishikawa on July DOI:.254/ T = T + z H + z H (T T + P r U 2 z 2 Cp H ( z (6.3 H p = constant ρ = p (6.4 RT where P r denotes the Prandtl number R is the gas constant and Cp is the specific heat capacity for constant pressure. In our numerical experiment H = 2 T =.8 and T =.85. The Mach number for the moving wall is taken as.. Reynold number is taken as with a constant viscosity µ =.. The computational domain is a rectangle ( x 2H z H. For this test case take ν r = 5. Weak Dirichlet boundary conditions are given for all the boudary faces. Hexahedral prismatic and hybrid mesh are used here with 5 successively refined grids. The second level meshes of each type are shown in Figure 6. The mesh convergence is shown in Figure 7. Errors in the temperature and its gradient in z direction converge at first-order for HNS(PP+P method. However a super convergence occurs when it comes to HNS(PP2+PP(VR and HNS(PP3+PP2(VR methods giving third- and forth- order of accuracy for temperature and its graident in z direction. VII. Conclusions A new formulation of hyperbolic Navier-Stokes method HNS2G suitable to the reconstructed discontinuous Galerkin framework has been developed. In this formulation the hyperbolic viscous system is established by introducing the gradients of density velocity and temperature as the additional variables. The high order derivatives of the gradient variables are obtained by performing a variational reconstruction to avoid the redundant variables. Three discretization schemes HNS(PP+P HNS(PP2+PP(VR and HNS(PP3+PP2(VR have been developed and the residual equations are solved by an implicit solver. The methods have been verified by the method of manufactured solutions and the designed orders of accuracy have bene confirmed for all the primitive variables and gradient variables including one-order higher-order of accuracy in the inviscid limit. In future work other high-order variants should be explored the methods will be applied to more realistic three-dimensional flow test cases on arbitrary grids and extended to unsteady problems. Acknowledgements This work has been partially funded by the U.S. Army Research Office under the contract/grant number W9NF with Dr. Mattew Munson as the program manager. 8 of 28
19 Y Y Z X Z X Y Y Z X Z X Downloaded by Hiroaki Nishikawa on July DOI:.254/ Figure. Second level of each set of mesh used for manufactured solution 9 of 28
20 Downloaded by Hiroaki Nishikawa on July DOI:.254/ resi_ resi_v resi_p resi_ xyz resi_v xyz resi_t xyz error_ error_v error_p error_ xyz error_v xyz error_t xyz 5 iter (HNS(PP+P 5 iter (HNS(PP2+PP_VR 5 5 iter (HNS(PP3+PP2_VR iter (HNS(PP+P iter (HNS(PP2+PP_VR iter (HNS(PP3+PP2_VR 2 of 28
21 Downloaded by Hiroaki Nishikawa on July DOI:.254/ iter (HNS(PP+P iter (HNS(PP2+PP_VR 2 4 iter (HNS(PP3+PP2_VR iter (HNS(PP+P iter (HNS(PP2+PP_VR iter (HNS(PP3+PP2_VR Figure 2. Residual and L2 error versus iteration manufactured solution Re =. 2 of 28
22 Downloaded by Hiroaki Nishikawa on July DOI:.254/ iter (HNS(PP+P 5 iter (HNS(PP2+PP_VR iter (HNS(PP3+PP2_VR 5 5 iter (HNS(PP+P iter (HNS(PP2+PP_VR iter (HNS(PP3+PP2_VR 22 of 28
23 Downloaded by Hiroaki Nishikawa on July DOI:.254/ iter (HNS(PP+P iter (HNS(PP2+PP_VR 2 4 iter (HNS(PP3+PP2_VR iter (HNS(PP+P iter (HNS(PP2+PP_VR iter (HNS(PP3+PP2_VR Figure 3. Residual and L2 error versus iteration for manufactured solution Re = of 28
24 3 3 error_ 5 7 HNS(PP+PRe=e hexa HNS(PP2+PP_VRRe=e hexa HNS(PP3+PP2_VRRe=e hexa HNS(PP+P Re=e regular HNS(PP2+PP_VR Re=e regular HNS(PP3+PP2_VR Re=e regular HNS(PP+P Re=e irregular HNS(PP2+PP_VR Re=e irregular HNS(PP3+PP2_VR Re=e irregular HNS(PP+P Re=e distorted HNS(PP2+PP_VR Re=e distorted HNS(PP3+PP2_VR Re=e distorted Slope Slope 2 Slope 3 Slope Log(h error_ xyz Log(h Downloaded by Hiroaki Nishikawa on July DOI:.254/ error_v error_p Log(h Log(h error_v xyz error_t xyz Log(h Log(h Figure 4. Mesh convergence for manufactured solution Re =. 24 of 28
25 3 3 error_ 5 7 HNS(PP+PRe=e5 hexa HNS(PP2+PP_VRRe=e5 hexa HNS(PP3+PP2_VRRe=e5 hexa HNS(PP+P Re=e5 regular HNS(PP2+PP_VR Re=e5 regular HNS(PP3+PP2_VR Re=e5 regular HNS(PP+P Re=e5 irregular HNS(PP2+PP_VR Re=e5 irregular HNS(PP3+PP2_VR Re=e5 irregular HNS(PP+P Re=e5 distorted HNS(PP2+PP_VR Re=e5 distorted HNS(PP3+PP2_VR Re=e5 distorted Slope Slope 2 Slope 3 Slope Log(h error_ xyz Log(h Downloaded by Hiroaki Nishikawa on July DOI:.254/ error_v error_p Log(h Log(h error_v xyz error_t xyz Log(h Log(h Figure 5. Mesh convergence for manufactured solution Re = of 28
26 Z Z Y X Y X Z Y X Downloaded by Hiroaki Nishikawa on July DOI:.254/ error_t Figure 6. Second level hexahedral prismatic and hybrid mesh for Couette flow HNS(PP+PRe=hexa HNS(PP2+PP_VRRe=hexa HNS(PP3+PP2_VRRe=hexa HNS(PP+PRe=pris HNS(PP2+PP_VRRe=pris HNS(PP3+PP2_VRRe=pris HNS(PP+PRe=hybrid HNS(PP2+PP_VRRe=hybrid HNS(PP3+PP2_VRRe=hybrid Slope Slope 2 Slope 3 Slope Log(h Log(h error_h z HNS(PP+PRe=hexa HNS(PP2+PP_VRRe=hexa HNS(PP3+PP2_VRRe=hexa HNS(PP+PRe=pris HNS(PP2+PP_VRRe=pris HNS(PP3+PP2_VRRe=pris HNS(PP+PRe=hybrid HNS(PP2+PP_VRRe=hybrid HNS(PP3+PP2_VRRe=hybrid Slope Slope 2 Slope 3 Slope 4 Figure 7. Mesh convergence for Couette flow Re =. 26 of 28
27 Downloaded by Hiroaki Nishikawa on July DOI:.254/ References [] Bassi F. and Rebay S. GMRES discontinuous Galerkin solution of the compressible Navier-Stokes equations Discontinuous Galerkin Methods: Theory Computations and Applications edited by B. C. G. Karniadakis and C. W. Shu Vol. Springer 2 pp [2] Guo Liu J. and Wang Shu C. A High-Order Discontinuous Galerkin Method for 2D Incompressible Flows J. Comput. Phys. Vol. 6 2 pp [3] Luo H. Luo L. Nourgaliev R. Mousseau V. and Ding N. A Reconstructed Discontinuous Galerkin Method for the Compressible Navier-Stokes Equations on Arbitrary Grids J. Comput. Phys. Vol. 229 No. 9 2 pp [4] Nishikawa H. A First-Order System Approach for Diffusion Equation. I: Second Order Residual Distribution Schemes J. Comput. Phys. Vol pp [5] Nishikawa H. A First-Order System Approach for Diffusion Equation. II: Unification of Advection and Diffusion J. Comput. Phys. Vol pp [6] Nishikawa H. First Second and Third Order Finite-Volume Schemes for Navier-Stokes Equations Proc. of 7th AIAA Theoretical Fluid Mechanics Conference AIAA Aviation and Aeronautics Forum and Exposition 24 AIAA Paper Atlanta GA 24. [7] Nishikawa H. Alternative Formulations for First- Second- and Third-Order Hyperbolic Navier-Stokes Schemes Proc. of 22nd AIAA Computational Fluid Dynamics Conference AIAA Paper Dallas TX 25. [8] Nishikawa H. New-Generation Hyperbolic Navier-Stokes Schemes: O(/h Speed-Up and Accurate Viscous/Heat Fluxes Proc. of 2th AIAA Computational Fluid Dynamics Conference AIAA Paper Honolulu Hawaii 2. [9] Liu Y. and Nishikawa H. Third-Order Inviscid and Second-Order Hyperbolic Navier-Stokes Solvers for Three-Dimensional Inviscid and Viscous Flows 46th AIAA Fluid Dynamics Conference AIAA Paper Washington D.C. 26. [] Liu Y. and Nishikawa H. Third-Order Edge-Based Hyperbolic Navier-Stokes Scheme for Three-Dimensional Viscous Flows 23rd AIAA Computational Fluid Dynamics Conference AIAA Paper Denver Colorado 27. [] Liu Y. and Nishikawa H. Third-Order Inviscid and Second-Order Hyperbolic Navier-Stokes Solvers for Three-Dimensional Unsteady Inviscid and Viscous Flows 55th AIAA Aerospace Sciences Meeting AIAA Paper Grapevine Texas 27. [2] Nakashima Y. Watanabe N. and Nishikawa H. Hyperbolic Navier-Stokes Solver for Three-Dimensional Flows 54th AIAA Aerospace Sciences Meeting AIAA Paper 26- San Diego CA 26. [3] Nagao T. Hashimoto A. and Sato T. Study on Time Evolution Method for Hyperbolic Navier-Stokes System 56th AIAA Aerospace Sciences Meeting AIAA Paper Kissimmee Florida 28. [4] Wu Z. and Wang D. The Development of a Hyperbolic RANS System for Analysing Turbomachinery Flow Field Proc. of Shanghai 27 Global Power and Propulsion Forum GPPS-7-98 Shanghai China 27. [5] Nishikawa H. and Liu Y. Hyperbolic Advection-Diffusion Schemes for High-Reynolds-Number Boundary-Layer Problems J. Comput. Phys. Vol pp [6] Nishikawa H. and Roe P. L. Third-Order Active-Flux Scheme for Advection Diffusion: Hyperbolic Diffusion Boundary Condition and Newton Solver Computers and Fluids Vol pp [7] Nishikawa H. and Liu Y. Hyperbolic Navier-Stokes Method for High-Reynolds-Number Boundary-Layer Flows 55th AIAA Aerospace Sciences Meeting AIAA Paper 27-8 Grapevine Texas 27. [8] Baty H. and Nishikawa H. Hyperbolic Method for Magnetic Reconnection Process in Steady State Magnetohydrodynamics Mon. Not. R. Astron. Soc. Vol pp [9] Kawashima R. Komurasaki K. and Schonherr T. A Flux-Splitting Method for Hyperbolic-Equation System of Magnetized Electron Fluids in Quasi-Neutral Plasmas J. Comput. Phys. Vol pp [2] Mazaheri A. Ricchiuto M. and Nishikawa H. A First-Order Hyperbolic System Approach for Dispersion J. Comput. Phys. Vol. 32 September 26 pp [2] Lou J. Liu X. Luo H. and Nishikawa H. Reconstructed Discontinuous Galerkin Methods for Hyperbolic Diffusion Equations on Unstructured Grids 55th AIAA Aerospace Sciences Meeting AIAA Paper 27-3 Grapevine Texas 27. [22] Nishikawa H. On Hyperbolic Merthod for Diffusion with Discontinous Coefficients J. Comput. Phys. Vol pp [23] Nishikawa H. First Second and Third Order Finite-Volume Schemes for Advection-Diffusion J. Comput. Phys. Vol pp [24] Nishikawa H. First- Second- and Third-Order Finite-Volume Schemes for Diffusion J. Comput. Phys. Vol pp of 28
28 Downloaded by Hiroaki Nishikawa on July DOI:.254/ [25] Lee E. Ahn H. T. and Luo H. Cell-Centered High-order Hyperbolic Finite Volume Method for Diffusion Equation on Unstructured Grids J. Comput. Phys. Vol pp [26] Chamarthi A. S. Wang Z. and Kawashima R. Weighted Nonlinear Schemes for Magnetized Electron Fluid in Quasineutral plasma 56th AIAA Aerospace Sciences Meeting AIAA Paper Kissimmee Florida 28. [27] Lou J. Li L. Luo H. and Nishikawa H. First-Order Hyperbolic System Based Reconstructed Discontinuous Galerkin Methods for Nonlinear Diffusion Equations on Unstructured Grids 56th AIAA Aerospace Sciences Meeting AIAA Paper Kissimmee Florida 28. [28] Lou J. Li L. Liu X. Luo H. and Nishikawa H. Reconstructed Discontinuous Galerkin Methods Based on First-Order Hyperbolic System for Advection-Diffusion Equations 23rd AIAA Computational Fluid Dynamics Conference AIAA Paper Denver Colorado 27. [29] Lou J. Luo H. and Nishikawa H. Discontinuous Galerkin Methods for Hyperbolic Advection-Diffusion Equation on Unstructured Grids Proc. of The 9th International Conference on Computational Fluid Dynamics Istanbul Turkey 26. [3] Lou J. Li L. Luo H. and Nishikawa H. Reconstructed discontinuous Galerkin methods for linear advectiondiffusion equations based on first-order hyperbolic system J. Comput. Phys. Vol pp [3] Nishikawa H. and Nakashima Y. Dimensional Scaling and Numerical Similarity in Hyperbolic Method for Diffusion J. Comput. Phys. Vol pp [32] Liu X. Xia Y. Cheng J. and Luo H. Development and assessment of a reconstructed discontinuous Galerkin method for the compressible turbulent flows on hybrid grids AIAA Paper [33] Liu X. Xuan L. Luo H. and Xia Y. A reconstructed discontinuous Galerkin method for the compressible Navier-Stokes equations on three-dimensional hybrid gridss Comput. Fluids Vol pp [34] Liu X. Xia Y. and Hong Luo L. X. A comparative study of Rosenbrock-type and implicit Runge-Kutta time integration for discontinuous Galerkin method for unsteady 3D compressible Navier-Stokes equations Commun. Compt. Phys. Vol. 2 No pp [35] Xia Y. Liu X. Luo H. and Nourgaliev R. A third-order implicit discontinuous Galerkin method based on a Hermite WENO reconstruction for time-accurate solution of the compressible NavierStokes equations Int. J. Numer. Meth. Fluids Vol. 79 No pp [36] Cheng J. Liu T. and Luo H. A hybrid reconstructed discontinuous Galerkin method for compressible flows on arbitrary grids Comput. Fluids Vol pp [37] Xia Y. Lou J. Hong L. Edwards J. and Mueller F. OpenACC acceleration of an unstructured CFD solver based on a reconstructed discontinuous Galerkin method for compressible flows Int. J. Numer. Meth. Fluids Vol. 78 No pp [38] Xia Y. Luo H. Frisbey M. and Nourgaliev R. A Set of Parallel Implicit Methods for a Reconstructed Discontinuous Galerkin Method for Compressible Flows on 3D Hybrid Grid Comput. Fluids Vol pp [39] Wang C. Cheng J. Berndt M. N.Carlson N. and HongLuo Application of Nonlinear Krylov Acceleration to a Reconstructed Discontinuous Galerkin Method for Compressible Flows Vol pp [4] Liu X. Xia Y. and Luo H. A Reconstructed Discontinuous Galerkin Method for Compressible Turbulent Flows on 3D curved Grids Comput. Fluids Vol pp [4] Li L. Liu X. Lou J. Luo H. Nishikawa H. and Ren Y. A Finite Volume Method Based on Variational Reconstruction for Compressible Flows on Arbitrary Grids 23rd AIAA Computational Fluid Dynamics Conference AIAA Paper Denver Colorado 27. [42] Li L. Liu X. Lou J. Luo H. Nishikawa H. and Ren Y. A Discontinuous Galerkin Method Based on Variational Reconstruction for Compressible Flows on Arbitrary Grids 56th AIAA Aerospace Sciences Meeting AIAA Paper Kissimmee Florida 28. [43] Toro E. F. Spruce M. and Speares W. Restoration of the Contact Surface in the HLL-Riemann Solver Shock Waves Volume 39 Issue 994 pp [44] R. Kannan and Z.J. Wang The direct discontinuous Galerkin (DDG viscous flux scheme for the high order spectral volume method Computers & Fluids Vol pp of 28
High-Order Hyperbolic Navier-Stokes Reconstructed Discontinuous Galerkin Method
AIAA SciTech Forum 7- January 29 San Diego California AIAA Scitech 29 Forum.254/6.29-5 High-Order Hyperbolic Navier-Stokes Reconstructed Discontinuous Galerkin Method Lingquan Li Jialin Lou and Hong Luo
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationExplicit Hyperbolic Reconstructed Discontinuous Galerkin Methods for Time-Dependent Problems
AIAA AVIATION Forum June 25-29 218 Atlanta Georgia 218 Fluid Dynamics Conference 1.2514/6.218-427 Explicit Hyperbolic Reconstructed Discontinuous Galerkin Metods for Time-Dependent Problems Jialin Lou
More informationActive Flux for Advection Diffusion
Active Flux for Advection Diffusion A Miracle in CFD Hiroaki Nishikawa National Institute of Aerospace! NIA CFD Seminar! August 25, 2015 In collaboration with the University of Michigan Supported by NASA
More informationAccuracy-Preserving Source Term Quadrature for Third-Order Edge-Based Discretization
Preprint accepted in Journal of Computational Physics. https://doi.org/0.06/j.jcp.07.04.075 Accuracy-Preserving Source Term Quadrature for Third-Order Edge-Based Discretization Hiroaki Nishikawa and Yi
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
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 informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationThird-Order Active-Flux Scheme for Advection Diffusion: Hyperbolic Diffusion, Boundary Condition, and Newton Solver
Preprint accepted for publication in Computers and Fluids, 16 http://d.doi.org/1.116/j.compfluid.15.1. Third-Order Active-Flu Scheme for Advection Diffusion: Hyperbolic Diffusion, Boundary Condition, and
More information1 PART1: Bratu problem
D9: Advanced Numerical Analysis: 3 Computer Assignment 3: Recursive Projection Method(RPM) Joakim Möller PART: Bratu problem. Background The Bratu problem in D is given by xx u + yy u + λe u =, u Γ = ()
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 informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
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 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 informationGoverning Equations of Fluid Dynamics
Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1
AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01
More informationJacobian-Free Newton Krylov Discontinuous Galerkin Method and Physics-Based Preconditioning for Nuclear Reactor Simulations
INL/CON-08-14243 PREPRINT Jacobian-Free Newton Krylov Discontinuous Galerkin Method and Physics-Based Preconditioning for Nuclear Reactor Simulations International Conference on Reactor Physics, Nuclear
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationDepartment of Mathematics University of California Santa Barbara, Santa Barbara, California,
The Ghost Fluid Method for Viscous Flows 1 Ronald P. Fedkiw Computer Science Department Stanford University, Stanford, California 9435 Email:fedkiw@cs.stanford.edu Xu-Dong Liu Department of Mathematics
More informationNumerical Solution of Partial Differential Equations governing compressible flows
Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
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 informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University
More informationThe hybridized DG methods for WS1, WS2, and CS2 test cases
The hybridized DG methods for WS1, WS2, and CS2 test cases P. Fernandez, N.C. Nguyen and J. Peraire Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics, MIT 5th High-Order
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More informationSlip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics
Slip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics J.J.W. van der Vegt and H. van der Ven Nationaal Lucht- en Ruimtevaartlaboratorium National
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 informationModule 2: Governing Equations and Hypersonic Relations
Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume
More informationApplication of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations. Xiaoyue Shen
Application of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations Xiaoyue Shen Supervisor: Prof. David W. Zingg A thesis submitted in conformity with the requirements
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
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 informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More information2. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations
. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations We need to review the governing equations of fluid mechanics before examining the methods of computational aerodynamics in detail.
More informationComputation of NACA0012 Airfoil Transonic Buffet Phenomenon with Unsteady Navier-Stokes Equations
5th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 9-2 January 22, Nashville, Tennessee AIAA 22-699 Computation of NACA2 Airfoil Transonic Buffet Phenomenon with
More informationAE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says
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 informationOptimizing Runge-Kutta smoothers for unsteady flow problems
Optimizing Runge-Kutta smoothers for unsteady flow problems Philipp Birken 1 November 24, 2011 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany. email:
More informationImproved Seventh-Order WENO Scheme
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Eposition 4-7 January 2, Orlando, Florida AIAA 2-45 Improved Seventh-Order WENO Scheme Yiqing Shen Gecheng Zha Dept.
More informationWALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID
WALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID Zia Ghiasi sghias@uic.edu Dongru Li dli@uic.edu Jonathan Komperda jonk@uic.edu Farzad
More informationON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY ON INTERIOR NODAL SETS
6th European Conference on Computational Mechanics (ECCM 6 7th European Conference on Computational Fluid Dynamics (ECFD 7 11 15 June 018, Glasgow, UK ON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value
More informationCompressible Navier-Stokes (Euler) Solver based on Deal.II Library
Compressible Navier-Stokes (Euler) Solver based on Deal.II Library Lei Qiao Northwestern Polytechnical University Xi an, China Texas A&M University College Station, Texas Fifth deal.ii Users and Developers
More informationHigh-Order Multidomain Spectral Difference Method for the Navier-Stokes Equations on Unstructured Hexahedral Grids
COMMUICATIOS I COMPUTATIOAL PHYSICS Vol. 2, o. 2, pp. 310-333 Commun. Comput. Phys. April 2007 High-Order Multidomain Spectral Difference Method for the avier-stokes Equations on Unstructured Hexahedral
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
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 informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationConvection Heat Transfer
Convection Heat Transfer Department of Chemical Eng., Isfahan University of Technology, Isfahan, Iran Seyed Gholamreza Etemad Winter 2013 Heat convection: Introduction Difference between the temperature
More informationImmersed boundary technique for compressible flow simulations on semi-structured meshes
Center for Turbulence Research Annual Research Briefs 2005 71 Immersed boundary technique for compressible flow simulations on semi-structured meshes By M. de Tullio AND G. Iaccarino 1. Motivation and
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 informationarxiv: v1 [physics.comp-ph] 21 Nov 2018
Implicit High-Order Flux Reconstruction Solver for High-Speed Compressible Flows Ray Vandenhoeck a,, Andrea Lani b a Department of Mechanical Engineering, KU Leuven, Celestijnenlaan 300, 3001, Leuven,
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationEffect of Applied Magnetic Field on Shock Boundary Layer Interaction
University of Kentucky UKnowledge Mechanical Engineering Faculty Publications Mechanical Engineering 1-2012 Effect of Applied Magnetic Field on Shock Boundary Layer Interaction Ovais U. Khan University
More informationA gas-kinetic theory based multidimensional high-order method for the compressible Navier Stokes solutions
Acta Mech. Sin. 2017) 334):733 741 DOI 10.1007/s10409-017-0695-2 RESEARCH PAPER A gas-kinetic theory based multidimensional high-order method for the compressible Navier Stokes solutions Xiaodong Ren 1
More informationThe Equivalent Differential Equation of the Gas-Kinetic BGK Scheme
The Equivalent Differential Equation of the Gas-Kinetic BGK Scheme Wei Gao Quanhua Sun Kun Xu Abstract The equivalent partial differential equation is derived for the gas-kinetic BGK scheme (GKS) [Xu J.
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 informationSTEADY AND UNSTEADY 2D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS FLOW. Radka Keslerová, Karel Kozel
Conference Applications of Mathematics 1 in honor of the th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 1 STEADY AND UNSTEADY D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationAssessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow
Assessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering - University
More informationAalto University School of Science and Technology CFD-group/ Department of Applied Mechanics. MEMO No CFD/MECHA DATE: March 15, 2012
Aalto University School of Science and Technology CFD-group/ Department of Applied Mechanics MEMO No CFD/MECHA-2-212 DATE: March 15, 212 TITLE The Effect of Free-Stream Turbulence Parameters on the SST
More informationControl Theory Approach to Aero Shape Optimization
Control Theory Approach to Aero Shape Optimization. Overview. Finite Volume Method 3. Formulation of Optimal Design Problem 4. Continuous Adjoint Approach 5. Discrete Adjoint Approach Overview of Adjoint
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 informationNumerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin
More informationStabilization of Explicit Flow Solvers Using a Proper Orthogonal Decomposition Technique
50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 09-12 January 2012, Nashville, Tennessee AIAA 2012-1096 Stabilization of Explicit Flow Solvers Using a Proper
More informationURANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver
URANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver *Yohan Choe 1), Hyeongjun Kim 1), Chongam Kim 2) 1), 2) Department of Aerospace Engineering, Seoul National
More informationThe Shallow Water Equations
The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow
More informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationA Speed-Up Strategy for Finite Volume WENO Schemes for Hyperbolic Conservation Laws
J Sci Comput (011) 46: 359 378 DOI 10.1007/s10915-010-9407-9 A Speed-Up Strategy for Finite Volume WENO Schemes for Hyperbolic Conservation Laws Fei Teng Li Yuan Tao Tang Received: 3 February 010 / Revised:
More informationFINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS
FINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS Lucia Catabriga Department of Computer Science, Federal University of Espírito Santo (UFES) Av. Fernando Ferrari,
More informationAeroelastic Analysis of Engine Nacelle Strake Considering Geometric Nonlinear Behavior
Aeroelastic Analysis of Engine Nacelle Strake Considering Geometric Nonlinear Behavior N. Manoj Abstract The aeroelastic behavior of engine nacelle strake when subjected to unsteady aerodynamic flows is
More informationarxiv: v4 [physics.flu-dyn] 23 Mar 2019
High-Order Localized Dissipation Weighted Compact Nonlinear Scheme for Shock- and Interface-Capturing in Compressible Flows Man Long Wong a, Sanjiva K. Lele a,b arxiv:7.895v4 [physics.flu-dyn] 3 Mar 9
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationPerformance tuning of Newton-GMRES methods for discontinuous Galerkin discretizations of the Navier-Stokes equations
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2685 Performance tuning of Newton-GMRES methods for discontinuous Galerkin
More informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More informationEigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations
NASA/CR-1998-208741 ICASE Report No. 98-51 Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations David L. Darmofal Massachusetts Institute of Technology, Cambridge,
More informationGT UNSTEADY SIMULATION OF A TWO-STAGE COOLED HIGH PRESSURE TURBINE USING AN EFFICIENT NON-LINEAR HARMONIC BALANCE METHOD
Proceedings of ASME Turbo Expo 213: Turbine Technical Conference and Exposition GT213 June 3-7, 213, San Antonio, Texas, USA GT213-94574 UNSTEADY SIMULATION OF A TWO-STAGE COOLED HIGH PRESSURE TURBINE
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
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 informationSimulation of unsteady muzzle flow of a small-caliber gun
Advances in Fluid Mechanics VI 165 Simulation of unsteady muzzle flow of a small-caliber gun Y. Dayan & D. Touati Department of Computational Mechanics & Ballistics, IMI, Ammunition Group, Israel Abstract
More information