Geotechnical characteristics of Belfast's estuarine clayey silt (sleech)

Bryan A McCabe; Barry M Lehane; Bryan A McCabe; Barry M Lehane

doi:10.3934/geosci.2025006

AIMS Geosciences

2025, Volume 11, Issue 1: 91-116. doi: 10.3934/geosci.2025006

Previous Article Next Article

Research article Special Issues

Geotechnical characteristics of Belfast's estuarine clayey silt (sleech)

Bryan A McCabe ^{1
,
,},
Barry M Lehane ²

1.
School of Engineering, University of Galway, Galway, Ireland
2.
School of Engineering, University of Western Australia, Perth, Australia

Received: 04 December 2024 Revised: 28 January 2025 Accepted: 26 February 2025 Published: 10 March 2025

Belfast, the capital city of Northern Ireland and the second largest city on the island of Ireland, has experienced significant construction activity in recent years, yet relatively little has been published on the problematic soft estuarine deposits known locally as sleech which underlie the city and its hinterland. Results of a detailed characterization of the sleech at a site near Holywood, Co. Down, 8 km northeast of the city centre, are presented in this paper. This characterization was carried out in conjunction with a unique suite of full-scale foundation load tests at the site, commencing in the late 1990s. In situ tests clearly identified distinct sandy and silty horizons within the deposit, the bulk of which is a high-plasticity lightly-overconsolidated organic clayey silt, with clay content increasing with depth and clay fraction dominated by illite and chlorite. Both the moderate organic content and presence of diatom microfossils explain the relatively high Atterberg limits, the high compression index and the high friction angle of this material. Undrained shear strengths in triaxial compression fall at the upper end of the expected range based on a widely-used correlation with overconsolidation ratio, but are broadly compatible with in situ shear vane strengths corrected for plasticity index. The constant volume friction angle is remarkably insensitive to the specimen's stress history and the particle size distribution. Despite the high silt content, permeabilities and coefficients of consolidation are more typical of a clay than a silt. While the sleech behaviour is shown to be broadly similar to the estuarine deposits at the well-characterised geotechnical test bed at Bothkennar in Scotland, the paper illustrates that low OCR clays can have their own distinguishing characteristics.

Keywords:

Citation: Bryan A McCabe, Barry M Lehane. Geotechnical characteristics of Belfast's estuarine clayey silt (sleech)[J]. AIMS Geosciences, 2025, 11(1): 91-116. doi: 10.3934/geosci.2025006

Related Papers:

[1]	Musawa Yahya Almusawa, Hassan Almusawa . Numerical analysis of the fractional nonlinear waves of fifth-order KdV and Kawahara equations under Caputo operator. AIMS Mathematics, 2024, 9(11): 31898-31925. doi: 10.3934/math.20241533
[2]	Aslı Alkan, Halil Anaç . The novel numerical solutions for time-fractional Fornberg-Whitham equation by using fractional natural transform decomposition method. AIMS Mathematics, 2024, 9(9): 25333-25359. doi: 10.3934/math.20241237
[3]	Hayman Thabet, Subhash Kendre, James Peters . Travelling wave solutions for fractional Korteweg-de Vries equations via an approximate-analytical method. AIMS Mathematics, 2019, 4(4): 1203-1222. doi: 10.3934/math.2019.4.1203
[4]	Muath Awadalla, Abdul Hamid Ganie, Dowlath Fathima, Adnan Khan, Jihan Alahmadi . A mathematical fractional model of waves on Shallow water surfaces: The Korteweg-de Vries equation. AIMS Mathematics, 2024, 9(5): 10561-10579. doi: 10.3934/math.2024516
[5]	M. Mossa Al-Sawalha, Rasool Shah, Adnan Khan, Osama Y. Ababneh, Thongchai Botmart . Fractional view analysis of Kersten-Krasil'shchik coupled KdV-mKdV systems with non-singular kernel derivatives. AIMS Mathematics, 2022, 7(10): 18334-18359. doi: 10.3934/math.20221010
[6]	Ihsan Ullah, Aman Ullah, Shabir Ahmad, Hijaz Ahmad, Taher A. Nofal . A survey of KdV-CDG equations via nonsingular fractional operators. AIMS Mathematics, 2023, 8(8): 18964-18981. doi: 10.3934/math.2023966
[7]	Muhammad Imran Liaqat, Sina Etemad, Shahram Rezapour, Choonkil Park . A novel analytical Aboodh residual power series method for solving linear and nonlinear time-fractional partial differential equations with variable coefficients. AIMS Mathematics, 2022, 7(9): 16917-16948. doi: 10.3934/math.2022929
[8]	Khalid Khan, Amir Ali, Manuel De la Sen, Muhammad Irfan . Localized modes in time-fractional modified coupled Korteweg-de Vries equation with singular and non-singular kernels. AIMS Mathematics, 2022, 7(2): 1580-1602. doi: 10.3934/math.2022092
[9]	Sayed Saifullah, Amir Ali, Zareen A. Khan . Analysis of nonlinear time-fractional Klein-Gordon equation with power law kernel. AIMS Mathematics, 2022, 7(4): 5275-5290. doi: 10.3934/math.2022293
[10]	Asif Khan, Tayyaba Akram, Arshad Khan, Shabir Ahmad, Kamsing Nonlaopon . Investigation of time fractional nonlinear KdV-Burgers equation under fractional operators with nonsingular kernels. AIMS Mathematics, 2023, 8(1): 1251-1268. doi: 10.3934/math.2023063

Abstract

1. Introduction

Currently, fractional-order calculus (FOC) plays a dynamic role in the development and applications of modern science and technology ^[1]. The FOC is a field of mathematics that deals with differentiation and integration under an arbitrary order of the operation; the order is not limited to the integer one but can also be any real or even complex number ^[2,3,4]. For a long time, it was ignored due to the integral complexity and freedom of integer order calculus, as well as the reality that FOC does not have entirely conventional physical or geometrical understandings ^[5]. Currently, there is a strong interest in generalising integer-order NPDEs to model complex engineering problems such as bio-engineering, electronics, viscoelasticity, control theory, thermal sciences, electrical networking, chemical, and fluid dynamics ^[6,7].

Korteweg-de Vries (KdV) equation has been extensively investigated in physics to describe the interaction and evaluation of non-linear waves. It was originally derived by D. J. Korteweg and G. de Vries, in an investigation to re-generate the solitary waves that govern a small amplitude in one dimension with lengthy propagating gravity waves spreading in a shallow channel of water ^[8]. Apart from solitary waves, the KdV equation has also been studied for numerous physical phenomena, including hydro-magnetic waves, lattice dynamics, physics of plasma such as electron-acoustic (EA) and ion-acoustic (IA) waves ^[9]. It is also used in aerodynamics, fluid dynamics, transportation of mass, boundary layer behavior, and modeling the mechanical behavior of materials modeled as a continuous mass for shock wave formations. The concept of excitation of the electron-acoustic (EA) mode was initially proposed by Fried and Gould in the reference ^[10]. They have noted that the Landau's dampening effect on the EA potentials weakens with increasing wave number. In further experiments, weak dampening of the EA in plasma containing both high- and low-temperature electrons was found to be caused by ^[11]. Such plasma circumstances had previously been noted in a number of contexts, including ^[12]. Iwamoto ^[13] has studied the development of the high frequency Langmuir mode and the electron-acoustic wave evolution in nonrelativistic electron-positron plasma. It has been demonstrated that, in compared to the Langmuir wave, the low frequency EA excitation Landau damped with a substantially higher growth rate. The propagation of electron acoustic waves (EAWs) in plasmas and in space environments has drawn a lot of interest due to its value in comprehending a variety of collective processes in laboratory equipment ^[14,15]. Irfan et al. ^[16] studied fast (Langmuir) and slow (electron-acoustic) modes in dense electron-positron-ions (EPI) plasma. They observed that there is a chance that the fast electron-acoustic (FEA) mode could occur. The slow electron-acoustic (SEA) solitons are not allowed to propagate because they have a negative phase dispersion impact. The EPI plasma may allow for the evolution of electron-holes at a relatively low positron concentration, leading to the formation of compressive FEA solitons. Broadband electrostatic noise (BEN) emissions in the auroral and other regions of the magnetosphere, such as the plasma sheet boundary layer (PSBL) and the polar cusp, have been detected by satellite data. The Korteweg-de Vries (KdV) equation is obtained from studies of small-amplitude EAWs in unmagnetized plasma using first-order reductive perturbation theory ^[17]. More than 20 percent of the solitary waves' amplitude is understated by the first order soliton solution. Consequently, higher-order fixes may be used to help solve this issue. Shewy et al. ^[18,19] investigated the contributions of higher-order corrections to the properties of electrostatic EAWs as well as the impacts of non-thermal distribution of hot electrons. The higher-order correction contribution is given by the time-fractional modified Korteweg-de Vries (TF-mKdV) which modulates the amplitude of the solitary wave.

The time-fractional modified KdV (TF-mKdV) equation is well-known for its contribution in the creation of the Lax pair and an unlimited quantity of conservation laws for the KdV equation ^[20]. The finite-gap integration approach and Whitham modulation theory are used to investigate the comprehensive classification of solutions to the defocusing complex modified KdV equation with step-like initial condition ^[21]. Wang et al. ^[22] employed the nonlinear steepest descent method of Deift and Zhou, the long-time asymptotics of the focusing Kundu-Eckhaus equation with nonzero boundary conditions at infinity. The TF-mKdV equation is obtained using perturbation expansions and the theory that the soliton width is modest in comparison to the plasma inhomogeneity scale length. In this case, the soliton retains its whole identity, as well as its amplitude, breadth, and speed. Further, Breathers and the unique soliton behavior of the TF-mKdV equation are well-known. Here, we study the localized solutions for time-fractional Korteweg-de Vries (TF-KdV) and time-fractional modified Korteweg-de Vries (TF-mKdV) with help of an analytical method called Aboodh transform decomposition method (ADM) given in ^[23].

From the conventional Fourier integral, the Aboodh transform is derived. The Aboodh transform was first proposed by Khalid Aboodh ^[24] to show how to solve various ordinary differential equations in the time domain. This was done because of the Aboodh transform's mathematical features and simplicity. Differential equations are typically solved using Fourier, Laplace, Elzaki, and Sumudu transforms as the primary mathematical methods. Motivated from these study we are interesting to apply the Aboodh transform together with decomposition method to find an approximate solution of our proposed model. The proposed models have been extensively used in modeling complex physical systems such as thermal as well as current flows in electric circuits. It also plays an important role in modeling particles vibrating in their lattices. It should be noted that the physical phenomena in the KdV equation will be considered as non-conservative to describe fractional differential equations. Recently, the time-fractional KdVB equation together with nonlinearities has been extensively investigated using various methods ^[25]. We consider the integer order of the time derivative term of the mKdV model given in ^[26] be in the arbitrary order of the following form

$\begin{eqnarray} &&\frac{\partial^\alpha u}{\partial t^\alpha}+a\, u^m\frac{\partial u}{\partial x}+b\, \frac{\partial^{3}u}{\partial x^3} = 0, \, \, \, \, 0 < \alpha\leq 1, \end{eqnarray}$

(1.1)

with

$\begin{eqnarray} u(x, 0) = f(x), \end{eqnarray}$

(1.2)

the numbers $a, b \in \mathbb{R}$ with $a$ is the non-linearity coefficient and $b$ is the dispersive coefficient of Eq (1.1) and the time fractional-order ( $0 < \alpha\leq 1$ ) derivative is supposed to be in Caputo's form. Equation (1.1) is known as time-fractional modified Korteweg-de Vries (TF-mKdV) equation and is used vastly in plasma physics to study solitary and shock waves ^[27]. When $m = 1$ , the TF-mKdV reduce to simple time fractional Korteweg-de Vries (TF-KdV) equation.

Recently, various descriptions for fractional-order operators are also broadly studied, for example, Caputo-Fabrizio (C-F), Liouville-Caputo, and Riemann-Liouville (R-L) ^[28]. The R-L derivatives involve the complexity of the assumed function with the power-law kernel, while, the LC derivatives include the complexity of the confined derivative of a considered problem with power-law. Recently, the fractional KdV type equations together with several other nonlinear systems involving the time-fractional order derivatives along various methods for their analytical solutions have been investigated. To investigate analytical as well as numerical solutions to FNPDEs/TF-mKdV equations, various techniques have been studied ^[29,30,31]. Here, we apply the Aboodh decomposition method (ADM) ^[23] to find an approximate solution of Eq (1.1) in Caputo's sense. The rest of the article is prepared as follows. The preliminaries section 2 includes the main definitions, remarks, and some important results regarding the proposed method. In section 3, the Aboodh transform and Aboodh decomposition method are discussed. In section 4, uniqueness and convergence analysis of the method is discussed. In section 5, we consider two examples in Caputo's form of TF-KdV and TF-mKdV equations with the application of the proposed method. This section also includes numerical solutions and discussion. In section 5, the accuracy and effectiveness of the ADM are confirmed by using the error analysis. We include the summary in section 6, while in section 7 we have given the future work and extension of our present work.

2. Preliminaries

In this section, we provide some basic definitions, theorems, lemmas, and remarks, which are essential for our proposed method. We also define some basic rules and definitions related to the Aboodh transform and decomposition methods with some related properties. The Aboodh transform (see the reference ^[32]) is defined for the functions of exponential order in the following set

$\begin{eqnarray} \mathcal{A} = \left\{u(t):\exists\, \, M, k_1, k_2, \, \, \left|u(t)\right|\le M\, e^{-vt}, \, \ k_1\leq v\leq k_2\right\}, \end{eqnarray}$

(2.1)

the constant $M$ must be a finite value for a given function $u(t)$ in the set $\mathcal{A}$ , while $k_1, k_2 \in \mathbb{R}$ . The Aboodh transform $\mathcal{A}$ $u(t)$ is given by

$\begin{eqnarray} \mathcal{A}(u(t))(v) = K(v) = \frac{1}{v}\int_{0}^{\infty}\, u(t)e^{-vt}dt, \quad\quad \, \, \, t\geq 0, \end{eqnarray}$

(2.2)

where $v \in [k_1, k_2]$ is the variable utilised to factor the variable $t$ in the function $u(.)$ .

Definition 1. ^[33] The Aboodh transform for two variables function $x\geq 0$ , $t\geq 0$ is defined by

$\begin{eqnarray} \mathcal{A}(u(x, t))(u, v) = K(u, v) = \frac{1}{uv}\int_{0}^{\infty}\int_{0}^{\infty} e^{-ut}\, e^{-vt}u(x, t)dxdt, \end{eqnarray}$

(2.3)

while the inverse Aboodh transform is give by

$\begin{equation*} u(x, t) = \mathcal{A}^{-1}\{K(u, v)\} = \frac{1}{2\pi\, i}\int\limits_{\beta-i\infty}^{\beta+i\infty}u\, e^{ut}\bigg[\frac{1}{2\pi\, i} \int\limits_{\gamma-i\infty}^{\gamma+i\infty}ve^{vx}K(u, v)\, dv\bigg]du, \end{equation*}$

where, $K(u, v)$ , is a function such that $u$ and $v$ defined such that $Re$ $(u)\geq \beta$ and $Re$ $(v)\geq \gamma$ . It should be noted that $\beta, \gamma \in \mathbb{R}$ to be considered accordingly.

Definition 2. ^[34] The Caputo's derivative for $\alpha > 0$ for $u(t)$ is defined by function $u : (0, \infty)\rightarrow {\rm I\!R}$ is given by

$\begin{eqnarray} ^cD^\alpha u(t) = \frac{1}{\Gamma{(n-\alpha)}}\int_{0}^{t}(t-s)^{n-\alpha-1}u^n(s)ds, \end{eqnarray}$

(2.4)

where $\alpha \in (n-1, n]$ , with $n = [\alpha] + 1, [\alpha]$ represent the integer part of a real number $\alpha$ , providing that the right-hand side of the above integral is continuously well-defined on $(0, \infty)$ .

Definition 3. ^[35] Applying Aboodh transform on function $u(x, t)$ in the Caputo's sense is given by

$\begin{equation} \mathcal{A}\{^cD^{\alpha}u(x, t)\} = v^\alpha K(u, v)-\sum\limits_{k = 0}^{n-1}\frac{u^{(k)}(x, 0)}{v^{2-\alpha+k}}. \end{equation}$

(2.5)

Definition 4. ^[35] The Aboodh transform of some other functions and derivatives are given as

$\begin{eqnarray*} \mathcal{A}\{t^n\} = \frac{n!}{v^{n+2}}, \quad n\in N, \, \, \, \, \mathcal{A}\{t^\alpha\} = \frac{\Gamma(\alpha+1)}{v^{\alpha+2}}, \quad \alpha \ge -1. \end{eqnarray*}$

3. Aboodh transform decomposition method

Here, we introduce ADM applied to attain a series of solutions of NODEs and NPDEs. It is a very much effective technique to find an approximate solution of dynamical systems. Writing Eq (1.1) in Caputo's sense, we obtain

$\begin{equation} \begin{aligned} ^cD_t^\alpha\, u+a\, u^m\frac{\partial u}{\partial x}+b\, \frac{\partial ^3u}{\partial x^3} = 0, \end{aligned} \end{equation}$

(3.1)

$0 < \alpha\leq 1$ and

$\begin{eqnarray} u(x, 0) = f(x). \end{eqnarray}$

(3.2)

Using ADM and the definitions given in the previous section, we obtain

$\begin{equation} \mathcal{A}\{^cD_t^\alpha\, u\}+a\mathcal{A}\{u^m\frac{\partial u}{\partial x}\}+b\mathcal{A}\{\frac{\partial^3 u}{\partial x^3}\} = 0. \end{equation}$

(3.3)

Applying the technique given in section 2 for ADM on fractional order, we obtain

$\begin{equation} \mathcal{A}\{u(x, t)\} = \frac{1}{v^2}\mathcal{A}\{u(x, 0)\}-a\frac{1}{v^\alpha}\mathcal{A} \{u^m\frac{\partial u}{\partial x}\}-b\frac{1}{v^\alpha}\mathcal{A}\{\frac{\partial^3 u}{\partial x^3}\}, \end{equation}$

(3.4)

where Eq (3.2) becomes

$\begin{eqnarray} \mathcal{A}\{u(x, 0)\} = K(u). \end{eqnarray}$

(3.5)

By considering the series solution

$\begin{equation} u = \sum\limits_{n = 0}^\infty u_n. \end{equation}$

(3.6)

The non-linearity can be obtained as

$\begin{eqnarray*} u^mu_x(x, t) = \sum\limits_{i = 0}^\infty A_n. \end{eqnarray*}$

The value of $A_n$ for the functions $u_i, \, i = 1, 2, 3, \dots$ , is given by ^[36]

$\begin{equation} A_n = \frac{1}{n!}\frac{d^n}{d\lambda^n}\left[\sum\limits_{k = 0}^n \lambda^ku^m_k(x, t)\sum\limits_{k = 0}^n \lambda^ku_{kx}(x, t)\right]_{\lambda = 0}. \end{equation}$

(3.7)

Applying an inverse Aboodh transform to above equation, we obtain

$\begin{equation} \begin{split} u_0& = \mathcal{A}^{-1}\{\frac{1}{v^2}K(u, 0)\} = u(x, 0).\\ u_1& = -a\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\{A_0\}\right\} -b\mathcal{A}^{-1}\left\{\frac{1}{v^\alpha}\mathcal{A}\{u_{0xxx}\}\right\}, \\ u_2& = -a\mathcal{A}^{-1}\left\{\frac{1}{v^\alpha}\mathcal{A}\{A_1\}\right\}- b\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\{u_{1xxx}\}\right\}, \\ u_3& = -a\mathcal{A}^{-1}\left\{\frac{1}{v^\alpha}\mathcal{A}\{A_2\}\right\}- b\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\{u_{2xxx}\}\right\}, \\ &\vdots\\ u_{n+1}& = -a\mathcal{A}^{-1}\left\{\frac{1}{v^\alpha}\mathcal{A}\{A_n\}\right\}- b\mathcal{A}^{-1}\left\{\frac{1}{v^\alpha}\mathcal{A}\{u_{nxxx}\}\right\}. \end{split} \end{equation}$

(3.8)

The general solution can be written as

$\begin{equation} u(x, t) = u_0(x, t)-a\mathcal{A}^{-1}\left\{\frac{1}{v^\alpha}\mathcal{A}\{A_n\}\right\}- b\mathcal{A}^{-1}\left\{\frac{1}{v^\alpha}\mathcal{A} \{u_{nxxx}\}\right\}. \end{equation}$

(3.9)

4. Uniqueness and convergence of the method

The necessary condition that ensures the presence of a unique solution is studied in this section, and the convergence of the solutions is discussed, we follow these theorems applied to the series solutions in ^[37,38].

Theorem 1. (Uniqueness theorem) For $0 < \gamma < 1$ , the solution (3.9) of Eq (3.1) is a unique solution. where

$\begin{eqnarray} \gamma = \frac{(\ell_1+\ell_2+\ell_3)t^{\alpha+1}}{\Gamma(\alpha)}, \end{eqnarray}$

(4.1)

where $\ell_1$ and $\ell_2$ , are Lipschitz constants, $0 < \alpha\leq 1$ and $\Gamma$ is the well known Gamma function can be defined as

$\Gamma(\alpha) = \int_0^\infty e^{-t} t^{\alpha-1} dt, \quad \mathbb{R} > 0.$

Proof. Define a mapping $\mathscr{D}$ : $\mathscr{B}$ $\rightarrow$ $\mathscr{B}$ with $\mathscr{B}$ = $(C [J], \left\lVert{.}\right\rVert)$ is the Banach space on $J = [0, T]$ . We can write Eq (3.9) as

$\begin{equation*} u_{n+1} = u(x, 0)+\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\mathscr{H}u_n(x, t)+\mathscr{M}u_n(x, t)+\mathscr{N}u_n(x, t)\bigg)\right\}, \end{equation*}$

where $\mathscr{H}u = \frac{\partial u}{\partial x}$ , $\mathscr{M}u = \frac{\partial^3u }{\partial x^3}$ and $\mathscr{N}u = u\frac{\partial u}{\partial x}$ , and suppose that $\mathscr{H}u$ , and $\mathscr{M}u$ are also Lipschitzian with

$\begin{eqnarray*} \left|\mathscr{H}u-\mathscr{H}\bar{u}\right| < \ell_1\left|u-\bar{u}\right|, \, \, \, \, \left|\mathscr{M}u-\mathscr{M}\bar{u}\right| < \ell_2\left|u-\bar{u}\right|, \end{eqnarray*}$

where $u, \bar{u}$ are the function's distinct values. Now we proceed as follow

$\begin{array}{l} \left\lVert{\mathscr{D}u\, -\mathscr{D}\bar{u}\, }\right\rVert = \max\limits_{t\in J}\left|\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\mathscr{H}u_n\, +\mathscr{M}u_n\, +\mathscr{N}u_n\, \bigg)\right\}- \mathcal{A}^{-1}\\ \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\mathscr{H}\bar{u_n}\, +\mathscr{M}\bar{u_n}\, +\mathscr{N}\bar{u_n}\, \bigg)\right\}\right|\\ \leq\max\limits_{t\in J}\Bigg|\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\mathscr{H}u_n\, -\mathscr{H}\bar{u}\, \bigg)\right\} +\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\mathscr{M}u\, -\mathscr{M}\bar{u}\, \bigg)\right\} \\ +\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\mathscr{N}u\, -\mathscr{N}\bar{u}\, \bigg)\right\}\Bigg| \\ \leq\max\limits_{t\in J}\Bigg|\ell_1\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(u\, -\bar{u}\, \bigg)\right\} +\ell_2\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(u\, -\bar{u}\, \bigg)\right\} \\ +\ell_3\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\mathscr{N}u\, -\bar{u}\, \bigg)\right\}\Bigg| \\ \leq\max\limits_{t\in J}(\ell_1+\ell_2+\ell_3) \Bigg|\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(u\, -\bar{u}\, \bigg)\right\}\Bigg| \\ = \frac{(\ell_1+\ell_2+\ell_3)t^{(\alpha-1)}}{\Gamma(\alpha)} \left\lVert{u\, -\bar{u}\, }\right\rVert. \end{array}$

The mapping is contraction under the condition $0 < \gamma < 1$ . As a result of the Banach fixed point theorem for contraction, Eq (3.1) has a unique solution.

Next, we discuss convergence analysis of the problems.

Theorem 2. (Convergence theorem) The solution of Eq (3.1) in general form will be convergence.

Proof. The Banach space of all continuous functions on the interval $J$ with the norm $\left\lVert{u(x, t)}\right\rVert$ = $\max_{t\in J}\left\lVert{u(x, t)}\right\rVert$ is denoted as $(C[J], \left\lVert{.}\right\rVert)$ . Define the $\left\{\mathscr{S} n\right\}$ sequence of partial sums that is $\mathscr{S}_n$ = $\sum_{j = 0}^{n}$ . With $n\geq m$ , let $\mathscr{S}_n$ and $\mathscr{S}_m$ be arbitrary partial sums. In this Banach space, we will show that $\mathscr{S}_n$ is a Cauchy sequence. This can obtained by employing a new formulation of Adomian polynomials.

$P(\mathscr{S}_n) = \bar{A}_{n}+\sum\limits_{r1 = 0}^{n-1}\bar{A}_{{r_1}}, \quad Q(\mathscr{S}_n) = \bar{A}_{n}+\sum\limits_{r2 = 0}^{n-1}\bar{A}_{{r_2}}.$

Now

$\begin{array}{l} \left\lVert{\mathscr{S}_n-\mathscr{S}_m}\right\rVert = \left|\sum\limits_{j = 0}^nu_j\, -\sum\limits_{k = 0}^mu_k\, \right| = \max\limits_{t\in J}\left|\sum\limits_{j = m+1}^nu_j\, \right|\\ \leq\max\limits_{t\in J}\Bigg|\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\sum\limits_{j = m+1}^n\mathscr{H}(u_{j-1}\, )\bigg)\right\} +\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\sum\limits_{j = m+1}^n\mathscr{M}(u_{j-1}\, )\bigg)\right\} \\ +\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\sum\limits_{j = m+1}^n(A_{j-1}\, )\bigg)\right\}\Bigg| \\ = \max\limits_{t\in J}\Bigg|\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\sum\limits_{j = m}^{n-1}\mathscr{H}(u_{j}\, )\bigg)\right\} +\\ \mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\sum\limits_{j = m}^{n-1}\mathscr{M}(u_{j}\, )\bigg)\right\} +\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\sum\limits_{j = m}^{n-1}(A_{j}\, )\bigg)\right\}\Bigg|\\ \leq\max\limits_{t\in J}\Bigg|\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\sum\limits_{j = m}^{n-1}\mathscr{H}(\mathscr{S}_{n-1})-\mathscr{H}(\mathscr{S}_{m-1})\bigg)\right\} \\+\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\sum\limits_{j = m}^{n-1}\mathscr{M}(\mathscr{S}_{n-1})-\mathscr{M}(\mathscr{S}_{m-1})\bigg)\right\} \\ +\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg(\sum\limits_{j = m}^{n-1}\mathscr{N}(\mathscr{S}_{n-1})-\mathscr{N}(\mathscr{S}_{m-1})\bigg)\right\}\Bigg| \end{array}$

$\begin{equation*} \begin{split} &\leq \ell_1\max\limits_{t\in J}\Bigg|\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg((\mathscr{S}_{n-1})-(\mathscr{S}_{m-1})\bigg)\right\}\Bigg| +\ell_2\max\limits_{t\in J}\Bigg|\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg((\mathscr{S}_{n-1})-(\mathscr{S}_{m-1})\bigg)\right\}\Bigg| \\&+\ell_3\max\limits_{t\in J}\Bigg|\mathcal{A}^{-1} \left\{\frac{1}{v^\alpha}\mathcal{A}\bigg((\mathscr{S}_{n-1})-(\mathscr{S}_{m-1})\bigg)\right\}\Bigg| = \frac{(\ell_1+\ell_2+\ell_3)t^{\alpha-1}}{\Gamma(\alpha)}\left\lVert{\mathscr{S}_{n-1}-\mathscr{S}_{m-1}}\right\rVert. \end{split} \end{equation*}$

Choosing $n = m+1$ then

$\left\lVert{\mathscr{S}_{m+1}-\mathscr{S}_{m}}\right\rVert\leq\gamma\left\lVert{\mathscr{S}_m-\mathscr{S}_{m-1}}\right\rVert\leq\gamma^2\left\lVert{\mathscr{S}_{m-1}-\mathscr{S}_{m-2}}\right\rVert\leq\cdots\cdots \leq\gamma^m\left\lVert{\mathscr{S}_1-\mathscr{S}_{0}}\right\rVert,$

with $\gamma = \frac{(\ell_1+\ell_2+\ell_3)t^{\alpha-1}}{\Gamma(\alpha)},$ by using the following triangular inequality

$\begin{equation*} \begin{split}\left\lVert{\mathscr{S}_n-\mathscr{S}_{m}}\right\rVert&\leq\left\lVert{\mathscr{S}_m-\mathscr{S}_{m+1}}\right\rVert\leq\left\lVert{\mathscr{S}_{m+1}-\mathscr{S}_{m+2}}\right\rVert\leq\cdots\cdots \leq\gamma^m\left\lVert{\mathscr{S}_1-\mathscr{S}_{0}}\right\rVert\\& \leq\left(\gamma^m+\gamma^{m+1}+\cdots+\gamma^n\right)\left\lVert{\mathscr{S}_1-\mathscr{S}_0}\right\rVert\\& \leq\gamma^m\left(1+\gamma+\gamma^2+\cdots+\gamma^{n-m-1}\right)\left\lVert{\mathscr{S}_1-\mathscr{S}_0}\right\rVert\\& \leq\gamma^m\left(\frac{1-\gamma^{n-m}}{1-\gamma}\right)\left\lVert{u_1}\right\rVert. \end{split} \end{equation*}$

Now by definition $0 < \gamma < 1$ , we have $1-\gamma^{n-m} < 1$ , thus we have

$\begin{equation} \left\lVert{\mathscr{S}_n-\mathscr{S}_m}\right\rVert\leq\frac{\gamma^m}{1-\gamma}\max\limits_{t\in J}\left\lVert{u_1}\right\rVert, \end{equation}$

(4.2)

and also as $|u| < \infty$ ( $u\,$ is bounded), therefore, $\left\lVert{\mathscr{S}_n-\mathscr{S}_m}\right\rVert\rightarrow 0$ , hence $\mathscr{S}_n$ is a Cauchy sequence in the Banach space $\mathscr{B}$ , hence $\sum_{j = 0}^{n}u_j\,$ is convergent.

5. Applications of ADM

In this section, we consider two specific examples of time-fractional Korteweg-de Vries (TF-KdV) and time-fractional modified Korteweg-de Vries (TF-mKdV) equations in the form of Eq (1.1) with some initial conditions and apply ADM to find their approximate solutions.

Next we consider two examples one for time-fractional order KdV and the other for time-fractional mKdV and use the propose method (ADM) to obtain an approximate solution.

Example 1. For $m = 1$ in Eq (3.1) we can take the following TF-KdV in Caputo's form ^[39]

$\begin{eqnarray} ^cD_t^\alpha u+au\frac{\partial u}{\partial x}+b\frac{\partial ^3u}{\partial x^3}& = &0, \ \ \, 0 < \alpha\leq 1, \end{eqnarray}$

(5.1)

with initial condition of the form

$\begin{eqnarray} u_0 = u(x, 0) = A \operatorname{sech}^2(kx). \end{eqnarray}$

(5.2)

The exact solution for $\alpha = 1$ of Eq (5.1) ^[39]

$\begin{equation} u\, = A\, \operatorname{sech}^2(kx-4k^3bt), \end{equation}$

(5.3)

where $A = \frac{12k^2b}{a}$ is amplitude and $4 b k^3$ is the speed of the wave function. Applying ADM to Eq (5.1) and decomposing the non-linear term $uu_x$ by using the Adomian formula (3.7), we obtain the following Adomian polynomials

$\begin{array}{l} A_0 = u_0u_{0x}, \quad A_1 = u_0u_{1x}+u_1u_{0x}, \quad A_2 = u_0u_{2x}+u_1u_{1x}+u_2u_{0x},\\ \quad A_3 = u_0u_{3x}+u_1u_{2x}+u_3u_{0x}+u_2u_{1x}, \end{array}$

putting the above values in Eq (3.8) with assumption $h_0(x) = u_0(x, 0)$ with $h_1(x) = h_2(x) = h_3(x) = 0$ , we obtain the following single soliton solution for the TF-KdV Eq (5.1) in the form

$\begin{eqnarray*} u_0 = h_0(x), \quad u_1 = \frac{t^\alpha}{\Gamma(\alpha+1)}f_1(x), \quad u_2 = -\frac{t^{2\alpha}}{\Gamma(2\alpha+1)}f_2(x), \quad u_3 = \frac{t^{3\alpha}}{\Gamma(3\alpha+1)}f_3(x), \end{eqnarray*}$

where

$\begin{eqnarray*} &f_1(x)& = ah_0h_{0x}+bh_{0xxx}, \quad f_2(x) = a(h_0f_{1x}+2h_0f_1h_{0x})+bf_{1xxx}, \\ &f_3(x)& = a\left[h_0f_{2x}+2h_0f_2h_{0x}+\frac{\Gamma(2\alpha+1)}{\Gamma(\alpha+1)}(2h_0f_1f_{1x}h_{0x}+f_{1}h_{0x})\right]+bf_{2xxx}. \end{eqnarray*}$

The final solution in the series form up to $\mathcal{O}(4)$ is given by

$\begin{equation} U(x, t) = h_0(x)+\frac{t^\alpha}{\Gamma(\alpha+1)}f_1(x)-\frac{t^{2\alpha}}{\Gamma(2\alpha+1)}f_2(x)+ \frac{t^{3\alpha}}{\Gamma(3\alpha+1)}f_3(x). \end{equation}$

(5.4)

For $\alpha = 1$ , we obtain

$\begin{equation*} u(x, t)\, = \frac{12k^2b}{a}\, \operatorname{sech}^2(kx-4k^3bt). \end{equation*}$

The ADM solution for different $\alpha$ are reported in , by considering $k = 0.1, a = 2, b = 0.2, t = 1$ .

Table 1. Evaluation between approximate (

$U(x, t))$ versus an exact solution (u(x, t) for different values a = 0.2, b = 0.4, k = 0.1,

$-10\leq x\leq 10$ and

$0\leq t\leq 1$ Exp.1.

t	$U(\alpha=0.3)$	$U(\alpha=0.5)$	$U(\alpha=0.7)$	$U(\alpha=0.9)$	$U(\alpha=1)$	Exact	MAE
x=-10
0.00	0.0050869	0.0050869	0.0050869	0.0050869	0.0050869	0.0050869
0.25	0.0050056	0.0050276	0.0050447	0.0050575	0.0050625	0.0050634
0.50	0.0049834	0.0049994	0.0050155	0.0050302	0.0050367	0.0050400	4.0619E-05
0.75	0.0049673	0.0049763	0.0049888	0.0050025	0.0050093	0.0050167
1.00	0.0049542	0.0049559	0.0049633	0.0049742	0.0049803	0.0049935
x=-5
0.00	0.030238	0.030238	0.030238	0.030238	0.030238	0.030238
0.25	0.029905	0.029984	0.030052	0.030105	0.030127	0.030128
0.50	0.029825	0.029877	0.029933	0.029989	0.030015	0.030018	1.2525E-04
0.75	0.02977	0.029793	0.029831	0.0050025	0.029902	0.029908
1.00	0.029726	0.029722	0.029737	0.0049742	0.029787	0.029798
x=0
0.00	0.0720000	0.0720000	0.0720000	0.0720000	0.0720000	0.0720000
0.25	0.071994	0.071997	0.071998	0.071999	0.0720000	0.0720000
0.50	0.07199	0.071993	0.071996	0.071998	0.0719987	0.0719987	6.2213E-06
0.75	0.071988	0.07199	0.071993	0.071995	0.071996	0.071996
1.00	0.071985	0.071987	0.071989	0.071992	0.071993	0.071993
x=5
0.00	0.030238	0.030238	0.030238	0.030238	0.030238	0.030238
0.25	0.030555	0.030483	0.03042	0.030369	0.030348	0.030349
0.50	0.030625	0.030582	00.030532	0.030481	0.030457	0.03046	1.0040E-04
0.75	0.030673	0.030657	0.030626	0.030587	0.030565	0.030571
1.00	0.030711	0.030719	0.03071	0.030687	0.030671	0.030682
x=10
0.00	0.0050869	0.0050869	0.0050869	0.0050869	0.0050869	0.0050869
0.25	0.005144	0.0051338	0.0050447	0.0051233	0.0051096	0.0051105
0.50	0.0051538	0.0051496	0.0050155	0.0051431	0.0051308	0.0051341	1.4575E-05
0.75	0.0051598	0.0051603	0.0049888	0.0051583	0.0051505	0.0051579
1.00	0.0051641	0.0051684	0.0049633	0.0051706	0.0051687	0.0051819

| Show Table

DownLoad: CSV

, displays the zeroth and higher order solitary wave solutions (5.3) and (5.4) versus the space variable $x$ respectively for time-fractional KdV Eq (5.1) for small time of $t$ . It infers that at $t > 0$ the higher order solitary wave solution is significantly modified. The wave three-dimensional profiles for solutions (5.3) and (5.4) in , also reveal that variations in temporal variables affect the amplitude and spatial extensions of the higher order wave solutions. We have plotted the approximate solution (5.3) against $t$ at $x = 0.5$ and $1$ respectively, with variation in the time-fractional index $\alpha$ [see , ]. A degree of enhancement in $\alpha$ reduces the spatial width of the solitary waves. Thus, the time-fractional index leads to a reduction in the wave dispersion that in turn localises the wave profile. The effect of the non-linearity coefficient "a" is shown in . We see that the amplitude of the solitary wave decreases regularly when the values of the non-linearity coefficient are increasing. This definitely occurs for wave equations because the non-linearity coefficient has an inverse proportionality with the solitary wave amplitude. On the other hand, the wave amplitude is gradually increasing with the increasing values of the dispersion coefficient "b" in . The dispersive coefficient is proportional to the amplitude of a solitary wave. We have depicted the higher order solution (5.4) for the TF-KdV Eq (5.1), verses $x$ at time $t = 3$ (solid curve), at 3.5 (dashed curve), and at 4 (dotted curve) [see ]. Obviously, the enhancement in $t$ gives rise to the pulse amplitude while increasing the spatial width. Moreover, the pulse shape solutions suffer from oscillations for $t > 0$ . The wave solution (5.4) presented in at different values of the temporal index ( $\alpha$ ) reveals strengthening of the wave with wider spatial extension at large $\alpha$ . One can see that the solitary wave solutions for the TF-KdV equation are observed. By employing an Aboodh decomposition (ADM) technique and with the addition of $t$ the Caputo operator for the solitary excitations are derived in Eq (5.4). Importantly, the wave solution significantly changes in time. The large amplitude electrostatic excitations in the auroral zone ^[27], associated with the intensified electric field. These may be described by the solutions (5.3) and (5.4) for the TF-KdV equation of Exp.1. Our results obtained in this paper are so important because the waves produced by the solutions of the TF-KdV may have very high amplitudes. The amplitudes of these waves can be reduced by using the fractional order of the wave equations. In other words, by taking the fractional order of the time derivative term, a better analysis of the waves can be achieved.

Figure 1. In the top panel, (a) shows exact solution of Eq (5.3), while (b) shows comparison of approximate (ADM) and exact solution Eq (5.4)

$k = 0.1, a = 0.2, b = 0.4$ , for small time

$t = 0.2$ . The bottom panel (c, d) are corresponding surface 3D plots for both solutions against the temporal variables

$(x, t)$ with parameters (

$\alpha = 1$ for exact, and

$0.7$ , for the approximate solution, with fixed

$k = 0.1, a = 2, b = 0.2, t = 1$ ).

DownLoad: Full-Size Img PowerPoint

Figure 2. The above plots show the dispersive behaviours of the solitary waves by varying time and the fractional order

$\alpha$ with fixed special variable (a)

$x = 0.5$ and (b)

$x = 1$ of the exact (

$\alpha = 1$ star black) solution Eq (5.3) and approximate (

$\alpha < 1$ blue) solution Eq (5.4) respectively.

DownLoad: Full-Size Img PowerPoint

Figure 3. Solitary waves profiles for different values of (a) the non-linearity coefficient "a" with b = 0.4 and (b) the dispersive coefficient "b" with a = 0.2, when the wave number k = 0.1, the fractional index

$\alpha = 0.5$ and t = 0.5 of Eq (5.4).

DownLoad: Full-Size Img PowerPoint

Figure 4. Plots of the approximate solution for (a) different values of t with

$\alpha = 0.9$ (b) different values of

$\alpha$ for larger time and the other parameters are

$a = 0.9, k = 0.5, b = 0.9$ of Eq (5.4).

DownLoad: Full-Size Img PowerPoint

Example 2. For $m = 2, a = 6, b = 1$ , the modified time-fractional KdV ^[26] in Caputo's sense is given by

$\begin{eqnarray} &^CD_t^\alpha\, v+6\, v^2\frac{\partial\, v}{\partial x}+\frac{\partial ^3\, v}{\partial x^3} = 0, \, \, \, 0 < \alpha\leq 1, \end{eqnarray}$

(5.5)

with

$\begin{eqnarray} v_0 = v(x, 0) = \sqrt{c} \operatorname{sech}^2(k+\sqrt{c}x). \end{eqnarray}$

(5.6)

When $\alpha = 1$ , Eq (5.5) gives ^[26]

$\begin{equation} v(x, t) = \sqrt{c}\, \operatorname{sech}(k+\sqrt{c}x-ct). \end{equation}$

(5.7)

Applying the ADM to Eq (5.5) and decompose the non-linear term $v^2v_x$ in Eq (5.5) by using the Adomian formula (3.7), we obtain

$\begin{array}{l} A_0 = v_0^2v_{0x}, \quad A_1 = v_0^2v_{1x}+v_1^2v_{0x}, \quad A_2 = v_0^2v_{2x}+v_1^2v_{1x}+v_2^2v_{0x}, \\ A_3 = v_0^2v_{3x}+v_1^2v_{2x}+v_3^2v_{0x}+v_2^2v_{1x}. \end{array}$

Putting $a = 6$ and $b = 1$ in Eq (3.8) and assume that $g_0(x) = v(x, 0)$ with $g_1(x) = g_2(x) = g_3(x)$ , we obtain

$\begin{eqnarray*} v_0 = g_0(x), \quad v_1 = \frac{t^\alpha}{\Gamma(\alpha+1)}q_1(x), \quad v_2 = -\frac{t^{2\alpha}}{\Gamma(2\alpha+1)}q_2(x), \quad v_3 = \frac{t^{3\alpha}}{\Gamma(3\alpha+1)}q_3(x), \end{eqnarray*}$

where

$\begin{eqnarray*} &q_1(x)& = 6g_0^2g_{0x}+g_{0xxx}, \quad q_2(x) = 6(g_0^2q_{1x}+2g_0q_1g_{0x})+q_{1xxx}, \\ &q_3(x)& = 6\left[g_0^2q_{2x}+2g_0q_2g_{0x}+\frac{\Gamma(2\alpha+1)}{\Gamma(\alpha+1)}(2g_0q_1q_{1x}g_{0x}+q_{1}^2g_{0x})\right]+q_{2xxx}. \end{eqnarray*}$

The solution up to $\mathcal{O}(4)$ can be written as

$\begin{equation} V(x, t) = g_0(x)+\frac{t^\alpha}{\Gamma(\alpha+1)}q_1(x)-\frac{t^{2\alpha}}{\Gamma(2\alpha+1)}q_2(x)+ \frac{t^{3\alpha}}{\Gamma(3\alpha+1)}q_3(x). \end{equation}$

(5.8)

For $\alpha = 1$ we found the following exact solution

$\begin{equation*} v(x, t) = \sqrt{c}\, \operatorname{sech}(k+\sqrt{c}x-ct). \end{equation*}$

We have displayed the zeroth (higher) order solitary wave solutions (5.7) and (5.8) in , respectively, for the modified time-fractional KdV Eq (5.5) verses x for small t. Notice that the wave solution of the TF-mKdV equation suffers only a relatively modest steeping effect with a reduction in the wave dispersion. Moreover, at $t > 0$ the higher-order solution illustrates noticeable modifications in the wave profiles. The three dimensional profiles Figure 5(c, d) for the time fractional (exact and approximate) solutions also confirm that both the solutions are quite similar. , display the solution (5.8) at $x = 0.5$ and 1 respectively with variation in $\alpha$ . Contrary to the TF-KdV solution, the wave solution for TF-mKdV admits reduced wave dispersion at $t > 0$ . We observed that, the dispersive approximate solution ( $\alpha < 1$ the blue curves) approaching to the exact solution $(\alpha = 1$ , black curve) when then $\alpha$ values are approaching to 1. In Figure 7(a), we depicted the approximate (ADM) solution against the special variable x by changing the wave speed "c = 0.1 (solid curve), 0.2 (dashed curve), and 0.3 (dotted curve)". We observed that when the speed "c" is increasing, the wave amplitude is increasing while its width is decreasing, which shows that the amplitude is directly proportional to the wave speed while the width of the wave is inversely related to the wave speed. Similarly, in , we displayed the approximate solution against "x" by changing the wave number "k = 0.5 (solid curve), 1 (dashed curve), and 1.5 (dotted curve)". We see that, for different values of the wave number and equal amplitudes, wave profiles are obtained. The wave solution (5.5) for TF-mKdV in shows at different times that the solitary waves suffer from oscillation due to the external perturbations. It also shows that the solitary solution for the TF-mKdV is relatively less stable against the perturbations. Moreover, the increase in $t$ and $\alpha$ modify the amplitude and width of the wave. Variation of the $\alpha$ value for the approximate solution is shown in against $x$ . We observed from the , when the temporal variable (t) and the fractional order $\alpha$ increasing the amplitudes of the wave profile getting larger values.

Figure 5. In the top panel, (a) shows exact solution for Eq (5.7) while (b) shows a comparison of exact and approximate solution (ADM) of Eq (5.8). The bottom panel (c, d) are the corresponding surface 3D plots against the temporal variables

$(x, t)$ with parameters (

$\alpha = 1$ for exact, and

$0.7$ , for the approximate solution, with fixed

$k = 0, c = 1, t = 0.2$ ).

DownLoad: Full-Size Img PowerPoint

Figure 6. Dispersive behaviours of the solitary waves by varying time and the fractional order

$\alpha$ with fixed special variable (a)

$x = 0.5$ and (b)

$x = 1$ for exact (

$\alpha = 1$ stars blue curve of Eq (5.7)) and approximate (

$\alpha < 1$ blue curves Eq (5.8)) respectively.

DownLoad: Full-Size Img PowerPoint

Figure 7. Solitary waves profiles for different values of (a) the speed "c" with k = 0.5 and (b) the wave number "k", with c = 0.2 when the fractional index

$\alpha = 0.7$ and t = 0.5 of Eq (5.8).

DownLoad: Full-Size Img PowerPoint

Figure 8. Plots of the approximate solution for (a) different values of

$t$ with

$\alpha = 0.9$ (b) different values of

$\alpha$ for larger time values and other parameters are

$c = 0.6, k = 0$ of Eq (5.8).

DownLoad: Full-Size Img PowerPoint

The mean absolute error (MAE) is a statistic that calculates the variance in errors amongst matching interpretations representing identical singularities. The relationships between an exact and an approximate solution, future time against beginning time, and one measuring method versus another quantity technique are all examples of $y$ versus $x$ . The MAE can be computed as

$\begin{eqnarray} MAE = \frac{1}{n}\sum\limits_{i = 1}^{n}|y_i-x_i|, \end{eqnarray}$

(5.9)

where, $x_i$ is the expected value and $y_i$ is the real value of the system. The following figures are obtained from the error analysis in the Tables 1–4.

Table 2. Comparison between the approximate solution (

$V(x, t)$ ), versus an exact solution (

$v(x, t)$ ) for

$-10\leq x\leq 10$ and

$0\leq t\leq 1$ with k = 0.5 and c = 0.82 of Exp.2.

t	$V(\alpha=0.3)$	$V(\alpha=0.5)$	$V(\alpha=0.7)$	$V(\alpha=0.9)$	$V(\alpha=1)$	Exact	MAE
x=-10
0.00	0.010216	0.010216	0.010216	0.010216	0.010216	0.010216
0.25	0.0095819	0.00972	0.0098439	0.0099468	0.0099898	0.0097175
0.50	0.0094455	0.0095258	0.0096209	0.0097203	0.0097689	0.0092437	3.6330E-04
0.75	0.0093537	0.0093809	0.0094359	0.0095107	0.0095529	0.008793
1.00	0.0092825	0.0092616	0.0092729	0.0093128	0.0093417	0.0083642
x=-5
0.00	0.094515	0.094515	0.094515	0.094515	0.094515	0.094515
0.25	0.088823	0.090001	0.091122	0.092066	0.09246	0.090002
0.50	0.087702	0.088282	0.089091	0.08999	0.090436	0.085696	3.5E-03
0.75	0.086993	0.087053	0.087436	0.08807	0.088455	0.081589
1.00	0.086476	0.08609	0.086022	0.086289	0.08653	0.077671
x=0
0.00	0.447210	0.447210	0.447210	0.447210	0.447210	0.447210
0.25	0.445160	0.446220	0.446770	0.447030	0.44710	0.446660
0.50	0.444000	0.445150	0.446010	0.446560	0.446740	0.444990	1.46E-03
0.75	0.443020	0.444030	0.445030	0.445830	0.446130	0.442230
1.00	0.442140	0.442860	0.443860	0.444830	0.445250	0.438420
x=5
0.00	0.094515	0.094515	0.094515	0.094515	0.094515	0.094515
0.25	0.099243	0.099243	0.099243	0.099243	0.099243	0.030349
0.50	0.104190	0.104190	0.10419	0.030481	0.10419	0.104190	4.446E-03
0.75	0.109380	0.109380	0.109380	0.0305870	0.10938	0.109380
1.00	0.114800	0.114800	0.114800	0.0306870	0.1148	0.030682
x=10
0.00	0.010216	0.010216	0.010216	0.010216	0.010216	0.010216
0.25	0.010928	0.010752	0.010606	0.010492	0.010447	0.010739
0.50	0.011106	0.01117	0.011083	0.010978	0.010683	0.01129	4.3146E-04
0.75	0.01133	0.011332	0.011289	0.011215	0.010924	0.011868
1.00	0.114800	0.11480	0.114800	0.030687	0.011171	0.012477

| Show Table

DownLoad: CSV

Table 3. Comparison of exact solution (u(x, t)) with approximate (U(x, t)) solution for

$\alpha = 0.5$ , a = 0.2, b = 0.4, and k = 1 of Exp.1.

$x$	t	$\xi$ = $kx-4bk^3t$	Exact solution	ADM Solution	Error
-50	0.00	-49.9434	-0.35995	-0.35995	2.5009 $\times 10^{-08}$
-40	0.25	-39.6454	-0.35965	-0.35965	1.9820 $\times 10^{-07}$
-30	0.50	-29.6453	-0.35745	-0.35746	2.1868 $\times 10^{-06}$
-20	0.75	-19.6423	-0.34165	-0.34170	5.1926 $\times 10^{-05}$
$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$
$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$
$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$
$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$
$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$
20	0.25	20.3634	0.11987	0.11988	-1.6363 $\times 10^{-05}$
30	0.50	30.3653	0.12000	0.12000	-3.3767 $\times 10^{-07}$
40	0.75	40.3634	0.12000	0.12000	-6.2851 $\times 10^{-09}$
50	1.00	50.3643	0.12000	0.12000	-1.1537 $\times 10^{-10}$

| Show Table

DownLoad: CSV

Table 4. Comparison of exact (v(x, t)) solution with approximate (V(x, t) solution for k = 0.5, and c = 0.82 of Exp.2.

$x$	t	$\xi$ = $k+\sqrt{c}x-ct$	Exact solution	ADM Solution	Error
-50	0.00	-45.5969	-0.35995	-0.35995	-2.8783 $\times 10^{-07}$
-40	0.25	-36.5415	-0.35966	-0.35966	-2.1145 $\times 10^{-06}$
-30	0.50	-27.4862	-0.35751	-0.35750	-1.4960 $\times 10^{-05}$
-20	0.75	-18.4308	-0.34207	-0.34199	-7.7286 $\times 10^{-05}$
$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$
$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$
$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$
$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$
$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$	$\cdots$
20	0.25	17.7908	0.11986	0.11987	-1.3770 $\times 10^{-05}$
30	0.50	26.8462	0.12000	0.12000	-2.8188 $\times 10^{-07}$
40	0.75	35.9015	0.12000	0.12000	-5.2413 $\times 10^{-09}$
50	1.00	44.9569	0.12000	0.12000	-9.6194 $\times 10^{-11}$

| Show Table

DownLoad: CSV

The accuracy and effectiveness of the ADM are confirmed by using the error analysis given in Tables 1 and 2 followed by its MAE Figure 9(a, b) respectively. The other two Tables 3 and 4 confirm the errors that occur between the exact and approximate solution of the model. Their plots Figure 10(a, b) are shown, which are the graphical conformations for the error analysis between the two solutions (exact and approximate).

Figure 9. (a, b) 3D plots represent the mean absolute error (MAE) for different values of x and t for Tables 1 and 2 respectively.

DownLoad: Full-Size Img PowerPoint

Figure 10. (a, b) 3D plots represent the mean absolute error (MAE) for different values of x and t for Tables 3 and 4 respectively.

DownLoad: Full-Size Img PowerPoint

6. Conclusions

We consider time-fractional Korteweg-de Vries (TF-KdV) and time-fractional modified Korteweg-de Vries (TF-mKdV) equations for solitary wave solutions with the help of an analytical method called the "Aboodh decomposition method (ADM)." The proposed method is a combination of the Aboodh transform and the decomposition method, which is a very authentic and valuable method for solving fractional-order non-linear problems. The approximate solution obtained with the help of ADM is compared with the exact solution of the models to verify the efficiency of the method. We observed from the graphical analysis in the manuscript the time-fractional order $\alpha$ can modify the wave profiles in such a way that when $\alpha$ values is increasing the amplitudes of the wave is reducing which is very much important to study the small amplitude's characteristics of such wave equations. The effect of the nonlinearity coefficient "a" and the dispersive coefficient "b" in Exp.1 are also discussed, form the graphical analysis we have shown that both of these coefficients have significant affect on the wave amplitudes as well as its widths. We have studied the higher-order series solutions (up to fourth order) for time-fractional KdV and modified KdV equations with the help of the ADM. The convergence analysis determined that the obtained results are similar to the exact solutions of the models. The uniqueness results confirmed that ADM is an effective and systematic scheme for solving non-linear dynamical problems. The analytical and numerical solutions are plotted by marginal changes in parameters $x$ , $t$ and $\alpha$ . It is observed that a significant variation in the index values causes a large variation in the amplitude of the wave profile. For small-time $t$ , the series solutions coincide with the exact solutions of the models, and the solutions are stable. However, for large-time $t$ the results show that the amplitudes of the wave profiles no longer remain constant at the given parameters. The results are tabulated for altered $x$ , $t$ and $\alpha$ . It is observed that mean absolute error (MAE) is decreasing for small-time and large values of $\alpha$ (fractional order). This means that for taking $\alpha$ values nearest to 1 (exact), the MAE values decrease, which means that the obtained results coincide with the solution of the considered models.

7. Future work

We studied the analytical solution of the TF-KdV and TF-mKdV equations with the help of ADM. The effectiveness of the proposed method is studied with the help of graphical analysis and some numerical error tables. The method is also essential for $K(m, n)$ equation obtained in ^[40]. Here, the author shows that $K(m, n)$ equations are only Lax integrable for specific values of the parameters $(\alpha, m, n)$ . For the integrable instances, nontrivial prolongation structures and Lax pairings are provided.

Acknowledgments

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R8). Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflict of interest

The authors declare that they have no conflicts of interest.

References

[1]	Hancock M (2019) Ground improvement: slaying the sleech, Ground Engineering. Available from: https://www.geplus.co.uk/features/ground-improvement-slaying-sleech-20-11-2019/.
[2]	Crooks JHA (1973) Laboratory studies of the Belfast estuarine deposits. PhD thesis, Queens University Belfast.
[3]	Crooks JHA, Graham J (1972) Stress-strain properties of Belfast estuarine clay. Eng Geol 6: 275–288. https://doi.org/10.1016/0013-7952(72)90012-9 doi: 10.1016/0013-7952(72)90012-9
[4]	Crooks JHA, Graham J (1976) Geotechnical properties of the Belfast estuarine deposits. Géotechnique 26: 293–315. https://doi.org/10.1680/geot.1976.26.2.293 doi: 10.1680/geot.1976.26.2.293
[5]	Bell A (1977) Laboratory Studies of the Belfast Estuarine Deposits. PhD Thesis, Queen's University of Belfast.
[6]	Gregory BJ, Bell AL (1991) Geotechnical properties of Quaternary deposits in the Belfast area. Geol Soc Eng Geol Spec Publ 7: 219–228. https://doi.org/10.1144/GSL.ENG.1991.007.01.20 doi: 10.1144/GSL.ENG.1991.007.01.20
[7]	McCabe BA (2002) Experimental investigations of driven pile group behaviour in Belfast soft clay. PhD thesis, Trinity College Dublin.
[8]	Hight DW, Paul MA, Barras BF, et al. (2003) The characterization of Bothkennar Clay. Charact Eng Prop Nat Soils 1: 543–597.
[9]	Lehane BM, Jardine RJ, McCabe BA (2003) One-way axial cyclic tension loading of driven piles in clay. Proceedings of the BGA International Conference on Foundations, 493–506.
[10]	McCabe BA, Lehane BM (2003) Stress changes associated with driving pile groups in clayey silt. Proceedings of the 13^th European Conference on Soil Mechanics and Geotechnical Engineering, National University of Ireland Galway, 2: 271–276.
[11]	McCabe BA, Lehane BM (2006) Behavior of axially loaded pile groups driven in clayey silt. J Geotech Geoenviron Eng 132: 401–410. https://doi.org/10.1061/(ASCE)1090-0241(2006)132:3(401) doi: 10.1061/(ASCE)1090-0241(2006)132:3(401)
[12]	McCabe BA, Lehane BM (2006) Prediction of pile group response using a simplified non-linear finite element model, Methods in Geotechnical Engineering, CRC Press, Graz, 589–594.
[13]	McCabe BA, Lehane BM (2008) Measured and predicted t-z curves for a driven single pile in lightly overconsolidated silt. Proceedings of the 3^rd International Conference on Site Characterization, Taipei, 487–491.
[14]	McCabe BA, Phillips DT (2008) Design lessons from full-scale foundation load tests. Proceedings of the 3^rd International Conference on Site Characterization, Taipei, 615–620.
[15]	McCabe BA, Gavin KG, Kennelly ME (2008) Installation of a reduced-scale pile group in silt. 2^nd BGA International Conference on Foundations, IHS BRE Press, Dundee, 1: 607–616.
[16]	Doherty P (2010) Factors affecting the capacity of piles in clay. PhD thesis, University College Dublin.
[17]	Gavin KG, Gallagher D, Doherty P, et al. (2010) Field investigation of the effect of installation method on the shaft resistance of piles in clay. Can Geotech J 47: 730–741. https://doi.org/10.1139/T09-146 doi: 10.1139/T09-146
[18]	Doherty P, Gavin KG (2013) Pile aging in cohesive soils. J Geotech Geoenviron Eng 139: 1620–1624. https://doi.org/10.1061/(ASCE)GT.1943-5606.0000884 doi: 10.1061/(ASCE)GT.1943-5606.0000884
[19]	McCabe BA, Sexton BG, Lehane BM (2013) Operational coefficient of consolidation around a pile group in clay/silt. Geotech Geol Eng 31: 183–197. https://doi.org/10.1007/s10706-012-9579-1 doi: 10.1007/s10706-012-9579-1
[20]	Lehane BM, Phillips DT, Paul TS, et al. (1999) Instrumented piles subjected to combined vertical and lateral loads. Proceedings of the 12^th European Conference on Soil Mechanics and Geotechnical Engineering, Amsterdam, 2: 1145–1150.
[21]	Phillips DT (2002) Field tests on single piles subjected to lateral and combined axial and lateral loads. PhD Thesis, Trinity College Dublin.
[22]	Lehane BM (2003) Vertically loaded shallow foundation on soft clayey silt: A case history. Geotech Eng 156: 17–26. https://doi.org/10.1680/geng.2003.156.1.17 doi: 10.1680/geng.2003.156.1.17
[23]	Timoney MJ (2015) Strength verification methods for stabilised soil-cement columns: a laboratory investigation of the PORT and PIRT. PhD thesis, University of Galway.
[24]	Timoney MJ, McCabe BA (2017) Strength verification of stabilised soil-cement columns: a laboratory investigation of the Push-In Resistance Test (PIRT). Can Geotech J 54: 789–805. https://doi.org/10.1139/cgj-2016-0230 doi: 10.1139/cgj-2016-0230
[25]	Bowen DQ, Rose J, McCabe AM, et al. (1986) Correlation of Quaternary glaciations in England, Ireland, Scotland and Wales. Quat Sci Rev 5: 299–340. https://doi.org/10.1016/0277-3791(86)90194-0 doi: 10.1016/0277-3791(86)90194-0
[26]	Wilson HE (1972) Recent Geology of Northern Ireland, HMSO, Belfast.
[27]	Manning PI, Robbie JA, Wilson HE (1970) Geology of Belfast and the Lagan Valley, HMSO, Belfast.
[28]	Hight DW, Bond AJ, Legge JD (1992) Characterization of the Bothkennar clay: an overview. Géotechnique 42: 303–347. https://doi.org/10.1680/geot.1992.42.2.303 doi: 10.1680/geot.1992.42.2.303
[29]	Mitchell JK (1976) Fundamentals of Soil Behaviour, 1st edition, London: Wiley.
[30]	Lutenegger AJ, Cerato A (2001) Surface area and engineering properties of fine-grained soils. Proceedings of the 15^th International Conference on Soil Mechanics and Geotechnical Engineering, 1–3: 603–606.
[31]	Xu S, Feng Q, Deng Y, et al. (2024) Effect of diatom inclusion on physical property identification of soft clay. J Soils Sediments 24: 603–614. https://doi.org/10.1007/s11368-023-03673-x doi: 10.1007/s11368-023-03673-x
[32]	Gundersen AS, Hansen RC, Lunne T, et al. (2019) Characterization and engineering properties of the NGTS Onsøy soft clay site. AIMS Geosci 5: 665–703. https://doi.org/10.3934/geosci.2019.3.665 doi: 10.3934/geosci.2019.3.665
[33]	Shibuya S, Tamrakar SB, Thermast N (2001) Geotechnical site characterization on engineering properties of Bangkok Clay. J Southeast Asian Geotech Soc 32: 139–151.
[34]	Low HE, Landon Maynard M, Randolph MF, et al. (2011) Geotechnical characterisation and engineering properties of Burswood clay. Géotechnique 61: 575–591. https://doi.org/10.1680/geot.9.P.035 doi: 10.1680/geot.9.P.035
[35]	Zhang X, Liu X, Xu Y, et al. (2024) Compressibility, permeability and microstructure of fine-grained soils containing diatom microfossils. Géotechnique 74: 661–675. https://doi.org/10.1680/jgeot.22.00155 doi: 10.1680/jgeot.22.00155
[36]	Robertson PK (1990) Soil classification using the cone penetration test. Can Geotech J 27: 151–158. https://doi.org/10.1139/t90-014 doi: 10.1139/t90-014
[37]	Suzuki Y, Lehane BM (2014) Field observations of CPT penetration rate effects in Burswood clay. Proceedings of 3^rd International Conference on Cone Penetration Testing, 403–410.
[38]	Long M, Colreavy C, Ward D, et al. (2014) Piezoball tests in soft Irish clays. Proceedings of the 3^rd International Symposium on Cone Penetration Testing, Las Vegas, Nevada, 467–475.
[39]	Houlsby GT, Teh CI (1988) Analysis of the Piezocone in Clay. Conference on Penetration Testing, Balkema, Rotterdam, 777–783.
[40]	Leroueil S, Lerat P, Hight DW, et al. (1992) Hydraulic Conductivity of a Recent Estuarine Silty Clay at Bothkennar, Scotland. Géotechnique 42: 275–288. https://doi.org/10.1680/geot.1992.42.2.275 doi: 10.1680/geot.1992.42.2.275
[41]	Shiwakoti DR, Tanaka H, Tanaka M, et al. (2002) Influences of diatom microfossils on engineering properties of soils. Soils Found 42: 1–17. https://doi.org/10.3208/sandf.42.3_1 doi: 10.3208/sandf.42.3_1
[42]	McCabe BA, Sheil BB, Buggy F, et al. (2014) Empirical correlations for the compression index of Irish soft soils. Proc Inst Civ Eng Geotech Eng 167: 510–517. https://doi.org/10.1680/geng.13.00116 doi: 10.1680/geng.13.00116
[43]	Burland JB, Eng F (1990) On the Compressibility and Shear Strength of Natural Clays. Géotechnique 40: 329–378. https://doi.org/10.1680/geot.1990.40.3.329 doi: 10.1680/geot.1990.40.3.329
[44]	Mesri G, Godlewski PM (1977) Time and stress-compressibility inter-relationship. J Geotech Eng Div 103: 417–430. https://doi.org/10.1061/AJGEB6.0000421 doi: 10.1061/AJGEB6.0000421
[45]	Glossop NH, Saville DR, Moore JS, et al. (1979) Geotechnical aspects of shallow tunnel construction in Belfast estuarine deposits. Tunnelling 79, London, 45–56.
[46]	DeGroot DJ, Landon ME, Poirier SE (2019) Geology and engineering properties of sensitive Boston Blue Clay at Newbury, Massachusetts. AIMS Geosci 5: 412–447. https://doi.org/10.3934/geosci.2019.3.412 doi: 10.3934/geosci.2019.3.412
[47]	Chen C, Wang Y, Sun Z, et al. (2024) Rate effects of cylindrical cavity expansion in fine grained soil. J Rock Mech Geotech Eng, in press. https://doi.org/10.1016/j.jrmge.2024.09.032
[48]	Davies JA, Humpheson C (1981) Comparison between the performance of two types of vertical drain beneath a trial embankment in Belfast. Géotechnique 31: 19–31. https://doi.org/10.1680/geot.1981.31.1.19 doi: 10.1680/geot.1981.31.1.19
[49]	Jaky J (1944) The coefficient of earth pressure at rest. J Soc Hung Archit Eng, 355–358.
[50]	Ladd CC, Foott R, Ishihara K, et al. (1977) Stress-Deformation and Strength Characteristics. Proceedings of the 9th International Conference on Soil Mechanics and Foundation Engineering, Tokyo, 2: 421–494.
[51]	Padfield CJ, Mair RJ (1983) Design of retaining walls embedded in soft clay, CIRIA report 104.
[52]	Jardine RJ, Chow FC, Matsumoto T, et al. (1998) A new design procedure for driven piles and its applications to two Japanese clays. Soils Found 38: 207–219. https://doi.org/10.3208/sandf.38.207 doi: 10.3208/sandf.38.207
[53]	Lupini JF, Skinner AE, Vaughan PR (1981) Drained residual strength of cohesive soils. Géotechnique 31: 181–213. https://doi.org/10.1680/geot.1981.31.2.181 doi: 10.1680/geot.1981.31.2.181
[54]	Strick van Linschoten CJ (2004) Driven pile capacity in organic clays with particular reference to square piles in Belfast sleech. MSc thesis, Imperial College London.
[55]	Doherty P, Gavin KG (2009) Experimental investigation of the effect of shearing rate on the capacity of piles in soft silt. Contemp Top Deep Found, 575–582. https://doi.org/10.1061/41021(335)72 doi: 10.1061/41021(335)72

This article has been cited by:

Ahmed A. Gaber, Abdul-Majid Wazwaz, Dynamic wave solutions for (2+1)-dimensional DJKM equation in plasma physics, 2024, 9, 2473-6988, 6060, 10.3934/math.2024296

Reader Comments

Your name:*

Email:*
© 2025 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)