Next Article in Journal
Weak Solution for a Fractional Langevin Inclusion with the Katugampola–Caputo Fractional Derivative
Next Article in Special Issue
An Application of the Distributed-Order Time- and Space-Fractional Diffusion-Wave Equation for Studying Anomalous Transport in Comb Structures
Previous Article in Journal
An Advanced Fractional Order Method for Temperature Control
Previous Article in Special Issue
Fourth-Order Numerical Solutions for a Fuzzy Time-Fractional Convection–Diffusion Equation under Caputo Generalized Hukuhara Derivative
 
 
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 General Discontinuous Kernel

by
Sharifah E. Alhazmi
1,* and
Mohamed A. Abdou
2
1
Mathematics Department, Al-Qunfudah University College, Umm Al-Qura University, Al-Qunfudhah 28821, Saudi Arabia
2
Department of Mathematics, Faculty of Education, Alexandria University, Alexandria P.O. Box 21526, Egypt
*
Author to whom correspondence should be addressed.
Submission received: 7 December 2022 / Revised: 30 January 2023 / Accepted: 31 January 2023 / Published: 9 February 2023

Abstract

:
In this study, a fractional nonlinear mixed integro-differential equation (Fr-NMIDE) is presented and has a general discontinuous kernel based on position and time space. Conditions of the existence and uniqueness of the solution is provided through the principal form of the integral equation, based on the Banach fixed point theorem. After applying the properties of a fractional integral, the Fr-NMIDE conformed to the Volterra–Hammerstein integral equation (V-HIE) of the second kind, with a general discontinuous kernel in position with the Hammerstein integral term and a continuous kernel in time to the Volterra term. Then, using a technique of the separating method, we obtained HIE, where its physical coefficients were variable in time. The Toeplitz matrix method (TMM) and its schemes were used to obtain a nonlinear algebraic system by studying the convergence of the system. The Maple 18 program was implemented to present the numerical results, along with corresponding errors.

1. Introduction

As integral equations, integro-differential equations and fractional integro-differential equations (IEs/IDEs/fIDEs) can be used to simulate a wide range of problems in the basic sciences, many scientists have focused their attention on presenting the solutions for these systems. These equations have played a significant role in finding solutions using diverse methods, which is in line with the rapid development in finding the answers to diverse problems originating from the basic sciences. Currently, several studies have concentrated on creating more sophisticated and effective techniques for solving the IEs/IDEs, such as the Rieman–Stieltjes integral conditions [1,2], Lerch polynomials method [3] and Legendre–Chebyshev spectral method [4], along with the numerical observations based on semi-analytical approaches, e.g., Adomian’s decomposition method [5] and HOBW method [6]. The linear/nonlinear equations (IEs/IDEs/ fIDEs) have various uses in fluid mechanics [7], Stokes flow [8], airfoil [9], quantum mechanics [10], integral models [11], mathematical engineering [12], nuclear physics [13] and the theory of laser [14]. The orthogonal polynomials method is considered one of the most significant operators used to solve various scientific problems. Alhazmi [15] used a new technique based on the separation of variables and the orthogonal polynomials method to obtain many spectral relationships based on the mixed integral equation using a generalized potential kernel. Nemati et al. [16] applied the Legendre polynomials scheme for the outcomes of a second-order two-dimensional (2D) Volterra integral model, together with a continuous kernel. Mirzaee and Samadyar [17] discussed the convergence of 2D-orthonormal Bernstein collocation method for solving 2D-mixed Volterra-Fredholm integral equations. Basseem and Alalyani [18] used Chebyshev polynomials to discuss the numerical solution of the quadratic integral equation with a logarithmic kernel. Katani [19] implemented a quadrature scheme for the numerical outcomes of the second kind of Fredholm integral model. Al-Bugami [20] used the Simpson and Trapezoidal schemes to perform numerical representations based on an integral model using 2D surface crack layers. Brezinski and Zalglia [21] used the extrapolation approach to achieve numerical computing results based on the second kind of nonlinear integral model that has a continuous kernel. Baksheesh [22] proposed using the Galerkin scheme to find the approximate results based on the Volterra integral equations of the second kind, which have discontinuous kernels. Alkan and Hatipoglu [23] applied the sinc-collocation method for solving the Volterra-Fredholm IDEs of fractional order. Mosa et al. [24] studied the semi-group scheme to assess uniqueness and existence based on the partial and fractional integro models of heat performance in the Banach space using the Adomian decomposition scheme. Bin Jebreen and Dassios [25] proposed an efficient algorithm to find an approximate solution via the wavelet collocation method for fractional Fredholm integro-differential equations. Akram et al. [26] interpreted the collocation approach to tackle the fractional partial integro-differential equation by employing the extended cubic B-spline Abdelkawy et al. [27] applied the Jacobi–Gauss collocation method after using the Riemann–Liouville fractional integral and derivative fractional to obtain the approximate solution for variable-order fractional integro-differential equations with a weakly singular kernel.
The initial value of the Fr-NMIDE is presented as:
μ α Φ ( x , t ) t α + v Φ ( x , t ) = λ Ω k ( | x y | ) Φ m ( y , t ) d y + g ( x , t ) , ( Φ ( x , 0 ) = ψ ( x ) ) .
Here, g ( x , t ) and Φ ( x , t ) are the known and unknown continuous functions, respectively, in L 2 [ Ω ] X C [ 0 , T ] . Ω is the integration domain and m = 1 , 2 , M . In addition, μ indicates the constant values of Equation (1) and λ and v are constants and have physical sense. The kernel k ( | x y | ) , in general, has a singular term. For reference, the essential properties and definitions have been stated using fractional calculus theory.
Definition 1.
For the function f : ( 0 , ) R , the fractional Riemann–Liouville integral with order α > 0 is shown as [28]:
I 0 + α l ( t ) = 1 Γ ( α ) 0 t ( t s ) α 1 l ( s ) d s .
In addition, we define the Caputo derivatives of order as α > 0
D 0 + α l ( t ) = 1 Γ ( n α ) 0 t l ( n ) ( s ) d s ( t s ) α n + 1 , ( n < α n + 1 ) .
The time Abel kernel ( t τ ) α 1 , t , τ [ 0 , T ] , 0 τ t T < 1 satisfies the following features: 0 t 1 t 2 t T < 1 for every continuous function h ( t ) . Integrals
0 t ( t τ ) α 1 h ( τ ) d τ , max 0 t T 0 t ( t τ ) α 1 d τ , t 1 t 2 ( t τ ) α 1 h ( τ ) d τ ,
are the continuous time function, i.e.,
0 t ( t τ ) α 1 h ( τ ) d τ M .
In this paper, Section 2 presents the fractional definition to obtain a NMIE. Then, the theorem-based Banach fixed point is discussed to prove the existence and uniqueness of the solution of NMIE. Section 3 presents the convergence of the solution. Section 4 indicates the technique of separation of the variables for the Hammerstein integral model in position and its coefficients. This scheme would help researchers choose the time known function in an easier way, enabling them to choose the necessary time to obtain the required results. Section 5 indicates the convergence analyses of the Hammerstein integral equation. Section 6 shows the Toeplitz matrix scheme to conform the Hammerstein integral model of the nonlinear algebraic system (NAS). The TMM is considered the best numerical method for solving singular integral equations, where the singular terms disappear and we have simple integrals. Section 7 represents the NAS convergence. Section 8 provides the convergence of the error using one of the famous theorems. Section 9 provides the numerical solutions through Maple 18, together with the kernel of the nonlinear integral equation that takes the logarithmic form, the Carleman function and Hilbert kernel. In addition, the corresponding errors are computed.

2. The Solution’s Existence and Uniqueness

The fundamental Caputo fractional integral is used to find the second order NMIE as:
μ Φ ( x , t ) + v Γ ( α ) 0 t ( t τ ) α 1 Φ ( x , τ ) d τ λ Γ ( α ) 0 t Ω ( t τ ) α 1 k ( | x y | ) Φ m ( y , τ ) d y d τ = f ( x , t ) , f ( x , t ) = 1 Γ ( α ) 0 t ( t τ ) α 1 g ( x , τ ) d τ + ψ ( x ) , 0 < α < 1 .
In Equation (2), the free term f ( x , t ) L 2 ( Ω ) × C [ 0 , T ] and the unknown function Φ ( x , t ) will be discussed in the same space, L 2 ( Ω ) × C [ 0 , T ] , along with the discontinuous kernel k | x y | L 2 [ Ω ] × [ Ω ) . The discontinuous kernel of time ( t s ) α 1 , 0 < α < 1 , is considered in class C [ 0 , T ] , where T < 1 .
To prove the existence of NMIE (2) based on the unique results, the integral operator form is given below:
U ¯ Φ ( x , t ) = λ μ U 2 Φ ( x , t ) v μ U 1 Φ ( x , t ) + 1 μ f ( x , t ) , U 1 Φ ( x , t ) = 1 Γ ( α ) 0 t ( t τ ) α 1 Φ ( x , y , τ ) d τ , U 2 Φ ( x , t ) = 1 Γ ( α ) 0 t Ω ( t τ ) 1 + α k ( | y x | ) Φ ( y , τ ) d y d τ .
The following conditions are presented as:
(i) (i-a) The position kernel k ( | x y | ) satisfies
Ω Ω k ( | x y | ) 2 d x d y 1 2 d τ = C , ( C constant ) .
(i-b) Therefore, the kernel of position and time ( t τ ) α 1 k ( | x y | ) in L 2 ( Ω ) × C [ 0 , T ] satisfies
max 0 t T 0 t Ω Ω ( t τ ) α 1 k ( | x y | ) 2 d x d y 1 2 d τ = T α C α ( 0 < α < 1 ) .
(ii) The continuous function f ( x , t ) L 2 ( Ω ) × C [ 0 , T ] and its norm is shown as f ( x , t ) L 2 ( Ω ) × C [ 0 , T ] = max 0 t T 0 t Ω | f ( x , τ ) | 2 d x 1 2 d τ = G and G is taken as a constant.
(iii) The decreasing function Q > P , constant Q > Q 1 and Φ m ( x , t ) are used as:
(iii-a) max 0 t T 0 t Ω Φ m ( x , τ ) 2 d x 1 2 d τ Q 1 Φ ( x , t ) L 2 ( Ω ) × C [ 0 , T ] .
(iii-b) Φ 1 m ( y , t ) Φ 2 m ( y , t ) N ( t , x ) Φ 1 ( x , t ) Φ 2 ( x , t ) , | N ( t , x ) | = p .
Theorem 1
(Banach Fixed Point [29]). Consider X = ( X , d ) to be a metric space, where X { Φ } , Φ is a null set. Suppose X is complete and T : X X is the X contraction, then T contains exactly a single fixed point.
To prove this, K shows the contraction operator on B with integral form K Ψ = Ψ , Ψ presents the unique form of the solution, the general kernel H | y x | , ( t τ ) = ( t τ ) α 1 × k ( | y x | ) satisfies in (i-b), L 2 ( Ω ) × C [ 0 , T ] is the domain of integration, with respect to position Ω and the time t , τ [ 0 , T ] , 0 τ t T < 1 .
Theorem 2.
Principal theorem: The NMIE (2) with the use of above conditions and the space L 2 ( Ω ) × C [ 0 , T ] takes the form of:
| μ | Γ ( α + 1 ) > ( v + C Q λ ) T α .
Lemmas (1) and (2) must be proven to satisfy the above theorem.
Lemma 1.
Under conditions (i) to (iii-a), the W ¯ operator maps the space L 2 ( Ω ) × C [ 0 , T ] onto itself:
Proof. 
Equation (3) is used to prove
U ¯ Φ ( x , t ) λ | μ | U 2 Φ ( x , t ) + 1 | μ | f ( x , t ) + | v | | μ | U 1 Φ ( x , t ) .
Using (i)–(iii-a) and the inequality of Cauchy–Schwarz, we have:
U ¯ Φ ( x , t ) Q | μ | + σ Φ ( x , t ) , σ = ( ν + λ C Q ) T α | μ | Γ ( α + 1 ) .
In the above statement, the U ¯ operator maps the S r ball as:
r = G r ( α + 1 ) | μ | Γ ( α + 1 ) ( v + λ C Q ) T α .
As r > 0 , G Γ ( α + 1 ) > 0 , therefore σ < 1 . Furthermore, the inequality (5) includes the boundedness operators U 1 , U 2 and U ¯ . □
Lemma 2.
If (i)–(iii-b) conditions are fulfilled, then U ¯ operator is a contractive-based Banach space L 2 ( Ω ) × C [ 0 , T ] .
Proof. 
For Φ 1 ( x , t ) and Φ 2 ( x , t ) functions using L 2 ( Ω ) × C [ 0 , T ] space, Formula (3) becomes:
U ¯ Φ 1 ( x , t ) Φ 2 ( x , t ) | v | | μ | U 1 Φ 1 ( x , t ) Φ 2 ( x , t ) + U 2 Φ 1 ( x , t ) Φ 2 ( x , t ) ,
The conditions (i), (ii) and (iii-b) have been applied to the Cauchy–Schwarz inequality as:
U ¯ Φ 1 ( x , t ) Φ 2 ( x , t ) σ Φ 1 ( x , t ) Φ 2 ( x , t ) .
Inequality (7) presents the operator U ¯ (contraction operator), which shows the continuity in the L 2 ( Ω ) × C [ 0 , T ] space. □

3. Convergence of the Solution

Consider the simple iteration Φ 1 ( x , y ) , , Φ n 1 ( x , t ) , Φ n ( x , t ) , Φ ( x , t ) , where the two functions Φ n 1 ( x , t ) , Φ n ( x , t ) are used to satisfy
μ Φ n ( x , t ) Φ n 1 ( x , t ) + v Γ ( α ) 0 t ( t τ ) α 1 Φ n 1 ( x , t ) Φ n 2 ( x , t ) d τ = λ Γ ( α ) 0 t Ω ( t τ ) α 1 k ( | x y | ) Φ n 1 m ( y , τ ) Φ n 2 m ( y , τ ) d y d τ
Consider
Φ n ( x ; t ) = i = 0 n Ψ i ( x ; t ) ,
where
Ψ n ( x , t ) = Φ n ( x , t ) Φ n 1 ( x , t ) ; ( n 1 ) , Ψ 0 ( x , t ) = f ( x , t ) .
Equation (8) is updated by using Equation (9):
μ Ψ n ( x , t ) v Γ ( α ) 0 t ( t τ ) α 1 Ψ n 1 ( x , t ) d τ + λ Γ ( α ) 0 t Ω ( t τ ) α 1 k ( | x y | ) Ψ n 1 m ( y , τ ) d y d τ .
Taking n = 1 , the above formula becomes:
Ψ 1 ( x , t ) σ G , σ = ( v + λ C Q ) T α | μ | Γ ( α + 1 ) ,
and we see
Ψ n ( x , t ) σ n G , σ < 1 .
Equation (10) shows the convergent sequence Ψ n ( x , t ) uniformly. Moreover, it provides the convergent solution of the sequence Φ n ( x , t ) . As Ψ i ( x ; t ) is continuous and Φ ( x , t ) = lim n Φ n ( x , t ) = lim n i = 0 n Ψ i ( x ; t ) , Φ ( x ; t ) is uniformly continuous with an infinite Φ n ( x , t ) n = 0 series. This proves the lemma.

4. Separation of Variables Scheme

In the problems of mathematical physics, we find that researchers are interested in finding the unidentified potential function, which is linked to time and position. A variety of methods can be used to obtain the unknown function. One of these methods is time division that turns the mixed integral equation into an algebraic system of integral equations. Researchers apply the separating variable method to solve the mixed integral equation using the coefficients of the space functions. Moreover, these time coefficients are in the form of an integral operator of the Volterra type (Jan [30,31]). This scheme helps researchers to choose the time known function in an easier way, which enables them to choose the necessary time to obtain the required results. The unidentified and known functions Φ ( x , t ) and g ( x , t ) are shown in the separation form as:
Φ ( x , t ) = X ( x ) Y ( t ) , g ( x , t ) = b ( x ) Y ( t ) , Y ( 0 ) 0 ,
where X ( x ) is an unknown function in a position that is to be determined, b ( x ) is the given function in a position and Y ( t ) shows the known function in time.
The time function is chosen in the form of a series based on the polynomial constants. This form helps the researcher to categorize the function of time based on the constants as a famous function or time representation in several other forms. It is noted that time series convergence is based on the premise of the experiment time being less than one and the start time is not equal to zero. Assume that
Y ( t ) = t α n = 0 a n t n , a 0 0 , t [ 0 , T ] , T < 1 .
Using Equations (11) and (12) in Equation (2), shown as
μ X ( x ) ρ ( t ) Ω k ( | y x | ) X m ( y ) d y = H ( x , t ) ,
where
ρ ( t ) = λ Q ( t ) 0 t ( t τ ) α 1 Y m ( τ ) d τ , m = 1 , 2 , , M , Q ( t ) = v 0 t ( t τ ) α 1 Y ( τ ) d τ + μ Γ ( α ) Y ( t ) 1 , H ( x , t ) = Q ( t ) Γ ( α ) ψ ( x ) + h ( x ) 0 t ( t τ ) α 1 Y ( τ ) d τ .
By using the updated form of Equation (13) by using Equation (14), the reader can use the separation of time scheme for the nonlinear mixed integral model in time and position, which leads to a nonlinear integral equation in a position with coefficients linked to t, where t [ 0 , T ] , T < 1 . Moreover, Equation (13) represents that the single solution condition on the nonlinear integral equation is an equivalence relationship between position and time, which is given as:
k ( | x y | ) δ , δ = μ Γ ( α ) Y ( t ) + v 0 t ( t τ ) α 1 Y ( τ ) d τ λ 0 t ( t τ ) α 1 Y m ( τ ) d τ 1 < 1 .
Equation (15) shows the relationship between the position represented in the kernel form (which represents the properties of matter) and the time required for the continuity of these properties, which is known under the condition that there is a single solution. It is noted that at a certain time, an increase that exceeds the standard value of the nucleus may occur, and this leads to instability.

5. Convergence Investigations Based on Nonlinear Integral Model

To check the nonlinear integral model (13) convergence, the solution sequence takes the form:
X ( x ) = X 0 ( x ) , X 1 ( x ) , , X m ( x ) , X ( x ) , ,
where ψ n and ψ m are two distinct arbitrary partial sums of sequence X j ( x ) : ( n > m ) , and
d ψ n , ψ = max n , m j X n X δ d X n , X 1 δ 2 d X n , X 2 δ d X n , X 0 ,
where δ is defined in the Equation (15). Finally, we have
d ψ n , ψ δ 1 δ d ψ 1 , ψ 0 ; 0 < δ < 1 .
As and for the fixed values, d ψ 1 , ψ 0 must approach zero. Hence, it is concluded that ψ n shows a Cauchy sequence throughout the metric space, which shows the convergent series.

6. Toeplitz Matrix Method (Abdou et al. [32])

Many numerical methods have been used to solve integral equations with continuous or unconnected kernels. In singular integral equations, the best way to solve them is the Toeplitz matrix method (TMM), due to the following reasons: The singular term directly disappears, being transformed into simple integrals that can be solved quickly, and then, forms a linear/nonlinear algebraic system of equations. The degree of convergence in the relative error is less than in the other methods.
There are some methods that give a quick approximation when studying the error, including a circulated preconditioned iterative scheme. Xian et al. [33] considered the conjugate gradient technique with a circulated preconditioned scheme for a linear discretized model. They obtained a matrix whose coefficient kept the structure of a symmetric Toeplitz-plus-tridiagonal. They also studied linear large-scale systems with coefficient matrices based on the non-Hermition Toeplitz and applied a product of a fast Toeplitz matrix-vector through iterative solvers for a linear discretized model (circulant preconditioned method). This numerical scheme provides fast results to reduce computational costs, using O m 3 to O ( m log m ) , and storage through O m 2 to O ( m ) deprived of loss compression, where m is the number of spatial grid nodes.
For TMM, consider Ω ( 1 , 1 ) , using the nonlinear integral term of ( 2 ) as:
1 1 k ( | y x | ) X m ( y ) d y = = N N 1 a = h ( + 1 ) h k ( | y x | ) X m ( y ) d y = = N N 1 A ( x ) X m ( h ) + B ( x ) X m ( h + h ) + R , h = 1 N .
The two functions A ( x ) and B ( x ) are
A ( x ) = 1 Δ ( h + h ) m I ( x ) J ( x ) , B ( x ) = 1 Δ J ( x ) ( h ) m I ( x ) ,
where
Δ = h m ζ = 1 m Γ ( m + 1 ) Γ ( ζ + 1 ) Γ ( m ζ + 1 ) ( m ζ ) , ( N N )
and
I ( x ) = h ( + 1 ) h k ( | y x | ) d y , J ( x ) = h ( + 1 ) h k ( | y x | ) y m d y .
The integral term of Equation (17) after using Equation (18) becomes:
1 1 k ( | x y | ) X m ( y ) d y = = N N D ( x ) X m ( h ) D ( x ) = A N ( x ) , = N A n ( x ) + B n 1 ( x ) , N < < N B N 1 ( x ) , = N .
Equations (18) and (20) are used in x = j h , ( N j N ) , to obtain
X ( j h ) = X j , A ( j h ) = A , j , B ( j h ) = B , j , D ( j h ) = D , j , H ( j h , t ) = H j ( t ) .
The nonlinear integral in Equation (15) shows the NAS of a (2N + 1) system
μ X j ρ ( t ) = N N D , j X m = H j ( t ) , N , j N .
The matrices D n , show the Toeplitz matrix as:
D , j = V n U n , V j = A , j + B 1 , j , ( N , N j ) , U l , j = B N 1 , j = N 0 N < < N A N , j = N .
Equation (23) shows the matrix of two types with a ( 2 N + 1 ) order, including V j (Toeplitz matrix) and U , j , which shows the zero elements, excluding the first and last columns or rows. The error R n , can be obtained based on the following formula:
R = max 0 j N h h + h y 2 m k ( | x y | ) d y G ( x ) ( h ) 2 m + H ( x ) ( h + h ) 2 m = O h 3 m , ( x = j h ) .

7. The Nonlinear Algebraic Toeplitz Matrix System

This section provides the existence based on NAS (22) using the Banach space × C [ 0 , T ] . The operator form is written as:
T ¯ X j = T X j + H j μ ,
where
T X j = ρ ( t ) μ = N N D , j X m ; ( M m M , 0 t T < 1 ) .
The following lemma is
Lemma 3.
If the position kernel satisfies the conditions below:
( i ) h h + h j h j h + h k 2 ( | x y | ) d x d y 1 2 C , ( i i ) lim x x k x , y k ( x , y ) L 2 = 0 , x , x ( 1 , 1 ) .
Then,
( a ) sup N j = N N D , j e x i s t s , ( b ) lim l N sup N j = N N D , j D , j = 0 .
Proof. 
From Equations (18) and (19), we have
A ( x ) 1 | Δ | ( h + h ) m h ( + 1 ) h k ( | y x | ) d y + h ( + 1 ) h k ( | y x | ) y m d y
Applying the Cauchy–Schwarz inequality and taking the sum from = N to = N , the above inequality yields
= N N A ( x ) 1 | Δ | k ( | x y | ) = N N ( h + h ) m + y m
Based on Equation (27), the function continuity y m in input ( 1 , 1 ) , a small constant E 1 exists, i.e., = N N A ( x ) E 1 , N . As each value of = N N A ( x ) is bounded, take x = j h as:
sup N = N N A ( j h ) E 1 .
Similarly, by taking a small value of the constant E 2 for Equations (18) and (19), we have
sup N n = N N B n ( m h ) E 2 .
Using (29) and (30), we have
sup N = N N D , j sup N = N N A ( j h ) + sup N = N N B ( j h ) E .
Hence, sup N = N N D , j exists.
To prove the second equation of (28) for x , x ( 1 , 1 ) , the Cauchy–Schwarz inequality is applied by taking the sum from = N to
sup N = N N A x A ( x ) 1 | Δ | k x , y k ( x , y ) L 2 sup N = N ( h + h ) m + y m
Using x = j h , x = j h and Equation (27), we obtain that as x x ,
lim j j sup N = N N A j h A ( j h ) = 0 .
Similarly, from (18) and (19), it is proved that
lim j m sup N t = N N B j h B ( j h ) = 0 .
Finally, we obtain
lim m m sup N n = N N D m n D m n = 0 .
Now, the principal theorem is proven based on the nonlinear algebraic systems.
Theorem 3.
The NAS (22) using the Banach × C [ 0 , T ] space shows the unique form of the solution, as follows:
sup j H j ( t ) H ¯ < ,
sup N = N N D , J E ¯ ,
where H ¯ and E ¯ are constants.
The functions X m ( j h ) , where m = 1 , 2 , , M for the constants Q ¯ > Q ¯ 1 , Q ¯ > P ¯ 1 satisfy the following:
sup j X m ( j h ) Q 1 ¯ X
sup j X m ( j h ) Z m ( j h ) P ¯ 1 X Z
where X = sup j X j , X ( j h ) = X j for each integer j.
The below lemmas must be proven for the above theorem.
Lemma 4.
If Equations (34)–(36) are tested, then the T ¯ operator is defined by using Equation (25), which maps the × C [ 0 , T ] space onto itself.
Proof. 
Suppose U shows all the functions set as X = X j in × C [ 0 , T ] , i.e., Φ × C [ 0 , T ] β ¯ , β ¯ is constant and the T ¯ Φ operator is based on the Banach × C [ 0 , T ] space:
T ¯ X × C [ 0 , T ] = sup j T ¯ X j , for all j
Using conditions (34) and (35), we have
T ¯ X j ρ ( t ) μ Q X sup j = N N D , j + sup j H ¯ μ .
For each integer j, the above inequality is shown as:
sup j T ¯ X j σ 1 X + H ¯ μ , σ 1 = ρ μ Q ¯ E ¯
The above inequality (38) represents that the T ¯ operator maps the U set, where
β ¯ = H ¯ ( | μ | | ρ | Q ¯ E ¯ ) .
Hence, σ 1 < 1 have been taken, whereas T and T ¯ operators are bounded. □
Lemma 5.
Under the conditions (34), (35) and (37), T ¯ shows a contraction operator in the × C [ 0 , T ] space
Proof. 
For X and Z functions in × C [ 0 , T ] , Formulas (25) and (26) become:
T ¯ X j T ¯ Z j ρ μ n = N N D , j sup j X j Z j .
Using conditions (35) and (37),
T ¯ X T ¯ Z × C [ 0 , T ] σ 1 X Z × C [ 0 , T ] , σ 1 = ρ μ Q ¯ E ¯ .
The above form indicates that the T ¯ operator is continuous in × C [ 0 , T ] space. T ¯ shows the contraction operator based on σ 1 < 1 . Hence, T ¯ presents a unique fixed point with specific solutions of the NAS using × C [ 0 , T ] space. □

8. The Error of the Toeplitz Matrix Method

In any practical use of the TMM, some estimation of the error size is involved. Hence, these two definitions were used to calculate the error of the TMM.
Definition 2.
A local error R j is used as:
X ( x ) X ξ ( x ) = = N N D , j X j m X j , ξ m + R j , ( x = j h ) ,
where X ξ ( x ) shows the approximate results of (2).
In the other formula, Equation (40) is used as:
R j = 1 1 k ( | y x | ) X m ( y ) d y = N N D , j X j m .
Definition 3.
The TMM shows a convergence of r order in input ( 1 , 1 ) ; conversely, by taking the large values of N , D ¯ > 0 exists based on N independently, that is,
X ( x ) X N ( x ) D ¯ N r .
Now, we present the theorem, which is provided using NAS, based on Equation (22), which has a unique solution.
Theorem 4.
The error R j is considered negligible as j
lim j R j = 0 .
Proof. 
Equation (40) is used as:
R ξ X j X j ξ + = N N D , j sup j X j m X j m ξ .
Using conditions (35) and (37), along with each integer ξ , we obtain
R ξ × C [ 0 , T ] ( 1 + E Q ) X X ξ × C [ 0 , T ] .
Since X X ξ × C [ 0 , T ] 0 as  ξ , t [ 0 , T ] , then R j 0 . □

9. Applications

For numerical results and tables, we used Maple 2022.1 software, Version 15, March 2022, Windows, 10, 8 G RAM, 64-bit. In this section, the NMIE of (1) in the following special form was considered
α Φ ( x , t ) t α + 0.5 Φ ( x , t ) = g ( x , t ) + 0.33 1 1 k ( | y x | ) Φ m ( y , t ) d y , Φ ( x , 0 ) = x 2
The true result of Equation (44) is Φ ( x , t ) = 0.5 t 0.5 + 0.25 t 1.5 + x 2 .
Example 1
(For logarithmic kernel).
k ( | x y | ) = ln ( | x y | ) .
After using (17)–(19) including the famous integral reference as Gradstein et al. [34]:
x m ln ( a + b x ) d x = 1 m + 1 x m + 1 ( a ) m + 1 b m + 1 ln ( a + b x ) + 1 m + 1 k = 1 m + 1 ( 1 ) k x m k + 2 a k 1 ( m k + 2 ) b k 1 ,
we have
A J ( h ) = h J m ( J + 1 ) m 1 ( m + 1 ) ( J + 1 ) m + 1 m + 1 ln | ( J 1 ) h | J m + 1 m + 1 ln | ( J ) h | + ( J + 1 ) m [ ( J 1 ) ln | ( J 1 ) h | ( J ) ln | ( J ) h | + 1 ] 1 ( m + 1 ) k = 1 m + 1 ( J + 1 ) m k + 2 J m k + 2 k 1 ( m k + 2 )
and
B J ( h ) = h ( J + 1 ) m J m 1 ( m + 1 ) ( J + 1 ) m + 1 m + 1 ln | ( J 1 ) h | J m + 1 m + 1 ln | ( J ) h | + h ( J + 1 ) m J m J m [ ( J 1 ) ln | ( J 1 ) h | ( J ) ln | ( J ) h | + 1 ] 1 ( m + 1 ) k = 1 m + 1 ( J + 1 ) m k + 2 J m k + 2 k 1 ( m k + 2 ) .
Using (47) and (48), the coefficients D , J of the nonlinear algebraic system (22) are presented as:
D , J = h J m ( 1 + J ) m ( 1 + m ) 1 ( 1 + J ) m + 1 m + 1 ln | ( 1 J ) h | J m + 1 m + 1 ln | ( J ) h | + ( J + 1 ) m [ ( J 1 ) ln | ( J 1 ) h | ( J ) ln | ( J 1 ) h | + 1 ] 1 ( m + 1 ) k = 1 m + 1 ( J + 1 ) m k + 2 J m k + 2 k 1 ( m k + 2 ) + h J m ( J 1 ) m 1 ( 1 + m ) J m + 1 m + 1 ln | ( J ) h | ( J 1 ) m + 1 m + 1 ln | ( J + 1 ) h | .
In addition,
R = a a + h ln x y φ i ( y ) d y A n ( x ) φ i ( a ) B n ( x ) φ i ( a + h ) .
The formula (50) takes the form:
R = C h 2 m + 1 , C = 2 m + 1 2 m + 1 m + 1 m + 1 ln h 2 m + 1 1 2 m + 1 m + 1 1 m + 1 ln h ( 1 ) k = 1 2 m + 1 k 1 ( 2 m k + 2 ) ( 2 m + 1 ) + k = 1 m + 1 k 1 ( m k + 2 ) ( m + 1 ) .
The linear case can be obtained from Equations (47)–(51) by letting m = 1 .
Example 2
(For Carleman kernel).
k ( | x y | ) = | x y | β ( 0 < β < 1 ) .
The significance based on the Carleman kernel was shown in the Arutiunian [35] work that presented the plane contact problem using nonlinear plasticity theory as its first calculation. It was reduced through the first-order Fredholm integral model with a Carleman kernel.
After using (17) to (19), we obtain
A J ( h ) = h 1 β J m ( 1 + J ) m k = 0 m m ! J m k | J | k + 1 β ( J + 1 ) m k | J 1 | k + 1 β ( m k ) ! ( 1 β ) ( 2 β ) ( k + 1 β ) + ( 1 + J ) m ( β + 1 ) [ J 1 1 m | J | 1 m
and
B J ( h ) = h 1 β J m + ( 1 + J ) m k = 0 m m ! J m k | J | k + 1 β ( J + 1 ) m k | J 1 | k + 1 β ( m k ) ! ( 1 β ) ( 2 β ) ( k + 1 β ) + J m ( 1 β ) [ J 1 1 β | J | 1 β
Therefore, the Toeplitz matrix D , J becomes:
D , J = h 1 β 1 J m ( J + 1 ) m k = 0 m m ! J m k | J | k + 1 β ( J + 1 ) m k | J 1 | k + 1 β ( m k ) ! ( 1 ) ( 2 β ) ( k + 1 β ) + ( J + 1 ) m ( 1 β ) | J 1 | 1 β | J | 1 β + 1 J m ( J 1 ) m k = 0 m m ! ( J 1 ) i k | J + 1 | k + 1 β J m k | J | k + 1 β ( m k ) ! ( 1 β ) ( 2 β ) ( k + 1 β ) + ( J 1 ) i ( 1 β ) | J | 1 β | J + 1 | 1 β
The error R is shown as:
| R | C h 2 m + 1 β , C = k = 0 m m ! k + 1 β 1 1 k + 1 β ( m k ) ! ( 1 β ) ( 2 β ) ( k + 1 β ) k = 0 2 m ( 2 m ) ! k + 1 β 1 1 k + 1 β ( 2 m k ) ! ( 1 β ) ( 2 β ) ( k + 1 β ) .
Example 3
(Suppose the Hilbert kernel).
k ( | x y | ) = cot x y 2 , φ ( ± π , t ) = 0 .
The exact output of Equation (57) is φ ( x , t ) = 0.5 t 0.5 + 0.25 t 1.5 sin x .
The integral equation based on the Hilbert kernel, together with the crack problem used in elasticity theory are discussed in [34].
n h ( n + 1 ) h x m cot x d x = s = 0 ( 1 ) s 2 s 2 s B 2 s ( m + 2 s ) ( 2 s ) ! x m + 2 s ( m 1 , | x | < π ) .
where Bernoulli numbers B 2 s is used as:
D J , = 2 ( J + 1 ) ln sin h ( J + 1 ) 2 2 ( J ) ln sin h ( J ) 2 + ( J 1 ) ln sin h ( J 1 ) 2 s = 0 ( 1 ) s 2 2 s B 2 s ( m + 2 s ) ( 2 s ) ! × ( J + 1 ) 1 + 2 s 2 ( J ) 1 + 2 s + ( J 1 ) 1 + 2 s , m = 1 , 2 , , M .
From the above, the results of using the logarithmic function and the effect of time in the first example, as well as the arithmetic error, are described numerically in Table 1 and Figure 1a,b in the nonlinear case ( m = 2 ) and linear case ( m = 1 ). In the second example, the results of using the Carleman function over time periods for the nonlinear and linear cases, as well as the resulting arithmetic error, were derived in Table 2 and Figure 2a,b, respectively. While Table 3 and Figure 3a,b represent the numerical results and arithmetic errors corresponding to the nonlinear and linear cases, respectively, for different values of the coefficient of the Carleman function. In the third example of the Hilbert kernel, the numerical error of the nonlinear case and the linear case has been calculated at t = 0.1 in Figure 4a,b, respectively. While for t = 0.4 the corresponding results of errors are described in Figure 5a,b.

10. Conclusions

The following conclusions were drawn:
1
In this paper, the existence of a unique solution is proven using the Banach fixed point theorem. In addition, the reader could use the successive approximate method (Picard method) to arrive at the same conclusion. In the homogeneous case of Equation (1), the successive approximate method fails to prove the existence of a unique solution. For this, we can only use the Banach fixed point theorem.
2
If the two conditions of (i) and (ii) are not satisfied, this means that we have at least one solution. In this case, we would use one of the following theorems: Brouwer fixed point theorem or Schauder fixed point theorem.
3
Using TMM, we have an NAS where the coefficient of the nonlinear term is a function of time. Hence, the existence of a unique solution for the NAS is discussed in the space × C [ 0 , T ] , T < 1 .
4
The fractional nonlinear mixed integro-differential Equation (1), under certain relations of μ and v, represents the nonlinear integral equation of the fractional phase-lag term
v Φ ( x , t + δ t ) = g ( x , t ) + λ Ω k ( | x y | ) Φ m ( y , t ) d y , μ = v ( δ t ) α Γ ( α ) .
The delaying or advancing of time reveals the natural phenomena, especially in the presence of thermoelectricity and magnetic media. Some of applications of fractional integro-differential equations are found in physics, chemistry, economics, and biology [12,29]. Equation (59) explains the physical meaning of the fractional equation of time as the first fractional approximation of the time lag equation, and this lag may be before or after real time.
5
The variable separation technique used in Equation (11) enabled the researchers to find the necessary time relationship between the nucleus and time, as presented in Equation (15).
6
The significance of the logarithmic kernel was approved from its derivatives f with these cases:
(a)
x k ( | y x | ) = 1 | y x | Cauchy kernel.
(b)
2 x 2 k ( | y x | ) = 1 | y x | 2 Strong singular kernel
(c)
The Carleman function was also established as:
ln | y x | = ( ln | y x | ) | y x | v U ( y , x ) | y x | v
where U ( y , x ) is a continuous function.
7
When the kernel of the equation was in the logarithmic function form k ( | x y | ) = ln ( | x y | ) , the relative error increased with increasing time. It was also noted that the error in the non-linear case was slightly larger than in the linear case.
8
In Example (2), when the kernel took the Carleman function k ( | x y | ) = x y β ( 0 < β < 1 ) , we noticed that the behavior of the error when increasing time was the same as that of the logarithmic function. However, by comparison, we found that at small times, the error in the logarithmic function was higher than in the Carleman function. With increasing time, we find that the relative error in the Carleman function is higher than its counterpart in the logarithmic function.
9
In Example (3), the error behavior of the Hilbert kernel k ( | x y | ) = cot x y 2 was the same as that of the logarithmic form and Carleman function.

11. Future Work

Future work will attempt to solve Equation (1) when the coefficients of the equation are variable. This will lead to solutions for many applications in the sciences related to nonlinear elasticity.

Author Contributions

Supervision, M.A.A.; Project administration, S.E.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the project number: 22UQU4282396DSR01, through the Deanship for Research & Innovation, Ministry of Education in Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Acknowledgments

The authors would like to thank the reviewers for their suggestions that helped improve the research. The authors thank the Deanship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number: 22UQU4282396DSR01.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. El-Sayed, A.M.A.; Gamal Ahmed, R. Infinite point and Riemann Stieltjes integral conditions for an integro differential equation. Nonlinear Anal. Model. Control 2019, 24, 733–754. [Google Scholar] [CrossRef]
  2. El-Sayed, A.M.A.; Gamal Ahmed, R. Solvability of a coupled system of functional integro-differential equations with infinite point and Riemann-Stieltjes integral conditions. Appl. Math. Comput. 2020, 370, 124918. [Google Scholar] [CrossRef]
  3. Mohamed, D.S. Application of Lerch polynomials to approximate solution of singular Fredholm integral equations with Cauchy kernel. Appl. Math. Inf. Sci. 2022, 16, 565–574. [Google Scholar]
  4. Hafez, R.M.; Youssri, Y.H. Spectral Legendre-Chebyshev Treatment of 2D linear and nonlinear mixed Volterra-Fredholm integral equation. Math. Sci. Lett. 2020, 9, 37–47. [Google Scholar]
  5. El-Sayed, S.M.; Abdel-Aziz, M.R. A comparison of Adomian’s decomposition method and wavelet-Galerkin method for solving integro-differential equations. Appl. Math. Comput. 2003, 136, 151–159. [Google Scholar] [CrossRef]
  6. Ali, M.R.; Mousa, M.M.; Ma, W.-Z. Solution of nonlinear Volterra integral equations with weakly singular kernel by using the HOBW method. Adv. Math. Phys. 2019, 2019, 1–10. [Google Scholar] [CrossRef]
  7. Klaseboer, E.; Sun, Q.; Chan, D.Y. Non-singular boundary integral methods for fluid mechanics applications. J. 525 Fluid Mech. 2012, 696, 468–478. [Google Scholar] [CrossRef]
  8. Ata, K.; Sahin, M. An integral equation approach for the solution of the Stokes flow with Hermite surfaces. Eng. Anal. Bound. Elem. 2018, 96, 14–22. [Google Scholar] [CrossRef]
  9. Kuzmina, K.; Marchevsky, I. The boundary integral equation solution in vortex methods with the airfoil surface line discretization into curvilinear panels. Proc. Top. Probl. Fluid Mech. 2019, 2019, 131–138. [Google Scholar]
  10. Lienert, M.; Tumulka, R. A new class of Volterra-type integral equations from relativistic quantum physics. J. Integral Equations Appl. 2019, 31, 535–569. [Google Scholar] [CrossRef]
  11. Sidorov, D. Integral Dynamical Models: Singularities, Signals and Control; World Scientific: Singapore, 2014. [Google Scholar]
  12. Abdou, M.A.; Basseem, M. Thermopotential function in position and time for a plate weakened by curvilinear hole. Arch. Appl. Mech. 2022, 92, 867–883. [Google Scholar] [CrossRef]
  13. Matoog, R.T. Treatments of probability potential function for nuclear integral equation. J. Phys. Math. 2017, 8, 2090-0902. [Google Scholar]
  14. Gao, J.; Condon, M.; Iserles, A. Spectral computation of highly oscillatory integral equations in laser theory. J. Comput. Phys. 2019, 395, 351–381. [Google Scholar] [CrossRef]
  15. Alhazmi, S.E.A. New Model for Solving Mixed Integral Equation of the First Kind with Generalized Potential Kernel. J. Math. Res. 2017, 9, 18–29. [Google Scholar] [CrossRef]
  16. Nemati, S.; Lima, P.M.; Ordokhani, Y. Numerical solution of a class of two-dimensional nonlinear Volterra integral equations using Legendre polynomials. J. Comput. Appl. Math. 2013, 242, 53–69. [Google Scholar] [CrossRef]
  17. Mirzaee, F.; Samadyar, N. Convergence of 2D-orthonormal Bernstein collocation method for solving 2D-mixed Volterra-Fredholm integral equations. Trans. A. Razmadze Math. Inst. 2018, 172, 631–641. [Google Scholar] [CrossRef]
  18. Basseem, M.; Alalyani, A. On the Solution of Quadratic Nonlinear Integral Equation with Different Singular Kernels. Math. Probl. Eng. 2020, 2020, 1–7. [Google Scholar] [CrossRef]
  19. Katani, R. Numerical solution of the Fredholm integral equations with a quadrature method. SeMA J. 2019, 76, 271276. [Google Scholar] [CrossRef]
  20. 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, 1–12. [Google Scholar] [CrossRef]
  21. 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]
  22. 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]
  23. 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]
  24. 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]
  25. Jebreen, B.; Dassios, I. On the Wavelet Collocation Method for Solving Fractional Fredholm Integro-Differential Equations. Mathematics 2022, 10, 1272. [Google Scholar] [CrossRef]
  26. 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]
  27. 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]
  28. Podlubny, I.; Chechkin, A.; Skovranek, T.; Chen, Y.; Jara, B.M.V. Matrix approach to discrete fractional calculus II: Partial fractional differential equations. J. Comput. Phys. 2009, 228, 3137–3153. [Google Scholar]
  29. Kreyszig, E. Introductory Functional Analysis with Applications; John Wiley & Sons Inc.: New York, NY, USA, 1978. [Google Scholar]
  30. Jan, A.R. An asymptotic model for solving mixed integral equation in position and time. J. Math. 2022, 2022, 8063971. [Google Scholar] [CrossRef]
  31. Jan, A.R. Solution of nonlinear mixed integral equation via collocation method basing on orthogonal polynomials. Heliyen 2022, 8, e1182. [Google Scholar] [CrossRef]
  32. Abdou, M.A.; El-Borai, M.M.; El-Kojok, M.M. Toeplitz matrix method and nonlinear integral equation of Hammerstein type. J. Comput. Appl. Math. 2009, 223, 765–776. [Google Scholar]
  33. Gu, X.M.; Huang, T.Z.; Zhao, X.L.; Xu, W.R.; Li, H.B.; Li, L. Circulant preconditioned iterative methods for peridynamic model simulation. Appl. Math. Comput. 2014, 248, 470–479. [Google Scholar] [CrossRef]
  34. Gradstein, I.S.; Ryzhik, I.M.; Romer, R.H. Tables of Integrals, Series, and Products; Academic Press: Cambridge, MA, USA, 1988. [Google Scholar]
  35. Arutiunian, N.K. Ploskaia kontaktnaia zadacha teorii polzuchesti (The plane contact problem of the theory of creep). J. Appl. Math. Mech. 1959, 23, 1283–1313. [Google Scholar] [CrossRef]
Figure 1. (a) Shows the corresponding error of NMIE with a logarithmic kernel at different times, whereas (b) illustrates the corresponding error for the linear case at the same times.
Figure 1. (a) Shows the corresponding error of NMIE with a logarithmic kernel at different times, whereas (b) illustrates the corresponding error for the linear case at the same times.
Fractalfract 07 00173 g001
Figure 2. (a) Shows the error for the nonlinear case at β = 0.01 and (b) shows the corresponding error for the linear case at β = 0.01 .
Figure 2. (a) Shows the error for the nonlinear case at β = 0.01 and (b) shows the corresponding error for the linear case at β = 0.01 .
Fractalfract 07 00173 g002
Figure 3. (a,b) Shows the error of Equation (44) respectively for the nonlinear and linear cases of Carleman coefficients at t = 0.2 and N = 20 .
Figure 3. (a,b) Shows the error of Equation (44) respectively for the nonlinear and linear cases of Carleman coefficients at t = 0.2 and N = 20 .
Fractalfract 07 00173 g003
Figure 4. (a,b) Shows the error of Equation (44) using the Hilbert kernel for the nonlinear and the linear cases at time t = 0.1 , N = 20 .
Figure 4. (a,b) Shows the error of Equation (44) using the Hilbert kernel for the nonlinear and the linear cases at time t = 0.1 , N = 20 .
Fractalfract 07 00173 g004
Figure 5. (a,b) Shows the error of Equation (44) with Hilbert kernel for the nonlinear and the linear cases at time t = 0.4 , N = 20 .
Figure 5. (a,b) Shows the error of Equation (44) with Hilbert kernel for the nonlinear and the linear cases at time t = 0.4 , N = 20 .
Fractalfract 07 00173 g005
Table 1. Solutions of the linear/nonlinear solutions of MIE (44) with a logarithmic kernel, along with the corresponding errors using TMM.
Table 1. Solutions of the linear/nonlinear solutions of MIE (44) with a logarithmic kernel, along with the corresponding errors using TMM.
TimeExact SolutionNonlinear Case m = 2, N = 20Linear Case m = 1, N = 20
t φ Exact φ Non . E Non . φ Lin . E Lin .
1.0001000 1.00010270 0.269 × 10 7 1.00010561 0.567 × 10 7
t = 0.01 0.9345000 0.93450708 0.362 × 10 7 0.93450758 0.384 × 10 7
0.4012111 0.40122142 0.496 × 10 7 0.40122136 0.187 × 10 7
0.0901000 0.09010238 0.976 × 10 7 0.09010146 0.705 × 10 7
1.0100000 1.01003928 0.286 × 10 6 1.01006525 0.525 × 10 6
t = 0.1 0.9444440 0.94444790 0.324 × 10 6 0.94444837 0.392 × 10 6
0.6500000 0.65000474 0.297 × 10 6 0.65000101 0.898 × 10 6
0.5977778 0.59779381 0.160 × 10 6 0.59778257 0.172 × 10 6
t = 0.4 1.1600000 1.60004729 0.292 × 10 5 1.16000215 0.215 × 10 5
1.09444440 1.09444173 0.434 × 10 5 1.09444627 0.182 × 10 5
0.7477778 0.74772763 0.288 × 10 5 0.74779495 0.982 × 10 6
0.6977778 0.69777818 0.341 × 10 6 0.69777785 0.231 × 10 5
t = 0.8 1.640000 1.64007850 0.211 × 10 5 1.64001375 0.137 × 10 5
1.177778 1.17784366 0.588 × 10 5 1.17781185 0.340 × 10 5
0.711111 0.71115465 0.559 × 10 5 0.71118519 0.740 × 10 5
Table 2. Describes the nonlinear and linear cases for MIE (44), with the Carleman kernel, time and β = 0.01 .
Table 2. Describes the nonlinear and linear cases for MIE (44), with the Carleman kernel, time and β = 0.01 .
TimeExact SolutionNonlinear Case m = 2, N = 20Linear Case m = 1, N = 20
t φ Exact φ Non . E Non . φ Lin . E Lin .
t = 0.01 1.0001000 1.0001000 0.545 × 10 8 1.0001000 0.567 × 10 8
0.6945000 0.6945000 0841 × 10 8 0.6945000 0.956 × 10 8
0.0178777 0.0178777 0.708 × 10 8 0.0178777 0344 × 10 8
0.0278777 0.0278777 0.710 × 10 8 0.0278777 0.547 × 10 8
t = 0.1 1.0100000 1.0100000 0.456 × 10 7 1.0100000 0.648 × 10 7
0.9444444 0.9444448 0.465 × 10 7 0.9444447 0.375 × 10 7
0.6500000 0.6500065 0.342 × 10 7 0.6500004 0.173 × 10 7
0.2600000 0.2600008 0.179 × 10 7 0.2600008 0.160 × 10 7
t = 0.4 1.1600000 1.1600875 0.700 × 10 5 1.1600075 0.750 × 10 5
0.8000000 0.8000373 0.753 × 10 5 0.8000290 0.251 × 10 5
0.6044444 0.604439 0.194 × 10 5 0.6044499 0.244 × 10 5
0.2711111 0.2711951 0.832 × 10 5 0.2711990 0.612 × 10 5
t = 0.8 1.6400000 1.6400332 0.331 × 10 4 1.6400542 0.542 × 10 4
1.0844444 1.0849735 0.549 × 10 4 1.0849810 0.416 × 10 4
0.7300000 0.7304171 0.416 × 10 4 0.73082523 0.311 × 10 4
0.6444444 0.6443994 0.053 × 10 4 0.64440364 0.170 × 10 5
Table 3. Shows the nonlinear and linear cases for Equation (44) for different values of Carleman coefficients at time t = 0.2 , N = 20 .
Table 3. Shows the nonlinear and linear cases for Equation (44) for different values of Carleman coefficients at time t = 0.2 , N = 20 .
Carleman CoefficientsExact Solution at t = 0.2 Nonlinear Case at m = 2Linear Case at m = 1
β φ Exact φ Non . E Non . φ Lin . E Lin .
0.0400000 0.0400099 0.519 × 10 6 0.04003999 0.199 × 10 6
β = 0.22 0.0096040 0.00960747 0.467 × 10 6 0.00965943 0.624 × 10 6
0.4000000 0.40000256 0.166 × 10 6 0.40000373 0.373 × 10 6
0.0400000 0.04000081 0.158 × 10 6 0.0400006 0.167 × 10 6
0.0400000 0.04000220 0.822 × 10 5 0.04003999 0.890 × 10 5
β = 0.32 0.0163840 0.01632257 0.798 × 10 5 0.01631233 0.466 × 10 5
0.0010240 0.00102484 0.364 × 10 5 0.00102149 0.260 × 10 5
0.0400000 0.04000523 0.220 × 10 5 0.04000061 0.111 × 10 5
0.0400000 0.04030011 0.434 × 10 4 0.04000210 0.210 × 10 5
β = 0.7 0.0010240 0.00100645 0.754 × 10 5 0.00100266 0.262 × 10 5
0.0163840 0.01638701 0.197 × 10 5 0.01638794 0.457 × 10 5
0.0400000 0.04003430 0.430 × 10 5 0.04999902 0.973 × 10 6
0.0400000 0.04000037 0.375 × 10 5 0.04000021 0.219 × 10 6
β = 0.8 0.0096040 0.00960777 0.237 × 10 5 0.00960137 0.974 × 10 6
0.0000640 0.00006437 0.243 × 10 5 0.00006407 0.292 × 10 5
0.0163840 0.016385043 0.112 × 10 5 0.01638486 0.214 × 10 5
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

Alhazmi, S.E.; Abdou, M.A. A Physical Phenomenon for the Fractional Nonlinear Mixed Integro-Differential Equation Using a General Discontinuous Kernel. Fractal Fract. 2023, 7, 173. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7020173

AMA Style

Alhazmi SE, Abdou MA. A Physical Phenomenon for the Fractional Nonlinear Mixed Integro-Differential Equation Using a General Discontinuous Kernel. Fractal and Fractional. 2023; 7(2):173. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7020173

Chicago/Turabian Style

Alhazmi, Sharifah E., and Mohamed A. Abdou. 2023. "A Physical Phenomenon for the Fractional Nonlinear Mixed Integro-Differential Equation Using a General Discontinuous Kernel" Fractal and Fractional 7, no. 2: 173. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract7020173

Article Metrics

Back to TopTop