A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
|
|
- Gladys Kelly
- 5 years ago
- Views:
Transcription
1 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 Engineering, University of Maryland - College Park January 8, 2018
2 Overview 1 WENO, CRWENO, and Central Schemes 2 Motivation 3 1D CCRWENO Scheme 4 2D Extension 5 Numerical Tests 6 Conclusion
3 WENO Reconstruction Weighted Essentially Non-Oscillatory Use a convex combination of linear subschemes: u j+1/2 = i ω i u (i) j+1/2
4 WENO Reconstruction Weighted Essentially Non-Oscillatory Use a convex combination of linear subschemes: u j+1/2 = i ω i u (i) j+1/2 (1) u (1) j+1/2 = ūj 2 6ūj 1+ 6 ūj ω 1 = 1 10 u (2) j+1/2 = 1 5 6ūj ω 2 = 3 6ūj 3ūj+1 5 u (3) j+1/2 = ω 3 = 3 3ūj 6ūj+1 6ūj+2 10
5 WENO Reconstruction Weighted Essentially Non-Oscillatory Use a convex combination of linear subschemes: u j+1/2 = i ω i u (i) j+1/2 (1) u (1) j+1/2 = ūj 2 6ūj 1+ 6 ūj ω 1 = 1 10 u (2) j+1/2 = 1 5 6ūj ω 2 = 3 6ūj 3ūj+1 5 u (3) j+1/2 = ω 3 = 3 3ūj 6ūj+1 6ūj+2 10 u is smooth in all stencils ω i ω i and Eq. (1) is O ( x 5).
6 WENO Reconstruction Weighted Essentially Non-Oscillatory Use a convex combination of linear subschemes: u j+1/2 = i ω i u (i) j+1/2 (1) u (1) j+1/2 = ūj 2 6ūj 1+ 6 ūj ω 1 = 1 10 u (2) j+1/2 = 1 5 6ūj ω 2 = 3 6ūj 3ūj+1 5 u (3) j+1/2 = ω 3 = 3 3ūj 6ūj+1 6ūj+2 10 u is smooth in all stencils ω i ω i and Eq. (1) is O ( x 5). u is not smooth in stencil i ω i 0.
7 Non-Oscillatory Weights Produce this behavior by scaling the ideal weights ω i according to smoothness indicators β i [Jiang & Shu, 1996]: β 1 = (u j 2 2u j 1 + u j ) (u j 2 4u j 1 + 3u j ) 2 β 2 = (u j 1 2u j + u j+1 ) (u j 1 u j+1 ) 2 β 3 = (u j 2u j+1 + u j+2 ) (3u j 4u j+1 + u j+2 ) 2 ω i α i = (ɛ + β i ) 2 ω i = α i m α m
8 CRWENO Schemes Ghosh and Baeder [Ghosh & Baeder, 2012] used compact subschemes to produce the Compact-Reconstruction WENO scheme: 2 3 u j 1/ u j+1/2 = 1 5 6ūj 1+ 6ūj ω 1 = u j 1/ u j+1/2 = ūj 6ūj+1 ω 2 = u j+1/ u j+3/2 = ūj 6ūj+1 ω 3 = 3 10 Solve a tridiagonal system for the interface values.
9 CRWENO Schemes Ghosh and Baeder [Ghosh & Baeder, 2012] used compact subschemes to produce the Compact-Reconstruction WENO scheme: 2 3 u j 1/ u j+1/2 = 1 5 6ūj 1+ 6ūj ω 1 = u j 1/ u j+1/2 = ūj 6ūj+1 ω 2 = u j+1/ u j+3/2 = ūj 6ūj+1 ω 3 = 3 10 Solve a tridiagonal system for the interface values. CRWENO has smaller truncation error than WENO.
10 CRWENO Schemes Ghosh and Baeder [Ghosh & Baeder, 2012] used compact subschemes to produce the Compact-Reconstruction WENO scheme: 2 3 u j 1/ u j+1/2 = 1 5 6ūj 1+ 6ūj ω 1 = u j 1/ u j+1/2 = ūj 6ūj+1 ω 2 = u j+1/ u j+3/2 = ūj 6ūj+1 ω 3 = 3 10 Solve a tridiagonal system for the interface values. CRWENO has smaller truncation error than WENO. Implicitness in space improves resolution of small-scale features.
11 Extension to Systems WENO and CRWENO reconstructions fit into the Godunov framework: reconstruct left- and right-biased interface values and use a Riemann solver. Oscillations can arise if reconstructions are not performed in characteristic variables ξ = Cu, where C is the matrix of left eigenvectors of the flux Jacobian.
12 Extension to Systems WENO and CRWENO reconstructions fit into the Godunov framework: reconstruct left- and right-biased interface values and use a Riemann solver. Oscillations can arise if reconstructions are not performed in characteristic variables ξ = Cu, where C is the matrix of left eigenvectors of the flux Jacobian. The tridiagonal system of CRWENO becomes a block-tridiagonal system.
13 Extension to Systems WENO and CRWENO reconstructions fit into the Godunov framework: reconstruct left- and right-biased interface values and use a Riemann solver. Oscillations can arise if reconstructions are not performed in characteristic variables ξ = Cu, where C is the matrix of left eigenvectors of the flux Jacobian. The tridiagonal system of CRWENO becomes a block-tridiagonal system. CRWENO with characteristic variables becomes less efficient than WENO with characteristic variables.
14 Extension to Systems WENO and CRWENO reconstructions fit into the Godunov framework: reconstruct left- and right-biased interface values and use a Riemann solver. Oscillations can arise if reconstructions are not performed in characteristic variables ξ = Cu, where C is the matrix of left eigenvectors of the flux Jacobian. The tridiagonal system of CRWENO becomes a block-tridiagonal system. CRWENO with characteristic variables becomes less efficient than WENO with characteristic variables. M1: Can we avoid using characteristic variables?
15 Extension to Systems WENO and CRWENO reconstructions fit into the Godunov framework: reconstruct left- and right-biased interface values and use a Riemann solver. Oscillations can arise if reconstructions are not performed in characteristic variables ξ = Cu, where C is the matrix of left eigenvectors of the flux Jacobian. The tridiagonal system of CRWENO becomes a block-tridiagonal system. CRWENO with characteristic variables becomes less efficient than WENO with characteristic variables. M1: Can we avoid using characteristic variables? M2: Can we eliminate the Riemann solver, and does doing so achieve M1?
16 Central Schemes A Riemann solver is needed to define a flux where the reconstruction is discontinuous. Central schemes avoid this by introducing a staggered grid: Staggered grid Main grid x j x j 1/2 x j+1/2
17 Central Schemes A Riemann solver is needed to define a flux where the reconstruction is discontinuous. Central schemes avoid this by introducing a staggered grid: Staggered grid Main grid x j x j 1/2 x j+1/2 Applying the integral form of a conservation law to a staggered cell gives: u (u) tn+1 + f = 0 ū n+1 t x j+1/2 = ūn j+1/2 f (u(x j+1 )) f (u(x j )) dx t n x
18 Time Advancement In Central Schemes Evaluate the time integrals by quadrature (Simpson s rule). Evolve the point values using Runge-Kutta with the differential form. Sample a natural continuous extension at quadrature nodes.
19 Time Advancement In Central Schemes RK4: Evaluate the time integrals by quadrature (Simpson s rule). Evolve the point values using Runge-Kutta with the differential form. Sample a natural continuous extension at quadrature nodes. u n+1 = u n + t 4 b k G k k=1
20 Time Advancement In Central Schemes RK4: Evaluate the time integrals by quadrature (Simpson s rule). Evolve the point values using Runge-Kutta with the differential form. Sample a natural continuous extension at quadrature nodes. u n+1 = u n + t 4 b k G k k=1 Natural Continuous Extension [Zennaro, 1986]: z(t n + θ t) = u n + t 4 b k (θ)g k = u(t) + O ( t 4) k=1 z(t n ) = u n and z(t n+1 ) = u n+1.
21 Central Scheme Steps 1 From cell averages on the main grid at time index n, reconstruct left-subcell averages. 2 Compute right-subcell averages by conservation. 3 Compute averages over staggered cells. 4 Reconstruct point values at staggered-cell interfaces (i.e., main-cell midpoints). 5 Compute fluxes and flux derivatives at staggered-cell interfaces. 6 Repeat step 5 as needed to calculate the time integrals in the update equation. 7 Compute staggered-cell averages at time index n + 1.
22 Central WENO Schemes Scheme Order Reconstruction Characteristic approach variables? [Levy, Puppo, & Russo 2002] 4th Polynomials No CWENO5 [Qiu & Shu, 2002] 5th FV formulas Yes, for SA
23 Central WENO Schemes Scheme Order Reconstruction Characteristic approach variables? [Levy, Puppo, & Russo 2002] 4th Polynomials No CWENO5 [Qiu & Shu, 2002] 5th FV formulas Yes, for SA Observation: The polynomial approach guarantees identical weights for the subcell-average and point-value reconstructions. Those weights do not match in CWENO5.
24 Central WENO Schemes Scheme Order Reconstruction Characteristic approach variables? [Levy, Puppo, & Russo 2002] 4th Polynomials No CWENO5 [Qiu & Shu, 2002] 5th FV formulas Yes, for SA Observation: The polynomial approach guarantees identical weights for the subcell-average and point-value reconstructions. Those weights do not match in CWENO5. M3: Would a central WENO scheme using the same weights for these reconstructions still require characteristic variables?
25 Motivation 1 Restore the advantage of compact reconstructions for systems by avoiding the need for characteristic variables.
26 Motivation 1 Restore the advantage of compact reconstructions for systems by avoiding the need for characteristic variables. 2 Eliminate the Riemann solver.
27 Motivation 1 Restore the advantage of compact reconstructions for systems by avoiding the need for characteristic variables. 2 Eliminate the Riemann solver. 3 Test the significance of identical weights for subcell-average and point-value reconstructions.
28 Constraints On The CCRWENO Subschemes
29 Constraints On The CCRWENO Subschemes 1 SA and PV reconstructions must use the same ideal weights.
30 Constraints On The CCRWENO Subschemes 1 SA and PV reconstructions must use the same ideal weights. Normalization: 2 Ideal weights must be positive and sum to 1. 3 LHS coefficients must sum to 1 for all weight combinations.
31 Constraints On The CCRWENO Subschemes 1 SA and PV reconstructions must use the same ideal weights. Normalization: 2 Ideal weights must be positive and sum to 1. 3 LHS coefficients must sum to 1 for all weight combinations. Accuracy: 4 SA and PV reconstructions must be at least fifth-order accurate.
32 Constraints On The CCRWENO Subschemes 1 SA and PV reconstructions must use the same ideal weights. Normalization: 2 Ideal weights must be positive and sum to 1. 3 LHS coefficients must sum to 1 for all weight combinations. Accuracy: 4 SA and PV reconstructions must be at least fifth-order accurate. Stability: 5 LHS coefficient matrix must be diagonally dominant for all weight combinations.
33 Constraints On The CCRWENO Subschemes 1 SA and PV reconstructions must use the same ideal weights. Normalization: 2 Ideal weights must be positive and sum to 1. 3 LHS coefficients must sum to 1 for all weight combinations. Accuracy: 4 SA and PV reconstructions must be at least fifth-order accurate. Stability: 5 LHS coefficient matrix must be diagonally dominant for all weight combinations. Symmetry: 6 Ideal weights must be symmetric: ω 1 = ω 3. 7 LHS of the combined scheme with ideal weights must be symmetric.
34 Identical Coefficient Matrices If we use the same smoothness indicators for SA and PV reconstructions, we can merge them into one linear system by designing their subschemes so that the two coefficient matrices are identical. 8 LHS coefficients of each SA subscheme must match those of the corresponding PV subscheme.
35 Identical Coefficient Matrices If we use the same smoothness indicators for SA and PV reconstructions, we can merge them into one linear system by designing their subschemes so that the two coefficient matrices are identical. 8 LHS coefficients of each SA subscheme must match those of the corresponding PV subscheme. A one-parameter family of such schemes exists that satisfy all eight conditions.
36 The CCRWENO Subschemes: SA and PV LHS LHS coefficients of subschemes: L = 1 d d d d 2 1 d d 1 1 d 1 2 2
37 The CCRWENO Subschemes: SA and PV LHS LHS coefficients of subschemes: L = 1 d d d d 2 1 d d 1 1 d d 1 > 0, d 2 = 5 + 8d d 1 ω 1 = 1 ω 2, ω 2 = d , ω 3 = 1 ω 2 2 d 1 + d 2 2
38 The CCRWENO Subschemes: SA and PV RHS ū L j 1 d d d 2 16 L ūj L = 0 ūj+1 L 0 0 u j d d 2 L u j = 0 24 u j d d d 2 5d d 1 1 3d 1 8 d d d d d ū j 2 ū j 1 ū j ū j+1 ū j+2 ū j 2 ū j 1 ū j ū j+1 ū j+2
39 The CCRWENO Subschemes: SA and PV RHS ū L j 1 d d d 2 16 L ūj L = 0 ūj+1 L 0 0 u j d d 2 L u j = 0 24 u j d d d 2 5d d 1 1 3d 1 8 d d d d d ū j 2 ū j 1 ū j ū j+1 ū j+2 ū j 2 ū j 1 ū j ū j+1 ū j+2 The PV reconstruction is unstable for d All results will use d 1 = 1.3.
40 The CCRWENO Subschemes: FD Reconstruction Each time step involves four FD reconstructions.
41 The CCRWENO Subschemes: FD Reconstruction Each time step involves four FD reconstructions. Cannot reuse SA/PV coefficient matrix because either the weights or the LHS coefficients must be different in order to satisfy the accuracy conditions.
42 The CCRWENO Subschemes: FD Reconstruction Each time step involves four FD reconstructions. Cannot reuse SA/PV coefficient matrix because either the weights or the LHS coefficients must be different in order to satisfy the accuracy conditions. No symmetric and diagonally dominant scheme with positive weights exceeds the fourth-order accuracy provided by: x(u x ) j = 1 2 u j 2 2u j u j ω 1 = 1 6 x(u x ) j = 1 2 u j 1 + 2u j u j+1 ω 2 = 2 3 x(u x ) j = 3 2 u j + 2u j u j+2 ω 3 = 1 6 We use these subschemes for the flux derivative reconstruction.
43 2D Staggered Grid Main grid Staggered grid Integrate fluxes over faces by quadrature.
44 Outer Product Extension Compose a 1D subscheme r in the x direction with a 1D subscheme s in the y direction. Order of accuracy is preserved. x: L r mūi+m L = m Rr mū i+m m y: L s nūj+n L = n Rs nū j+n n m,n L r ml s nū LL i+m,j+n = m,n R r mr s nū i+m,j+n
45 Outer Product Extension Compose a 1D subscheme r in the x direction with a 1D subscheme s in the y direction. Order of accuracy is preserved. x: L r mūi+m L = m Rr mū i+m m y: L s nūj+n L = n Rs nū j+n n m,n L r ml s nū LL i+m,j+n = m,n R r mr s nū i+m,j+n j + 1 L r 1 Ls 1 L r 0 Ls 1 L r 1 Ls 1 j L r 1 Ls 0 L r 0 Ls 0 L r 1 Ls 0 ω r,s = ω r ω s j 1 L r 1 Ls 1 i 1 L r 0 Ls 1 i L r 1 Ls 1 i + 1
46 Linear Advection: Convergence Linear advection of sin 2 (πx) sin 2 (πy) with velocity (u, v) = (1, 1).
47 Linear Advection: Efficiency Compared to [Levy, Puppo, & Russo 2002] and the outer-product extension of CWENO5 in [Qiu & Shu, 2002].
48 Isentropic Vortex Advection: Convergence Smooth solution of the Euler equations.
49 Isentropic Vortex Advection: Efficiency
50 Lax Problem
51 Shu-Osher Problem
52 2D Riemann Problem: Configuration 5 Density contours for Configuration 5 in [Lax & Liu, 1998] computed with CCRWENO cells
53 2D Riemann Problem: Configuration 16 Density contours for Configuration 16 in [Lax & Liu, 1998] computed by CCRWENO cells
54 2D Riemann Problem: Configuration 16 A cross-section of the CCRWENO solution at x = 0.9 shows the small amplitude of the oscillations.
55 2D Riemann Problem: Configuration 16 Density contours computed by CWENO cells
56 Advection of High-Frequency Waves Euler advection of ρ(x, y) = 2 + sin(2πkx) sin(2πky) with velocity (u, v) = (1, 0) for a range of wavenumbers k.
57 Conclusion Original motivations: Avoid characteristic variables and preserve advantages of compact reconstructions.
58 Conclusion Original motivations: Avoid characteristic variables and preserve advantages of compact reconstructions. Small oscillations appear in the Lax problem and Configuration 16. CCRWENO does resolve fine features more efficiently.
59 Conclusion Original motivations: Avoid characteristic variables and preserve advantages of compact reconstructions. Small oscillations appear in the Lax problem and Configuration 16. CCRWENO does resolve fine features more efficiently. Avoid Riemann solvers. 2D physics captured without special treatment.
60 Conclusion Original motivations: Avoid characteristic variables and preserve advantages of compact reconstructions. Small oscillations appear in the Lax problem and Configuration 16. CCRWENO does resolve fine features more efficiently. Avoid Riemann solvers. 2D physics captured without special treatment. Avoid oscillations by using identical SA and PV weights. Fewer oscillations in Configuration 16.
61 Conclusion Original motivations: Avoid characteristic variables and preserve advantages of compact reconstructions. Small oscillations appear in the Lax problem and Configuration 16. CCRWENO does resolve fine features more efficiently. Avoid Riemann solvers. 2D physics captured without special treatment. Avoid oscillations by using identical SA and PV weights. Fewer oscillations in Configuration 16. CCRWENO is more computationally efficient than the Levy scheme or CWENO5.
62 Future Work Remove oscillations in the Lax problem and Configuration 16. Implement more general boundary conditions. Extend to Navier-Stokes equations.
63 References I Jiang, G.-S. and Shu, C.-W. (1996) Efficient Implementation of Weighted ENO Schemes Journal of Computational Physics 126 (), Ghosh, D. and Baeder, J. (2012) Compact reconstruction schemes with weighted ENO limiting for hyperbolic conservation laws SIAM Journal on Scientific Computing 34 (3), A1678 A1706 Levy, D. and Puppo, G. and Russo, G. (2002) A fourth-order central WENO scheme for multidimensional hyperbolic systems of conservation laws SIAM Journal on Scientific Computing 24 (2), Qiu, J. and Shu, C.-W. (2002) On the construction, comparison, and local characteristic decomposition for high-order central WENO schemes Journal of Computational Physics 183 (1), Zennaro, M. (1986) Natural continuous extensions of Runge-Kutta methods Mathematics of Computation 46 (173), Lax, P. and Liu, X.-D. (1998) Solution of two-dimensional Riemann problems of gas dynamics by positive schemes SIAM Journal on Scientific Computing 19 (2), Shi, J. and Hu, C. and Shu, C.-W. (2002) A technique of treating negative weights in WENO schemes Journal of Computational Physics 175 (1),
64 References II Liu, X.-D. and Osher, S. and Chan, T. (1994) Weighted essentially non-oscillatory schemes Journal of Computational Physics 115 (1), Lele, S. (1992) Compact finite difference schemes with spectral-like resolution Journal of Computational Physics 103 (1), Levy, D. and Puppo, G. and Russo, G. (1999) Central WENO schemes for hyperbolic systems of conservation laws ESAIM: Mathematical Modelling and Numerical Analysis 33 (3), Bianco, F. and Puppo, G. and Russo, G. (1999) High-order central schemes for hyperbolic systems of conservation laws SIAM Journal on Scientific Computing 21 (1), Levy, D. and Puppo, G. and Russo, G. (2000) A third order central WENO scheme for 2D conservation laws Applied Numerical Mathematics 33 (1-4), Levy, D. and Puppo, G. and Russo, G. (2000) Compact central WENO schemes for multidimensional conservation laws SIAM Journal on Scientific Computing 22 (2),
A 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 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 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 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 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 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 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 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 informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
More informationA Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes
A Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory SIAM Annual Meeting
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 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 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 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 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 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 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 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 informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
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 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 informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
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 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 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 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 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 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 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 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 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 informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More 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 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 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 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 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 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 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 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 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 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 informationThird Order Reconstruction of the KP Scheme for Model of River Tinnelva
Modeling, Identification and Control, Vol. 38, No., 07, pp. 33 50, ISSN 890 38 Third Order Reconstruction of the KP Scheme for Model of River Tinnelva Susantha Dissanayake Roshan Sharma Bernt Lie Department
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
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 informationICES 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 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 informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
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 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 informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationA minimum entropy principle of high order schemes for gas dynamics. equations 1. Abstract
A minimum entropy principle of high order schemes for gas dynamics equations iangxiong Zhang and Chi-Wang Shu 3 Abstract The entropy solutions of the compressible Euler equations satisfy a minimum principle
More informationAn 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 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 informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
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 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 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 informationSMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES
Proceedings of ALGORITMY 2016 pp. 155 164 SMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES TAMER H. M. A. KASEM AND FRANÇOIS G. SCHMITT Abstract. The weighted essentially non-oscillatory
More informationFinite difference approximation
Finite difference approximation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen January
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 informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationFinite-Volume-Particle Methods for Models of Transport of Pollutant in Shallow Water
Journal of Scientific Computing ( 006) DOI: 0.007/s095-005-9060-x Finite-Volume-Particle Methods for Models of Transport of Pollutant in Shallow Water Alina Chertock, Alexander Kurganov, and Guergana Petrova
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 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 informationEulerian Vortex Motion in Two and Three Dimensions
Eulerian Vortex Motion in Two and Three Dimensions Igor Yanovsky June 2006 1 Objective An Eulerian approach is used to solve the motion of an incompressible fluid, in two and three dimensions, in which
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 NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationAn Assessment of Semi-Discrete Central Schemes for Hyperbolic Conservation Laws
SANDIA REPORT SAND2003-3238 Unlimited Release Printed September 2003 An Assessment of Semi-Discrete Central Schemes for Hyperbolic Conservation Laws Mark A. Christon David I. Ketcheson Allen C. Robinson
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 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationFully implicit higher-order schemes applied to polymer flooding. K.-A. Lie(SINTEF), T.S. Mykkeltvedt(IRIS), X. Raynaud (SINTEF)
Fully implicit higher-order schemes applied to polymer flooding K.-A. Lie(SINTEF), T.S. Mykkeltvedt(IRIS), X. Raynaud (SINTEF) Motivation Challenging and complex to simulate water-based EOR methods, unresolved
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 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 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 informationADER schemes and high order coupling on networks of hyperbolic conservation laws
ADER schemes and high order coupling on networks of hyperbolic conservation laws Raul Borsche a,, Jochen Kall a a Erwin Schrödinger Straße, TU Kaiserslautern, Building 48, 67663 Kaiserslautern, Germany
More informationSpatial Discretization
Spatial Discretization Shengtai Li and Hui Li Los Alamos National Laboratory 1 Overview Finite Difference Method Finite Difference for Linear Advection Equation 3 Conservation Laws Modern Finite Difference
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 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 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 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 informationA new type of Hermite WENO schemes for Hamilton-Jacobi equations
A new type of Hermite WENO schemes for Hamilton-Jacobi equations Jun Zhu and Jianxian Qiu Abstract In this paper, a new type of Hermite weighted essentially non-oscillatory (HWENO) schemes is designed
More informationA gas-kinetic theory based multidimensional high-order method for the compressible Navier Stokes solutions
Acta Mech. Sin. 2017) 334):733 741 DOI 10.1007/s10409-017-0695-2 RESEARCH PAPER A gas-kinetic theory based multidimensional high-order method for the compressible Navier Stokes solutions Xiaodong Ren 1
More informationAdaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems
EPJ Web of Conferences 108, 0019 (016) DOI: 10.1051/ epjconf/ 0161080019 C Owned by the authors, published by EDP Sciences, 016 Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear
More informationRESEARCH HIGHLIGHTS. WAF: Weighted Average Flux Method
RESEARCH HIGHLIGHTS (Last update: 3 rd April 2013) Here I briefly describe my contributions to research on numerical methods for hyperbolic balance laws that, in my view, have made an impact in the scientific
More 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 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 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 informationA Hybrid WENO Scheme for Simulation of Shock Wave-Boundary Layer Interaction
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 29, Orlando, Florida AIAA 29-36 A Hybrid WENO Scheme for Simulation of Shock Wave-Boundary Layer
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
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 informationOn the Comparison of the Finite Volume and Discontinuous Galerkin Methods
Diploma Thesis Institute for Numerical Simulation, TUHH On the Comparison of the Finite Volume and Discontinuous Galerkin Methods Corrected version Katja Baumbach August 17, 2006 Supervisor: Prof. Dr.
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationarxiv: v2 [math.na] 1 Sep 2015
Space-time adaptive ADER discontinuous Galerkin finite element schemes with a posteriori sub-cell finite volume limiting Olindo Zanotti a, Francesco Fambri a, Michael Dumbser a, Arturo Hidalgo b a Department
More information