Abstract
In this study, the Fredholm hypersingular integral equation of the first kind with a singular right-hand function on the interval is solved. The discontinuous solution on the domain is approximated by a piecewise polynomial, and a collocation method is introduced to evaluate the unknown coefficients. This method, which can be applied to both linear and nonlinear integral equations, is very simple and straightforward. The presented illustrations relate that the results are very accurate compared to the other methods in the literature.
1. Introduction
Hypersingular integral equations (HSIEs) are widely presented in different fields of physics and engineering such as crack problems in fracture mechanics [1], radiation problems involving thin-submerged plates [2–4], antenna theory [5], aerodynamics [6, 7], elliptic wing case in Prandtl’s equation [8, 9], electrodynamics, and quantum physics [6].
Many two-dimensional boundary value problems in mathematical physics can be formulated as HSIEs. It is investigated that how the problems such as two-dimensional flow post a rigid plate in an infinite fluid, acoustic scattering by a hard plate, water wave interaction with thin impermeable barriers, stress fields around cracks, certain specific problems involving cracks in an elastic medium and scattering of surface water waves by obstacles in the form of thin barriers, and some other issues give rise to HSIEs (see [7]). There are two strategies to solve these problems. The direct methods in which the problem is solved using the method of solving differential equations are finite difference methods [10] and Runge-Kutta-embedded methods [11]. Another robust tool is indirect methods in which the differential equation is transformed into an equivalent equation such as a HSIE. The numerical solutions for the second category is often less complicated.
In certain cases, there exist some analytical methods for solving HSIEs, but in general, there is no known analytical method. Thus, different numerical methods are introduced and developed for solving HSIEs [8, 9, 12–26].
The following first kind HSIE with various justifications on the boundaries was studied in [21] by using different types of Chebyshev polynomials. In Equation (1), the right-hand function is considered a smooth function, and the hypersingular integral is understood as Hadamard finite part as follows:
Recently, a special class of Equation (1) was introduced under the assumption that the right-hand function possesses singularities at some points of the interval [5, 6, 27]. In this case, the solution is not smooth in the integration domain. It was shown in [21] that if the right-hand function in Equation (1) belongs to a Hölder class, then, under some additional conditions on the boundaries, the solution can be represented as or , where is a smooth function. While, for the singular right-hand function , the class of solutions for Equation (1) is not determined yet.
Authors in [28, 29] introduced a variant of the generalized Euler-Maclaurin summation formula for the product integration rule and applied it for approximating weakly singular integrals. Also, with the expansion of the theory of fractional calculus in recent years, many attempts have been made to use fractional calculus to solve singular integral equations. These methods are suitable for solving weakly singular integral equations, but they cannot be used for a stronger singular integral equations such as HSIEs.
Recently, an iterative projection method is presented for solving linear and nonlinear HSIEs with non-Riemann integrable functions on the right-hand side [25]. They applied the continuous method and presented the convergence of the proposed method based on the Lyapunov stability theory. Not many papers have been published to solve the Equation (1) with singular right-hand function yet. To avoid the complexity of the method presented in [25] in the present study, we consider a simple collocation method for the numerical solution of HSIE (Equation (1)) with the singular right-hand function. The most important advantage of our method is its simplicity and adaptability to different states of the right-hand function. The rest of this paper is organized as follows: in Section 2, the piecewise polynomial collocation method is introduced and applied to the proposed linear and nonlinear hypersingular integral equation. The numerical illustrations are presented in Section 3, and finally, a brief conclusion appeared in Section 4.
2. Method of Solution
In this section, we consider a new collocation method for solving the following linear and nonlinear HSIE of the first kind: where the right-hand function has different types of singularities at the points of the interval . The singularities of may have the following structure: or any other singular type. In Equation (4), and , are considered smooth functions. To introduce the new method, we need the following definitions.
Definition 1. We define the Cauchy and Hadamard integrals, respectively, as follows
Theorem 2. The Cauchy and Hadamard integrals defined in Equations (5) and (6) are calculated as follows:
Proof. By the definition of Cauchy integral (Equation (5)), we have Also, by the definition of Hadamard integral (Equation (6)), we have
Now, we intend to solve Equation (3) by representing the solution as a piecewise continuous polynomial. To do this, we set and , and then, we seek a solution for Equation (3) on the interval . According to the structure of the function on the right hand, the solution may be discontinuous at the points , . Therefore, we consider the solution in the following form: where is a given integer.
2.1. Application to Linear HSIE
Consider Equation (3) with . We rewrite Equation (3) as follows:
Substituting Equation (11) in Equation (3), we get which gives
We select distinct points on the th subinterval , .
Now, we collocate (14) as follows:
Equation (15) represents a system of linear equations of unknowns , , , which can be converted to matrix format as where and is a square matrix. The elements of are as follows:
Remark 3. Due to Definition 1, the Hadamard integrals , are linearly independent functions. Also, the points are distinct and . Then, from Haar property for independent functions, we conclude that . This means that the linear system (Equation (16)) has a unique solution.
2.2. Application to Nonlinear HSIE
Let us consider Equation (3) with . For solving the nonlinear HSIE (Equation (3)), we rewrite it as follows:
Substituting Equation (11) in Equation (19), we get after some simplification we get which gives
We note that in Equation (21) is the corresponding coefficient of in the expansion of , which is calculated using Adomian polynomials as follows:
If we collocate Equation (22) at the points , , we get the following nonlinear system of equations:
3. Numerical Illustrations
In this section, we present the effectiveness of the proposed method by solving some HSIEs. All the simulations have been done using Maple software with an accuracy of 20 digits. The first few are as follows (see Equation (8)).
Example 4. Suppose the following linear HSIE of the first kind [25], where
Here, we set , , and , where has singularity. Then, we apply the proposed method with and get which gives the exact solution. In Figure 1, we plotted the graph of exact and approximate solutions for and .
Example 5. Consider the following linear HSIE of the first kind, where
Here, we set , , , and . Then, we apply the proposed method with and get which is the exact solution to the problem. The graph for the exact and approximate solution for and is plotted in Figure 2.
Example 6. Suppose the following nonlinear HSIE of the first kind [25], where
The right-hand function is singular at three points , , and . Then, we apply the proposed method with , and we get which is the exact solution to the problem. The corresponding graphs for the exact and approximate solution are plotted in Figure 3.
Example 7. For the next example, we suppose the following nonlinear HSIE of the first kind, where
The right-hand function is singular at three points , , and , while the exact solution is not known. The graph for approximate solutions for different values of is plotted in Figure 4.
Example 8. Suppose the following nonlinear HSIE of the first kind, where the right-hand function is singular at . We solved this problem with different values of . The graphs for approximate solutions are plotted in Figure 5. Note that the exact solution is not known here.
Remark 9. According to the illustrated examples, we can see that when the solution of the HSIE is polynomial, the present method produces the exact solution of the problem. In case the exact solution is not available, the numerical approximations converge to each other by increasing (see Examples 7 and 8). Comparison of numerical results with the results of [25] show that the present method is very efficient.
4. Conclusion
In this study, we have developed a novel collocation method for solving the first kind of hypersingular integral equations of the first kind with a singular right-hand function. The solution is considered as a piecewise polynomial on the integration domain, and the unknown coefficients are determined by solving an algebraic system of equations. The new method applies to both linear and nonlinear integral equations. This method covers a wide range of singular integral equations with different types of singularities on the right-hand function. It can also be applied for supersingular integral equations with singular right-hand function.
For further research in this field, it is suggested that the class of the solution space of the equation be analyzed analytically based on the type of singularity of the right-hand function.
Data Availability
No underlying data was collected or produced in this study.
Conflicts of Interest
The authors declare that they have no conflict of interest.