Next Article in Journal
Collocation-Based Approximation for a Time-Fractional Sub-Diffusion Model
Next Article in Special Issue
A Study on Generalized Degenerate Form of 2D Appell Polynomials via Fractional Operators
Previous Article in Journal
Global Dynamics of Fractional-Order Discrete Maps
Previous Article in Special Issue
Dynamical Behaviour, Control, and Boundedness of a Fractional-Order Chaotic System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Physical Phenomenon for the Fractional Nonlinear Mixed Integro-Differential Equation Using a Quadrature Nystrom Method

1
Department of Mathematics, Faculty of Applied Science, Umm Al–Qura University, Makkah P.O. Box 715, Saudi Arabia
2
Department of Mathematics, Faculty of Education, Alexandria University, Alexandria 21544, Egypt
3
Department of Mathematics, Faculty of Engineering, Sinai University, El Arish 16020, Egypt
*
Author to whom correspondence should be addressed.
Submission received: 19 July 2023 / Revised: 22 August 2023 / Accepted: 23 August 2023 / Published: 31 August 2023
(This article belongs to the Special Issue Advances in Fractional Modeling and Computation)

Abstract

:
In this work, the existence and uniqueness solution of the fractional nonlinear mixed integro-differential equation (FrNMIoDE) is guaranteed with a general discontinuous kernel based on position and time-space   L 2 Ω × C 0 , T ,   T < 1 . The FrNMIoDE conformed to the Volterra-Hammerstein integral equation (V-HIE) of the second kind, after applying the characteristics of a fractional integral, with a general discontinuous kernel in position for the Hammerstein integral term and a continuous kernel in time to the Volterra integral (VI) term. Then, using a separation technique methodology, we developed HIE, whose physical coefficients were time-variable. By examining the system’s convergence, the product Nystrom technique (PNT) and associated schemes were employed to create a nonlinear algebraic system (NAS).

1. Introduction

Due to the vast number of applications that can be found for fractional nonlinear/linear integro-differential equations FrNIoEs/FrLIoEs containing time-dependent coefficients in physics, engineering, and other scientific domains, their significance has been growing steadily over the past several decades. These equations are ideal for accurately representing a variety of events that occur in the real world because they capture both the non-local and local behavior of a large number of complex systems. Hermann [1] introduced some applications of fractional calculus to the field of physics. Oyedepo et al. [2] presented a numerical solution to the linear FrLIoDE problem by employing the method of standard least squares. Bernstein piecewise polynomials were exploited by Osama and Sarmad [3] in order to find an approximate solution to the FrLIoDE problem. In Daşcıoğlu and Bayram [4], approximate solutions to FrLIoDEs were found by using Laguerre polynomials. Mohammed [5] utilized the approach of least squares in conjunction with a shifted Chebyshev polynomial in order to solve the FrLIoDE problem. In [6], Mahdy et al. investigated the numerical solution of FrLIoDE by employing the least squares approach and supplementing it with a shifted Laguerre polynomial. In order to locate the numerical solution of FrLIoDE with the Caputo derivative, Nanware et al. [7] used the Bernstein polynomial to solve it. The least square technique and the homotopy perturbation method were both proposed by Oyedepo et al. [8] to discuss the solution to the FrIDE problem. Several techniques are used in efficient ways to solve IE and FrIE. For instance, Basseem and Alayani [9] solved a nonlinear quadratic mixed IE of the second kind with a singular kernel by employing the Toeplitz matrix method in conjunction with the PNT. A quadrature scheme was implemented by I. Katani [10] for the numerical outcomes of the second type of the Fredholm integral (FI) model. Al-Bugami [11] employed the Simpson and Trapezoidal methods to perform numerical representations based on an integral model that utilized 2D surface crack layers. In order to obtain numerical computing results for the second kind of nonlinear integral model that has a continuous kernel, Brezinski and Zalglia [12] employed the extrapolation method. Baksheesh [13] suggested making use of the Galerkin scheme in order to find the approximate solution for the VIEs of the second kind. The sinc-collocation method was utilized by Alkan and Hatipoglu [14] to find solutions of V-FIDEs of fractional order. Mosa et al. [15] researched the semigroup scheme to evaluate uniqueness and existence based on the partial and fractional integro models of heat performance in the Banach space using the Adomian decomposition scheme. An effective methodology for finding an approximate solution using the wavelet collocation method to FrFIDEs was proposed by Bin Jebreen and Dassios [16]. Using the extended cubic B-spline, Akram et al. [17] interpreted the collocation strategy in order to solve the fractional partial integro-differential problem. After employing the Riemann–Liouville fractional integral and fractional derivative, Abdelkawy et al. [18] used the Jacobi–Gauss collocation method to achieve an approximate solution for a variable-order of FrLIoDE with a weakly singular kernel. This was accomplished by applying the Jacobi–Gauss collocation method. Khalil [19] used a Jacobi polynomial for a solution of coupled system of fractional differential equations. Some different numerical methods for NLIoEs are introduced in Abdou [20].
This paper is divided into 12 sections. In Section 2, we derive the fractional mixed integral equation from the phase lag mixed integral equation. In Section 3, using a mixed NIE of the second kind with local conditions, we establish FrNIoDE, using the Caputo-fractional integral. In addition, the existence of a unique solution is guaranteed. In Section 4, the convergence of this solution is established and proved. In Section 5, the technique of separation variables method is applied to change the problem to HNIE of the second kind where its coefficients are parameters of time. PNT is employed in Section 6 to obtain the NAS. Then, we discuss the existence and uniqueness solution of NAS in Section 7. While the convergence system is considered in Section 8, the estimated error of PNT is discussed in Section 9. Illustrative numerical examples are involved to demonstrate the propriety and effectuality of the technique and some conclusions are stated in Section 10 and Section 11. In Section 12, we propose some parameters for our future work.

2. Time Fractional and Phase Lag Integral Equation

The integral equations play an important role in the phase-lag problems with local conditions. Consider, in the time fractional calculus, the phase lag integral equation:
φ x , t + q = f x , t + λ 0 t Ω k x y G t , τ γ y , τ , φ y , τ d y d τ , 0 q 1 ,
with conditions
φ x , 0 = V 1 x , t φ x , t t = 0 = V 2 x ,
where f ( x ,   t ) is a known function in the space L 2 Ω × C 0 , T ,   T < 1 represents the free term of the problem, λ is a constant that depends on the kind of material (in applied mathematics and has many physical meanings). The function G t , τ is a smooth kernel in time, while k x y is a singular kernel in position, which will be taken as a logarithmic form and Carleman function, γ x , t , φ x , t is a known nonlinear function of the unknown function   φ x , t , V 1 x and V 2 x are two given initial position functions, the constant q is a small quantity that represents the delay of time.
Here, the aim of this research is to predict the near future, by studying fractional derivatives, and using the initial conditions, where it is known that the differential derivatives express the breaking of the ionic bond between the particles of the substance, and that the use of fractional time enables the researcher to deepen this study. The past time is studied when q is negative.
Abdou and Raad [21] and Mosal et al. [22] discussed the solution of mixed IE with nonlocal conditions. However, in this research, the above is developed by studying the fractional delay over time for a local phase-lag problem of a FrNMIoDE with continuous kernel in time and singular kernel in position.
Using Taylor’s expansion, in the fractional calculus, to have
φ x , t + q φ x , t + q α Γ α + 1 α t α φ x , t + q α + 1 Γ α + 2 α + 1 t α + 1 φ x , t + q α + 2 Γ α + 3 α + 1 t α + 1 φ x , t + , ,       n 1 < α < n .
In this work, we have focused on n = 2
φ x , t + q = φ x , t + q α Γ α + 1 α t α φ x , t + q α + 1 Γ α + 2 α + 1 t α + 1 φ x , t .
Using the basic formula of the Caputo-fractional integral
I a α f x = 1 Γ α a x x t α 1 f t d t ,
I a α I a β f x = I a α + β f x = 1 Γ α + β a x x t α + β 1 f t d t ,
and
0 t 0 τ n 1 0 τ 2 0 τ 1 f τ d τ d τ 1 d τ n 2 d τ n 1 = 1 Γ ( n ) 0 t t τ n 1 f τ d τ .
In the view of Equations (1) and (4), we have
q α + 1 Γ α + 2 α t α φ x , t t φ x , 0 + q α Γ α + 1 I 1 α φ x , t + 0 t φ x , τ d τ λ 0 t 0 τ 1 Ω G t , τ k x y γ y , τ , φ y , τ d y d τ 1 d τ = 0 t f x , τ d τ .
Applying Equations (5)–(7) in Equation (8), we obtain
μ 1 φ x , t + μ 2 0 t 1 + α q α t τ α φ x , τ d τ + μ 3 0 t Ω t τ α + 1 G t , τ k x y γ y , τ , φ y , τ d y d τ = F x , t ,
where
μ 1 = q α + 1 Γ α + 2 , μ 2 = q α Γ α + 1 , μ 3 = λ 2 Γ α F x , t = 1 Γ α 0 t t τ α f x , τ d τ + q α + 1 V 1 x Γ α + 2 Γ α + 1 t α + q α + 1 Γ α + 2 V 2 x .
In order to guarantee the existence of a unique solution of the considered problem of Equation (1) or its equivalent in Equation (9), we assume the following conditions:
(i)
The unknown function   φ x , t and its derivatives are in the space L 2 Ω × C 0 , T and its norm is defined as
φ x , t = max 0 t T 0 t Ω φ 2 x , τ d x 1 / 2 d τ .
(ii)
For the constant Q , the known function γ x , t , φ satisfies the following conditions:
γ x , t , φ Q φ x , t . γ x , t , φ 1 γ x , t , φ 2 Q φ 1 x , t φ 2 x , t
(iii)
The given function f ( x , t ) satisfies
f x , t = max 0 t T 0 t Ω f 2 x , τ d x 1 / 2 d τ B , B   i s   a   c o n s t a n t .
(iv)
The two functions V i x for the constants D i   ,   i = 1 ,   2 satisfy the following:
V 1 x D 1   and   V 2 x D 2 .
(v)
The position kernel in the space L 2 Ω satisfies
k x y = Ω Ω k 2 x y d x d y 1 / 2 = C , C   is a constant .
(vi)
The continuous function G ( t , τ ) in time satisfies
m a x 0 t T G t , τ = E ,     E   is a constant .

3. Existence and Uniqueness

To prove the existence and uniqueness of Equation (9), it can be written in the following integral operator form
χ φ x , t = χ 1 φ x , t χ 2 φ x , t + V 2 x + V 1 x Γ α 0 t t τ 1 d τ + α α + 1 q α + 1 0 t t τ α f x , τ d τ . χ 1 φ x , t = α + 1 q 0 t φ x , τ d τ + α α + 1 q α + 1 0 t t τ α φ x , τ d τ χ 2 φ x , t = λ α α + 1 2 q α + 1 0 t Ω t τ α + 1 G t , τ k x y γ y , τ , φ y , τ d y d τ .
Theorem 1.
The solution of Equation (9) exists and is unique under the condition:
α + 1 T q α + α T α + 1 + λ α α + 1 T α + 2 2 α + 2 E C Q < q α + 1 .
The following lemmas must be proved to satisfy the above theorem.
Lemma 1.
The operator χ   maps the space L 2 Ω × C 0 , T onto itself under the conditions (i)–(vi).
Proof. 
From Equation (10), we obtain
χ φ x , t χ 1 φ x , t + χ 2 φ x , t + V 2 x + V 1 x Γ α 0 t t τ 1 d τ + α α + 1 q α + 1 0 t t τ α f x , τ d τ
Using conditions (i)(v) and Cauchy–Schwartz inequality, we have
χ φ x , t δ φ x , t + D 1 T α Γ α + 1 + D 2 + α B T α + 1 q α + 1
where
δ = α + 1 T q α + α T α + 1 + λ α α + 1 T α + 2 2 α + 2 E C Q q α + 1 .
It is obvious that the operator χ maps the ball B r     L 2 1 ,   1 × C 0 ,   T Š onto itself where
r = σ 1 δ ,   σ = D 1 T α Γ α + 1 + D 2 + α B T α + 1 q α + 1 .
The inequality (13) involves the boundedness of the operator χ under the condition δ < 1. □
In the previous Lemma, we considered that the discontinuous kernel t τ 1 , f o r   a l l   t , τ   [ 0 , T ] , satisfies for every continuous function h t , τ , h t , τ c o n s t a n t s , t h e   f o l l o w i n g 0 t t τ α 1 h t , τ d τ or t 1 t 2 t τ α 1 h t , τ d τ , 0 < t 1 t 2 t , max t 0 t t τ α 1 d τ , e x i s t s .
Lemma 2.
If conditions (i)–(vi) are satisfied, then χ   is a contraction operator in Banach space L 2 1 ,   1 × C 0 ,   T .
Proof. 
Let two functions φ 1 x , t and φ 2 x , t be two solutions of (9), and then, the Formula (10) leads to the
χ φ 1 x , t χ φ 2 x , t χ 1 φ 1 x , t φ 2 x , t + χ 2 φ 1 x , t φ 2 x , t .
Using conditions (i)–(iv) and Cauchy–Schwarz inequality, we deduce that
χ φ 1 x , t χ φ 2 x , t δ φ 1 x , t φ 2 x , t
It follows that for δ < 1, χ is a contraction operator of Equation (10). Hence, there exists a unique solution in L 2 Ω by a Banach fixed point theorem for every     t     C 0 ,   T ,   T < 1 . □

4. Convergence of Solution

For this aim, take the straightforward iteration φ 1 x ,   t ,   φ 2 x ,   t ,   . . . ,   φ n x ,   t ,   . . .   φ x ,   t . Then, use Equation (9), to have
q α + 1 Γ α + 2 φ n x , t φ n 1 x , t + q α Γ α + 1 0 t φ n 1 x , τ φ n 2 x , τ d τ   + 1 Γ α 0 t t τ α φ n 1 x , τ φ n 2 x , τ d τ = λ 2 Γ α 0 t Ω t τ α + 1 G t , τ k x y γ n y , τ , φ n 1 y , τ                                                                                           γ n y , τ , φ n 2 y , τ d y d τ ,
let
φ n x , t = i = 0 n ψ i x , t ,
where
ψ n x , t = φ n x , t φ n 1 x , t , ψ 0 x , t = 1 Γ α 0 t t τ α f x , τ d τ . n 1 .
Lemma 3.
A sequence φ n x ,   t of Equation (17) is uniformly convergent under the condition δ < 1 .
Proof. 
By applying Cauchy–Schwarz and using (17) in (16), we obtain
q α + 1 Γ α + 2 ψ n x , t q α Γ α + 1 0 t d τ ψ n 1 x , t + 1 Γ α 0 t t τ α d τ × ψ n 1 x , t + λ N 2 Γ α 0 t Ω t τ α + 1 G t , τ k x y d y d τ ψ n 1 x , t
Taking n = 1 , the above formula becomes
ψ n x , t δ T α + 1 B Γ α + 1 ,
and then
ψ n x , t δ n T α + 1 B Γ α + 1 , δ < 1 .
The sequence ψ n x , t is uniformly convergent by Equation (19). Additionally, it gives the sequence’s φ n x , t = i = 0 n ψ i x , t convergent solution.
As n ,   φ n x , t φ x , t , hence the solution φ x , t is uniformly convergent under the condition δ < 1 . This demonstrates the lemma. □

5. Separation of Variables Technique

We discover that researchers are drawn to the unknown potential function, which is connected to time and place, in the issues of mathematical physics. The unknown function can be obtained using a variety of approaches. Time division is one of these techniques, converting the mixed integral problem into an algebraic system of integral equations.
Some researchers use the separation technique method which is a powerful mathematical tool that allows us to transform FrNIoE with time-dependence into a class of integral equations with coefficients on time only. This technique simplifies the problem by separating the time-dependent part from the integral part, enabling us to handle the integral equations more efficiently. The unknown and well-known functions are shown in the separation form as
                              φ x , t = N t ψ x , f x , t = g x M t , γ x , t , φ x , t = γ 1 t , N ( t ) γ 2 x , ψ ( x ) .
Hence, after using (20), the Formula (9) yields,
μ 1 t ψ x μ 2 t Ω k x y γ 2 ψ y , y d y = μ 3 t g x + q α + 1 V 1 x Γ α 0 t t τ 1 d τ + q α + 1 V 2 x ,
where
μ 1 t = q α + 1 N t + α + 1 q α 0 t N τ d τ + α α + 1 0 t t τ α N τ d τ . μ 2 t = λ 2 α α + 1 0 t t τ + 1 G t , τ γ 1 τ , N ( τ ) d τ , μ 3 t = α α + 1 0 t t τ M τ d τ .
Equation (21) considers HIE of the second kind with coefficients specifying the time domain μ i t ,   i = 1 ,   2 ,   3 . Here, μ 1 t and μ 2 t indicate the time term of the unknown function   φ x , t . Meanwhile, μ 3 t describes the time in the free term   f x , t .

6. Product Nystrom Technique (PNT)

To solve integral equations with continuous or disconnected kernels, several numerical approaches have been utilized. The PNT is the best approach for solving singular integral equations for the following reasons: the singular term vanishes instantly, turned into simple integrals that can be solved rapidly, and then creates a NAS. The relative error approach has a lower degree of convergence than the other methods.
Here, by using the product integration, we approximate the integral part of Equation (21) by a suitable Lagrange interpolation polynomial. For this, let x = x m and the integral part can be written as
Ω k x m y γ 2 ψ y , y d y = n = 0 N I m , n γ 2 ψ y n , y n = n = 0 N 2 2 y 2 n y 2 n + 2 k x m y γ 2 ψ y , y d y .
Without sacrificing generality, we take Ω = 1 ,   1 , x m = y m = 1 + m h , m = 1,2 , 3 , , N where h = 2 N and N is an even number.
The Lagrange interpolation polynomial is used to approximate the nonsingular component of the integral throughout each interval [ y 2 n , y 2 n + 2 ] at the points 2n, 2n + 1 and 2n + 2. As a result, the integral term of (23) becomes
1 1 k x m y γ 2 ψ y , y d y = n = 0 N 2 2 y 2 n y 2 n + 2 k x m y y 2 n + 1 y y 2 n + 2 y y 2 n + 1 y 2 n y 2 n + 2 y 2 n γ 2 ψ y 2 n , y 2 n + y 2 n y y 2 n + 2 y y 2 n y 2 n + 1 y 2 n y 2 n + 1 γ 2 ψ y 2 n + 1 , y 2 n + 1 + y 2 n y y 2 n + 1 y y 2 n y 2 n + 2 y 2 n + 1 y 2 n + 2 γ 2 ψ y 2 n + 2 , y 2 n + 2 d y .
Comparing Equations (23) and (24), we deduce
I m , 0 = 1 2 h 2 y 0 y 2 k x m y y 1 y y 2 y d y , I m , 2 n + 1 = 1 h 2 y 2 n y 2 n + 2 k x m y y y 2 n y 2 n + 2 y d y , I m , 2 n = 1 2 h 2 y 2 n y 2 n + 2 k x m y y 2 n + 1 y y 2 n + 2 y + y 2 n 2 y 2 n k x m y y y 2 n 1 y y 2 n 2 d y , I m , N = 1 2 h 2 y N 2 y N k x m y y y N 2 y y N 1 d y .
Introduce the following notations
α n y m = 1 2 h 2 y 2 n 2 y 2 n k x m y y y 2 n 2 y y 2 n 1 d y ,
β n y m = 1 2 h 2 y 2 n 2 y 2 n k x m y y 2 n 1 y y 2 n y d y ,
and
ζ n y m = 1 2 h 2 y 2 n 2 y 2 n k x m y y y 2 n 2 y 2 n y d y ,
then
I m , 0 = β 1 x m , I m , 2 n + 1 = 2 ζ n + 1 x m , I m , 2 n = α n x m + β n + 1 ( x m ) , I m , N = α N / 2 x m .
By substituting in Equation (21), we obtain
μ 1 t ψ x m μ 2 t n = 0 N I m , n γ 2 ψ y n , y n = μ 3 t g x m                                                     + q α + 1 Γ α + 1 { V 1 x m t α + V 2 x m Γ α + 1 } .
Equation (28) represents the NAS which gives an approximate solution of Equation (4) in a certain domain of time.

7. The Existence of a Unique Solution of NAS

To prove the existence of a unique solution of Equation (28), we write it in the following operator form
χ ¯ ψ x m = μ 2 t μ 1 t n = 0 N I m , n γ 2 ψ y n , y n + μ 3 t μ 1 t g x m + q α + 1 μ 1 t Γ α + 1 { V 1 x m t α + V 2 x m Γ α + 1 } ,
where the following assumptions are held:
(a)
The parameters μ i t satisfy max t μ i t A i t 0 , T , T < 1 , i = 1 , 2 , 3 where A i are constants.
(b)
g x m l 2 = m = 0 g x m 2 1 2 B ¯ .
(c)
sup m V i x m D i ¯ , i = 1 , 2 where D i ¯ are constants.
(d)
I m , n l 2 = n = 0 m = 0 I 2 x m x n 1 2 C ¯ .
(e)
The unknown function   ψ x m is in the space l 2 and its norm is defined as
ψ x m = m = 0 ψ x m 2 1 2
(f)
For the constant Q ¯ , the known function γ 2 x m , ψ x m satisfies the following conditions:
γ 2 x m , ψ x m Q ¯ ψ x m
and
γ 2 x m , ψ 1 x m γ 2 x m , ψ 2 x m Q ¯ ψ 1 x m ψ 2 x m .
Theorem 2.
The approximate solution of the NAS of Equation (28) exists and unique under the condition
A 2 C ¯ Q ¯ < A 1 .
The following lemmas must be proved to satisfy the above theorem.
Lemma 4.
The operator χ ¯   maps the space l 2 onto itself under the conditions (a)–(f).
Proof. 
From Equation (29), applying Cauchy–Minkowski inequality, we obtain
χ ¯ ψ x m μ 2 t μ 1 t γ 2 ψ y m , y m n = 0 m = 0 I 2 x m x n 1 2 + μ 3 t μ 1 t g x m + q α + 1 μ 1 t Γ α + 1 { V 1 x m T α + Γ α + 1 V 2 x m ,
using the conditions (a)–(f), we obtain
χ ¯ ψ x m A 2 A 1 Q ¯ ψ x m C ¯ + A 3 A 1 B + q α + 1 A 1 Γ α + 1 [ D 1 T α + D 2 Γ α + 1 ] ,
so, we have
χ ¯ ψ x m δ ψ x m + σ
where
δ = A 2 A 1 Q ¯ C ¯ .
It is obvious that the operator χ ¯ maps the ball B r *     l 2 × C[0, T] onto itself where
r = σ 1 δ ,   σ = A 3 A 1 B + q α + 1 A 1 Γ α + 1 [ D 1 T α + D 2 Γ α + 1 ] .
The inequality (33) involves the boundedness of the operator χ ¯ under the assumption δ * < 1 . □
Lemma 5.
If conditions (a)–(f) are satisfied, then χ ¯   is a contraction operator in Banach space l 2 .
Proof. 
Let two functions ψ 1 x m and ψ 2 x m be two solutions of (29), and then, the Formula (31) leads to
χ ¯ ψ 1 x m χ ¯ ψ 2 x m μ 2 t μ 1 t n = 0 N I m , n γ 2 x m , ψ 1 x m γ 2 x m , ψ 2 x m .
Using conditions (a)–(f) and Cauchy–Minkowski inequality, we obtain
χ ¯ ψ 1 x m χ ¯ ψ 2 x m δ ψ 1 x m ψ 2 x m .
It follows that for δ * < 1, χ ¯ is a contraction operator of a system (36). Hence, there exists a unique solution in l 2 by a Banach fixed point theorem for every   t     C 0 ,   T .

8. Convergence of the Approximate Solution of NAS

To discuss the convergence of the system (28), we state the following theorem.
Theorem 3.
The NAS (28) for all values of time t 0 ,   T ,   T < 1 is convergent in the Banach space l 2 under the condition δ * < 1 .
Proof. 
We construct the sequence ψ z x m , and using Equation (29), we have
μ 1 t ψ z x m ψ z 1 x m = μ 2 t n = 0 N I m , n γ 2 ψ z 1 y n , y n γ 2 ψ z 2 y n , y n ,
consider
ψ z x m = s = 0 z η s x m ,
where
η z x m = ψ z x m ψ z 1 x m , η 0 x m = μ 3 t g x m . z 1 .
By applying Cauchy–Schwarz and using (38) in (37), we obtain
A 1 η z x m A 2 C ¯ Q ¯ ψ n 1 x , t .
Taking n = 1 , the above formula becomes
η 1 x m δ A 3 B ¯ ,
and then
η z x m δ z A 3 B ¯ , δ < 1 .
The sequence η z x m is uniformly convergent by Equation (39). Additionally, it gives the sequence’s ψ z x m = s = 0 z η z x m convergent solution.
As z ,   ψ z x m ψ x m , hence the solution ψ x m is uniformly convergent under the condition δ * < 1 . This demonstrates the lemma. □

9. The Error of the Product Nystrom Technique

The following two definitions are used to calculate the error of this technique:
Definition 1.
The Nystrom method is said to be convergent of order r in the interval [−1, 1], if and only if, for N sufficient large, there exists a constant K > 0 independent of N such that
ψ x ψ N ( x ) K N r
Definition 2.
The estimated error of this method can be calculated in the form
R N =   1 1 k x m y γ 2 ψ y , y d y n = 0 N I m , n γ 2 ψ y n , y n ,
As N , R N 0 . In this case, the approximate solution of (28) is equivalent to the exact solution of (21) in the space L 2 1 , 1 × C 0 , T ,   T < 1 .

10. Numerical Results

In this section, some numerical applications are considered to show the accuracy and applicable of the proposed methods.
Example 1.
φ x , t + 0.03 = f x , t + 0.01 0 t 1 1 k ( x y ) τ 3 t 3 φ 2 y , τ d y d τ ,
with conditions
φ x , 0 = x , φ x , 0 t = 0.5 x 2 .
Then, by taking
φ x , t + 0.03 = φ x , t + 0.03 0.75 Γ 1.75 0.75 t 0.75 φ x , t + 0.03 1.75 Γ 2.75 1.75 t 1.75 φ x , t ,
and k x y = x y υ , we obtain
0.03 1.75 Γ 2.75 φ x , t 0.5 x 2 0.03 1.75 x Γ 2.75 Γ 0.75 0 t t τ 0.25 d τ + 0.03 0.75 Γ 1.75 0 t φ x , τ d τ + 1 Γ 0.75 0 t t τ α φ x , τ d τ λ 2 Γ 0.75 0 t 1 1 t τ 1.75 x y υ τ 3 t 3 φ 2 y , τ d y d τ = 1 Γ 0.75 0 t t τ α f x , τ d τ ,
where f ( x , t ) is given by putting φ x , t = X x N t , X = x 2 , N ( t ) = t + 0.5 as an exact value. We have the following two cases:
Case (i): if f ( x , t ) and the unknown function φ x , t have the same function of time where M ( t ) = t + 0.5 .
Case (ii): if f ( x , t ) and the unknown function φ x , t have different time function in which M t = 0.7 e t .
The absolute errors in case (i) with some values of x 1 ,   1 are shown in Table 1 and the approximate solution φ x , t = ( t + 0.5 ) x 2 is represented by Figure 1 while the errors for different T are shown by Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6. The errors in case (ii) are shown in Table 2 and Figure 7 and Figure 8.
The error increases with time and decreases by increasing the number of iterations, see Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8.
Example 2.
Consider Equation (41) with k x y = ln x y , we obtain
0.03 1.75 Γ 2.75 φ x , t 0.5 x 2 0.03 1.75 x Γ 2.75 Γ 0.75 0 t t τ 0.25 d τ + 0.03 0.75 Γ 1.75 0 t φ x , τ d τ + 1 Γ 0.75 0 t t τ α φ x , τ d τ λ 2 Γ 0.75 0 t 1 1 t τ 1.75 ln x y τ 3 t 3 φ 2 y , τ d y d τ = 1 Γ 0.75 0 t t τ α f x , τ d τ .
Here, f x , t = 0.1 + sin t   g ( x ) is given by putting φ x , t = X x N t , X = x 3 , N t = 0.3 cos t as an exact value.
The rate of errors is evaluated using the following formula: R a t e = log 2 E r r o r 2 N E r r o r N .
By increasing N, the error decreases and the rate of convergence is given below, see Table 3, Table 4, Table 5 and Table 6.
The approximate solution is shown in Figure 9, and the absolute errors with different time are shown in Figure 10 and Figure 11.

11. Conclusions

Fractional calculus has proved to be a valuable tool for modeling and analyzing numerous phenomena with non-local and memory-dependent characteristics. In this research, we focus on the study of fractional nonlinear integral equations with a time-dependent coefficient using the separation technique method. By employing this approach, we aim to transform the integral equation into a class of integral equations with coefficients on time, which can be subsequently solved using the Nystrom method. From the previous work and discussion, we can establish the following:
  • The separation technique method is a powerful mathematical tool that allows us to transform fractional nonlinear integral equations with time-dependent coefficients into a class of integral equations with coefficients on time only. This technique simplifies the problem by separating the time-dependent part from the integral part, enabling us to handle the integral equations more efficiently.
  • We present the Nystrom method as a numerical scheme to solve the sub-equations derived from the separation process. The Nystrom method efficiently discretizes the integral equation using a set of points and exploits the smoothness of the kernel functions to achieve high accuracy and computational efficiency. We discuss the convergence analysis and error estimates associated with the Nystrom method for the proposed fractional integral equations.
  • In the Carleman kernel, when decreasing   ν , the solution becomes better, see Table 1.
  • The error increases with time, see Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8.
  • The error decreases by increasing the number of iterations, see Table 3, Table 4, Table 5 and Table 6.
  • CPU Time for Intel® Core(TM) i7 CPU M 620 @ 2.67 GHz (64-bit Operating system) in example 2 takes the following:
NTime in SecondsMemory
410.459 M.
814.34219.19 M.
1634.32215.9 M.
  • The error is the same to an extreme degree between positive and negative values for x in the region of integration. We note that the maximum error occurs at the ends of x values while it takes the minimum in the middle, see Figure 2.

12. Future Work

We will attempt to solve Equation (1) as a phase-lag problem with a history memory function in the sense of general form of fractional calculus and differ between the three following cases of time: terrestrial ( T < 1 ), eternal ( T > 1 ), and isthmus lives ( T = 1 ).

Author Contributions

Conceptualization, M.A.A.; Methodology, A.R.J. and M.B.; Software, M.B.; Validation, M.A.A.; Formal analysis, A.R.J., M.A.A. and M.B.; Writing—original draft, A.R.J. and M.B.; Writing—review & editing, M.A.A. and M.B.; Supervision, M.A.A.; Funding acquisition, A.R.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No instructional records have been created or comfortable with the ongoing assessment data in the original copy.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Herrmann, R. Fractional Calculus: An Introduction for Physicists; World Scientific: Singapore, 2014. [Google Scholar]
  2. Oyedepo, T.; Taiwo, O.A. Numerical Studies for Solving Linear Fractional Integro-differential Equations Based on Constructed Orthogonal Polynomials. ATBU J. Sci. Technol. Educ. 2019, 7, 1–13. [Google Scholar]
  3. Mohammed, O.H.; Altaie, S.A. Approximate solution of Fractional Integro-Differential Equations by using Bernstein polynomials. Eng. Technol. J. 2012, 30, 1362–1373. [Google Scholar] [CrossRef]
  4. Daşcıoğlu, A.; Bayram, D.V. Solving Fractional Fredholm Integro-Differential Equations by Laguerre Polynomials. Sains Malays. 2019, 48, 251–257. [Google Scholar] [CrossRef]
  5. Mohammed, D.S. Numerical solution of Fractional Integro-differential Equations by least squares method and shifted Chebyshev polynomial. Comp. Methods Eng. Sci. 2014, 2014, 431965. [Google Scholar] [CrossRef]
  6. Mahdy, A.M.S.; Mohamed, E.M.H.; Marai, G.M.A. Numerical solution of fractional integro-differential Equations by least squares method and shifted Chebyshev polynomials of the third kind method. Theor. Math. Appl. 2016, 6, 87–101. [Google Scholar]
  7. Nanware, J.A.; Goud, P.M.; Holambe, T.L. Solution of Fractional Integro-differential Equations by Bernstein Polynomials. Malaya. J. Mat. 2020, 1, 581–586. [Google Scholar] [CrossRef]
  8. Oyedepo, T.; Taiwo, O.A.; Adebisi, A.F.; Ishola, C.Y.; Faniyi, O.E. Least Squares Method and Homotopy Perturbation Method for Solving Fractional Integro-Differential Equations. Pac. J. Sci. Technol. 2019, 20, 86–95. [Google Scholar]
  9. Basseem, M.; Alalyani, A. On the Solution of Quadratic Nonlinear Integral Equation with Different Singular Kernels. Math. Probl. Eng. 2020, 2020, 7856207. [Google Scholar] [CrossRef]
  10. Katani, R. Numerical solution of the Fredholm integral Equations with a quadrature method. SeMA J. 2018, 76, 271–276. [Google Scholar] [CrossRef]
  11. Al-Bugami, A.M. Numerical Treating of Mixed Integral Equation Two-Dimensional in Surface Cracks in Finite Layers of Materials. Adv. Math. Phys. 2022, 2022, 3398175. [Google Scholar] [CrossRef]
  12. Brezinski, C.; Redivo-Zaglia, M. Extrapolation methods for the numerical solution of nonlinear Fredholm integral Equations. J. Integral Equ. Appl. 2019, 31, 29–57. [Google Scholar] [CrossRef]
  13. Bakhshayesh, S.J. Discontinuous Galerkin Approximations for Volterra Integral Equations of the First Kind with Convolution Kernel. Indian J. Sci. Technol. 2015, 8, 33. [Google Scholar] [CrossRef]
  14. Alkan, S.; Hatipoglu, V.F. Approximate solutions of Volterra-Fredholm integro-differential Equations of fractional order. Tbil. Math. J. 2017, 10, 1–13. [Google Scholar] [CrossRef]
  15. Mosa, G.A.; Abdou, M.A.; Gawish, F.A.; Abdalla, M.H. On the behavior solutions of fractional and partial integro differential heat Equations and its numerical solutions. Math. Slovak 2022, 72, 397–410. [Google Scholar] [CrossRef]
  16. Jebreen, H.B.; Dassios, I. On the Wavelet Collocation Method for Solving Fractional Fredholm Integro-Differential Equations. Mathematics 2022, 10, 1272. [Google Scholar] [CrossRef]
  17. Akram, T.; Ali, Z.; Rabiei, F.; Shah, K.; Kumam, P. A Numerical Study of Nonlinear Fractional Order Partial Integro-Differential Equation with a Weakly Singular Kernel. Fractal Fract. 2021, 5, 85. [Google Scholar] [CrossRef]
  18. Abdelkawy, M.A.; Amin, A.Z.M.; Lopes, A.M.; Hashim, I.; Babatin, M.M. Shifted Fractional-Order Jacobi Collocation Method for Solving Variable-Order Fractional Integro-Differential. Fractal Fract. 2021, 6, 19. [Google Scholar] [CrossRef]
  19. Khalil, H. The use of Jacobi polynomials in the numerical solution of coupled system of fractional differential Equation. Int. J. Comput. Math. 2014, 92, 1452–1472. [Google Scholar] [CrossRef]
  20. Abdou, M.A. On the solution of linear and nonlinear integral Equation. J. App. Math. Comp. 2014, 146, 857–871. [Google Scholar] [CrossRef]
  21. Abdou, M.A.; Raad, S.A. Nonlocal solution of a nonlinear partial differential Equation and its equivalent of nonlinear integral Equation. Inter. J. Comp. Math. 2014, 92, 1452–1472. [Google Scholar] [CrossRef]
  22. Mosal, G.A.; Abdou, M.A.; Rahby, A.S. Numerical solutions for nonlinear Volterra-Fredholm integral equations of the second kind with a phase lag. AIMS Math 2021, 6, 8525–8543. [Google Scholar] [CrossRef]
Figure 1. Approximate solution   φ x , t = ( t + 0.5 ) x 2 where   T = 0.9 .
Figure 1. Approximate solution   φ x , t = ( t + 0.5 ) x 2 where   T = 0.9 .
Fractalfract 07 00656 g001
Figure 2. The error of example 1, case (i), where N = 4 ,     T = 0.9 .
Figure 2. The error of example 1, case (i), where N = 4 ,     T = 0.9 .
Fractalfract 07 00656 g002
Figure 3. The error of example 1, case (i), where N = 16 ,   T = 0.1 .
Figure 3. The error of example 1, case (i), where N = 16 ,   T = 0.1 .
Fractalfract 07 00656 g003
Figure 4. The error of example 1, case (i), where N = 8 ,     T = 0.5 .
Figure 4. The error of example 1, case (i), where N = 8 ,     T = 0.5 .
Fractalfract 07 00656 g004
Figure 5. The error of example 1, case (i), where N = 16 ,   T = 0.5 .
Figure 5. The error of example 1, case (i), where N = 16 ,   T = 0.5 .
Fractalfract 07 00656 g005
Figure 6. The error of example 1, case (i), where N = 16 ,   T = 0.9 .
Figure 6. The error of example 1, case (i), where N = 16 ,   T = 0.9 .
Fractalfract 07 00656 g006
Figure 7. The error of example 1, case (ii), where N = 16 ,     T = 0.27 .
Figure 7. The error of example 1, case (ii), where N = 16 ,     T = 0.27 .
Fractalfract 07 00656 g007
Figure 8. The error of example 1, case (ii), where N = 16 ,     T = 0.75 .
Figure 8. The error of example 1, case (ii), where N = 16 ,     T = 0.75 .
Fractalfract 07 00656 g008
Figure 9. Approximate solution   φ x , t = 0.3 cos t X .
Figure 9. Approximate solution   φ x , t = 0.3 cos t X .
Fractalfract 07 00656 g009
Figure 10. The error of example 2, where T = 0.9 .
Figure 10. The error of example 2, where T = 0.9 .
Fractalfract 07 00656 g010
Figure 11. The error of example 2, where T = 0.03 .
Figure 11. The error of example 2, where T = 0.03 .
Fractalfract 07 00656 g011
Table 1. The error for different.   ν . of the approximate solution in case (i) in which X ( x ) = [ 0.5625 ,   0.2500 ,   0.0625 ,   0.2500 ,   0.5625 ] at the given points above reduced to 10 5 .
Table 1. The error for different.   ν . of the approximate solution in case (i) in which X ( x ) = [ 0.5625 ,   0.2500 ,   0.0625 ,   0.2500 ,   0.5625 ] at the given points above reduced to 10 5 .
TX Error   ν = 0.01 Error   ν = 0.47
0.0–0.752.000 × 10−102.000 × 10−10
−0.51.000 × 10−101.000 × 10−10
0.251.000 × 10−111.000 × 10−11
0.50.0000.000
0.750.0000.000
0.1−0.750.0000.000
−0.51.000 × 10−101.000 × 10−10
0.251.000 × 10−111.000 × 10−11
0.51.000 × 10−101.000 × 10−10
0.750.0000.000
0.5−0.758.000 × 10−101.200 × 10−9
−0.59.000 × 10−101.300 × 10−9
0.258.400 × 10−101.400 × 10−10
0.59.000 × 10−101.200 × 10−9
0.759.000 × 10−101.200 × 10−9
0.9−0.758.770 × 10−81.294 × 10−7
−0.58.460 × 10−81.326 × 10−7
0.258.533 × 10−81.337 × 10−8
0.58.480 × 10−81.327 × 10−7
0.758.750 × 10−81.293 × 10−7
Table 2. The error for different.   ν . of the approximate solution in case (ii).
Table 2. The error for different.   ν . of the approximate solution in case (ii).
TX Error   ν   =   0.01 Error   ν   =   0.47
0.0–0.750.0000.000
−0.51.000 × 10−101.000 × 10−10
0.250.0000.000
0.50.0000.000
0.751.000 × 10−101.000 × 10−10
0.1−0.752.000 × 10−102.000 × 10−10
−0.51.000 × 10−101.000 × 10−10
0.251.000 × 10−111.000 × 10−11
0.50.0000.000
0.752.000 × 10−102.000 × 10–10
0.5−0.759.000 × 10−101.300 × 10−9
−0.59.000 × 10−101.400 × 10−9
0.258.400 × 10−101.300 × 10−10
0.51.000 × 10−91.200 × 10−9
0.759.000 × 10−101.300 × 10−9
0.9−0.758.770 × 10−81.293 × 10−7
−0.58.460 × 10−81.326 × 10−7
0.258.532 × 10−81.337 × 10−8
0.58.480 × 10−81.326 × 10−7
0.758.740 × 10−81.293 × 10−7
Table 3. Convergence rate in both methods with fixed time T = 0.0 .
Table 3. Convergence rate in both methods with fixed time T = 0.0 .
NMean ErrorRate
42.4 × 10−100.49
81.7 × 10−100.79
169.8 × 10−110.16
328.8 × 10−11---
Table 4. Convergence rate in both methods with fixed time T = 0.3 .
Table 4. Convergence rate in both methods with fixed time T = 0.3 .
NMean ErrorRate
44.8 × 10−100.94
82.5 × 10−101.32
161.0 × 10−100.12
329.2 × 10−11---
Table 5. Convergence rate in both methods with fixed time T = 0.6 .
Table 5. Convergence rate in both methods with fixed time T = 0.6 .
NMean ErrorRate
45.8 × 10−80.16
86.5 × 10−83.92
164.3 × 10−92.10
321.0 × 10−9---
Table 6. Convergence rate in both methods with fixed time T = 0.9 .
Table 6. Convergence rate in both methods with fixed time T = 0.9 .
NMean ErrorRate
48.9 × 10−72.89
81.2 × 10−73.79
168.7 × 10−91.92
322.3 × 10−9---
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jan, A.R.; Abdou, M.A.; Basseem, M. A Physical Phenomenon for the Fractional Nonlinear Mixed Integro-Differential Equation Using a Quadrature Nystrom Method. Fractal Fract. 2023, 7, 656. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7090656

AMA Style

Jan AR, Abdou MA, Basseem M. A Physical Phenomenon for the Fractional Nonlinear Mixed Integro-Differential Equation Using a Quadrature Nystrom Method. Fractal and Fractional. 2023; 7(9):656. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7090656

Chicago/Turabian Style

Jan, A. R., M. A. Abdou, and M. Basseem. 2023. "A Physical Phenomenon for the Fractional Nonlinear Mixed Integro-Differential Equation Using a Quadrature Nystrom Method" Fractal and Fractional 7, no. 9: 656. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7090656

Article Metrics

Back to TopTop