Next Article in Journal
A Robust Zn-Hydroxamate Metal–Organic Framework Constructed from an Unsymmetrical Ligand for Iodine Capture
Next Article in Special Issue
Criteria of Oscillation for Second-Order Mixed Nonlinearities in Dynamic Equations
Previous Article in Journal
Anti-Fatigue-Damage-Oriented Through-Life Optimization and Control of High-Power IGCT Converters in Wind Energy Systems
Previous Article in Special Issue
New Monotonic Properties for Solutions of Odd-Order Advanced Nonlinear Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis of the Solution for the Mixed Integral Equation with Symmetric Kernel in Position and Time with Its Applications

by
Faizah M. Alharbi
Mathematics Department, Faculty of Sciences, Umm Al-Quraa University, Makkah 24227, Saudi Arabia
Symmetry 2024, 16(8), 1048; https://doi.org/10.3390/sym16081048
Submission received: 16 July 2024 / Revised: 8 August 2024 / Accepted: 10 August 2024 / Published: 14 August 2024

Abstract

:
Under certain assumptions, the existence of a unique solution of mixed integral equation (MIE) of the second type with a symmetric kernel is discussed, in L 2 [ Ω ] × C 0 , T , T < 1 , Ω is the position domain of integration and T is the time. The convergence error and the stability error are considered. Then, after using the separation technique, the MIE transforms into a system of Hammerstein integral equations (SHIEs) with time-varying coefficients. The nonlinear algebraic system (NAS) is obtained after using the degenerate method. New and special cases are derived from this work. Moreover, numerical results are computed using MATLAB R2023a software.

1. Introduction

It is known that integral equations are the structural origin of differential equations. Therefore, it is natural and known that to solve a differential equation or transform it into an integral equation, the solution outcome depends on the given conditions. If the given conditions are initial, the result represents the Volterra integral equation. If the conditions are marginal, the result represents the Fredholm integral equation. If the conditions are mixed, we obtain the Volterra–Fredholm integral equation. We find these types of problems frequently in contact problems in elasticity science or mixed problems in communication mechanics and other applied sciences.
The versatility of integral equations (IEs), with their various forms, permits the simulation of a broad array of issues in the realm of fundamental sciences, prompting a considerable effort by researchers to present the solutions to these systems. Scientists have used many different ways to solve these problems, which shows that they are getting better quickly.
The orthogonal polynomial approach is well recognized as a crucial tool for solving diverse scientific challenges encountered in numerous fields of study. Alhazmi [1] employed a novel method to extract several spectral connections from the mixed integral equation. He did this by combining the generalized potential kernel with variable separation and orthogonal polynomials. Nemati et al.’s study [2] used a novel method that included a continuous kernel to examine the output of the two-dimensional (2D) Volterra integral model. Hafez and Youssri [3] examined the numerical solutions of the two-dimensional integral model, which is founded on the linear Volterra–Fredholm equation. They employed the collocation method using the Legendre–Chebyshev polynomials. In their study, Abdou et al. [4] employed Chebyshev polynomials to effectively address quadratic integral equations that incorporate a phase-lag factor in the time domain.
Presently, numerous studies have focused on developing advanced and efficient methods for solving the integral equations (IEs) and integral differential equations (Io-DEs). These include the utilization of Riemann–Stieltjes integral conditions [5,6], the application of the Lerch polynomial method [7,8], and the implementation of the Spectral Legendre–Chebyshev method [9]. Additionally, numerical observations based on semi-analytical approaches in [10] have been employed to solve the VIEs. The wavelet collocation method was used [11] for solving fractional integro-differential equations (FrIo-DEs). In [12], the authors applied the extended cubic B-spline to interpret the collocation strategy to solve the FrIo-DEs. Zhang and Li, in [13], establish the basic structure of the exponential Euler difference form for Caputo–Fabrizio fractional-order differential equations with multiple lags. This kind of study provides a deep understanding of behavior for modifying systems.
The linear and nonlinear integral equations (LIEs/NIEs) are utilized in airfoil analysis, quantum mechanics, astrophysics, and lasers; see [14,15,16,17,18,19,20].
Consider the NIE of the second kind of type Volterra–Hammerstein in (n + 1) dimensional
μ U ( u , t ) = H ( u , t ) + λ 0 t Ω f t , τ p u , y Ζ τ , w , U w , τ d w d τ , ( μ , λ c o n s t a n t s ) ( u = u ¯ ( u 1 , u 2 , . . . , u n ) , w = w ¯ ( w 1 , w 2 , . . . , w n ) )
Here, in Equation (1), the given functions H u , t , Ζ ( t , u , U ( u , t ) ) belong to L 2 ( Ω ) × C 0 , T —space. The integration domain Ω is a closed bounded set and depends on the vector of position, while the time t 0 , T . The kernel p ( u , w ) L 2 ( [ Ω ] × [ Ω ] ) , while f ( t , τ ) C ( [ 0 , T ] × [ 0 , T ] ) ; t , τ 0 , T , T < 1 . The constant μ describes the type of the integral equation, and λ carries a physical significance. Finally, U ( u , t ) is the unknown function and will be discussed in the space of the given function H ( u , t ) .
Differentiating Equation (1) with respect to the variable t, we obtain
μ U u , t t = H u , t t + λ f t , t Ω p u , w Ζ t , w , U w , t d w + λ 0 t Ω f ( t , τ ) t p ( x , y ) Z ( τ , w , U ( w , τ ) ) d w d τ .
The integro-differential Equation (Io-DE) (2) is equivalent to IE (1). Therefore, the two equations have the same solution after neglecting the constant term.
In this article, the author is interested in studying the effect of time on the integral equation, so that the position part is linked to boundary conditions, while the time part is linked to an initial condition. There is little research that has linked position and time in integral equations. Therefore, the interest in this research is to study a nonlinear integral equation linked between position and time. In addition, this work examines and establishes the existence of a solitary solution for NIE (1), subject to specific circumstances. Furthermore, the analysis includes the examination of convergence and the assessment of error stability. A numerical technique was applied to derive a set of location-wise Hammerstein integral equations (SHIEs). In addition, the degenerate kernel method is employed to calculate the numerical solution of the SHIEs, resulting in a nonlinear algebraic system (NAS) of equations. A unique solution is obtained. Ultimately, a multitude of situations with varying kernels was resolved.

2. The Existence Solution of the NV-HIE

A fixed-point theorem and the successive approximations approach will be used to prove the existence of a single solution to Equation (1) and then converge this solution as well as the error. Therefore, we represent it in the structure of an integral operator,
W ¯ U ( u , t ) = 1 μ H ( u , t ) + 1 μ W U ( u , t ) , ( μ 0 )
W U ( u , t ) = λ 0 t Ω f ( t , τ ) p ( u , w ) Z ( τ , w , U ( w , τ ) ) d w d τ .
Then, assume the assumptions:
(i) Position kernel p ( u , w ) , u = u ¯ ( u 1 , u 2 , . . . , u n ) , w = w ¯ ( w 1 , w 2 , . . . , w n ) , check the offline condition, Ω Ω p ( u , w ) 2 d u d w 1 2 = M , (M-constant).
(ii) Time kernel f ( t , τ ) C 0 , T check status f ( t , τ ) f ¯ , ( f ¯ —constant).
(iii) The given function H ( u , t ) belongs to L 2 ( Ω ) × C 0 , T —space, and its norm is
H ( u , t ) L 2 ( Ω ) × C 0 , T = m a x 0 t T 0 t Ω H 2 ( u , τ ) d u 1 2 d τ = G , ( G constant )
(iv) The function Z t , u , U u , t , for the constants Q > L , Q > Q 1 , satisfies the following conditions:
( a ) m a x 0 t T 0 t Ω Z ( τ , u , U ( u , τ ) ) 2 d u 1 2 d τ Q 1 U ( u , t ) L 2 ( Ω ) × C 0 , T , U ( u , t ) = m a x 0 t T 0 t Ω U ( u , τ ) 2 d u 1 2 d τ .
b Z t , u , U 1 u , t Z t , u , U 2 u , t N t , u U 1 u , t U 2 u , t , N ( t , u ) L 2 ( Ω ) × C 0 , T = m a x 0 t T 0 t Ω N 2 ( τ , u ) d u 1 2 d τ = L , ( L c o n s t a n t ) .
Theorem 1.
If the assumptions (i)–(iv) are verified, then Equation (1) has a single solution in L 2 ( Ω ) × C 0 , T , under the condition
μ > λ f ¯ M Q T
In order to establish the validity of this theorem, it is necessary to examine the following lemmas:
Lemma 1.
By the assumptions (i)–(iv-a), W ¯ transforms L 2 ( Ω ) × C 0 , T into itself.
Proof. 
Based on Equations (3) and (4), we obtain
W ¯ U ( u , t ) 1 μ H ( u , t ) + λ μ 0 t Ω f ( t , τ ) P ( u , w ) Z ( τ , w , U ( w , τ ) ) d w d τ .
By utilizing the assumptions (ii)–(iii) and subsequently applying the Cauchy Schwarz inequality, we obtain,
W ¯ U ( u , t ) G μ + λ μ f ¯ (   m a x 0 t d t 0 t T   )   Ω Ω p ( u , w ) 2 d u d w 1 2 . m a x 0 t T 0 t Ω Z ( τ , w , U ( w , τ ) 2 d w 1 2 d τ .
Considering the given constraints (i), (iv-a), we have,
W ¯ U ( u , t ) G μ + σ U ( u , t ) , ( σ = λ μ   f ¯ M Q T ) .
Hence, W ¯ maps the ball S ρ into itself where
ρ = G [ μ λ f ¯ M Q T ] .
the positive values of ρ , G , lead to μ λ f ¯ M Q T > 0 . This leads us to say that σ < 1 .
Furthermore, W is bounded where
W U ( u , t ) σ U ( u , t ) .
Lemma 2.
Under (i), (ii), and (iv-b), W ¯ is contraction operator in L 2 ( Ω ) × C 0 , T .
Proof. 
Regarding the two functions { U 1 ( u , t ) , U 2 ( u , t ) } L 2 ( Ω ) × C 0 , T , then from Equations (3) and (4), we find
W ¯ ( U 1 u , t U 2 x , t )   λ μ 0 t Ω f ( t , τ ) p ( u , w ) Z ( τ , w , U 1 ( w , τ ) ) Z ( τ , w , U 2 ( w , τ ) ) d w d τ .
With the aid of (i), (ii), and (iv-b), and Cauchy–Schwarz inequality, we obtain
W ¯ ( U 1 ( u , t ) U 2 ( u , t ) ) σ U 1 ( u , t ) U 2 ( u , t ) .
Hence, W ¯ is continuous in L 2 ( Ω ) × C 0 , T . Moreover, under the condition σ < 1 ,   W ¯ is a contraction operator. Hence, W ¯ has a single fixed point, which is the unique solution of Equation (1). □

3. Results the Convergence and the Error Stability of Solution

The successive approximation method was used to prove the convergence and error stability of the solution. For this, we assume U n u , t U i u , t } i = 0 , to have
μ U n u , t = H u , t                         + λ 0 t Ω f ( t , τ ) p ( u , w ) Z ( τ , w , U n 1 ( w , τ ) ) d w d τ , ( n 1 , U 0 ( u , t ) = H ( u , t ) )
set,
Ψ n ( u , t ) = U n ( u , t ) U n 1 ( u , t )
to have
U n ( u , t ) = i = 1 n Ψ i ( u , t ) , Ψ 0 ( u , t ) = H ( x , t ) .
Equation (1), yields
μ U n ( u , t ) U n 1 ( u , t )                                             λ 0 t Ω f ( t , τ ) p ( u , w ) | | Z ( τ , w , U n 1 ( w , τ ) )                                             Z ( τ , w , U n 2 ( w , τ ) ) d w d τ .
Using assumptions (ii)–(iv-a), to have
μ U n ( u , t ) U n 1 ( u , t ) λ f ¯ 0 t Ω p ( u , w ) N ( τ , w ) U n 1 ( w , τ ) U n 2 ( w , τ ) ) d w d τ .
By utilizing the Cauchy–Schwarz inequality on the Hammerstein integral term, and considering the assumptions (i) and (iv-b), the inequality is simplified to
Ψ n ( u , t ) σ Ψ n 1 ( u , t ) , ( n 1 ) .
When n = 1 , the inequality (12), with condition (ii) takes the form
Ψ 1 ( u , t ) σ G .
By induction, we write
Ψ n ( u , t ) σ n G , n = 0,1 , . . . ,
Since Equation (14), for n = 0,1 , . . . , is obviously true then the sequence U n ( u , t ) converges, and
U ( u , t ) = i = 0 Ψ i ( u , t ) .
The series (15) exhibits uniform convergence due to the uniform behavior of its terms. Ψ i ( u , t ) are dominated by σ i .
To discuss the error stability, we assume U n ( u , t ) is the numerical solution of Equation (1); hence, we have
μ R u , t = H u , t H n ( u , t )                                                 + λ 0 t Ω f t , τ p u , w Z τ , w , U w , τ Z n ( τ , w , U n ( w , τ ) d w d τ
Following the same previous way, we have
R ( u , t ) H ( u , t ) H n ( u , t ) μ λ f ¯ M Q T .
The inequality (17) proved that the error is stable and has a unique representation under the condition,
μ > λ f ¯ M Q T .

4. System of Hammerstein Integral Equations (SHIEs) in Position

In this part, Equation (1) employs a numerical approach to derive SHIEs in position. In order to achieve this objective, we divide 0 , T as t = t k , k = 0,1 , 2 , . . . , N , then from Equation (1), we have (see [2,20])
0 t k Ω f ( t k , τ ) p ( x , y ) Z ( τ , w , U ( w , τ ) ) d w d τ = j = 0 k ω j f ( t k , t j ) Ω p ( u , w ) Z ( t j , w , U ( w , t j ) ) d w + O ( k p + 1 ) , P > 0 , , k = m a x 0 j k h j , h j = t j + 1 t j , ω k = 1 2 h k , ω j = h j , j = 0 , k .
Using Equation (18) in Equation (1), and neglecting O ( k p + 1 ) , we have
μ U k ( u ) = H k ( u ) + λ j = 0 k ω j f j , k Ω p ( u , w ) Z j ( w , U j ( w ) ) d w .
Here, we used the notations,
U u , t k = U k u , H u , t k = H k u , f t k , t j = f j , k , Z ( t j , u , U ( u , t j ) ) = Z j ( u , U j ( u ) )
Equation (19) represents a set of n-dimensional HIEs, and its solution relies on the specified function H k ( u ) , k = 0,1 , . . . , N , the kernel p (u, w), and the known function Z j ( u , U j ( u ) ) .

5. Special Cases

(i) As an important special case let, in (19), k = 0 to have
μ U 0 ( u ) = H 0 ( u ) + λ ω 0 f 0,0 Ω p u , w Z 0 w , U 0 w d w .
Equation (20) represents an integral equation of Hammerstein type and its solution depends on the kind of kernel.
(ii) Let u = u ¯ u 1 , w = w ¯ w 1 , Z 0 w , U 0 w = U 0 w , Ω = 1,1 , t o h a v e
μ U 0 u = H 0 u + ζ 1 1 p u , w U 0 w d w , ζ = λ ω 0 f 0,0
Formula (21) represents a Fredholm integral equation of the second kind. The importance of the kernel when it takes different cases of a singular term like logarithmic kernel, Carleman kernel, Cauchy kernel, and strong singular kernel,
As an important special case when in (21) p u , w = l n u v (logarithmic kernel) hence, we have
μ U 0 u = H 0 u + ζ 1 1 l n u v U 0 w d w .
Differentiating the above equation with respect to u, we obtain
μ d U 0 u d u = d H 0 u d u + ζ 1 1 1 u v U 0 w d w
Using the substation u = 2 x 1 , v = 2 y 1 , we have
d θ x d x = g x + ζ ¯ 0 1 1 u v U 0 w d w , { g x = d H 0 2 x 1 2 μ d x , ζ ¯ = ζ μ }
The integral Equation (22) has appeared in both combined infrared gaseous radiation and molecular conduction, see [20].
(iii) Three-dimensional integral equation
Let u = u ¯ u 1 , u 2 , u 3 , w = w ¯ w 1 , w 2 , w 3 , Z 0 w , U 0 w = U 0 w , we have
μ U 0 u , w = H 0 u , w + β Ω p u ζ , w η U 0 ζ , η d ζ d η , β = λ ω 0 f 0,0
here Ω   Ω = { u , w , v ϵ Ω : u 2 + w 3 2 a , v = 0 .
If p u ζ , w η = [ u ζ 2 + w + η 2 ] υ , 0 ν < 1 , the integral Equation (22) was nvestigated from the semi-symmetric Hertz problem for two different elastic materials in three dimensions when the modules of elasticity change according to the power law σ i = K 0 ε i ν   ( 0 ν 1 ) ,where σ i and ε i , i = 1,2 , 3 are the stress and strain rate intensities, respectively, while K 0 , ν are physical constants.

6. The Degenerate Kernel Method and the Nonlinear Algebraic System

The straightforwardness of utilizing the Degenerate method to solve the nonlinear algebraic integral system (19) compels one to contemplate substituting the provided kernel p ( u , w ) approximately by a kernel p n ( u , w ) ; that is
p n ( u , w ) = i = 1 n B i ( u ) C i ( w ) .
The set of functions {Bi(u)} and {Ci(w)} are linearly independent, such that
Ω Ω p ( u , w ) p n ( u , w ) 2 d u d w 1 2 0   a s   n
hence, the solution of Equation (19) associated to the kernel p n ( u , w ) takes the form
μ U n , k ( u ) = H k ( u ) + λ j = 0 k ω j f j , k Ω p n ( u , w ) Z j ( w , U n , j ( w ) ) d w .
Using Equation (20) in Equation (22), we have
μ U n . k ( u ) = H k ( u ) + λ i = 1 n j = 0 k ω j f ( j , k ) A i , j B i ( u ) , ( μ 0 ) , k = 0,1 , 2 , . . . , N .
where
A i , j = Ω C i ( w ) Z j ( w , U n , j ( w ) ) d w , j = 0,1 , 2 , . . . , k
here, A i , j s are constants to be determined from the following formula
A m , j = Ω C m ( w ) Z j ( w , 1 μ H j ( w ) + λ μ i = 1 n r = 0 j ω r f r , j A i , j B i ( w ) ) d w , ( m = 1,2 , . . . , n ) .
Define
Ε m , j ( A 1 , j , A 2 , j , . . . , A n , j ) = Ω C m ( w ) Z j ( w , 1 μ H j ( w ) + λ μ i = 1 n r = 0 j ω r f r , j A i , j B i ( w ) ) d w ,
Equation (28) represents a system of NAEs that can be written as a matrix equation
A 1 , j A 2 , j A 3 , j . . . A n , j = Ε 1 , j ( A 1 , j , A 2 , j , . . . , A n , j ) Ε 2 , j ( A 1 , j , A 2 , j , . . . , A n , j ) Ε 3 , j ( A 1 , j , A 2 , j , . . . , A n , j ) . . . Ε n , j ( A 1 , j , A 2 , j , . . . , A n , j ) ,
the nonlinear algebraic system (29) can be solved numerically.

The Existence of a Single Solution of the Nonlinear Algebraic System

This section will provide an exposition on the existence and validation of a distinct NAS solution (24). To accomplish this, the following theorem will be formulated.
Theorem 2.
Assume that the known continuous functions Z j ( w , ψ ( w , A i , j ) ) in Equation (27) satisfy the following conditions,
Ω Z j ( w , ψ ( w , A i , j ) ) 2 d w 1 2 L i = 1 n A i , j 2 1 2 , ( L   i s   c o n s t a n t   )
and,
Z j ( w , ψ ( w , A i , j ) ) Z j ( w , ψ ( w , D i , j ) ) M 1 ψ w , A i , j ψ w , D i , j , ( M 1   i s   c o n s t a n t )
then, the NAS (29) has a single solution   A ¯ j , and U n , k ( u ) is the single solution of Equation (27) in l 2 -space.
To demonstrate the validity of this theorem, it is necessary to examine the following two lemmas.
Lemma 3.
By the aid of the condition (30), the operator E ¯ j of Equation (28) maps l 2 -space into itself.
Proof. 
Let V be the set of functions Ξ = ξ i in l 2 such that
Ξ l 2 = i = 1 ξ i 2 1 2 β , ( β constant ) .
From Equation (27), we have
E m , j ( A 1 , j , A 2 , j , . . . , A n , j ) Ω C m ( w ) w , 1 μ H j ( w ) + λ μ i = 1 n r = 0 j ω r f r , j A i , j B i ( w ) ) ) d w .
Hence, after applying Cauchy–Schwarz inequality, and using condition (30), we have
m = 1 n E m , j ( A 1 , j , A 2 , j , . . . , A n , j ) 2 1 2 M 2 i = 1 n A i , j 2 1 2 , ( M 2 L m = 1 n Ω C m ( w ) 2 d w 1 2 ) constant ) ,
as n , the last inequality yields
E ¯ j ( A ¯ j ) l 2 M 2 A ¯ j l 2 .
Hence, H ¯ i is a bounded operator that maps the set U into itself, where
β = M 2 A ¯ ( j ) l 2 .
Lemma 4.
Under the condition (31), E ¯ j is a contraction mapping in l 2 .
For the functions A ¯ j = ( A 1 , j , A 2 , j , . . . , A n , j ) and D ¯ j = ( D 1 , j , D 2 , j , . . . , D n , j ) in l 2 , Equation (24) leads to
E m , j ( A 1 , j , A 2 , j , . . . , A n , j ) E m , j ( D 1 , j , D 2 , j , . . . , D n , j ) Ω C m ( w ) Z j ( w , 1 μ H j ( w ) + λ μ i = 1 n r = 0 j ω r f r , j A i , j B i ( w ) ) Z j ( w , 1 μ H j ( w ) + λ μ i = 1 n r = 0 j ω r f r , j D i , j B i ( w ) ) d w .
Introducing condition (31), then applying Cauchy–Schwarz inequality three times, the above inequality takes the form
m = 1 n H m , j ( A 1 , j , A 2 , j , . . . , A n , j ) H m , j ( D 1 , j , D 2 , j , . . . , D n , j ) 2 1 2 M 3 i = 1 n A i , j D i , j 2 1 2 , M 3 = λ μ M 1 r = 0 j u ( r ) 2 1 2 r = 0 j F ( r , j ) 2 1 2 m = 1 n Ω C m ( y ) 2 d y 1 2 i = 1 n Ω B i ( y ) 2 d y 1 2 .
As n , the previous inequality can be reduced to
E ¯ j ( A ¯ j ) E ¯ j ( D ¯ j ) l 2 M 3 A ¯ j D ¯ j l 2 .
Thus, E ¯ j is a continuous operator in the space l 2 . If M 3 < 1 , then E ¯ j is a contraction operator. Hence, E ¯ j has a unique fixed point A ¯ j , which is the unique solution of the algebraic system (27). In view of Theorem 2, the algebraic system (27) has a single solution U n , k ( u ) in l 2 .

7. Examples and Numerical Results

Example 1.
For the NV-HIE
U ( u , t ) 0 t 0 1 τ 3 ( 1 + u w ) U l ( w , τ ) d w d τ = H ( u , t ) , l = 1,2 , . . . , N , ( U ( u , t ) = u t ) .
If we divide the interval [0, T] as 0 = t 0 < t 1 < t 2 < t 3 = t , t = t k ; k = 0,1 , 2,3 , and using the degenerate method, we obtain
U k ( u ) j = 0 k m = 0 n ω j t j 3 u m m ! A j , m = H k ( u ) , ( k = 0,1 , 2,3 ) ,
where
A k , m = 0 1 y m H k ( w ) + j = 0 k m = 0 n ω j t j 3 u m m ! A j , m 2 d w , H k ( u ) = H ( u , t k ) .
The solution of Equation (35) gives the following results,
U 0 ( u ) = 0 , U 1 ( u ) = u t 3 , U 2 ( u ) = 2 u t 3 , U 3 ( u ) = u t .
Here,   U 3 ( u ) is the exact solution of Equation (34).
Example 2.
For the NV-HIE
U ( u , t ) 0 t 0 1 τ 3 e u w U l ( w , τ ) d w d τ = H ( u , t ) , l = 1,2 , . . . , N , ( U ( u , t ) = u t ) .
Approximate the kernel of Equation (36) p ( u , w ) = e u w in the form,
p ( u , w ) p n ( u , w ) = m = 0 n ( u w ) m m ! ,
where
0 1 0 1 p u , w p n u , w 2 d u d w 1 2 0   a s   n .
Using the degenerate kernel, we obtain
U k ( u ) j = 0 k m = 0 n ω j t j 3 u m m ! A j , m = H k ( u ) , A k , m = 0 1 w m H k w + j = 0 k m = 0 n ω j t j 3 u m m ! A j ( m ) 2 d w , ( k = 0,1 , 2,3 ) .
The solution of Equation (38), gives U 0 ( u ) = 0 , U 1 ( u ) = u t 3 , U 2 ( u ) = 2 u t 3 , U 3 ( u ) = u t . Here, U 3 ( u ) is the exact solution of Equation (36).
Example 3.
Assume the kernel in the Legendre or Chebyshev polynomials forms. Then, we have two cases as follows,
Example 3.1.
Consider the MIE with the kernel in Legendre polynomial form,
U ( u , t ) = H ( x , t ) + 0 t 1 1 τ 2 n = 0 20 P n ( u ) P n ( w ) U l ( w , τ ) d w d τ ;
The exact solution is U ( u , t ) = u 2 + t 2 .
In the third example, the kernel generation method was applied when the kernel is in the form of a Legendre polynomial, considering the difference in time. The equation was also applied when it was linear and when it was nonlinear. From Table 1, Table 2 and Table 3, it was noted that the error in the linear case is slightly higher than the error in the nonlinear case. Also, as time increases, the cumulative error increases.
Example 3.2.
Consider the MIE with the kernel in the Chebyshev polynomial form,
U u , t = H u , t + 0 t 1 1 τ 2 n = 0 20 T n u T n w U l w , τ d w d τ ;
The exact solution is U ( u , t ) = u 2 + t 2 .
In this example, the kernel generation method was also applied in the linear and nonlinear integral equation at different times for the continuous time kernel. The position kernel was also imposed in the form of a Chebyshev polynomial function of the first kind. The resulting error was also studied at different times, and this is clear in Table 4, Table 5 and Table 6.

8. Conclusions

The study successfully proves the existence of a solitary solution for NMIE (1) under specified conditions. The research includes an analysis of convergence and an evaluation of error stability. To solve the problem, a numerical technique is utilized to transform the problem into a set of location-wise Hammerstein integral equations (SHIEs). The degenerate kernel method is a powerful technique to compute the numerical solutions of these SHIEs, leading to a nonlinear algebraic system (NAS) of equations that yields a unique solution. The study concludes by resolving numerous scenarios involving different kernels, including the kernels in the form of Legendre and Chebyshev polynomials, respectively, and calculating the numerical solution in the linear and nonlinear cases. We deduce that the error estimates associated with the nonlinear case are less than the linear ones; see Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6 and Figure 1, Figure 2, Figure 3 and Figure 4.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. 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]
  2. 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]
  3. 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]
  4. Abdou, M.A.; Soliman, A.; Abdel-Aty, M. Analytical results for quadratic integral equations with phase-clag term. J. Appl. Anal. Comput. 2020, 10, 1588–1598. [Google Scholar] [CrossRef] [PubMed]
  5. El-Sayed, A.M.A.; Ahmed, G.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]
  6. El-Sayed, A.M.A.; Ahmed, R.G. Infinite point and Riemann Stieltjes integral conditions for an integro differential equation. Nonlinear Anal. Model. Control 2019, 24, 733–754. [Google Scholar] [CrossRef]
  7. 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]
  8. Cayan, S.; Sezer, M.A. Novel study based on Lerch polynomials for approximate solutions of pure neumann problem. Int. J. Appl. Comput. Math. 2022, 8, 8. [Google Scholar] [CrossRef]
  9. Ali, M.R.; Mousa, M.M.; Ma, W.-X. Solution of nonlinear Volterra integral equations with weakly singular kernel by using the HOBW method. Adv. Math. Phys. 2019, 2019, 1705651. [Google Scholar] [CrossRef]
  10. 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. Comp. 2003, 136, 151–159. [Google Scholar] [CrossRef]
  11. Bin Jebreen, H.; Dassios, I. On the Wavelet Collocation Method for Solving Fractional Fredholm Integro-Differential Equations. Mathematics 2022, 10, 1272. [Google Scholar] [CrossRef]
  12. Mahdy, A.M.S.; Abdou, M.A.; Mohamed, D.S. Computational methods for solving higher-order (1 + 1) dimensional mixed difference integro-differential equations with variable coefficients. Mathematics 2023, 11, 2045. [Google Scholar] [CrossRef]
  13. Zhang, T.; Li, Y. Exponential Euler scheme of multi-delay Caputo–Fabrizio fractional-order differential equations. Appl. Math. Lett. 2022, 124, 107709. [Google Scholar] [CrossRef]
  14. 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]
  15. 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]
  16. Kuzmina, K.; Marchevsky, I. The boundary integral equation solution in vortex methods with the airfoil surface line discretization into curvilinear panels. In Proceedings of the Topical Problems of Fluid Mechanics, Prague, Czech Republic, 20–22 February 2019; pp. 131–138. [Google Scholar]
  17. Lienert, M.; Tumulka, R. A new class of Volterra-type integral equations from relativistic quantum physics. J. Integral Equ. Appl. 2019, 31, 535–569. [Google Scholar] [CrossRef]
  18. 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]
  19. 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]
  20. Matoog, R.T. Treatments of probability potential function for nuclear integral equation. J. Phys. Math. 2017, 8, 1000226. [Google Scholar]
Figure 1. Comparison between the exact solution and numerical solution in the linear case of the integral equation at T = 0.01.
Figure 1. Comparison between the exact solution and numerical solution in the linear case of the integral equation at T = 0.01.
Symmetry 16 01048 g001
Figure 2. Comparison between the exact solution and numerical solution in the nonlinear case of the integral equation at T = 0.01.
Figure 2. Comparison between the exact solution and numerical solution in the nonlinear case of the integral equation at T = 0.01.
Symmetry 16 01048 g002
Figure 3. Comparison between the exact solution and numerical solution in the linear case of the integral equation at T = 0.8.
Figure 3. Comparison between the exact solution and numerical solution in the linear case of the integral equation at T = 0.8.
Symmetry 16 01048 g003
Figure 4. Comparison between the exact solution and numerical solution in the nonlinear case of the integral equation at T = 0.8.
Figure 4. Comparison between the exact solution and numerical solution in the nonlinear case of the integral equation at T = 0.8.
Symmetry 16 01048 g004
Table 1. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.01 , in the linear case l = 1 , and nonlinear case.
Table 1. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.01 , in the linear case l = 1 , and nonlinear case.
u E x a c t U l i n e a r U n o n l i n e a r E e r r o r l i n e a r E r r o r n o n l i n e a r
−11.00011.0000273011.000156796 7.2699 × 10 5 5.6796 × 10 5
−0.50.25010.2505587940.2501511588 4.58794 × 10 4 5.511588 × 10 5
00.00010.000322780.000093462 2.2278 × 10 4 6.538 × 10 6
0.50.25010.250603190.2503023022 5.0319 × 10 4 2.023022 × 10 4
11.00011.000270341.0000273 1.7034 × 10 4 7.27 × 10 5
Table 2. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.4 , in the linear case l = 1 , and nonlinear case.
Table 2. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.4 , in the linear case l = 1 , and nonlinear case.
u E x a c t U l i n e a r U n o n l i n e a r E e r r o r l i n e a r E r r o r n o n l i n e a r
−11.161.2094647291.160215981 0.04946429 2.15981 × 10 4
−0.50.410.44195804770.4099824773 0.0319580477 1.75227 × 10 5
00.160.15535721230.1599678394 4.6427877 × 10 3 3.21606 × 10 5
0.50.410.39778821590.4099281713 0.0122117841 7.18287 × 10 5
11.161.2094647291.160085598 00.049464729 8.5598 × 10 5
Table 3. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.8 , in the linear case l = 1 , and nonlinear case.
Table 3. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.8 , in the linear case l = 1 , and nonlinear case.
u E x a c t U l i n e a r U n o n l i n e a r E e r r o r l i n e a r E r r o r n o n l i n e a r
−11.641.85102868851.640901375 0.211028688 5.4191 × 10 4
−0.50.890.83715834410.8896882523 0.0528416559 3.117477 × 10 4
00.640.61044123540.6396523326 0.0295587646 3.476674 × 10 4
0.50.890.8394536210.8896515856 0.050546379 3.484144 × 10 4
11.641.8510286861.640038297 0.211028686 3.8297 × 10 5
Table 4. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.01 , in the linear case l = 1 , and nonlinear case.
Table 4. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.01 , in the linear case l = 1 , and nonlinear case.
u E x a c t U l i n e a r U n o n l i n e a r E e r r o r l i n e a r E r r o r n o n l i n e a r
−11.00011.0006455521.000156796 5.45552 × 10 4 5.6796 × 10 4
−0.50.25010.24992990370.2508547062 7.547062 × 10 4 1.700963 × 10 4
00.00010.00060810.0001345009 5.081 × 10 4 3.45009 × 10 4
0.50.25010.2492580.2498976982 8.42 × 10 4 2.23018 × 10 4
11.00011.0006455521.0002703450 5.45552 × 10 4 1.70345 × 10 4
Table 5. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.4 in the linear case l = 1 , and nonlinear case.
Table 5. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.4 in the linear case l = 1 , and nonlinear case.
u E x a c t U l i n e a r U n o n l i n e a r E e r r o r l i n e a r E r r o r n o n l i n e a r
−11.161.1096579131.159812825 0.050342087 1.87175 × 10 4
−0.50.410.3944369430.4101825043 0.015563057 1.825043 × 10 4
00.160.15098163820.1601612092 9.0183618 × 10 3 1.612092 × 10 4
0.50.410.40032110.40981749 9.6789 × 10 3 1.8251 × 10 4
11.161.210342081.160056829 0.05034208 5.6829 × 10 5
Table 6. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.8 , in the linear case l = 1 , and nonlinear case.
Table 6. Discusses the numerical results and estimating errors at μ = 1 , λ = 1 , T = 0.8 , in the linear case l = 1 , and nonlinear case.
u E x a c t U l i n e a r U n o n l i n e a r E e r r o r l i n e a r E r r o r n o n l i n e a r
−11.641.852486331.64054191 0.21248633 5.4191 × 10 4
−0.50.890.782164420.890210195 0.10783558 2.10195 × 10 4
00.640.61044123540.6403476674 0.0295587646 3.476674 × 10 4
0.50.890.8394536210.88965158 0.050546379 3.4842 × 10 4
11.641.8510286861.640038297 0.211028686 3.8297 × 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

Alharbi, F.M. Stability Analysis of the Solution for the Mixed Integral Equation with Symmetric Kernel in Position and Time with Its Applications. Symmetry 2024, 16, 1048. https://doi.org/10.3390/sym16081048

AMA Style

Alharbi FM. Stability Analysis of the Solution for the Mixed Integral Equation with Symmetric Kernel in Position and Time with Its Applications. Symmetry. 2024; 16(8):1048. https://doi.org/10.3390/sym16081048

Chicago/Turabian Style

Alharbi, Faizah M. 2024. "Stability Analysis of the Solution for the Mixed Integral Equation with Symmetric Kernel in Position and Time with Its Applications" Symmetry 16, no. 8: 1048. https://doi.org/10.3390/sym16081048

APA Style

Alharbi, F. M. (2024). Stability Analysis of the Solution for the Mixed Integral Equation with Symmetric Kernel in Position and Time with Its Applications. Symmetry, 16(8), 1048. https://doi.org/10.3390/sym16081048

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop