Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
|
|
- Michael Whitehead
- 6 years ago
- Views:
Transcription
1 Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University Collaborators: J. Zhu (NUAA), C.-W. Shu (Brown), H. Zhu (NUPT), X. Zhong (MSU), G. Li (Qingdao), Y. Cheng (Baidu Company, Beijing), M. Dumbser (Trento)
2 Outline Introduction Numerical results Conclusions
3 Introduction 1 Introduction We consider hyperbolic conservation laws: { u t + f(u) = 0, u(x, 0) = u 0 (x). Hyperbolic conservation laws and convection dominated PDEs play an important role arise in applications, such as gas dynamics, modeling of shallow waters,... There are special difficulties associated with solving these systems both on mathematical and numerical methods, for discontinuous may appear in the solutions for nonlinear equations, even though the initial conditions are smooth enough. 1
4 Introduction This is why devising robust, accurate and efficient methods for numerically solving these problems is of considerable importance and as expected, has attracted the interest of many researchers and practitioners. Within recent decades, many high-order numerical methods have been developed to solve these problem. Among them, we would like to mention Discontinuous Galerkin (DG) method and Weighted essentially non-oscillatory (WENO) scheme. DG method is a high order finite element method. WENO scheme is finite difference or finite volume scheme. Both DG and WENO are very important numerical methods for the Convection Dominated PDEs. 2
5 Introduction The first DG was presented by Reed and Hill in 1973, in the framework of neutron transport (steady state linear hyperbolic equations). From 1987, a major development of the DG method was carried out by Cockburn, Shu et al. in a series of papers. They established a framework to easily solve nonlinear time dependent hyperbolic conservation laws using explicit, nonlinearly stable high order Runge- Kutta time discretization and DG discretization in space. These methods are termed RKDG methods. DG employs useful features from high resolution finite volume schemes, such as the exact or approximate Riemann solvers serving as numerical fluxes, and limiters. 3
6 Introduction Limiter is an important component of RKDG methods for solving convection dominated problems with strong shocks in the solutions, which is applied to detect discontinuities and control spurious oscillations near such discontinuities. Many such limiters have been used in the literature on RKDG methods such as the minmod type TVB limiter by Coukburn and Shu et al., the moment based limiter developed by Flaherty et al.. Limiters have been an extensively studied subject for the DG methods, however it is still a challenge to find limiters which are robust, maintaining high order accuracy in smooth regions including at smooth extrema, and yielding sharp, non-oscillatory discontinuity transitions. 4
7 Introduction WENO schemes have following advantages: Uniform high order accuracy in smooth regions including at smooth extrema Sharp and essentially non-oscillatory (to the eyes) shock transition. Robust for many physical systems with strong shocks. Especially suitable for simulating solutions containing both discontinuities and complicated smooth solution structure, such as shock interaction with vortices. The limiters used to control spurious oscillations in the presence of strong shocks are less robust than the strategies of WENO finite volume and finite difference methods. In this presentation, we would like to show the design of a robust limiter for the RKDG methods based on WENO methods. 5
8 2 s We consider one dimensional conservation laws: u t + f(u) x = 0. Let x i are the centers of the cells I i = [x i 1 2, x i+ 1 2 ], x i = x i+ 1 2 x i 1 2, h = sup i x i. The solution and the test function space: V k h = {p : p Ii P k (I i )}. A local orthogonal basis over I i, v (i) 0 (x) = 1, v(i) 1 (x) = x x ( ) 2 i, v (i) x 2 x (x) = xi 1 i x i 12, 6
9 The numerical solution u h (x, t): u h (x, t) = k l=0 u (l) i (t)v (i) l (x), for x I i The degrees of freedom u (l) i (t) are the moments: u (l) i (t) = 1 u h (x, t)v (i) l (x)dx, l = 0, 1,, k a l I i where a l = I i (v (i) l (x)) 2 dx 7
10 In order to evolve the degrees of freedom u (l) i (t), we time equation u t +f(u) x = 0 with basis v (i) l (x), and integrate it on cell I i, using integration by part, we obtain: d dt u(l) i (t) + 1 ( f(u h (x, t)) d a l I i f(u h (x i 1/2, t))v (i) l (x i 1/2 ) ) dx v(i) l (x)dx + f(u h (x i+1/2, t))v (i) l (x i+1/2 ) = 0, l = 0, 1,, k However, the boundary terms f(u i+1/2 ) and v i+1/2 etc. are not well defined when u and v are in this space, as they are discontinuous at the cell interfaces. 8
11 From the conservation and stability (upwinding) considerations, we take A single valued monotone numerical flux to replace f(u i+1/2 ): ˆf i+1/2 = ˆf(u i+1/2, u+ i+1/2 ) where ˆf(u; u) = f(u) (consistency); ˆf(, ) (monotonicity) and ˆf is Lipschitz continuous with respect to both arguments. Values from inside I i for the test function v: v (i) l (x i+1/2 ), v(i) l (x + i 1/2 ) We get semi-discretization scheme: d dt u(l) i (t) + 1 ( f(u h (x, t)) d a l I i dx v(i) l (x)dx + ˆf(u i+1/2, u+ i+1/2 )v(i) l (x i+1/2 ) ˆf(u i 1/2, u i+1/2 )v(i) l (x + i 1/2 ) ) = 0, l = 0, 1,, k. ( ) 9
12 Using explicit, nonlinearly stable high order Runge-Kutta time discretizations. [Shu and Osher, JCP, 1988] The semidiscrete scheme ( ) is written as: u t = L(u) is discretized in time by a nonlinearly stable Runge-Kutta time discretization, e.g. the third order version. u (1) = u n + tl(u n ) u (2) = 3 4 un u(1) tl(u(1) ) u n+1 = 1 3 un u(2) tl(u(2) ). 10
13 Lax problem. t = cells. Density. Left: k = 1. Right: k = 2. k = 3 code blows up. For Blast Wave problem, code blows up for any k. 11
14 Limiters Many limiters have been used in the literature, such as: The minmod based TVB limiter.( Cockburn and Shu, Math. Comp. 1989) Moment limiter. (Biswas, Devine and Flaherty, Appl. Numer. Math, 1994) A modification of moment limiter.(burbean, Sagaut and Brunean, JCP, 2001) The monotonicity preserving (MP) limiter. (Suresh and Huynh, JCP, 1997) A modification of the MP limiter. (Rider and Margolin, JCP, 2001) 12
15 These limiters tend to degrade accuracy when mistakenly used in smooth regions of the solution. Burgers equation, initial condition u(x, 0) = sin(π(2x 10)), with periodic boundary condition, RKDG with TVB limiter, t=0.05. Cockburn and Shu, JSC (2001) 13
16 Burbean, Sagaut and Brunean, JCP, (2001) 14
17 WENO Type limiter In order to overcome the drawback of these limiters, from 2003, with my colleagues, we have studied using WENO as limiter for RKDG methods, with the goal of obtaining a robust and high order limiting procedure to simultaneously obtain uniform high order accuracy and sharp, non-oscillatory shock transition for RKDG methods. We separate limiter procedure into two parts: Identify the troubled cells, namely those cells which might need the limiting procedure; Reconstruct polynomials in troubled cells using WENO reconstruction which only maintain the original cell averages (conservation). 15
18 For the first part, we can use the following troubled-cell indicators: TVB: based on the TVB minmod function BDF: moment limiter of Biswas, Devine and Flaherty BSB: modified moment limiter of Burbeau et al. MP: monotonicity-preserving limiter MMP: modified monotonicity-preserving limiter KXRCF: A shock detector of Krivodonova et al., Applied Numer. Math (2004) Harten: Discontinuous detection technique based on Harten s subcell resolution, (Qiu and Shu, SISC, 2005). 16
19 TVB indicator Let ũ i = u h (x i+1/2 ) u(0) i, ũ i = u h (x + i 1/2 ) + u(0) i. These are modified by the modified minmod function where m is given by ũ (mod) i = m(ũ i, u (0) i+1 u(0) i, u (0) i u (0) i 1 ), ũ (mod) i = m(ũ i, u (0) i+1 u(0) i, u (0) i u (0) i 1 ), m(a 1, a 2,..., a n ) { a1 if a 1 M( x) 2, = m(a 1, a 2,..., a n ) otherwise. 17
20 The minmod function m is given by m(a 1, a 2,..., a n ) { s min1 j n a = i if sign(a 1 ) = = sign(a n ) = s, 0 otherwise. If ũ (mod) i ũ i or ũ (mod) i ũ i, we declare the cell I i as a troubled cell. 18
21 KXRCF indicator Partition the boundary of a cell I i into two portions Ii (inflow, v n < 0) and I i + (outflow, v n > 0). The cell I i is identified as a troubled cell, if I (u h i Ii u h Ini )ds > 1, I uh i Ii h k+1 2 i here h i is the radius of the circumscribed circle in the element I i. I ni is the neighbor of I i on the side of Ii and the norm is based on an element average in onedimensional case. 19
22 WENO reconstruction Reconstruct polynomials in troubled cells using WENO reconstruction which only maintain the original cell averages (conservation). x G is Gauss or Gauss-Lobatto quadrature point 20
23 S 0 : 1 x S 1 : 1 x S k : 1 x T : 1 x I i+l q 0 (x)dx = u (0) i+l, l = k,, 0; I i+l q 1 (x)dx = u (0) i+l, l = k + 1,, 1; I i+l q k (x)dx = u (0) i+l, l = 0,, k; I i+l Q(x)dx = u (0) i+l, l = k,, k; 21 First Prev Next Last Go Back Full Screen Close Quit
24 We find the combination coefficients, also called linear weights γ j, j = 0, 1,, k satisfying: A : I i Q(x)v (i) l (x)dx = k γ j j=0 B : Q(x G ) = I i q j (x)v (i) l (x)dx, l = 1,, k k γ j q j (x G ). j=0 22
25 We compute the smoothness indicator, denoted as β j for each stencil S j, which measures how smooth the function q j (x) on cell I i, β j = k l=1 xi+1/2 x i 1/2 ( x) 2l 1 (q (l) j )2 dx, where q (l) j is the lth-derivative of q j (x). We compute the nonlinear weight ω j based on the smoothness indicator ω j = α j k l=0 α l, with α j = γ j (ε + β j ) 2, j = 0, 1,..., k, where ε > 0 is a small number to avoid the denominator to become 0. 23
26 The final WENO approximation is then given by: A : u (l) i = 1 a l k ω j j=0 B : u(x G ) = I i q j (x)v (i) l (x)dx, l = 1,, k; k ω j q j (x G ). Reconstruction of moments based on the reconstructed point values: j=0 Remark I: u (l) i = x a l G w G u(x G )v (i) l (x G ), l = 1,, k. For procedure A, there are not the linear weights for P 3 case. 24
27 For procedure B: For the P 1 case, we use the two-point Gauss quadrature points. For the P 2 case, we use either the four-point Gauss-Lobatto quadrature points or three-point Gauss quadrature points. But there are negative linear weights when three-point Gauss quadrature points are used. For the P 3 case, we use the four-point Gauss quadrature points. Remark II: WENO limiters work well in all our numerical test cases, including 1D, 2D and 3D, structure and unstructured meshes, but for P 2 and P 3 cases, the compactness of DG is destroyed. 25
28 Hermite WENO (HWENO) reconstruction Reconstruct polynomials which maintain the original cell averages (conservation). 26
29 For P 2 case, we obtain following reconstructed polynomials: q 0 (x)dx = u (0) i+j a 0, j = 1, 0; q 0 (x)v (i 1) 1 (x)dx = u (1) i 1 a 1 I i+j I i 1 q 1 (x)dx = u (0) i+j a 0, j = 0, 1; q 1 (x)v (i+1) 1 (x)dx = u (1) i+1 a 1 I i+j I i+1 q 2 (x)dx = u (0) i+j a 0, j = 1, 0, 1 I i+j Q(x)dx = u (0) i+j a 0, j = 1, 0, 1; Q(x)v (i+j) 1 (x)dx = u (1) i+j a 1, j = 1, 1. I i+j I i+j Follow the routine A of WENO reconstruction, we can obtain new moment u (1) i. 27
30 To reconstruct u (2) i : q 0 (x)dx = u (0) i+j a 0, I i+j q 1 (x)dx = u (0) i+j a 0, I i+j q 2 (x)dx = u (0) I i+j I i+j Q(x)dx = u (0) i+j a 0, I i+j q 0 (x)v (i+j) 1 (x)dx = u (1) i+j a 1, j = 1, 0 q 1 (x)v (i+j) 1 (x)dx = u (1) i+j a 1, j = 0, 1 I i+j i+j a 0, j = 1, 0, 1; q 2 (x)v (i) 1 dx = u(1) i a 1 I i I i+j Q(x)v (i+j) 1 (x)dx = u (1) i+j a 1, j = 1, 0, 1, Follow the routine A of WENO reconstruction, we can obtain new moment u (2) i. Remark III: For P 3 case, we should extend stencil. 28
31 New HWENO type reconstruction I i is a troubled cell, we use stencil S = {I i 1, I i, I i+1 }. Denote the solutions of the DG method on these three cells as polynomials q 0 (x), q 1 (x) and q 2 (x), respectively. We would like to modify q 1 (x) to q new 1 (x). Procedure by Zhong and Shu, JCP (2013): In order to make sure that the reconstructed polynomial maintains the original cell average of q 1 in the target cell I i, the following modifications are taken: q 0 = 1 x i q 0 = q 0 q 0 + q 1, q 1 = q 1, q 2 = q 0 q 2 + q 1 I i q 0 (x)dx, q 1 = 1 x i I i q 1 (x)dx, q 2 = 1 x i I i q 2 (x)dx, 29
32 The final nonlinear WENO reconstruction polynomial q new 1 (x) is now defined by a convex combination of these modified polynomials: q new 1 (x) = ω 0 q 0 (x) + ω 1 q 1 (x) + ω 2 q 2 (x) If ω 0 + ω 1 + ω 2 = 1, then q new 1 has the same cell average and order of accuracy as q 1. Computational formula of ω 0, ω 1, and ω 2 are same as in WENO reconstruction. The linear weights can be chosen to be any set of positive numbers adding up to one. Since for smooth solutions the central cell is usually the best one, a larger linear weight is put on the central cell than on the neighboring cells, i.e. γ 0 < γ 1 and γ 1 > γ 2. 30
33 In Zhong and Shu, JCP (2013), they take: γ 0 = 0.001, γ 1 = γ 2 = which can maintain the original high order in smooth regions and can keep essentially non-oscillatory shock transitions in all their numerical examples. 31
34 Procedure by Zhu, Zhong, Shu and Q., 2014: In order to make sure that the reconstructed polynomial maintains the original cell average of q 1 in the target cell I i, the following modifications are taken: ( q 0 (x) q 0 (x)) 2 dx = min (φ(x) q 0 (x)) 2 dx I i 1 I i 1 ( q 2 (x) q 2 (x)) 2 dx = min (φ(x) q 2 (x)) 2 dx I i+1 I i+1 for φ(x) P k with I i φ(x)dx = I i q 1 (x)dx For notational consistency we also denote q 1 (x) = q 1 (x). Then we follow the routine of Zhong and Shu JCP (2013), and obtain the final nonlinear WENO reconstruction polynomial q new 1 (x). 32
35 For two dimensional case, we select the HWENO reconstruction stencil as S = {I i 1,j, I i,j 1, I i+1,j, I i,j+1, I ij } for simplify, we renumber these cells as I l, l = 0,, 4, and denote the DG solutions on these five cells to be p l (x, y), respectively. 33
36 { ( p l (x, y) p l (x, y)) 2 dxdy = min (φ(x, y) p l (x, y)) 2 dxdy I l I l ( 2 ( ) } 2 + (φ(x, y) p l (x, y))dxdy) + (φ(x, y) p l (x, y))dxdy, I l for φ(x, y) P k with I 4 φ(x, y)dxdy = I 4 p 4 (x, y)dxdy, where l = mod(l 1, 4) and l = mod(l + 1, 4). For notational consistency we also denote p 4 (x, y) = p 4 (x, y). Then we follow the routine of WENO reconstruction: Take linear weights: γ 0 = γ 1 = γ 2 = γ 3 = 0.001, γ 4 = I l 34
37 We compute the smoothness indicators, β l = k α =1 I ij α 1 where α = (α 1, α 2 ) and α = α 1 + α 2. I ij ( ) α 2 x α p 1 y α l (x, y) dxdy, l = 0,, 4, 2 We compute the non-linear weights based on the smoothness indicators. The final nonlinear HWENO reconstruction polynomial p new 4 (x, y) is defined by a convex combination of the (modified) polynomials in the stencil: p new p new 4 (x, y) = 4 ω l p l (x, y). l=0 4 (x, y) has the same cell average and order of accuracy as the original one p 4 (x, y) on condition that 4 l=0 ω l = 1. 35
38 Numerical results 3 Numerical results We show the the numerical results of one- and twodimensional cases to illustrate the performance of the WENO type limiters. Accuracy test Test cases with shock 36
39 Numerical results Nonuniform Meshes 37
40 Numerical results WENO limiter on unstructured meshes 38
41 Numerical results 2D Euler equation, HWENO limiter unstructured meshes 39
42 Numerical results 2D Euler quation, New HWENO limiter 40
43 Numerical results Lax problem: Euler equations with initial condition (ρ, v, p) = { (0.445, 0.698, 3.528) if x 0, (0.5, 0, 0.571) if x > 0. 41
44 Numerical results 42
45 Numerical results New HWENO limiter, P 1, P 2 and P 3 from left to right using KXRCF troubled-cell indicator. 43
46 Numerical results Blast wave problem: Euler equations with initial condition (1, 0, 1000) if 0 x < 0.1, (ρ, v, p) = (1, 0, 0.01) if 0.1 x < 0.9, (1, 0, 100) if 0.9 x First Prev Next Last Go Back Full Screen Close Quit
47 Numerical results 45
48 Numerical results New HWENO limiter, P 1, P 2 and P 3 from left to right using KXRCF troubled-cell indicator. 46
49 Numerical results Two-dimensional Euler equations The PDEs are ρ ρu ρu ρu ρv p ρuv + E u(e + p) t x ρv ρuv ρv 2 + p v(e + p) y = 0. 47
50 Numerical results Double mach reflection problem The computational domain for this problem is [0, 4] [0, 1]. The reflecting wall lies at the bottom, starting from x = 1 6. Initially a right-moving Mach 10 shock is positioned at x = 1 6, y = 0 and makes a 60 angle with the x-axis. For the bottom boundary, the exact post-shock condition is imposed for the part from x = 0 to x = 1 6 and a reflective boundary condition is used for the rest. At the top boundary, the flow values are set to describe the exact motion of a Mach 10 shock. We compute the solution up to t =
51 Numerical results Form top to bottom, WENO limiter ( cells), HWENO limiter ( cells) and New HWENO limiter ( cells). 49
52 Numerical results Forward step problem A Mach 3 wind tunnel with a step. The wind tunnel is 1 length unit wide and 3 length units long. The step is 0.2 length units high and is located 0.6 length units from the left-hand end of the tunnel. The problem is initialized by a right-going Mach 3 flow. Reflective boundary conditions are applied along the wall of the tunnel and in/out flow boundary conditions are applied at the entrance/exit. We compute the solution up to t = 4. 50
53 Numerical results Form top to bottom, WENO limiter ( cells), HWENO limiter ( cells) and New HWENO limiter ( cells). 51
54 Conclusions 4 Conclusions We have developed a new limiter for the RKDG methods solving hyperbolic conservation laws using finite volume high order WENO and HWENO reconstructions. First identify troubled cells by troubled cell indicator. Then reconstruct the polynomial solution inside the troubled cells by WENO type reconstruction using the cell averages and moments of neighboring cells, while maintaining the original cell averages of the troubled cells. Numerical results show that the methods are stable, accurate and robust in maintaining accuracy. Troubled cell indicator was used as discontinuous indicator for h-adaptive and r- adaptive methods (H. Zhu Q. JCP 2009, ACM 2013) and hybrid WENO methods (G. Li Q, JCP 2010, JSC 2012, ACM 2013 ). 52
55 Refereces 1. J. Qiu and C.-W. Shu: J. Comput. Phys., 193 (2004) J. Qiu and C.-W. Shu: Computers & Fluids, 34 (2005) J. Qiu and C.-W. Shu: SIAM J. Sci. Comput., 26 (2005), J. Qiu and C.-W. Shu: SIAM J. Sci. Comput., 27 (2005), J. Zhu, J. Qiu, C.-W. Shu and M. Dumbser: J. Comput. Phys., 227 (2008) J. Zhu and J. Qiu: J. Sci. Comput., 39 (2009), H. Zhu and J. Qiu: J. Comput. Phys., 228 (2009), G. Li and J. Qiu: J. Comput. Phys., 229 (2010), J. Zhu and J. Qiu: J. Comput. Phys., 230 (2011), J. Zhu and J. Qiu: Commun.Comput. Phys., 11 (2012), G. Li, C. Lu and J. Qiu: J. Sci. Comput., 51 (2012), J. Zhu and J. Qiu: J. Sci. Comput., 55 (2013), H. Zhu, Y. Cheng and J. Qiu: Adv. Appl. Math. Mech., 5 (2013), J. Zhu, X. Zhong, C.-W. Shu and J. Qiu: J. Comput. Phys., 248 (2013), H. Zhu and J. Qiu: Adv. Comput. Math., 39(2013), G. Li and J. Qiu: Adv. Comput. Math., 40 (2014), J. Zhu, X. Zhong, C.-W. Shu and J. Qiu: Runge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter, submitted to Commun.Comput. Phys. 53
56 THANK YOU!
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationarxiv: v2 [math.na] 11 Jul 2014
Runge-Kutta Discontinuous Galerkin Method for Traffic Flow Model on Networks Suncica Canic 1, Benedetto Piccoli, Jing-Mei Qiu 3, Tan Ren 4 arxiv:1403.3750v [math.na] 11 Jul 014 Abstract. We propose a bound-preserving
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 informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
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 informationDiscontinuous Galerkin method for convection-dominated time-dependent PDEs Chi-Wang Shu
Discontinuous Galerkin method for convection-dominated time-dependent PDEs Chi-Wang Shu Division of Applied Mathematics Brown University Outline Hyperbolic equations Convection-diffusion equations Higher
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationRunge Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems
Journal of Scientific Computing, Vol. 16, No. 3, September 2001 ( 2001) Review Article Runge Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems Bernardo Cockburn 1 and Chi-Wang Shu
More informationAdaptive TVD-RK Discontinuous Galerkin Algorithms for Shallow Water Equations
Adaptive TVD-RK Discontinuous Galerkin Algorithms for Shallow Water Equations Thida Pongsanguansin, Khamron Mekchay and Montri Maleewong Abstract The adaptive Discontinuous Galerkin DG method for solving
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 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 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 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 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 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 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 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 Dept. of
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 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 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 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 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 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 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 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 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 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 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 informationIntroduction In this paper we review existing and develop new local discontinuous Galerkin methods for solving time dependent partial differential equ
Local Discontinuous Galerkin Methods for Partial Differential Equations with Higher Order Derivatives Jue Yan and Chi-Wang Shu 3 Division of Applied Mathematics Brown University Providence, Rhode Island
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 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 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 informationStrong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators
Journal of Scientific Computing, Vol. 8, No., February 3 ( 3) Strong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators Sigal Gottlieb
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 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 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 informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationStrong Stability Preserving Time Discretizations
AJ80 Strong Stability Preserving Time Discretizations Sigal Gottlieb University of Massachusetts Dartmouth Center for Scientific Computing and Visualization Research November 20, 2014 November 20, 2014
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 informationHigh order positivity-preserving discontinuous Galerkin methods for radiative transfer equations
High order positivity-preserving discontinuous Galerkin methods for radiative transfer equations Daming Yuan 1, Juan Cheng 2 and Chi-Wang Shu 3 Abstract The positivity-preserving property is an important
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 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 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 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 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 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 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 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 informationarxiv: v1 [math.na] 22 Nov 2018
Asymptotic preserving Deferred Correction Residual Distribution schemes Rémi Abgrall and Davide Torlo arxiv:1811.09284v1 [math.na] 22 Nov 2018 Abstract This work aims to extend the residual distribution
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 informationSpectral Difference Method for Unstructured Grids II: Extension to the Euler Equations
Journal of Scientific Computing, Vol. 32, No. 1, July 2007 ( 2006) DOI: 10.1007/s10915-006-9113-9 Spectral Difference Method for Unstructured Grids II: Extension to the Euler Equations Z. J. Wang, 1,4
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 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 information