Nothing Special   »   [go: up one dir, main page]

skip to main content
article

Multi-domain Fourier-continuation/WENO hybrid solver for conservation laws

Published: 01 October 2011 Publication History

Abstract

We introduce a multi-domain Fourier-continuation/WENO hybrid method (FC-WENO) that enables high-order and non-oscillatory solution of systems of nonlinear conservation laws, and which enjoys essentially dispersionless, spectral character away from discontinuities, as well as mild CFL constraints (comparable to those of finite difference methods). The hybrid scheme employs the expensive, shock-capturing WENO method in small regions containing discontinuities and the efficient FC method in the rest of the computational domain, yielding a highly effective overall scheme for applications with a mix of discontinuities and complex smooth structures. The smooth and discontinuous solution regions are distinguished using the multi-resolution procedure of Harten [J. Comput. Phys. 115 (1994) 319-338]. We consider WENO schemes of formal orders five and nine and a FC method of order five. The accuracy, stability and efficiency of the new hybrid method for conservation laws is investigated for problems with both smooth and non-smooth solutions. In the latter case, we solve the Euler equations for gas dynamics for the standard test case of a Mach three shock wave interacting with an entropy wave, as well as a shock wave (with Mach 1.25, three or six) interacting with a very small entropy wave and evaluate the efficiency of the hybrid FC-WENO method as compared to a purely WENO-based approach as well as alternative hybrid based techniques. We demonstrate considerable computational advantages of the new FC-based method, suggesting a potential of an order of magnitude acceleration over alternatives when extended to fully three-dimensional problems.

References

[1]
Albin, N. and Bruno, O.P., A spectral FC solver for the compressible Navier-Stokes equations in general domains I: explicit time-stepping. J. Comput. Phys. v230. 6248-6270.
[2]
Boyd, J.P., A comparison of numerical algorithms for Fourier extension of the first, second, and third kinds. J. Comput. Phys. v178. 118-160.
[3]
Boyd, J.P. and Ong, J.R., Exponentially-convergent strategies for defeating the Runge phenomenon for the approximation of non-periodic functions, Part I: single-interval schemes. Commun. Comput. Phys. v5. 484-497.
[4]
Boyd, J.P. and Ong, J.R., Exponentially-convergent strategies for defeating the Runge phenomenon for the approximation of non-periodic functions, part two: multi-interval polynomial schemes and multidomain Chebyshev interpolation. J. Appl. Numer. Math. v61. 460-472.
[5]
Bruno, O.P., Han, Y. and Pohlman, M.M., Accurate, high-order representation of complex three-dimensional surfaces via Fourier continuation analysis. J. Comput. Phys. v227. 1094-1125.
[6]
Bruno, O.P. and Lyon, M., High-order unconditionally stable FC-AD solvers for general smooth domains I. Basic elements. J. Comput. Phys. v229. 2009-2033.
[7]
Lyon, M. and Bruno, O.P., High-order unconditionally stable FC-AD solvers for general smooth domains II. Elliptic, parabolic and hyperbolic PDEs; theoretical considerations. J. Comput. Phys. v229. 3358-3381.
[8]
L.N. Trefethen, D. Bau III, Numerical Linear Algebra, Soc. for Industr. and Appl. Math., Philadelphia, 1997.
[9]
Harten, A., Adaptive multiresolution schemes for shock computations. J. Comput. Phys. v115. 319-338.
[10]
Costa, B. and Don, W.S., Multi-domain hybrid spectral-WENO methods for hyperbolic conservation laws. J. Comput. Phys. v224. 3358-3381.
[11]
Costa, B. and Don, W.S., High order hybrid central-WENO finite difference scheme for conservation laws. J. Comput. Appl. Math. v204. 209-218.
[12]
Costa, B., Don, W.S., Gottlieb, D. and Sendersky, R., Two-dimensional multi-domain hybrid spectral-WENO methods for conservation laws. Cummun. Comput. Phys. v1. 550-577.
[13]
Jiang, G.S. and Shu, C.W., Efficient implementation of weighted ENO schemes. J. Comput. Phys. v126. 202-228.
[14]
Shu, C.-W., Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws. In: Quarteroni, A. (Ed.), Lecture Notes in Mathematics, vol. 1697. Springer, Cetraro, Italy. pp. 325
[15]
Shu, C.W., High order weighted essentially nonoscillatory schemes for convection dominated problems. SIAM Rev. v51. 82-126.
[16]
Balsara, D.S. and Shu, C.W., Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy. J. Comput. Phys. v160. 405-452.
[17]
Shu, C.-W. and Osher, S., Efficient implementation of essentially non-oscillatory shock capturing schemes. J. Comput. Phys. v77. 439-471.
[18]
Liu, X-D., Osher, S. and Chan, T., Weighted essentially non-oscillatory schemes. J. Comput. Phys. v115. 200-212.
[19]
Borges, R., Carmona a, M., Costa, B. and Don, W.S., An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws. J. Comput. Phys. v227. 3191-3211.
[20]
Gottlieb, D. and Hesthaven, J.S., Spectral methods for hyperbolic problems. J. Comput. Appl. Math. v128. 83-131.
[21]
Hesthaven, J.S., Gottlieb, S. and Gottlieb, D., . 2007. Cambridge Monographs on Applied and Computational Mathematics 21, 2007.Cambridge University Press, Cambridge, UK.
[22]
Morkovin, M.K., Note on the assessment of flow disturbances at a blunt body traveling at supersonic speeds owing to flow disturbance in free stream. J. Appl. Mech. v27. 223-229.
[23]
Shi, J., Zhang, Y. and Shu, C.W., Resolution of high order WENO schemes for complicated flow structures. J. Comput. Phys. v186. 690-696.
[24]
Latini, M., Schilling, O. and Don, W.S., Effects of WENO flux reconstruction order and spatial resolution on reshocked two-dimensional Richtmyer-Meshkov instability. J. Comput. Phys. v221. 805-836.
[25]
Schilling, O. and Latini, M., High-order WENO simulations of three-dimensional reshocked Richtmyer-Meshkov instability to late times: dynamics, dependence on initial conditions, and comparisons to experimental data. Acta Math. Sci. v30. 595-620.
[26]
Babuska, I.A. and Sauter, S.A., Is the pollution the FEM avoidable for the Helmholtz equation considering high wave numbers?. SIAM J. Numer. Anal. v34. 2392-2423.
[27]
Jameson, L., High order schemes for resolving waves: number of points per wavelength. J. Sci. Comput. v15. 417-433.
[28]
Gottlieb, S. and Shu, C.W., Total variation diminishing Runge-Kutta schemes. Math. Comp. v67. 73-85.
[29]
Dahlquist, G. and Bjork, A., Numerical Methods. 1974. Prentice Hall, Englewood Clis.
[30]
Gustafsson, B., High Order Difference Methods for Time Dependent PDE. 2008. Springer.
[31]
Canuto, C.G., Hussaini, M.Y., Quarteroni, A. and Zang, T.A., Spectral Methods: Evolution to Complex Geometries and Applications to Fluid Dynamics. 2007. Springer.
[32]
Kopriva, D.A. and Kolias, J.H., A conservative staggered-grid Chebyshev multidomain method for compressible flows. J. Comput. Phys. v125. 244-261.
[33]
Funaro, D. and Gottlieb, D., A new method of imposing boundary conditions in pseudospectral approximations of hyperbolic equations. Math. Comput. v51. 599-613.
[34]
Gottlieb, D., Hussaini, M.Y. and Orszag, S.A., Theory and application of spectral methods. In: Voigt, R.G., Gottlieb, D., Hussaini, M.Y. (Eds.), Spectral Methods for Partial Differential Equations, SIAM, Philadelphia.
[35]
Castro, M., Costa, B. and Don, W.S., High order weighted essentially non-oscillatory WENO-Z schemes for hyperbolic conservation laws. J. Comput. Physic. v230. 1766-1792.

Cited By

View all
  • (2018)A Domain Decomposition Fourier Continuation Method for Enhanced $$L_1$$L1 Regularization Using Sparsity of Edges in Reconstructing Fourier DataJournal of Scientific Computing10.1007/s10915-017-0467-y74:2(851-871)Online publication date: 1-Feb-2018
  • (2017)Robust second-order scheme for multi-phase flow computationsJournal of Computational Physics10.1016/j.jcp.2017.03.025339:C(163-178)Online publication date: 15-Jun-2017
  • (2017)A numerical study of the local monotone polynomial edge detection for the hybrid WENO methodJournal of Computational and Applied Mathematics10.1016/j.cam.2017.02.029321:C(232-245)Online publication date: 1-Sep-2017
  • Show More Cited By
  1. Multi-domain Fourier-continuation/WENO hybrid solver for conservation laws

      Recommendations

      Comments

      Please enable JavaScript to view thecomments powered by Disqus.

      Information & Contributors

      Information

      Published In

      cover image Journal of Computational Physics
      Journal of Computational Physics  Volume 230, Issue 24
      October, 2011
      189 pages

      Publisher

      Academic Press Professional, Inc.

      United States

      Publication History

      Published: 01 October 2011

      Author Tags

      1. Conservation laws
      2. Fourier continuation methods
      3. High-order WENO methods
      4. Multi-resolution methods
      5. Shock waves

      Qualifiers

      • Article

      Contributors

      Other Metrics

      Bibliometrics & Citations

      Bibliometrics

      Article Metrics

      • Downloads (Last 12 months)0
      • Downloads (Last 6 weeks)0
      Reflects downloads up to 19 Dec 2024

      Other Metrics

      Citations

      Cited By

      View all
      • (2018)A Domain Decomposition Fourier Continuation Method for Enhanced $$L_1$$L1 Regularization Using Sparsity of Edges in Reconstructing Fourier DataJournal of Scientific Computing10.1007/s10915-017-0467-y74:2(851-871)Online publication date: 1-Feb-2018
      • (2017)Robust second-order scheme for multi-phase flow computationsJournal of Computational Physics10.1016/j.jcp.2017.03.025339:C(163-178)Online publication date: 15-Jun-2017
      • (2017)A numerical study of the local monotone polynomial edge detection for the hybrid WENO methodJournal of Computational and Applied Mathematics10.1016/j.cam.2017.02.029321:C(232-245)Online publication date: 1-Sep-2017
      • (2016)An FC-based spectral solver for elastodynamic problems in general three-dimensional domainsJournal of Computational Physics10.1016/j.jcp.2015.11.060307:C(333-354)Online publication date: 15-Feb-2016
      • (2015)Hybrid Fourier-Continuation Method and Weighted Essentially Non-oscillatory Finite Difference Scheme for Hyperbolic Conservation Laws in a Single-Domain FrameworkJournal of Scientific Computing10.1007/s10915-014-9913-264:3(670-695)Online publication date: 1-Sep-2015
      • (2014)The Fourier approximation of smooth but non-periodic functions from unevenly spaced dataAdvances in Computational Mathematics10.1007/s10444-014-9342-740:5-6(1073-1092)Online publication date: 1-Dec-2014
      • (2013)Multi-dimensional hybrid Fourier continuation-WENO solvers for conservation lawsJournal of Computational Physics10.5555/2743135.2743305253:C(209-225)Online publication date: 15-Nov-2013

      View Options

      View options

      Media

      Figures

      Other

      Tables

      Share

      Share

      Share this Publication link

      Share on social media