Abstract
In this paper, the Mohand transform-based homotopy perturbation method is proposed to solve two-dimensional linear and non-linear shallow water wave equations. This approach has been proved suitable for a broad variety of non-linear differential equations in science and engineering. The variation trend of the water surface elevation at different time levels and depths are given by some graphs. Moreover, the obtained solutions are compared with the existing results, which show higher efficiency and fewer computations than other approaches studied in the literature.
1 Introduction
Water waves that have a horizontal scale much larger than the depth of the fluid are considered shallow water waves (SWWs). SWWs describe the evolution of incompressible flow, neglecting density change along the depth, which are widely used to simulate the propagation of tsunami waves, tidal currents, storm floods, and shock waves [–]. Karunakar and Chakraverty applied the homotopy perturbation method (HPM) to solve linear and non-linear one-dimensional shallow water wave equations (SWWEs) [], while the two-dimensional equations were solved in []. Sahoo and Chakraverty [] used the Sawi transform-based homotopy perturbation method to solve the one-dimensional SWWEs. The application of the variational iteration method and homotopy perturbation method (HPM) was given in [] for solving one-dimensional shallow water equations with crisp and fuzzy uncertain initial conditions. Noeiaghdam and Sidorov [] used the homotopy analysis method to solve the non-linear shallow water wave equation. Safari [] gave solutions to two extended model equations for shallow water waves by He’s variational iteration method. In [], the Adomian decomposition method was used to solve two extended model equations for SWWEs. In addition, Bekir and Aksoy [] applied the Exp-function method to construct exact solutions to non-linear wave equations.
Mohand transform (MT) [] is an integral transformation, which is obtained from the classical Fourier integral, and it was initiated by Mohand Mahgoub to solve some integral and differential equations (–). Khandelwal [] combined the Adomian decomposition method and Mohand transform for solving a differential equation of mixing layers that arise in a viscous incompressible form. The comparison of two integral transforms between Mohand and Laplace is given in []; the results show that the two transforms are closely related. Aggarwal, Sharma, and Chauhan [] gave the duality relations of Kamal transform with Laplace, Laplace−Carson, Aboodh, Sumudu, Elzaki, Mohand, and Sawi transforms []. They used an effective scheme based on Mohand transform and the homotopy perturbation method to find numerical solutions to non-linear fractional shock wave equations. Althobaiti, Dubey, and Prasad [] presented the local fractional Mohand transform with the Adomian decomposition method to solve the local fractional generalized Fokker–Planck equation. Shah, Khan, and Farooq et al. [] employed the Mohand transform to provide the analytical solution to the one-dimensional time fractional system of PDEs.
The homotopy perturbation method (HPM) [] has been widely used for solving differential and integral equations, which was first proposed and developed by the Chinese mathematician Ji-Huan He [–]. Laplace transformation coupled with the homotopy perturbation method is called the He−Laplace method. Madani and Mishra [, ] adopted this method to solve non-linear ordinary and partial differential equations. It is proved that this method needs fewer computations and is much easier than other methods. He, Moatimid, and Zekry [] used the two-scale transforms and He−Laplace method to get an analytic approximate solution to a fractal−differential model []. They presented the combination of the enhanced perturbation method and the parameter expansion technology based on Li−He’s approach, aiming to find an approximate solution to the given second-order non-linear ordinary differential equation []. They applied the basic idea of the perturbation method to construct a homotopy equation with a higher order; the results show a highly effective obtained solution. The approximation of the damped non-linear Klein–Gordon equation was given [] by the coupling of the exponential decay parameter with the homotopy perturbation method, which overcomes the shortcomings of the classical method []. They used the reducing rank method with the homotopy perturbation method to solve the third-order Duffing equation.
The main purpose of the present work is to solve two-dimensional linear and non-linear shallow water wave equations by coupling Mohand transform with the homotopy perturbation method. The remainder of the paper is organized as follows: in Section 2, we present some preliminaries about shallow water wave equations (SWWEs), Mohand transform (MT), and the homotopy perturbation method (HPM). In Section 3, the Mohand transform (MT)-based homotopy perturbation method (HPM) is applied to solve coupled linear and non-linear SWWEs in two dimensions. In Section 4, the variations in the water surface elevation have been demonstrated at different time levels and depths. Finally, conclusions have been drawn in Section 5.
2 Preliminaries
2.1 Shallow water wave equations (SWWEs)
The linear SWWEs in two dimensions [, ] may be given as follows:As such, the non-linear SWWEs [, ] may be written aswhere is the water surface elevation, is the time, and are the depth-averaged volume flux in and directions, is the acceleration due to gravity , and is the water depth.
2.2 Mohand transform (MT) and some of its properties
The Mohand transform of the function is given by []
Mohand transform of some useful fundamental functions is given in Table 1.
TABLE 1
| 1 | 1 | |
| 2 | 1 | |
| 3 | ||
| 4 | ||
| 5 |
MT of frequently encountered functions.
The linearity property of MT is given as follows []: if , then , where and are the arbitrary constants.
Let ; to obtain MT of a partial derivative, we use integration by parts, and then we derive the following []:
The function is called the inverse Mohand transform of the function if it has the following property: . The inverse Mohand transform of some useful fundamental functions [] is given in Table2.
TABLE 2
| S. No. | ||
|---|---|---|
| 1 | 1 | |
| 2 | 1 | |
| 3 | ||
| 4 | ||
| 5 |
Inverse MT of frequently encountered functions.
2.3 Homotopy perturbation method (HPM)
To illustrate the basic ideas of the HPM [], we consider the following non-linear differential equation:with the following boundary conditions:Here, is a general differential operator, is a boundary operator, is a known analytic function, and is the boundary of the domain .
The operator can be divided into two parts and , where is linear and is non-linear; Eq. 5 can be rewritten as follows:
Using the homotopy technique, one may construct a homotopy , which satisfies the following:orwhere is an embedding parameter and is an initial approximation of Eq. 5, which satisfies the boundary conditions. Evidently, from Eqs 8, 9, we obtain the following:The changing process of from zero to unity is same as that of from to . In topology, this is called deformation, and and are called homotopic functions.
We consider the solution to Eq. 8 as a power series in , given as follows:
The approximate solution to Eq. 5 will be obtained by substituting in Eq. 11.
3 Approximate solutions to shallow water wave equations using the Mohand transform-based homotopy perturbation method (MHPM)
3.1 Linear shallow water wave equations
In this section, we obtain the series solution to linear two-dimensional shallow water Eq. 1 by the MHPM.
Let us consider the Gaussian initial conditions for and Q, given as follows:
Applying Mohand transform in Eq. 1, we obtain
Using the derivative characteristics of Mohand transform, Eq. 14 may be written as
Taking the inverse of Mohand transform on both sides of Eq. 17, we may obtain the following:
Similarly, on solving Eqs 15, 16 by applying the same aforementioned procedure, we may obtain
Now using the HPM technique, we may construct a homotopy for Eqs 18, 19, 20 in the following form:
By simplifying the aforementioned formulas, we may obtain
Let us consider the solution to Eq. 1 in the power series form as:we may get the following results by applying Eq. 23 in Eq. 22 and comparing the coefficients of :
Comparing the coefficients of , we getSimilarly, comparing the coefficients of and , we get
In the same way, by correlating the coefficients of , we obtain the following generic pattern:Finally, by performing MT on the right hand side of Eq. 25 and then taking inverse MT, at a constant depth , we obtain
We can also get by Eq. 28; the solution of Eq. 1 may be found in the following form:
3.2 Non-linear shallow water wave equations
In this section, we obtain the series solution to the non-linear shallow water wave equation in two-dimensional Eq. 2 using the MHPM with the same initial conditions in Eq. 24.
Applying the MT and HPM in Eq. 2, we may obtain
Here, we again consider the power series solution to Eq. 2 as in (Eq. 23). We may get the following results by substituting Eq. 23 in Eq. 31 and comparing the coefficients of :Comparing the coefficients of , we getSimilarly, comparing the coefficients of , we get
By continuing to compare the coefficients of , we getBy this way, comparing the coefficients of , we getwhere ; calculating Eqs 33–36 with respect to , we get . The solution to Eq. 2 may be obtained in the series form by combining given as follows:
4 Numerical results and discussion
In this section, the results obtained for both linear and non-linear SWWEs have been presented, and the graphs of solutions obtained by the MHPM for various time levels (), water depths (), and distances () are investigated. Figure 1 represents the term-wise solutions to linear SWWEs in Eq. 1, which are obtained by the MHPM at the fixed basin depth of and time .
FIGURE 1
The term-wise solutions in and at the same time and depth are presented in Table 3. One may observe that the convergent solution is obtained as the number of terms increasing. It can also be seen from Figure 1 that the images possess a symmetrical property, and the images of non-linear SWWEs in Eq. 2 also have this property. Figures 2–7 show the plane images for the variation in the water surface elevation (WSE) (η) with respect to , when is fixed at . Depending on the symmetry of the images, the same results will be obtained when x has a fixed value.
TABLE 3
| Number of terms considered | Water surface elevation () |
|---|---|
| 3-term solution | 23.83 |
| 4-term solution | 23.82 |
| 5-term solution | 24.97 |
| 6-term solution | 24.97 |
| 7-term solution | 24.98 |
| 8-term solution | 24.98 |
MHPM at , and .
FIGURE 2
The water surface elevation (WSE) () of linear SWWEs are depicted at various time levels in Figure 2, at the fixed water depth of . Figure 2A shows the WSE vs. the global coordinate at three distinct time levels (, and ), whereas Figures 2B, C represent time levels and , respectively.
Diagrams in Figure 3 depict changes in the amplitude of linear SWWEs as the basin depth varies from to at two different time levels and . It can be seen in Table 4 that the WSE increases when the water depth () increases from .1 to .2 at a fixed time . Similar observations can be found for . Additionally, when time is increased from to , the WSE increases at a fixed water depth . Similarly, we get comparable results for .
FIGURE 3
TABLE 4
| x | ||||
|---|---|---|---|---|
| .1 | .2 | .1 | .2 | |
| 23.0 | 26.8973 | 186.4153 | 3.2798 × 103 | 2.4953 × 104 |
| 23.2 | 32.2187 | 242.9376 | 4.5436 × 103 | 3.7379 × 104 |
| 23.4 | 35.5048 | 276.6909 | 5.3663 × 103 | 4.5951 × 104 |
| 23.6 | 36.8368 | 288.2283 | 5.7489 × 103 | 5.0608 × 104 |
| 23.8 | 36.4372 | 279.7550 | 5.7251 × 103 | 5.1581 × 104 |
| 24.0 | 34.6307 | 254.7271 | 5.3546 × 103 | 4.9340 × 104 |
| 24.2 | 31.8003 | 217.3891 | 4.7140 × 103 | 4.4520 × 104 |
| 24.4 | 28.3449 | 172.3047 | 3.8886 × 103 | 3.7848 × 104 |
| 24.6 | 24.6413 | 123.9262 | 2.9639 × 103 | 3.0068 × 104 |
| 24.8 | 21.0149 | 76.2404 | 2.0189 × 103 | 2.1883 × 104 |
| 25.0 | 17.7191 | 32.5119 | 1.1204 × 103 | 1.3901 × 104 |
WSE of linear SWWEs at and with different values.
Figure 4 represents the results of WSE of non-linear SWWEs at various time levels with . Also in Figure 4A, the WSE is plotted against the spatial coordinate x at three different time levels , and , while Figures 4B, C include results at time levels and , respectively.
FIGURE 4
Figure 5 demonstrates changes in non-linear SWWEs as the water depth changes from to over two distinct time periods and , respectively. Table 5 shows that at a fixed time , the WSE increases when the water depth () increases from .1 to 0.5. Similar results can be found for . Additionally, at a fixed water depth , the WSE increases when time is increased from to . We obtain same results when .
FIGURE 5
TABLE 5
| x | ||||
|---|---|---|---|---|
| 0.1 | 0.5 | 0.1 | 0.5 | |
| 17.0 | 0.3977 × 104 | 0.2515 × 105 | 0.2576 × 106 | 1.6124 × 106 |
| 17.2 | 1.1673 × 104 | 0.3689 × 105 | 0.7504 × 106 | 2.3641 |
| 17.4 | 2.0858 × 104 | 0.4985 × 105 | 1.3385 × 106 | 3.1944 × 106 |
| 17.6 | 3.1173 × 104 | 0.6351 × 105 | 1.9989 × 106 | 4.0684 × 106 |
| 17.8 | 4.2052 × 104 | 0.7709 × 105 | 2.6953 × 106 | 4.9382 × 106 |
| 18.0 | 5.2745 × 104 | 0.8969 × 105 | 3.3799 × 106 | 5.7447 × 106 |
| 18.2 | 6.2365 × 104 | 1.0025 × 105 | 3.9958 × 106 | 6.4214 × 106 |
| 18.4 | 6.9961 × 104 | 1.0772 × 105 | 4.4821 × 106 | 6.8999 × 106 |
| 18.6 | 7.4621 × 104 | 1.1111 × 105 | 4.7804 × 106 | 7.1171 × 106 |
| 18.8 | 7.5578 × 104 | 1.0962 × 105 | 4.8417 × 106 | 7.0220 × 106 |
| 19.0 | 7.2314 × 104 | 1.0277 × 105 | 4.6328 × 106 | 6.5834 × 106 |
WSE of nonlinear SWWEs in and with different h.
Figure 6 shows a comparison graph between the MHPM and HPM for the water surface elevation of the linear SWWEs at the fixed water depth of at various time levels (). The comparison for the WSE of non-linear SWWEs at fixed and by the MHPM and HPM is given in Figure 7, respectively.
FIGURE 6
FIGURE 7
From Figures 2, 4, we may conclude that as time increases, the water surface elevations also increase both in linear and non-linear cases. It can be seen in Figures 3, 5 that the amplitudes of the linear and non-linear shallow water waves increase when the water depth () increases at a fixed time. The results obtained are consistent with the characteristics of shallow water waves. By comparing Figures 6, 7, the solutions to both linear and non-linear SWWEs obtained by the MHPM are found to be in excellent agreement with the results obtained by the HPM in [].
5 Conclusion
The main goal of this work is to solve shallow water wave equations using the Mohand transform-based homotopy perturbation method. The variation in the water surface height at different time levels and depths are given in this paper; the results are consistent with the characteristics of shallow water waves. It is proved that this method is a very good tool for solving SWWEs, which can easily be applied in finding out the approximate analytic solutions. The main advantage of this method over the HPM is that it is a powerful and efficient method to determine the analytical solution of the wave equation. In future research, we can also explore the effectiveness of this method in solving other problems, so as to improve the problem-solving efficiency.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding authors.
Author contributions
Conceptualization, YL, YZ, and JP; methodology, YL, YZ, and JP; software, YL and YZ; validation, YL and JP; writing—original draft preparation, YL, YZ, and JP; writing—review and editing, YL, YZ, and JP; supervision, JP; funding acquisition, JP. All authors have read and agreed to the published version of the manuscript.
Funding
This research was supported by the National Natural Science Foundation of China (Grant No. 10561151), the Basic Science Research Fund in the universities directly under the Inner Mongolia Autonomous Region (Grant No. JY20220003), and the Scientific Research Project of Hetao College (Grant No. HYZQ202122).
Acknowledgments
The authors would like to express their gratitude to the reviewers and editors for their insightful remarks and ideas, which helped them improve the paper’s quality.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
PerumandlaKChakravertyS. Comparison of solutions of linear and non-linear shallow water wave equations using homotopy perturbation method. Int J Numer Methods Heat Fluid flow (2017) 27(9):2015–29. 10.1108/HFF-09-2016-0329
2
LiuYShiYYuenD ASevreEODYuanXXingHL. Comparison of linear and nonlinear shallow wave water equations applied to tsunami waves over the China Sea. Acta Geotechnica (2009) 4(2):129–37. 10.1007/s11440-008-0073-0
3.
ChoYSSohnDHLeeSO. Practical modified scheme of linear shallow-water equations for distant propagation of tsunamis. Ocean Eng (2007) 34(11):1769–77. 10.1016/j.oceaneng.2006.08.014
4.
TandelPPatelHPatelT. Tsunami wave propagation model: A fractional approach. J Ocean Eng Sci (2021) 7:509–20. 10.1016/j.joes.2021.10.004
5.
WuNJChenCTsayTK. Application of weighted-least-square local polynomial approximation to 2D shallow water equation problems. Eng Anal Boundary Elem (2016) 68:124–34. 10.1016/j.enganabound.2016.04.010
6.
KarunakarPChakravertyS. 2-D shallow water wave equations with fuzzy parameters. Germany: Springer (2018).
7.
SahooMChakravertyS. Sawi transform based homotopy perturbation method for solving shallow water wave equations in fuzzy environment. Mathematics (2022) 10:2900. 10.3390/math10162900
8.
PerumandlaKChakravertyS. Solving shallow water equations with crisp and uncertain initial conditions. Int J Numer Methods Heat Fluid Flow (2018) 28:2801–15. 10.1108/hff-09-2017-0351
9.
NoeiaghdamLSidorovDNNoeiaghdamS. Dynamical control on the homotopy analysis method for solving nonlinear shallow water wave equation. Dynamic Syst Comp Sci Theor Appl (Dysc) (2020) 1847:012010. 10.1088/1742-6596/1847/1/012010
10.
SafariMSafariM. Analytical solution of two extended model equations for shallow water waves by He's variational iteration method. Am J Comput Math (2011) 1(4):235–9. 10.4236/ajcm.2011.14027
11.
SafariM. Analytical solution of two extended model equations for shallow water waves by adomian's decomposition method. Adv Pure Math (2011) 1(4):238–42. 10.4236/apm.2011.14042
12.
BekirAAksoyE. Exact solutions of extended shallow water wave equations by exp-function method. Int J Numer Methods Heat Fluid Flow (2013) 23(2):305–19. 10.1108/09615531311293489
13.
MohandMMahgoubA. The new integral transform Mohand transform. Appl Math Sci (2017) 12(2):113–20.
14.
AggarwaSSharmaNChauhanR. Solution of linear volterra integral equations of second kind using Mohand transform. Int J Res Advent Tech (2018) 6(11):3098–102.
15.
AggarwalSVyasASharmaSD. Mohand transform for handling convolution type volterra integro-differential equation of first kind. Int J Latest Tech Eng Manag Appl Sci (2020) IX(VII):78–84.
16.
NadeemMHeJHIslamA. The homotopy perturbation method for fractional differential equations: Part 1 Mohand transform. Int J Numer Methods Heat Fluid Flow (2021) 31(11):3490–504. 10.1108/hff-11-2020-0703
17.
KhandelwalRKhandelwalY. Solution of blasius equation concerning with Mohand transform. Int J Appl Comput Math (2020) 6(5):128. 10.1007/s40819-020-00871-w
18.
AggarwalSChaudharyR. A comparative study of Mohand and Laplace transforms. J Emerging Tech Innovative Res (2019) 6(2):230–40.
19.
AggarwalSSharmaNChauhanR. Duality relations of kamal transform with Laplace, laplace-carson, Aboodh, Sumudu, Elzaki, Mohand and Sawi transforms. SN Appl Sci (2020) 2(1):135–8. 10.1007/s42452-019-1896-z
20.
FangJHNadeemMHabibMAkgülA. Numerical investigation of nonlinear shock wave equations with fractional order in propagating disturbance. Symmetry (2022) 14:1179. 10.3390/sym14061179
21.
AlthobaitiSDubeyRSPrasadJG. Solution of local fractional generalized fokker-planck equation using local fractional Mohand adomian decomposition method. Fractals (2022) 30(01):2240028. 10.1142/s0218348x2240028x
22.
ShahRKhanHFarooqUBaleanuDKumamPArifM. A new analytical technique to solve system of fractional-order partial differential equations. IEEE Access (2019) 7:150037–50. 10.1109/access.2019.2946946
23.
HeJH. Homotopy perturbation technique. Comp Methods Appl Mech Eng (1999) 178(3-4):257–62. 10.1016/s0045-7825(99)00018-3
24.
AnjumNHeJH. Higher-order homotopy perturbation method for conservative nonlinear oscillators generally and microelectromechanical systems' oscillators particularly. Int J Mod Phys B (2020) 34(32):2050313. 10.1142/s0217979220503130
25.
YuDNHeJHGarcaAG. Homotopy perturbation method with an auxiliary parameter for nonlinear oscillators. J Low Frequency Noise, Vibration Active Control (2019) 38(3-4):1540–54. 10.1177/1461348418811028
26.
HeJHEl-DibYO. Homotopy perturbation method with three expansions. J Math Chem (2021) 59:1139–50. 10.1007/s10910-021-01237-3
27.
HeJH. Application of homotopy perturbation method to nonlinear wave equations. Chaos Solitons and Fractals (2005) 26(3):695–700. 10.1016/j.chaos.2005.03.006
28.
HemedaAA. Homotopy perturbation method for solving systems of nonlinear coupled equations. Appl Math Sci (2012) 6(96):4787–800.
29.
MadaniMFathizadehMKhanYYildirimA. On the coupling of the homotopy perturbation method and Laplace transformation. Math Comp Model (2011) 53:1937–45. 10.1016/j.mcm.2011.01.023
30.
MishraHK. A comparative study of variational iteration method and He-Laplace method. Appl Math (2012) 3:1193–201. 10.4236/am.2012.310174
31.
HeJHMoatimidGMZekryMH. Forced nonlinear oscillator in a fractal space. Facta Universitatis Series: Mech Eng (2022) 20(1):001–20. 10.22190/fume220118004h
32.
AnjumNHeJHAinQTTianD. LI-HE’S modified homotopy perturbation method for doubly-clamped electrically actuated microbeams-based microelectromechanical system. Ser Mech Eng (2021) 19(4):601–12. 10.22190/fume210112025a
33.
LiXXHeCH. Homotopy perturbation method coupled with the enhanced perturbation method. J Low Frequency Noise Vibration Active Control (2018) 0(0):1399–403. 10.1177/1461348418800554
34.
HeJHEl-DibYO. The enhanced homotoyp perturbation method for axxial vibration of strings. Facta Universitatis Ser Mech Eng (2021) 19(4):735–50. 10.22190/fume210125033h
35.
HeJEl-DibYO. The reducing rank method to solve third‐order Duffing equation with the homotopy perturbation. Numer Methods Partial Differential Equations (2020) 37:1800–8. 10.1002/num.22609
Summary
Keywords
approximate solutions, shallow water wave equations, Mohand transform, homotopy perturbation method, water surface elevation
Citation
Liu Y, Zhang Y and Pang J (2023) Approximate solutions to shallow water wave equations by the homotopy perturbation method coupled with Mohand transform. Front. Phys. 10:1118898. doi: 10.3389/fphy.2022.1118898
Received
08 December 2022
Accepted
27 December 2022
Published
13 January 2023
Volume
10 - 2022
Edited by
Chun-Hui He, Xi’an University of Architecture and Technology, China
Reviewed by
Dragan Marinkovic, Technical University of Berlin, Germany
Chaolu Temuer, Shanghai Maritime University, China
Dan Tian, Xi’an University of Architecture and Technology, China
Updates
Copyright
© 2023 Liu, Zhang and Pang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yanni Zhang, 641992269@qq.com; Jing Pang, pang_j@imut.edu.cn
This article was submitted to Interdisciplinary Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.