Double controlled $M$-metric spaces and some fixed point results

1.   Introduction

First, Peregrine ^[1] promulgated the regularized long-wave (RLW) equation to describe the propagation of unidirectional weakly nonlinear dispersive water waves. Furthermore, the authors of ^[2,3] engaged this equation to illuminate an enormous class of real-world problems as a substitute of the well-known Korteweg–De Vries (KdV) equation. These studies revealed that the RLW equation is more impressive than the latter one. The RLW equation takes part as a fundamental role in the study of the non-linear dispersive waves that have a lot of norms in various precise areas, e.g., magnetohydrodynamic waves as well as ion acoustic plasma waves, longitudinal dispersive and pressure waves in elastic rods and liquid-gas bubble mixtures, and rotating flow down a tube. Bona et al. ^[4] proposed an integer-ordered formulation of the RLW equation for describing the surface water wave's propagation in a channel.

The RLW equation has been studied by means of numerous procedures. For instance, this equation has been approximated, numerically, by the Galerkin finite element method (FEM) ^[15,16,34], Petrov-Galerkin FEM ^[17], least squares FEM ^[18], CBS and least squares CBS finite element methods (FEMs) ^[19,25], respectively, least squares quadratic B-spline (QdBS) FEM ^[20], splitting methods with CBS and QdBS FEMs ^[21,22], respectively, quintic B-spline Galerkin finite element method (QBS-GFEM) ^[23], linearized implicit finite difference method (FDM) ^[24], splitting-up technique with CBS and QdBS ^[26], quartic B-spline, QBS, and fourth-order CBS collocation techniques ^[27,28,29], respectively, CBS differential quadrature method ^[30], and lumped Galerkin QdBS FEM ^[33].

It is generally recognized that the trajectory's characteristic of the fractional derivatives is non-local as the remembrance outcome ^[5]. Many researchers prove that fractional differential equations (FDEs) are more appropriate than integer-order ones, as fractional derivatives demonstrate the memory and inherited possessions of several materials and processes ^[6,7,8,9].

Furthermore, the time-fractional partial differential equations (TFDEs) have generated further consideration for a number of real-life applications such as signal processing, electrical network systems, optics, financial estimation and forecast, mathematical biology, electromagnetic control theory, fluid flows in multi-dimension, material science, acoustics, biological systems associated with predator-prey models, etc. ^{[10,11,12,13]}. The application of fractional models is rising for enhanced precision in real-life models, and points out substantial necessities for improved fractional mathematical models. In ^[14], the author implemented Caputo's fractional derivative for dynamical investigation of a generalized tumor model. This derivative is being used for modeling of biological systems, comprising tumor growth. In biomedical research, tumor growth models have been expansively used to examine the dynamics of tumor expansion and estimating possible treatments.

Recently, the TFRLW equation was approximated by some analytical and numerical methods. For instant, the authors of ^[8] applied a method based on the q-homotopy analysis transform for approximating the TFRLW equation, while in ^[9], they presented a new fractional extension of the RLW equation. Besides, they used the fixed-point theorem to prove the existence and uniqueness of the solutions. Nikan et al. ^[35] obtained the traveling-wave solutions of the TFRLW equation using the radial basis function (RBF) collocation technique. Maarouf et al. ^[38] systematically examined the Lie group analysis technique of the TFRLW equation with the Riemann-Liouville fractional derivative. Naeem et al. ^[39] developed numerical methodologies that use the Yang transform, the homotopy perturbation method (HPM), and the Adomian decomposition method (ADM) to analyze this equation.

The TFRLW equation is one of the most substantial nonlinear evolution equations used to model various physical phenomena such as ion-acoustic plasma waves, shallow water waves, and longitudinal waves in elastic rods. Hossain et al. ^[40] used a modified simple equation integral technique in the TFRLW equation to create kink waves, anti-kink waves, brilliant and dark bell waves, double periodic waves, and combinations of solitons and periodic waves. The fractional RLW equations were used to mathematically model the nonlinear waves in the ocean, and similarly, the fractional RLW equations are used to describe the huge ocean waves known as tsunamis ^[41]. According to ^[42], the TFRLW equation can be used to study many phenomena such as plasma waves in complex media, water wave propagation in shallow water, and long-wave occupancy dynamics in the ocean, including tsunamis and tidal waves. Results in ^[43] can aid in understanding ion-acoustic waves in plasma, shallow water waves in oceans, and the development of a three-dimensional wave packet with finite depth on water under weak nonlinearity using the TFRLW. In ^[44], the authors obtained the soliton and periodic wave solutions for the TFRLW, which is the first step toward understanding ocean models' structural and physical behavior and coastal and harbor regions of the oceans. The Kudryashov approach was used to investigate the TFRLW problem in ^[45], which has prospective applications in applied science, nonlinear dynamics, mathematical physics, and engineering and is also important in biosciences, neurosciences, plasma physics, geochemistry, and fluid mechanics.

The most important part of this work is that the Caputo fractional derivative is used in the RLW equation for analyzing the nature of the displacement of shallow-water waves and ion acoustic plasma waves. It takes a broad view of the RLW equation for interpretation of the water waves. In interpretation of the excessive significance of fractional derivatives, we consider a TFRLW equation emerging in ion acoustic plasma waves (1). We use the cubic CBS collocation procedure to discretize the spatial derivatives. The Caputo's definition is used for time-fractional derivative.

with the initial and boundary conditions

where $p$ is a positive integer, $\hat \gamma$, $\hat \beta$ and dissipative term $\hat \mu$ are positive constants, $w(\zeta , \tau )$ represents the vertical displacement of the water surface, and the function $f(\zeta , \tau )$ denotes a source term. The notation $\frac{{{\partial ^\alpha }w}}{{\partial {\tau ^\alpha }}}$ indicates the Caputo's time-fractional derivative as

Definition 1.1. Caputo's integral and derivative of $g(\tau ) \in \mathbb{R}$ of order $\alpha \geqslant 0$ are, respectively, defined by

and

2.   Discretization of the TFRLW equation

This part implements the discretization process of the TFRLW equation by means of the CBS collocation procedure. First, we fix an identical partition of $\left[ {0, {\text{ }}T} \right]$ with size $\Delta \tau = \frac{T}{N}$, where $N$ is the partition's number of the time variable. Now, we discretize the fractional derivative $\frac{{{\partial ^\alpha }w}}{{\partial {\tau ^\alpha }}}$ for $0 < \alpha < 1$ at $\tau = {\tau _{j + 1}}$ by the $L1$ formula ^[20,30,31] as follows:

where ${a_0} = \frac{{\Delta {\tau ^{ - \alpha }}}}{{\left| \!{\overline {\, {2 - \alpha } \, }} \right. }}$, ${\chi _l} = {\left( {l + 1} \right)^{1 - \alpha }} - {l^{1 - \alpha }}$, $j = 0, 1, \ldots , N$, and ${\hat \delta ^{j + 1}}$ is the truncate error termed by ${\hat \delta ^{j + 1}} \leqslant {L_w}\Delta {\tau ^{2 - \alpha }}$, where the constant ${L_w}$ is associated to $w$.

Lemma 2.1. The element ${\chi _k}$, arising in Eq (5), fulfils the following possessions:

Next, the CBS collocation procedure is used to discretize derivatives of the spatial variable. The domain [a, b] is apportioned consistently with $h = \Delta \zeta = \frac{{b - a}}{M}$ by ${\zeta _i} = a + ih$, $i = 0, 1, \ldots , M$, such that $a = {\zeta _0} < {\zeta _1} < {\zeta _2} < \ldots < {\zeta _M} = b$. Now, we define CBS functions ${\Upsilon _i}(x)$ for $i = - 1, {\text{ }}0, {\text{ }}..., {\text{ }}M + 1$ as:

where $\left\{ {{\Upsilon _0}, {\Upsilon _1}, ..., {\Upsilon _M}, {\Upsilon _{M + 1}}} \right\}$ are preferred so that they form a basis over $\left[ {a, b} \right]$. The ${\Upsilon _i}(\zeta )$, ${\Upsilon _i}^\prime (\zeta )$, and ${\Upsilon _i}^{\prime \prime }(\zeta )$ at knot points are valued by the subsequent table (see Table 1).

where ${\sigma _1} = 1$, ${\sigma _2} = 4$, ${\sigma _3} = \frac{1}{h}$, and ${\sigma _4} = \frac{6}{{{h^2}}}$. We define the approximate solutions as

where ${C_i}\left( {{\tau _j}} \right)$ are unknown extents. The variation of $w\left( {\zeta , {\tau _j}} \right)$ is defined by

Using Eq (8), we approximate $w$ and its first-and second-order derivatives ${w_\zeta }$ and ${w_{\zeta \zeta }}$, respectively, with respect to $\zeta$ as

and

At $\tau = {\tau _{j + 1}}$, using the Eq (5) for $\frac{{{\partial ^\alpha }w}}{{\partial {\tau ^\alpha }}}$ and the $\theta$ scheme, we discretize problem (1) as

Now, to linearize the nonlinear term $\left( {{w^p}{w_\zeta }} \right)_i^{j + 1}$, we use the Rubin-Graves procedure as:

Taking $\theta = \frac{1}{2}$ and using Eq (13) in (12) with some manipulation, we have

where

Next, using the CBS collocation technique, we get

where

The Eq (15) forms a linear system with $M + 1$ equations and $M + 3$ unknowns. For making it uniquely solvable, we use the boundary conditions $w\left( {a, \tau } \right) = {\varphi _1}\left( \tau \right)$ and $w\left( {b, \tau } \right) = {\varphi _2}\left( \tau \right)$ as

From Eqs (16) and (17), we have

For $i = 0$ and $i = M$, inverting the Eq (18) in (15), we get

and

where

Equations (19), (15) and (20) form the following system of linear equations:

where $\tilde A_0^j = \frac{{{\sigma _2}{\sigma _3}}}{{{\sigma _1}}}B_0^j - \frac{{{\sigma _2}}}{{{\sigma _1}}}\bar A_0^j + \bar B_0^j$, $\tilde B_i^j = {\sigma _1}A_i^j + {\sigma _4}D$, $\tilde D_i^j = {\sigma _2}A_i^j - 2{\sigma _4}D$, and $\tilde A_M^j = - \frac{{{\sigma _2}{\sigma _3}}}{{{\sigma _1}}}B_M^j - \frac{{{\sigma _2}}}{{{\sigma _1}}}\bar A_M^j + \bar B_M^j$.

To solve the system (21), it is necessary to define the initial vector $\left( {C_0^0, C_1^0, ..., C_{M - 1}^0, C_M^0} \right)$ from $w(\zeta , 0) = \vartheta (\zeta )$ which provides $M + 1$ equations with $M + 3$ unknowns. To take out $C_{ - 1}^0$ and $C_{M + 1}^0$, we use ${w_\zeta }(a, 0) = {\vartheta _\zeta }(a)$ and ${w_\zeta }(b, 0) = {\vartheta _\zeta }(b)$ which gives

Now using Eq (22) and the initial condition, we have the subsequent system of linear equations:

3.   Stability analysis

This section establishes the stability for the discretized system of the TFRLW equation using the von Neumann scheme ^[32]. According to Duhamels' principle ^[36], the stability of an inhomogeneous system is the same as the stability of the corresponding homogeneous system. Therefore, we choose $f = 0$, and taking ${\left( {{w_\zeta }} \right)^p} = {\hat k_1}^p$ as locally constant to linearize ${w^p}{w_\zeta }$, and $\theta = \frac{1}{2}$, the Eq (12) can be written as

With the help of Eqs (9)–(11), we get

where ${A^*} = {a_0}{\sigma _1} - {\sigma _4}D, {\text{ }}{B^*} = {a_0}{\sigma _2} + 2{\sigma _4}D, {\text{ }}D = \frac{{\hat \mu }}{{\Delta \tau }}, {\text{ and }}{E^*} = \frac{1}{2}\left( {\hat \gamma + \hat \beta {{\hat k}_1}^p} \right)$.

Now, using the Fourier mode's growth factor $C_i^j = {\xi ^j}{e^{li\varepsilon h}}$, where $l = \sqrt { - 1}$, $\xi$ is the constraint depending on time, and we have

Now, we define $\xi _{\max }^j = \mathop {\max }\limits_{0 \leqslant i \leqslant j} \left| {{\xi ^i}} \right|$.

Using it in the Eq (26), and by means of the property $\sum\limits_{k = 0}^{j - 1} {\left( {\left( {{\chi _k} - {\chi _{k + 1}}} \right) + {\chi _j}} \right)} = 1$, we have

where ${S_1} = {\left( {2{a_0}{\sigma _1}\cos \varepsilon h + {a_0}{\sigma _2} + 2{\sigma _4}D - 2{\sigma _4}D\cos \varepsilon h} \right)^2} + 4{E^*}^2{\sigma _3}^2{\sin ^2}\varepsilon h$, and ${S_2} = \left( {2{A^*}\cos \varepsilon h + } \right.$ ${\left. {{B^*}} \right)^2} + 4{E^*}^2{\sin ^2}\varepsilon h$. Using the values ${A^*}, {B^*}, D, {\text{ }}{E^*}$, ${\sigma _1}$, ${\sigma _3}$, ${\sigma _4}$, and simplifying terms, we have

Hence, we conclude that $|\xi | \leqslant 1$. So, the discretized system of the TFRLW equation is unconditionally stable.

4.   Result and discussion

This division provides an example of the TFRLW equation to investigate the efficacy and validation of the projected technique. For this purpose, we use

and approximate error = $\frac{{|w({\zeta _j}, {\tau _{N + 1}}) - w({\zeta _i}, {\tau _N})|}}{{|w({\zeta _j}, {\tau _{N + 1}})|}}$,

where $W$ represents the exact solution. The ROC is analyzed by ${\text{ROC}} = \frac{{\ln \left( {err({h_1})/err({h_2})} \right)}}{{\ln \left( {{{{h_1}} / {{h_2}}}} \right)}}$, where the terms $err({h_1})$ and $err({h_2})$ represents errors with ${h_1}$ and ${h_2}$, in that order. The conservation possessions belonging to the TFRLW equation are measured by calculating quantities analogous to mass, momentum, and energy, respectively, as follows:

Now, we consider the TFRLW equation (1) with $\hat \gamma = 1 = \hat \beta = \hat \mu = p$ together with initial and boundary conditions $w(\zeta , 0) = 3\rho \sec {h^2}(\eta \zeta )$ and $w(a, \tau ) = w(b, \tau ) = 0$. Here, $3\rho$ is the amplitude and $\eta = \frac{1}{2}\sqrt {\frac{\rho }{{1 + \rho }}}$. When $\alpha$ = 1, the TFRLW equation has the subsequent single solitary wave solution $w(\zeta , \tau ) = 3\rho \sec {h^2}(\eta \zeta - \varpi \tau + {\zeta _0})$, where $\varpi = \frac{1}{2}\sqrt {\rho (1 + \rho )}$, ${\zeta _0}$ is an arbitrary constant, and $\eta$ and ${\varpi \mathord{\left/ {\vphantom {\varpi \eta }} \right. } \eta }$ represent the width and velocity, respectively. For all calculations, we have chosen ${\zeta _0}$ = 0.

Figure 1 signifies the estimated solution $w(\zeta , \tau )$ with admiration of the time $\tau$ for several values of $\rho$. From this figure, it can be revealed that the estimated solution $w(\zeta , \tau )$ increases as the value of $\rho$ increases. The approximate solutions with $h =$0.4, $\Delta \tau =$0.01, $\rho$ = 0.03 at times $\tau$ = 5, 7, 10 and 20 are demonstrated in Figure 2 for time-fractional orders $\alpha$ = 0.3, 0.5, 0.7, and 0.8. The figures show the influence of the Caputo order $\alpha$ of the fractional derivative on the evolution of the obtained solutions over time. An apparent dependence of $\alpha$ on the solutions can be seen clearly when the time is large. Table 2 shows the approximate errors together with an ROC for $\alpha =$0.9 with $\rho =$0.1, $h =$0.2, and $\tau = 0.1$ with respect to various time intervals. It can be perceived that the errors are very small and the projected method is linearly convergent with respect to the time variable. Table 3 shows the approximate errors for $\alpha =$0.4 with $\rho =$0.03, $h =$0.2, $\zeta =$2 and 4 for various time intervals at $\tau =$1 while Table 4 illustrates the approximate errors with $\rho =$0.1, $h =$0.2 for fractional orders $\alpha =$0.5, and 0.7 at times $\tau =$5 and 10. It can be noticed from these tables that the approximate errors are small which confirms the accuracy of the proposed technique.

Table 5 shows the ROC with respect to the space variable including errors in invariants for $\alpha =$1 with $\rho =$0.1, $\Delta \tau =$0.01 at $\tau =$1. It can be noticed from this table that the projected method is second-order convergent in space as well as that the small difference among the numerical and analytical values of ${I_1}$, ${I_2}$, and ${I_3}$ that extends in the invariants remains almost inconsistent for the duration of the computer run.

Table 6 demonstrates a comparison between the projected method and those available in refs. ^{[18,19,33,34]} in terms of ${L_2}$ and ${L_\infty }$ errors. The values of the single solitary wave's invariants are also compared for $\alpha =$1 with $\rho =$0.1, $h =$0.125, $\Delta \tau$ = 0.1, and $\zeta \in {\text{ }}[ - 40, 60]$ at various times. It is observed from Table 6 that the magnitudes in the invariants keep almost insistent in the course of the computer run. At $\tau$ = 16, the difference among the numerical and analytical values of the conservation constants are $\Delta {I_1}$ = 4.815941e-05, $\Delta {I_2}$ = 1.856193e-06, $\Delta {I_3}$ = 2.651635e-08. It is obvious from the table that the ${L_\infty }$ error norms at each time achieved by the projected method are much lower than those given in refs. ^{[18,19,33,34]}, However, the ${L_2}$ error norm is only higher than in ^[18] and is lower than the others. Also, the ${L_2}$ and ${L_\infty }$ errors in ^[33] are slightly smaller than those achieved by the projected method.

Table 7 compares the projected method and existing methods refs. ^{[19,23,33,34,37]} in terms of ${L_2}$ and ${L_\infty }$ errors as well as invariants for $\alpha =$1 with $\rho =$0.1, $h =$0.125, $\Delta \tau$ = 0.1, and $\zeta \in [ - 40, 60]$ at time $\tau$ = 20. It can be perceived from this table that the ${L_2}$ and ${L_\infty }$ error norms achieved by the projected method are very much smaller than those obtained in ^{[19,23,33,34,37]}, whereas, the errors obtained by the QBGM1 are almost similar to the projected method. The magnitudes in the invariants keep on nearly consistent in the course of the computer run. It is found that the difference among the numerical and analytical values of ${I_1}$, ${I_2}$, and ${I_3}$ are $\Delta {I_1}$ = 2.496862e-05, $\Delta {I_2}$ = 2.642822e-06, and $\Delta {I_3}$ = 3.886086e-08. Table 8 compares the invariants obtained by the projected method with ref. ^[37] and analytical quantities for $\alpha =$0.5 with $\rho =$0.03, $h =$0.1, $\Delta \tau$ = 0.0001 at various times $\tau$. It can be remarked from this table that the obtained invariant quantities are very close to analytical values and are much better than what is presented in ref. ^[37]. Table 9 shows the absolute errors in the invariants obtained by the projected method and analytical quantities for $\alpha =$0.6 with $\rho =$0.03, $h =$0.2, $\Delta \tau$ = 0.001 at various times $\tau$. It can be seen that the invariant quantities are nearly ${10^{ - 3}}$ accurate.

Figure 3 demonstrates the plots of the estimate solution $w(\zeta , \tau )$ contrasted with spatial as well as time variables $\zeta$ and $\tau$, respectively, for the values of $\alpha =$0.5 and $\alpha =$ 0.75 showing that the appearances of this figure are stable with ref. ^[35] (Figures 2 and 3). Figure 4 illustrates the approximate errors for $\alpha =$0.5 with $\rho =$0.1, $h =$0.1 for various time interval sizes $\Delta \tau$ at $\tau =$0.1. It can be seen from this figure that the approximate errors are decreasing on increasing $\Delta \tau$. Also, it is observed that the approximate errors are less than ${10^{ - 5}}$ which shows the accuracy of the projected method. The 3D plot of the approximate errors for $\alpha =$0.9 with $\rho =$0.1, $h =$0.2, $\Delta \tau$ = 0.001, and $\tau \in [0, 0.1]$ is depicted in Figure 5. The depiction of single solitary wave solutions with absolute errors by assuming $\alpha =$1, $h =$0.3, $\Delta \tau$ = 0.1 for $\rho =$0.1 and $\rho =$0.03 at $\tau = 1$ is described in Figure 6.

5.   Conclusions

The traveling-wave solutions are obtained for the TFRLW equation via a CBS collocation technique. The spatial derivatives are discretized by the aforesaid technique while the time-fractional derivative is discretized through Caputo's definition. The nonlinear term is commenced by the Rubin-Graves linearization procedure. The von-Neumann analysis confirms that the discretized structure of the TFRLW equation is enthusiastically stable. It is also established that the technique is second-order convergent in the spatial variable while linearly convergent in time. Three invariant capacities corresponding to mass, momentum, and energy are assessed for further justification. It is demonstrated that these invariants remain almost inconsistent for the duration of the computer run, and absolute errors are very small, approximately $\approx {10^{ - 8}}$ to ${10^{ - 5}}$. It is also observed that the obtained results by the projected technique are much better than the existing ones in ^{[18,19,23,33,34,37]}.

Author contributions

All authors of this article have been contributed equally. All authors have read and approved the final version of the manuscript for publication.

Use of Generative-AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

The authors gratefully acknowledge the funding of the Deanship of Graduate Studies and Scientific Research, Jazan University, Saudi Arabia, through project number RG24-S011.

Conflict of interest

There is no competing interest among the authors regarding the publication of the article.