Next Article in Journal
A Ferroelectric Memristor-Based Transient Chaotic Neural Network for Solving Combinatorial Optimization Problems
Previous Article in Journal
Analysis of a Productive Topological Index Correlated to Fullerenes’ Physical Properties
Previous Article in Special Issue
A Class of Adaptive Exponentially Fitted Rosenbrock Methods with Variable Coefficients for Symmetric Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Error Estimates of a Symmetric Spectral Method for a Linear Volterra Integral Equation

College of Mathematics and Statistics, Hanshan Normal University, Chaozhou 521041, China
*
Author to whom correspondence should be addressed.
Submission received: 16 June 2022 / Revised: 16 July 2022 / Accepted: 22 July 2022 / Published: 26 December 2022

Abstract

:
A symmetric spectral method is applied to investigate the two-dimensional Volterra integral equation with weakly singular kernels and delays. In this work, the solution of the equation we considered is assumed to be sufficiently smooth so that the spectral method can be applied naturally. Employing three couples of variable transformations, we apply the two-dimensional Gauss quadrature rule to approximate the weakly singular integral with delays and obtain the spectral discretization. Then we derive the convergence results of the proposed approximation scheme. We show that the errors of solution decay exponentially in both the infinity norm and weighted square norm. In the end, we carry out numerical experiments to verify the theoretical results.

1. Introduction

Volterra equations with weakly singular kernel and/or delay frequently appear in many areas of science and engineering, such as the reaction–diffusion study [1], memory problems [2], population growth models [3], the competitive system [4], a model for tumor growth [5] and so on. In view of the extensive applications, Volterra integral equations with weakly singular kernel and/or delay have drawn much attention in recent years.
The literature [6] is the pioneer work, in which the authors proposed a spectral Jacobi-collocation approximation for the linear Volterra integral equations of the second kind with weakly singular kernels under the sufficiently smooth solutions. While the underlying solution is not sufficiently smooth, this case was considered in [7] by the same authors one year later. In [8], D. Li et al. discussed the Volterra integro-differential equations which are derived by combining a time memory term in the flux and a delay parameter in the reaction term. It is well known that spectral methods enjoy exponential convergence and high-order accuracy. As a result, there are many references on numerically investigating Volterra equations by using spectral methods. The literature [9] investigated the convergence of the spectral methods for the weakly singular Volterra integro-differential equations with smooth solutions. After that, weakly singular Volterra integral equations with pantograph delays were studied in [10]. In [11], a kind of Legendre spectral scheme was put forward to solve nonlinear systems of Volterra integral equations. The literature [12] dealt with a fractional order collocation method for second kind Volterra integral equations with weakly singular kernels. Fractional pantograph equations were taken into account for their long time numerical behaviors in [13]. G. Deng et al. [14] developed high accurate pseudo-spectral Galerkin scheme for pantograph type Volterra integro-differential equations with singular kernels, but the underlying solution they considered was sufficiently smooth. Based on a simple change of variable, a novel numerical approach had been devoted to the time-fractional parabolic equations with non-smooth solutions in [15]. Liang and Brunner established the analogous convergence analysis for Volterra integral equations with weakly singular kernels for both uniform and graded meshes [16]. For a more detailed description of this subject, we refer the readers to the recent papers [17,18,19,20,21,22,23,24] and the references therein. For the two-dimensional Volterra equations, there are only a few introductions in [25,26]. To the best of our knowledge, none of the two papers refer to the two-dimensional Volterra integral equations with weakly singular kernels or delays. In this paper, we will use the spectral method to numerically solve the two-dimensional Volterra integral equation with weakly singular kernels and delays. The novelty of our work lies in applying the two-dimensional Gauss quadrature rule to approximate the weakly singular integral with delays and employing the spectral error analysis for the Jacobi spectral collocation method.
The two-dimensional Volterra integral equation with weakly singular kernels and delays we shall study is of the form:
0 p X 0 q Y ( p X S 1 ) μ ( q Y ξ 1 ) ρ K ˜ ( X , Y , S 1 , ξ 1 ) U ( S 1 , ξ 1 ) d ξ 1 d S 1 = U ( X , Y ) + G ( X , Y ) , X [ 0 , T ] , Y [ 0 , Γ ] ,
where U ( X , Y ) is the unknown function, K ˜ , G are given functions and 0 < μ , ρ , p , q 1 are given constants. S 1 and ξ 1 are two integral variables. All functions in the above equation are assumed to be sufficiently smooth. As we know, this kind of equation not only arises from biological systems, but also found applications in physics, chemistry and memory material. Without losing generality, we discuss the last model, where U ( X , Y ) stands for the initial shape. A temporary shape that is fixed is well illustrated by the double integral, in which the weak terms express the memory characteristic and two parameters p and q depend on the degree of deformation of the specific material in the X-axis and Y-axis, respectively. The deformation of the specific material is described by K ˜ . In the end, the function G ( x , y ) represents the external force.
In order to apply the theory of orthogonal polynomials, we need to employ some linear transformations. First, let S 1 = p S , ξ 1 = q ξ and (1) become:
p q 0 X 0 Y ( p X p S ) μ ( q Y q ξ ) ρ K ˜ ( X , Y , p S , q ξ ) U ( p S , q ξ ) d ξ d S = U ( X , Y ) + G ( X , Y ) .
Then we make the changes of variables X = T 2 ( 1 + x ) and Y = Γ 2 ( 1 + y ) , where ( x , y ) Ω ¯ ( Ω = ( 1 , 1 ) × ( 1 , 1 ) ), and the above equation turns into:
0 T 2 ( 1 + x ) 0 Γ 2 ( 1 + y ) p T 2 ( 1 + x ) p S μ q Γ 2 ( 1 + y ) q ξ ρ K ¯ ( x , y , p S , q ξ ) U ( p S , q ξ ) d ξ d S = u ( x , y ) + g ( x , y ) ,
where
u ( x , y ) = U T 2 ( 1 + x ) , Γ 2 ( 1 + y ) , K ¯ ( x , y , p S , q ξ ) = p q K ˜ T 2 ( 1 + x ) , Γ 2 ( 1 + y ) , p S , q ξ , g ( x , y ) = G T 2 ( 1 + x ) , Γ 2 ( 1 + y ) .
Later, two linear transformations S = T 2 ( 1 + s ) and ξ = Γ 2 ( 1 + τ ) ( s , τ [ 1 , 1 ] ) are employed and (3) changes into:
1 x 1 y ( x s ) μ ( y τ ) ρ K ^ ( x , y , s , τ ) u ( p s + p 1 , q τ + q 1 ) d τ d s = u ( x , y ) + g ( x , y ) ,
where
K ^ ( x , y , s , τ ) = p T 2 1 μ q Γ 2 1 ρ K ˜ T 2 ( 1 + x ) , Γ 2 ( 1 + y ) , p T 2 ( 1 + s ) , q Γ 2 ( 1 + τ ) .
Spectral discretization will be carried out below. We first introduce the two-dimensional collocation points. Then we apply some linear transformations to obtain an equivalent problem to get the high-order accuracy. At last, the Gauss integration formula and Lagrange interpolation polynomial are employed to gain the high precision numerical format.
Firstly, we denote the collocation points in two-dimensional space by { ( x i , y j ) } i , j = 0 N , where { x i } i = 0 N and { y j } j = 0 N are the Gauss points in one dimension, equipped with the weights ω i x and ω j y , respectively. Equation (4) holds at the two-dimensional collocation points, namely
1 x i 1 y j ( x i s ) μ ( y j τ ) ρ K ^ ( x i , y j , s , τ ) u ( p s + p 1 , q τ + q 1 ) d τ d s = u ( x i , y j ) + g ( x i , y j ) .
To apply Gauss quadrature rule, we need to transfer the two integral intervals [ 1 , x i ] and [ 1 , y j ] into [ 1 , 1 ] , respectively, by using the following changes of variables:
s = s ( θ ) = x i + 1 2 θ + x i 1 2 , τ = τ ( η ) = y j + 1 2 η + y j 1 2 ,
where i , j = 1 , 2 , , N , θ , η [ 1 , 1 ] , then (5) becomes:
1 1 1 1 ( 1 θ ) μ ( 1 η ) ρ k ^ ( x i , y j , s ( θ ) , τ ( η ) ) u ( p s ( θ ) + p 1 , q τ ( η ) + q 1 ) d η d θ = u ( x i , y j ) + g ( x i , y j ) ,
where
k ^ ( x i , y j , s ( θ ) , τ ( η ) ) = x i + 1 2 1 μ y j + 1 2 1 ρ K ^ ( x i , y j , s ( θ ) , τ ( η ) ) .
The Gauss integration formula tells that the integration term in (7) can be approximated by:
1 1 1 1 ( 1 θ ) μ ( 1 η ) ρ k ^ ( x i , y j , s ( θ ) , τ ( η ) ) u ( p s ( θ ) + p 1 , q τ ( η ) + q 1 ) d η d θ k , l = 0 N k ^ ( x i , y j , s ( θ k ) , τ ( η l ) ) u ( p s ( θ k ) + p 1 , q τ ( η l ) + q 1 ) ω k x ω l y .
I N u ( x , y ) denotes the Lagrange interpolation polynomial, i.e.,
I N u ( x , y ) = i , j = 0 N u ( x i , y j ) F i j ( x , y ) ,
where F i j ( x , y ) = F i ( x ) F j ( y ) , F i and F j stand for two Lagrange interpolation basis functions which are based on the collocation points { x i } i = 0 N and { y j } j = 0 N , respectively. I N u satisfies
I N u ( x i , y j ) = u ( x i , y j ) .
Denote u i j u ( x i , y j ) and u N ( x , y ) u ( x , y ) , where
u N ( x , y ) = i , j = 0 N u i j F i j ( x , y ) .
Then the spectral collocation method tells that all of the u i j , i , j = 0 , . . . , N can be found by solving the following linear system
i , j = 0 N u i j k , l = 0 N k ^ ( x i , y j , s ( θ k ) , τ ( η l ) ) F i ( p s ( θ k ) + p 1 ) F j ( q τ ( θ l ) + q 1 ) ω k x ω l y = u i j + g ( x i , y j ) .
Define U = ( u i j ) T and G = ( g ( x i , y j ) ) T . We rewrite (10) as the following matrix form:
A U = U + G ,
where
A i j = k , l = 0 N k ^ ( x i , y j , s ( θ k ) , τ ( η l ) ) F i ( p s ( θ k ) + p 1 ) F j ( q τ ( θ l ) + q 1 ) ω k x ω l y .
For any given k ^ and G, we can obtain the values of u i j ( i , j = 0 , , N ) by solving the matrix Equation (11). Meanwhile, the approximate solution u N ( x , y ) can also be obtained by the expression of (9).
Now we give the arrangement of the following sections. Three kinds of norms and seven lemmas are described in Section 2. Error estimates are proposed in Section 3. Section 4 contains numerical examples to confirm the theoretical analysis. At last, we end the paper with some conclusions.

2. Some Preliminaries

In this section, some elementary lemmas will be provided, which are important for establishing the error estimates in the subsequent section. Before giving the lemmas, we introduce three kinds of norms one by one.
The first kind of norm v denotes the infinite norm in the usual sense, given by
v = e s s sup ( x , y ) Ω | v ( x , y ) | .
The second kind of norm u ω α , β represents the weighted square norm which is endowed in the space:
L ω α , β 2 ( Ω ) = u : 1 1 1 1 | u ( x , y ) | 2 ω α , β ( x , y ) d y d x < + ,
whose norm is defined as:
u ω α , β = 1 1 1 1 | u ( x , y ) | 2 ω α , β ( x , y ) d y d x 1 2 ,
and the inner product is given by:
( u , v ) ω α , β = 1 1 1 1 u ( x , y ) v ( x , y ) ω α , β ( x , y ) d y d x ,
where ω α , β ( x , y ) = ( 1 x ) α ( 1 + x ) β ( 1 y ) α ( 1 + y ) β , α , β ( 1 , 0 ) .
The third kind of norm is denoted as v ω α , β m = k = 0 m k v ω α , β 2 1 2 , which is equipped in the space
H ω α , β m ( Ω ) = v : k v L ω α , β 2 ( Ω ) , 0 k m ,
along with the semi-norm
| v | ω α , β m ; N = k = min ( m , N + 1 ) m k v ω α , β 2 1 2 .
Then, seven important lemmas will be introduced. They will play key roles in the sharp error estimates for the fully discrete scheme in the next section. Throughout the paper, C denotes a generic positive constant, not necessarily the same at different places, which is dependent on the given data and the solution, but independent of the number of collocation points.
Lemma 1
([27] Gronwall inequality). Suppose a nonnegative integrable function E ( x , y ) satisfies
E ( x , y ) M 1 x 1 y ( x s ) μ ( y τ ) ρ E ( s , τ ) d τ d s + G ( x , y ) , ( x , y ) Ω ,
where M is a non-negative constant and G ( x , y ) is a non-negative integrable function too. Then we have E C G and E ω α , β C G ω α , β .
Lemma 2
([28]). Define | | I N | | = max ( x , y ) Ω ¯ i , j = 0 N | F i ( x ) F j ( y ) | and it holds that:
| | I N | | = O log 2 N , 1 < α , β 1 2 , O N 2 max ( α , β ) + 1 , otherwise .
Lemma 3
([29]). Assume P N denotes a space in which the degrees of all polynomials not exceed N and the product v ϕ is integrated by the Gauss quadrature formula, where v H ω α , β m ( Ω ) for some m > 1 and ϕ P N . Then
( v , ϕ ) ω α , β ( v , ϕ ) N C N m | v | ω α , β m ; N ϕ ω α , β ,
where
( v , ϕ ) N = i , j = 0 N v ( θ i , η j ) ϕ ( θ i , η j ) ω i x ω j y .
Lemma 4
([29,30]). Assume that u ( x , y ) H ω α , β m ( Ω ) for m > 1 . Then we have the following estimates:
u I N u C N 4 m | u | ω α , β m ; N ,
u I N u ω α , β C N m | u | ω α , β m ; N .
Lemma 5
([31]). Assume v ( x , y ) is a bounded function and there is:
sup N i , j = 0 N v ( x i , y j ) F i j ( x , y ) ω α , β C max ( x , y ) Ω ¯ | v ( x , y ) | ,
where C is a constant independent of v ( x , y ) .
Lemma 6
([32,33]). For a nonnegative integer r and κ ( 0 , 1 ) , there exists a constant C r , κ > 0 such that for any function v C r , κ ( Ω ¯ ) , there is a polynomial function J N v P N such that
v J N v C r , κ N ( r + κ ) v C r , κ ( Ω ¯ ) ,
where
v C r , κ ( Ω ¯ ) = max α r max ( x , y ) Ω ¯ | α v ( x , y ) | + max α r sup ( x , y ) ( x , y ) Ω ¯ α v ( x , y ) α v ( x , y ) [ ( x x ) 2 + ( y y ) 2 ] κ / 2 .
It is obvious that J N is a linear operator from C r , κ ( Ω ¯ ) into P N .
Lemma 7
([34]). For any continuous function v defined in C ( [ 1 , 1 ] × [ 1 , 1 ] ) , assume M v is defined by:
( M v ) ( x , y ) = 1 x 1 y ( x s ) μ ( y τ ) ρ k ^ ( x , y , s , τ ) v ( s , τ ) d τ d s .
For any different points ( x 1 , y 1 ) , ( x 2 , y 2 ) [ 1 , 1 ] × [ 1 , 1 ] and δ ( 0 , min { 1 μ , 1 ρ } ) , it is valid that:
| ( M v ) ( x , y ) ( M v ) ( x , y ) | [ ( x x ) 2 + ( y y ) 2 ] κ / 2 C max ( x , y ) Ω ¯ | v ( x , y ) | .
That is to say:
M v C 0 , κ ( Ω ¯ ) C max ( x , y ) Ω ¯ | v ( x , y ) | .

3. Error Estimates

We present the error estimates for the numerical scheme (10) and derive that the proposed approximation possesses the desired exponential rate of convergence. Firstly, we carry out the error estimate in the sense of the L norm.
Theorem 1.
Suppose u ( x , y ) is the exact solution of (4), which is supposed to be sufficiently smooth, u N ( x , y ) is the approximate solution of u ( x , y ) and the error function is denoted as e ( x , y ) = u ( x , y ) u N ( x , y ) . If m > 4 , then we claim that:
e C N m log 2 N max ( x , y ) Ω ¯ | k ^ ( x , y , s , τ ) | ω α , β m ; N u ω α , β + N 4 | u | ω α , β m ; N , α , β ( 1 , 1 2 ] , N 2 max ( α , β ) + 1 max ( x , y ) Ω ¯ | k ^ ( x , y , s , τ ) | ω α , β m ; N u ω α , β + N 4 | u | ω α , β ( Ω ) m ; N , otherwise .
Proof. 
Subtracting (10) from (7), we get that:
1 1 1 1 ( 1 θ ) μ ( 1 η ) ρ k ^ ( x i , y j , s ( θ ) , τ ( η ) ) u ( p s ( θ ) + p 1 , q τ ( η ) + q 1 ) d η d θ i , j = 0 N u i j k , l = 0 N k ^ ( x i , y j , s ( θ k ) , τ ( η l ) ) L i ( p s ( θ k ) + p 1 ) L j ( q τ ( η l ) + q 1 ) ω k x ω l y = u ( x i , y j ) u i j ,
and we can rewrite (25) into the following form:
1 1 1 1 ( 1 θ ) μ ( 1 η ) ρ k ^ ( x i , y j , s ( θ ) , τ ( η ) ) e ( p s ( θ ) + p 1 , q τ ( η ) + q 1 ) d η d θ + B ( x i , y j ) = u ( x i , y j ) u i j ,
by defining
B ( x i , y j ) = 1 1 1 1 ( 1 θ ) μ ( 1 η ) ρ k ^ ( x i , y j , s ( θ ) , τ ( η ) ) u N ( p s ( θ ) + p 1 , q τ ( η ) + q 1 ) d η d θ i , j = 0 N u i j k , l = 0 N k ^ ( x i , y j , s ( θ k ) , τ ( η l ) ) L i ( p s ( θ k ) + p 1 ) L j ( q τ ( θ l ) + q 1 ) ω k x ω l y .
Noting (6), we have that
1 x i 1 y j ( x i s ) μ ( y j τ ) ρ K ^ ( x i , y j , s , τ ) e ( p s + p 1 , q τ + q 1 ) d τ d s + B ( x i , y j ) = u ( x i , y j ) u i j .
Multiplying F i ( x ) F j ( y ) on both sides of (27), then summing up the indexes i and j from 0 to N, respectively, and the error equation can be obtained as:
e ( x , y ) = 1 x 1 y ( x s ) μ ( y τ ) ρ K ^ ( x , y , s , τ ) e ( p s + p 1 , q τ + q 1 ) d τ d s + m = 1 3 B m ( x , y ) ,
here
B 1 ( x , y ) = i , j = 0 N B ( x i , y j ) F i ( x ) F j ( y ) , B 2 ( x , y ) = u ( x , y ) I N u ( x , y ) , B 3 ( x , y ) = I N 1 x 1 y ( x s ) μ ( y τ ) ρ K ^ ( x , y , s , τ ) e ( p s + p 1 , q τ + q 1 ) d τ d s 1 x 1 y ( x s ) μ ( y τ ) ρ K ^ ( x , y , s , τ ) e ( p s + p 1 , q τ + q 1 ) d τ d s .
As a result, we obtain that
| e ( x , y ) | M 1 x 1 y ( x s ) μ ( y τ ) ρ | e ( p s + p 1 , q τ + q 1 ) | d τ d s + m = 1 3 | B m ( x , y ) | ,
along with
M = max 1 s x 1 , 1 τ y 1 | K ^ ( x , y , s , τ ) | .
Using Lemma 1, we get that
| | e | | C | | B 1 | | + | | B 2 | | + | | B 3 | | .
Applying Lemmas 2 and 3 yields:
| | B 1 | | C | | I N | | max 0 i , j N | B ( x i , y j ) | C | | I N | | N m max 0 i , j N | k ^ ( x i , y j , s , τ ) | ω α , β m ; N ( u ω α , β + e ) C N m max 0 i , j N | k ^ ( x i , y j , s , τ ) | ω α , β m ; N ( u ω α , β + e ) log 2 N , α , β ( 1 , 1 2 ] , N 2 max ( α , β ) + 1 , otherwise .
By using the estimate (17) in Lemma 4, we have that:
| | B 2 | | C N 4 m | u | ω α , β m ; N .
Following Lemmas 6 and 7, we arrive at:
B 3 = ( I I N ) M e = ( I I N ) M e ( I I N ) J N M e = ( I I N ) ( M e J N M e ) ( 1 + | | I N | | ) M e J N M e C I N N κ M e C 0 , κ ( Ω ) ¯ C I N N κ e ,
where in the second step of (32) we notice the fact that ( I I N ) J N M e = J N M e J N M e = 0 .
Combining (29)–(32), we can get the estimate for | | e | | immediately. □
Next, we carry out the error analysis in L ω α , β 2 space.
Theorem 2.
If the hypotheses given in Theorem 1 hold, then:
| | e | | ω α , β C N m max ( x , y ) Ω ¯ | k ^ ( x , y , s , τ ) | ω α , β m ; N u ω α , β + | u | ω α , β m ; N + N κ e .
Proof. 
With the help of (28) and Gronwall inequality, we have that:
| | e | | ω α , β C | | B 1 | | ω α , β + | | B 2 | | ω α , β + | | B 3 | | ω α , β .
By virtue of Lemmas 1 and 5, we get that:
| | B 1 | | ω α , β C max ( x , y ) Ω ¯ | B ( x , y ) | C N m max 0 i , j N | k ^ ( x i , y j , s , τ ) | ω α , β m ; N u ω α , β + e ω α , β .
Applying Lemma 4 results in:
| | B 2 | | ω α , β C N m | u | ω α , β m ; N .
By using Lemmas 5–7, we find that:
| | B 3 | | ω α , β = | | ( I I N ) M e | | ω α , β = | | ( I I N ) ( M e J N M e ) | | ω α , β = | | I ( M e J N M e ) | | ω α , β + | | I N ( M e J N M e ) | | ω α , β | | M e J N M e | | ω α , β + | | I N ( M e J N M e ) | | ω α , β C | | M e J N M e | | C N κ | | M e | | C r , κ ( Ω ¯ ) C N κ | | e | | .
Together with (34)–(37), and we can get the desired estimate. □

4. Error Tests

In this section, two examples will be presented to confirm the theoretical analysis. We perform the computations by using Matlab software.
Example 1.
Let p = q = 1 , then (4) has the following form:
1 x 1 y ( x s ) μ ( y τ ) ρ K ^ ( x , y , s , τ ) u ( s , τ ) d τ d s = u ( x , y ) + g ( x , y ) .
Here we choose u ( x , y ) = 2 e cos ( x y ) as the exact solution and select
K ^ ( x , y , s , τ ) = ( x s ) 2 ( y τ ) e cos ( s τ ) ,
to this end
g ( x , y ) = 2 ( x + 1 ) 3 μ ( y + 1 ) 2 ρ ( 3 μ ) ( 2 ρ ) 2 e cos ( x y ) .
Let μ = 0.5 , ρ = 1 / 3 , and we solve the matrix Equation (11), and the infinity error in L norm and weighted square error in L ω α , β 2 norm for different N which denotes the number of collocation points are listed in Table 1. To show the exponential rate of convergence more visually, Figure 1 exhibits the errors in the L and L ω α , β 2 norms versus N. From Table 1 and Figure 1, one can see that the spectral method achieves the exponential rate of convergence, illustrating a fact that the numerical simulations highly agree with the theoretical results.
Example 2.
We consider a normal situation. We take:
u ( x , y ) = x y
and
K ^ ( x , y , s , τ ) = 1
into Equation (4) and it becomes:
1 x 1 y ( x s ) μ ( y τ ) ρ ( p s + p 1 ) ( q τ + q 1 ) d τ d s = x y + g ( x , y ) .
By calculation, g ( x , y ) is determined by:
g ( x , y ) = p ( x + 1 ) 2 μ ( 1 μ ) ( 2 μ ) ( x + 1 ) 1 μ 1 μ q ( y + 1 ) 2 ρ ( 1 ρ ) ( 2 ρ ) ( y + 1 ) 1 ρ 1 ρ x y .
As the theoretical analysis shown, if four parameters μ , ρ , p , q are given, we can carry out the spectral numerical test to yield the approximate solution. Now we let μ = 0.8 , ρ = 0.1 , p = 0.65 , q = 1 / 3 . Then we carry out the computation and get the errors in L and L ω α , β 2 norms. More information please see Table 2 and Figure 2. Again, we conclude that the spectral method exhibits the exponential rate of convergence. These simulations further confirm our theoretical results.

5. Conclusions

In this paper, we have presented the error estimates for a weakly singular Volterra integral equation with delays in two-dimensional space. The main contribution of this work is that we have demonstrated rigorously that the errors of spectral approximations decay exponentially in both the L norm and the L ω α , β 2 norm, which is a desired feature for a spectral method. Two numerical examples verify the theoretical results and exhibit an improvement of accuracy over the existing method [35]. Next we will extend our work to two-dimensional Volterra integro-differential equations.

Author Contributions

Software, W.Z.; Writing—original draft, Y.C.; Writing—review and editing, D.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by Research Projects of Guangdong Provincial Education Department (2021KTSCX071, 2022KTSCX077, QD202212, HSGDJG21356-372) and a Project of Hanshan Normal University (521036).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dixon, J.A. A nonlinear weakly singular Volterra integro-differential equation arising from reaction-diffusion study in a small cell. J. Comput. Appl. Math. 1987, 18, 289–305. [Google Scholar] [CrossRef]
  2. Shaw, S.; Warby, M.K.; Whiteman, J.R. Error estimates with sharp constants for a fading memory Volterra problem in linear solid viscoelasticity. SIAM J. Numer. Anal. 1997, 34, 1237–1254. [Google Scholar] [CrossRef]
  3. Al-Khaled, K. Numerical approximations for population growth models. Appl. Math. Comput. 2005, 160, 865–873. [Google Scholar] [CrossRef]
  4. Zhang, C.; Yan, X. Positive solution bifurcating from zero solution in a Lotka-Volterra competitive system with cross-diffusion effects. Appl. Math. J. Chin. Univ. 2011, 26, 342–352. [Google Scholar] [CrossRef]
  5. Villasana, M.; Radunskaya, A. A delay differential equation model for tumor growth. J. Math. Biol. 2003, 47, 270–294. [Google Scholar] [CrossRef]
  6. Chen, Y.; Tang, T. Spectral methods for weakly singular Volterra integral equations with smooth solutions. J. Comput. Appl. Math. 2009, 233, 938–950. [Google Scholar] [CrossRef] [Green Version]
  7. Chen, Y.; Tang, T. Convergence analysis of the Jacobi spectral-collocation methods for Volterra integral equations with a weakly singular kernel. Math. Comput. 2010, 79, 147–167. [Google Scholar] [CrossRef] [Green Version]
  8. Li, D.; Zhang, C.; Wang, W. Long time behavior of non-Fickian delay reaction-diffusion equations. Nonlinear Anal-Real. 2012, 13, 1401–1415. [Google Scholar] [CrossRef]
  9. Wei, Y.; Chen, Y. Convergence analysis of the spectral methods for weakly singular Volterra integro-differential equations with smooth solutions. Adv. Appl. Math. Mech. 2012, 4, 1–20. [Google Scholar] [CrossRef]
  10. Zhang, R.; Zhu, B.; Xie, H. Spectral methods for weakly singular Volterra integral equations with pantograph delays. Front. Math. China 2013, 8, 281–299. [Google Scholar] [CrossRef]
  11. Tohidi, E.; Navid, O.R.; Shateyi, S. Convergence analysis of Legendre Pseudospectral scheme for solving nonlinear systems of Volterra integra equations. Adv. Math. Phys. 2014, 2014, 307907. [Google Scholar] [CrossRef] [Green Version]
  12. Cai, H.; Chen, Y. A fractional order collocation method for second kind Volterra integral equations with weakly singular kernels. J. Sci. Comput. 2018, 75, 970–992. [Google Scholar] [CrossRef]
  13. Li, D.; Zhang, C. Long time numerical behaviors of fractional pantograph equations. Math. Comput. Simul. 2020, 172, 244–257. [Google Scholar] [CrossRef]
  14. Deng, G.; Yang, Y.; Tohidi, E. High accurate pseudo-spectral Galerkin scheme for pantograph type Volterra integro-differential equations with singular kernels. Appl. Math. Comput. 2021, 396, 125866. [Google Scholar] [CrossRef]
  15. Li, D.; Sun, W.; Wu, C. A novel numerical approach to time-fractional parabolic equations with nonsmooth solutions. Numer. Math. Theor. Meth. Appl. 2021, 14, 355–376. [Google Scholar]
  16. Liang, H.; Brunner, H. The convergence of collocation solutions in continuous piecewise polynomial spaces for weakly singular Volterra integral equations. SIAM J. Numer. Anal. 2019, 57, 1875–1896. [Google Scholar] [CrossRef]
  17. Shokri, A. A symmetric P-stable hybrid Obrechkoff methods for the numerical solution of second order IVPS. TWMS J. Pure Appl. Math. 2014, 5, 28–35. [Google Scholar]
  18. Shokri, A.; Saadat, H. P-stability, TF and VSDPL technique in Obrechkoff methods for the numerical solution of the Schrodinger equation. Bull. Iran. Math. Soc. 2016, 42, 687–706. [Google Scholar]
  19. Gu, Z. Chebyshev spectral collocation method for system of nonlinear volterra integral equations. Numer. Algorithms 2020, 83, 243–263. [Google Scholar] [CrossRef]
  20. Yaghoobnia, A.R.; Khodabin, M.; Ezzati, R. Numerical solution of stochastic It o ^ -Volterra integral equations based on Bernstein multi-scaling polynomials. Appl. Math. J. China 2021, 36, 317–329. [Google Scholar] [CrossRef]
  21. Elkot, N.A.; Zaky, M.A.; Doha, E.H.; Ameen, I.G. On the rate of convergence of the Legendre spectral collocation method for multidimensional nonlinear Volterra-Fredholm integral equations. Theor. Phys. China 2021, 73, 025002. [Google Scholar]
  22. Agram, N.; Labed, S.B.; Øksendal, B.; Yakhlef, S. Singular control of stochastic Volterra integral equations. Acta Math. Sci. 2022, 42, 1003–1017. [Google Scholar] [CrossRef]
  23. Chen, X.; Wei, W.; Luo, A. Spectral distribution and numerical methods for rational eigenvalue problems. Symmetry 2022, 14, 1270. [Google Scholar] [CrossRef]
  24. Wu, N.; Zheng, W.; Gao, W. Symmetric spectral collocation method for a kind of nonlinear Volterra integral equation. Symmetry 2022, 14, 1091. [Google Scholar] [CrossRef]
  25. Tang, T.; Xu, X.; Cheng, J. On spectral methods for Volterra integral equations and the convergence analysis. J. Comput. Math. 2008, 26, 825–837. [Google Scholar]
  26. Yang, Y.; Chen, Y.; Huang, Y.; Yang, W. Convergence analysis of Legendre-collocation methods for nonlinear Volterra type integro equations. Adv. Appl. Math. Mech. 2015, 7, 74–88. [Google Scholar] [CrossRef]
  27. Headley, V.B. A multidimensional nonlinear Gronwall inequality. J. Math. Anal. Appl. 1974, 47, 250–255. [Google Scholar] [CrossRef] [Green Version]
  28. Fedotov, A.I. Lebesgue constant estimation in multidimensional Sobolev space. J. Math. 2004, 14, 25–32. [Google Scholar]
  29. Canuto, C.; Hussaini, M.Y.; Quarteroni, A.; Zang, T.A. Spectral Methods: Fundamentals in Single Domains; Springer: Berlin, Germany, 2006. [Google Scholar]
  30. Wei, Y.; Chen, Y.; Shi, X.; Zhang, Y. Spectral method for multidimensional Volterra integral equation with regular kernel. Front. Math. China 2019, 14, 435–448. [Google Scholar] [CrossRef]
  31. Nevai, P. Mean convergence of Lagrange interpolation. Trans. Am. Math. Soc. 1984, 282, 669–698. [Google Scholar] [CrossRef]
  32. Ragozin, D.L. Polynomial approximation on compact manifolds and homogeneous spaces. Trans. Am. Math. Soc. 1970, 150, 41–53. [Google Scholar] [CrossRef]
  33. Ragozin, D.L. Constructive polynomial approximation on spheres and projective spaces. Trans. Am. Math. Soc. 1971, 162, 157–170. [Google Scholar]
  34. Wei, Y.; Chen, Y. A Jacobi spectral method of multidimensional linear Volterra integral equation of the second kind. J. Sci. Comput. 2019, 79, 1801–1813. [Google Scholar] [CrossRef]
  35. Boykov, I.V.; Tynda, A.N. Numerical methods of optimal accuracy for weakly singular Volterra integral equations. Ann. Funct. Anal. 2015, 6, 114–133. [Google Scholar] [CrossRef]
Figure 1. The errors in L and L ω α , β 2 norms for Example 1.
Figure 1. The errors in L and L ω α , β 2 norms for Example 1.
Symmetry 15 00060 g001
Figure 2. The errors in L and L ω α , β 2 norms for Example 2.
Figure 2. The errors in L and L ω α , β 2 norms for Example 2.
Symmetry 15 00060 g002
Table 1. The infinity error and weighted square error for Example 1.
Table 1. The infinity error and weighted square error for Example 1.
N2468
L -error0.8907540.8131497.16 × 10 7 4.65 × 10 10
L ω α , β 2 -error1.0643460.8881065.91 × 10 7 5.31 × 10 10
N10121416
L -error1.43 × 10 11 7.40 × 10 13 7.82 × 10 13 7.35 × 10 13
L ω α , β 2 -error1.67 × 10 11 6.63 × 10 13 6.92 × 10 13 6.45 × 10 13
Table 2. The infinity error and weighted square error for Example 2.
Table 2. The infinity error and weighted square error for Example 2.
N2468
L -error1.1268581.8473363.18 × 10 3 1.01 × 10 5
L ω α , β 2 -error1.7771182.4216844.59 × 10 3 7.04 × 10 6
N10121416
L -error1.34 × 10 7 4.30 × 10 10 1.85 × 10 12 1.34 × 10 12
L ω α , β 2 -error5.19 × 10 8 8.30 × 10 11 2.02 × 10 12 2.84 × 10 12
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

Wu, D.; Zheng, W.; Chen, Y. Error Estimates of a Symmetric Spectral Method for a Linear Volterra Integral Equation. Symmetry 2023, 15, 60. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15010060

AMA Style

Wu D, Zheng W, Chen Y. Error Estimates of a Symmetric Spectral Method for a Linear Volterra Integral Equation. Symmetry. 2023; 15(1):60. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15010060

Chicago/Turabian Style

Wu, Danna, Weishan Zheng, and Yanfeng Chen. 2023. "Error Estimates of a Symmetric Spectral Method for a Linear Volterra Integral Equation" Symmetry 15, no. 1: 60. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15010060

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