ICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
|
|
- Kristina Horton
- 6 years ago
- Views:
Transcription
1 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 University of Texas at Austin Austin, Texas 787 Reference: Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao, "A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws," ICES REPORT 7-, The Institute for Computational Engineering and Sciences, The University of Texas at Austin, August 7.
2 A MULTILEVEL-WENO TECHNIQUE FOR SOLVING NONLINEAR CONSERVATION LAWS TODD ARBOGAST, CHIEH-SEN HUANG, AND XIKAI ZHAO Abstract. We present a multi-level WENO polynomial reconstruction technique. Polynomials of various orders defined over various size and shifted stencils may be combined by using novel weighting strategy. Numerical results show that the technique provides some improvement over standard WENO reconstructions. Weighted essentially non-oscillatory, polynomial reconstruction, hyperbolic con- Key words. servation law AMS subect classifications. 65M8, 65M6, 76M, 76M. Introduction. We consider finite volume and finite difference methods for systems of nonlinear, hyperbolic conservation laws of the form u t + f(x) x =, t >, x R, u R d, d. (.) (For simplicity, we do not treat boundary conditions, although we might also consider a bounded space domain if we impose periodic boundary conditions.) Weighted essentially non-oscillatory (WENO) schemes are a popular choice for solving such systems of equations. They allow one to reconstruct a high order version of the solution merely from approximations of cell averages (in finite volume schemes) or point values (in finite difference schemes). The key is to average approximations defined on various stencils, and to weight them so as to avoid stencils containing a shock discontinuity in the solution. We present a novel multi-level WENO polynomial reconstruction technique. Polynomials of various orders defined over various size and shifted stencils may be combined by using novel weighting strategy. Numerical results in -D for scalar problems and Euler systems show that the technique provides some improvement over standard WENO reconstructions.. Standard WENO framework. To establish the context for our new multilevel WENO technique and to fix notation, we review the standard theory. For simplicity of exposition, we consider the finite volume framework. The finite difference framework is similar. Partition time as = t < t < < t n. In the domain of (.), take arbitrary grid points < x < x < x <, and define the cell I i = [x i, x i+ ], and the length of the cell x i = x i+ x i. The mid point of I i is denoted as x i+ (x i + x i+ ). Let ū n i be the numerically approximated average of u on the cell I i at This work was supported by the U.S. National Science Foundation under grant DMS-875. Department of Mathematics, University of Texas, 55 Speedway, C, Austin, TX 787- and Institute for Computational Engineering and Sciences, University of Texas, East th St., C, Austin, TX (arbogast@ices.utexas.edu) Department of Applied Mathematics and National Center for Theoretical Sciences, National Sun Yat-sen University, Kaohsiung 8, Taiwan, R.O.C. (huangcs@math.nsysu.edu.tw) Department of Mathematics, University of Texas, 55 Speedway, C, Austin, TX 787- and Institute for Computational Engineering and Sciences, University of Texas, East th St., C, Austin, TX (xzhao@math.utexas.edu) =
3 ARBOGAST, HUANG, AND ZHAO time t n, so ū n i u(x, t n ) dx. x i I i Now assume u is smooth. For a kth order approximation of u on the given cell I i, we consider the ordered stencil S i, which contains r cells to the left and s cells to the right of I i with r + s + = k, i.e., S i = {I i r,..., I i,..., I i+s }, r + s + = k, (.) from which we obtain the kth order reconstruction polynomial P i by imposing the interpolation conditions P i (x) dx = ū m, for all I m S i. x m I m.. Background on WENO Reconstruction. Given the cell I i, suppose we are interested in a (k )st general WENO reconstruction. From now and for the remainder of the paper, we fix a value of i and drop it from the notation. First consider all the small stencils S k containing I i and with k cells (k ) S k = {I i k++,..., I i+ }, k. From each S k we can reconstruct P k of degree k. Moreover, we can define the large stencil S k = Sk and reconstruct a higher order polynomial P k of degree (k ). At a fixed point x, the polynomial P k can usually be written as a convex combination of P k, so P k (x) = α P k (x), α =. (.) where we drop the subscript i. We refer to α as a linear weight. When there are discontinuities in the data over stencil S i, we want to make use of the relatively small stencils on which u is smooth in order to achieve the essentially non-oscillatory property. The standard WENO reconstruction is where α = R(x) = ˆα m ˆα, ˆα = m α P k (x), α α =, (.) (ɛ + σ P k, for all, (.) ) τ for some smoothness indicator σ P k and exponent τ (usually τ = and ɛ = 6 > avoids any possibility of division by zero). We refer to α as a nonlinear weight. The smoothness indicator σ defined by Jiang and Shu in [] is normally used to measure the smoothness of the reconstruction polynomials on the cell I i. It is given by k ( xi+ d m σ P k = ( x i ) m P k(x) ) dx m dx, (.5) m= x i
4 A Multilevel-WENO Technique where again k is the degree of the polynomial P k. If the grid is uniform in the cell size x = h, then in the regions where u is smooth, the Taylor expansion of (.5) gives σ P k = (u h) ( + O(h )), which implies α = α + O(h ). If there are discontinuities in u within some of the stencils S k, then σ P = O(). So the weights of the smooth and non-smooth k stencils are O() and O(h τ ), respectively. This leads to an oder of accuracy of the reconstruction R decreasing to k... Levy-Puppo-Russo Reconstructions. In [], Levy, Puppo and Russo describe a third order compact, central WENO scheme. They use a somewhat different WENO reconstruction than the standard one (.), because the linear weights in (.) fail to exist when x = x i+. For the given cell I i, reconstruct the optimal quadratic polynomial P using the large stencil S = {I i, I i, I i+ } and two linear polynomials P and P from the stencils S = {I i, I i } and S = {I i, I i+ }, respectively. The idea is to introduce a centered quadratic polynomial p c such that P = α P + α P + α c p c, α =, α, {, c, }, p c = α c (P α P α P ), where α are the linear weights, subect only to the symmetry α = α. Then the WENO reconstruction is R = α P + α P + α c p c, α =, α, {, c, }, and the nonlinear weights α are computed by (.) with τ = and ɛ =. Here ɛ needs to satisfy: ) ɛ σ in smooth region, ) ɛ σ near discontinuities.. Multilevel-WENO. In this section, we present the multilevel-weno scheme in one space dimension. We proceed as if u were smooth. For the kth order reconstruction (k ), consider the full stencil S i as in (.). We will have hierarchically smaller stencils. Consider the small stencils at level l, where l k. We can construct an (l + )st order accurate reconstruction using the stencils S l+ with (l + ) cells contained in S i. Define S l+ = {I i l+,..., I i+ } S, max(, l r) min(l, s), from which we obtain the (l + )st order reconstruction polynomial P l. Then the kth order multilevel-weno reconstruction is given by the following algorithm, where the nonlinear weights are computed by (.) with corresponding linear weights and smoothness indicators, to be defined after presenting the algorithm. Step. For l =, reconstruct P from the stencil S, and P from the stencil S. Then combine these two polynomials to obtain the second order reconstruction R R = γ P + γ P, γ + γ =, γ, γ, with linear weights γ and smoothness indicators σ P. Perhaps the unbiased choice γ = γ = is most reasonable.
5 ARBOGAST, HUANG, AND ZHAO Step. For l, first construct the (l + )th order polynomial P l from the stencil S l+ for all. Then following Levy-Puppo-Russo s idea [], the reconstruction R l will be a convex combination of Rl and the centered polynomial p l,c = β l (P l αr l l ), that is, R l = α l R l + β l p l,c, α l + β l =, α l, β l, all, with linear weights α l, βl and using the modified smoothness indicators given below. We generally take α l = 5 and βl = 5. Step. If l = k, then S k is the same as the entire stencil S i, so we stop here. If l < k, then again, using the WENO methodology, we combine all the reconstructions R l to get the final (l + )th reconstruction Rl at level l: R l = γ l R l, γ l =, γ l, all, with linear weights γ and smoothness indicators σ P l... Smoothness indicators. We now define the smoothness indicators. For l k. (a) If = max(, l r), then { σ l,r = σ R l σ P l σ P l + P l, σ l,p = σ P l σ P l + P l. (.) (b) If max(, l r) < < min(l, s), then { σ l,r = σ R l σ P l σp l P σ l P+ l P, l σ,p l = σ P l σp l P σ l P+ l P. (.) l (c) If = min(l, s), then { σ l,r = σ R l σ P l σ P l P l, σ l,p = σ P l σ P l P l. (.). For l = k. { σ k,r = σ R k σ P k, σ,p k = σ P. k (.).. Level. Notice that for the above algorithm, we started from l =. So to add in the level l =, we can first obtain the reconstruction R and then implement Step to Step 5 recursively starting from l =. We will describe the method for reconstructing R below. When l =, we only have one stencil S = {I i } from which we can obtain a first order reconstruction polynomial P = ū i. We will reconstruct R in two steps.
6 A Multilevel-WENO Technique 5. First, reconstruct P from the stencil S and P from the stencil S. Then combine P and P to get R R = α P + β P, =,, where the nonlinear weights are defined by (.) with smoothness indicators σ P and σ P. If we use (.5) to compute the smoothness indicator σ P, then σ P =. This will result in the nonlinear weights being highly biased to α, and therefore we lose accuracy. We want to define σ P such that it has order O(h ) ust like the other smoothness indicators in the smooth region, while maintaining the higher order polynomials when the data fits. So define σ P = σ P +P = (ū i+ ū i ). (.5). Next, combine the reconstructions R and R to get R = γ R + γ R, and the nonlinear weights are computed by (.) with linear weight γ and smoothness indicator σ P. Remark. Suppose instead of using the above method, we let R = P. Then in Step we have R = α R + β p, =,, where we use the smoothness indicator σ R = σ P given by (.5) and σ P to compute nonlinear weights. Near the shock, in theory we expect this modified method to cause oscillations. For instance, when there are shocks on one side of the cell I i, say σ P = O(h ), σ P = O(), then σ R may be much larger than σ P so that the nonlinear weights bias to β and R p. This will cause overshooting or undershooting near the shocks. However, this method seems to give slightly better results than the original method given above for the examples we tested in this paper... Reduction. When implement the algorithm, we can choose which levels to use rather than using all the levels. For example, for the fifth order scheme, we can choose level to level simply by setting α =, for all, or only use level and level by setting α = and β =.. Numerical results in One Dimensional Space. In our one dimensional tests, the L and L errors are computed by the formulas u(x, t n ) dx ū n i x i i x i (L error), (.) I i max i u(x, t n ) dx ū n i x i (L error). (.) I i In all of the tests in this paper, for both third order and fifth order schemes, we use α l = 5, βl = 5 for all l and γ = γ =. For the fifth order scheme, use γ = γ = and γ =. We will not include level in our tests.
7 6 ARBOGAST, HUANG, AND ZHAO.. Linear advection. We first test u t + u x = with periodic boundary conditions on [, ] and initial condition u (x) = sin(πx), at time T = on both uniform and nonuniform grids using the third order ML-WENO scheme with ɛ =, see Tables.., and the fifth order ML-WENO scheme with ɛ = 6, see Tables... We show the results of Shu s linear test in Figures... Finally, we test the initial condition with a sawtooth pattern on a uniform grid on (, ) with cells; that is, ū i = if i is odd and i = if i is even. The results are shown in Figure.. N L error L order L error L order e e e-6..79e e e-7. Table.: Linear advection. ML-WENO error and convergence rate on the uniform grid at T =, with ɛ =, t =.5 x. N L error L order L error L order e e e-6..5e e e Table.: Linear advection. ML-WENO error and convergence rate on the nonuniform grid at T =, with ɛ =. N L error L order L error L order e e e e e e e e e e e e-.896 Table.: Linear advection. ML-WENO5 error and convergence rate on the uniform grid at T =, with ɛ = 6, t =. x.
8 A Multilevel-WENO Technique 7 N L error L order L error L order e e e e e e e e e e e e-.977 Table.: Linear advection. ML-WENO5 error and convergence rate on the nonuniform grid at T =, with ɛ = (a) time step. N = (b) time steps. N =. Fig..: Shu s linear test. t = x/. The plots are WENO5-JS (black line) and ML-WENO5 (blue line) (a) time step. N = (b) time steps. N =. Fig..: Shu s linear test. t = x/. The plots are WENO5-JS (black line) and ML-WENO5 with levels - (red line) (a) time step (b) time steps. Fig..: Sawtooth. t = x/. The plots are WENO5-JS (black line) and ML- WENO5 (blue line).
9 8 ARBOGAST, HUANG, AND ZHAO.. Burger s equation. We solve the Burger s equation u t + (u /) x = with the initial condition u (x) = sin(πx) on [, ] at time T =. The results are shown in Tables N L error L order L error L order e e e e e Table.5: Burger s equation. ML-WENO error and convergence rate on the uniform grid at T =, with ɛ =, t =. x. N L error L order L error L order e e e e e-6.99 Table.6: Burger s equation. ML-WENO error and convergence rate on the nonuniform grid at T =, with ɛ =. N L error L order L error L order e e e e e e e-8.97 Table.7: Burger s equation. ML-WENO5 error and convergence rate on the uniform grid at T =, with ɛ = 6, t =. x... Buckley-Leverett equation. In this example we test the Buckley-Leverett equation, which has the nonconvex flux f(u) = u u + ( u).
10 A Multilevel-WENO Technique 9 N L error L order L error L order e e e e e e e Table.8: Burger s equation. ML-WENO5 error and convergence rate on the nonuniform grid at T =, with ɛ = 6. We use the initial condition x, for x.5, u (x) =.5, for.5 x.,, otherwise. The results are shown in Figures. and.5 on uniform grids using N = cells at t =., t =., t =.6 and dt =. x. We take ɛ = 6 for both of the tests (a) t =.. (b) t =.. (c) t =.6. Fig..: Buckley-Leverett equation. N =, t =. x. The plots are WENO5- JS (black line) and ML-WENO5 (blue line) (a) t =.. (b) t =.. (c) t =.6. Fig..5: Buckley-Leverett equation. N =, t =. x. The plots are WENO5- JS (black line) and ML-WENO5 with levels - (red line).
11 ARBOGAST, HUANG, AND ZHAO.. D Euler system. The one dimensional Euler system of gas dynamics is given by ρ m + t x E m ρu + p =, (.) u(e + p) where m = ρu, E = p/(γ ) + ρu / and ρ, u, m, p, E are the density, velocity, momentum, pressure and energy, respectively, and γ =.. For all the tests below, we implemented ML-WENO with ɛ = 6, ML-WENO5 with ɛ = and WENO- JS with ɛ = 6, WENO5-JS with ɛ = Accuracy test. We first test the order of convergence on the Euler equations on a smooth problem. The initial condition is given by ρ(x, ) = +. sin(πx), u(x, ) =, p(x, ) =, with periodic boundary conditions. The exact solution is See Tables.9.. ρ(x, t) = +.sin(πx πt), u =, p =. N L error L order L error L order e-5..75e e e-6..8e e e e e-9..5e Table.9: Euler s equations. ML-WENO error and convergence rate on the nonuniform grid at T =.5 and ɛ =. N L error L order L error L order e e e-7.6.8e e e e-9.7.7e-9.97 Table.: Euler s equations. ML-WENO error and convergence rate on the nonuniform grid at T =.5 and ɛ = 6.
12 A Multilevel-WENO Technique N L error L order L error L order e e e e e e e e e e- 5.5 Table.: Euler s equations. ML-WENO5 error and convergence rate on the nonuniform grid at T = Shock tube test with Sod s and Lax s initial conditions. Sod s initial condition is { ρ l =, m l =, E l =.5, for x <.5, (ρ, m, E) = ρ r =.5, m r =, E r =.5, for x >.5, and Lax s initial condition is { ρ l =.5, m l =., E l = 8.98, for x <.5, (ρ, m, E) = ρ r =.5, m r =, E r =.75. for x >.5. The results are shown in Figures.6.8 using cells and t =. x. The ML- WENO5 result is somewhat better than WENO5-JS (a) Sod. N =. Third order scheme (b) Sod. N =. Zoom (c) Lax. N =. Third order scheme (d) Lax. N =. Zoom. Fig..6: Sod and Lax density at T =.6 using N = cells. The plots are the reference solution (green), WENO-JS (black) and ML-WENO (blue).
13 ARBOGAST, HUANG, AND ZHAO (a) Sod. N =. Fifth order scheme (b) Sod. N =. Zoom (c) Lax. N =. Fifth order scheme (d) Lax. N =. Zoom. Fig..7: Sod and Lax density at T =.6 using N = cells. The plots are the reference solution (green), WENO5-JS (black) and ML-WENO5 (blue) (a) Sod. N =. Fifth order scheme (b) Sod. N =. Zoom (c) Lax. N =. Fifth order scheme (d) Lax. N =. Zoom. Fig..8: Sod and Lax density at T =.6 using N = cells. The plots are the reference solution (green), WENO5-JS (black) and ML-WENO5 with levels - (red).
14 A Multilevel-WENO Technique... Shu and Osher s shock interaction with entropy waves. Next we solve the shock interaction with entropy waves problem given in [], which has a moving Mach shock interacting with sine waves in density. The initial condition is { ρ l =.857, u l =.6996, p l =., for < x <., (ρ, u, p) = ρ r = +.sin(5(x 5)), u r =, p r =, for. < x <. We compute the density at T =.6 using t =. x and N = cells. The results are shown in Figures (a) N =. Third order scheme (b) N =. Zoom. Fig..9: Shu and Osher s shock interaction with entropy waves density at T =.6 using N = cells. The plots are WENO-JS (black) and ML-WENO (blue) (a) N =. Fifth order scheme (b) N =. Zoom. Fig..: Shu and Osher s shock interaction with entropy waves density at T =.6 using N = cells. The plots are WENO5-JS (black) and ML-WENO5 (blue) (a) N =. Fifth order scheme (b) N =. Zoom. Fig..: Shu and Osher s shock interaction with entropy waves density at T =.6 using N = cells. The plots are WENO5-JS (black) and ML-WENO5 with levels - (red).
15 ARBOGAST, HUANG, AND ZHAO... Woodward and Colella s double blast test. This test uses the initial condition ρ l =, m l =, E l = /(γ ), for < x <., (ρ, m, E) = ρ m =, m m =, E m =./(γ ), for. < x <.9, ρ r =, m r =, E r = /(γ ), for.9 < x <. See Figures. and (a) N = 99. Third order scheme (b) N = 99. Zoom Fig..: Woodward and Colellas double blast test density at T =.8 using N = 99 cells. The plots are the reference solution (green), WENO-JS (black) and ML-WENO (blue) (a) N = 99. Fifth order scheme (b) N = 99. Zoom Fig..: Woodward and Colellas double blast test density at T =.8 using N = 99 cells. The plots are the reference solution (green), WENO5-JS (black), ML-WENO5 (blue) and ML-WENO5 with levels - (red). REFERENCES [] G.-S. Jiang and C.-W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys., 6 (996), pp. 8. [] R. J. Leveque, High-resolution conservative algorithms for advection in incompressible flow, SIAM J. Numer. Anal., (996), pp [] D. Levy, G. Puppo, and G. Russo, Compact central WENO schemes for multidimensional conservation laws, SIAM J. Sci. Comput., (), pp [] C.-W. Shu and S. Osher, Efficient implementation of essentially non-oscillatory shockcapturing schemes, ii, Journal of Computational Physics, 8 (989), pp. 78.
ICES REPORT Accuracy of WENO and Adaptive Order WENO Reconstructions for Solving Conservation Laws
ICES REPORT 17-29 October 2017 Accuracy of WENO and Adaptive Order WENO Reconstructions for Solving Conservation Laws by Todd Arbogast, Chieh-Sen Huang, Xikai Zhao The Institute for Computational Engineering
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 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 informationAn Eulerian-Lagrangian WENO Scheme for Nonlinear Conservation Laws
An Eulerian-Lagrangian WENO Scheme for Nonlinear Conservation Laws Chieh-Sen Huang a,, Todd Arbogast b, a Department of Applied Mathematics, National Sun Yat-sen University, Kaohsiung 84, Taiwan, R.O.C.
More informationA semi-lagrangian finite difference WENO scheme for scalar nonlinear conservation laws
A semi-lagrangian finite difference WENO scheme for scalar nonlinear conservation laws Chieh-Sen Huang a,, Todd Arbogast b,, Chen-Hui Hung c a Department of Applied Mathematics, National Sun Yat-sen University,
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 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 informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationA Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws
A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University
More informationAn 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationA minimum entropy principle of high order schemes for gas dynamics. equations 1. Abstract
A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang and Chi-Wang Shu 3 Abstract The entropy solutions of the compressible Euler equations satisfy a minimum principle
More informationA 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 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 informationExtremum-Preserving Limiters for MUSCL and PPM
arxiv:0903.400v [physics.comp-ph] 7 Mar 009 Extremum-Preserving Limiters for MUSCL and PPM Michael Sekora Program in Applied and Computational Mathematics, Princeton University Princeton, NJ 08540, USA
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 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 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 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 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 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 informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
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 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 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 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 informationTodd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
CONSERVATIVE CHARACTERISTIC METHODS FOR LINEAR TRANSPORT PROBLEMS Todd Arbogast Department of Mathematics and, (ICES) The University of Texas at Austin Chieh-Sen (Jason) Huang Department of Applied Mathematics
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 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 informationThe WENO Method for Non-Equidistant Meshes
The WENO Method for Non-Equidistant Meshes Philip Rupp September 11, 01, Berlin Contents 1 Introduction 1.1 Settings and Conditions...................... The WENO Schemes 4.1 The Interpolation Problem.....................
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 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 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 informationLinear high-resolution schemes for hyperbolic conservation laws: TVB numerical evidence 1
Linear high-resolution schemes for hyperbolic conservation laws: TVB numerical evidence C. Bona a,2, C. Bona-Casas a,3, J. Terradas b,4 arxiv:885v2 [gr-qc] 2 Dec 28 a Departament de Física. Universitat
More informationEfficient Implementation of Weighted ENO Schemes*
JOURNAL OF COMPUTATIONAL PHYSICS 126, 202 228 (1996) ARTICLE NO. 0130 Efficient Implementation of Weighted ENO Schemes* GUANG-SHAN JIANG AND CHI-WANG SHU Division of Applied Mathematics, Brown University,
More informationA 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 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 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 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 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 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 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 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 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 informationEntropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing
Entropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing Siddhartha Mishra Seminar for Applied Mathematics ETH Zurich Goal Find a numerical scheme for conservation
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 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 informationarxiv: v3 [math.na] 31 Jan 2017
Embedded WENO: a design strategy to improve existing WENO schemes arxiv:1607.05019v3 [math.na] 31 Jan 017 Bart S. van Lith a,1, Jan H.M. ten Thije Boonamp a, Wilbert L. IJzerman a,b a Department of Mathematics
More informationOn Convergence of Minmod-Type Schemes
On Convergence of Minmod-Type Schemes Sergei Konyagin, Boan Popov, and Ognian Trifonov April 16, 2004 Abstract A class of non-oscillatory numerical methods for solving nonlinear scalar conservation laws
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 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 informationFOURIER PADÉ APPROXIMATIONS AND FILTERING FOR SPECTRAL SIMULATIONS OF AN INCOMPRESSIBLE BOUSSINESQ CONVECTION PROBLEM
MATHEMATICS OF COMPUTATION Volume 7, Number 9, July 7, Pages 7 9 S -78(7)8- Article electronically published on February, 7 FOURIER PADÉ APPROXIMATIONS AND FILTERING FOR SPECTRAL SIMULATIONS OF AN INCOMPRESSIBLE
More informationDetonation Simulations with a Fifth-Order TENO Scheme
Commun. Comput. Phys. doi: 10.4208/cicp.OA-2018-0008 Vol. x, No. x, pp. 1-37 xxx 201x Detonation Simulations with a Fifth-Order TENO Scheme Haibo Dong 1, Lin Fu 2, Fan Zhang 3,, Yu Liu 4 and Jun Liu 1
More informationOn the Reduction of Numerical Dissipation in Central-Upwind Schemes
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol., No., pp. 4-63 Commun. Comput. Phys. February 7 On the Reduction of Numerical Dissipation in Central-Upwind Schemes Alexander Kurganov, and Chi-Tien Lin Department
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 informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More informationSemi-discrete central-upwind schemes with reduced dissipation for Hamilton Jacobi equations
IMA Journal of Numerical Analysis 005 5, 113 138 doi: 10.1093/imanum/drh015 Semi-discrete central-upwind schemes with reduced dissipation for Hamilton Jacobi equations STEVE BRYSON Programme in Scientific
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 informationUnderstanding the Shu Osher Conservative Finite Difference Form
Journal of Scientific Computing, Vol. 9, Nos. 3, December 003 ( 003) Understanding the Shu Osher Conservative Finite Difference Form Barry Merriman Received June 4, 00; accepted (in revised form) October
More informationarxiv: v1 [math.na] 11 Jun 2018
High-order residual distribution scheme for the time-dependent Euler equations of fluid dynamics arxiv:1806.03986v1 [math.na] 11 Jun 2018 R. Abgrall, P. Bacigaluppi, S. Tokareva Institute of Mathematics,
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 informationENO RECONSTRUCTION AND ENO INTERPOLATION ARE STABLE. Communicated by Nira Dyn
ENO RECONSTRUCTION AND ENO INTERPOLATION ARE STABLE ULRIK S. FJORDHOLM, SIDDHARTHA MISHRA, AND EITAN TADMOR Communicated by Nira Dyn ABSTRACT. We prove that the ENO reconstruction and ENO interpolation
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
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 informationSemi-Discrete Central-Upwind Schemes with Reduced Dissipation for Hamilton-Jacobi Equations
Semi-Discrete Central-Upwind Schemes with Reduced Dissipation for Hamilton-Jacobi Equations Steve Bryson, Aleander Kurganov, Doron Levy and Guergana Petrova Abstract We introduce a new family of Godunov-type
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 information1. Introduction. We are interested in the scalar hyperbolic conservation law
SIAM J. NUMER. ANAL. Vol. 4, No. 5, pp. 1978 1997 c 005 Society for Industrial and Applied Mathematics ON CONVERGENCE OF MINMOD-TYPE SCHEMES SERGEI KONYAGIN, BOJAN POPOV, AND OGNIAN TRIFONOV Abstract.
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationAn Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term
An Accurate Deterministic Projection Method for Hyperbolic Systems with Stiff Source Term Alexander Kurganov Department of Mathematics, Tulane University, 683 Saint Charles Avenue, New Orleans, LA 78,
More informationDevice Benchmark Comparisons via Kinetic, Hydrodynamic, and High-Field Models
Device Benchmark Comparisons via Kinetic, Hydrodynamic, and High-Field Models Carlo Cercignani, Irene M. Gamba, Joseph W. Jerome, and Chi-Wang Shu Abstract This paper describes benchmark comparisons for
More informationBurgers equation - a first look at fluid mechanics and non-linear partial differential equations
Burgers equation - a first look at fluid mechanics and non-linear partial differential equations In this assignment you will solve Burgers equation, which is useo model for example gas dynamics anraffic
More informationA Central Discontinuous Galerkin Method for Hamilton-Jacobi Equations
J Sci Comput (00) 45: 404 48 DOI 0.007/s095-009-9340-y A Central Discontinuous Galerkin Method for Hamilton-Jacobi Equations Fengyan Li Sergey Yakovlev Received: 6 August 009 / Revised: November 009 /
More 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 informationThe Penultimate Scheme for Systems of Conservation Laws: Finite Difference ENO with Marquina s Flux Splitting 1
The Penultimate Scheme for Systems of Conservation Laws: Finite Difference ENO with Marquina s Flux Splitting 1 Ronald P. Fedkiw Computer Science Department Stanford University, Stanford, California 94305
More informationApplication of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms
Future Generation Computer Systems () 65 7 Application of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms R. Naidoo a,b, S. Baboolal b, a Department
More information