21 results on '"Chi-Wang Shu"'
Search Results
2. Review Article:Review of Entropy Stable Discontinuous Galerkin Methods for Systems of Conservation Laws on Unstructured Simplex Meshes
- Author
-
Tianheng Chen Chi-Wang Shu
- Subjects
Conservation law ,Simplex ,Discontinuous Galerkin method ,Applied mathematics ,Polygon mesh ,General Medicine ,Mathematics - Published
- 2020
3. On New Strategies to Control the Accuracy of WENO Algorithm Close to Discontinuities II: Cell Averages and Multiresolution
- Author
-
Sergio Amat, Chi-Wang Shu, and Juan Ruiz
- Subjects
Computational Mathematics ,Signal processing ,Computer science ,Classification of discontinuities ,Control (linguistics) ,Algorithm - Published
- 2020
4. Runge-Kutta Discontinuous Galerkin Method with a Simple and Compact Hermite WENO Limiter on Unstructured Meshes
- Author
-
Jun Zhu, Jianxian Qiu, Chi-Wang Shu, and Xinghui Zhong
- Subjects
Physics::Computational Physics ,Polynomial ,Hermite polynomials ,Physics and Astronomy (miscellaneous) ,Computer science ,Mathematical analysis ,010103 numerical & computational mathematics ,01 natural sciences ,Mathematics::Numerical Analysis ,010101 applied mathematics ,Nonlinear system ,Runge–Kutta methods ,Discontinuous Galerkin method ,Limiter ,Applied mathematics ,Polygon mesh ,Convex combination ,0101 mathematics - Abstract
In this paper we generalize a new type of compact Hermite weighted essentially non-oscillatory (HWENO) limiter for the Runge-Kutta discontinuous Galerkin (RKDG) method, which was recently developed in [38] for structured meshes, to two dimensional unstructured meshes. The main idea of this HWENO limiter is to reconstruct the new polynomial by the usage of the entire polynomials of the DG solution from the target cell and its neighboring cells in a least squares fashion [11] while maintaining the conservative property, then use the classical WENO methodology to form a convex combination of these reconstructed polynomials based on the smoothness indicators and associated nonlinear weights. The main advantage of this new HWENO limiter is the robustness for very strong shocks and simplicity in implementation especially for the unstructured meshes considered in this paper, since only information from the target cell and its immediate neighbors is needed. Numerical results for both scalar and system equations are provided to test and verify the good performance of this new limiter.
- Published
- 2017
5. High Order Fixed-Point Sweeping WENO Methods for Steady State of Hyperbolic Conservation Laws and Its Convergence Study
- Author
-
Shuhai Zhang, Chi-Wang Shu, Yong-Tao Zhang, and Liang Wu
- Subjects
Conservation law ,Mathematical optimization ,Partial differential equation ,Physics and Astronomy (miscellaneous) ,MathematicsofComputing_NUMERICALANALYSIS ,010103 numerical & computational mathematics ,Fixed point ,01 natural sciences ,010101 applied mathematics ,symbols.namesake ,Nonlinear system ,Monotone polygon ,Rate of convergence ,Convergence (routing) ,Euler's formula ,symbols ,Applied mathematics ,0101 mathematics ,Mathematics - Abstract
Fixed-point iterative sweeping methods were developed in the literature to efficiently solve static Hamilton-Jacobi equations. This class of methods utilizes the Gauss-Seidel iterations and alternating sweeping strategy to achieve fast convergence rate. They take advantage of the properties of hyperbolic partial differential equations (PDEs) and try to cover a family of characteristics of the corresponding Hamilton-Jacobi equation in a certain direction simultaneously in each sweeping order. Different from other fast sweeping methods, fixed-point iterative sweeping methods have the advantages such as that they have explicit forms and do not involve inverse operation of nonlinear local systems. In principle, it can be applied in solving very general equations using any monotone numerical fluxes and high order approximations easily. In this paper, based on the recently developed fifth order WENO schemes which improve the convergence of the classical WENO schemes by removing slight post-shock oscillations, we design fifth order fixed-point sweeping WENO methods for efficient computation of steady state solution of hyperbolic conservation laws. Especially, we show that although the methods do not have linear computational complexity, they converge to steady state solutions much faster than regular time-marching approach by stability improvement for high order schemes with a forward Euler time-marching.
- Published
- 2016
6. Runge-Kutta Discontinuous Galerkin Method with a Simple and Compact Hermite WENO Limiter
- Author
-
Chi-Wang Shu, Jun Zhu, Jianxian Qiu, and Xinghui Zhong
- Subjects
Physics::Computational Physics ,Polynomial ,Conservation law ,Hermite polynomials ,Physics and Astronomy (miscellaneous) ,Compact stencil ,Mathematical analysis ,010103 numerical & computational mathematics ,01 natural sciences ,Mathematics::Numerical Analysis ,Euler equations ,010101 applied mathematics ,Runge–Kutta methods ,symbols.namesake ,Discontinuous Galerkin method ,symbols ,Limiter ,0101 mathematics ,Mathematics - Abstract
In this paper, we propose a new type of weighted essentially non-oscillatory (WENO) limiter, which belongs to the class of Hermite WENO (HWENO) limiters, for the Runge-Kutta discontinuous Galerkin (RKDG) methods solving hyperbolic conservation laws. This new HWENO limiter is a modification of the simple WENO limiter proposed recently by Zhong and Shu [29]. Both limiters use information of the DG solutions only from the target cell and its immediate neighboring cells, thus maintaining the original compactness of the DG scheme. The goal of both limiters is to obtain high order accuracy and non-oscillatory properties simultaneously. The main novelty of the new HWENO limiter in this paper is to reconstruct the polynomial on the target cell in a least square fashion [8] while the simple WENO limiter [29] is to use the entire polynomial of the original DG solutions in the neighboring cells with an addition of a constant for conservation. The modification in this paper improves the robustness in the computation of problems with strong shocks or contact discontinuities, without changing the compact stencil of the DG scheme. Numerical results for both one and two dimensional equations including Euler equations of compressible gas dynamics are provided to illustrate the viability of this modified limiter.
- Published
- 2016
7. Analysis of Sharp Superconvergence of Local Discontinuous Galerkin Method for One-Dimensional Linear Parabolic Equations
- Author
-
Yang Yang and Chi-Wang Shu
- Subjects
Computational Mathematics ,Discretization ,Discontinuous Galerkin method ,Mathematical analysis ,Piecewise ,Boundary value problem ,Superconvergence ,Computer Science::Numerical Analysis ,Parabolic partial differential equation ,Projection (linear algebra) ,Finite element method ,Mathematics::Numerical Analysis ,Mathematics - Abstract
In this paper, we study the superconvergence of the error for the local discontinuous Galerkin (LDG) finite element method for one-dimensional linear parabolic equations when the alternating flux is used. We prove that if we apply piecewise k-th degree polynomials, the error between the LDG solution and the exact solution is (k+2)-th order superconvergent at the Radau points with suitable initial discretization. Moreover, we also prove the LDG solution is (k + 2)-th order superconvergent for the error to a particular projection of the exact solution. Even though we only consider periodic boundary condition, this boundary condition is not essential, since we do not use Fourier analysis. Our analysis is valid for arbitrary regular meshes and for P k polynomials with arbitrary k � 1. We perform numerical experiments to demonstrate that the superconvergence rates proved in this paper are sharp.
- Published
- 2015
8. High Order Finite Difference Methods with Subcell Resolution for Stiff Multispecies Discontinuity Capturing
- Author
-
Chi-Wang Shu, Wei Wang, Helen C. Yee, Björn Sjögreen, and D. V. Kotov
- Subjects
Discontinuity (linguistics) ,Physics and Astronomy (miscellaneous) ,Flow (mathematics) ,Shock (fluid dynamics) ,Reaction step ,Inviscid flow ,Mathematical analysis ,Finite difference method ,Geometry ,Euler system ,Classification of discontinuities ,Mathematics - Abstract
In this paper, we extend the high order finite-difference method with sub- cell resolution (SR) in (34) for two-species stiff one-reaction models to multispecies and multireaction inviscid chemical reactive flows, which are significantly more difficult because of the multiple scales generated by different reactions. For reaction problems, when the reaction time scale is very small, the reaction zone scale is also small and the governing equations become very stiff. Wrong propagation speed of discontinu- ity may occur due to the underresolved numerical solution in both space and time. The present SR method for reactive Euler system is a fractional step method. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with certain computed flow variables in the shock region modified by the Harten subcell resolution idea. Several numerical examples of multispecies and multireaction reactive flows are performed in both one and two dimensions. Studies demonstrate that the SR method can capture the correct propaga- tion speed of discontinuities in very coarse meshes. AMS subject classifications: 65M06, 76V05, 80A32
- Published
- 2015
9. A Survey of High Order Schemes for the Shallow Water Equations
- Author
-
Chi-Wang Shu and Yulong Xing
- Subjects
Finite volume method ,Dimension (vector space) ,Discontinuous Galerkin method ,Mathematical analysis ,Finite difference ,Polygon mesh ,Classification of discontinuities ,Shallow water equations ,Finite element method ,Mathematics - Abstract
In this paper, we survey our recent work on designing high order positivity- preserving well-balanced finite difference and finite volume WENO (weighted essen- tially non-oscillatory) schemes, and discontinuous Galerkin finite element schemes for solving the shallow water equations with a non-flat bottom topography. These schemes are genuinely high order accurate in smooth regions for general solutions, are essentially non-oscillatory for general solutions with discontinuities, and at the same time they preserve exactly the water at rest or the more general moving water steady state solutions. A simple positivity-preserving limiter, valid under suitable CFL condi- tion, has been introduced in one dimension and reformulated to two dimensions with triangular meshes, and we prove that the resulting schemes guarantee the positivity of the water depth.
- Published
- 2014
10. A Conservative Lagrangian Scheme for Solving Compressible Fluid Flows with Multiple Internal Energy Equations
- Author
-
Qinghong Zeng, Juan Cheng, and Chi-Wang Shu
- Subjects
Mathematical optimization ,Physics and Astronomy (miscellaneous) ,Internal energy ,Compressible flow ,Euler equations ,symbols.namesake ,Flow (mathematics) ,Scheme (mathematics) ,symbols ,Applied mathematics ,Inertial confinement fusion ,Compressible fluid flow ,Lagrangian ,Mathematics - Abstract
Lagrangian methods are widely used in many fields for multi-material compressible flow simulations such as in astrophysics and inertial confinement fusion (ICF), due to their distinguished advantage in capturing material interfaces automatically. In some of these applications, multiple internal energy equations such as those for electron, ion and radiation are involved. In the past decades, several staggered-grid based Lagrangian schemes have been developed which are designed to solve the internal energy equation directly. These schemes can be easily extended to solve problems with multiple internal energy equations. However such schemes are typically not conservative for the total energy. Recently, significant progress has been made in developing cell-centered Lagrangian schemes which have several good properties such as conservation for all the conserved variables and easiness for remapping. However, these schemes are commonly designed to solve the Euler equations in the form of the total energy, therefore they cannot be directly applied to the solution of either the single internal energy equation or the multiple internal energy equations without significant modifications. Such modifications, if not designed carefully, may lead to the loss of some of the nice properties of the original schemes such as conservation of the total energy. In this paper, we establish an equivalency relationship between the cell-centered discretizations of the Euler equations in the forms of the total energy and of the internal energy. By a carefully designed modification in the implementation, the cell-centered Lagrangian scheme can be used to solve the compressible fluid flow with one or multiple internal energy equations and meanwhile it does not lose its total energy conservation property. An advantage of this approach is that it can be easily applied to many existing large application codes which are based on the framework of solving multiple internal energy equations. Several two dimensional numerical examples for both Euler equations and three-temperature hydrodynamic equations in cylindrical coordinates are presented to demonstrate the performance of the scheme in terms of symmetry preserving, accuracy and non-oscillatory performance.
- Published
- 2012
11. Improvement on Spherical Symmetry in Two-Dimensional Cylindrical Coordinates for a Class of Control Volume Lagrangian Schemes
- Author
-
Chi-Wang Shu and Juan Cheng
- Subjects
Class (set theory) ,symbols.namesake ,Classical mechanics ,Generalized coordinates ,Physics and Astronomy (miscellaneous) ,Log-polar coordinates ,symbols ,Circular symmetry ,Cylindrical coordinate system ,Volume element ,Control volume ,Lagrangian ,Mathematics - Abstract
In, Maire developed a class of cell-centered Lagrangian schemes for solving Euler equations of compressible gas dynamics in cylindrical coordinates. These schemes use a node-based discretization of the numerical fluxes. The control volume version has several distinguished properties, including the conservation of mass, momentum and total energy and compatibility with the geometric conservation law (GCL). However it also has a limitation in that it cannot preserve spherical symmetry for one-dimensional spherical flow. An alternative is also given to use the first order area-weighted approach which can ensure spherical symmetry, at the price of sacrificing conservation of momentum. In this paper, we apply the methodology proposed in our recent work to the first order control volume scheme of Maire in to obtain the spherical symmetry property. The modified scheme can preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid, and meanwhile it maintains its original good properties such as conservation and GCL. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the good performance of the scheme in terms of symmetry, non-oscillation and robustness properties.
- Published
- 2012
12. Conservative Semi-Lagrangian Finite Difference WENO Formulations with Applications to the Vlasov Equation
- Author
-
Chi-Wang Shu and Jing-Mei Qiu
- Subjects
Strang splitting ,Finite volume method ,Physics and Astronomy (miscellaneous) ,Advection ,Differential form ,Mathematical analysis ,Compressibility ,Vlasov equation ,Finite difference ,Landau damping ,Mathematics - Abstract
In this paper, we propose a new conservative semi-Lagrangian (SL) finite difference (FD) WENO scheme for linear advection equations, which can serve as a base scheme for the Vlasov equation by Strang splitting [4]. The reconstruction procedure in the proposed SL FD scheme is the same as the one used in the SL finite volume (FV) WENO scheme [3]. However, instead of inputting cell averages and approximate the integral form of the equation in a FV scheme, we input point values and approximate the differential form of equation in a FD spirit, yet retaining very high order (fifth order in our experiment) spatial accuracy. The advantage of using point values, rather than cell averages, is to avoid the second order spatial error, due to the shearing in velocity (v) and electrical field (E) over a cell when performing the Strang splitting to the Vlasov equation. As a result, the proposed scheme has very high spatial accuracy, compared with second order spatial accuracy for Strang split SL FV scheme for solving the Vlasov-Poisson (VP) system. We perform numerical experiments on linear advection, rigid body rotation problem; and on the Landau damping and two-stream instabilities by solving the VP system. For comparison, we also apply (1) the conservative SL FD WENO scheme, proposed in [22] for incompressible advection problem, (2) the conservative SL FD WENO scheme proposed in [21] and (3) the non-conservative version of the SL FD WENO scheme in [3] to the same test problems. The performances of different schemes are compared by the error table, solution resolution of sharp interface, and by tracking the conservation of physical norms, energies and entropies, which should be physically preserved.
- Published
- 2011
13. A Conservative Modification to the Ghost Fluid Method for Compressible Multiphase Flows
- Author
-
Wei Liu, Chi-Wang Shu, and Li Yuan
- Subjects
Momentum ,Third order ,Level set method ,Physics and Astronomy (miscellaneous) ,Discretization ,Computer science ,Multiphase flow ,Compressibility ,Finite difference ,Applied mathematics ,Thermodynamics ,Conservation of mass - Abstract
A conservative modification to the ghost fluid method (GFM) is developed for compressible multiphase flows. The motivation is to eliminate or reduce the conservation error of the GFM without affecting its performance. We track the conservative variables near the material interface and use this information to modify the numerical solution for an interfacing cell when the interface has passed the cell. The modification procedure can be used on the GFM with any base schemes. In this paper we use the fifth order finite difference WENO scheme for the spatial discretization and the third order TVD Runge-Kutta method for the time discretization. The level set method is used to capture the interface. Numerical experiments show that the method is at least mass and momentum conservative and is in general comparable in numerical resolution with the original GFM.
- Published
- 2011
14. Local Discontinuous Galerkin Methods for the Degasperis-Procesi Equation
- Author
-
Chi-Wang Shu and Yan Xu
- Subjects
Nonlinear system ,Physics and Astronomy (miscellaneous) ,Computer simulation ,Discontinuous Galerkin method ,Piecewise ,Applied mathematics ,Degasperis–Procesi equation ,Constant (mathematics) ,Galerkin method ,Stability (probability) ,Mathematics::Numerical Analysis ,Mathematics - Abstract
In this paper, we develop, analyze and test local discontinuous Galerkin (LDG) methods for solving the Degasperis-Procesi equation which contains nonlinear high order derivatives, and possibly discontinuous or sharp transition solutions. The LDG method has the flexibility for arbitrary h and p adaptivity. We prove the L2 stability for general solutions. The proof of the total variation stability of the schemes for the piecewise constant P0 case is also given. The numerical simulation results for different types of solutions of the nonlinear Degasperis-Procesi equation are provided to illustrate the accuracy and capability of the LDG method.
- Published
- 2011
15. A High Order Adaptive Finite Element Method for Solving Nonlinear Hyperbolic Conservation Laws
- Author
-
Jinchao Xu, Chi-Wang Shu, and Zhengfu Xu
- Subjects
Computational Mathematics ,Nonlinear system ,Conservation law ,Order (business) ,Applied mathematics ,Finite element method ,Mathematics - Published
- 2011
16. On the Order of Accuracy and Numerical Performance of Two Classes of Finite Volume WENO Schemes
- Author
-
Chi-Wang Shu, Mengping Zhang, and Rui Zhang
- Subjects
Mathematical optimization ,Finite volume method ,Physics and Astronomy (miscellaneous) ,Computer simulation ,Linear system ,Order of accuracy ,law.invention ,Nonlinear system ,law ,Applied mathematics ,Order (group theory) ,Polygon mesh ,Cartesian coordinate system ,Mathematics - Abstract
In this paper we consider two commonly used classes of finite volume weighted essentially non-oscillatory (WENO) schemes in two dimensional Cartesian meshes. We compare them in terms of accuracy, performance for smooth and shocked solutions, and efficiency in CPU timing. For linear systems both schemes are high order accurate, however for nonlinear systems, analysis and numerical simulation results verify that one of them (Class A) is only second order accurate, while the other (Class B) is high order accurate. The WENO scheme in Class A is easier to implement and costs less than that in Class B. Numerical experiments indicate that the resolution for shocked problems is often comparable for schemes in both classes for the same building blocks and meshes, despite of the difference in their formal order of accuracy. The results in this paper may give some guidance in the application of high order finite volume schemes for simulating shocked flows.
- Published
- 2011
17. Dissipative Numerical Methods for the Hunter-Saxton Equation
- Author
-
Yan Xu and Chi-Wang Shu
- Subjects
Computational Mathematics - Published
- 2010
18. STRONG STABILITY PRESERVING PROPERTY OF THE DEFERRED CORRECTION TIME DISCRETIZATION.
- Author
-
Yuan Liu, Chi-Wang Shu, and Mengping Zhang
- Subjects
- *
MESHFREE methods , *APPROXIMATION theory , *PARTIAL differential equations , *DIFFERENTIAL equations , *MATHEMATICS - Abstract
In this paper, we study the strong stability preserving (SSP) property of a class of deferred correction time discretization methods, for solving the method-of-fines schemes approximating hyperbolic partial differential equations. [ABSTRACT FROM AUTHOR]
- Published
- 2008
19. NUMERICAL BOUNDARY CONDITIONS FOR THE FAST SWEEPING HIGH ORDER WENO METHODS FOR SOLVING THE EIKONAL EQUATION.
- Author
-
Ling Huang, Chi-Wang Shu, and Mengping Zhang
- Subjects
- *
BOUNDARY value problems , *HAMILTON-Jacobi equations , *EIKONAL equation , *DIFFERENTIAL equations , *MATHEMATICS - Abstract
High order fast sweeping methods have been developed recently in the literature to solve static Hamilton-Jacobi equations efficiently. Comparing with the first order fast sweeping methods, the high order fast sweeping methods are more accurate, but they often require additional numerical boundary treatment for several grid points near the boundary because of the wider numerical stencil It is particularly important to treat the points near the inflow boundary accurately, as the information would flow into the computational domain and would affect global accuracy. In the literature, the numerical solution at these boundary points are either fixed with the exact solution, which is not always feasible, or computed with a first order discretization, which could reduce the global accuracy. In this paper, we discuss two strategies to handle the inflow boundary conditions. One is based on the numerical solutions of a first order fast sweeping method with several different mesh sizes near the boundary and a Richardson extrapolation, the other is based on a Lax-Wendroff type procedure to repeatedly utilizing the PDE to write the normal spatial derivatives to the inflow boundary in terms of the tangential derivatives, thereby obtaining high order solution values at the grid points near the inflow boundary. We explore these two approaches using the fast sweeping high order WENO scheme in [18] for solving the static Eikonal equation as a representative example. Numerical examples are given to demonstrate the performance of these two approaches. [ABSTRACT FROM AUTHOR]
- Published
- 2008
20. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT.
- Author
-
Zhengfu Xu and Chi-Wang Shu
- Subjects
- *
FINITE differences , *POLLUTANTS , *FLUX (Metallurgy) , *NUMERICAL analysis , *MATHEMATICAL analysis - Abstract
In this paper we further explore and apply our recent anti-diffusive flux corrected high order finite difference WENO schemes for conservation laws [18] to compute the Saint-Venant system of shallow water equations with pollutant propagation, which is described by a transport equation. The motivation is that the high order anti-diffusive WENO scheme for conservation laws produces sharp resolution of contact discontinuities while keeping high order accuracy for the approximation in the smooth region of the solution. The application of the anti-diffusive high order WENO scheme to the Saint-Venant system of shallow water equations with transport of pollutant achieves high resolution [ABSTRACT FROM AUTHOR]
- Published
- 2006
21. LOCAL DISCONTINUOUS GALERKIN METHODS FOR THREE CLASSES OF NONLINEAR WAVE EQUATIONS.
- Author
-
Yan Xu and Chi-wang Shu
- Subjects
- *
GALERKIN methods , *NONLINEAR wave equations , *KORTEWEG-de Vries equation , *NONLINEAR theories , *PARTIAL differential equations - Abstract
In this paper, we further develop the local discontinuous Galerkin method to solve three classes of nonlinear wave equations formulated by the genera] KdV-Burgers type equations, the general fifth-order KdV type equations and the fully nonlinear K(n, n, n) equations, and prove their stability for these genera] classes of nonlinear equations. The schemes we present extend the previous work of Yau and Shu [30, 31] and of Levy, Shu and Yah [24] on local discontinuous Galerkin method solving partial differential equations with higher spatial derivatives. Numerical examples for nonlinear problems am showa to illustrate the accuracy and capability of the methods. The numerical experiments include stationary solitons, soliton interactions and oscillatory solitary wave solutions, The numerical experiments also include the compacton solutions of a generalized fifthorder KdV equation in which the highest order derivative term is nonlinear and the fully nonlinear K(n, n, n) equations. [ABSTRACT FROM AUTHOR]
- Published
- 2004
Catalog
Discovery Service for Jio Institute Digital Library
For full access to our library's resources, please sign in.