A High-Order Dirac-Delta Regularization with Optimal Scaling in the Spectral Solution of One-Dimensional Singular Hyperbolic Conservation Laws
Pages A1831 - A1849
Abstract
A regularization technique based on a class of high-order, compactly supported piecewise polynomials is developed that regularizes the time-dependent, singular Dirac-delta sources in spectral approximations of hyperbolic conservation laws. The regularization technique provides higher-order accuracy away from the singularity. A theoretical criterion that establishes a lower bound on the support length (optimal scaling) has to be satisfied to achieve optimal order of convergence. The optimal scaling parameter has been shown to be, instead of a fixed constant value, a function of the smoothness of the compactly supported piecewise polynomial and the number of vanishing moment of the polynomial when integrated with respect to mononomial of degree up to some order (moment). The effectiveness of the criterion is illustrated in the solutions of a linear and a nonlinear (Burgers) scalar hyperbolic conservation law with a singular source, as well as the nonlinear Euler equations with singular sources, a system of hyperbolic conservation laws governing compressible fluid dynamics with shocks and particles. A Chebyshev collocation method (spectral) discretizes the spatial derivatives in the scalar equation tests. A multidomain hybrid spectral-WENO method discretizes the Euler equations. The WENO scheme captures shocks and high gradients in the gas flow, whereas the spectral discretization applies in all other subdomains including the regularization of the singular source. In the advection and Burgers problems, the convergence order of the numerical solution follows the asymptotic behavior of the singular source approximation. The multidomain hybrid scheme is in excellent agreement with the computations that are based on a WENO discretization only.
References
[1]
H. Abe, S. Natsuhiko, and R. Itatani, High-order spline interpolations in the particle simulation, J. Comput. Phys., 63 (1986), pp. 247--267.
[2]
C. K. Birdsall and A. B. Langdon, Plasma Physics via Computer Simulation, McGraw-Hill, New York, 1985.
[3]
R. Borges, M. Carmona, B. Costa, and W. S. Don, An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws, J. Comput. Phys., 227 (2008), pp. 3101--3211.
[4]
A. Bressan, Hyperbolic Systems of Conservation Laws, Oxford University Press, New York, 2000.
[5]
M. Castro, B. Costa, and W. S. Don, High order weighted essentially non-oscillatory WENO-Z schemes for hyperbolic conservation laws, J. Comput. Phys., 230 (2011), pp. 1766--1792.
[6]
B. Costa and W. S. Don, Multi-domain hybrid spectral-WENO methods for hyperbolic conservation laws, J. Comput. Phys., 224 (2007), pp. 970--991.
[7]
B. Costa, W. S. Don, D. Gottlieb, and R. Sendersky, Two-dimensional multi-domain hybrid spectral-WENO methods for conservation laws, Commun. Comput. Phys., 1 (2006), pp. 550--577.
[8]
S. L. Davis, T. B. Dittman, G. Jacobs, and W. S. Don, Disperson of a cloud of particles by a moving shock: Effects of shape, angle of rotation and aspect ratio, J. Appl. Mech. Tech. Phys., 54 (2013), pp. 900--912.
[9]
W. S. Don, A. De Gregorio, J. P. Suarez, and G. Jacobs, Assessing the performance of a three dimensional hybrid central-WENO finite difference scheme with computation of a sonic injector in supersonic cross flow, Adv. Appl. Math. Mech., 4 (2012), pp. 719--736.
[10]
B. Engquist, A. K. Tornberg, and R. Tsai, Discretization of Dirac delta functions in level set methods, J. Comput. Phys., 207 (2005), pp. 28--51.
[11]
E. Godlewski and P. A. Raviart, Numerical Approximation of Hyperbolic Systems of Conservation Laws, Springer-Verlag, New York, 1996.
[12]
S. Gottlieb and C.-W. Shu, Total variation diminishing Runge--Kutta schemes, Math. Comp., 67 (1998), pp. 73--85.
[13]
A. Harten, High resolution schemes for hyperbolic conservation laws, J. Comput. Phys., 49 (1983), pp. 357--393.
[14]
A. Harten, Adaptative multiresolution schemes for shock computations, J. Comput. Phys., 115 (1994), pp. 319--338.
[15]
J. S. Hesthaven, S. Gottlieb, and D. Gottlieb, Spectral Methods for Time-Dependent Problems, Vol. 21, Cambridge University Press, Cambridge, 2007.
[16]
R. W. Hockney and J. W. Eastwood, Computer Simulation Using Particles, Taylor & Francis, New York, 1988.
[17]
G. Jacobs and W. S. Don, A high-order WENO-Z finite difference scheme based particle-source-in-cell method for computation of particle-laden flows with shocks, J. Comput. Phys., 228 (2009), pp. 1365--1379.
[18]
G. Jacobs and J. S. Hesthaven, High-order nodal discontinuous Galerkin particle-in-cell method on unstructured grids, J. Comput. Phys., 214 (2006), pp. 96--121.
[19]
J. H. Jung, A note on the spectral collocation approximation of some differential equations with singular source terms, J. Sci. Comput., 39 (2009), pp. 49--66.
[20]
J. H. Jung and W. S. Don, Collocation methods for hyperbolic partial differential equations with singular sources, Adv. Appl. Math. Mech., 1 (2009), pp. 769--780.
[21]
R. P. Kanwal, Generalized Functions: Theory and Technique, 2nd ed., Birkhäuser Boston, Cambridge, MA, 1998.
[22]
R. J. Leveque, Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, Cambridge, 2002.
[23]
J. J. Monaghan, Shock simulation by the particle method SPH, J. Comput. Phys., 52 (1983), pp. 374--389.
[24]
J. J. Monaghan, Extrapolating B splines for interpolation, J. Comput. Phys., 60 (1985), pp. 253--262.
[25]
P. J. Montgomery and T. B. Moodie, Jump conditions for hyperbolic systems of forced conservation laws with an application to gravity currents, Stud. Appl. Math., 106 (2001).
[26]
A. Quarteroni, R. Sacco, and F. Saleri, Numerical Mathematics, Springer-Verlag, New York, 2000.
[27]
J. Ryan and C.-W. Shu, On a one-sided post-processing technique for the discontinuous Galerkin methods, Methods Appl. Anal., 10 (2003), pp. 295--308.
[28]
C.-W. Shu, Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, Lecture Notes in Math. 1697, 1998, pp. 325--432.
[29]
A. K. Tornberg, Multi-dimensional quadrature of singular and discontinuous functions, BIT, 42 (2002), pp. 644--669.
[30]
L. N. Trefethen, Is Gauss quadrature better than Clenshaw--Curtis?, SIAM Rev., 50 (2008), pp. 67--87.
[31]
Y. Yang and C.-W. Shu, Discontinuous Galerkin method for hyperbolic equations involving $\delta$-singularities: Negative-order norm error estimates and applications, Numer. Math., 124 (2013), pp. 753--781.
Index Terms
- A High-Order Dirac-Delta Regularization with Optimal Scaling in the Spectral Solution of One-Dimensional Singular Hyperbolic Conservation Laws
Index terms have been assigned to the content through auto-classification.
Recommendations
Regularization of Singularities in the Weighted Summation of Dirac-Delta Functions for the Spectral Solution of Hyperbolic Conservation Laws
Singular source terms expressed as weighted summations of Dirac-delta functions are regularized through approximation theory with convolution operators. We consider the numerical solution of scalar and one-dimensional hyperbolic conservation laws with ...
A Moving Mesh Method for One-dimensional Hyperbolic Conservation Laws
We develop an adaptive method for solving one-dimensional systems of hyperbolic conservation laws that employs a high resolution Godunov-type scheme for the physical equations, in conjunction with a moving mesh PDE governing the motion of the spatial grid ...
On maximum-principle-satisfying high order schemes for scalar conservation laws
We construct uniformly high order accurate schemes satisfying a strict maximum principle for scalar conservation laws. A general framework (for arbitrary order of accuracy) is established to construct a limiter for finite volume schemes (e.g. essentially ...
Comments
Please enable JavaScript to view thecomments powered by Disqus.Information & Contributors
Information
Published In
© 2014, Society for Industrial and Applied Mathematics.
Publisher
Society for Industrial and Applied Mathematics
United States
Publication History
Published: 01 January 2014
Author Tags
Author Tags
Qualifiers
- Research-article
Contributors
Other Metrics
Bibliometrics & Citations
Bibliometrics
Article Metrics
- 0Total Citations
- 0Total Downloads
- Downloads (Last 12 months)0
- Downloads (Last 6 weeks)0
Reflects downloads up to 12 Jan 2025