Next Article in Journal
Mindfulness Model Using Polariton Oscillation in Plasmonic Circuit for Human Performance Management
Next Article in Special Issue
Eigenfunction Families and Solution Bounds for Multiplicatively Advanced Differential Equations
Previous Article in Journal / Special Issue
On the Periodicity of General Class of Difference Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Approximate Methods for Solving Linear and Nonlinear Hypersingular Integral Equations

1
Department of Mathematics, The Penza State University, 40, Krasnaya Str., 440026 Penza, Russia
2
Department of Computational Physics, Saint Petersburg State University, 7/9 Universitetskaya Emb., 199034 Saint Petersburg, Russia
*
Author to whom correspondence should be addressed.
Axioms 2020, 9(3), 74; https://doi.org/10.3390/axioms9030074
Submission received: 2 June 2020 / Revised: 24 June 2020 / Accepted: 25 June 2020 / Published: 1 July 2020

Abstract

:
We propose an iterative projection method for solving linear and nonlinear hypersingular integral equations with non-Riemann integrable functions on the right-hand sides. We investigate hypersingular integral equations with second order singularities. Today, hypersingular integral equations of this type are widely used in physics and technology. The convergence of the proposed method is based on the Lyapunov stability theory of solutions of ordinary differential equation systems. The advantage of the method for linear equations is in simplicity of unique solvability verification for the approximate equations system in terms of the operator logarithmic norm. This makes it possible to estimate the norm of the inverse matrix for an approximating system. The advantage of the method for nonlinear equations is that neither the existence or reversibility of the nonlinear operator derivative is required. Examples are given illustrating the effectiveness of the proposed method.

1. Introduction

The importance of developing analytical and numerical methods for solving hypersingular integral equations is determined by a variety of fields of mathematics and by applications that use hypersingular integral equations.
Hadamard introduced the concept of a finite part of an integral, or the hypersingular integral in modern terminology, when studying hyperbolic equations. The Riemann boundary problem leads to hypersingular integral equations in exceptional cases. The boundary integral equations method reduces the dimensions of partial differential equations; that leads to hypersingular integral equations.
Hypersingular integral equations, singular integral equations and Riemann boundary problem are widely used in aerodynamics, electrodynamics, quantum physics, antennae theory and many other fields of physics and engineering [1,2,3,4,5].
Analytical methods for solving singular and hypersingular integral equations are known only for certain particular types of equations [6,7,8]. Thus, the importance of constructing numerical solutions is clear.
Developing numerical methods for solving singular integral equations began in the middle of the last century. By now, exhaustive results have been obtained for many types of equations. A detailed account of numerical methods for solving singular integral equations as well as numerous bibliography references can be found in [9,10,11,12,13,14].
Numerical methods for solving hypersingular integral equations have been developed to a much lesser extent. Mostly numerical methods to solve hypersingular integral equations of the first kind have been developed. Numerical methods for solving hypersingular integral equations of the second kind have been much less developed. Apparently, hypersingular integral equations of the first kind are more common. Naturally, the equations of the first kind are widely used in aerodynamics (one-dimensional [15] and multi-dimensional [5,16] Prandtl equation), electrodynamics, antennae theory, etc.
The following methods are used in solving hypersingular integral equations of the first kind.
Collocations, mechanical quadratures and Galerkin methods were employed to solve equations with p = 2 singularity [6,17,18,19].
Approximate methods for solving hypersingular integral equations having singularities of order p = 2 , 3 , , and defined on closed smooth integration contours are constructed in [20].
In [21,22] spline-collocation methods for solving hypersingular and polyhypersingular integral equations of the second kind with odd and even singularities have been developed and justified. The spline-collocation methods for solving nonlinear hypersingular and polyhypersingular integral equations have been developed and justified in [23].
An iterative projection method for solving linear and nonlinear hypersingular integral equations, and polyhyperpersingular and multidimensional hypersingular equations, was proposed in [24].
In [22] the unique solvability of hypersingular integral equations with even singularities (p = 2, 4, …) was proven. Meanwhile the convergence of approximate solution to the exact one was not justified. In [24] a unique solvability of the spline-collocation method was proven. In addition, for hypersingular integral equations with bounded right-hand sides the convergence of an approximate solution sequence to the exact solution was proven under certain additional conditions.
The iterative projection method proposed here overcomes these limitations. It was shown that if the exact equation has a solution for large enough N, where N is the dimension of an approximate system of equations, an approximate solution converges to the exact one.
Hypersingular integral equations with bounded right-hand sides are a small subset of the hypersingular integral equations family. Therefore, the problem arises of constructing and justifying approximate methods for solutions for hypersingular integral equations with non-Riemann integrable functions on the right-hand sides. This paper is devoted to those issues.
A large number of works are devoted to approximate methods for solving hypersingular integral equations of the first kind
1 1 x ( τ ) d τ ( τ t ) 2 + 1 1 h ( t , τ ) d τ = f ( t ) .
To solve the Equation (1), collocation and mechanical quadrature methods [17,18], the method of orthogonal polynomials [25], the method of discrete vortices [19], the method of homotopy [26] and others are used.
In the works [27,28,29] computational schemes for the approximate solution of the Equation (1) are constructed and their justification is carried out under the assumption that the solution has the forms x ( t ) = ( 1 t 2 ) ± 1 / 2 ω ( t ) or x ( t ) = ( ( 1 t ) / ( 1 + t ) ) ± 1 / 2 ω ( t ) , where ω ( t ) is a smooth function.
The hypersingular integral equations
1 π 1 1 x ( τ ) d τ ( τ t ) 2 = f ( t ) , 1 < t < 1 ,
are widely used in aerodinamical problems and in the theory of antennae [30,31]. In the works [30,31] the Equation (2) is investigated under the assumption that the right-hand side has the form f ( t ) = 1 / ( t c ) or f ( t ) = δ ( t c ) , where δ ( t ) is the delta-function. An analytical solution of the Equation (2) with the indicated right-hand sides is obtained under the assumption that it has the form x ( t ) = 1 t 2 φ ( t ) .
A fairly detailed review of analytical and numerical methods for solving hypersingular integral equations is given in [32].
In this paper, we propose an approach to solving linear and nonlinear hypersingular integral equations, the right parts of which contain functions with power features.
In particular, the right-hand sides of the form
f ( t ) = g ( t ) 1 t c 1 1 t c 2 1 t c l , l = 1 , 2 , , 1 < c 1 < < c l < 1 ,
are considered. Here g ( t ) is a smooth function.
Below, for simplicity of notation, we put l = 1 in (3).
Remark 1.
It can be shown that if in the hypersingular integral Equation (1) of the first kind the right side f ( t ) H , H is a Holder class, then the solution to this equation has the form x ( t ) = ( 1 t 2 ) ± 1 / 2 or x ( t ) = ( ( 1 + t ) / ( 1 t ) ) ± 1 / 2 . For singular right-hand sides, the classes of solutions of (1) are unknown.
Below, when constructing and justifying the computational method, we assume that the Equation (1) with a given right-hand side has a unique solution.
The proposed method has the following advantages:
(1)
It allows us to extend collocations and mechanical quadratures methods to hypersingular integral equations with non-Riemann integrable right sides;
(2)
For linear hypersingular integral equations, it allows one to verify the inverse operator existence and estimate its norm quite easily;
(3)
The method is stable with respect to the operator and right hand side perturbations;
(4)
The method does not require the existence and reversibility of the nonlinear operator derivative.
The paper is organized as follows. The continuous method for linear and nonlinear operator equations is explained in Section 2. The numerical method for solving hypersingular integral equations is presented in Section 2.

2. Continuous Method and Its Convergence Properties

The method we employ in the next section for solving hypersingular integral equations is based on the continuous method introduced in [33].

Continuous Method for Solving Operator Equations

The continuous method for solving operator equations is based on the Lyapunov theory of stability.
Let x ( t ) be a solution of the differential equation in a Banach space B
d x d t = F ( t , x )
which is defined for all t t 0 . The solution x ( t ) is said to be stable if (i) for each ε > 0 there is a corresponding δ = δ ( ε ) > 0 such that any solution x ˜ ( t ) of (4) which satisfies the inequality | x ˜ ( t 0 ) x ( t 0 ) | < δ exists and satisfies the inequality | x ˜ ( t ) x ( t ) | < ε for all t t 0 .
It is said to be asymptotically stable if in addition (ii) | x ˜ ( t ) x ( t ) | 0 if t whenever | x ˜ ( t 0 ) x ( t 0 ) | is sufficiently small.
We will use the following notation:
B ( a , r ) = { z B : z a r } , S ( a , r ) = { z B : z a = r } , R e ( K ) = ( K ) = ( K + K ) / 2 Λ ( K ) = lim h 0 ( I + h K 1 ) h 1
.
Here B is a Banach space, a B , K is a linear operator on B, Λ ( K ) is the logarithmic norm [34] of the operator K, K is the conjugate operator to K and I is the identity operator.
The analytical expressions for logarithmic norms are known for operators in many spaces. We restrict ourselves to a description of the three norms.
Let A = { a i j } , i , j = 1 , 2 , , n , be a matrix.
In the n-dimensional space R n of vectors x = ( x 1 , , x n ) the following norms are often used:
octahedral - x 1 = i = 1 n | x i | ; cubic - x 2 = max 1 i n | x i | ; spherical   ( Euclidean ) - x 3 = ( i = 1 n x i 2 ) 1 / 2 .
Here are analytical expressions of the logarithmic norm of n × n matrix A = ( a i j ) , due to the above norms of the vectors:
octahedral logarithmic norm Λ 1
Λ 1 ( A ) = max 1 j n a j j + i j | a i j | ;
cubic logarithmic norm Λ 2
Λ 2 ( A ) = max 1 i n a i i + j i | a i j | ;
spherical (Euclidean) logarithmic norm Λ 3
Λ 3 ( A ) = λ max A + A 2 ,
where A is the conjugate matrix for A .
Note that the logarithmic norm of the same matrix can be positive in one space and negative in another.
The logarithmic norm has the some properties which are very useful for numerical mathematics.
Let A , B be n × n matrices with complex elements; and x = ( x 1 , , x n ) , y = ( y 1 , , y n ) , ξ = ( ξ 1 , , ξ n ) and η = ( η 1 , , η n ) are n-dimensional vectors with complex components. Let the systems of algebraic equations A x = ξ and B y = η be given. The norm of a vector and its subordinate operator norm of the matrix are agreed upon; the logarithmic norm Λ ( A ) corresponds to the operator norm.
Theorem 1
([35]). If Λ ( A ) < 0 , the matrix A is non-singular and A 1 1 / | Λ ( A ) | .
Theorem 2
([35]). Let A x = ξ , B y = η and Λ ( A ) < 0 , Λ ( B ) < 0 . Then
x y ξ η | Λ ( B ) | + A B | Λ ( A ) Λ ( B ) | .
Some properties of the logarithmic norm in a Banach space, which are useful in numerical mathematics, are given in [34].
Let us consider in a Banach space B, the Cauchy problem
d x ( t ) d t = A ( x ( t ) ) ,
x ( 0 ) = x 0 .
Let us assume that the nonlinear operator A has a Frechet derivative and A ( 0 ) = 0 .
The sufficiently satisfying conditions of asymptotically stability for the solution of the Cauchy problem (5), (6) were obtained in [36,37].
Theorem 3.
Let the integral 0 t Λ ( A ( φ ( τ ) ) ) d τ be non-positive (respectively, be negative and satisfy lim t + 1 t 0 t Λ ( A ( φ ( τ ) ) ) d τ α φ , α φ > 0 ) for any differentiable curve φ ( t ) lying in a ball B ( 0 , r ) of some radius r. Then the trivial solution of Equation (5) is stable (respectively, asymptotically stable).
Remark 2.
Additionally, the Theorem is valid under r = .
Let us consider in a Banach space B a nonlinear operator equation
A ( x ) f = 0 ,
where operator A acts from B into B.
We associate Equation (7) with the Cauchy problem
d x ( t ) d t = A ( x ( t ) ) f ,
x ( 0 ) = x 0 .
Let x be a the solution of Equation (7). Let us make the change of variable x = x + v in Equation (8). This change reduces the Cauchy problem (8), (9) to the form
d v ( t ) d t = A ( x + v ( t ) ) A ( x ) ,
v ( 0 ) = x 0 x .
It is easy to see that if the trivial solution of Equation (10) is globally asymptotically stable, then lim t + v ( t ) 0 . So, for any initial value the solution of the Cauchy problem (8), (9) tends to x . It follows from the next assertions which were proven in [33].
Theorem 4.
Let E q u a t i o n ( 7 ) have a solution x , and let inequality
lim t + 1 t 0 t Λ ( A ( g ( τ ) ) ) d τ α g , α g > 0 ,
be true on each differentiable curve g ( t ) lying in the Banach space B. Then the solution of the Cauchy problem (8), (9) converges to the solution x of E q u a t i o n ( 7 ) for any initial value.
Theorem 5.
Let Equation (7) have a solution x , and let the following conditions be satisfied on any differentiable curve g ( t ) lying in the ball B ( x , r ) .
1.
The inequality
0 t Λ ( A ( g ( τ ) ) ) d τ 0
holds for all t ( t > 0 ) .
2.
Inequality (12) is satisfied.
Then the solution of the Cauchy problem (8), (9) converges to the solution x of Equation (7).
Remark 3.
The sufficient condition for convergence of the Cauchy problem (8), (9) solution to the solution of the operator Equation (7) is given above. It was obtained by analysing Lyapunov stability. One of the first basic results in accretive operator theory was a relation between the solution of operator equation A u = 0 , where A is a locally Lipschitzian and accretive operator, and the differential equation d u d t = A u was obtained in [38].
Later, accretive operator theory and its applications for finding fixed points and constructing iterative procedures were studied by many authors. Basic results and a detailed bibliography devoted to the subject may be found in [39,40,41,42].

3. An Solution of Hypersingular Integral Equations with the Continuous Method

Let us consider the method of mechanical quadrature for solving hypersingular integral equation of the types
a ( t ) x ( t ) + 1 1 h ( t , τ ) x ( τ ) d τ ( τ t ) 2 = f ( t ) .
and
a ( t ) x ( t ) + 1 1 h ( t , τ , x ( τ ) ) d τ ( τ t ) 2 = f ( t ) .
It is assumed that in the Equations (13) and (14) the right-hand sides have features of the following types
f ( t ) = i = 1 l g i ( t ) 1 t c i , f ( t ) = g ( t ) i = 1 l 1 t c i ,
where 1 < c i < 1 , i = 1 , 2 , , l , l = 1 , 2 , ; g ( t ) , g i ( t ) , i = 1 , 2 , , l , —are continuous functions.
In what follows, without loss of generality, we set l = 1 .
Let us recall the Hadamard definition of hypersingular integrals [43].
Definition 1
([43]). The integral of the type
a b A ( x ) d x ( b x ) p + α
for an integer p and 0 < α < 1 , is defined as
lim x b a x A ( t ) d t ( b t ) p + α + B ( x ) ( b x ) p + α 1 ,
if A ( x ) has p derivatives in the neighborhood of the point b. Here B ( x ) is any function that satisfies the following two conditions:
(i) 
The above limit exists;
(ii) 
B ( x ) has at least p derivatives in the neighborhood of the point x = b .
It is easy to see [43], that the conditions (i) and (ii) are sufficient for the existence of the limit.
Chikin in [44] introduced the following definition.
Definition 2
([44]). The Cauchy–Hadamard principal value of the integral
a b φ ( τ ) d τ ( τ c ) p , a < c < b ,
is defined as
a b φ ( τ ) d τ ( τ c ) p = lim v 0 a c v φ ( τ ) d τ ( τ c ) p + c + v b φ ( τ ) d τ ( τ c ) p + ξ ( v ) v p 1 ,
where ξ ( v ) is a function constructed so that the limit exists.

3.1. An Approximate Solution of Linear Hypersingular Integral Equations with Second Order Singularity

Consider a one-dimensional hypersingular integral equation of the type
K x a ( t ) x ( t ) + 1 1 h ( t , τ ) x ( τ ) d τ ( τ t ) 2 = f ( t ) ,
where f ( t ) = g ( t ) / ( t c ) or f ( t ) = g ( t ) / ( ( 1 t 2 ) ( t c ) ) , 1 < c < 1 , g ( t ) C [ 1 , 1 ] .
Divide the interval [ 1 , 1 ] into two subintervals [ 1 , c ] , [ c , 1 ] .
Let us fix a positive integer N 0 . Put h = 2 / N 0 , N 1 = ( 1 + c ) / h , N 2 = ( 1 c ) / h , N = N 1 + N 2 .
Divide the interval [ 1 , c ] into N 1 subintervals at the points t k = 1 + ( c + 1 ) k / N 1 , k = 0 , 1 , , N 1 .
Divide the interval [ c , 1 ] into N 2 subintervals at the points τ k = c + ( 1 c ) k / N 2 , k = 0 , 1 , , N 2 .
Let us introduce the nodes t ¯ 0 = t 0 + 1 / 2 ( N 1 ) 2 , t ¯ k = t k , k = 1 , 2 , , N 1 1 , t ¯ N 1 = t N 1 1 / 2 ( N 1 ) 2 ; τ ¯ 0 = τ 0 + 1 / 2 ( N 2 ) 2 , τ ¯ k = τ k , k = 1 , 2 , , N 2 1 , τ ¯ N 2 = 1 1 / 2 ( N 2 ) 2 .
As an approximate solution of (16), we shall seek in the form of a continuous function
x N ( t ) = k = 0 N 1 α k φ k ( t ) + k = 0 N 2 β k ψ k ( t ) ,
where φ k ( t ) , k = 0 , 1 , , N 1 , ψ k ( t ) , k = 0 , 1 , , N 2 is a family of basis functions.
For nodes t k , k = 1 , , N 1 1 , the corresponding basis elements are determined by
φ k ( t ) = 0 , t k 1 t t k 1 + 1 N 1 2 , N 1 2 ( 1 + c ) N 1 2 ( t t k 1 ) 1 ( 1 + c ) N 1 2 , t k 1 + 1 N 1 2 t t k 1 N 1 2 , 1 , t k 1 N 1 2 t t k + 1 N 1 2 , N 1 2 ( 1 + c ) N 1 2 ( t t k 1 N 1 2 ) + 1 , t k + 1 N 1 2 t t k + 1 1 N 1 2 , 0 , t k + 1 1 N 1 2 t t k + 1 , 0 , t [ 1 , 1 ] \ [ t k 1 , t k + 1 ] .
For boundary nodes t k , k = 0 and k = N 1 the corresponding basis elements are defined as
φ 0 ( t ) = 1 , 1 t 1 + 1 N 1 2 , N 1 2 ( 1 + c ) N 1 2 ( t + 1 1 N 1 2 ) + 1 , 1 + 1 N 1 2 t t 1 1 N 1 2 , 0 , t 1 1 N 1 2 t t 1 , 0 , [ 1 , 1 ] \ [ t 0 , t 1 ] ;
and
φ N 1 ( t ) = 0 , 1 t t N 1 1 + 1 N 1 2 , N 1 2 ( 1 + c ) N 1 2 ( t t N 1 1 ) 1 ( 1 + c ) N 1 2 , t N 1 1 + 1 N 1 2 t c 1 N 1 2 , 1 , c 1 N 1 2 t c .
For nodes τ k , k = 0 , 1 , , N 2 , the corresponding basis elements ψ k , k = 0 , 1 , , N 2 , are determined in a the similar way: For nodes τ k , k = 1 , , N 2 1 , the corresponding basis elements are determined by
ψ k ( t ) = 0 , τ k 1 t τ k 1 + 1 N 2 2 , N 2 2 ( 1 c ) N 2 2 ( t τ k 1 ) 1 ( 1 c ) N 2 2 , τ k 1 + 1 N 2 2 t τ k 1 N 2 2 , 1 , τ k 1 N 2 2 t τ k + 1 N 2 2 , N 2 2 ( 1 c ) N 2 2 ( t τ k 1 N 2 2 ) + 1 , τ k + 1 N 2 2 t τ k + 1 1 N 2 2 , 0 , τ k + 1 1 N 2 2 t τ k + 1 , 0 , t [ 1 , 1 ] \ [ τ k 1 , τ k + 1 ] .
For boundary nodes τ k , k = 0 and k = N 2 the corresponding basis elements are defined as
ψ 0 ( t ) = 1 , c t c + 1 N 1 2 , N 2 2 ( 1 c ) N 2 2 ( t c 1 N 2 2 ) + 1 , c + 1 N 2 2 t τ 1 1 N 2 2 , 0 , τ 1 1 N 2 2 t τ 1 , 0 , [ 1 , 1 ] \ [ c , τ 1 ] ;
and
ψ N 2 ( t ) = 0 , 1 t τ N 2 1 + 1 N 2 2 , N 2 2 ( 1 c ) N 2 2 ( t τ N 2 1 ) 1 ( 1 c ) N 2 2 , τ N 2 1 + 1 N 2 2 t 1 1 N 2 2 , 1 , 1 1 N 2 2 t 1 .
To simplify the description of computational scheme, we introduce the following notation:
(1)
Unite the nodes t k , k = 0 , 1 , , N 1 and τ l , l = 0 , 1 , , N 2 , denoting them by v i , i = 0 , 1 , , N , N = N 1 + N 2 ;
(2)
Unite the nodes t ¯ k , k = 0 , 1 , , N 1 and τ ¯ l , l = 0 , 1 , , N 2 , denoting them by v ¯ i , i = 0 , 1 , , N + 1 ;
(3)
Denote the family of basis functions { φ k } , k = 0 , 1 , , N 1 , { ψ l } , l = 0 , 1 , , N 2 by { ζ j } , j = 0 , 1 , , N + 1 ;
(4)
Denote by { γ k } , k = 0 , 1 , , N + 1 , unknowns { α i } , i = 0 , 1 , , N 1 , { β j } , j = 0 , 1 , , N 2 .
Here v i = t i , i = 0 , 1 , , N 1 , v N 1 + i = τ i , i = 1 , 2 , , N 2 ,
γ i = α i , i = 0 , 1 , , N 1 , γ N 1 + 1 + i = β i , i = 0 , 1 , , N 2 ,
ζ i = φ i , i = 0 , 1 , , N 1 , ζ N 1 + 1 + i = ψ i , i = 0 , 1 , , N 2 .
Let us recall that the points t N 1 and τ 0 coincide.
Applying the collocation method on the knots v ¯ k , k = 0 , 1 , , N + 1 to the Equation (16), we obtain the following system of algebraic equations for finding unknown coefficients { γ k } of the polygon (17)
a ( v ¯ k ) γ k + l = 0 N + 1 h ( v ¯ k , v l ) γ l 1 1 ζ l ( τ ) ( τ v ¯ k ) 2 d τ = f ( v ¯ k ) ,
k = 0 , 1 , , N + 1 .
Using the definition of hypersingular integrals, we receive:
v k 1 v k + 1 ζ k ( τ ) d τ ( τ v ¯ k ) 2 = 2 N 1 2 ( 1 + c ) N 1 2 ln ( ( 1 + c ) N 1 1 ) , k = 1 , 2 , , N 1 1 ;
v k 1 v k + 1 ζ k + 1 ( τ ) d τ ( τ v ¯ k ) 2 = 2 N 2 2 ( 1 c ) N 2 2 ln ( ( 1 c ) N 2 1 ) , k = N 1 + 2 , , N 1 ;
1 v 1 ζ 0 ( τ ) d τ ( τ + 1 1 2 N 1 2 ) 2 = 2 N 1 2 N 1 2 ln ( 2 ( c + 1 ) N 1 3 ) ( c + 1 ) N 1 2 ,
v N 1 1 v N 1 ζ N 1 ( τ ) d τ ( τ v N 1 + 1 2 N 1 2 ) 2 = 2 N 1 2 N 1 2 ln ( 2 ( c + 1 ) N 1 3 ) ( c + 1 ) N 1 2 ,
v N 1 v N 1 + 1 ζ N 1 + 1 ( τ ) d τ ( τ v N 1 1 2 N 2 2 ) 2 = 2 N 2 2 N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ,
v N 1 1 ζ N + 1 ( τ ) d τ ( τ 1 + 1 2 N 2 ) 2 = 2 N 2 2 N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ,
1 1 l = 1 N + 1 ζ l ( τ ) d τ ( τ + 1 1 2 N 1 2 ) 2 = 2 N 1 2 4 N 1 2 1 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ,
1 1 l = 0 N ζ l ( τ ) d τ ( τ 1 + 1 2 N 2 2 ) 2 = 2 N 2 2 4 N 2 2 1 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ,
1 1 l = 0 N + 1 ζ l ( τ ) d τ ( τ ( c 1 2 N 1 2 ) ) 2 = 2 N 1 2 4 N 1 2 1 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ,
1 1 l = 0 N + 1 ζ l ( τ ) d τ ( τ ( c + 1 2 N 2 2 ) ) 2 = 2 N 2 2 4 N 2 2 1 + N 2 2 ln ( 2 ( 1 + c ) N 2 3 ) ( 1 + c ) N 2 2 ,
1 1 l = 0 N + 1 φ l ( τ ) d τ ( τ v k ) 2 = N 1 ( c + 1 ) k N 1 2 N 1 ( c + 1 ) k + 2 N 1 2 ( 1 + c ) N 1 2 ln ( ( c + 1 ) N 1 1 ) ;
1 1 l = 0 N + 1 φ l ( τ ) d τ ( τ v k ) 2 = N 2 ( 1 c ) k N 2 2 N 2 ( 1 c ) k + 2 N 2 2 ( 1 c ) N 2 2 ln ( ( 1 c ) N 2 1 ) .
Here l , l , l , l indicates a summation over l N 1 , l N 1 + 1 , l k ( 1 k N 1 1 ) , l k ( N 1 + 2 k N 1 ) , respectively. Detailed calculations are given in [23].
We can rewrite the system (24) as
a ( v ¯ k ) γ k h ( v ¯ k , v ¯ k ) 2 N 1 2 ln ( N 1 1 ) ( 1 + c ) N 1 2 γ k + l = 0 N + 1 γ l h ( v ¯ k , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ k ) 2 = f ( v ¯ k ) , k = 1 , , N 1 1 ;
a ( v ¯ k ) γ k h ( v ¯ k , v ¯ k ) 2 N 2 2 ln ( N 2 1 ) ( 1 c ) N 2 2 γ k + l = 0 N + 1 γ l h ( v ¯ k , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ k ) 2 = f ( v ¯ k ) , k = N 1 + 2 , , N ;
a ( v ¯ 0 ) γ 0 h ( v ¯ 0 , v ¯ 0 ) ( 2 N 1 2 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ) γ 0 + l = 1 N + 1 γ l h ( v ¯ k , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ 0 ) 2 = f ( v ¯ 0 ) ;
a ( v ¯ N 1 ) γ N 1 h ( v ¯ N 1 , v ¯ N 1 ) ( 2 N 1 2 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ) γ N 1 + l = 0 N + 1 γ l h ( v ¯ N 1 , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 ) 2 = f ( v ¯ N 1 ) ;
a ( v ¯ N 1 + 1 ) γ N 1 + 1 h ( v ¯ N 1 + 1 , v ¯ N 1 + 1 ) ( 2 N 2 2 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ) γ N 1 + 1 + l = 0 N + 1 γ l h ( v ¯ N 1 + 1 , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 + 1 ) 2 = f ( v ¯ N 1 ) ;
a ( v ¯ N + 1 ) γ N + 1 h ( v ¯ N + 1 , v ¯ N + 1 ) ( 2 N 2 2 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ) γ N + 1 + l = 0 N γ l h ( v ¯ N + 1 , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 ) 2 = f ( v ¯ N + 1 ) .
Here , , indicates a summation over l k , l N 1 , l N 1 + 1 , respectively.
The system (37)–(42) is equivalent to the system
( s g n h ( v k , v k ) ) a ( v ¯ k ) γ k h ( v ¯ k , v ¯ k ) 2 N 1 2 ln ( N 1 1 ) ( 1 + c ) N 1 2 γ k + l = 0 N + 1 γ l h ( v ¯ k , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ k ) 2 = ( s g n h ( t k , t k ) ) f ( v k ) , k = 1 , , N 1 1 ;
( s g n h ( v k , v k ) ) a ( v ¯ k ) γ k h ( v ¯ k , v ¯ k ) 2 N 2 2 ln ( N 2 1 ) ( 1 c ) N 2 2 γ k + l = 0 N + 1 γ l h ( v ¯ k , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ k ) 2 = ( s g n h ( t k , t k ) ) f ( v ¯ k ) , k = N 1 + 2 , , N ;
( s g n h ( v 0 , v 0 ) ) a ( v ¯ 0 ) γ 0 h ( v ¯ 0 , v ¯ 0 ) ( 2 N 1 2 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ) γ 0 + l = 1 N + 1 γ l h ( v ¯ k , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ 0 ) 2 = ( s g n h ( v 0 , v 0 ) ) f ( v ¯ 0 ) ;
( s g n h ( v N 1 , v N 1 ) ) a ( v ¯ N 1 ) γ N 1 h ( v ¯ N 1 , v ¯ N 1 ) ( 2 N 1 2 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ) γ N 1 + l = 0 N + 1 γ l h ( v ¯ N 1 , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 ) 2 = ( s g n h ( v N 1 , v N 1 ) ) f ( v ¯ N 1 ) ;
( s g n h ( v N 1 + 1 , v N 1 + 1 ) ) a ( v ¯ N 1 + 1 ) γ N 1 + 1 h ( v ¯ N 1 + 1 , v ¯ N 1 + 1 ) ( 2 N 2 2 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ) γ N 1 + 1 + l = 0 N + 1 γ l h ( v ¯ N 1 , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 + 1 ) 2 = ( s g n h ( v N 1 , v N 1 ) ) f ( v ¯ N 1 ) ;
( s g n h ( v N + 1 , v N + 1 ) ) a ( v ¯ N + 1 ) γ N + 1 h ( v ¯ N + 1 , v ¯ N + 1 ) ( 2 N 2 2 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ) γ N + 1 + l = 0 N γ l h ( v ¯ N 1 , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ N + 1 ) 2 = ( s g n h ( v N + 1 , v N + 1 ) ) f ( v ¯ N + 1 ) .
Here , , indicates a summation over l k , l N 1 , l N 1 + 1 , respectively.
Let us write the system (43)–(48) in a matrix form
D X = F ,
where D = { d k l } , k , l = 0 , 1 , , N + 1 , X = ( x 0 , x 1 , , x N + 1 ) , F = ( f 0 , f 1 , , f N + 1 ) . The values { d k l } , { x k } , and { f k } are obvious.
The diagonal elements in the left-hand side of the system of Equations (43)–(48) have the following forms
d k k = s g n h ( v ¯ k , v ¯ k ) a ( v ¯ k ) h ( v ¯ k , v ¯ k ) 2 N 1 2 ln ( N 1 1 ) ( 1 + c ) N 1 2
k = 1 , 2 , , N 1 1 ,
d k k = s g n h ( v ¯ k , v ¯ k ) a ( v ¯ k ) h ( v ¯ k , v ¯ k ) 2 N 1 2 ln ( N 1 1 ) ( 1 + c ) N 1 2 ,
k = N 1 + 2 , , N ,
d 00 = s g n h ( v ¯ 0 , v ¯ 0 ) a ( v ¯ 0 ) h ( v ¯ 0 , v ¯ 0 ) ( 2 N 1 2 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ,
d N 1 , N 1 = ( s g n h ( v N 1 , v N 1 ) ) a ( v ¯ N 1 ) h ( v ¯ N 1 , v ¯ N 1 ) ( 2 N 1 2 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ) ,
d N 1 + 1 , N 1 + 1 = ( s g n h ( v N 1 + 1 , v N 1 + 1 ) ) a ( v ¯ N 1 + 1 ) h ( v ¯ N 1 + 1 , v ¯ N 1 + 1 ) ( 2 N 2 2 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ) ,
d N + 1 , N + 1 = ( s g n h ( v N + 1 , v N + 1 ) ) a ( v ¯ N + 1 ) h ( v ¯ N + 1 , v ¯ N + 1 ) ( ( 2 N 2 2 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ) ) .
The cubic logarithmic norm of the matrix D is equal to
Λ 2 ( D ) = max max 1 k N 1 1 d k k + l = 0 N + 1 | h ( v ¯ k , v ¯ l ) | 1 1 ζ l ( τ ) d τ ( τ v ¯ k ) 2 , max N 1 + 2 k N d k k + l = 0 N + 1 | h ( v ¯ k , v ¯ l ) | 1 1 ζ l ( τ ) d τ ( τ v ¯ k ) 2 , d 00 + l = 1 N + 1 | h ( v ¯ 0 , v ¯ l ) | 1 1 ζ l ( τ ) d τ ( τ v ¯ 0 ) 2 , d N 1 N 1 + l = 0 N + 1 | h ( v ¯ N 1 , v ¯ l ) | 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 ) 2 , d N 1 + 1 , N 1 + 1 + l = 0 N + 1 | h ( v ¯ N 1 + 1 , v ¯ l ) | 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 + 1 ) 2 , d N + 1 , N + 1 + l = 0 N | h ( v ¯ N + 1 , v ¯ l ) | 1 1 ζ l ( τ ) d τ ( τ v ¯ N + 1 ) 2 .
From (25)–(36) it follows that for sufficiently large N Λ 2 ( D ) < 0 occurs. By Theorem 2 it is clear that the system (43)–(48) (and (37)–(42)) has a unique solution x N ( t ) and D 1 1 / | Λ 2 ( D ) | .
Let x ( t ) and x N be solutions of (16 ) and (37)–(42), respectivety.
We recall the following definitions.
Definition 3.
The class W r ( M , [ a , b ] ) , r = 1 , 2 , , consists of all functions f C ( [ a , b ] ) , which have an absolutely continuous derivative f ( r 1 ) ( x ) and piecewise derivative f ( r ) ( x ) with | f ( r ) ( x ) | M .
Definition 4.
Denote by W r ( f : f 1 , f 2 ; M , c ) , r = 1 , 2 , , a set of functions f ( x ) , x [ a , b ] , such that f ( x ) = f 1 ( x ) , x [ a , c ) , f ( x ) = f 2 ( x ) , x ( c , b ] , where f 1 ( x ) W r ( M , [ a , c ] ) , f 2 ( x ) W r ( M , [ c , b ] ) , f 1 ( c ) f 2 ( c ) , c ( a , b ) .
Repeating the proof presented in [24] we see that the approximation of f ( t ) W 1 ( ( f : f 1 , f 2 ; M , c ) ) by piecewise linear functions constructed on the basis ζ k ( t ) , k = 0 , 1 , , N + 1 , has the error C N max ( ω ( f 1 ( 1 ) , 1 N ) , ω ( f 1 ( 1 ) , 1 N ) ) for f ( t ) W 1 ( ( f : f 1 , f 2 ; M , c ) ) , and C N 2 for f ( t ) W 2 ( ( f : f 1 , f 2 ; M , c ) ) .
In this paper, we denote the constants that do not depend on N by C.
Let x ( t ) W 2 ( ( x : x 1 , x 2 ; M , c ) ) , and x 1 ( 1 ) ( t ) C M 1 , t [ a , c ] , x 2 ( 1 ) ( t ) C M 2 , t [ c , b ] , M = m a x ( M 1 , M 2 ) , 0 < M < , where M is a bounded constant.
Repeating the arguments given in [24], we arrive at the following statement.
Theorem 6.
Let the following conditions be fulfilled:
(1) 
Equation (16) has the unique solution x ( t ) W 2 ( x 1 , x 2 ; M , c ) , 1 < c < 1 , M = c o n s t .
(2) 
For all t [ 1 , 1 ] the function h ( t , t ) 0 .
(3) 
Λ 2 ( D ) < 0 .
Then the system of Equations (37)–(42) has a unique solution x N ( t ) and the following estimate holds: | | x x N | | 1 C N 1 ln N .

3.2. Nonlinear Hypersingular Integral Equations

Consider the nonlinear hypersingular integral equation:
a ( t ) x ( t ) + 1 1 h ( t , τ , x ( τ ) ) d τ ( τ t ) 2 = f ( t ) .
The approximate solution of the Equation (49) we shall seek as a continuous function (17) with the coefficients γ k . The coefficients γ k are determined by the following system of nonlinear algebraic equations
a ( v ¯ k ) γ k + l = 0 N + 1 h ( v ¯ k , v l , γ l ) 1 1 ζ l ( τ ) ( τ v ¯ k ) 2 d τ = f ( v ¯ k ) , k = 0 , 1 , , N + 1 .
Remark 4.
Note that the set γ k , k = 0 , 1 , , N + 1 , is union of sets α k , k = 0 , 1 , , N 1 , and β k , k = 0 , 1 , , N 2 .
By computing the hypersingular integrals in (50), we can rewrite the system of Equation (50) as
a ( v ¯ k ) γ k h ( v ¯ k , v ¯ k , γ k ) 2 N 1 2 ln ( N 1 1 ) ( 1 + c ) N 1 2 + l = 0 N + 1 h ( v ¯ k , v ¯ l , γ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ k ) 2 = f ( v ¯ k ) , k = 1 , , N 1 1 ;
a ( v ¯ k ) γ k h ( v ¯ k , v ¯ k , γ k ) 2 N 2 2 ln ( N 2 1 ) ( 1 c ) N 2 2 + l = 0 N + 1 γ l h ( v ¯ k , v ¯ l , γ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ k ) 2 = f ( v ¯ k ) , k = N 1 + 2 , , N ;
a ( v ¯ 0 ) γ 0 h ( v ¯ 0 , v ¯ 0 , γ 0 ) ( 2 N 1 2 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ) + l = 1 N + 1 h ( v ¯ k , v ¯ l , γ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ 0 ) 2 = f ( v ¯ 0 ) ;
a ( v ¯ N 1 ) γ N 1 h ( v ¯ N 1 , v ¯ N 1 , γ N 1 ) ( 2 N 1 2 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ) + l = 0 N + 1 h ( v ¯ N 1 , v ¯ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 ) 2 = f ( v ¯ N 1 ) ;
a ( v ¯ N 1 + 1 ) γ N 1 + 1 h ( v ¯ N 1 + 1 , v ¯ N 1 + 1 , γ N 1 + 1 ) ( 2 N 2 2 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ) + l = 0 N + 1 h ( v ¯ N 1 + 1 , v ¯ l , γ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 + 1 ) 2 = f ( v ¯ N 1 + 1 ) ;
a ( v ¯ N + 1 ) γ N + 1 h ( v ¯ N + 1 , v ¯ N + 1 , γ N + 1 ) ( 2 N 2 2 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ) + l = 0 N h ( v ¯ N + 1 , v ¯ l , γ l ) 1 1 ζ l ( τ ) d τ ( τ v ¯ N + 1 ) 2 = f ( v ¯ N + 1 ) .
Here , , indicate summations over l k , l N 1 , l N 1 + 1 , respectively.
The Frechet derivative on a vector ( α ¯ 0 , α ¯ 1 , , α ¯ N + 1 ) in the space R N + 1 is equal to
a ( v ¯ k ) γ k h 3 ( v ¯ k , v ¯ k , γ ¯ k ) 2 N 1 2 ln ( N 1 1 ) ( 1 + c ) N 1 2 γ k + l = 0 N + 1 h 3 ( v ¯ k , v ¯ l , γ ¯ l ) γ l 1 1 ζ l ( τ ) d τ ( τ v ¯ k ) 2 , k = 1 , , N 1 1 ; a ( v ¯ k ) γ k h 3 ( v ¯ k , v ¯ k , γ ¯ k ) 2 N 2 2 ln ( N 2 1 ) ( 1 c ) N 2 2 γ k + l = 0 N + 1 h 3 ( v ¯ k , v ¯ l , γ ¯ l ) γ l 1 1 ζ l ( τ ) d τ ( τ v ¯ k ) 2 , k = N 1 + 2 , , N ; a ( v ¯ 0 ) γ 0 h 3 ( v ¯ 0 , v ¯ 0 , γ ¯ 0 ) γ 0 ( 2 N 1 2 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ) + l = 1 N + 1 h 3 ( v ¯ 0 , v ¯ l , γ ¯ l ) γ l 1 1 ζ l ( τ ) d τ ( τ v ¯ 0 ) 2 ; a ( v ¯ N 1 ) γ N 1 h 3 ( v ¯ N 1 , v ¯ N 1 , γ ¯ N 1 ) γ N 1 ( 2 N 1 2 + N 1 2 ln ( 2 ( 1 + c ) N 1 3 ) ( 1 + c ) N 1 2 ) + l = 0 N + 1 h 3 ( v ¯ N 1 , v ¯ l , γ ¯ l ) γ l 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 ) 2 ; a ( v ¯ N 1 + 1 ) γ N 1 + 1 h 3 ( v ¯ N 1 + 1 , v ¯ N 1 + 1 , γ ¯ N 1 + 1 ) γ N 1 + 1 ( 2 N 2 2 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ) + l = 0 N + 1 h 3 ( v ¯ N 1 + 1 , v ¯ l , γ ¯ l ) γ l 1 1 ζ l ( τ ) d τ ( τ v ¯ N 1 + 1 ) 2 ; a ( v ¯ N + 1 ) γ N + 1 h 3 ( v ¯ N + 1 , v ¯ N + 1 , γ ¯ N + 1 ) γ N + 1 ( 2 N 2 2 + N 2 2 ln ( 2 ( 1 c ) N 2 3 ) ( 1 c ) N 2 2 ) + l = 0 N h 3 ( v ¯ N + 1 , v ¯ l , γ ¯ l ) γ l 1 1 ζ l ( τ ) d τ ( τ v ¯ N + 1 ) 2 .
Here , indicate summations over l k , l N 1 , l N 1 + 1 , respectively.
The notation h 3 ( t , τ , u ) = δ h ( t , τ , u ) δ u is used here.
Let the Equation (49) has the unique solution x ( t ) inside the ball B ( x , δ ) . We shall assume that the Frechet derivative (57) in the ball R N + 1 ( x , δ ) satisfies the conditions of Theorem 5. Thus, according to statements of the Theorem 5, the solution of the system of differential equations
d α l ( σ ) d σ = a ( t ¯ l ) α l ( σ ) k = 0 N + 1 h ( t ¯ l , t ¯ k , α k ( σ ) ) N 1 2 l 2 k + 1 1 2 l 2 k 1 f ( t ¯ l ) , l = 0 , 1 , , N + 1 ,
converges to the solution of the Equation (49).
Thus, we have proven the following statement.
Theorem 7.
Let the following conditions hold:
(1) 
Equation (49) has a unique solution x ( t ) inside some ball B ( x , δ ) , x W 2 ( x : x 1 , x 2 ; M , c ) ;
(2) 
The Frechet derivative (57) in the ball R N + 1 ( x , δ ) satisfies the conditions of Theorem 5.
Then the system of Equations (51)–(56) has a unique solution inside the ball B ( x , δ ) , and the solution of Equation (58) converges to this solution.
The effectiveness of the presented algorithms is illustrated by solving two hypersingular integral equations modeling aerodynamics problems.
Example 1.
Let us illustrate the effectiveness of continuous method by solving the following linear hypersingular equation
1 1 x ( τ ) ( τ t ) 2 d τ = f ( γ 1 , γ 2 , t ) ,
where f ( γ 1 , γ 2 , t ) is the given right-hand side of the equation:
f ( γ 1 , γ 2 ; t ) = γ 1 γ 2 + ( a 1 a 2 ) 1 t ( a 1 + γ 1 t ) 1 1 + t ( a 2 + γ 2 t ) 1 1 t + γ 1 ln | t 1 + t | + γ 2 ln | 1 t t | .
The exact solution of the equation is x ( t ) = ( x 1 ( t ) , x 2 ( t ) ) ; x i ( t ) = a i + γ i t , i = 1 , 2 .
To solve the Equation (59) numerically we use the continuous method for solving operator equations and arrive to the following evolution equation
d α k ( σ ) d σ = l = 0 N + 1 N α l ( σ ) ( 1 2 k + 2 l 1 1 2 k + 2 l + 1 ) f ( γ 1 , γ 2 ; v ¯ k ) , k = 0 , 1 , , N + 1 .
Nodes v k , v ¯ k , k = 0 , 1 , , N + 1 , have been entered above.
In Figure 1 we show the trajectories of the exact solution of the Equation (59); its approximate solution, received with continuous method; and values of error.
Here a 1 = 1 , a 2 = 1.5 , γ 1 = 0.5 , γ 2 = 0.3 .
Example 2.
Let us illustrate the effectiveness of the continuous method for the solutions of nonlinear hypersingular equations
1 1 x 2 ( τ ) ( τ t ) 2 d τ = f ( γ 1 , γ 2 , t )
where f ( γ 1 , γ 2 , t ) is the given right-hand side of the equation:
f ( γ 1 , γ 2 ; t ) = γ 1 2 + 2 γ 1 a 1 + γ 1 2 t + a 1 2 t + ( 2 γ 1 a 1 + 2 γ 1 2 t ) ln | t 1 + t | ( a 1 2 + 2 γ 1 a 1 t + γ 1 2 t 2 ) 1 1 + t + γ 2 2 2 a 2 γ 2 γ 2 2 t a 2 2 t ( a 2 2 + 2 a 2 γ 2 t + γ 2 2 t 2 ) 1 1 t + ( 2 a 2 γ 2 + 2 γ 2 2 t ) ln | 1 t t | .
The exact solution of the equation is x ( t ) = ( x 1 ( t ) , x 2 ( t ) ) ; x i ( t ) = a i + γ i t , i = 1 , 2 .
It easy to see that, if x ( t ) is a solution of the Equation (60), then functions x ( t ) , | x ( t ) | and | x ( t ) | are solutions of this equation too.
To solve the Equation (60) numerically we use the continuous method and receive the following evolution equation
d α k ( σ ) d σ = l = 0 N + 1 N α l 2 ( σ ) ( 1 2 k + 2 l 1 1 2 k + 2 l + 1 ) f ( γ 1 , γ 2 , v ¯ l ) ,
k = 0 , 1 , , N + 1 .
At first, we take α k ( 0 ) = 0.0 as an initial condition in order to demonstrate applicability of our method in cases of the Newton–Kantorovich method, the minimal residual method and other numerical methods; using in their construction the derivative of nonlinear operator is not applicable. Indeed, in this case the Frechet derivative (57) is not only degenerate—and, therefore, not invertable—but is an identical zero.
In Figure 2 we put a 1 = 1 , a 2 = 1.4 , γ 1 = 0.5 , γ 2 = 0.4 .
In Figure 2 we show the trajectories of the exact solution of the Equation (60), its approximate solution, received with continuous method and values of error.
The exact solution at t = 0 has a jump discontinuity of h = 0.4 . The slopes of the exact solution also change at t = 0 . In Figure 2 we demonstrate that the numerical solution approximates the exact one at [ 1 , 0 ) well. At t = 0 the approximate solution has a jump h ˜ = 0.15 .

4. Summary and Discussion

An iterative projection method for solving linear and nonlinear hypersingular integral equations has been proposed. The method is based on the use of sufficient conditions for asymptotic stability of ODE systems. Stability conditions are expressed in terms of the logarithmic norms of the corresponding matrices. In a number of spaces often used in computational mathematics, the calculation of logarithmic norms does not cause difficulties, even for large-dimensional matrices.
What are the advantages of the presented method?
(1)
The method is applicable for solving linear and nonlinear hypersingular integral equations, whose right-hand sides contain non-Riemann integrable functions.
(2)
In Section 3.1 the continuous method is applied to linear hypersingular integral equations with the singularities of the second order. The conditions for the unique solvability of the constructed computing scheme are obtained and the convergence of the sequence of approximate solutions to the exact one is proven. It is shown that for linear hypersingular integral equations, the method converges for sufficiently large N and for b ( t ) 0 , t [ 1 , 1 ] .
(3)
In Section 3.2 the continuous method is applied to nonlinear hypersingular integral equations with the singularities of the second order. Conditions are given for the convergence of the constructed iterative spline-collocation method to the solution of a nonlinear hypersingular integral equation. It should be noted that the method is applicable to hypersingular integral equations of the first and second kinds.
The detailed bibliography of approximation methods of hypersingular integral equations of the first and the second kinds is given in [32]. The bibliography on solving hypersingular integral equations of the first kind is presented in [45].
Mostly, papers devoted to hypersingular integral equations of the first kind focused to seek solutions in the class of functions 1 t 2 ϕ ( t ) , where ϕ ( t ) is a smooth function. The presented method provides solutions in a general form.
The theoretical justification of the method is based on Lyapunov stability theory. It connects convergence of the method to the sign of the approximate system matrix logarithmic norm.
Said justification has advantages that allow us
1.
To obtain a set of convergence conditions owing to logarithmic norm values in various spaces;
2.
To determine the norm of the inverse matrix of an approximate system;
3.
To determine stability boundaries for solutions with respect to variations of kernels and right-hand sides of the equations.
The major advantage of the method for nonlinear equations is as follows.
The Newton–Kantorovich method requires the Frechet derivative reversibility at each iteration step. Similar conditions are required when using other iteration methods. Our method lacks such a deficiency. It does not put any restrictions on the Frechet derivative of the nonlinear operator.

Author Contributions

Conceptualization, I.B.; Data curation, V.R.; Formal analysis, I.B. and V.R.; Funding acquisition, V.R.; Investigation, I.B. and V.R.; Methodology, V.R.; Project administration, I.B.; Resources, V.R.; Software, V.R. and A.B.; Supervision, I.B.; Visualization, A.B.; Writing—original draft, I.B.; Writing—review and editing, V.R. and A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bogolubov, N.N.; Mesheryakov, V.A.; Tavkhelidze, A.N. An application of Muskhelishvili’s methods in the theory of elementary particles. In Proceedings of the Symposium on Solid Mechanics and Related Problems of Analysis, Metsniereba, Tbilisi; 1971; Volume 1, pp. 5–11. (In Russian). [Google Scholar]
  2. Brown, D.D.; Jackson, A.D. The Nucleon-Nucleon Interaction; North-Holland: Amsterdam, The Netherlands, 1976. [Google Scholar]
  3. Faddeev, L.; Takhtajan, L. Hamiltonian Approach to Solution Theory; Springer: Berlin, Germany, 1986. [Google Scholar]
  4. Ioakimidis, N.I. Two methods for the numerical solution of Bueckher’s singular integral equation for plane elasticity crack problems, Comput. Methods Appl. Mech. Engrg. 1982, 2, 169–177. [Google Scholar] [CrossRef]
  5. Lifanov, I.K.; Poltavskii, L.N.; Vainikko, G.M. Hypersingular Integral Equations and their Applications; Chapman Hall/CRC: Boca Raton, FL, USA; CRC Press Company: London, UK; New York, NY, USA; Washington, DC, USA, 2004. [Google Scholar]
  6. Mandal, B.N.; Chakrabarti, A. Applied Singular Integral Equations; CRC Press: Boca Raton, FL, USA, 2011; 260p. [Google Scholar]
  7. Boykov, I.V.; Boykova, A.I. Analytical methods for solution of hypersingular and polyhypersingular integral equations. arXiv 2019, arXiv:1901.04880v1. [Google Scholar]
  8. Gakhov, F.D. Boundary Value Problems; Dover Publication: Mineola, NY, USA, 1990; 561p. [Google Scholar]
  9. Ivanov, V.V. The Theory of Approximate Methods and their Application to the Numerical Solution of Singular Integral Equations; Noordhoff International Publishing: Leiden, The Netherlands, 1976. [Google Scholar]
  10. Gohberg, I.C.; Fel’dman, I.A. Convolution Equation and Projection Methods for Their Solution; Nauka: Moscow, Russia, 1971; English Translation: Translations of Mathematical Monographss Volume 41; American Mathematical Society: Providence, RI, USA, 1974. [Google Scholar]
  11. Golberg, M.A. Introduction to the numerical solution of Cauchy singular integral equations. In Mathematical Concepts and Methods in Science and Engineering; Plenum Press: New York, NY, USA, 1990; Volume 42. [Google Scholar]
  12. Michlin, S.G.; Prossdorf, S. Singular Integral Operatoren; Acad.-Verl.: Berlin, Germany, 1980. [Google Scholar]
  13. Boykov, I.V. Approximate Methods of Solution of Singular Integral Equations; The Penza State University: Penza, Russia, 2004. (In Russian) [Google Scholar]
  14. Boykov, I.V. Numerical methods for solutions of singular integral equations. arXiv 1973, arXiv:1610.09611. [Google Scholar]
  15. Capobiano, M.R.; Criscuolo, G.; Junghanns, P. On the numerical solution of a nonlinear integral equation of Prandtl’s type. In Operator Theory: Advances and Applications; Birkhauser Verlag: Basel, Switzerland, 2005; Volume 160, pp. 53–79. [Google Scholar]
  16. Oseledets, I.V.; Tyrtyshnicov, E.E. Approximate invention of matrices in the process of solving hypersingular integral equation. Comput. Math. Math. Phys. 2005, 45, 302–313. [Google Scholar]
  17. Golberg, M.A. The convergence of several algorithms for solving integral equations with finite-part integrals, I. Integral Equ. 1983, 5, 329–340. [Google Scholar] [CrossRef]
  18. Golberg, M.A. The convergence of several algorithms for solving integral equations with finite-part integrals, II. Integral Equ. 1985, 9, 267–275. [Google Scholar] [CrossRef]
  19. Lifanov, I.K. Singular Integral Equations and Discrete Vortices; VSP: Utrecht, The Netherlands, 1996. [Google Scholar]
  20. Boykov, I.V.; Zakharova Yu, F. An approximate solution to hypersingular integro-differential equations. Univ. Proc. Volga Reg. Phys. Math. Sci. Math. 2010, 1, 80–90. (In Russian) [Google Scholar]
  21. Boykov, I.V.; Boykova, A.I. An approximate solution of hypersingular integral equations with odd singularities of integer order. Univ. Proc. Volga Reg. Phys. Math. Sci. Math. 2010, 3, 15–27. (In Russian) [Google Scholar]
  22. Boykov, I.V.; Ventsel, E.S.; Boykova, A.I. An approximate solution of hypersingular integral equations. Appl. Numer. Math. 2010, 60, 607–628. [Google Scholar] [CrossRef]
  23. Boykov, I.V.; Ventsel, E.S.; Roudnev, V.A.; Boykova, A.I. An approximate solution of nonlinear hypersingular integral equations. Appl. Numer. Math. 2014, 86, 1–21. [Google Scholar] [CrossRef]
  24. Boykov, I.V.; Roudnev, V.A.; Boykova, A.I.; Baulina, O.A. New iterative method for solving linear and nonlinear hypersingular integral equations. Appl. Numer. Math. 2018, 127, 280–305. [Google Scholar] [CrossRef]
  25. Kaya, A.C.; Erdogan, F. On The Solution of Integral Equations with Strongly Singular Kernels. Quart. Appl. Math. 1987, 95, 105–122. [Google Scholar] [CrossRef] [Green Version]
  26. Eshkuvatov, Z.K.; Zulkarnain, F.S.; Nik Long, N.M.A.; Muminov, Z. Modified homotopy perturbation method for solving hypersingular integral equations of the first kind. SpringerPlus 2016, 5, 1473. [Google Scholar] [CrossRef] [Green Version]
  27. Boykov, I.V.; Boykova, A.I.; Syomov, M.A. An approximate solution of hypersingular integral equations of the first kind. Univ. Proc. Volga Reg. Phys. Math. Sci. Math. 2015, 3, 11–27. (In Russian) [Google Scholar]
  28. Boykov, I.V.; Boykova, A.I. Approximate methods for solving hypersingular integral equations of the first kind with second-order singularities on classes of functions with weights (1-t2)-1/2. Univ. Proc. Volga Reg. Phys. Math. Sci. Math. 2017, 2, 79–90. (In Russian) [Google Scholar]
  29. Boykov, I.V.; Boykova, A.I. An approximate solution of hypersingular integral equations of the first kind with singularities of the second order on the class of functions with weight ((1-t)/(1+t))±1/2. Univ. Proc. Volga Reg. Phys. Math. Sci. Math. 2019, 3, 76–92. (In Russian) [Google Scholar]
  30. Lifanov, I.K.; Nenashev, A.S. Hypersingular integral equations and the theory of wire antennas. Diff. Equ. 2005, 41, 126–145. [Google Scholar] [CrossRef]
  31. Lifanov, I.K.; Nenashev, A.S. Analysis of Some Computational Schemes for a Hypersingular Integral Equation on an Interval. Diff. Equ. 2005, 41, 1343–1348. [Google Scholar] [CrossRef]
  32. Boykov, I.V. Analytical and numerical methods for solving hypersingular integral equations. Dyn. Syst. 2019, 9, 244–272. (In Russian) [Google Scholar]
  33. Boikov, I.V. On a continuous method for solving nonlinear operator equations. Differ. Equ. 2012, 48, 1308–1314. [Google Scholar] [CrossRef]
  34. Daletskii Yu, L.; Krein, M.G. Stability of Solutions of Differential Equations in Banach Space; Nauka: Moscow, Russia, 1970; 536p, English Translation: Translations of Mathematical Monographss, Volume 43; American Mathematical Society: Providence, RI, USA, 1974. (In Russian) [Google Scholar]
  35. Lozinskii, S.M. Note on a paper by V.S. Godlevskii. Ussr Comput. Math. Math. Phys. 1973, 13, 232–234. [Google Scholar] [CrossRef]
  36. Boikov, I.V. On the stability of solutions of differential and difference equations in critical cases. Soviet Math. Dokl. 1991, 42, 630–632. [Google Scholar]
  37. Boikov, I.V. Stability of Solutions of Differential Equations; Publishing House of Penza State University: Penza, Russia, 2008; 24p. (In Russian) [Google Scholar]
  38. Browder, F.E. Nonlinear mappings of the nonexpansive and accretive type in Banach spaces. Bull. Am. Math. Soc. 1967, 73, 875–882. [Google Scholar] [CrossRef] [Green Version]
  39. Malkowski, E.; Rakoćević, V. Advanced Functional Analysis; CRS Press, Taylor and Francis Group: Boca Raton, FL, USA, 2019. [Google Scholar]
  40. Ćirić, L.; Rafiq, A.; Radenović, S.; Rajović, M. Jeong Sheok Ume On Mann implicit iterations for strongly accreative and strongly pseudo-contractive mappings. Appl. Math. Comput. 2008, 198, 128–137. [Google Scholar]
  41. Todoriević, V. Harmonic Quasiconformal Mappings and Hyperbolic Type Metrics; Springer Nature: Cham, Switzerland, 2019. [Google Scholar]
  42. Todoriević, V. Subharmonic behavior and quasiconformal mappings. Anal. Math. Phys. 2019, 9, 1211–1225. [Google Scholar] [CrossRef]
  43. Hadamard, J. Lectures on Cauchy’s Problem in Linear Partial Differential Equations; Dover Publ. Inc.: New York, NY, USA, 1952. [Google Scholar]
  44. Chikin, L.A. Special cases of the Riemann boundary value problems and singular integral equations. Sci. Notes Kazan State Univ. 1953, 113, 53–105. (In Russian) [Google Scholar]
  45. Chan, Y.-S.; Fannjiang, A.C.; Paulino, G.H. Integral equations with hypersingular kernels-theory and applications to fracture mechanics. Int. J. Eng. Sci. 2003, 41, 683–720. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Numerical solutions for the linear hypersingular equation with a discontinuous right-hand side example.
Figure 1. Numerical solutions for the linear hypersingular equation with a discontinuous right-hand side example.
Axioms 09 00074 g001
Figure 2. Numerical solution for the nonlinear hypersingular equation with a discontinuous right-hand side example.
Figure 2. Numerical solution for the nonlinear hypersingular equation with a discontinuous right-hand side example.
Axioms 09 00074 g002

Share and Cite

MDPI and ACS Style

Boykov, I.; Roudnev, V.; Boykova, A. Approximate Methods for Solving Linear and Nonlinear Hypersingular Integral Equations. Axioms 2020, 9, 74. https://doi.org/10.3390/axioms9030074

AMA Style

Boykov I, Roudnev V, Boykova A. Approximate Methods for Solving Linear and Nonlinear Hypersingular Integral Equations. Axioms. 2020; 9(3):74. https://doi.org/10.3390/axioms9030074

Chicago/Turabian Style

Boykov, Ilya, Vladimir Roudnev, and Alla Boykova. 2020. "Approximate Methods for Solving Linear and Nonlinear Hypersingular Integral Equations" Axioms 9, no. 3: 74. https://doi.org/10.3390/axioms9030074

APA Style

Boykov, I., Roudnev, V., & Boykova, A. (2020). Approximate Methods for Solving Linear and Nonlinear Hypersingular Integral Equations. Axioms, 9(3), 74. https://doi.org/10.3390/axioms9030074

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