CMES CMES CMES Computer Modeling in Engineering & Sciences 1526-1506 1526-1492 Tech Science Press USA 21332 10.32604/cmes.2022.021332 Article The Fractional Investigation of Fornberg-Whitham Equation Using an Efficient Technique The Fractional Investigation of Fornberg-Whitham Equation Using an Efficient Technique The Fractional Investigation of Fornberg-Whitham Equation Using an Efficient Technique KhanHassan1 2 KumamPoom3 4poom.kum@kmutt.ac.th NawazAsif1 KhanQasim1 KhanShahbaz1 Department of Mathematics, Abdul Wali Khan University, Mardan, 23200, Pakistan Department of Mathematics, Near East University TRNC, Mersin, 10, Turkey Center of Excellence in Theoretical and Computational Science (TaCS-CoE) & KMUTTFixed Point Research Laboratory, Room SCL 802 Fixed Point Laboratory, Science Laboratory Building, Departments of Mathematics, Faculty of Science, King Mongkut's University of Technology Thonburi (KMUTT), Bangkok, 10140, Thailand Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, 40402, Taiwan Corresponding Author: Poom Kumam. Email: poom.kum@kmutt.ac.th 27 09 2022 135 1 259 273 0912022 1752022 © 2023 Khan et al. 2023 Khan et al. This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In the last few decades, it has become increasingly clear that fractional calculus always plays a very significant role in various branches of applied sciences. For this reason, fractional partial differential equations (FPDEs) are of more importance to model the different physical processes in nature more accurately. Therefore, the analytical or numerical solutions to these problems are taken into serious consideration and several techniques or algorithms have been developed for their solution. In the current work, the idea of fractional calculus has been used, and fractional Fornberg Whitham equation (FFWE) is represented in its fractional view analysis. A well-known method which is residual power series method (RPSM), is then implemented to solve FFWE. The RPSM results are discussed through graphs and tables which conform to the higher accuracy of the proposed technique. The solutions at different fractional orders are obtained and shown to be convergent toward an integer-order solution. Because the RPSM procedure is simple and straightforward, it can be extended to solve other FPDEs and their systems.

Caputo derivative fractional partial differential equations fornberg-whitham residual power series method
Introduction

In mathematical physics, the Fornberg-Whitham equation is a fundamental mathematical model. The Fornberg-Whitham equation [1,2] is written as: uξuξ+u=uuuu+3uu,ξ>0.

This equation was introduced to investigate how non-linear dispersive water waves break. The Fornberg-Whitham equation is shown to allow peakon solutions, as well as the occurrence of wave breaking, as a mathematical model for waves of limiting height. Fractional calculus (FC) is now widely used and accepted, owing to its well-established applications in a wide range of seemingly disparate domains of science and engineering . Many scholars, including Gupta et al. , Merdan et al. , Singh et al. , have examined the fractional extension of the Fornberg-Whitham equation relevant to the Caputo fractional derivative.

FC has the potential to explain various difficult phenomena like memory and heredity. In recent years, researchers have taken a keen interest in the subject of fractional differential equations (FDEs), such as viscoelasticity, fluid mechanics, nanotechnology, electrochemistry, modelling for shape memory polymers, biological population models, optics and signal processing, modelling control theory, the damping behaviour of materials, economics and chemistry, signal processing, creeping and relaxation for viscoelastic materials and diffusion and reaction processes . Accurate modelling of time-fractional two-mode coupled burgers equation is done with the help of fractional derivatives (FD) .

Analytical and numerical techniques are frequently used for the solution of FPDEs and their systems. The fractional problems that have been modelled by using FPDEs are found in various disciplines, because the mathematical modelling of real-life phenomena is usually modelled accurately by using FPDEs.

In this connection, the important fractional mathematical models are solved by using various techniques such as Chun-hui He’s algorithm , reproducing kernel method , Runge-Kutta method for time discretization and Fourier transform for spatial discretization , Fourier spectral method , He-Laplace method , Fractional variational iteration method , operational matrix of fractional Riemann-Liouville integration with Legendre basis and zeros of Chebyshev polynomials . Khan et al.  have solved nonlinear fractional differential equations using an efficient approach. The non-linear differential equations are solved by using perturbation transform method and He’s polynomials . The one-dimensional non-homogeneous partial differential equations with a variable coefficient are solved by using homotopy perturbation method and Laplace transformation . Linear and non-linear differential equations arising in circuit analysis are investigated by using Maclaurin series method . The analysis of Caputo fractional-order dynamics of Middle East Lungs Coronavirus (MERS-CoV) model are discussed in . He’s fractional derivative and fractional complex transform are implemented for the time fractional Camassa-Holm equation . The fractional complex transform is utilised to solve time-fractional Schro¨dinger equation . A two-scale fractal theory is suggested for the population dynamics . The idea of the two-scale fractal is utilised to model the time-fractional tsunami wave traveling on an unsmooth surface . A study of projectile motion in a quadratic resistant medium via fractional differential transform method . The design of a variable fractional delay (VFD) FIR filter is handled using least square method . Homotopy analysis Sumudu transform method is applied to solve delay Fractional Bagley-Torvik equation (HASTM) . Coupled fractional Navier-Stokes equation is computed numerically using the proposed q-homotopy analysis transform method (q-HATM) . The solutions of the fifth-order dispersive equations with porous medium type non-linearity is investigated by the classical Riccati equations method . The analytical solution of time-fractional Navier-Stokes equation is calculated by using Adomian decomposition method  and similarly the Residual power series method (RPSM) was used to analyze the solution of fractional partial differential equations .

Many researchers have worked hard to find the solutions to FPDEs by using RPSM and other novel techniques that have been used for the solutions of FPDEs like: Senol et al.  solved the time-fractional nonlinear coupled Jaulent-Miodek system with energy-dependent Schrodinger Potential using RPSM in 2019, Korpinar et al.  have analysed the solution of the fractional cancer model by RPSM, Kurt  has implemented RPSM to obtain the solution of fractional Bogoyavlensky-Konopelchenko equation, Xu et al.  have implemented RPSM to obtain the solution of fractional Boussinesq equations, Freihet et al.  have found the solution of fractional Burgers-Huxley equations in by using RPSM. Jena et al.  have found the solution of the fractional model of the vibration equation of large membranes using RPSM.

In the current work, the solution of FFWE is investigated by RPSM. The RPSM was found to be a very effective technique for finding the analytical solution of FFWE . In the current work, we have applied the RPSM technique to solve different FPDEs and obtained the series form solutions. The closed form series solution is achieved for FFWE by using the proposed analytical technique. The accuracy of RPSM is represented by using graphs and tables. The solutions at different fractional orders are of greater interest, which show the useful information about the actual dynamics of the given problems. The graphs show the convergence of fractional solutions of the targeted problem towards integer order solution. The proposed method is also applied to other physical problems in applied sciences and engineering. The article layout is organized as follows; the fundamental concepts regarding FC are described in Section 2; the basic methodology is discussed in Section 3; the application of RPSM in the Fornberg-Whitham equation involving the Caputo fractional derivative is contained in Section 4; results and discussion are in Section 5 and the conclusion is given in Section 6.

Basic Definitions

In this section, we discussed some preliminaries and definitions.

Definition

The integral operator of Reimann-Liouville having order ϑ is given as [45,46]: (Iζϑh)(ζ)={1Γ(ϑ)0ζh(τ)(ζτ)1ϑdτ,ϑ>0,h(ζ),ϑ=0.

Its fractional derivative for ϑ0 is defined as: (Dζϑh)(ζ)=(ddζ)m(Imϑh)(ζ),(ϑ>0,m<ϑ<m1),where m is an integer.

Definition

The Caputo FD operator of the fractional order σ is given as : Danσf(ξ)=Ianσf(n)(ξ)=1Γ(nσ)a(τ)nσ1fn(τ)dτ.n1<σn,n.

The operators Dσ and Iσ satisfy the following properties:

DaσIaσf(ξ)=f(ξ).

Daσc=0,foranyconstantc.

Iaσξϑ=Γ(ϑ+1)Γ(ϑ+1σ)ξϑσ.

IaσDaσf(ξ)=f(ξ)j=0n1f(j)(0)ξjj!, ξ>0.

Definition

A Power Series (PS) expansion is defined as : n=0fn()(ξξo)nσ=fo()+f1()(ξξo)σ+f2(ξξo)()2σ+f3()(ξξo)3σ+,0n1<σn,ξξo.

Theorem

Assume that u(,ξ) has a multiple fractional PS representation at ξ=ξo of the form , μ(,ξ)=n=0fn()(ξξo)nσ.0n1<σn,I,ξoξ<ξo+.

If Dξnσ are continuous on I×(ξo,ξo+), n=0,1,, then the coefficients of fn() are given as: fn()=Dξnσu(,ξo)Γ(1+nσ),n=0,1,2,.

RPSM Methodology

To the understand the procedure of RPSM , let consider FPDEs of the form Dξσu(,ξ)=L(u)+N(u),ξ>0,0<σ1,having initial condition, u(,0)=f(),where Dξσ is Caputo type fractional derivative of the function u(,ξ), L(u) is linear term and N(u) is the non-linear term. Initially, RPSM is used by taking the fractional PS expansion about a fix point ξ = ξo, that is, f(,ξ)=n=0Dξnσfn(,ξ)Γ(nσ+1)(ξξo)nσ,0<σ1,[a,b],0ξ<,the 0th RPSM approximate solutions of u(,ξ) is given below as: uo(,0)=f().

Let uk(, ξ) represent the kth term, uk(,ξ)=n=0kfn()ξnσΓ(1+nσ),exchange the kth truncated series uk(,ξ) and of Eq. (3) in Eq. (2) which is the kth residual function and is given as: Resuk(,ξ)=Dξσu(,ξ)L(u)N(u),ξ>0,0<σ1,

Also obtain Resu(,ξ)=limkResuk(,ξ). It is observed that Resu(,ξ)=0, values of [a,b]. This Res u(,ξ) is infinitely differentiable at =a and dkdk1Resu(,0)=dkdk1Resuk(,0) which provide the basic rule of RPSM technique, allowing it to solve a broad range FPDEs. To determine the kth RPSM solution, we review Eq. (3) and then differentiate with respect to independent variables and ξ and finally replace ξ=0, we obtain f and constant parameters. By replacing these constant parameters in uk(,ξ), we obtain the kth truncated series, which we then assign to Eq. (4) to obtain an approximate solution. This method can be repeated using RPSM solutions for other arbitrary coefficient orders of Eq. (4).

Solution of Time-Fractional F-W Equation

Let us consider the time fractional F-W equation, uξσuξ+u=uuuu+3uu,ξ>0,0<σ1,with initial condition, u(,0)=e12.

Let uk(,ξ) denote the kth truncated series of u(,ξ), we have uk(,ξ)=n=0kfn()ξnσΓ(1+nσ).

We define the kth residual function of FW equation as: Resuk(,ξ)=Dξσuk(,ξ)(uk)ξ(,ξ)(uk)(,ξ)(uk)(,ξ)3(uk)(,ξ)(uk)(,ξ)+(uk)(,ξ)(uk)(,ξ)+(uk)(,ξ).

As Resu(,ξ)=0 and Resu(,ξ)=limkResuk(,ξ),I and ξ0, DξnσResu(,ξ)=0, the FD of a constant is 0 in the Caputo sense. Also, the FD DξnσResu(,ξ)=0 and Resuk(,ξ) and assign at ξ=0 for each n=0,1,2,3,,k. We first replace uk(,ξ) into Eq. (7), and find the FD formula Dξ(k1)αResu(,0)=0 for k=1,2,3

For the first step, k=1 put in FW equation Eq. (8), in the form as: Resu1(,ξ)=Dξσu1(,ξ)(u1)ξ(,ξ)(u1)(,ξ)(u1)(,ξ)3(u1)(,ξ)(u1)(,ξ)+(u1)(,ξ)(u1)(,ξ)+(u1)(,ξ)

Now we can write Eq. (7), for k=1, as: u1(,ξ)=f()+f1()ξσΓ(1+σ).

Differentiating of Eq. (10), and put Eq. (9), we have Resu1(,ξ)=f1()f1(xξσ1Γ(1+σ)[(f()+f1()ξσΓ(1+σ))(f()+f1()ξσΓ(1+σ))]3[(f()+f1()ξσΓ(1+σ))(f()+f1()ξσΓ(1+σ))]+[(f()+f1()ξσΓ(1+σ))(f()+f1()ξσΓ(1+σ))]+f()+f1()ξσΓ(1+σ),from Eq. (11), and ξ=0, we have Resu1(x,0)=f1()f()f()3f()f()+f()f()+f(),therefore f1()=f()f()+3f()f()f()f()f(),then from Eq. (7), and Eq. (12), we get u1(,ξ)=f()+ξσΓ(1+σ)[f()f()+3f()f()f()f()f(x],using Eq. (6), put in Eq. (13), then we have u1(,ξ)=e1223ξσΓ(1+σ)(e2).

Step 2. for k=2, put in Eq. (8), can be written as: Resu2(,ξ)=Dξσu2(,ξ)(u2)ξ(,ξ)(u2)(,ξ)(u2)(,ξ)3(u2)(,ξ)(u2)(,ξ)+(u2)(,ξ)(u2)(,ξ)+(u2)(,ξ).

Now we can written Eq. (7), for k=2, as: u2(,ξ)=f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ).

Similarly, Resu2(,ξ)=f1()+f2()ξσΓ(1+σ)[σf1()ξσ1Γ(1+σ)+2σf2()ξσ1Γ(1+2σ)][(f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ))(f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ))]3[(f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ))×((f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ))]+[(f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ))((f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ))]×f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ),applying Dξσ on both sides of Eq. (16) and then put ξ=0, yields DξσResu2(,0)=f2()[f()f1()+f1()f()]3[f()f1()+f1()f()]+f()f1()+f1()f()+f1(), f2()=f()f1()+f1()f()+3f()f1()+3f1()f()f()f1()f1()f()f1().

Now using Eq. (12) and Eq. (17), in Eq. (15), we get u2(,ξ)=f()+ξσΓ(1+σ)[f()f()+3f()f()f()f()f()]×ξ2σΓ(1+2σ)[f()f1()+f1()f()+3f()f1()+3f1()f()f()f1()f1()f()f1()].

Putting Eq. (6), and Eq. (12), in Eq. (18), we have u2(,ξ)=e1223ξσΓ(1+σ)(e2)+13ξ2σΓ(1+2σ)(e2).

For the third step, putting k=3, in FW Eq. (8), as: Resu3(,ξ)=Dξσu3(,ξ)(u3)ξ(,ξ)(u3)(,ξ)(u3)(,ξ)3(u3)(,ξ)(u3)(,ξ)+(u3)(,ξ)(u3)(,ξ)+(u3)(,ξ).

From Eq. (7), put k=3, the truncated series as: u3(,ξ)=f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ)+f3()ξ3σΓ(1+3σ),same as the above procedure Resu3(,ξ)=f1()+f2()ξσΓ(1+σ)+f3()ξ2σΓ(1+2σ)[σf1()ξσ1Γ(1+σ)+2σf2()ξ2σ1Γ(1+2σ)+3σf3()ξ3σ1Γ(1+3σ)][(f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ)+f3()ξ3σΓ(1+3σ))×(f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ)+f3()ξ3σΓ(1+3σ))]3[(f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ)+f3()ξ3σΓ(1+3σ))×(f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ)+f3()ξ3σΓ(1+3σ))]+[(f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ)+f3()ξ3σΓ(1+3σ))×(f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ)+f3()ξ3σΓ(1+3σ))]+f()+f1()ξσΓ(1+σ)+f2()ξ2σΓ(1+2σ)+f3()ξ3σΓ(1+3σ).

Applying Dξ2σ on both sides of Eq. (21) and then put ξ=0, yields Dξ2σResu3(,0)=f3()[f()f2()+2f1()f1()+f2()f()]3[f()f2()+2f1()f1()+f2()f()]+f()f2()+2f1()f1()f2()f1()+f2(),therefore f3()=f()f2()+2f1()f1()+f2()f()+3f()f2()+6f1()f1()+3f2()f()f()f2()2f1()f1()f2()f1()f2(), u3(,ξ)=f()+ξσΓ(1+σ)[f()f()+3f()f()f()f()f()]×ξ2σΓ(1+2σ)[f()f1()+f1()f()+3f()f1()+3f1()f()f()f1()f1()f()f1()]+ξ3σΓ(1+3σ)[f()f2()+2f1()f1()+f2()f()+3f()f2()+6f1()f1()+3f2()f()f()f2()2f1()f1()f2()f1()f2()].

Using Eq. (6), Eq. (12), and Eq. (18), put in Eq. (22), we have u3(,ξ)=(e2)23ξσΓ(1+σ)(e2)+13ξ2σΓ(1+2σ)(e2)16ξ3σΓ(1+3σ)(e2).

Equivalence of the 3rd-order approximate solution of u3(,ξ) with the exact solution is required to demonstrate the capability of the RPSM. Even when a low-order approximate solution, u3(,ξ), is used, the numerical results show that RPSM performs admirably for the FW equation. The accuracy can be improved by using a higher-order RPSM solution. Eq. (23) provides the 3rd approximate solution of the FW equation. The results are discussed through graphs and tables in the next section of the paper.

Results and Discussion

In Fig. 1, the exact and RPSM solutions at σ=1 are presented, respectively. The RPSM solution is calculated for k=3, that is up to the third iteration. From Fig. 1, it is very clear that the exact and RPSM solutions are very close to each other. The comparison also confirmed the greater accuracy of the present method. The domain of the given problem is 44, while 0ξ4. Fig. 2 shows the RPSM solution up to the 3rd iteration at σ=0.9 and σ=0.8. The solution graph at the fractional order σ=0.9 is slightly different as compared to the integer order solution, which shows another useful dynamics of the given physical problem. Similarly, another solution plot of the targeted problem is discussed at σ=0.8 is slightly different as compared to the fractional order σ=0.9. The plots at σ=0.9 and σ=0.8 are vary closed to each other. In Fig. 3, the solution graphs of various mathematical models with different fractional orders of the FW equation are displayed. A very consistent relationship is observed between the solutions graphs at various fractional orders σ. The overall graphical analysis of the problem has provided the opportunity to select that mathematical model of the fractional order which has a close relationship with the actual data of the problem. It is concluded that we have upgraded the existing model with the fractional order derivative to represent the dynamics of any physical problem in a rather sophisticated manner as compared to an integer model. In Fig. 3, the 2D solution plot is presented which shows the exact solution of the FW equation. Fig. 4 shows the approximate solution of the given problem at σ=1, which is nearly equal to the exact solution. The solutions graphs obtained by exact and approximate solutions are compared with each other. The 2D plots in Figs. 3 and 4 are very close to each other in the same fractional order. Another plot in Fig. 4 shows us the combined graphs of RPSM solutions at different fractional orders σ.

Comparison plots of exact and approximate solutions at the fractional order σ = 1Comparison plots of approximate solutions at the fractional order σ = 0.9 and σ = 0.8Combine 3D graphs for different fractional orders σ and exact solution of 2D plot2D plot and combine plot for various fractional order σ

Table 1 lists the result of RPSM u3(,ξ) solution at various fractional orders σ. The numerical results have shown that RPSM works excellently to solve FPDEs, even when using low-order approximate solution. The accuracy can be improve by using higher-order approximate solutions of FFWE. The Numerical results for various particular cases of σ are presented in Figs. 13.

Comparison table of exact and RPSM solutions of <inline-formula id="ieqn-89"><mml:math id="mml-ieqn-89"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi mathvariant="normal">ℑ</mml:mi><mml:mo>,</mml:mo><mml:mi>ξ</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>
ξ Exact σ=1 σ=0.9 σ=0.8 σ=0.7 σ=0.6 σ=0.5
−4 0.01 0.1792480649 0.1795470609 0.1789669892 0.1780335597 0.1765509395 0.1742355152 0.1707005599
−3 0.01 0.2955300975 0.2960230583 0.2950666817 0.2935277167 0.2910832893 0.2872658000 0.2814376440
−2 0.01 0.4872467579 0.4880595131 0.4864827145 0.4839453902 0.4799152106 0.4736212348 0.4640122301
−1 0.01 0.8033340940 0.8046741005 0.8020743991 0.7978910586 0.7912464160 0.7808694041 0.7650268336
0 0.01 1.3244740080 1.3266833050 1.3223971220 1.3154999600 1.3045447970 1.2874359960 1.2613160130
1 0.01 2.1836884690 2.1873309860 2.1802642650 2.1688927670 2.1508307550 2.1226231120 2.0795585410
2 0.01 3.6002936290 3.6062991210 3.5946480680 3.5758996370 3.5461204150 3.4996138740 3.4286123980
3 0.01 5.9358806850 5.9457820710 5.9265727310 5.8956617930 5.8465641560 5.7698878340 5.6528261920
4 0.01 9.7866127450 9.8029373700 9.7712665230 9.7203030030 9.6393546840 9.5129368000 9.3199347810

In Table 2, the solutions comparison of RPSM with He’s HPM is discussed. It is analyzed that the RPSM solution up to three terms is quite closed to the He’s HPM solutions which is calculated up to five terms. The tables and graphs have confirmed the higher degree of accuracy of RPSM. In Table 3, some nomenclatures are given which are frequently used in this paper.

Comparison table of Absolute error (AE) of RPSM (upto third term) and He’s HPM (upto fifth term) [<xref ref-type="bibr" rid="ref-55">55</xref>] solutions
ξ RPSM3 at σ=1 He’s HPM5 at σ=1
0.3 0.2 0.002395942 0.000190752
0.4 0.1 0.000653162 0.000195353
0.8 0.5 0.017409892 0.000332255
0.6 0.9 0.0459783993 0.0006452894
Nomenclature
 FPDEs Fractional partial differential equations FDEs Fractional differential equations DEs Differential equations RPSM Residual power series method FFWE Fractional fornberg whitham equations
Conclusion

In this work, we have implemented an efficient analytical technique, which is known as RPSM, to get an approximate series solution of FFWE with initial conditions. The suggested problems are first converted into their fractional form of the derivative and then incorporated the Caputo definition into the given problem to define FD. The general formulation for the proposed problem is discussed and then implemented for the solutions of FFWE. The proposed technique is applied to both fractional and integer orders of the suggested problem. It is observed that the procedure of the present technique is very effective and straight-forward. For verification and a better understanding of the obtained solutions, the graphical and tabular scenarios are presented. In Figs. 13, the exact and approximate solutions of the problem are presented, respectively. The solution graphs are presented at different fractional orders of the derivatives and show the various useful dynamics of the problem. It is investigated through numerical results and graphs that the fractional solutions are convergent to the integer-order solution. The analysis has shown the best contact between the RPSM and the exact solutions to the problem. For this purpose, the 2D graphs are discussed through Graphs 3 and 4. The 2D graphs clearly indicate that the proposed technique provides sufficient accuracy for problems with fractional orders. On the basis of the above analysis, it can be concluded that the proposed technique can be expanded easily to solve other problems in science and engineering. The values of the exact and approximate solutions of various fractional orders of the proposed method are given in Table 1, respectively. The exact and RPSM solution analyses given in Table 1 confirmed the higher accuracy.

Authors Contribution: Hassan Khan (Supervision), Poom Kumam (Funding, Draft Writing), Asif Nawaz (Methodology), Qasim Khan (Methodology, Investigation), Shahbaz Khan (Draft Writing).

Availability of Data and Material: Not applicable.

Funding Statement: The authors acknowledge the financial support provided by the Center of Excellence in Theoretical and Computational Science (TaCS-CoE), KMUTT. This research project is supported by Thailand Science Research and Innovation (TSRI) Basic Research Fund: Fiscal year 2022 under Project No. FRB650048/0164.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.