Hybrid Fourier-Continuation Method and Weighted Essentially Non-oscillatory Finite Difference Scheme for Hyperbolic Conservation Laws in a Single-Domain Framework

Peng Li1, Zhen Gao2, Wai-Sun Don3,2, Shusen Xie2
1State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing, China
2School of Mathematical Sciences, Ocean University of China, Qingdao, China
3Division of Applied Mathematics, Brown University, Providence, USA

Tóm tắt

We investigate a hybrid Fourier-Continuation (FC) method (Bruno and Lyon, J Comput Phys 229:2009–2033, 2010) and fifth order characteristic-wise weighted essentially non-oscillatory (WENO) finite difference scheme for solving system of hyperbolic conservation laws on a uniformly discretized Cartesian domain. The smoothness of the solution is measured by the high order multi-resolution algorithm by Harten (J Comput Phys 49:357–393, 1983) at each grid point in a single-domain framework (Costa and Don, J Comput Appl Math 204(2):209–218, 2007) (Hybrid), as opposed to each subdomain in a multi-domain framework (Costa et al., J Comput Phys 224(2):970–991, 2007; Shahbazi et al., J Comput Phys 230:8779–8796, 2011). The Hybrid scheme conjugates a high order shock-capturing WENO-Z5 (nonlinear) scheme (Borges et al., J Comput Phys 227:3101–3211, 2008) in non-smooth WENO stencils with an essentially non-dissipative and non-dispersive FC (linear) method in smooth FC stencils, yielding a high fidelity scheme for applications containing both discontinuous and complex smooth structures. Several critical and unique numerical issues in an accurate and efficient implementation (such as reasonable choice of parameters, singular value decomposition, fast Fourier transform, symmetry preservation, and overlap zone) of the FC method, due to a dynamic spatial and temporal change in the size of data length in smooth FC stencils in a single-domain framework, will be illustrated and addressed. The accuracy and efficiency of the Hybrid scheme in solving one and two dimensional system of hyperbolic conservation laws is demonstrated with several classical examples of shocked flow, such as the one dimensional Riemann initial value problems (123, Sod and Lax), the Mach 3 shock–entropy wave interaction problem with a small entropy sinusoidal perturbation, the Mach 3 shock–density wave interaction problem, and the two dimensional Mach 10 double Mach reflection problem. For a sufficiently large problem size, a factor of almost two has been observed in the speedup of the Hybrid scheme over the WENO-Z5 scheme.

Từ khóa


Tài liệu tham khảo

Adams, N., Shariff, K.: High-resolution hybrid compact-ENO scheme for shock-turbulence interaction problems. J. Comput. Phys. 127, 27–51 (1996) Albin, N., Bruno, O.P.: A spectral FC solver for the compressible Navier–Stokes equations in general domains I: explicit time-stepping. J. Comput. Phys. 230(16), 2009–2033 (2010) Balsara, D., Shu, C.W.: Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy. J. Comput. Phys. 160, 405–452 (2000) Borges, R., Carmona, M., Costa, B., Don, W.S.: An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws. J. Comput. Phys. 227, 3101–3211 (2008) Bruno, O.P., Lyon, M.: High-order unconditionally stable FC-AD solvers for general smooth domains I. Basic elements. J. Comput. Phys. 229, 2009–2033 (2010) Castro, M., Costa, B., Don, W.S.: High order weighted essentially non-oscillatory WENO-Z schemes for hyperbolic conservation laws. J. Comput. Phys. 230, 1766–1792 (2011) Costa, B., Don, W.S.: On the computation of high order pseudospectral derivatives. Applied Numer. Math. 33, 151–159 (2000) Costa, B., Don, W.S.: High order Hybrid Central-WENO finite difference scheme for conservation laws. J. Comput. Appl. Math. 204(2), 209–218 (2007) Costa, B., Don, W.S.: Multi-domain hybrid spectral-WENO methods for hyperbolic conservation laws. J. Comput. Phys. 224(2), 970–991 (2007) Costa, B., Don, W.S., Gottlieb, D., Sendersky, R.: Two-dimensional multi-domain hybrid spectral-WENO methods for conservation laws. Commun. Comput. Phys. 1(3), 548–574 (2006) Dahlquist, G., Björk, A.: Numerical Methods. Prentice Hall, Englewood Clis (1974) Don, W.S., Borges, R.: Accuracy of the weighted essentially non-oscillatory conservative finite difference schemes. J. Comput. Phys. 250, 347–372 (2013) Don, W.S.: HOPEpack. http://www.math.hkbu.edu.hk/~wsdon/ Gottlieb, D., Orzag, S.: Numerical analysis of spectral methods: theory and applications, SIAM (1977). Gao, Z., Don, W.S.: Mapped hybrid central-WENO finite difference scheme for detonation waves simulations. J. Sci. Comput. 55, 351–371 (2013) Harten, A.: High resolution schemes for hyperbolic conservation laws. J. Comput. Phys. 49, 357–393 (1983) Henrick, A.K., Aslam, T.D., Powers, J.M.: Mapped weighted essentially non-oscillatory schemes: achieving optimal order near critical points. J. Comput. Phys. 207, 542–567 (2005) Hill, D., Pullin, D.: Hybrid tuned center-difference-WENO method for large eddy simulations in the presence of strong shocks. J. Comput. Phys. 194, 435–450 (2004) Isaacson, E., Keller, H.B.: Analysis of Numerical Methods. Wiley, New York (1966) Jacobs, G.B., Don, W.S.: A high-order WENO-Z finite difference based particle-source-in-cell method for computation of particle-laden flows with shocks. J. Comput. Phys. 228, 1365–1379 (2008) Jacobs, G.B., Don, W.S., Dittmann, T.: High-order resolution eulerian-lagrangian simulations of particle dispersion in the accelerated flow behind a moving shock. Theor. Comput. Fluid Dyn. 26, 37–50 (2012) Jiang, G.S., Shu, C.W.: Efficient implementation of weighted ENO schemes. J. Comput. Phys. 126, 202–228 (1996) Loan, C.V.: Computational frameworks for the fast fourier transform, SIAM (1992) Lyon, M., 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. 229, 3358–3381 (2010) Pirozzoli, S.: Conservative hybrid compact-WENO schemes for shock–turbulence interaction. J. Comput. Phys. 178(1), 81–117 (2002) Shahbazi, K., Albin, N., Bruno, O., Hesthaven, J.: Multi-domain fourier-continuation/WENO hybrid solver for conservation laws. J. Comput. Phys. 230, 8779–8796 (2011) Shahbazi, K., Hesthaven, J., Zhu, X.: Multi-dimensional hybrid Fourier continuation-CWENO solvers for conservation laws. J. Comput. Phys. 253, 209–225 (2013) Shu, C.W.: High order weighted essentially nonoscillatory schemes for convection dominated problems. SIAM Rev. 51(1), 82–126 (2009) Woodward, P., Collela, P.: The numerical simulation of two dimensional fluid flow with strong shocks. J. Comput. Phys. 54, 115–173 (1984)