A class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
|
|
- Sabrina Watts
- 5 years ago
- Views:
Transcription
1 Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, math.scichina.com A class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes ZHU Jun 1, & QIU JianXian 1 1 Department of Mathematics, Nanjing University, Nanjing 10093, China College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing 10016, China ( zhujun@nuaa.edu.cn, jxqiu@nju.edu.cn) Abstract In this paper, we developed a class of the fourth order accurate finite volume Hermite weighted essentially non-oscillatory (HWENO) schemes based on the work (Computers & Fluids, 34: (005)) by Qiu and Shu, with Total Variation Diminishing Runge-Kutta time discretization method for the two-dimensional hyperbolic conservation laws. The key idea of HWENO is to evolve both with the solution and its derivative, which allows for using Hermite interpolation in the reconstruction phase, resulting in a more compact stencil at the expense of the additional work. The main difference between this work and the formal one is the procedure to reconstruct the derivative terms. Comparing with the original HWENO schemes of Qiu and Shu, one major advantage of new HWENO schemes is its robust in computation of problem with strong shocks. Extensive numerical experiments are performed to illustrate the capability of the method. Keywords: finite volume HWENO scheme, conservation laws, Hermite polynomial, TVD Runge- Kutta time discretization method MSC(000): 65M06, 65M99, 35L65 1 Introduction In recent decades, numerical schemes for solving hyperbolic conservation laws have been developed rapidly. These schemes improve the first-order methods of Godunov [1] to arbitrary order of accuracy. In order to achieve uniform high order accuracy, Harten and Osher [] gave a weaker version of the Total Variation Diminishing (TVD) [3] criterion, and on which they established the framework for the reconstruction of high order essentially non-oscillatory (ENO) type schemes. Then Harten, et al. [4] developed ENO schemes to solve one dimensional problems. The key idea of ENO schemes is applying the most smooth stencil among all candidate stencils to approximate the variables at cell boundaries to a high order of accuracy and avoid oscillations near discontinuities. In [5], Harten extended the finite volume ENO scheme to the two-dimensional hyperbolic conservation laws. Following the work of Harten, Casper [6] andwithatkins [7] considered the finite volume approach in developing multi-dimensional, high order accurate ENO schemes, and Abgrall [8] extended the method on unstructured meshes. Received March 1, 008; accepted May 8, 008 DOI: /s Corresponding author This work was partially supported by the National Natural Science Foundation of China (Grant No ), the European project ADIGMA on the development of innovative solution algorithms for aerodynamic simulations, Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry and the Natural Science Foundation of Jiangsu Province (Grant No. BK006511)
2 1550 ZHU Jun & QIU JianXian In 1994, Liu, Osher and Chan [9] proposed a Weighted ENO (WENO) scheme that was constructed from the r-th order ENO schemes to obtain (r+1)-th order accuracy. In 1996, Jiang and Shu [10] proposed the framework to construct finite difference WENO schemes from the r-th order (in L 1 norm sense) ENO schemes to get (r 1)-th order accuracy, gave a new way of measuring the smoothness indicators, and emulated the ideas of minimizing the total variation of the approximation. The construction of finite volume WENO methods on unstructured meshes was presented by Friedrichs [11]. Instead of construction of two-dimensional finite volume WENO schemes used dimensional by dimensional methods, Hu and Shu [1] proposed a full dimensional reconstruction methodology for the third order WENO schemes by using a combination of two-dimensional linear polynomials and the fourth order WENO schemes by using a combination of two dimensional quadratic polynomials, and presented a new way of measuring two dimensional smoothness of numerical solutions. In [13, 14], Qiu and Shu developed the methodology, termed as Hermite WENO (HWENO) schemes and applied them as limiters to Runge-Kutta discontinuous Galerkin methods. The key idea of HWENO is to evolve both with the solution and its derivative. This allows for using Hermite interpolation in the reconstruction phase, resulting in a more compact stencil at the expense of additional work. In this paper, following the ideas of [14], we construct a kind of two-dimensional finite volume HWENO schemes which has fourth order accuracy, the main difference between this work and [14] is the procedure to reconstruct the derivative terms. Comparing with the original HWENO schemes of Qiu and Shu [14], one major advantage of new HWENO schemes is its robust in the computation of problem with strong shocks. The organization of this paper is as follows: In Section we review and construct two-dimensional finite volume HWENO scheme in details and present extensive numerical results in Section 3 to verify the accuracy and stability of this approach. Concluding remarks are given in Section 4. Reconstruction of HWENO scheme In this section we consider two-dimensional conservation laws (1): u t + f(u) x + g(u) y =0, u(x, y, 0) = u 0 (x, y). (1) For simplicity of presentation, we assume that the mesh is uniform with the cell size x i+1/ x i 1/ =Δx, y j+1/ y j 1/ =Δy and cell centers (x i,y j )=( 1 (x i+1/ + x i 1/ ), 1 (y j+1/ + y j 1/ )). We also denote the cells by I ij =[x i 1/,x i+1/ ] [y j 1/,y j+1/ ]. Let v = u x, w = u y. Taking the derivatives of (1), we obtain v t + h x + r y =0, v(x, y, 0) = u 0(x, y) (), x w t + q x + s y =0, w(x, y, 0) = u 0(x, y) (3), y where h(u, v) =f (u) v, r(u, v) =g (u) v, q(u, w) =f (u) w, s(u, w) =g (u) w. We integrate the equations (1), () and (3) on a control volume I ij to obtain the semi-discrete
3 A class of fourth order finite volume hermite weighted essentially non-oscillatory schemes 1551 finite volume scheme as: d dt u ij = 1 F nds, (4) ΔxΔy I ij d dt v ij = 1 H nds, (5) Δy I ij d dt w ij = 1 Q nds, (6) Δx I ij where u ij = 1 udxdy, v ij = 1 u ΔxΔy I ij Δy I ij x dxdy, w ij = 1 u Δx I ij y dxdy and F =(f,g) T,H =(h, r) T,Q=(q, s) T. The line integrals in (4) (6) are discretized by a q-point Gaussian integration formula q F nds I ij ω l F (u(g l,t)) n, (7) I ij l=1 q H nds I ij ω l H(u(G l,t),v(g l,t)) n, (8) I ij l=1 q Q nds I ij ω l Q(u(G l,t),w(g l,t)) n. (9) I ij l=1 Since we are constructing schemes up to fourth-order accuracy, two-point Gaussian will be used in each line integration, and F (u(g l,t)) n, H(u(G l,t),v(g l,t)) n, Q(u(G l,t),w(g l,t)) n are replaced by numerical fluxes such as the Lax-Friedrichs fluxes: f(u(g l,t)) 1 [f(u (G l,t)) + f(u + (G l,t)) α(u + (G l,t) u (G l,t))], (10) h(u(g l,t),v(g l,t)) 1 [h(u (G l,t),v (G l,t)) + h(u + (G l,t),v + (G l,t)) α(v + (G l,t) v (G l,t))], (11) q(u(g l,t),w(g l,t)) 1 [q(u (G l,t),w (G l,t)) + q(u + (G l,t),w + (G l,t)) α(w + (G l,t) w (G l,t))], (1) for G l =(x i±1/,y j± 3/6 ), u ± (G l,t), v ± (G l,t), w ± (G l,t) are the left and right limits of the solutions u, v, w at the cell interface G l respectively; and g(u(g l,t)) 1 [g(u (G l,t)) + g(u + (G l,t)) α(u + (G l,t) u (G l,t))], (13) r(u(g l,t),v(g l,t)) 1 [r(u (G l,t),v (G l,t)) + r(u + (G l,t),v + (G l,t)) α(v + (G l,t) v (G l,t))], (14) s(u(g l,t),w(g l,t)) 1 [s(u (G l,t),w (G l,t)) + s(u + (G l,t),w + (G l,t)) α(w + (G l,t) w (G l,t))], (15) for G l =(x i± 3/6,y j±1/ ), u ± (G l,t), v ± (G l,t), w ± (G l,t) are the bottom and top limits of the solutions u, v, w at the cell interface G l respectively. What we want to do is to reconstruct
4 155 ZHU Jun & QIU JianXian u ± (G l,t), v ± (G l,t), w ± (G l,t)from{u ij, v ij, w ij }. The procedure of reconstruction of u ± (G l,t) is the same as that in [14], we do not repeat it here for saving space. Reconstruction of v (G l,t)atg l =(x i+ 1,y j+ 1 ). For simplicity we relabel the cell I ij and 3 its neighboring cells as I 1,...,I 9 as shown in Figure 1, where I ij is relabeled as I 5. Figure 1 The big stencil First, we choose the following stencils: S 1 = {I 1,I,I 4,I 5 }, S = {I,I 3,I 5,I 6 }, S 3 = {I 4,I 5,I 7,I 8 }, S 4 = {I 5,I 6,I 8,I 9 }, S 5 = {I 1,I,I 3,I 4,I 5,I 7 }, S 6 = {I 1,I,I 3,I 5,I 6,I 9 }, S 7 = {I 1,I 4,I 5,I 7,I 8,I 9 }, S 8 = {I 3,I 5,I 6,I 7,I 8,I 9 }. We construct the polynomials p n (x, y)(n =1,...,4) in the small stencils to approximate function u such that 1 1 p n (x, y) p n (x, y)dxdy = u k, dxdy = v kx, ΔxΔy I k Δy I kx x 1 p n (x, y) dxdy = w ky, Δx I ky y for n =1,k=1,, 4, 5, k x =1, 4, 5, k y =1,, 5; n =,k=, 3, 5, 6, k x =3, 5, 6, k y =, 3, 5; n =3,k=4, 5, 7, 8, k x =4, 5, 7, k y =5, 7, 8; n =4,k=5, 6, 8, 9, k x =5, 6, 9, k y =5, 8, 9; and p n (x, y)(n =5,...,8) in the small stencils to approximate v, the first derivative of u, such 1 that Δy I kx p n (x, y)dxdy = v kx for n =5,k x =1,, 3, 4, 5, 7; n =6,k x =1,, 3, 5, 6, 9; n =7,k x =1, 4, 5, 7, 8, 9; n =8,k x =3, 5, 6, 7, 8, 9. For simplicity, we rewrite u, v,w as u,v,w in formulae (16) (3) and (6) (33). At point G l =(x i+ 1,y j+ 1 )wehave 3 x p 1(x i+ 1,y j+ 1 )= 7 3 v v 4 4u 5 +4u u1 3u 3u4 + 3u v w1 1 3w 1 3v4, (16)
5 A class of fourth order finite volume hermite weighted essentially non-oscillatory schemes 1553 x p (x i+ 1,y j+ 1 )= 1 3u6 1 3w u3 u 5 +u 6 1 v 5 1 v u 3 1 3u w, (17) 3 6 x p 3(x i+ 1,y j+ 1 )= 4u 5 +4u v v u4 3u v u u8 1 3v w7 1 3w8, (18) x p 4(x i+ 1,y j+ 1 )= 1 3u6 u 5 +u v 5 1 v u9 1 3w u u w8, (19) 3 6 p 5 (x i+ 1,y j+ 1 )= v 1 3 v v 3 1 v v v1 1 3v 1 3v v v7, (0) 1 p 6 (x i+ 1,y j+ 1 )= v v 1 6 v v v 6 1 3v + 1 3v5 1 3v v9, (1) 1 p 7 (x i+ 1,y j+ 1 )= 1 3 v v v 7 3 v v v v4 1 3v v v8, () 6 4 p 8 (x i+ 1,y j+ 1 )= 1 3 v v v v v v3 1 3v v v8. (3) 1 In order to obtain the 4th order approximation to function u, the approximation to the derivatives of u can be only the 3rd order. From definition of approximation polynomials of p n (x, y), it is easy to see that the x p n(x i+ 1,y j+ 1 )(n =1,...,4) and p n (x 3 i+ 1,y j+ 1 )(n = 3 5,...,8) approximates to v = u x at point G l up to the 3rd order when u is smooth enough. In this paper, we just choose the linear weights to be γ l = 1 8 (l =1,...,8). As in [1,15], we use the following smoothness indicator: β l = β l = ( I 5 α 1 D α ) I 5 x p l(x, y) dxdy, l =1,...,4, α =1 α =1 I 5 I 5 α 1 (D α p l (x, y)) dxdy, l =5,...,8. In the actual numerical implementation the smoothness indicator β l (l =1,...,8) can be written explicitly as follows: β 1 = 16 3 w w 1w w + 04u 4 v 5 +18u 1 36u 1 u 36u 1 u 4 +36u 1 u u 1v u 1v 4 +18u +36u u 4 36u u u v u v u 4 40u 4 u u 4v u 4v u u 5v u 5v 4 04u 5 v 5 +55v 5
6 1554 ZHU Jun & QIU JianXian +94v 5 v 4 +v 1 w 1 v 1 w w 1 v 4 +w v u 5w u 5w v v 1v v u 1w u 1w 44 3 u w u w 44 3 u 4w u 4w, (6) β = 10u u 6v w 40u 5 u u 5v u 5v 3 +18u 36u u u u 5 v 5 +55v u 5w u w 36u u 3 +36u u u v u v 6 +18u 3 +36u 3 u 5 36u 3 u u 3v u 3v u 6v v v 3v v u w u 3w u 3w u 5w u 6w u 6w w w 3w 04u 6 v 5 +94v 5 v 6 v 3 w 3 +v 3 w +w 3 v 6 w v 6, (7) β 3 = 10u 4 40u 4 u 5 36u 4 u 7 +36u 4 u u 4v u 4v 4 +36u 5 u 7 36u 5 u u 5v u 5v 4 +18u 7 36u 7 u u 7v u 7v 4 +18u u 8v u 8v v v 7v v u 4w u 4w u 5w u 5w u 7w u 7w u 8w u 8w w 7 v 7 w 7 +v 7 w w 7w u 4 v w 8 v 4 w u 5 04u 5 v 5 +55v 5 +94v 5 v 4 +w 7 v 4, (8) β 4 = 36u 5 u 8 +36u 5 u u 5v u 5v u 6 +36u 6 u 8 36u 6 u u 6v u 6v 6 +18u 8 36u 8 u u 8v u 8v 6 +18u u 9v u 9v v v 9v u 5 40u 5 u u 5 v 5 04u 6 v 5 +55v 5 +94v 5 v 6 +v 9 w 9 v 9 w 8 +v 6 w 8 w 9 v v u 5w u 5w u 6w u 6w u 8w u 8w u 9w u 9w w w 9w w 8, (9) β 5 = 4 3 v 7 + v 5 v v 4v v 4v v 1v 6 3 v 1v v 1v v v v v v 1v v 1v v v v v v v v v v 7 v 3 v 4 + v 3 v 5, (30) β 6 = 16 3 v 1v + v 1 v v v v 1v v v v 3 v 1 v 6 v v 9 + v 5 v v v v v 3v v v v 3v v 3v v 5v v v 6v v 9, (31) β 7 = 3 6 v v 5v v 4v v 4v v 1v 4 + v 1 v v 1v 7 + v 5 v v v v 1 v 4 v 9 v 1 v v 4v v 5v v 7v v v 7v v 8v v 9, (3)
7 A class of fourth order finite volume hermite weighted essentially non-oscillatory schemes 1555 β 8 = v 5 v v 5v v 5v v 5v v v 6v v 6v v 8v v v v v 3v v 3v v 7v v 7v v 7 + v 3 v 5 v 3 v 8 v 6 v v 5. (33) Then we compute the nonlinear weight based on the linear weight and smoothness indicator [15] : ω l = v i+ 1,j+ 1 3 ω l 8 k=1 ω, ω l = k l=1 γ l 8 k=1 (ε + β, l =1,...,8, (34) k) where γ l is the linear weight determined in the above step, and ε is a small positive number to avoid the denominator to become 0. We use ε =10 6 in all the computations in this paper. The final approximation is given by 4 8 ω l x p l(x i+ 1,y j+ 1 )+ ω l p l (x 3 i+ 1,y j+ 1 ). 3 Remarks. The main difference between the procedure of reconstruction of v and that in [14] is that in this work, more derivative terms are used in the reconstruction, and full -degree polynomials are constructed for the derivatives in the small stencils. The reconstructions to v + (G l )andw ± (G l ) are similar to the above procedure. For systems of conservation laws, such as the Euler equations of gas dynamics, all of the reconstructions are performed in the local characteristic directions to avoid oscillations. The semidiscrete scheme (4) (6), written as u t = L(u) is then discretized in time by a TVD Runge-Kutta method [16], for example the third order version given by u (1) = u n +ΔtL(u n ), u () = 3 4 un u(1) ΔtL(u(1) ), (35) u n+1 = 1 3 un + 3 u() + 3 ΔtL(u() ). 3 Numerical tests In this section we present the results of numerical tests for the scheme in the previous section. Example 3.1. We solve the following nonlinear scalar Burgers equation in two dimensions: ( ) ( ) u u u t + + = 0 (36) x y with the initial condition u(x, y, 0) = 0.5+ sin(π(x + y)/) and periodic boundary conditions in two directions. We compute the solution up to t =0.5/π. When t =0.5/π, the solution is still smooth. The errors and numerical orders of accuracy by the HWENO scheme are shown in Table 1. We can see that the designed order is obtained. Table 1 u t +( u )x +(u)y =0. u(x, y, 0) = 0.5+sin(π(x + y)/). Periodic boundary conditions in two directions. T =0.5/π. L 1 and L errors cells L 1 error L 1 order L error L order e 3.79e e e e e e e e e l=5
8 1556 ZHU Jun & QIU JianXian Example 3.. We solve the Euler equations ρ ρu ρv ρu + ρu + p + ρuv =0. (37) t ρv x ρuv y ρv + p E u(e + p) v(e + p) In which ρ is density, u is x-direction velocity, v is y-direction velocity, E is total energy, p is pressure. The initial conditions are: ρ(x, y, 0) = sin(π(x + y)), u(x, y, 0) = 0.7, v(x, y, 0) = 0.3, p(x, y, 0) = 1, periodic boundary conditions in two directions. We compute the density solution up to t =.0. The exact solution is ρ(x, y, t) = sin(π(x + y t)). In Table, we show the errors and numerical orders of accuracy by the HWENO scheme and we also can see that the designed order is obtained. Table D Euler equations: initial data ρ(x, y, 0) = sin(π(x + y)), u(x, y, 0) = 0.7, v(x, y, 0) = 0.3, p(x, y, 0) = 1. Periodic boundary conditions in two directions. T =.0. L 1 and L errors. cells L 1 error L 1 order L error L order e.13e e e e e e e e e Example 3.3. We solve the same nonlinear Burgers equation (36) for Riemann problem in a computational domain of [ 1, 3] [ 1, 3] and set the initial conditions as: 0.1, x>0,y >0, (1) u(x, y, 0) =.5, x<0,y >0, (38) 1.1, x<0,y <0, 1.5, x>0,y <0, 1.1, x>0,y >0, () u(x, y, 0) = 3.1, x<0,y >0, (39).1, x<0,y <0, 0.1, x>0,y <0. Figure Burgers equation for Riemann problem. T = cells. Left: initial condition (1), 15 equally spaced contours from 0.49 to.351; right: initial condition (), 15 equally spaced contours from 0.86 to.913.
9 A class of fourth order finite volume hermite weighted essentially non-oscillatory schemes 1557 We plot the results at t =0.8. The solutions are shown in Figure. We can see that the scheme gives non-oscillatory shock transitions for this problem. Example 3.4. D Euler equations for Riemann problem [17,18]. We solve the Euler equations (37) in a computational domain of [0, 1] [0, 1] and set the initial conditions as: (1) (0.5313, 0, 0, 0.4) T, x>0.5,y>0.5, (ρ, u, v, p) T (1, 0.776, 0, 1) T, x<0.5,y>0.5, = (0.8, 0, 0, 1) T, x<0.5,y<0.5, (1, 0, 0.776, 1) T, x>0.5,y<0.5, (40) () (1.1, 0, 0, 1.1) T, x>0.5,y >0.5, (ρ, u, v, p) T (0.5065, , 0, 0.35) T, x<0.5,y >0.5, = (1.1, , , 1.1) T, x<0.5,y <0.5, (0.5065, 0, , 0.35) T, x>0.5,y <0.5, (41) (3) (1, 0.1, 0, 1) T, x>0.5,y >0.5, (ρ, u, v, p) T (0.5313, 0.876, 0, 0.4) T, x<0.5,y >0.5, = (0.8, 0.1, 0, 0.4) T, x<0.5,y <0.5, (0.5313, 0.1, 0.776, 0.4) T, x>0.5,y <0.5, (4) (4) (0.5313, 0.1, 0.1, 0.4) T, x>0.5,y >0.5, (ρ, u, v, p) T (1.0, , 0.1, 1) T, x<0.5,y >0.5, = (0.8, 0.1, 0.1, 1) T, x<0.5,y <0.5, (1, 0.1, 0.876, 1) T, x>0.5,y <0.5, (43) (5) (1, 0.75, 0.5, 1) T, x>0.5,y>0.5, (ρ, u, v, p) T (, 0.75, 0.5, 1) T, x<0.5,y>0.5, = (1, 0.75, 0.5, 1) T, x<0.5,y<0.5, (3, 0.75, 0.5, 1) T, x>0.5,y<0.5. (44) In Figure 3, we show the computational results for density at (1) t =0.5, () t =0.5, (3) t =0.3, (4) t =0., (5) t =0.3, respectively. Example 3.5. Double mach reflection problem. We solve the Euler equations (37) in a computational domain of [0, 4] [0, 1]. A reflection wall lies at the bottom of the domain starting from x = 1 6, y=0,makinga60 angle with the x-axis. The reflection boundary condition is used at the wall, which for the rest of the bottom boundary (the part from x =0 to x = 1 6 ), the exact post-shock condition is imposed. At the top boundary is the exact motion of the mach 10 shock. The results shown are at t =0.. We present both the pictures of region [0, 3] [0, 1] and the blow-up region around the double mach stems in Figures 4 and 5 respectively. All pictures are the density contours with 30 equal spaced contour lines from 1.5 to.7. Remark. The HWENO schemes developed in [14] failed in this test case.
10 1558 ZHU Jun & QIU JianXian Figure 3 D Euler equations for Riemann problem with cells. From left to right and top to bottom. Initial condition (1) T = equally spaced density contours from 0.54 to 1.70; Initial condition () T = equally spaced density contours from 0.5 to 1.9; Initial condition (3) T = equally spaced density from contours 0.55 to 1.1; Initial condition (4) T = equally spaced density contours from 0.53 to 0.99; Initial condition (5) T = equally spaced density contours from 0.5 to Figure 4 Double Mach refection problem, 30 equally spaced density contours from 1.5 to.7. From top to bottom: cells, cells, cells. Example 3.6. A mach 3 wind tunnel with a step. The setup of the problem is as follows: The wind tunnel is 1 length unit wide and 3 length units long. The step is 0. length units high and is located 0.6 length units from the left end of the tunnel. Initially, a right going mach 3 flow is used. Reflective boundary conditions are applied alone the walls of the tunnel and in flow and out flow boundary conditions are applied at the entrance and the exit. The results are shown at t = 4. We present the pictures of whole region [0, 3] [0, 1] in Figure 6. All pictures are the density contours with 30 equal spaced contour lines from 0.3 to Remark. The HWENO schemes developed in [14] failed in this test case.
11 A class of fourth order finite volume hermite weighted essentially non-oscillatory schemes 1559 Figure 5 Double Mach refection problem zoom in, 30 equally spaced density contours from 1.5 to.7. From left to right and top to bottom: cells, cells, cells. Figure 6 A mach 3 wind tunnel with a step problem, 30 equally spaced density contours from 0.3 to From top to bottom: cells, cells and cells. 4 Concluding remarks In this paper, we constructed a class of fourth order accurate finite volume HWENO schemes based on the procedure of [14], for two dimensional hyperbolic conservation laws. The key
12 1560 ZHU Jun & QIU JianXian idea of HWENO is to evolve both with the solution and its derivative, this allows for using Hermite interpolation in the reconstruction phase, resulting in a more compact stencil at the expense of additional work. We use the procedure presented in [14] to reconstruct the values of solution, and the main difference between the procedure of reconstruction of derivatives in this paper and that in [14] is that more derivative terms are used in the reconstruction, and full -degree polynomials are constructed for the derivatives in the small stencils in this paper. Comparing with the original HWENO schemes in [14], one major advantage of new HWENO schemes is its robust in computation of problem with strong shocks, such as the double Mach reflection problem and forward step problem. Numerical experiments for two dimensional Burgers equation and Euler equations of compressible gas dynamics are presented to show the effectiveness of the methods. References 1 Godunov S K. A finite-difference method for the numerical computation of discontinuous solutions of the equations of fluid dynamics. Mat Sb, 47(3): (1959) Harten A, Osher S. Uniformly high-order accurate non-oscillatory schemes. IMRC Technical Summary Rept. 83, Univ. of Wisconsin, Madison, WI, May, Harten A. High resolution schemes for hyperbolic conservation laws. J Comput Phys, 49(3): (1983) 4 Harten A, Engquist B, Osher S, Chakravarthy S. Uniformly high order accurate essentially non-oscillatory schemes III. J Comput Phys, 71(): (1987) 5 Harten A. Preliminary results on the extension of ENO schemes to two-dimensional problems. In: C. Carasso et al. eds. Proceedings of International Conference on Nonlinear Hyperbolic Problems, Saint- Etienne, 1986, Lecture Notes in Mathematics, Berlin: Springer-Verlag, Casper J. Finite-volume implementation of high-order essentially nonoscillatory schemes in two dimensions. AIAA Journal, 30(1): (199) 7 Casper J, Atkins H L. A finite-volume high-order ENO scheme for two-dimensional hyperbolic systems. J Comput Phys, 106: 6 76 (1993) 8 Abgrall R. On essentially non-oscillatory schemes on unstructured meshes: Analysis and implementation. J Comput Phys, 114: (1994) 9 Liu X D, Osher S, Chan T. Weighted essentially non-oscillatory schemes. J Comput Phys, 115: 00 1 (1994) 10 Jiang G S, Shu C W. Efficient implementation of weighted ENO schemes. J Comput Phys, 16: 0 8 (1996) 11 Friedrichs O. Weighted essentially non-oscillatory schemes for the interpolation of mean values on unstructured grids. J Comput Phys, 144: (1998) 1 Hu C Q, Shu C W. Weighted essentially non-oscillatory schemes on triangular meshes. J Comput Phys, 150: (1999) 13 Qiu J, Shu C W. Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method: one-dimensional case. J Comput Phys, 193: (004) 14 Qiu J, Shu C W. Hermite WENO schemes and their application as limiters for Runge-Kutta discontiuous Galerkin method II: two dimensional case. Computers & Fluids, 34: (005) 15 Shu C W. Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws. In: Advanced Numerical Approximation of Nonlinear Hyperbolic Equations. A. Lecture Notes in Mathematics, CIME subseries Berlin-New York: Springer-Verlag, ICASE Report Shu C W, Osher S. Efficient implementation of essentially non-oscillatory shock capturing schemes. J Comput Phys, 77: (1988) 17 Lax P D, Liu X D. Solution of two dimensional Riemann problems of gas dynamics by positive schemes. SIAM J Sci Comput, 19(): (1998) 18 Brio M, Zakharian A R, Webb G M. Two dimensional Riemann solver for Euler equations of gas dynamics. J Comput Phys, 167: (001)
Weighted 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 informationRunge-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 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 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 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 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 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 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 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 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 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 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 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 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 informationRunge-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, College of Science, Nanjing University of Aeronautics and
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 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 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 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 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 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 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 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 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 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 informationARTICLE IN PRESS Mathematical and Computer Modelling ( )
Mathematical and Computer Modelling Contents lists available at ScienceDirect Mathematical and Computer Modelling ournal homepage: wwwelseviercom/locate/mcm Total variation diminishing nonstandard finite
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationAdaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems
EPJ Web of Conferences 108, 0019 (016) DOI: 10.1051/ epjconf/ 0161080019 C Owned by the authors, published by EDP Sciences, 016 Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear
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 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 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 informationHFVS: An Arbitrary High Order Flux Vector Splitting Method
HFVS: An Arbitrary High Order Flu Vector Splitting Method Yibing Chen, Song Jiang and Na Liu Institute of Applied Physics and Computational Mathematics, P.O. Bo 8009, Beijing 00088, P.R. China E-mail:
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 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 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 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 informationLecture 16. Wednesday, May 25, φ t + V φ = 0 (1) In upwinding, we choose the discretization of the spatial derivative based on the sign of u.
Lecture 16 Wednesday, May 25, 2005 Consider the linear advection equation φ t + V φ = 0 (1) or in 1D φ t + uφ x = 0. In upwinding, we choose the discretization of the spatial derivative based on the sign
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 informationTHIRD ORDER NON-OSCILLATORY CENTRAL SCHEME FOR MULTIDIMENSIONAL HYPERBOLIC CONSERVATION LAWS
Palestine Journal of Mathematics Vol. 5(Special Issue: ) (6), 8 37 Palestine Polytechnic University-PPU 6 THIRD ORDER NON-OSCILLATORY CENTRAL SCHEME FOR MULTIDIMENSIONAL HYPERBOLIC CONSERVATION LAWS A.
More informationConjugate lter approach for shock capturing
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2003; 19:99 110 (DOI: 10.1002/cnm.573) Conjugate lter approach for shock capturing Yun Gu 1 1; 2; ; and G. W. Wei 1 Department
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 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 FOURTH-ORDER CENTRAL WENO SCHEME FOR MULTIDIMENSIONAL HYPERBOLIC SYSTEMS OF CONSERVATION LAWS
SIAM J. SCI. COMPUT. Vol. 4, No., pp. 48 56 c Society for Industrial and Applied Mathematics A FOURTH-ORDER CENTRAL WENO SCHEME FOR MULTIDIMENSIONAL HYPERBOLIC SYSTEMS OF CONSERVATION LAWS DORON LEVY,
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 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 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 informationA method for avoiding the acoustic time step restriction in compressible flow
A method for avoiding the acoustic time step restriction in compressible flow Nipun Kwatra Jonathan Su Jón T. Grétarsson Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationKrylov single-step implicit integration factor WENO methods for advection-diffusion-reaction equations
Accepted Manuscript Krylov single-step implicit integration factor WENO methods for advection diffusion reaction equations Tian Jiang, Yong-Tao Zhang PII: S0021-9991(16)00029-2 DOI: http://dx.doi.org/10.1016/j.jcp.2016.01.021
More informationApplied Mathematics and Mechanics (English Edition) Transversal effects of high order numerical schemes for compressible fluid flows
Appl. Math. Mech. -Engl. Ed., 403), 343 354 09) Applied Mathematics and Mechanics English Edition) https://doi.org/0.007/s0483-09-444-6 Transversal effects of high order numerical schemes for compressible
More informationInvestigation of Godunov Flux Against Lax Friedrichs' Flux for the RKDG Methods on the Scalar Nonlinear Conservation Laws Using Smoothness Indicator
American Review of Mathematics and Statistics December 2014, Vol. 2, No. 2, pp. 43-53 ISSN: 2374-2348 (Print), 2374-2356 (Online) Copyright The Author(s). 2014. All Rights Reserved. Published by American
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 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 informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationA New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods
A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods Raymond J. Spiteri Steven J. Ruuth Technical Report CS-- May 6, Faculty of Computer Science 65 University Ave.,
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 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 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 informationArt Checklist. Journal Code: Article No: 6959
Art Checklist Journal Code: JCPH Article No: 6959 Disk Recd Disk Disk Usable Art # Y/N Format Y/N Remarks Fig. 1 Y PS Y Fig. 2 Y PS Y Fig. 3 Y PS Y Fig. 4 Y PS Y Fig. 5 Y PS Y Fig. 6 Y PS Y Fig. 7 Y PS
More informationA re-averaged WENO reconstruction and a third order CWENO scheme for hyperbolic conservation laws
A re-averaged WENO reconstruction and a third order CWENO scheme for hyperbolic conservation laws Chieh-Sen Huang a,, Todd Arbogast b,, Chen-Hui Hung c,3 a Department of Applied Mathematics and National
More informationOrder of Convergence of Second Order Schemes Based on the Minmod Limiter
Order of Convergence of Second Order Schemes Based on the Minmod Limiter Boan Popov and Ognian Trifonov July 5, 005 Abstract Many second order accurate non-oscillatory schemes are based on the Minmod limiter,
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 informationWell-Balanced Schemes for the Euler Equations with Gravity
Well-Balanced Schemes for the Euler Equations with Gravity Alina Chertock North Carolina State University chertock@math.ncsu.edu joint work with S. Cui, A. Kurganov, S.N. Özcan and E. Tadmor supported
More informationFifth order multi-moment WENO schemes for hyperbolic conservation laws
J. Sci. Comput. manuscript No. (will be inserted by the editor) Fifth order multi-moment WENO schemes for hyperbolic conservation laws Chieh-Sen Huang Feng Xiao Todd Arbogast Received: May 24 / Accepted:
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
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 informationNUMERICAL SOLUTION TO NONLINEAR TRICOMI EQUATION USING WENO SCHEMES
Eighth Mississippi State - UAB Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conf. 9 ), pp. 35 44. ISSN: 7-669. URL: http://ejde.math.txstate.edu
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 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 informationA Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes )
A Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes ) Takahiro MIYOSHI and Kanya KUSANO 1) Hiroshima University, Higashi-Hiroshima 739-856, Japan 1) Nagoya University,
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
More informationTHE ALTERNATING EVOLUTION METHODS FOR FIRST ORDER NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS
THE ALTERNATING EVOLUTION METHODS FOR FIRST ORDER NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS HAILIANG LIU TO MARSHALL SLEMROD ON HIS 7TH BIRTHDAY WITH FRIENDSHIP AND APPRECIATION. Abstract. In this paper,
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 information