Hostname: page-component-78c5997874-ndw9j Total loading time: 0 Render date: 2024-11-10T09:04:06.675Z Has data issue: false hasContentIssue false

Development of a Combined Compact Difference Scheme to Simulate Soliton Collision in a Shallow Water Equation

Published online by Cambridge University Press:  16 March 2016

Ching-Hao Yu
Affiliation:
Ocean College, Zhejiang University, 866 Yuhangtang Road, Hangzhou, Zhejiang, China
Tony Wen-Hann Sheu*
Affiliation:
Department of Engineering Science and Ocean Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei, Taiwan Institute of Applied Mathematical Sciences, National Taiwan University, Taipei, Taiwan Center of Advanced Study in Theoretical Sciences (CASTS), National Taiwan University, Taipei, Taiwan
*
*Corresponding author. Email addresses: chyu@zju.edu.cn (C.-H. Yu), twhsheu@ntu.edu.tw (T. W.-H. Sheu)
Get access

Abstract

In this paper a three-step scheme is applied to solve the Camassa-Holm (CH) shallow water equation. The differential order of the CH equation has been reduced in order to facilitate development of numerical scheme in a comparatively smaller grid stencil. Here a three-point seventh-order spatially accurate upwinding combined compact difference (CCD) scheme is proposed to approximate the firstorder derivative term. We conduct modified equation analysis on the CCD scheme and eliminate the leading discretization error terms for accurately predicting unidirectional wave propagation. The Fourier analysis is carried out as well on the proposed numerical scheme to minimize the dispersive error. For preserving Hamiltonians in Camassa- Holm equation, a symplecticity conserving time integrator has been employed. The other main emphasis of the present study is the use of uPα formulation to get nondissipative CH solution for peakon-antipeakon and soliton-anticuspon head-on wave collision problems.

Type
Research Article
Copyright
Copyright © Global-Science Press 2016 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Xu, Y. and Shu, C.W.. Local discontinuous Galerkin methods for high-order time-dependent partial differential equation. Commun. Comput. Phys., 7:146, 2010.Google Scholar
[2]Camassa, R. and Holm, D.. An integrable shallow water equation with peaked solitons. Phys. Rev. Letter, 7:16611664, 1993.CrossRefGoogle Scholar
[3]Constantin, A.. Existence of permanent and breaking waves for a shallow water equation: a geometric approach. Ann. Inst. Fourier, Grenoble, 50:321362, 2000.CrossRefGoogle Scholar
[4]Parker, A.. Binary cuspon-soliton interactions. Chaos Soliton Fract., 41:15311549, 2009.CrossRefGoogle Scholar
[5]Xu, Y. and Shu, C.W.. A local discontinuous Galerkin method for the Camassa-Holm equation. SIAM J. Numer. Anal., 46(4):19982021, 2008.CrossRefGoogle Scholar
[6]Liang, X., Li, X., Fu, D., Ma, Y.. Complex transition of double-diffusive convection in a rectangular enclosure with hieght-to-length ratio equal to 4: Part I. Commun. Comput. Phys., 6(2):247268, 2009.CrossRefGoogle Scholar
[7]Zhang, Y. T. andShu, C. W.. Third order WENO schemes on three dimensional tetrahedral meshes. Commun. Comput. Phys., 5:863–848, 2009.Google Scholar
[8]Feng, B. F., Maruno, K. and Ohta, Y.. A self-adaptive mesh method for the Camassa-Holm equation. J. Comput. Appl. Math., 235(1):229424, 2010.CrossRefGoogle Scholar
[9]Kalisch, H. and Lenells, J..Numerical study of traveling-wave solutions for the Camassa-Holm equation. Chaos Soliton Fract., 25:287298, 2005.CrossRefGoogle Scholar
[10]Artebrant, R. and Schroll, H. J.. Numerical simulation of Camassa-Holm peakons by adaptive upwinding. Appl. Numer. Math., 56:695711, 2006.CrossRefGoogle Scholar
[11]Cohen, D., Owren, B. and Raynaud, X.. Multi-symplectic integration of the Camassa-Holm equation. J. Comput. Phys., 227:54925512, 2008.CrossRefGoogle Scholar
[12]Matsuo, T. and Yamaguchi, H.. An energy-conserving Galerkin scheme for a class of nonlinear dispersive equations. J. Comput. Phys., 228:43464358, 2009.CrossRefGoogle Scholar
[13]Matsuo, T.. A Hamiltonian-conserving Galerkin scheme for the Camassa Holm equation.. J. Comput. Appl. Math., 234:12581266, 2010.CrossRefGoogle Scholar
[14]Holm, D. D. and Staley, M. F..Wave structure and nonlinear balances in a family of evolutionary PDEs. SIAM J. Appl. Dyn. Syst., 2:323380, 2003.CrossRefGoogle Scholar
[15]Xu, Y. and Shu, C.W.. Local discontinuous Galerkin methods for the Degasperis-Procesi equation. Commun. Comput. Phys., 10:474508, 2011.CrossRefGoogle Scholar
[16]Feng, B.F. and Liu, Y.. An operator splitting method for the Degasperis-Procesi equation. J. Comput. Phys., 228:78057820, 2009.CrossRefGoogle Scholar
[17]Bridges, T. J. and Reich, S.. Multi-symplectic integrators: numerical schemes for Hamiltonian PDE that conserve symplecticity. Phys. Letter A, 284(4-5):184193, 2001.CrossRefGoogle Scholar
[18]Stanislavova, M. and Stefanov, A.. On global finite energy solutions of the Camassa-Holm equation. J. Fourier Anal. Appl., 11(5):511–53, 2005.CrossRefGoogle Scholar
[19]Chiu, P. H., Lee, L. and Sheu, T. W. H.. A dispersion-relation-preserving algorithm for a non-linear shallow-water wave equation. J. Comput. Phys., 228:80348052, 2009.CrossRefGoogle Scholar
[20]Sheu, T. W. H., Chiu, P. H. and Yu, C. H.. On the development of a high-order compact scheme for exhibiting the switching and dissipative solution natures in the Camassa-Holm equation. J. Comput. Phys., 230:53995416, 2011.CrossRefGoogle Scholar
[21]Oevel, W. and Sofroniou, W.. Symplectic Runge-Kutta schemes II: classification of symplectic method. Univ. of Paderborn, Germany, Preprint, 1997.Google Scholar
[22]Sengupta, T. K., Rajpoot, M. K. and Bhumkar, Y. G.. Space-time discretizing optimal DRP schemes for flow and wave propagation problems. Comput. Fluids, 47(1):144154, 2011.CrossRefGoogle Scholar
[23]Sengupta, T. K., Vijay, V. V. S. N. and Bhaumik, S.. Further improvement and analysis of CCD scheme: Dissipation discretization and de-aliasing properties. J. Comput. Phys., 228:61506168, 2009.CrossRefGoogle Scholar
[24]Nihei, T. and Ishii, K.. A fast solver of the shallow water equations on a sphere using a combined compact difference scheme. J. Comput. Phys., 187:639659, 2003.CrossRefGoogle Scholar
[25]Tam, C. K. W. and Webb, J. C.. Dispersion-relation-preserving finite difference schemes for computational acoustics. J. Comput. Phys., 107:262281, 1993.CrossRefGoogle Scholar
[26]Chu, P. C. and Fan, C.. A three-point combined compact difference scheme. J. Comput. Phys., 140:370399, 1998.CrossRefGoogle Scholar
[27]Parker, A. and Matsuno, Y.. The peakon limits of soliton solutions of the Camassa-Holm equation. J. Phys. Soc. Japan, 75(12):124001, 2006.CrossRefGoogle Scholar
[28]Lenells, J.. Traveling wave solutions of the Camassa-Holm eqaution. J. Differen. Equat., 217:393430, 2005.CrossRefGoogle Scholar
[29]Camassa, R., Holm, D. D. and Hyman, J. M.. A new integrable shallow water equation. Adv. Appl. Mech., 31:133, 1994.CrossRefGoogle Scholar
[30]Holden, H. and Raynaud, X.. Global dissipative multipeakon solutions of the Camassa-Holm equation. Comm. Partial Differen. Equat., 33(11):20402063, 2008.CrossRefGoogle Scholar
[31]Jiang, G.S. and Peng, D.. Weighted ENO schemes for Hamilton-Jacobi equations. SIAM J. Sci. Comput., 21:21262143, 2000.CrossRefGoogle Scholar
[32]Ohta, Y., Maruno, K I. and Feng, B. F.. An integrable semi-discretization of the Camassa-Holm equation and its determinant solution. J. Phys. A: Math. Theor., 41:355205, 2008.CrossRefGoogle Scholar