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
Here, in Equation (1), the given functions belong to —space. The integration domain is a closed bounded set and depends on the vector of position, while the time . The kernel while The constant describes the type of the integral equation, and carries a physical significance. Finally, is the unknown function and will be discussed in the space of the given function .
Differentiating Equation (1) with respect to the variable t, we obtain
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,
Then, assume the assumptions:
(i) Position kernel check the offline condition, (M-constant).
(ii) Time kernel check status —constant).
(iii) The given function
belongs to
—space, and its norm is
(iv) The function
for the constants
,
, satisfies the following conditions:
Theorem 1. If the assumptions (i)–(iv) are verified, then Equation (1) has a single solution in , under the condition 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), transforms into itself.
Proof. Based on Equations (3) and (4), we obtain
By utilizing the assumptions (ii)–(iii) and subsequently applying the Cauchy Schwarz inequality, we obtain,
Considering the given constraints (i), (iv-a), we have,
Hence,
maps the ball
into itself where
the positive values of
, lead to
. This leads us to say that
Furthermore, W is bounded where
□
Lemma 2. Under (i), (ii), and (iv-b), is contraction operator in .
Proof. Regarding the two functions
then from Equations (3) and (4), we find
With the aid of (i), (ii), and (iv-b), and Cauchy–Schwarz inequality, we obtain
Hence, is continuous in . Moreover, under the condition is a contraction operator. Hence, 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
, to have
set,
to have
Using assumptions (ii)–(iv-a), to have
By utilizing the Cauchy–Schwarz inequality on the Hammerstein integral term, and considering the assumptions (i) and (iv-b), the inequality is simplified to
When
, the inequality (12), with condition (ii) takes the form
Since Equation (14), for
is obviously true then the sequence
converges, and
The series (15) exhibits uniform convergence due to the uniform behavior of its terms. are dominated by .
To discuss the error stability, we assume
is the numerical solution of Equation (1); hence, we have
Following the same previous way, we have
The inequality (17) proved that the error is stable and has a unique representation under the condition,
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
as
,
, then from Equation (1), we have (see [
2,
20])
Using Equation (18) in Equation (1), and neglecting
, we have
Here, we used the notations,
Equation (19) represents a set of n-dimensional HIEs, and its solution relies on the specified function the kernel p (u, w), and the known function
5. Special Cases
(i) As an important special case let, in (19),
to have
Equation (20) represents an integral equation of Hammerstein type and its solution depends on the kind of kernel.
(ii) Let
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)
(logarithmic kernel) hence, we have
Differentiating the above equation with respect to u, we obtain
Using the substation
, we have
The integral Equation (22) has appeared in both combined infrared gaseous radiation and molecular conduction, see [
20].
(iii) Three-dimensional integral equation
Let
, we have
here
.
If , 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 ,where and , are the stress and strain rate intensities, respectively, while , 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
approximately by a kernel
; that is
The set of functions {Bi(u)} and {Ci(w)} are linearly independent, such that
hence, the solution of Equation (19) associated to the kernel
takes the form
Using Equation (20) in Equation (22), we have
where
here,
are constants to be determined from the following formula
Equation (28) represents a system of NAEs that can be written as a matrix equation
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 in Equation (27) satisfy the following conditions,and,then, the NAS (29) has a single solution , and is the single solution of Equation (27) in -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 of Equation (28) maps -space into itself.
Proof. Let V be the set of functions
in
such that
From Equation (27), we have
Hence, after applying Cauchy–Schwarz inequality, and using condition (30), we have
as
, the last inequality yields
Hence,
is a bounded operator that maps the set U into itself, where
□
Lemma 4. Under the condition (31), is a contraction mapping in .
For the functions
and
in
, Equation (24) leads to
Introducing condition (31), then applying Cauchy–Schwarz inequality three times, the above inequality takes the form
As
, the previous inequality can be reduced to
Thus, is a continuous operator in the space . If , then is a contraction operator. Hence, has a unique fixed point , which is the unique solution of the algebraic system (27). In view of Theorem 2, the algebraic system (27) has a single solution in
7. Examples and Numerical Results
If we divide the interval [0, T] as
, and using the degenerate method, we obtain
where
The solution of Equation (35) gives the following results,
Here, is the exact solution of Equation (34).
Approximate the kernel of Equation (36)
in the form,
where
Using the degenerate kernel, we obtain
The solution of Equation (38), gives Here, 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,
The exact solution is .
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,
The exact solution is .
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.