Runge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter
|
|
- Easter Greene
- 5 years ago
- Views:
Transcription
1 Runge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter Jun Zhu, inghui Zhong, Chi-Wang Shu 3 and Jianxian Qiu, College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing, Jiangsu 6, P.R. China. Department of Mathematics, Michigan State University, East Lansing, MI 88, USA. 3 Division of Applied Mathematics, Brown University, Providence, RI 9, USA. School of Mathematical Sciences and Fujian Provincial Key Laboratory of Mathematical Modeling and High-Performance Scientific Computation, iamen University, iamen, Fujian 365, P.R. China. Abstract. In this paper, we propose a new type of weighted essentially non-oscillatory (WENO) limiter, which belongs to the class of Hermite WENO (HWENO) limiters, for the Runge-Kutta discontinuous Galerkin (RKDG) methods solving hyperbolic conservation laws. his new HWENO limiter is a modification of the simple WENO limiter proposed recently by Zhong and Shu [9]. Both limiters use information of the DG solutions only from the target cell and its immediate neighboring cells, thus maintaining the original compactness of the DG scheme. he goal of both limiters is to obtain high order accuracy and non-oscillatory properties simultaneously. he main novelty of the new HWENO limiter in this paper is to reconstruct the polynomial on the target cell in a least square fashion [8] while the simple WENO limiter [9] is to use the entire polynomial of the original DG solutions in the neighboring cells with an addition of a constant for conservation. he modification in this paper improves the robustness in the computation of problems with strong shocks or contact discontinuities, without changing the compact stencil of the DG scheme. Numerical results for both one and two dimensional equations including Euler equations of compressible gas dynamics are provided to illustrate the viability of this modified limiter. AMS subject classifications: 65M6, 35L65 Key words: Runge-Kutta discontinuous Galerkin method; HWENO limiter; conservation law. Corresponding author. addresses: zhujun@nuaa.edu.cn (J. Zhu), zhongxh@math.msu.edu (. Zhong), shu@dam.brown.edu (C.-W. Shu), jxqiu@xmu.edu.cn (J. Qiu) Global Science Preprint
2 Introduction In this paper, we are interested in solving the hyperbolic conservation law { ut + f(u) x =, u(x,)=u (x), (.) and its two-dimensional version { ut + f(u) x +g(u) y =, u(x,y,)=u (x,y), (.) using the Runge-Kutta discontinuous Galerkin (RKDG) methods [ 7], where u, f(u) and g(u) can be either scalars or vectors. It is not an easy task to solve (.) and (.) since solutions may contain discontinuities even if the initial conditions are smooth. Discontinuous Galerkin (DG) methods can capture weak shocks and other discontinuities without further modification. However, for problems with strong discontinuous solutions, the numerical solution has significant oscillations near discontinuities, especially for high order methods. A common strategy to control these spurious oscillations is to apply a nonlinear limiter. One type of limiters is based on the slope methodology, such as the minmod type limiters [ 7], the moment based limiter [] and an improved moment limiter [3]. hese limiters do control the oscillations well, however they may degrade accuracy when mistakenly used in smooth regions of the solution. Another type of limiters is based on the weighted essentially nonoscillatory (WENO) methodology [ 3,7], which can achieve both high order accuracy and non-oscillatory properties. he WENO limiters introduced in [8,, 3] and the Hermite WENO limiters in [9, ] belong to this type. hese limiters are designed based on the WENO finite volume methodology which require a wider stencil for higher order schemes. herefore, it is difficult to implement them for multi-dimensional problems, especially on unstructured meshes. An alternative family of DG limiters which serves at the same time as a new PDE-based limiter, as well as a troubled cells indicator, was introduced by Dumbser et al. []. More recently, a particularly simple and compact WENO limiter, which utilizes fully the advantage of DG schemes in that a complete polynomial is available in each cell without the need of reconstruction, is designed for RKDG schemes in [9]. he two major advantages of this simple WENO limiter are the compactness of its stencil, which contains only immediate neighboring cells, and the simplicity in implementation, especially for unstructured meshes [3]. However, it was observed in [9] that the limiter might not be robust enough for problems containing very strong shocks or low pressure problem, especially for higher order polynomials, for example the blast wave problems [3,8] and the double rarefaction wave problem [6], making it necessary to apply additional positivity-preserving limiters [7] in such situation. In order to overcome this difficulty, without compromising the advantages of compact stencil and simplicity of linear weights, we present a modification of the limiter in the step of preprocessing the
3 3 polynomials in the immediate neighboring cells before applying the WENO reconstruction procedure. his preprocessing is necessary to maintain strict conservation, and is designed in [9] to be a simple addition of a constant to make the cell average of the preprocessed neighboring cell polynomial in the target cell matching the original cell average. In this paper, a more involved least square process [8] is used in this step. he objective is to achieve strict conservation while maintaining more information of the original neighboring cell polynomial before applying the WENO procedure. Numerical experiments indicate that this modification does improve the robustness of the limiter. his paper is organized as follows: In Section, we provide a brief review of RKDG for one dimensional and two dimensional cases. In Section 3, we provide the details of the new HWENO limiter for one dimensional scalar and system cases. In Section, we provide the details of the HWENO limiter for two dimensional scalar and system cases. We demonstrate the performance of our HWENO limiter with one and two dimensional numerical examples including Euler equations of compressible gas dynamics in Section 5. Concluding remarks are given in Section 6. Review of RKDG methods In this subsection, we give a brief review of the RKDG methods for solving one and two dimensional conservation laws. One dimensional case. Given a partition of the computational domain consisting of cells I j =[x j,x j+ ], j=,,n, denote the cell center by x j = (x j +x j+ ), and the cell size by x j = x j+ x j. he DG method has its solution as well as the test function space given by Vh k={v(x):v(x) I j P k (I j )}, where P k (I j ) denotes the set of polynomials of degree at most k defined on I j. he semi-discrete DG method for solving (.) is defined as follows: find the unique function u h (,t) Vh k such that for j=,,n, (u h ) t vdx f(u h )v x dx+ ˆf j+ v(x I j I j+ j ) ˆf j v(x + )= (.) j holds for all test functions v Vh k, where u± = u j+ h (x ±,t) are the left and right limits of j+ the discontinuous solution u h at the interface x j+ and ˆf j+ = ˆf(u,u + is a monotone j+ j+ ) flux for the scalar case and an exact or approximate Riemann solver for the system case. wo dimensional case. Given a partition of the computational domain consisting of rectangular mesh consisting of cells I ij =[x i,x i+ ] [y j,y j+ ], for i=,,n x and j=,,n y with the cell sizes x i+ x i = x i, y j+ y j = y j and cell centers(x i,y j )= ( (x i+ +x i ), (y j+ +y j )). We now give the new test function space W k h ={p : p I ij P k (I ij )} as the polynomial spaces of degree of at most k on the cell I ij. he semi-discrete DG method for solving (.) is defined as follows: find the unique function u W k h such
4 that, for all i N x and j N y, (u h ) t vdxdy= I i,j f(u h )v x dxdy I i,j + g(u h )v y dxdy I i,j,y)dy+ ˆfi I j ˆfi+ I j (y)v(x i+ ĝ j+ (x)v(x,y )dx+ I j+ i (y)v(x +,y)dy i ĝ i (x)v(x,y + )dx, I j i (.) hold for all the test function v Wh k, where the hat terms are again numerical fluxes. he semi-discrete schemes (.) and (.) can be written as u t = L(u) where L(u) is the spatial discretization operator. hey can be discretized in time by a non-linearly stable Runge-Kutta time discretization [5], e.g. the third-order version: u () = u n + tl(u n ), u () = 3 un + u() + tl(u() ), (.3) u n+ = 3 un + 3 u() + 3 tl(u() ). As in [9], to apply a nonlinear limiter for the RKDG methods, for simplicity, we take the forward Euler time discretization of the semi-discrete scheme (.) as an example. Starting from a solution u n h Vk h at time level n (u h is taken as the L projection of the analytical initial condition into Vh k ). We would like to limit it to obtain a new function u n,new before advancing it to next time level. hat is: find u n+ h Vh k, such that, for j=,,n, Ij u n+ h u n,new h t vdx f(u n,new h )v x dx+ I j n,new ˆf v(x ) j+ j+ n,new ˆf v(x + (.) j j )= holds for all test functions v Vh k. he limiting procedure to go from un h discussed in the following sections. to un,new h will be 3 New HWENO limiter in one dimension In this section, we describe the details of the modified HWENO reconstruction procedure as a limiter for the RKDG method in the one dimensional scalar and system cases.
5 5 3. he troubled cell indicator in one dimension An important component of the simple WENO limiter in [9] is the identification of troubled cells, which are cells that may contain discontinuities and in which the WENO limiter is applied. We will use the KRCF shock detection technique developed in [5] to detect troubled cells. We divide the boundary of the target cell I j into two parts: I j and I + j, where the flow is into (v n<, n is the normal vector to I j and v is the velocity) and out of (v n>) I j, respectively. Here, we define v, taking its value from inside the cell I j as f (u) and take u as the indicator variable for the scalar case. he target cell I j is identified as a troubled cell when I (u h (x,t) Ij u h (x,t) Il )ds j > C x k+ k, (3.) j I j u h (x,t) Ij where C k is a constant, usually, we take C k = as in [5]. Here I l, for l = j or j+, denotes the neighboring cell sharing the end point I j with I j. u h is the numerical solution corresponding to the indicator variable(s) and u h (x,t) Ij is the standard L norm in the cell I j. 3. HWENO limiting procedure in one dimension: scalar case In this subsection, we present the details of the HWENO limiting procedure for the scalar case. he idea of this new HWENO limiter is to reconstruct a new polynomial on the troubled cell I j which is a convex combination of three polynomials: the DG solution polynomial on this cell and the modified DG solution polynomials on its two immediate neighboring cells. he modification procedure is in a least square fashion [8]. he construction of the nonlinear weights in the convex combination coefficients follows the classical WENO procedure. Assume I j is a troubled cell. Step.. Denote the DG solution polynomials of u h on I j,i j+,i j as polynomials p (x), p (x) and p (x), respectively. Now we want to find the modified version of p (x), denoted as p (x) on the cell I j in a least square fashion [8]. he modification procedure is defined as follows: p (x) is the solution of the minimization problem { } min (φ(x) p (x)) dx, (3.) φ(x) P k (I j ) I j subject to φ= p, where φ= φ(x)dx, p = p (x)dx. x j I j x j I j Here and below denotes the cell average of the function on the target cell.
6 6 Similarly, p (x) is the solution of the minimization problem { } min (φ(x) p (x)) dx φ(x) P k (I j+ ) I j+ (3.3) subject to φ= p. For notational consistency we denote p (x)= p (x). he final nonlinear HWENO reconstruction polynomial p new (x) is now defined by a convex combination of these modified polynomials: p new (x)=ω p (x)+ω p (x)+ω p (x). (3.) According to the modification procedure, it is easy to prove that p new has the same cell average and order of accuracy as p if the weights satisfy ω +ω +ω =. he convex combination coefficients ω l,l=,, follow the classical WENO procedure [,, 3]. We discuss it in the following steps. Step.. We choose the linear weights denoted by γ,γ,γ. As in [9], we have used complete information of the three polynomials p (x), p (x), p (x) in the three cells I j, I j+, I j, hence we do not have extra requirements on the linear weights in order to maintain the original high order accuracy. he linear weights can be chosen to be any set of positive numbers adding up to one. he choice of these linear weights is then solely based on the consideration of a balance between accuracy and ability to achieve essentially nonoscillatory shock transitions. In all of our numerical tests, following the practice in [9, 9], we take γ =.998 and γ = γ =.. Step.3. We compute the smoothness indicators, denoted by β l, l=,,, which measure how smooth the functions p l (x), l=,,, are on the target cell I j. he smaller these smoothness indicators are, the smoother the functions are on the target cell. We use the similar recipe for the smoothness indicators as in [3]: β l = k l= ( x l j I j l! d l ) dx l p l(x) dx, l=,,. (3.5) Step.. We compute the non-linear weights based on the smoothness indicators: ω ω i = i l= ω, ω l = γ l l (ε+β l ). (3.6) Here ε is a small positive number to avoid the denominator to become zero. We take ε= 6 in our computation. Step.5. he final nonlinear HWENO reconstruction polynomial is given by (3.), i.e. u new h Ij = p new (x)=ω p (x)+ω p (x)+ω p (x). It is easy to verify that p new (x) has the same cell average and order of accuracy as the original one p (x), on the condition that l= ω l=.
7 7 3.3 HWENO limiting procedure in one dimension: system case In this subsection, we present the details of the HWENO limiting procedure for one dimensional systems. Consider equation (.) where u and f(u) are vectors with m components. In order to achieve better non-oscillatory qualities, the HWENO reconstruction limiter is used with a local characteristic field decomposition. In this paper, we consider the following Euler system with m=3. t ρ ρµ E + x ρµ ρµ + p µ(e+ p) =, (3.7) where ρ is the density, µ is the x-direction velocity, E is the total energy, p is the pressure and γ=. in our test cases. We denote the Jacobian matrix as f (u), where u=(ρ,ρµ,e). We then give the left and right eigenvector matrices of such Jacobian matrix as: and L j (u)= R j (u)= B +µ/c B µ+/c B B B µ B B µ/c B µ /c µ c µ µ+c H cµ µ / H+cµ B, (3.8), (3.9) where c= γp/ρ, B =(γ )/c, B = B µ / and H=(E+ p)/ρ. Assume the troubled cell I j is detected by the KRCF technique [5] by using (3.), where v=µ is the velocity and again v takes its value from inside the cell I j. We take both the density ρ and the total energy E as the indicator variables. Denote p, p, p to be the DG polynomial vectors, corresponding to I j s two immediate neighbors and itself. hen we perform the characteristic-wise HWENO limiting procedure as follows: Step. Compute L j = L j (ū j ) and R j = R j (ū j ) as defined in (3.8) and (3.9), where ū j is the cell average of u on the cell I j. Step.. Project the polynomial vectors p, p, p into the characteristic fields p l = L j p l, l=,,, each of them being a 3-component vector and each component of the vector is a k-th degree polynomial. Step.3. Perform Step. to Step.5 of the HWENO limiting procedure that has been specified for the scalar case, to obtain a new 3-component vector on the troubled cell I j as p new. Step... Project p new back into the physical space to get the reconstruction polyno- Ij = R j p new mial, i.e. u new h.
8 8 New HWENO limiter in two dimension. he troubled cell indicator in two dimension We use the KRCF shock detection technique developed in [5] to detect troubled cells in two dimensions. We divide the boundary of the target cell I ij into two parts: I ij and I + ij, where the flow is into (v n<, n is the normal vector to I ij) and out of (v n> ) I ij, respectively. Here we define v, taking its value from inside the cell I ij, as the vector ( f (u),g (u)) and take u as the indicator variable for the scalar case. For the Euler system (.), v, again taking its value from inside the cell I ij, is (µ,ν) where µ is the x-direction velocity and ν is the y-direction velocity, and we take both the density ρ and the total energy E as the indicator variables. he target cell I ij is identified as a troubled cell when I (u h (x,y,t) Iij u h (x,y,t) Il )ds ij > C h k+ k, (.) ij I ij u h(x,y,t) Iij where C k is a constant, usually, we take C k = as in [5]. Here we choose h ij as the radius of the circumscribed circle in I ij, and I l,l=(i,j ); (i,j); (i+,j); (i,j+), denote the neighboring cells sharing the edge(s) in I ij. u h is the numerical solution corresponding to the indicator variable(s) and u h (x,y,t) Iij is the standard L norm in the cell I ij.. HWENO limiting procedure in two dimensions: scalar case In this subsection, we give details of the HWENO limiter for the two dimensional scalar case. he idea is similar to the one dimensional case, i.e. to reconstruct a new polynomial on the troubled cell I ij which is a convex combination of the following polynomials: DG solution polynomial on this cell and the modified DG solution polynomials on its immediate neighboring cells. he nonlinear weights in the convex combination coefficients follow the classical WENO procedure. o achieve better non-oscillatory property, the modification procedure now is in a least square fashion [8] with necessary adjustment. Assume I ij is a troubled cell. Step 3.. We select the HWENO reconstruction stencil as S={I i,j, I i,j, I i+,j, I i,j+, I ij }, for simplicity, we renumber these cells as I l,l=,,, and denote the DG solutions on these five cells to be p l (x,y), respectively. Now we want to find the modified version of p l (x,y), denoted as p l (x,y), in a least square fashion [8] with necessary adjustment. he modification procedure not only use the complete information of the polynomial p l (x,y), but also use the cell averages of the polynomials from its other neighbors. For p (x,y), the modification procedure is defined as follows: p (x,y) is the solution
9 9 of the minimization problem { ( min φ(x,y) P k (I ) subject to φ= p and ) ( ) } (φ(x,y) p (x,y)) dxdy + (φ(x,y) p l (x,y))dxdy, I l L I l L ={,,3} {l : p l p <max( p p, p p, p 3 p )}, where p l = p I l l (x,y)dxdy is the cell average of the polynomial p l (x,y) on the cell I l I l and I l is the area of I l. Here and below denotes the cell average of the function on its own associated cell. For this modification, we try to find the polynomial p, which has the same cell average as the polynomial on the troubled cell p, to optimize the distance to p (x,y) and to the cell averages of those useful polynomials on the other neighboring cells. By comparing the distance between the cell averages of the polynomials on the other neighboring cells and the cell average of p on the target cell, if one is not the farthest, then this polynomial is considered useful. Remark: he one dimensional algorithm is consistent with the two dimensional one, because in one dimension we only have two immediate neighbors and hence L= is the empty set. In the extreme case (for example, if p p = p p = p 3 p ), this two dimensional algorithm could degenerate to the one dimensional case and L could also be. Similarly, p (x,y) is the solution of the minimization problem { ( ) ( ) } min (φ(x,y) p (x,y)) dxdy + (φ(x,y) p l (x,y))dxdy, φ(x,y) P k (I ) I l L I l subject to φ= p, where L ={,,3} {l : p l p <max( p p, p p, p 3 p )}. p (x,y) is the solution of the minimization problem { ( min φ(x,y) P k (I ) subject to φ= p, where ) ( ) } (φ(x,y) p (x,y)) dxdy + (φ(x,y) p l (x,y))dxdy, I l L I l L ={,,3} {l : p l p <max( p p, p p, p 3 p )}. p 3 (x,y) is the solution of the minimization problem { ( min φ(x,y) P k (I 3 ) ) ( ) } (φ(x,y) p 3 (x,y)) dxdy + (φ(x,y) p l (x,y))dxdy, I 3 l L I 3 l
10 subject to φ= p, where L 3 ={,,} {l : p l p <max( p p, p p, p p )}. We also define p (x,y)= p (x,y). Step 3.. We choose the linear weights denoted by γ,γ,γ,γ 3,γ. Similar as in the one dimensional case, we put a larger linear weight on the troubled cell and the neighboring cells get smaller linear weights. In all of our numerical tests, following the practice in [9, 9], we take γ =.996 and γ = γ = γ = γ 3 =.. Step 3.3. We compute the smoothness indicators, denoted by β l, l=,...,, which measure how smooth the functions p l (x,y), l=,...,, are on the target cell I ij. We use the similar recipe for the smoothness indicators as in [, 3, ]: β l = k α = I ij α I ij ( ) α α! x α p y α l (x,y) dxdy, l=,,, (.) where α=(α,α ) and α =α +α. Step 3.. We compute the nonlinear weights based on the smoothness indicators. Step 3.5. he final nonlinear HWENO reconstruction polynomial p new (x,y) is defined by a convex combination of the (modified) polynomials in the stencil: u new h Ii,j = p new (x,y) = l= ω l p l (x,y). (.3) It is easy to verify that p new (x,y) has the same cell average and order of accuracy as the original one p (x,y) on the condition that l= ω l=..3 HWENO limiting procedure in two dimensions: system case In this subsection, we present the details of the HWENO limiting procedure for two dimensional systems. Consider equation (.) where u, f(u) and g(u) are vectors with m components. In order to achieve better non-oscillatory qualities, the HWENO reconstruction limiter is used with a local characteristic field decomposition. In this paper, we only consider the following Euler system with m =. t ρ ρµ ρν E + x ρµ ρµ + p ρµν µ(e+ p) + y ρν ρµν ρν + p ν(e+ p) = (.) where ρ is the density, µ is the x-direction velocity, ν is the y-direction velocity, E is the total energy, p is the pressure and γ =. in our test cases. We then give the left and right
11 eigenvector matrices of Jacobian matrices f (u) and g (u) as: B +µ/c B µ+/c B ν Lij x (u)= B ν B B µ B ν B B µ/c B µ /c B ν B, (.5) µ c µ µ+c Rij x (u)= ν ν ν, (.6) µ +ν H cµ ν H+cµ L y ij (u)= B +ν/c B µ B ν+/c B µ B B µ B ν B B ν/c B µ B ν /c B, (.7) µ µ µ R y ij (u)= ν c ν ν+c, (.8) H cν µ µ +ν H+cν where c= γp/ρ, B =(γ )/c, B =B (µ +ν )/ and H=(E+p)/ρ. he troubled cell I ij is detected by the KRCF technique [5] using (.). Denote p l, l=,,, to be the DG polynomial vectors, corresponding to I ij s four immediate neighbors and itself. he HWENO limiting procedure is then performed as follows: Step.. We reconstruct the new polynomial vectors p x,new and p y,new by using the characteristic-wise HWENO limiting procedure: Step.. Compute Lij x = Lx ij (ū ij), L y ij = Ly ij (ū ij), Rij x = Rx ij (ū ij) and R y ij = Rx ij (ū ij) as defined in (.5)-(.8), where again ū ij is the cell average of u on the cell I ij.
12 Step... Project the polynomial vectors p l, l=,,, into the characteristic fields p l x= Lx ij p l and p y l = Ly ij p l, l=,,, each of them being a -component vector and each component of the vector is a k-th degree polynomial. Step..3. Perform Step 3. to Step 3.5 of the HWENO limiting procedure that has been specified for the scalar case, to obtain the new -component vectors on the troubled cell I ij such as p x,new and p y,new. Step... Project p x,new and p y,new into the physical space p x,new = Rij x p x,new and p y,new = R y ij p y,new, respectively. Step.. he final new -component vector on the troubled cell I ij is defined as u new h Iij = px,new + p y,new. 5 Numerical results In this section, we provide numerical results to demonstrate the performance of the HWENO reconstruction limiters for the RKDG methods described in previous sections. We would like to remark that, for some cases with fourth order RKDG methods, the usage of the WENO limiters in [9] might not work well or might even break down, while the new HWENO limiters in this paper could obtain good resolutions. he simple Lax-Friedrichs flux is used in all of our numerical tests. We first test the accuracy of the schemes in one and two dimensional problems. We adjust the constant C k in (3.) or (.) to be a small number. from Example 5. to Example 5., for the purpose of artificially generating a larger percentage of troubled cells in order to test accuracy when the HWENO reconstruction procedure is enacted in more cells. In such test cases, we assume the coarse mesh as Γ h, which has cell size h, and then divide each cell into two (in D) or four (in D) equal smaller cells and denote the associated fine mesh as Γ h, which has smaller mesh size h. he error in the target cell I i is set as e h i. he L and L convergence errors in the computational field are defined as: e h L = N N i= eh i and e h L = max i N e h i, where N is the number of cells in the computational domain. he numerical order of accuracy are given by: O L = log ( eh L ), O L = log ( eh L ). e h L e h L he CFL number is set to be.3 for the second order (P ),.8 for the third order (P ) and. for the fourth order (P 3 ) RKDG methods both in one and two dimensions. Example 5.. We solve the following scalar Burgers equation: ( ) µ µ t + =, x [,], (5.) x with the initial condition µ(x,) =.5+sin(πx) and periodic boundary conditions. We compute the solution up to t=.5/π, when the solution is still smooth. he errors and numerical orders of accuracy for the RKDG method with the HWENO limiter comparing
13 3 ( ) µ able : µ t + =. µ(x,)=.5+sin(πx). Periodic boundary condition. t=.5/π. x L and L errors. RKDG with the HWENO limiter compared to RKDG without limiter. DG with HWENO limiter DG without limiter cells L error order L error order L error order L error order.e-.3e-.38e-.e- 5.6E E E E-.67 P.E-3..6E E-.3 9.E E-..5E-3.3.E-..53E E E E E E-5..69E-.99.3E-5..68E-.99.7E-3.83E-.7E-3.8E-.6E E-3.7.E E-3.7 P.73E E E E E E E E E E E E E E E E E-.E-3.76E-.35E-3.69E E-.3.66E E-.5 P 3 9.6E E E E E E E-8..E E E E-9..3E E-.3 9.E E-. 9.E with the original RKDG method without a limiter are shown in able. We can see that the HWENO limiter keeps the designed order of accuracy, even when a large percentage of good cells are artificially identified as troubled cells. Example 5.. We solve the D Euler equations (3.7). he initial conditions are: ρ(x,) = +.sin(πx), µ(x,)=. and p(x,)=. he computing domain is x [,]. Periodic boundary conditions are applied in this test. he exact solution is ρ(x,t)=+.sin(π(x t)). We compute the solution up to t=. he errors and numerical orders of accuracy of the density for the RKDG method with the HWENO limiter comparing with the original RKDG method without a limiter are shown in able. Similar to the previous example, we can see that the HWENO limiter keeps the designed order of accuracy. Example 5.3. We solve the following scalar Burgers equation in two dimensions: ( ) ( ) µ µ µ t + + =, (x,y) [,] [,], (5.) x y with the initial condition µ(x,y,) =.5+sin(π(x+y)/) and periodic boundary condi-
14 able : D-Euler equations: initial data ρ(x,)=+.sin(πx), µ(x,)=. and p(x,)=. Periodic boundary condition. t=.. L and L errors. RKDG with the HWENO limiter compared to RKDG without limiter. DG with HWENO limiter DG without limiter cells L error order L error order L error order L error order.55e- 3.53E- 3.6E-3 8.9E-3.73E E E-.5.68E-3.63 P 3.8E E-3..73E-. 7.5E E E-.78.3E E E E-5.9.E E E-6..3E-5..9E-6..7E E- 8.E-.3E- 7.86E-.77E E-.99.6E-5 3..E-.9 P.9E E E E E E E E E E E-8 3..E E E E E E-5.83E- 5.E-6 3.6E-5.7E-6..63E E-7..9E P 3.9E-8 7..E E-8..E E E-9..E E E E E E E E-..73E-. 3.5E-.
15 5 ( ( ) µ able 3: µ t + )x + µ =. µ(x,y,)=.5+sin(π(x+y)/). Periodic boundary condition. t=.5/π. L y and L errors. RKDG with the HWENO limiter compared to RKDG without limiter. DG with HWENO limiter DG without limiter cells L error order L error order L error order L error order 3.99E- 3.7E- 3.9E- 3.E- 9.78E-3.3.E E-3..5E-.68 P.3E E E E E-. 9.E E-. 9.6E E-..38E E-..E E-3.77E- 5.E-3.8E- 8.6E-.6.3E E-.6.5E-. P.E E-3.77.E E E-5.96.E-3.6.E-5.96.E E E-.86.8E E-.86.9E-3 8.E-.9E-3 8.9E-.5E E E E P E E E E E E E E E E E E tions in both directions. We compute the solution up to t=.5/π, when the solution is still smooth. he errors and numerical orders of accuracy for the RKDG method with the HWENO limiter comparing with the original RKDG method without limiter are shown in able 3. We again observe good results as in the one-dimensional case. Example 5.. We solve the Euler equations (.). he initial conditions are: ρ(x,y,) = +.sin(π(x+y)), µ(x,y,)=.7, ν(x,y,)=.3, p(x,y,)=. he computational domain is (x,y) [,] [,]. Periodic boundary conditions are applied in both directions. he exact solution is ρ(x,y,t) = +.sin(π(x+y t)). We compute the solution up to t =. he errors and numerical orders of accuracy of the density for the RKDG method with the HWENO limiter comparing with the original RKDG method without a limiter are shown in able. he proposed HWENO limiter again keeps the designed order of accuracy. We now test the performance of the RKDG method with the HWENO limiters for problems containing shocks. From now on, we reset the constant to C k =. For comparison with the RKDG method using the minmod VB limiter, we refer to the results in [,7]. For comparison with the RKDG method using the previous versions of HWENO type limiters, we refer to the results in [8, 3]. Example 5.5. We consider the D Euler equations (3.7) with a Riemann initial con-
16 6 able : D-Euler equations: initial data ρ(x,y,) = +.sin(π(x+y)), µ(x,y,) =.7, ν(x,y,) =.3, and p(x,y,)=. Periodic boundarycondition. t=.. L and L errors. RKDG with the HWENO limiter compared to RKDG without limiter. DG with HWENO limiter DG without limiter cells L error order L error order L error order L error order.6e- 7.98E-.55E-.E- 5.6E E E E-3.5 P 7.83E E E-.8.5E E E-. 9.E E E-5..63E-..89E-5.7.3E E- 5.7E E- 5.5E-3.6E E-.63.E E-.65 P.3E E-.69.5E E E E E E E E E E-6.97.E- 9.86E-.86E-5 6.7E-.85E E E-6.8.E-5.5 P 3.75E-7..9E-6..75E-7.3.9E E-8..55E-7..8E-8..55E E-. 9.7E E-. 9.7E-9.
17 7 dition for the Lax problem: (ρ,µ,p) =(.5,.698,3.58) for x [ 5,); (ρ,µ,p) = (.5,,.57) for x [,5]. he computed density ρ is plotted at t =.3 against the exact solution in Figure and the time history of the troubled cells is shown in Figure Figure : he Lax problem. RKDG with the HWENO limiter. Solid line: the exact solution; squares: numerical solution. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells: Figure : he Lax problem. RKDG with the HWENO limiter. roubled cells. Squares denote cells which are identified as troubled cells subject to the HWENO limiting. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells:. Example 5.6. A higher order scheme would show its advantage when the solution contains both shocks and complex smooth region structures. A typical example for this is the problem of shock interaction with entropy waves [6]. We solve the Euler equations (3.7) with a moving Mach = 3 shock interacting with sine waves in density: (ρ,µ,p) = (3.8573,.69369, ) for x [ 5, ); (ρ,µ,p) =(+.sin(5x),,) for x [,5]. he computed density ρ is plotted at t =.8 against the referenced exact solution which is a converged solution computed by the fifth order finite difference WENO scheme [3] with grid points in Figure 3 and the time history of the troubled cells is shown in Figure.
18 Figure 3: he shock density wave interaction problem. RKDG with the HWENO limiter. Solid line: the exact solution; squares: numerical solution. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells: Figure : he shock density wave interaction problem. RKDG with the HWENO limiter. roubled cells. Squares denote cells which are identified as troubled cells subject to the HWENO limiting. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells:.
19 9 Example 5.7. We consider the interaction of blast waves of Euler equation (3.7) with the initial conditions: (ρ,µ,p) =(,,) for x [,.); (ρ,µ,p) =(,,.) for x [.,.9); (ρ,µ,p) =(,,) for x [.9,.]. he computed density ρ is plotted at t =.38 against the reference exact solution which is a converged solution computed by the fifth order finite difference WENO scheme [3] with grid points in Figure 5 and the time history of the troubled cells is shown in Figure 6. We observe that the new RKDG methods do a better job than that in [9] Figure 5: he blast wave problem. RKDG with the HWENO limiter. Solid line: the exact solution; squares: numerical solution. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells: Figure 6: he blast wave problem. RKDG with the HWENO limiter. roubled cells. Squares denote cells which are identified as troubled cells subject to the HWENO limiting. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells:. Remark. For the P 3 case, the RKDG methods with the WENO limiters specified in [9] used the positivity-preserving limiter to avoid negative density or negative pressure during the time evolution, but the resolution was affected, especially for the smearing of contact discontinuities. he RKDG methods with the new HWENO limiters in this paper can work well without the help from the positivity-preserving limiter.
20 Example 5.8. he Sedov blast wave problem. his problem contains very low density with strong shocks. he exact solution is specified in detail in [, 3]. he computational domain is [,] and initial conditions are: ρ=, µ= and E= everywhere except that the energy in the center cell is the constant 3 x (emulating a δ-function at the center). he final computational time is t =.. he computational results including the density, velocity and pressure pictures are shown in Figure 7 to Figure 9. he time history of the troubled cells is shown in Figure. he RKDG methods with the original WENO limiters in [9] would break down for this problem without the help of positivitypreserving limiters, yet the new limiter in this paper works well Figure 7: he Sedov blast wave problem. RKDG with the HWENO limiter.. Solid line: the exact solution; squares: numerical solution. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells: Velocity Velocity Velocity Figure 8: he Sedov blast wave problem. RKDG with the HWENO limiter. Velocity. Solid line: the exact solution; squares: numerical solution. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells:. Example 5.9. he double rarefaction wave problem [6]. his test case has low pressure and low density regions and is difficult to simulate. he initial conditions are: (ρ,µ,p) =
21 Pressure Pressure Pressure Figure 9: he Sedov blast wave problem. RKDG with the HWENO limiter. Pressure. Solid line: the exact solution; squares: numerical solution. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells: Figure : he Sedov blast wave problem. RKDG with the HWENO limiter. roubled cells. Squares denote cells which are identified as troubled cells subject to the HWENO limiting. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells:.
22 (7,,.) for x [,); (ρ,µ,p) =(7,,.) for x [,]. he final computational time is t =.6. he computational results including the density, velocity and pressure pictures are shown in Figure to Figure 3. he time history of the troubled cells is shown in Figure. For the P 3 case, the RKDG methods with the WENO limiters in [9] do not work for this problem without the help of positivity-preserving limiters, however the same RKDG methods with the new HWENO limiters in this paper produce good results Figure : he double rarefaction wave problem. RKDG with the HWENO limiter.. Solid line: the exact solution; squares: numerical solution. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells: Velocity Velocity Velocity Figure : he double rarefaction wave problem. RKDG with the HWENO limiter. Velocity. Solid line: the exact solution; squares: numerical solution. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells:. Example 5.. Double Mach reflection problem. his model problem is originally from [8]. We solve the Euler equations (.) in a computational domain of [,] [,]. he reflection boundary condition is used at the wall, which for the rest of the bottom boundary (the part from x= to x= 6 ), the exact post-shock condition is imposed. At the top
23 Pressure. Pressure. Pressure Figure 3: he double rarefaction wave problem. RKDG with the HWENO limiter. Pressure. Solid line: the exact solution; squares: numerical solution. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells: Figure : he double rarefaction wave problem. RKDG with the HWENO limiter. roubled cells. Squares denote cells which are identified as troubled cells subject to the HWENO limiting. Left: second order (P ); middle: third order (P ); right: fourth order (P 3 ). Cells:.
24 boundary is the exact motion of the Mach shock. he results shown are at t=.. hree different orders of accuracy for the RKDG with HWENO limiters, k=, k= and k=3 (second order, third order and fourth order), are used in the numerical experiments. he simulation results are shown in Figure 5. he zoomed-in pictures around the double Mach stem to show more details are given in Figure 6. he troubled cells identified at the last time step are shown in Figure 7. Clearly, the resolution improves with an increasing k on the same mesh. Y.5 3 Y.5 3 Y.5 3 Figure 5: Double Mach refection problem. RKDG with the HWENO limiter. 3 equally spaced density contours from.5 to.5. op: second order (P ); middle: third order (P ); bottom: fourth order (P 3 ). Cells: 8.
25 5.5.5 Y Y Y.5.5 Figure 6: Double Mach refection problem. RKDG with HWENO limiter. Zoom-in pictures around the Mach stem. 3 equally spaced density contours from.5 to.5. Left: second order (P ); right: third order (P ); bottom: fourth order (P 3 ). Cells: 8.
26 6 Y.5 3 Y.5 3 Y.5 3 Figure 7: Double Mach refection problem. RKDG with the HWENO limiter. roubled cells. Squares denote cells which are identified as troubled cells subject to the HWENO limiting. op: second order (P ); middle: third order (P ); bottom: fourth order (P 3 ). Cells: 8.
27 7 6 Concluding remarks We have constructed a class of Hermite weighted essentially non-oscillatory (WENO) limiters, based on the procedure of [9], for the Runge-Kutta discontinuous Galerkin (RKDG) methods. he general framework of such HWENO limiters for RKDG methods, namely first identifying troubled cells subject to the HWENO limiting, then reconstructing the polynomial solution inside the troubled cell by the freedoms of the solutions of the DG method on the target cell and its adjacent neighboring cells by a HWENO procedure in a least square sense [8], is followed in this paper. he main novelty of this paper is the HWENO reconstruction procedure, which uses only the information from the troubled cell and its immediate neighbors (two cells in D and four cells in D), without any other extensive usage of geometric information of the meshes, and with simple positive linear weights in the reconstruction procedure. he methodology in this paper is more sophisticated than the original one specified in [9], but it leads to better resolutions for some examples for the P 3 case, such as those in Example 5.7 to Example 5., without the help of positivity-preserving limiters, in one and two dimensions. Numerical results are provided to demonstrate good results, both in accuracy and in non-oscillatory performance, comparable with those in earlier literature with more complicated WENO [, ] or HWENO [8, 9, ] limiters. In the future we will extend this limiter to unstructured meshes and to higher order polynomials (p-version DG methods). Acknowledgments he research was partially supported by NSFC grant 375, 93 and 38 and DOE grant DE-FG-8ER5863 and NSF grant DMS-875. References [] D. S. Balsara and C.-W. Shu, Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy, Journal of Computational Physics, 6 (), 5-5. [] R. Biswas, K.D. Devine and J. Flaherty, Parallel, adaptive finite element methods for conservation laws, Applied Numerical Mathematics, (99), [3] A. Burbeau, P. Sagaut and C.H. Bruneau, A problem-independent limiter for high-order Runge- Kutta discontinuous Galerkin methods, Journal of Computational Physics, 69 (), -5. [] B. Cockburn, S. Hou and C.-W. Shu, he Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws IV: the multidimensional case, Mathematics of Computation, 5 (99), [5] B. Cockburn, S.-Y. Lin and C.-W. Shu, VB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one dimensional systems, Journal of Computational Physics, 8 (989), 9-3. [6] B. Cockburn and C.-W. Shu, VB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws II: general framework, Mathematics of Computation, 5 (989), -35.
28 8 [7] B. Cockburn and C.-W. Shu, he Runge-Kutta discontinuous Galerkin method for conservation laws V: multidimensional systems, Journal of Computational Physics, (998), 99-. [8] M. Dumbser, D.S. Balsara, E.F. oro and C.D. Munz, A unified framework for the construction of one-step finite-volume and discontinuous Galerkin schemes on unstructured meshes, Journal of Computational Physics, 7 (8), [9] M. Dumbser and M. Käser, Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems, Journal of Computational Physics, (7), [] M. Dumbser, O. Zanotti, R. Loubère and S. Diot, A posteriori subcell limiting of the discontinuous Galerkin finite element method for hyperbolic conservation laws, Journal of Computational Physics, 78 (), [] O. Friedrichs, Weighted essentially non-oscillatory schemes for the interpolation of mean values on unstructured grids, Journal of Computational Physics, (998), 9-. [] C. Hu and C.-W. Shu, Weighted essentially non-oscillatory schemes on triangular meshes, Journal of Computational Physics, 5 (999), [3] G. Jiang and C.-W. Shu, Efficient implementation of weighted ENO schemes, Journal of Computational Physics, 6 (996), -8. [] V.P. Korobeinikov, Problems of Point-Blast heory, American Institute of Physics, 99. [5] L. Krivodonova, J. in, J.-F. Remacle, N. Chevaugeon and J.E. Flaherty, Shock detection and limiting with discontinuous Galerkin methods for hyperbolic conservation laws, Applied Numerical Mathematics, 8 (), [6]. Linde, P.L. Roe, Robust Euler codes, in: 3th Computational Fluid Dynamics Conference, AIAA Paper [7]. Liu, S. Osher and. Chan, Weighted essentially non-oscillatory schemes, Journal of Computational Physics, 5 (99), -. [8] H. Luo, J.D. Baum and R. Lohner, A Hermite WENO-based limiter for discontinuous Galerkin method on unstructured grids, Journal of Computational Physics, 5 (7), [9] J. Qiu and C.-W. Shu, Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method: one dimensional case, Journal of Computational Physics, 93 (3), [] J. Qiu and C.-W. Shu, Runge-Kutta discontinuous Galerkin method using WENO limiters, SIAM Journal on Scientific Computing, 6 (5), [] J. Qiu and C.-W. Shu, A comparison of troubled-cell indicators for Runge-Kutta discontinuous Galerkin methods using weighted essentially nonoscillatory limiters, SIAM Journal on Scientific Computing, 7 (5), [] J. Qiu and C.-W. Shu, Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method II: two dimensional case, Computers and Fluids, 3 (5), [3] L.I. Sedov, Similarity and Dimensional Methods in Mechanics, Academic Press, New York, 959. [] C.-W. Shu, Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, In Advanced Numerical Approximation of Nonlinear Hyperbolic Equations, B. Cockburn, C. Johnson, C.-W. Shu and E. admor (Editor: A. Quarteroni), Lecture Notes in Mathematics, volume 697, Springer, 998, [5] C.-W. Shu and S. Osher, Efficient implementation of essentially non-oscillatory shock-capturing schemes, Journal of Computational Physics, 77 (988), [6] C.-W. Shu and S. Osher, Efficient implementation of essentially non-oscillatory shock capturing
29 schemes II, Journal of Computational Physics, 83 (989), [7] C. Wang,. Zhang, C.-W. Shu and J. Ning, Robust high order discontinuous Galerkin schemes for two-dimensional gaseous detonations, Journal of Computational Physics, 3 (), [8] P. Woodward and P. Colella, he numerical simulation of two-dimensional fluid flow with strong shocks, Journal of Computational Physics, 5 (98), [9]. Zhong and C.-W. Shu, A simple weighted essentially nonoscillatory limiter for Runge-Kutta discontinuous Galerkin methods, Journal of Computational Physics, 3 (3), [3] J. Zhu, J. Qiu, C.-W. Shu and M. Dumbser, Runge-Kutta discontinuous Galerkin method using WENO limiters II: Unstructured meshes, Journal of Computational Physics, 7 (8), [3] J. Zhu,. Zhong, C.-W. Shu and J.. Qiu, Runge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured meshes, Journal of Computational Physics, 8 (3), -. 9
Runge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter
Runge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter Jun Zhu, inghui Zhong, Chi-Wang Shu 3 and Jianxian Qiu 4 Abstract In this paper, we propose a new type of weighted
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
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 informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
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 informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
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 informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationHigh order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes
High order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes Feng Zheng, Chi-Wang Shu and Jianxian Qiu 3 Abstract In this paper, a new type of high order Hermite
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More informationFinite difference Hermite WENO schemes for the Hamilton-Jacobi equations 1
Finite difference Hermite WENO schemes for the Hamilton-Jacobi equations Feng Zheng, Chi-Wang Shu 3 and Jianian Qiu 4 Abstract In this paper, a new type of finite difference Hermite weighted essentially
More informationEntropy stable high order discontinuous Galerkin methods. for hyperbolic conservation laws
Entropy stable high order discontinuous Galerkin methods for hyperbolic conservation laws Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Tianheng Chen, and with Yong Liu
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationA new type of Hermite WENO schemes for Hamilton-Jacobi equations
A new type of Hermite WENO schemes for Hamilton-Jacobi equations Jun Zhu and Jianxian Qiu Abstract In this paper, a new type of Hermite weighted essentially non-oscillatory (HWENO) schemes is designed
More informationA minimum entropy principle of high order schemes for gas dynamics. equations 1. Abstract
A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang and Chi-Wang Shu 3 Abstract The entropy solutions of the compressible Euler equations satisfy a minimum principle
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationDISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS
DISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS CHI-WANG SHU Abstract. In these lectures we give a general survey on discontinuous Galerkin methods for solving
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationNumerische Mathematik
Numer. Math. ) :545 563 DOI.7/s--443-7 Numerische Mathematik A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang Chi-Wang Shu Received: 7 July / Revised: 5 September
More informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More informationDiscontinuous Galerkin method for hyperbolic equations with singularities
Discontinuous Galerkin method for hyperbolic equations with singularities Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Qiang Zhang, Yang Yang and Dongming Wei Outline Introduction
More informationA central discontinuous Galerkin method for Hamilton-Jacobi equations
A central discontinuous Galerkin method for Hamilton-Jacobi equations Fengyan Li and Sergey Yakovlev Abstract In this paper, a central discontinuous Galerkin method is proposed to solve for the viscosity
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
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 informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationA Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws
A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University
More informationAn interface treating technique for compressible multi-medium flow with. Runge-Kutta discontinuous Galerkin method. Abstract
An interface treating technique for compressible multi-medium flow with Runge-Kutta discontinuous Galerkin method Chunwu Wang 1 and Chi-Wang Shu 2 Abstract The high-order accurate Runge-Kutta discontinuous
More informationALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS
ALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS HAILIANG LIU AND MICHAEL POLLACK Abstract. In this work, we propose a high resolution Alternating Evolution Discontinuous
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 informationTOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES
MATHEMATICS OF COMPUTATION Volume 67 Number 221 January 1998 Pages 73 85 S 0025-5718(98)00913-2 TOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES SIGAL GOTTLIEB AND CHI-WANG SHU Abstract. In this paper we
More informationHIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS
June 6, 7 :7 WSPC/Lecture Notes Series: 9in x 6in chapter HIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS Chi-Wang Shu Division of Applied Mathematics, Brown University Providence,
More informationA general well-balanced finite volume scheme for Euler equations with gravity
A general well-balanced finite volume scheme for Euler equations with gravity Jonas P. Berberich, Praveen Chandrashekar, Christian Klingenberg Abstract We present a second order well-balanced Godunov-type
More informationA High Order WENO Scheme for a Hierarchical Size-Structured Model. Abstract
A High Order WENO Scheme for a Hierarchical Size-Structured Model Jun Shen 1, Chi-Wang Shu 2 and Mengping Zhang 3 Abstract In this paper we develop a high order explicit finite difference weighted essentially
More informationMaximum-principle-satisfying and positivity-preserving high order schemes for conservation laws: Survey and new developments
Maximum-principle-satisfying and positivity-preserving high order schemes for conservation laws: Survey and new developments By Xiangxiong Zhang 1 and Chi-Wang Shu 1 Department of Mathematics, Brown University,
More informationFUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC PROBLEMS WITH DISCONTINUITIES
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 018, Glasgow, UK FUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
More informationReceived 6 August 2005; Accepted (in revised version) 22 September 2005
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol., No., pp. -34 Commun. Comput. Phys. February 6 A New Approach of High OrderWell-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a
More informationarxiv: v3 [math.na] 10 Mar 2015
A Posteriori Subcell Limiting of the Discontinuous Galerkin Finite Element Method for Hyperbolic Conservation Laws Michael Dumbser a, Olindo Zanotti a, Raphaël Loubère b, Steven Diot c a Department of
More informationThe Discontinuous Galerkin (dg) method is a popular numerical method that has received much attention
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2595 Discontinuous Galerkin method for multifluid Euler equations Marc
More informationSimulations of Shallow Water Equations with Finite Difference Lax-Wendroff Weighted Essentially Non-oscillatory Schemes
J Sci Comput (2011) 47: 281 302 DOI 10.1007/s10915-010-9437-3 Simulations of Shallow Water Equations with Finite Difference Lax-Wendroff Weighted Essentially Non-oscillatory Schemes Changna Lu Jianxian
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
More informationThe Runge Kutta Discontinuous Galerkin Method for Conservation Laws V
JOURNAL OF COMPUTATIONAL PHYSICS 141, 199 224 1998) ARTICLE NO. CP985892 The Runge Kutta Discontinuous Galerkin Method for Conservation Laws V Multidimensional Systems Bernardo Cockburn,1 and Chi-Wang
More informationTHE RUNGE KUTTA DISCONTINUOUS GALERKIN METHOD FOR CONSERVATION LAWS V: MULTIDIMENSIONAL SYSTEMS. Bernardo Cockburn 1. School of Mathematics
THE RUNGE KUTTA DISCONTINUOUS GALERKIN METHOD FOR CONSERVATION LAWS V: MULTIDIMENSIONAL SYSTEMS Bernardo Cockburn School of Mathematics University of Minnesota Minneapolis, Minnesota 55455 and Chi-Wang
More informationA genuinely high order total variation diminishing scheme for one-dimensional. scalar conservation laws 1. Abstract
A genuinely high order total variation diminishing scheme for one-dimensional scalar conservation laws Xiangxiong Zhang and Chi-Wang Shu 3 Abstract It is well known that finite difference or finite volume
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 informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationA Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
More information30 crete maximum principle, which all imply the bound-preserving property. But most
3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical
More informationADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws
ADER Schemes on Adaptive Triangular Meshes for Scalar Conservation Laws Martin Käser and Armin Iske Abstract. ADER schemes are recent finite volume methods for hyperbolic conservation laws, which can be
More informationSolution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers
Solution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers Alexander Kurganov, 1, * Eitan Tadmor 2 1 Department of Mathematics, University of Michigan, Ann Arbor, Michigan
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationStrong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
More informationAn efficient implementation of the divergence free constraint in a discontinuous Galerkin method for magnetohydrodynamics on unstructured meshes
An efficient implementation of the divergence free constraint in a discontinuous Galerkin method for magnetohydrodynamics on unstructured meshes Christian Klingenberg, Frank Pörner, Yinhua Xia Abstract
More informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More informationRESEARCH HIGHLIGHTS. WAF: Weighted Average Flux Method
RESEARCH HIGHLIGHTS (Last update: 3 rd April 2013) Here I briefly describe my contributions to research on numerical methods for hyperbolic balance laws that, in my view, have made an impact in the scientific
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationHermite WENO schemes and their application as limiters for Runge Kutta discontinuous Galerkin method: one-dimensional case
Journal of Computational Physics 93 (3) 5 35 www.elsevier.com/locate/jcp Hermite WENO schemes and their application as limiters for Runge Kutta discontinuous Galerkin method: one-dimensional case Jianian
More informationConvex ENO Schemes for Hamilton-Jacobi Equations
Convex ENO Schemes for Hamilton-Jacobi Equations Chi-Tien Lin Dedicated to our friend, Xu-Dong Liu, notre Xu-Dong. Abstract. In one dimension, viscosit solutions of Hamilton-Jacobi (HJ equations can be
More informationCompact Central WENO Schemes for Multidimensional Conservation Laws
arxiv:math/9989v [math.na] 3 Nov 999 Compact Central WENO Schemes for Multidimensional Conservation Laws Doron Levy Gabriella Puppo Giovanni Russo Abstract We present a new third-order central scheme for
More informationOn the Comparison of the Finite Volume and Discontinuous Galerkin Methods
Diploma Thesis Institute for Numerical Simulation, TUHH On the Comparison of the Finite Volume and Discontinuous Galerkin Methods Corrected version Katja Baumbach August 17, 2006 Supervisor: Prof. Dr.
More informationEssentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws
NASA/CR-97-0653 ICASE Report No. 97-65 Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws Chi-Wang Shu Brown University Institute for Computer
More informationSMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES
Proceedings of ALGORITMY 2016 pp. 155 164 SMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES TAMER H. M. A. KASEM AND FRANÇOIS G. SCHMITT Abstract. The weighted essentially non-oscillatory
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationH. L. Atkins* NASA Langley Research Center. Hampton, VA either limiters or added dissipation when applied to
Local Analysis of Shock Capturing Using Discontinuous Galerkin Methodology H. L. Atkins* NASA Langley Research Center Hampton, A 68- Abstract The compact form of the discontinuous Galerkin method allows
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 informationA Central Discontinuous Galerkin Method for Hamilton-Jacobi Equations
J Sci Comput (00) 45: 404 48 DOI 0.007/s095-009-9340-y A Central Discontinuous Galerkin Method for Hamilton-Jacobi Equations Fengyan Li Sergey Yakovlev Received: 6 August 009 / Revised: November 009 /
More informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
More informationA Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang 3, and Hongkai Zhao 4 ABSTRACT In this paper, we construct a second order fast
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 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 resolution of discontinuous Galerkin methods for time dependent. wave equations 1. Abstract
Numerical resolution of discontinuous Galerkin methods for time dependent wave equations Xinghui Zhong 2 and Chi-Wang Shu Abstract The discontinuous Galerkin DG method is known to provide good wave resolution
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 informationAccepted Manuscript. A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
Accepted Manuscript A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang, Hongkai Zhao PII: S-()- DOI:./j.jcp... Reference: YJCPH To
More informationEmerging Numerical Methods for Atmospheric Modeling
Chapter 9 Emerging Numerical Methods for Atmospheric Modeling Ramachandran D. Nair, Michael N. Levy and Peter H. Lauritzen Abstract This chapter discusses the development of discontinuous Galerkin (DG)
More informationHigh order finite difference WENO schemes with positivity-preserving. limiter for correlated random walk with density-dependent turning rates
High order finite difference WENO schemes with positivity-preserving limiter for correlated random walk with density-dependent turning rates Yan Jiang, Chi-Wang Shu and Mengping Zhang 3 Abstract In this
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationA simple robust and accurate a posteriori subcell finite
A simple robust and accurate a posteriori subcell finite volume limiter for the discontinuous Galerkin method Michael Dumbser ERC Grant STiMulUs FET-HPC Grant ExaHyPE 1/ 69 Séminaire LJLL, Paris, France,
More informationEfficient Implementation of Weighted ENO Schemes*
JOURNAL OF COMPUTATIONAL PHYSICS 126, 202 228 (1996) ARTICLE NO. 0130 Efficient Implementation of Weighted ENO Schemes* GUANG-SHAN JIANG AND CHI-WANG SHU Division of Applied Mathematics, Brown University,
More informationA weighted essentially non-oscillatory numerical scheme for a multi-class traffic flow model on an inhomogeneous highway
Manuscript A weighted essentially non-oscillatory numerical scheme for a multi-class traffic flow model on an inhomogeneous highway Peng Zhang, S.C. Wong and Chi-Wang Shu. Department of Civil Engineering,
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 informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationA positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations
A positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations Sashank Srinivasan a, Jonathan Poggie a, Xiangxiong Zhang b, a School of Aeronautics and Astronautics,
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More 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 informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
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 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More information