41 results on '"Chi-Wang Shu"'
Search Results
2. L$^2$ Error Estimate to Smooth Solutions of High Order Runge--Kutta Discontinuous Galerkin Method for Scalar Nonlinear Conservation Laws with and without Sonic Points
- Author
-
Jingqi Ai, Yuan Xu, Chi-Wang Shu, and Qiang Zhang
- Subjects
Numerical Analysis ,Computational Mathematics ,Applied Mathematics - Published
- 2022
3. An Oscillation-free Discontinuous Galerkin Method for Scalar Hyperbolic Conservation Laws
- Author
-
Yong Liu, Chi-Wang Shu, and Jianfang Lu
- Subjects
Numerical Analysis ,Computational Mathematics ,Conservation law ,Discontinuous Galerkin method ,Oscillation ,Applied Mathematics ,Scalar (mathematics) ,Mathematical analysis ,Superconvergence ,High order ,Spurious oscillations ,Mathematics - Abstract
In this paper, we propose a novel discontinuous Galerkin (DG) method to control the spurious oscillations when solving the scalar hyperbolic conservation laws. Usually, the high order linear numeri...
- Published
- 2021
4. Error Estimate of the Fourth-Order Runge--Kutta Discontinuous Galerkin Methods for Linear Hyperbolic Equations
- Author
-
Chi-Wang Shu, Yuan Xu, and Qiang Zhang
- Subjects
Physics::Computational Physics ,Numerical Analysis ,Applied Mathematics ,010103 numerical & computational mathematics ,Computer Science::Numerical Analysis ,01 natural sciences ,Energy analysis ,Mathematics::Numerical Analysis ,Computational Mathematics ,Runge–Kutta methods ,Fourth order ,Discontinuous Galerkin method ,Applied mathematics ,Condensed Matter::Strongly Correlated Electrons ,0101 mathematics ,Hyperbolic partial differential equation ,Mathematics - Abstract
In this paper we consider the Runge--Kutta discontinuous Galerkin (RKDG) method to solve linear constant-coefficient hyperbolic equations, where the fourth-order explicit Runge--Kutta time-marching...
- Published
- 2020
5. On New Strategies to Control the Accuracy of WENO Algorithms Close to Discontinuities
- Author
-
Chi-Wang Shu, Juan Ruiz, Sergio Amat, Fundación Séneca, Ministerio de Economía y Competitividad, and National Science Foundation (NSF)
- Subjects
Signal processing ,Numerical Analysis ,12 Matemáticas ,Applied Mathematics ,Order of accuracy ,Matemática Aplicada ,010103 numerical & computational mathematics ,Classification of discontinuities ,01 natural sciences ,Raising (metalworking) ,Computational Mathematics ,Nonlinear system ,Improved adaption to discontinuities ,New optimal weights ,WENO schemes ,0101 mathematics ,Control (linguistics) ,Algorithm ,Interpolation ,Mathematics - Abstract
This paper is devoted to the construction and analysis of new nonlinear optimal weights for weighted ENO (WENO) interpolation capable of raising the order of accuracy close to discontinuities. The new nonlinear optimal weights are constructed using a strategy inspired by the original WENO algorithm, and they work very well for corner or jump singularities, leading to optimal theoretical accuracy. This is the first part of a series of two papers. In this first part we analyze the performance of the new algorithms proposed for univariate function approximation in the point values (interpolation problem). In the second part, we will extend the analysis to univariate function approximation in the cell averages (reconstruction problem). Our aim is twofold: to raise the order of accuracy of the WENO type interpolation schemes both near discontinuities and in the interval which contains the singularity. The first problem can be solved using the new nonlinear optimal weights, but the second one requires a new strategy that locates the position of the singularity inside the cell in order to attain adaption. This new strategy is inspired by the ENO-SR schemes proposed by Harten [J. Comput. Phys., 83 (1989), pp. 148--184]. Thus, we will introduce two different algorithms in the point values. The first one can deal with corner singularities and jump discontinuities for intervals not containing the singularity. The second algorithm can also deal with intervals containing corner singularities, as they can be detected from the point values, but jump discontinuities cannot, as the information of their position is lost during the discretization process. As mentioned before, the second part of this work will be devoted to the cell averages and, in this context, it will be possible to work with jump discontinuities as well. The work of the authors was supported by the Programa de Apoyo a la Investigatión de la Fundación Séneca-Agencia de Ciencia y Tecnología de la Región de Murcia 20928/PI/18, by the national research project MTM2015-64382-P (MINECO/FEDER), and by National Science Foundation grant DMS-1719410.
- Published
- 2019
6. Strong Stability of Explicit Runge--Kutta Time Discretizations
- Author
-
Zheng Sun and Chi-Wang Shu
- Subjects
Physics::Computational Physics ,Numerical Analysis ,Conservation law ,Applied Mathematics ,Numerical Analysis (math.NA) ,010103 numerical & computational mathematics ,Computer Science::Numerical Analysis ,01 natural sciences ,Stability (probability) ,Mathematics::Numerical Analysis ,Computational Mathematics ,Runge–Kutta methods ,FOS: Mathematics ,Energy method ,Applied mathematics ,Mathematics - Numerical Analysis ,0101 mathematics ,Mathematics - Abstract
Motivated by studies on fully discrete numerical schemes for linear hyperbolic conservation laws, we present a framework on analyzing the strong stability of explicit Runge-Kutta (RK) time discretizations for semi-negative autonomous linear systems. The analysis is based on the energy method and can be performed with the aid of a computer. Strong stability of various RK methods, including a sixteen-stage embedded pair of order nine and eight, has been examined under this framework. Based on numerous numerical observations, we further characterize the features of strongly stable schemes. A both necessary and sufficient condition is given for the strong stability of RK methods of odd linear order.
- Published
- 2019
7. Superconvergence of Discontinuous Galerkin Method for Scalar Nonlinear Hyperbolic Equations
- Author
-
Waixiang Cao, Yang Yang, Chi-Wang Shu, and Zhimin Zhang
- Subjects
Numerical Analysis ,Applied Mathematics ,Mathematical analysis ,Scalar (mathematics) ,010103 numerical & computational mathematics ,Superconvergence ,Computer Science::Numerical Analysis ,01 natural sciences ,Mathematics::Numerical Analysis ,010101 applied mathematics ,Computational Mathematics ,Nonlinear system ,Exact solutions in general relativity ,Discontinuous Galerkin method ,Bounded function ,Piecewise ,0101 mathematics ,Hyperbolic partial differential equation ,Mathematics - Abstract
In this paper, we study the superconvergence behavior of the semi-discrete discontinuous Galerkin (DG) method for scalar nonlinear hyperbolic equations in one spatial dimension. Superconvergence results for problems with fixed and alternating wind directions are established. On the one hand, we prove that, if the wind direction is fixed (i.e., the derivative of the flux function is bounded away from zero), both the cell average error and numerical flux error at cell interfaces converge at a rate of $2k+1$ when upwind fluxes and piecewise polynomials of degree $k$ are used. Moreover, we also prove that the function value approximation of the DG solution is superconvergent at interior right Radau points, and the derivative value approximation is superconvergent at interior left Radau points, with an order of k+2 and k+1, respectively. As a byproduct, we show a (k+2)th order superconvergence of the DG solution towards the Gauss--Radau projection of the exact solution. On the other hand, superconvergence resu...
- Published
- 2018
8. Optimal Error Estimates of the Semidiscrete Central Discontinuous Galerkin Methods for Linear Hyperbolic Equations
- Author
-
Yong Liu, Chi-Wang Shu, and Mengping Zhang
- Subjects
Numerical Analysis ,Conservation law ,Degree (graph theory) ,Applied Mathematics ,010103 numerical & computational mathematics ,01 natural sciences ,Mathematics::Numerical Analysis ,law.invention ,010101 applied mathematics ,Computational Mathematics ,Dimension (vector space) ,Discontinuous Galerkin method ,law ,Piecewise ,Applied mathematics ,Polygon mesh ,Cartesian coordinate system ,0101 mathematics ,Hyperbolic partial differential equation ,Mathematics - Abstract
We analyze the central discontinuous Galerkin method for time-dependent linear conservation laws. In one dimension, optimal a priori $L^2$ error estimates of order $k+1$ are obtained for the semidiscrete scheme when piecewise polynomials of degree at most $k$ ($k\geq0$) are used on overlapping uniform meshes. We then extend the analysis to multidimensions on uniform Cartesian meshes when piecewise tensor-product polynomials are used on overlapping meshes. Numerical experiments are given to demonstrate the theoretical results.
- Published
- 2018
9. AN OSCILLATION-FREE DISCONTINUOUS GALERKIN METHOD FOR SCALAR HYPERBOLIC CONSERVATION LAWS.
- Author
-
JIANFANG LU, YONG LIU, and CHI-WANG SHU
- Subjects
GALERKIN methods ,CONSERVATION laws (Mathematics) ,CONSERVATION laws (Physics) ,LINEAR orderings ,DATA structures ,OSCILLATIONS - Abstract
In this paper, we propose a novel discontinuous Galerkin (DG) method to control the spurious oscillations when solving the scalar hyperbolic conservation laws. Usually, the high order linear numerical schemes would generate spurious oscillations when the solution of the hyperbolic conservation laws contains discontinuities. The spurious oscillations may be harmful to the numerical simulation, as it not only generates some artificial structures not belonging to the problems, but also causes many overshoots and undershoots that make the numerical scheme less robust. To overcome this difficulty, in this paper we introduce a numerical damping term to control spurious oscillations based on the classic DG formulation. In comparison to the classic DG method, the proposed DG method still maintains many good properties, such as the extremely local data structure, conservation, L²-boundedness, optimal error estimates, and superconvergence. We also provide some numerical examples to show the good performance of the proposed DG scheme and verify our theoretical results. [ABSTRACT FROM AUTHOR]
- Published
- 2021
- Full Text
- View/download PDF
10. Stability and Error Estimates of Local Discontinuous Galerkin Methods with Implicit-Explicit Time-Marching for Advection-Diffusion Problems
- Author
-
Chi-Wang Shu, Haijin Wang, and Qiang Zhang
- Subjects
Numerical Analysis ,Discretization ,Advection ,Applied Mathematics ,Mathematical analysis ,Stability (probability) ,Square (algebra) ,Mathematics::Numerical Analysis ,Computational Mathematics ,Discontinuous Galerkin method ,Diffusion (business) ,Constant (mathematics) ,Convection–diffusion equation ,Mathematics - Abstract
The main purpose of this paper is to analyze the stability and error estimates of the local discontinuous Galerkin (LDG) methods coupled with carefully chosen implicit-explicit (IMEX) Runge--Kutta time discretization up to third order accuracy for solving one-dimensional linear advection-diffusion equations. In the time discretization the advection term is treated explicitly and the diffusion term implicitly. There are three highlights of this work. The first is that we establish an important relationship between the gradient and interface jump of the numerical solution with the independent numerical solution of the gradient in the LDG methods. The second is that, by aid of the aforementioned relationship and the energy method, we show that the IMEX LDG schemes are unconditionally stable for the linear problems in the sense that the time-step $\tau$ is only required to be upper-bounded by a constant which depends on the ratio of the diffusion and the square of the advection coefficients and is independent...
- Published
- 2015
11. A NEW WENO-2r ALGORITHM WITH PROGRESSIVE ORDER OF ACCURACY CLOSE TO DISCONTINUITIES.
- Author
-
AMAT, SERGIO P., RUIZ, JUAN, CHI-WANG SHU, and YÁÑEZ, DIONISIO F.
- Subjects
ALGORITHMS ,PROBLEM solving - Abstract
In this article a modification of the algorithm for data discretized in the point values introduced in [S. Amat, J. Ruiz, and C.-W. Shu, Appl. Math. Lett., 105 (2020), pp. 106-298] is presented. In the aforementioned work, the authors managed to obtain an algorithm that reaches a progressive and optimal order of accuracy close to discontinuities for WENO-6. For higher orders, i.e., WENO-8, WENO-10, etc., it turns out that the previous algorithm presents some shadows in the detection of discontinuities, meaning that the order of accuracy is better than the one attained by WENO of the same order, but not optimal. In this article a modification of the smoothness indicators used in the original algorithm is presented. It is oriented to solve this problem and to attain a WENO-2r algorithm with progressive order of accuracy close to the discontinuities. We also show proofs for the accuracy and explicit formulas for all the weights used for any order 2r of the algorithm. [ABSTRACT FROM AUTHOR]
- Published
- 2020
- Full Text
- View/download PDF
12. ERROR ESTIMATE OF THE FOURTH-ORDER RUNGE-KUTTA DISCONTINUOUS GALERKIN METHODS FOR LINEAR HYPERBOLIC EQUATIONS.
- Author
-
YUAN XU, CHI-WANG SHU, and QIANG ZHANG
- Subjects
- *
GALERKIN methods , *LINEAR equations , *HYPERBOLIC differential equations , *TRANSFER matrix , *ESTIMATES , *EQUATIONS - Abstract
In this paper we consider the Runge-Kutta discontinuous Galerkin (RKDG) method to solve linear constant-coefficient hyperbolic equations, where the fourth-order explicit Runge- Kutta time-marching is used. By the aid of the equivalent evolution representation with temporal differences of stage solutions, we make a detailed investigation on the matrix transferring process about the energy equations and then present a sufficient condition to ensure the L²-norm stability under the standard Courant-Friedrichs-Lewy condition. If the source term is equal to zero, we achieve the strong (boundedness) stability without the matrix transferring process to multiple-steps time-marching of the RKDG method. By carefully introducing the reference functions and their projections, we obtain the optimal (or suboptimal) error estimate under a mild smoothness assumption on the exact solution, which is independent of the stage number of the RKDG method. Some numerical experiments are also given to verify our conclusions. [ABSTRACT FROM AUTHOR]
- Published
- 2020
- Full Text
- View/download PDF
13. Optimal Error Estimates of the Semidiscrete Local Discontinuous Galerkin Methods for High Order Wave Equations
- Author
-
Chi-Wang Shu and Yan Xu
- Subjects
Numerical Analysis ,Discretization ,Applied Mathematics ,Mathematical analysis ,Order of accuracy ,Wave equation ,Mathematics::Numerical Analysis ,Schrödinger equation ,Computational Mathematics ,symbols.namesake ,Third order ,Discontinuous Galerkin method ,symbols ,Jump ,Order (group theory) ,Mathematics - Abstract
In this paper, we introduce a general approach for proving optimal $L^2$ error estimates for the semidiscrete local discontinuous Galerkin (LDG) methods solving linear high order wave equations. The optimal order of error estimates holds not only for the solution itself but also for the auxiliary variables in the LDG method approximating the various order derivatives of the solution. Examples including the one-dimensional third order wave equation, one-dimensional fifth order wave equation, and multidimensional Schrodinger equation are explored to demonstrate this approach. The main idea is to derive energy stability for the various auxiliary variables in the LDG discretization by using the scheme and its time derivatives with different test functions. Special projections are utilized to eliminate the jump terms at the cell boundaries in the error estimate in order to achieve the optimal order of accuracy.
- Published
- 2012
14. Analysis of Optimal Superconvergence of Discontinuous Galerkin Method for Linear Hyperbolic Equations
- Author
-
Chi-Wang Shu and Yang Yang
- Subjects
Numerical Analysis ,Conservation law ,Discretization ,Applied Mathematics ,Mathematical analysis ,Superconvergence ,Computer Science::Numerical Analysis ,Projection (linear algebra) ,Finite element method ,Mathematics::Numerical Analysis ,Computational Mathematics ,Discontinuous Galerkin method ,Piecewise ,Hyperbolic partial differential equation ,Mathematics - Abstract
In this paper, we study the superconvergence of the error for the discontinuous Galerkin (DG) finite element method for linear conservation laws when upwind fluxes are used. We prove that if we apply piecewise $k$th degree polynomials, the error between the DG solution and the exact solution is ($k+2$)th order superconvergent at the downwind-biased Radau points with suitable initial discretization. Moreover, we also prove the DG solution is ($k+2$)th order superconvergent both for the cell averages and for the error to a particular projection of the exact solution. The superconvergence result in this paper leads to a new a posteriori error estimate. Our analysis is valid for arbitrary regular meshes and for $\mathcal{P}^k$ polynomials with arbitrary $k\geq1$, and for both periodic boundary conditions and for initial-boundary value problems. We perform numerical experiments to demonstrate that the superconvergence rate proved in this paper is optimal.
- Published
- 2012
15. Superconvergence of Discontinuous Galerkin Methods for Scalar Nonlinear Conservation Laws in One Space Dimension
- Author
-
Chi-Wang Shu, Qiang Zhang, Xiong Meng, and Boying Wu
- Subjects
Numerical Analysis ,Conservation law ,Applied Mathematics ,Scalar (mathematics) ,Mathematical analysis ,Superconvergence ,Upper and lower bounds ,Mathematics::Numerical Analysis ,Computational Mathematics ,Nonlinear system ,Exact solutions in general relativity ,Discontinuous Galerkin method ,Piecewise ,Mathematics - Abstract
In this paper, an analysis of the superconvergence property of the semidiscrete discontinuous Galerkin (DG) method applied to one-dimensional time-dependent nonlinear scalar conservation laws is carried out. We prove that the error between the DG solution and a particular projection of the exact solution achieves $\left(k + \frac32\right)$th order superconvergence when upwind fluxes are used. The results hold true for arbitrary nonuniform regular meshes and for piecewise polynomials of degree $k$ ($k \ge 1$), under the condition that $|f'(u)|$ possesses a uniform positive lower bound. Numerical experiments are provided to show that the superconvergence property actually holds true for nonlinear conservation laws with general flux functions, indicating that the restriction on $f(u)$ is artificial.
- Published
- 2012
16. Superconvergence of Discontinuous Galerkin and Local Discontinuous Galerkin Schemes for Linear Hyperbolic and Convection-Diffusion Equations in One Space Dimension
- Author
-
Chi-Wang Shu and Yingda Cheng
- Subjects
Numerical Analysis ,Conservation law ,Applied Mathematics ,Numerical analysis ,Mathematical analysis ,Superconvergence ,Mathematics::Numerical Analysis ,Computational Mathematics ,Discontinuous Galerkin method ,Piecewise ,Galerkin method ,Hyperbolic partial differential equation ,Linear equation ,Mathematics - Abstract
In this paper, we study the superconvergence property for the discontinuous Galerkin (DG) and the local discontinuous Galerkin (LDG) methods for solving one-dimensional time dependent linear conservation laws and convection-diffusion equations. We prove superconvergence towards a particular projection of the exact solution when the upwind flux is used for conservation laws and when the alternating flux is used for convection-diffusion equations. The order of superconvergence for both cases is proved to be $k+\frac{3}{2}$ when piecewise $P^k$ polynomials with $k\geq1$ are used. The proof is valid for arbitrary nonuniform regular meshes and for piecewise $P^k$ polynomials with arbitrary $k\geq1$, improving upon the results in [Y. Cheng and C.-W. Shu, J. Comput. Phys., 227 (2008), pp. 9612-9627], [Y. Cheng and C.-W. Shu, Computers and Structures, 87 (2009), pp. 630-641] in which the proof based on Fourier analysis was given only for uniform meshes with periodic boundary condition and piecewise $P^1$ polynomials.
- Published
- 2010
17. Stability Analysis and A Priori Error Estimates of the Third Order Explicit Runge–Kutta Discontinuous Galerkin Method for Scalar Conservation Laws
- Author
-
Chi-Wang Shu and Qiang Zhang
- Subjects
Numerical Analysis ,Conservation law ,Discretization ,Applied Mathematics ,Courant–Friedrichs–Lewy condition ,Mathematical analysis ,Scalar (mathematics) ,Computer Science::Numerical Analysis ,Mathematics::Numerical Analysis ,Computational Mathematics ,Runge–Kutta methods ,Discontinuous Galerkin method ,Mathematics ,Numerical stability ,Linear stability - Abstract
In this paper we present an analysis of the Runge-Kutta discontinuous Galerkin method for solving scalar conservation laws, where the time discretization is the third order explicit total variation diminishing Runge-Kutta method. We use an energy technique to prove the $\mathrm{L}^2$-norm stability for scalar linear conservation laws and to obtain a priori error estimates for smooth solutions of scalar nonlinear conservation laws. Quasi-optimal order is obtained for general numerical fluxes, and optimal order is given for upwind fluxes. The theoretical results are obtained for piecewise polynomials with any degree $k\geq1$ under the standard temporal-spatial CFL condition $\tau\leq\gamma h$, where $h$ and $\tau$ are the element length and time step, respectively, and the positive constant $\gamma$ is independent of $h$ and $\tau$.
- Published
- 2010
18. A Genuinely High Order Total Variation Diminishing Scheme for One-Dimensional Scalar Conservation Laws
- Author
-
Chi-Wang Shu and Xiangxiong Zhang
- Subjects
Numerical Analysis ,Conservation law ,Finite volume method ,Applied Mathematics ,Numerical analysis ,Finite difference ,Computational Mathematics ,Third order ,Norm (mathematics) ,Total variation diminishing ,Calculus ,Applied mathematics ,Flux limiter ,Mathematics - Abstract
It is well known that finite difference or finite volume total variation diminishing (TVD) schemes solving one-dimensional scalar conservation laws degenerate to first order accuracy at smooth extrema [S. Osher and S. Chakravarthy, SIAM J. Numer. Anal., 21 (1984), pp. 955-984], thus TVD schemes are at most second order accurate in the $L^1$ norm for general smooth and nonmonotone solutions. However, Sanders [Math. Comp., 51 (1988), pp. 535-558] introduced a third order accurate finite volume scheme which is TVD, where the total variation is defined by measuring the variation of the reconstructed polynomials rather than the traditional way of measuring the variation of the grid values. By adopting the definition of the total variation for the numerical solutions as in [R. Sanders, Math. Comp., 51 (1988), pp. 535-558], it is possible to design genuinely high order accurate TVD schemes. In this paper, we construct a finite volume scheme which is TVD in this sense with high order accuracy (up to sixth order) in the $L^1$ norm. Numerical tests for a fifth order accurate TVD scheme will be reported, which include test cases from traffic flow models.
- Published
- 2010
19. A Local Discontinuous Galerkin Method for the Camassa–Holm Equation
- Author
-
Yan Xu and Chi-Wang Shu
- Subjects
Numerical Analysis ,Camassa–Holm equation ,Computer simulation ,Applied Mathematics ,Numerical analysis ,Mathematical analysis ,Stability (probability) ,Mathematics::Numerical Analysis ,Computational Mathematics ,Nonlinear system ,Nonlinear Sciences::Exactly Solvable and Integrable Systems ,Discontinuous Galerkin method ,Galerkin method ,Mathematics ,Numerical stability - Abstract
In this paper, we develop, analyze, and test a local discontinuous Galerkin (LDG) method for solving the Camassa–Holm equation which contains nonlinear high-order derivatives. The LDG method has the flexibility for arbitrary h and p adaptivity. We prove the $L^2$ stability for general solutions and give a detailed error estimate for smooth solutions, and provide numerical simulation results for different types of solutions of the nonlinear Camassa–Holm equation to illustrate the accuracy and capability of the LDG method.
- Published
- 2008
20. Central Discontinuous Galerkin Methods on Overlapping Cells with a Nonoscillatory Hierarchical Reconstruction
- Author
-
Eitan Tadmor, Chi-Wang Shu, Mengping Zhang, and Yingjie Liu
- Subjects
Numerical Analysis ,Conservation law ,Finite volume method ,Applied Mathematics ,Mathematical analysis ,Godunov's scheme ,Order of accuracy ,Computational Mathematics ,symbols.namesake ,Riemann problem ,Discontinuous Galerkin method ,symbols ,Applied mathematics ,MUSCL scheme ,Galerkin method ,Mathematics - Abstract
The central scheme of Nessyahu and Tadmor [J. Comput. Phys., 87 (1990), pp. 408-463] solves hyperbolic conservation laws on a staggered mesh and avoids solving Riemann problems across cell boundaries. To overcome the difficulty of excessive numerical dissipation for small time steps, the recent work of Kurganov and Tadmor [J. Comput. Phys., 160 (2000), pp. 241-282] employs a variable control volume, which in turn yields a semidiscrete nonstaggered central scheme. Another approach, which we advocate here, is to view the staggered meshes as a collection of overlapping cells and to realize the computed solution by its overlapping cell averages. This leads to a simple technique to avoid the excessive numerical dissipation for small time steps [Y. Liu, J. Comput. Phys., 209 (2005), pp. 82-104]. At the heart of the proposed approach is the evolution of two pieces of information per cell, instead of one cell average which characterizes all central and upwind Godunov-type finite volume schemes. Overlapping cells lend themselves to the development of a central-type discontinuous Galerkin (DG) method, following the series of works by Cockburn and Shu [J. Comput. Phys., 141 (1998), pp. 199-224] and the references therein. In this paper we develop a central DG technique for hyperbolic conservation laws, where we take advantage of the redundant representation of the solution on overlapping cells. The use of redundant overlapping cells opens new possibilities beyond those of Godunov-type schemes. In particular, the central DG is coupled with a novel reconstruction procedure which removes spurious oscillations in the presence of shocks. This reconstruction is motivated by the moments limiter of Biswas, Devine, and Flaherty [Appl. Numer. Math., 14 (1994), pp. 255-283] but is otherwise different in its hierarchical approach. The new hierarchical reconstruction involves a MUSCL or a second order ENO reconstruction in each stage of a multilayer reconstruction process without characteristic decomposition. It is compact, easy to implement over arbitrary meshes, and retains the overall preprocessed order of accuracy while effectively removing spurious oscillations around shocks.
- Published
- 2007
21. High Resolution Schemes for a Hierarchical Size‐Structured Model
- Author
-
Jun Shen, Chi-Wang Shu, and Mengping Zhang
- Subjects
Numerical Analysis ,Applied Mathematics ,Numerical analysis ,Finite difference method ,Stability (learning theory) ,Upwind scheme ,Computational Mathematics ,Convergence (routing) ,Calculus ,Applied mathematics ,Flux limiter ,High-resolution scheme ,Mathematics ,Numerical stability - Abstract
In this paper we discuss two explicit finite difference schemes, namely a first order upwind scheme and a second order high resolution scheme, for solving a hierarchical size-structured population model with nonlinear growth, mortality, and reproduction rates. We prove stability and convergence for both schemes and provide numerical examples to demonstrate their capability in solving smooth and discontinuous solutions.
- Published
- 2007
22. Error Estimates to Smooth Solutions of Runge–Kutta Discontinuous Galerkin Method for Symmetrizable Systems of Conservation Laws
- Author
-
Chi-Wang Shu and Qiang Zhang
- Subjects
Numerical Analysis ,Discretization ,Applied Mathematics ,Courant–Friedrichs–Lewy condition ,Mathematical analysis ,Computer Science::Numerical Analysis ,Finite element method ,Mathematics::Numerical Analysis ,Piecewise linear function ,Computational Mathematics ,Runge–Kutta methods ,Discontinuous Galerkin method ,Piecewise ,Galerkin method ,Mathematics - Abstract
In this paper we study the error estimates to sufficiently smooth solutions of symmetrizable systems of conservation laws for the Runge–Kutta discontinuous Galerkin (RKDG) method. Time discretization is the second‐order explicit TVD (total variation diminishing) Runge–Kutta method, and the $\mathbb{P}^k$ (piecewise polynomial) finite element is used. When $k=1$ (piecewise linear finite element), the error estimate is obtained under the usual CFL condition $\dt\leq \beta h$ for nonlinear systems in one dimension and for linear systems in multiple space dimensions. Here, h is the maximum element length, τ is the time step, and β is a positive constant independent of h and τ. Error estimates for $\mathbb{P}^k$ finite elements with $k>1$ are obtained under a more restrictive CFL condition.
- Published
- 2006
23. THE L2 -NORM STABILITY ANALYSIS OF RUNGE--KUTTA DISCONTINUOUS GALERKIN METHODS FOR LINEAR HYPERBOLIC EQUATIONS.
- Author
-
YUAN XU, QIANG ZHANG, CHI-WANG SHU, and HAIJIN WANG
- Subjects
GALERKIN methods ,LINEAR equations ,TRANSFER matrix ,POLYNOMIALS - Abstract
In this paper we propose a simple and unified framework to investigate the L
2 -norm stability of the explicit Runge--Kutta discontinuous Galerkin (RKDG) methods when solving the linear constant-coefficient hyperbolic equations. Two key ingredients in the energy analysis are the temporal differences of numerical solutions in different Runge--Kutta stages and a matrix transferring process. Many popular schemes, including the fourth order RKDG schemes, are discussed in this paper to show that the presented technique is flexible and useful. Different performances in the L2 - norm stability of different RKDG schemes are carefully investigated. For some lower-degree piecewise polynomials, the monotonicity stability is proved if the stability mechanism can be provided by the upwind-biased numerical fluxes. Some numerical examples are also given. [ABSTRACT FROM AUTHOR]- Published
- 2019
- Full Text
- View/download PDF
24. ON NEW STRATEGIES TO CONTROL THE ACCURACY OF WENO ALGORITHMS CLOSE TO DISCONTINUITIES.
- Author
-
AMAT, SERGIO, RUIZ, JUAN, and CHI-WANG SHU
- Subjects
NONLINEAR analysis ,ALGORITHMS ,INTERPOLATION ,CELL physiology ,UNIVARIATE analysis - Abstract
This paper is devoted to the construction and analysis of new nonlinear optimal weights for weighted ENO (WENO) interpolation capable of raising the order of accuracy close to discontinuities. The new nonlinear optimal weights are constructed using a strategy inspired by the original WENO algorithm, and they work very well for corner or jump singularities, leading to optimal theoretical accuracy. This is the first part of a series of two papers. In this first part we analyze the performance of the new algorithms proposed for univariate function approximation in the point values (interpolation problem). In the second part, we will extend the analysis to univariate function approximation in the cell averages (reconstruction problem). Our aim is twofold: to raise the order of accuracy of the WENO type interpolation schemes both near discontinuities and in the interval which contains the singularity. The first problem can be solved using the new nonlinear optimal weights, but the second one requires a new strategy that locates the position of the singularity inside the cell in order to attain adaption. This new strategy is inspired by the ENO-SR schemes proposed by Harten [J. Comput. Phys., 83 (1989), pp. 148-184]. Thus, we will introduce two different algorithms in the point values. The first one can deal with corner singularities and jump discontinuities for intervals not containing the singularity. The second algorithm can also deal with intervals containing corner singularities, as they can be detected from the point values, but jump discontinuities cannot, as the information of their position is lost during the discretization process. As mentioned before, the second part of this work will be devoted to the cell averages and, in this context, it will be possible to work with jump discontinuities as well. [ABSTRACT FROM AUTHOR]
- Published
- 2019
- Full Text
- View/download PDF
25. Error Estimates to Smooth Solutions of Runge--Kutta Discontinuous Galerkin Methods for Scalar Conservation Laws
- Author
-
Qiang Zhang and Chi-Wang Shu
- Subjects
Numerical Analysis ,Conservation law ,Discretization ,Applied Mathematics ,Courant–Friedrichs–Lewy condition ,Mathematical analysis ,Scalar (mathematics) ,Mathematics::Numerical Analysis ,Piecewise linear function ,Computational Mathematics ,Runge–Kutta methods ,Monotone polygon ,Discontinuous Galerkin method ,Mathematics - Abstract
In this paper we study the error estimates to sufficiently smooth solutions of scalar conservation laws for Runge--Kutta discontinuous Galerkin (RKDG) methods, where the time discretization is the second order explicit total variation diminishing (TVD) Runge--Kutta method. Error estimates for the $\mathbb{P}^1$ (piecewise linear) elements are obtained under the usual CFL condition $\tau\leq \gamma h$ for general nonlinear conservation laws in one dimension and for linear conservation laws in multiple space dimensions, where h and $\tau$ are the maximum element lengths and time steps, respectively, and the positive constant $\gamma$ is independent of $h$ and $\tau$. However, error estimates for higher order $\mathbb{P}^k(k\geq 2)$ elements need a more restrictive time step $\tau\leq \gamma h^{4/3}$. We remark that this stronger condition is indeed necessary, as the method is linearly unstable under the usual CFL condition $\tau\leq\gamma h$ for the $\mathbb{P}^k$ elements of degree $k\geq 2$. Error estimates of $O(h^{k+1/2}+\tau^2)$ are obtained for general monotone numerical fluxes, and optimal error estimates of $O(h^{k+1}+\tau^2)$ are obtained for upwind numerical fluxes.
- Published
- 2004
26. A Local Discontinuous Galerkin Method for KdV Type Equations
- Author
-
Jue Yan and Chi-Wang Shu
- Subjects
010101 applied mathematics ,Numerical Analysis ,Computational Mathematics ,Applied Mathematics ,010103 numerical & computational mathematics ,0101 mathematics ,01 natural sciences - Published
- 2002
27. The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems
- Author
-
Chi-Wang Shu and Bernardo Cockburn
- Subjects
Numerical Analysis ,Partial differential equation ,Applied Mathematics ,Mathematical analysis ,Triangulation (social science) ,Mathematics::Numerical Analysis ,Computational Mathematics ,Nonlinear system ,Runge–Kutta methods ,Rate of convergence ,Discontinuous Galerkin method ,Convection–diffusion equation ,Galerkin method ,Mathematics - Abstract
In this paper, we study the local discontinuous Galerkin (LDG) methods for nonlinear, time-dependent convection-diffusion systems. These methods are an extension of the Runge--Kutta discontinuous Galerkin (RKDG) methods for purely hyperbolic systems to convection-diffusion systems and share with those methods their high parallelizability, high-order formal accuracy, and easy handling of complicated geometries for convection-dominated problems. It is proven that for scalar equations, the LDG methods are L2-stable in the nonlinear case. Moreover, in the linear case, it is shown that if polynomials of degree k are used, the methods are kth order accurate for general triangulations; although this order of convergence is suboptimal, it is sharp for the LDG methods. Preliminary numerical examples displaying the performance of the method are shown.
- Published
- 1998
28. On the Gibbs Phenomenon III: Recovering Exponential Accuracy in a Sub-Interval From a Spectral Partial Sum of a Pecewise Analytic Function
- Author
-
David Gottlieb and Chi-Wang Shu
- Subjects
Numerical Analysis ,Gegenbauer polynomials ,Applied Mathematics ,Spectral element method ,Mathematical analysis ,Trigonometric polynomial ,Exponential function ,Gibbs phenomenon ,Computational Mathematics ,symbols.namesake ,symbols ,Piecewise ,Legendre polynomials ,Mathematics ,Analytic function - Abstract
We continue the investigation of overcoming the Gibbs phenomenon, i.e., obtaining exponential accuracy at all points, including at the discontinuities themselves, from the knowledge of a spectral partial sum of a discontinuous but piecewise analytic function. We show that if we are given the first N expansion coefficients of an $L^2 $ function $f(x)$ in terms of either the trigonometric polynomials or the Legendre polynomials, we can construct an exponentially convergent approximation to the point values of $f(x)$ in any sub-interval in which it is analytic.
- Published
- 1996
29. Nonlinearly Stable Compact Schemes for Shock Calculations
- Author
-
Chi-Wang Shu and Bernardo Cockburn
- Subjects
Numerical Analysis ,Computational Mathematics ,Nonlinear system ,Conservation law ,Runge–Kutta methods ,Applied Mathematics ,Mathematical analysis ,Scalar (mathematics) ,Compact finite difference ,Galerkin method ,Finite element method ,Numerical stability ,Mathematics - Abstract
The applications of high-order, compact finite difference methods in shock calculations are discussed. The main concern is to define a local mean which will serve as a reference for introducing a local nonlinear limiting to control spurious numerical oscillations while maintaining the formal accuracy of the scheme. For scalar conservation laws, the resulting schemes can be proven total-variation stable in one space dimension and maximum-norm stable in multiple space dimensions. Numerical examples are shown to verify accuracy and stability of such schemes for problems containing shocks. These ideas can also be applied to other implicit schemes such as the continuous Galerkin finite element methods.
- Published
- 1994
30. On One-Sided Filters for Spectral Fourier Approximations of Discontinuous Functions
- Author
-
David Gottlieb, Chi-Wang Shu, and Wei Cai
- Subjects
Numerical Analysis ,Applied Mathematics ,Mathematical analysis ,Existence theorem ,Spectral line ,Periodic function ,Gibbs phenomenon ,Computational Mathematics ,symbols.namesake ,Discontinuity (linguistics) ,Fourier transform ,Fourier analysis ,symbols ,Spectral method ,Mathematics - Abstract
The existence of one-sided filters, for spectral Fourier approximations of discontinuous functions, which can recover spectral accuracy up to discontinuity from one side, was proved. A least square procedure was also used to construct such a filter and test it on several discontinuous functions numerically.
- Published
- 1992
31. High-Order Essentially Nonoscillatory Schemes for Hamilton–Jacobi Equations
- Author
-
Chi-Wang Shu and Stanley Osher
- Subjects
Numerical Analysis ,Computational Mathematics ,Discontinuity (linguistics) ,Conservation law ,Applied Mathematics ,Mathematical analysis ,Initial value problem ,Differential (infinitesimal) ,Classification of discontinuities ,Optimal control ,Hyperbolic partial differential equation ,Hamilton–Jacobi equation ,Mathematics - Abstract
Hamilton–Jacobi (H–J) equations are frequently encountered in applications, e.g., in control theory and differential games. H–J equations are closely related to hyperbolic conservation laws—in one space dimension the former is simply the integrated version of the latter. Similarity also exists for the multidimensional case, and this is helpful in the design of difference approximations. In this paper high-order essentially nonoscillatory (ENO) schemes for H–J equations are investigated, which yield uniform high-order accuracy in smooth regions and sharply resolve discontinuities in the derivatives. The ENO scheme construction procedure is adapted from that for hyperbolic conservation laws. The schemes are numerically tested on a variety of one-dimensional and two-dimensional problems, including a problem related to control optimization, and high-order accuracy in smooth regions, good resolution of discontinuities in the derivatives, and convergence to viscosity solutions are observed.
- Published
- 1991
32. SUPERCONVERGENCE OF DISCONTINUOUS GALERKIN METHODS FOR TWO-DIMENSIONAL HYPERBOLIC EQUATIONS.
- Author
-
WAIXIANG CAO, CHI-WANG SHU, YANG YANG, and ZHIMIN ZHANG
- Subjects
- *
DISCONTINUOUS functions , *GALERKIN methods , *HYPERBOLIC differential equations , *INITIAL value problems , *BOUNDARY value problems , *DISCRETIZATION methods - Abstract
This paper is concerned with superconvergence properties of discontinuous Galerkin (DG) methods for two-dimensional linear hyperbolic conservation laws over rectangular meshes when upwind fluxes are used. We prove, under some suitable initial and boundary discretizations, the (2k + 1)th order superconvergence rate of the DG approximation at the downwind points and for the cell averages, when piecewise tensor-product polynomials of degree k are used. Moreover, we prove that the gradient of the DG solution is superconvergent with a rate of (k + 1)th order at all interior left Radau points; and the function value approximation is superconvergent at all right Radau points with a rate of (k+2)th order. Numerical experiments indicate that the aforementioned superconvergence rates are sharp. [ABSTRACT FROM AUTHOR]
- Published
- 2015
- Full Text
- View/download PDF
33. STABILITY AND ERROR ESTIMATES OF LOCAL DISCONTINUOUS GALERKIN METHODS WITH IMPLICIT-EXPLICIT TIME-MARCHING FOR ADVECTION-DIFFUSION PROBLEMS.
- Author
-
HAIJIN WANG, CHI-WANG SHU, and QIANG ZHANG
- Subjects
- *
ADVECTION-diffusion equations , *GALERKIN methods , *ERROR analysis in mathematics , *IMPLICIT functions , *RUNGE-Kutta formulas , *STABILITY theory - Abstract
The main purpose of this paper is to analyze the stability and error estimates of the local discontinuous Galerkin (LDG) methods coupled with carefully chosen implicit-explicit (IMEX) Runge-Kutta time discretization up to third order accuracy for solving one-dimensional linear advection-diffusion equations. In the time discretization the advection term is treated explicitly and the diffusion term implicitly. There are three highlights of this work. The first is that we establish an important relationship between the gradient and interface jump of the numerical solution with the independent numerical solution of the gradient in the LDG methods. The second is that, by aid of the aforementioned relationship and the energy method, we show that the IMEX LDG schemes are unconditionally stable for the linear problems in the sense that the time-step t is only required to be upper-bounded by a constant which depends on the ratio of the diffusion and the square of the advection coefficients and is independent of the mesh-size h, even though the advection term is treated explicitly. The last is that under this time step condition, we obtain optimal error estimates in both space and time for the third order IMEX Runge-Kutta time-marching coupled with LDG spatial discretization. Numerical experiments are also given to verify the main results. [ABSTRACT FROM AUTHOR]
- Published
- 2015
- Full Text
- View/download PDF
34. ANALYSIS OF OPTIMAL SUPERCONVERGENCE OF DISCONTINUOUS GALERKIN METHOD FOR LINEAR HYPERBOLIC EQUATIONS.
- Author
-
YANG YANG and CHI-WANG SHU
- Subjects
- *
STOCHASTIC convergence , *GALERKIN methods , *POLYNOMIALS , *DISCRETIZATION methods , *BOUNDARY value problems - Abstract
In this paper, we study the superconvergence of the error for the discontinuous Galerkin (DG) finite element method for linear conservation laws when upwind fluxes are used. We prove that if we apply piecewise kth degree polynomials, the error between the DG solution and the exact solution is (k + 2)th order superconvergent at the downwind-biased Radau points with suitable initial discretization. Moreover, we also prove the DG solution is (k + 2)th order superconvergent both for the cell averages and for the error to a particular projection of the exact solution. The superconvergence result in this paper leads to a new a posteriori error estimate. Our analysis is valid for arbitrary regular meshes and for Pk polynomials with arbitrary k ≥ 1, and for both periodic boundary conditions and for initial-boundary value problems. We perform numerical experiments to demonstrate that the superconvergence rate proved in this paper is optimal. [ABSTRACT FROM AUTHOR]
- Published
- 2012
- Full Text
- View/download PDF
35. SUPERCONVERGENCE OF DISCONTINUOUS GALERKIN METHODS FOR SCALAR NONLINEAR CONSERVATION LAWS IN ONE SPACE DIMENSION.
- Author
-
XIONG MENG, CHI-WANG SHU, QIANG ZHANG, and BOYING WU
- Subjects
- *
SUPERCONVERGENT methods , *GALERKIN methods , *DISCONTINUOUS functions , *POLYNOMIALS , *HYPERBOLIC differential equations - Abstract
In this paper, an analysis of the superconvergence property of the semidiscrete discontinuous Galerkin (DG) method applied to one-dimensional time-dependent nonlinear scalar conservation laws is carried out. We prove that the error between the DG solution and a particular projection of the exact solution achieves (k + 3/2)th order superconvergence when upwind fluxes are used. The results hold true for arbitrary nonuniform regular meshes and for piecewise polynomials of degree k (k ≥ 1), under the condition that |f'(u)| possesses a uniform positive lower bound. Numerical experiments are provided to show that the superconvergence property actually holds true for nonlinear conservation laws with general flux functions, indicating that the restriction on f(u)is artificial. [ABSTRACT FROM AUTHOR]
- Published
- 2012
- Full Text
- View/download PDF
36. A GENUINELY HIGH ORDER TOTAL VARIATION DIMINISHING SCHEME FOR ONE-DIMENSIONAL SCALAR CONSERVATION LAWS.
- Author
-
Xiangxiong Zhang and Chi-Wang Shu
- Subjects
- *
CONSERVATION laws (Mathematics) , *POLYNOMIALS , *FINITE differences , *FINITE element method , *NUMERICAL analysis - Abstract
It is well known that finite difference or finite volume total variation diminishing (TVD) schemes solving one-dimensional scalar conservation laws degenerate to first order accuracy at smooth extrema [S. Osher and S. Chakravarthy, SIAM J. Numer. Anal., 21 (1984), pp. 955- 984], thus TVD schemes are at most second order accurate in the L1 norm for general smooth and nonmonotone solutions. However, Sanders [Math. Comp., 51 (1988), pp. 535-558] introduced a third order accurate finite volume scheme which is TVD, where the total variation is defined by measuring the variation of the reconstructed polynomials rather than the traditional way of measuring the variation of the grid values. By adopting the definition of the total variation for the numerical solutions as in [R. Sanders, Math. Comp., 51 (1988), pp. 535-558], it is possible to design genuinely high order accurate TVD schemes. In this paper, we construct a finite volume scheme which is TVD in this sense with high order accuracy (up to sixth order) in the L1 norm. Numerical tests for a fifth order accurate TVD scheme will be reported, which include test cases from traffic flow models. [ABSTRACT FROM AUTHOR]
- Published
- 2010
- Full Text
- View/download PDF
37. ANALYSIS OF A LOCAL DISCONTINUOUS GALERKIN METHOD FOR LINEAR TIME-DEPENDENT FOURTH-ORDER PROBLEMS.
- Author
-
BO DONG and CHI-WANG SHU
- Subjects
- *
NUMERICAL analysis , *GALERKIN methods , *MATHEMATICAL analysis , *MULTIDIMENSIONAL scaling , *EXPERIMENTS , *EQUATIONS - Abstract
We analyze a local discontinuous Galerkin method for fourth-order time-dependent problems. Optimal error estimates are obtained in one dimension and in multidimensions for Cartesian and triangular meshes. We extend the analysis to higher even-order equations and the linearized Cahn-Hilliard type equations. Numerical experiments are displayed to verify the theoretical results. [ABSTRACT FROM AUTHOR]
- Published
- 2009
- Full Text
- View/download PDF
38. CONVERGENCE OF GODUNOV-TYPE SCHEMES FOR SCALAR CONSERVATION LAWS UNDER LARGE TIME STEPS.
- Author
-
Jing-Mei Qiu and Chi-Wang Shu
- Subjects
- *
MATHEMATICS education , *NUMERICAL solutions to conservation laws , *ACCELERATION of convergence in numerical analysis , *NUMERICAL analysis , *HYPERBOLIC differential equations , *QUANTUM entropy , *MATHEMATICAL analysis , *STOCHASTIC convergence , *CONVEX domains - Abstract
In this paper, we consider convergence of classical high-order Godunov-type schemes towards entropy solutions for scalar conservation laws. It is well known that sufficient conditions for such convergence include total variation boundedness of the reconstruction and cell or wavewise entropy inequalities. We prove that under large time steps, we only need total variation boundedness of the reconstruction to guarantee such convergence. We discuss high-order total variation bounded reconstructions to fulfill this sufficient condition and provide numerical examples on one-dimensional convex conservation laws to assess the performance of such large time step Godunov-type methods. To demonstrate the generality of this approach, we also prove convergence and give numerical examples for a large time step Godunov-like scheme involving Sanders' third-order total variation diminishing reconstruction using both cell averages and point values at cell boundaries. [ABSTRACT FROM AUTHOR]
- Published
- 2008
- Full Text
- View/download PDF
39. CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NONOSCILLATORY HIERARCHICAL RECONSTRUCTION.
- Author
-
Yingjie Liu, Chi-Wang Shu, Tadmor, Eitan, and Mengping Zhang
- Subjects
- *
GALERKIN methods , *NUMERICAL analysis , *DISCONTINUOUS functions , *CONSERVATION laws (Mathematics) , *RUNGE-Kutta formulas - Abstract
The central scheme of Nessyahu and Tadmor [J. Comput. Phys., 87 (1990), pp. 408-463] solves hyperbolic conservation laws on a staggered mesh and avoids solving Riemann problems across cell boundaries. To overcome the difficulty of excessive numerical dissipation for small time steps, the recent work of Kurganov and Tadmor [J. Comput. Phys., 160 (2000), pp. 241-282] employs a variable control volume, which in turn yields a semidiscrete nonstaggered central scheme. Another approach, which we advocate here, is to view the staggered meshes as a collection of overlapping cells and to realize the computed solution by its overlapping cell averages. This leads to a simple technique to avoid the excessive numerical dissipation for small time steps [Y. Liu, J. Comput. Phys., 209 (2005), pp. 82-104]. At the heart of the proposed approach is the evolution of two pieces of information per cell, instead of one cell average which characterizes all central and upwind Godunov-type finite volume schemes. Overlapping cells lend themselves to the development of a central-type discontinuous Galerkin (DG) method, following the series of works by Cockburn and Shu [J. Comput. Phys., 141 (1998), pp. 199-224] and the references therein. In this paper we develop a central DG technique for hyperbolic conservation laws, where we take advantage of the redundant representation of the solution on overlapping cells. The use of redundant overlapping cells opens new possibilities beyond those of Godunov-type schemes. In particular, the central DG is coupled with a novel reconstruction procedure which removes spurious oscillations in the presence of shocks. This reconstruction is motivated by the moments limiter of Biswas, Devine, and Flaherty [Appl. Numer. Math., 14 (1994), pp. 255-283] but is otherwise different in its hierarchical approach. The new hierarchical reconstruction involves a MUSCL or a second order ENO reconstruction in each stage of a multilayer reconstruction process without characteristic decomposition. It is compact, easy to implement over arbitrary meshes, and retains the overall preprocessed order of accuracy while effectively removing spurious oscillations around shocks. [ABSTRACT FROM AUTHOR]
- Published
- 2007
- Full Text
- View/download PDF
40. ERROR ESTIMATES TO SMOOTH SOLUTIONS OF RUNGE-KUTTA DISCONTINUOUS GALERKIN METHOD FOR SYMMETRIZABLE SYSTEMS OF CONSERVATION LAWS.
- Author
-
Qiang Zhang and Chi-Wang Shu
- Subjects
- *
NUMERICAL solutions to differential equations , *RUNGE-Kutta formulas , *GALERKIN methods , *CONSERVATION laws (Mathematics) , *FINITE element method , *NUMERICAL analysis - Abstract
In this paper we study the error estimates to sufficiently smooth solutions of symmetrizable systems of conservation laws for the Runge-Kutta discontinuous Galerkin (RKDG) method. Time discretization is the second-order explicit TVD (total variation diminishing) Runge-Kutta method, and the ℙκ (piecewise polynomial) finite element is used. When κ = 1 (piecewise linear finite element), the error estimate is obtained under the usual CFL condition τ ≤ βh for nonlinear systems in one dimension and for linear systems in multiple space dimensions. Here, h is the maximum element length, τ is the time step, and β is a positive constant independent of h and τ. Error estimates for ℙκ finite elements with κ > 1 are obtained under a more restrictive CFL condition. [ABSTRACT FROM AUTHOR]
- Published
- 2006
- Full Text
- View/download PDF
41. THE LOCAL DISCONTINUOUS GALERKIN METHOD FOR TIME-DEPENDENT CONVECTION-DIFFUSION SYSTEMS.
- Author
-
Cockburn, Bernardo and Chi-Wang Shu
- Subjects
- *
GALERKIN methods , *NUMERICAL analysis , *HYPERBOLIC differential equations , *NONLINEAR statistical models , *GEOMETRY , *MATHEMATICS - Abstract
In this paper, we study the local discontinuous Galerkin (LDG) methods for nonlinear, time-dependent convection-diffusion systems. These methods are an extension of the RungeKutta discontinuous Galerkin (RKDG) methods for purely hyperbolic systems to convection-diffusion systems and share with those methods their high parallelizability, high-order formal accuracy, and easy handling of complicated geometries for convection-dominated problems. It is proven that for scalar equations, the LDG methods are L2-stable in the nonlinear case. Moreover, in the linear case, it is shown that if polynomials of degree k are used, the methods are kth order accurate for general triangulations; although this order of convergence is suboptimal, it is sharp for the LDG methods. Preliminary numerical examples displaying the performance of the method are shown. [ABSTRACT FROM AUTHOR]
- Published
- 1998
- Full Text
- View/download PDF
Catalog
Discovery Service for Jio Institute Digital Library
For full access to our library's resources, please sign in.