Next Article in Journal
Degenerate Fubini-Type Polynomials and Numbers, Degenerate Apostol–Bernoulli Polynomials and Numbers, and Degenerate Apostol–Euler Polynomials and Numbers
Next Article in Special Issue
A Space-Time Legendre-Petrov-Galerkin Method for Third-Order Differential Equations
Previous Article in Journal
Stability for a Class of Differential Set-Valued Inverse Variational Inequalities in Finite Dimensional Spaces
Previous Article in Special Issue
Stability of a Nonlinear ML-Nonsingular Kernel Fractional Langevin System with Distributed Lags and Integral Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A High-Order Approximate Solution for the Nonlinear 3D Volterra Integral Equations with Uniform Accuracy

School of Data Science and Information Engineering, Guizhou Minzu University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(9), 476; https://doi.org/10.3390/axioms11090476
Submission received: 9 August 2022 / Revised: 13 September 2022 / Accepted: 14 September 2022 / Published: 16 September 2022
(This article belongs to the Special Issue Numerical Methods for Fractional and Integer PDEs)

Abstract

:
In this paper, we present a high-order approximate solution with uniform accuracy for nonlinear 3D Volterra integral equations. This numerical scheme is constructed based on the three-dimensional block cubic Lagrangian interpolation method. At the same time, we give the local truncation error analysis of the numerical scheme based on Taylor’s theorem. Through theoretical analysis, we reach the conclusion that the optimal convergence order of this high-order numerical scheme is 4. Finally, we verify the effectiveness and applicability of the method through four numerical examples.

1. Introduction

Nonlinear Volterra integral equations (VIEs) are widely used in many practical problems, such as electromagnetism and plasmon progeny physics, hydrodynamics, oscillation theory, polymer rheology, chemical dynamics, biomechanics, and control theory, etc. [1]. However, nonlinear VIEs are equations with nonlinear and integral terms, and only a small part of them can obtain accurate solutions directly. In the process of rapid development in mathematical physics, engineering, medicine, and other fields, more and more researchers have begun to study the numerical calculation methods of nonlinear VIEs [2].
In this paper, we will consider the following nonlinear 3D-VIEs of the form
μ ( x , y , z ) = g ( x , y , z ) + 0 x 0 y 0 z K ( x , y , z , w , t , ν , μ ( w , t , ν ) ) d ν d t d w ,
where ( x , y , z ) Θ , g ( x , y , z ) and K ( x , y , z , w , t , ν , μ ( w , t , ν ) ) are given continuous functions defined on Θ = [ 0 , g ¯ ] 3 , Φ = Θ × Θ × ( , + ) , respectively, and μ ( x , y , z ) is unknown on Θ .
The rigid-body collisions with friction can be described by a linear VIE considering the impulses and sliding velocities as functions of the direction of the sliding velocity in [3]. The numerical solution for solving the system of VIEs without special starting procedures was given by the block-by-block method (BLBM) in [4]. In [5], the authors give a numerical solution for the second kind of VIE under rigorous error analysis by spectral methods. In [6], the second kind of VIE, with smooth or weakly singular kernels, was solved by the Taylor-series expansion method. In [7], the Laplace transform can be used to solve VIEs by using a new transformation. In [8], the Urysohn-type nonlinear VIEs were solved by the Euler and trapezoidal discretization methods. An ordinary method for constructing the numerical schemes of the fractional ordinary differential equations (FODEs) was introduced by using the idea of BLBM without solving the coupling unknown solutions at each block step in [9]. In [10], they introduced a high-order numerical scheme for FODEs with the Caputo derivative. This numerical scheme is implemented by dividing the interval into a great quantity of subintervals and using quadratic interpolation on them. At the same time, it is shown that the numerical scheme is proved to be unconditionally stable for general nonlinear equations with the uniform sharp convergence order. In [11], the authors propose a high-order numerical scheme for an integro-differential equation with a fractional Caputo derivative based on the shifted Legendre polynomials by coupling the idea of the Gauss–Legendre quadrature rule and spectral collocation method. The second kind of VIE was solved by the Bernstein’s approximation in [12]. Based on the row and column sweep parallel algorithms, a parallel solution of the second kind of linear VIE was given with the suitable quadrature technique in [13]. Some efficient numerical solutions for the linear system of Volterra–Fredholm integral equations (VFIEs) were constructed based on the collocation method with the help of the Bernstein polynomial in [14]. In [15], the numerical solution for the separable kernels VIEs were given by the method of differential transform. By the radial basis functions with scattered points, the VIEs were solved with the analysis of the error estimation and the convergence rate estimation in [16]. In [17], the authors studied the numerical solution of the nonlinear 2D VIE by using the collocation method and iterated collocation method. In [18], the Euler-type numerical solution for the first kind 2D-VIEs was introduced with its convergence analysis. It gave a numerical scheme for solving the 2D differential transform for double integrals in [19]. The numerical solution for 2D nonlinear fractional VFIEs was constructed based on the collocation method by using the block-pulse functions of 2D and the shifted Legendre polynomials of two variables in [20]. In [21], the second 2D nonlinear VIE was solved by the Galerkin method with the help of the moving least squares method. The paper proposed a numerical method to solve 2D-VIEs with fractional order, weakly singular kernels based on 2D Euler polynomials combined with the Gauss–Jacobi quadrature formula in [22]. In [23], an efficient method was presented for solving the second kind of 3D-VFIEs based on 3D Bernstein polynomials. The spectrally accurate collacation method was introduced to solve the second integral equation of a weakly singular kernel by using multivarate Jacobi approximating in [24]. In [25], the results of 2D-VIE was extended to solve 3D-VIE by using the reduced differential transform method. The dimensional differential transform method was applied for solving nonlinear 3D-VIEs in [26]. The authors present an optimal homotopy asymptotic method for solving 3D-VIEs of the second kind in [27]. A new numerical method for solving 3D VFIEs was presented based upon 3D block-pulse functions approximation in [28].
The study of a numerical shceme for nonlinear 3D-VIEs is an essential research topic. Thus far, scholars have limited research on higher-order numerical schemes for nonlinear 3D-VIEs with uniform accuracy, which further strengthens the significance of our research on this topic. In this article, we divide the region into many subdomains and construct a new numerical scheme for nonlinear 3D-VIE by using three-dimensional block cubic Lagrangian interpolation. The scheme has a uniform accuracy and an optimal convergence order of 4. At the same time, we can obtain the convergence analysis of the numerical scheme based on the numerical analysis method, the Granwall inequality and the idea of [18].
The arrangement of the paper is as follows. The existence and uniqueness of the analytical solution is proposed based on the compressed mapping in Section 2. High-order approximate solutions are constructed based on the idea of Lagrangian interpolation methods in Section 3. In Section 4, the local truncation error of the constructed high-order approximate solution is given based on Taylor’s theorem. In Section 5, the convergence analysis of the high-order approximate solution is studied by using the Gronwall inequality. Four numerical examples are given to verify the correctness of convergence theoretical analysis in Section 6. In Section 7, some conclusions are given to discuss the computational efficiency of the high-order approximate solution, and further studies are discussed.

2. Existence and Uniqueness of Solution

For our numerical solution to be meaningful, we will first prove that Equation (1) has a solution and only one solution. Given this condition, a numerical scheme constructed based on Equation (1) makes sense. In this section, to prove this condition, we will establish the existence and uniqueness of the solution in Equation (1) based on the compressive mapping method, and give the corresponding theorem as follows. In this theorem, we define the following symbols: X = ( x , y , z ) and W = ( w , t , ν ) for narrative convenience.
Theorem 1.
Let K ( X , W , μ ) satisfy the following condition with respect to variable μ,
| K ( X , W , μ 1 ) K ( X , W , μ 2 ) | L | μ 1 μ 2 | , L > 0 ,
and then Equation (1) has a unique solution in Θ = [ 0 , g ¯ ] 3 .
Proof. 
Set φ = max X Θ e ( x + y + z ) φ ( X ) 2 , where · 2 is the L 2 norm. Suppose S is a continuous space of function on [ 0 , g ¯ ] 3 R 3 is defined as S = { ψ ( X ) : | | ψ ( X ) | | < + , ψ ( X ) , X [ 0 , g ¯ ] 3 } . Let A : S S be an operator as follows:
A ( φ ) ( X ) = g ( X ) + 0 x 0 y 0 z K ( X , W , φ ( W ) ) d ν d t d w .
In the following, we will prove that A is a compressed mapping. For φ , ψ S , by a directly calculate we can obtain
A ( φ ) A ( ψ ) = max X Θ e ( x + y + z ) 0 x 0 y 0 z K ( X , W , φ ( W ) ) d ν d t d w 0 x 0 y 0 z K ( X , W , ψ ( W ) ) d ν d t d w 2 max X Θ e ( x + y + z ) 0 x 0 y 0 z L e ( w + t + ν ) e ( w + t + ν ) φ ( W ) ψ ( W ) 2 d ν d t d w max X Θ e ( x + y + z ) 0 x 0 y 0 z L e ( w + t + ν ) d ν d t d w φ ψ max X Θ L ( 1 e x ) ( 1 e y ) ( 1 e z ) 3 φ ψ L 3 φ ψ .
Choosing > L 1 / 3 > 0 , it is easy to show that A is a compressed mapping. Therefore, according to the principle of a compressed fixed point, there is a unique solution μ S that satisfies Equation (1). The proof is then completed. □

3. Construction of the High-Order Approximate Solution

Now we construct the higher-order approximate solution of (1). Let the domain [ 0 , g ¯ ] 3 be divided into N 3 equal subdomains, and set x i = y i = z i = i h , i = 0 , 1 , , N , here h = g ¯ N . In the following, denoting by μ i , j q is the numerical solution at point ( x i , y j , z q ) . Meanwhile, according to this definition and formula (1), we have μ i , 0 0 = g ( x i , 0 , 0 ) , μ 0 , j 0 = g ( 0 , y j , 0 ) , μ 0 , 0 q = g ( 0 , 0 , z q ) , μ i , j 0 = g ( x i , y j , 0 ) , μ i , 0 q = g ( x i , 0 , z q ) and μ 0 , j q = g ( 0 , y j , z q ) , i , j , q = 0 , 1 , , N . For the sake of simplicity, we set K ( x i , y j , z q , w , t , ν , μ ( w , t , ν ) ) = K i , j , q ( w , t , ν , μ ( w , t , ν ) ) and g ( x i , y j , z q ) = g i , j , q in the following parts.
Let ψ i , j ( s ) , i = 0 , 1 , 2 , 3 ; j N be the four basic functions of cubic interpolation at points s j , s j + 1 , s j + 2 , s j + 3 , respectively, and s j = j h , the definition of ψ i , j ( s ) is
ψ 0 , j ( s ) = ( s ( j + 1 ) h ) ( s ( j + 2 ) h ) ( s ( j + 3 ) h ) 6 h 3 , ψ 1 , j ( s ) = ( s j h ) ( s ( j + 2 ) h ) ( s ( j + 3 ) h ) 2 h 3 , ψ 2 , j ( s ) = ( s j h ) ( s ( j + 1 ) h ) ( s ( j + 3 ) h ) 2 h 3 , ψ 3 , j ( s ) = ( s j h ) ( s ( j + 1 ) h ) ( s ( j + 2 ) h ) 6 h 3 .
First, we construct an approximate solution of μ ( x 1 , y 1 , z 1 ) , which has
μ ( x 1 , y 1 , z 1 ) = g 1 , 1 , 1 + 0 x 1 0 y 1 0 z 1 K 1 , 1 , 1 ( w , t , ν , μ ( w , t , ν ) ) d ν d t d w g 1 , 1 , 1 + 0 x 1 0 y 1 0 z 1 a , b , c = 0 3 ψ a , 0 ( w ) ψ b , 0 ( t ) ψ c , 0 ( ν ) K 1 , 1 , 1 ( x a , y b , z c , μ ( x a , y b , z c ) ) d ν d t d w = g 1 , 1 , 1 + a , b , c = 0 3 λ a , 0 1 λ ˜ b , 0 1 λ ¯ c , 0 1 K 1 , 1 , 1 ( x a , y b , z c , μ ( x a , y b , z c ) ) ,
where
λ a , 0 1 = 0 x 1 ψ a , 0 ( w ) d w , λ ˜ b , 0 1 = 0 y 1 ψ b , 0 ( t ) d t , λ ¯ c , 0 1 = 0 z 1 ψ c , 0 ( ν ) d ν , a , b , c = 0 , 1 , 2 , 3 .
Secondly, we compute an approximate solution at point ( x 2 , y 1 , z 1 ) .
μ ( x 2 , y 1 , z 1 ) = g 2 , 1 , 1 + 0 x 2 0 y 1 0 z 1 K 2 , 1 , 1 ( w , t , ν , μ ( w , t , ν ) ) d ν d t d w g 2 , 1 , 1 + 0 x 2 0 y 1 0 z 1 a , b , c = 0 3 ψ a , 0 ( w ) ψ b , 0 ( t ) ψ c , 0 ( ν ) K 2 , 1 , 1 ( x a , y b , z c , μ ( x a , y b , z c ) ) d ν d t d w = g 2 , 1 , 1 + a , b , c = 0 3 λ a , 0 2 λ ˜ b , 0 1 λ ¯ c , 0 1 K 2 , 1 , 1 ( x a , y b , z c , μ ( x a , y b , z c ) ) ,
with
λ a , 0 2 = 0 x 2 ψ a , 0 ( w ) d w , a = 0 , 1 , 2 , 3 .
Because the approximate solution can be easily obtained by using the calculation similar to (4) and (5) for the other cases, the calculation process is omitted here. Therefore, the approximate solution for μ ( x k , y l , z m ) ,   k , l , m = 1 , 2 , 3 is given as follows:
μ ( x k , y l , z m ) g k , l , m + a , b , c = 0 3 λ a , 0 k λ ˜ b , 0 l λ ¯ c , 0 m K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) , k , l , m = 1 , 2 , 3 ,
where
λ a , 0 k = 0 x k ψ a , 0 ( w ) d w , λ ˜ b , 0 l = 0 y l ψ b , 0 ( t ) d t , λ ¯ c , 0 m = 0 z m ψ c , 0 ( ν ) d ν , a , b , c = 0 , 1 , 2 , 3 ; k , l , m = 1 , 2 , 3 .
From (6), we can get the high-order numerical scheme for μ ( x k , y l , z m ) , k , l , m = 1 , 2 , 3 as follows:
μ k , l m = g k , l , m + a , b , c = 0 3 λ a , 0 k λ ˜ b , 0 l λ ¯ c , 0 m K k , l , m ( x a , y b , z c , μ a , b c ) , k , l , m = 1 , 2 , 3 .
Finally, for μ ( x k , y l , z m ) , k , l , m 4 , we can obtain an approximate solution of the form
μ ( x k , y l , z m ) = g k , l , m + 0 x k 0 y l 0 z m K k , l , m ( w , t , ν , μ ( w , t , ν ) ) d ν d t d w = g k , l , m + 0 x 3 0 y 3 0 z 3 K k , l , m ( w , t , ν , μ ( w , t , ν ) ) d ν d t d w + d = 4 k e = 4 l f = 4 m x d 1 x d y e 1 y e z f 1 z f K k , l , m ( w , t , ν , μ ( w , t , ν ) ) d ν d t d w g k , l , m + a , b , c = 0 3 0 x 3 0 y 3 0 z 3 ψ a , 0 ( w ) ψ b , 0 ( t ) ψ c , 0 ( ν ) K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) d ν d t d w + d = 4 k e = 4 l f = 4 m a , b , c = 0 3 x d 1 x d y e 1 y e z f 1 z f ψ a , d 3 ( w ) ψ b , e 3 ( t ) ψ c , f 3 ( ν ) d ν d t d w × K k , l , m ( x d 3 + a , y e 3 + b , z f 3 + c , μ ( x d 3 + a , y e 3 + b , z f 3 + c ) ) = g k , l , m + a , b , c = 0 3 λ a , 0 3 λ ˜ b , 0 3 λ ¯ c , 0 3 K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) + d = 4 k e = 4 l f = 4 m a , b , c = 0 3 τ a , d 3 d τ ˜ b , e 3 e τ ¯ c , f 3 f K k , l , m ( x d 3 + a , y e 3 + b , z f 3 + c , μ ( x d 3 + a , y e 3 + b , z f 3 + c ) ) ,
where ψ a , 0 ( w ) , ψ a , d 3 ( w ) , ψ b , 0 ( t ) , ψ b , e 3 ( t ) , ψ c , 0 ( ν ) , ψ c , f 3 ( ν ) , and λ a , 0 3 , λ ˜ b , 0 3 , λ ¯ c , 0 3 are defined in (3) and (7), respectively, and
τ a , d 3 d = x d 1 x d ψ a , d 3 ( w ) d w , a = 0 , 1 , 2 , 3 ; d = 4 , 5 , , k , τ ˜ b , e 3 e = y e 1 y e ψ b , e 3 ( t ) d t , b = 0 , 1 , 2 , 3 ; e = 4 , 5 , , l , τ ¯ c , f 3 f = z f 1 z f ψ c , f 3 ( ν ) d ν , c = 0 , 1 , 2 , 3 ; f = 4 , 5 , , m .
According to the estimated formula of (9), we can easily obtain the numerical scheme at point ( x k , y l , z m ) , k , l , m = 4 , 5 , , N as follows:
μ k , l m = g k , l , m + a , b , c = 0 3 λ a , 0 3 λ ˜ b , 0 3 λ ¯ c , 0 3 K k , l , m ( x a , y b , z c , μ a , b c ) + d = 4 k e = 4 l f = 4 m a , b , c = 0 3 τ a , d 3 d τ ˜ b , e 3 e τ ¯ c , f 3 f K k , l , m ( x d 3 + a , y e 3 + b , z f 3 + c , μ d 3 + a , e 3 + b f 3 + c ) .
To summarize, by combining (8) and (10), we can get the higher-order numerical scheme of problem (1) as follows:
μ k , l m = g k , l , m + a , b , c = 0 3 λ a , 0 k λ ˜ b , 0 l λ ¯ c , 0 m K k , l , m ( x a , y b , z c , μ a , b c ) , k , l , m = 1 , 2 , 3 . μ k , l m = g k , l , m + a , b , c = 0 3 λ a , 0 3 λ ˜ b , 0 3 λ ¯ c , 0 3 K k , l , m ( x a , y b , z c , μ a , b c ) + d = 4 k e = 4 l f = 4 m a , b , c = 0 3 τ a , d 3 d τ ˜ b , e 3 e τ ¯ c , f 3 f K k , l , m ( x d 3 + a , y e 3 + b , z f 3 + c , μ d 3 + a , e 3 + b f 3 + c ) , k , l , m = 4 , 5 , , N .

4. Estimation of the Truncation Error

This section mainly deals with the truncation error analysis of the numerical scheme (11) of problem (1). Here we give the following definition:
E k , l , m μ ( x k , y l , z m ) μ ¯ ( x k , y l , z m ) ,
where E k , l , m represents the truncation error at point ( x k , y l , z m ) produced by the numerical scheme proposed in Section 3, μ ( x k , y l , z m ) is the exact solution of problem (1) in this paper at point ( x k , y l , z m ) , and μ ¯ ( x k , y l , z m ) is an approximate value. We replace the numerical solution μ k , l m on the right side of Equation (11) with the value obtained after the corresponding exact solution μ ( x k , y l , z m ) , which is expressed as μ ¯ ( x k , y l , z m ) ; for example,
μ ¯ ( x k , y l , z m ) = g k , l , m + a , b , c = 0 3 λ a , 0 k λ ˜ b , 0 l λ ¯ c , 0 m K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) , k , l , m = 1 , 2 , 3 . μ ¯ ( x k , y l , z m ) = g k , l , m + a , b , c = 0 3 λ a , 0 3 λ ˜ b , 0 3 λ ¯ c , 0 3 K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) + d = 4 k e = 4 l f = 4 m a , b , c = 0 3 τ a , d 3 d τ ˜ b , e 3 e τ ¯ c , f 3 f K k , l , m ( x d 3 + a , y e 3 + b , z f 3 + c , μ ( x d 3 + a , y e 3 + b , z f 3 + c ) ) , k , l , m = 4 , 5 , , N .
Based on Taylor’s theorem, we will introduce the following lemma to estimate the truncation error of the numerical scheme (11) proposed in Section 3. For the convenience of description, we set 4 K s 4 = s 4 K .
Lemma 1.
Suppose the function K C 4 ( [ 0 , g ¯ ] 3 ) , and E k , l , m represents the truncation error defined in (12), and then it holds that
| E k , l , m | C h 4 , k , l , m = 1 , 2 , , N ,
where C is a positive constant that only depends on g ¯ , M , with
M = max [ 0 , g ¯ ] 3 ( | w 4 K ( x , y , z , w , t , ν , μ ( w , t , ν ) ) | , | t 4 K ( x , y , z , w , t , ν , μ ( w , t , ν ) ) | , | ν 4 K ( x , y , z , w , t , ν , μ ( w , t , ν ) ) | ) .
Proof. 
We first estimate E k , l , m , k , l , m = 1 , 2 , 3 and according to the definition of truncation error, we bring (1) and (13) into (12) to obtain
E k , l , m = 0 x k 0 y l 0 z m K k , l , m ( w , t , ν , μ ( w , t , ν ) ) d ν d t d w 0 x k 0 y l 0 z m a , b , c = 0 3 ψ a , 0 ( w ) ψ b , 0 ( t ) ψ c , 0 ( ν ) K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) d ν d t d w = 0 x k 0 y l 0 z m ( K k , l , m ( w , t , ν , μ ( w , t , ν ) ) a , b , c = 0 3 ψ a , 0 ( w ) ψ b , 0 ( t ) ψ c , 0 ( ν ) K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) ) d ν d t d w 0 x k 0 y l 0 z m γ 1 d ν d t d w .
According to Taylor’s theorem, we can know that for all ( w , t , ν ) [ 0 , x k ] × [ 0 , y l ] × [ 0 , z m ] , k , l , m = 1 , 2 , 3 , there exists ( ξ k ( w ) , ϑ l ( t ) , ϱ m ( ν ) ) [ 0 , x k ] × [ 0 , y l ] × [ 0 , z m ] , which satisfies
γ 1 = 1 4 ! w 4 K k , l , m ( ξ k ( w ) , t , ν , μ ( ξ k ( w ) , t , ν ) ) a = 0 3 ( w x a ) + a = 0 3 ψ a , 0 ( w ) 1 4 ! t 4 K k , l , m ( x a , ϑ l ( t ) , ν , μ ( x a , ϑ l ( t ) , ν ) ) b = 0 3 ( t y b ) + a = 0 3 b = 0 3 ψ a , 0 ( w ) ψ b , 0 ( t ) 1 4 ! ν 4 K k , l , m ( x a , y b , ϱ m ( ν ) , μ ( x a , y b , ϱ m ( ν ) ) ) c = 0 3 ( ν z c ) .
Through the definition of the interpolation basis functions ψ a , 0 ( w ) , ψ b , 0 ( t ) , with w ( 0 , x k ) , t ( 0 , y l ) , a , b = 0 , 1 , 2 , 3 ; k , l = 1 , 2 , 3 , we can know that | ψ a , 0 ( w ) | 1 , | ψ b , 0 ( t ) | 1 . Therefore, we can get the conclusion that | a = 0 3 ψ a , 0 ( w ) | 4 and | a = 0 3 b = 0 3 ψ a , 0 ( w ) ψ b , 0 ( t ) | 16 .
Combining (16) and (17) for k , l , m = 1 , we directly get
| E 1 , 1 , 1 | | 0 x 1 0 y 1 0 z 1 1 4 ! w 4 K 1 , 1 , 1 ( ξ 1 ( w ) , t , ν , μ ( ξ 1 ( w ) , t , ν ) ) a = 0 3 ( w x a ) d ν d t d w | + | 0 x 1 0 y 1 0 z 1 a = 0 3 ψ a , 0 ( w ) 1 4 ! t 4 K 1 , 1 , 1 ( x a , ϑ 1 ( t ) , ν , μ ( x a , ϑ 1 ( t ) , ν ) ) b = 0 3 ( t y b ) d ν d t d w | + | 0 x 1 0 y 1 0 z 1 a = 0 3 b = 0 3 ψ a , 0 ( w ) ψ b , 0 ( t ) 1 4 ! ν 4 K 1 , 1 , 1 ( x a , y b , ϱ 1 ( ν ) , μ ( x a , y b , ϱ 1 ( ν ) ) ) c = 0 3 ( ν z c ) d ν d t d w | γ 1 1 + γ 1 2 + γ 1 3 .
Now, we analyze the terms on the right side of Equation (18) one by one, and we get
γ 1 1 0 x 1 0 y 1 0 z 1 1 4 ! M × | a = 0 3 ( w x a ) | d ν d t d w
M h 4 0 x 1 0 y 1 0 z 1 1 d ν d t d w = M h 7 ,
γ 1 2 4 4 ! M h 4 0 x 1 0 y 1 0 z 1 1 d ν d t d w M h 7 ,
γ 1 3 16 4 ! M h 4 0 x 1 0 y 1 0 z 1 1 d ν d t d w M h 7 ,
where M is defined in (15).
Through (18)–(21), we can directly get the conclusion of
| E 1 , 1 , 1 | 3 M h 7 .
By using the same method, we can estimate the error as
| E k , l , m | 3 M h 7 , k , l , m = 1 , 2 , 3 .
Next, we will analyze the case of k , l , m 4 , combining (1), (12), and (13), we have
| E k , l , m | = | 0 x k 0 y l 0 z m K k , l , m ( w , t , ν , μ ( w , t , ν ) ) d ν d t d w a , b , c = 0 3 λ a , 0 3 λ ˜ b , 0 3 λ ¯ c , 0 3 K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) d = 4 k e = 4 l f = 4 m a , b , c = 0 3 τ a , d 3 d τ ˜ b , e 3 e τ ¯ c , f 3 f K k , l , m ( x d 3 + a , y e 3 + b , z f 3 + c , μ ( x d 3 + a , y e 3 + b , z f 3 + c ) ) | | 0 x 3 0 y 3 0 z 3 ( K k , l , m ( w , t , ν , μ ( w , t , ν ) ) a , b , c = 0 3 ψ a , 0 ( w ) ψ b , 0 ( t ) ψ c , 0 ( ν ) K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) ) d ν d t d w | + | d = 4 k e = 4 l f = 4 m x d 1 x d y e 1 y e z f 1 z f ( K k , l , m ( w , t , ν , μ ( w , t , ν ) ) a , b , c = 0 3 ψ a , d 3 ( w ) ψ b , e 3 ( t ) ψ c , f 3 ( ν ) K k , l , m ( x d 3 + a , y e 3 + b , z f 3 + c , μ ( x d 3 + a , y e 3 + b , z f 3 + c ) ) d ν d t d w | | 0 x 3 0 y 3 0 z 3 γ 2 d ν d t d w | + | d = 4 k e = 4 l f = 4 m x d 1 x d y e 1 y e z f 1 z f γ 3 d ν d t d w | .
First, we analyze the first term on the right side of (23). Similar to the estimation of (18), we have
| 0 x 3 0 y 3 0 z 3 γ 2 d ν d t d w | 3 M h 7 .
Next, we estimate the second term on the right side of (23), and for all ( w , t , ν ) [ x d 1 , x d ] × [ y e 1 , y e ] × [ z f 1 , z f ] , there are ( ξ d ( w ) , ϑ e ( t ) , ϱ f ( ν ) ) [ x d 1 , x d ] × [ y e 1 , y e ] × [ z f 1 , z f ] , which satisfies
γ 3 = 1 4 ! w 4 K k , l , m ( ξ d ( w ) , t , ν , μ ( ξ d ( w ) , t , ν ) ) a = 0 3 ( w x d 3 + a ) + a = 0 3 ψ a , d 3 ( w ) 1 4 ! t 4 K k , l , m ( x d 3 + a , ϑ e ( t ) , ν , μ ( x d 3 + a , ϑ e ( t ) , ν ) ) b = 0 3 ( t y e 3 + b ) + a = 0 3 b = 0 3 ψ a , d 3 ( w ) ψ b , e 3 ( t ) 1 4 ! ν 4 K k , l , m ( x d 3 + a , y e 3 + b , ϱ f ( ν ) , μ ( x d 3 + a , y e 3 + b , ϱ f ( ν ) ) ) c = 0 3 ( ν z f 3 + c ) .
Therefore, we obtain that
| d = 4 k e = 4 l f = 4 m x d 1 x d y e 1 y e z f 1 z f γ 3 d ν d t d w | | d = 4 k e = 4 l f = 4 m x d 1 x d y e 1 y e z f 1 z f 1 4 ! w 4 K k , l , m ( ξ d ( w ) , t , ν , μ ( ξ d ( w ) , t , ν ) ) a = 0 3 ( w x d 3 + a ) d ν d t d w | + | d = 4 k e = 4 l f = 4 m x d 1 x d y e 1 y e z f 1 z f a = 0 3 ψ a , d 3 ( w ) 1 4 ! t 4 K k , l , m ( x d 3 + a , ϑ e ( t ) , ν , μ ( x d 3 + a , ϑ e ( t ) , ν ) ) × b = 0 3 ( t y e 3 + b ) d ν d t d w | + | d = 4 k e = 4 l f = 4 m x d 1 x d y e 1 y e z f 1 z f a = 0 3 b = 0 3 ψ a , d 3 ( w ) ψ b , e 3 ( t ) × 1 4 ! ν 4 K k , l , m ( x d 3 + a , y e 3 + b , ϱ f ( ν ) , μ ( x d 3 + a , y e 3 + b , ϱ f ( ν ) ) ) × c = 0 3 ( ν z f 3 + c ) d ν d t d w | γ 3 1 + γ 3 2 + γ 3 3 .
Next, we analyze each item in (26) one by one, and we have
γ 3 1 M h 4 d = 4 k e = 4 l f = 4 m x d 1 x d y e 1 y e z f 1 z f 1 d ν d t d w
= M h 4 x 3 x k y 3 y l z 3 z m 1 d ν d t d w g ¯ 3 M h 4 ,
γ 3 2 M h 4 d = 4 k e = 4 l f = 4 m x d 1 x d y e 1 y e z f 1 z f 1 d ν d t d w g ¯ 3 M h 4 ,
γ 3 3 M h 4 d = 4 k e = 4 l f = 4 m x d 1 x d y e 1 y e z f 1 z f 1 d ν d t d w g ¯ 3 M h 4 .
Consequently, combining (23), (24) and (26)–(29), we can get the following result
| E k , l , m | 3 M h 7 + 3 g ¯ 3 M h 4 , k , l , m = 4 , 5 , , N .
From the conclusions of above (22) and (30), we can obtain the truncation error of the proposed higher-order numerical scheme (11) as follows:
| E k , l , m | C h 4 , k , l , m = 1 , 2 , , N ,
where C is a positive constant that only depends on g ¯ , M . □

5. Convergence Analysis

For the coefficients appearing in the higher-order numerical scheme (11) proposed in Section 3, we can get the following conclusions through accurate calculation.
| λ a , 0 j | C h , | λ ˜ a , 0 j | C h , | λ ¯ a , 0 j | C h , a = 0 , 1 , 2 , 3 ; j = 1 , 2 , 3 ,
| τ a , j 3 j | C h , | τ ˜ a , j 3 j | C h , | τ ¯ a , j 3 j | C h , a = 0 , 1 , 2 , 3 ; j = 4 , 5 , , N ,
where C is a constant independent of h.
Theorem 2.
Let the function K C 4 ( [ 0 , g ¯ ] 3 ) satisfy (2). Let μ ( x k , y l , z m ) and μ k , l m be the solutions of (1) and (11) at points ( x k , y l , z m ) , respectively. If h is small enough, then
| μ ( x k , y l , z m ) μ k , l m | C h 4 , k , l , m = 1 , 2 , , N ,
where C is independent of h.
Proof. 
We first define the error, let ε k , l , m = μ ( x k , y l , z m ) μ k , l m , k , l , m = 0 , 1 , , N . Combining the definitions of ε k , l , m and Equation (1), we can easily get
ε k , 0 , 0 = ε 0 , l , 0 = ε 0 , 0 , m = ε k , l , 0 = ε k , 0 , m = ε 0 , l , m = 0 , k , l , m = 0 , 1 , , N .
Now, we begin to analyze ε k , l , m , k , l , m = 1 , 2 , 3 , and ε k , l , m satisfies
ε k , l , m = μ ( x k , y l , z m ) μ k , l m = μ ( x k , y l , z m ) μ ¯ ( x k , y l , z m ) + μ ¯ ( x k , y l , z m ) μ k , l m = E k , l , m + μ ¯ ( x k , y l , z m ) μ k , l m .
According to (11), (13), Lemma 1 and (32), we have
| ε k , l , m | | E k , l , m | + | C h 3 a , b , c = 0 3 [ K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) K k , l , m ( x a , y b , z c , μ a , b c ) ] | C h 4 + C L h 3 a , b , c = 0 3 | μ ( x a , y b , z c ) μ a , b c | = C h 4 + C L h 3 a , b , c = 0 3 | ε a , b , c | , k , l , m = 1 , 2 , 3 .
Simultaneously with the inequality in (36), we easily get
| ε k , l , m | C h 4 , k , l , m = 1 , 2 , 3 .
Next, we analyze ε k , l , m , k , l , m = 4 , 5 , , N . According to the definition of ε k , l , m , combined with (32), (33), and Lemma 1, we can obtain
| ε k , l , m | | C h 3 a = 0 k b = 0 l c = 0 m ( K k , l , m ( x a , y b , z c , μ ( x a , y b , z c ) ) K k , l , m ( x a , y b , z c , μ a , b c ) ) | + | E k , l , m | C L h 3 a = 0 k b = 0 l c = 0 m | ε a , b , c | + C h 4 , k , l , m = 4 , 5 , , N .
At the same time, we can get
| ε k , l , m | C L h 3 ( a = 0 k 1 b = 0 l 1 c = 0 m 1 | ε a , b , c | + b = 0 l 1 c = 0 m 1 | ε k , b , c | + a = 0 k 1 c = 0 m 1 | ε a , l , c | + a = 0 k 1 b = 0 l 1 | ε a , b , m | + c = 0 m 1 | ε k , l , c | + b = 0 l 1 | ε k , b , m | + a = 0 k 1 | ε a , l , m | + | ε k , l , m | ) + C h 4 .
Take | | ε ^ a , b | | = max 0 c N | ε a , b , c | , a = 0 , 1 , , k ; b = 0 , 1 , 2 , l ; c = 0 , 1 , , m . Then the above inequality can be equivalent to
( 1 C L g ¯ h 2 C L h 3 ) | | ε ^ k , l | | ( C L g ¯ h 2 + C L h 3 ) a = 0 k 1 b = 0 l 1 | | ε ^ a , b | | + ( C L g ¯ h 2 + C L h 3 ) b = 0 l 1 | | ε ^ k , b | | + ( C L g ¯ h 2 + C L h 3 ) a = 0 k 1 | | ε ^ a , l | | + C h 4 .
For the small-enough h and by using the Gronwall inequality [18], we have
| | ε ^ k , l | | C h 4 ,
and then
| ε k , l , m | C h 4 , k , l , m = 4 , 5 , , N .
Combining (37) and (41), we completed the proof of Theorem 2. □

6. Numerical Examples

In this section, we use the higher-order numerical scheme (11) to solve the 3D-VIEs. The following four numerical examples are used to demonstrate the convergence order of the numerical solution. In the following numerical examples, we implement higher-order numerical scheme with h = 1 N , N = 5 , 10 , 20 , 40 , 80 and error E h = max i , j , q | μ ( x i , y j , z q ) μ i , j q | , i , j , q = 0 , 1 , , N . The convergence order is calculated by log 2 E 2 h E h . At the same time, without losing generality, we all choose the g ¯ = 1 to do the following analysis.
Example 1.
We consider the linear 3D-VIEs as follows:
μ ( x , y , z ) = g ( x , y , z ) + 0 x 0 y 0 z ( x y z w 2 t ν ) μ ( w , t , ν ) d ν d t d w , ( x , y , z ) [ 0 , 1 ] 3 ,
with the function g defined as
g ( x , y , z ) = x 2 y 5 z 3 1 175 x 6 y 8 z 6 .
It is straightforward to show that μ ( x , y , z ) = x 2 y 5 z 3 is the exact solution of the above equation.
First, we perform numerical example according to the higher-order numerical solution in Section 3. According to different step lengths h, the calculated maximum absolute error and its convergence order are shown in Table 1.
The theoretical analysis results of Theorem 2 in Section 5 show that our optimal convergence order is 4 for our proposed higher-order numerical scheme. At the same time, through Table 1, we can know that when the value of h = 1 5 , 1 10 , 1 20 , 1 40 , 1 80 gradually becomes smaller, the corresponding maximum absolute error also becomes smaller from 10 4 to 10 9 , and the reduction is extremely obvious. In this process, we can know that our convergence order is gradually approaching 4.
Next, we show the distribution of function value for the numerical solution and the exact solution. The function value distribution of N = 80 in the following Figure 1, where it is a comparison of the numerical solution and the exact solution corresponding to the 3D-VIEs. Figure 1 is the three-dimensional surface mesh map of the corresponding μ ( x , y , 1 ) = x 2 y 5 when z = 1 .
As can be seen from Figure 1, the numerical solution obtained by our calculation according to the numerical scheme is a good agreement with the exact solution of this example, and the values on the corresponding nodes are almost identical. Combining the data in Table 1 and Figure 1 in this example, we know that our proposed higher-order scheme has good convergence for linear 3D-VIE, which is consistent with the theoretical results.
Example 2.
Let us consider the following linear 3D-VIEs:
μ ( x , y , z ) = g ( x , y , z ) + 0 x 0 y 0 z sin ( x w ) sin 2 ( y + t ) cos ( 2 z + ν ) μ ( w , t , ν ) d ν d t d w ,
where ( x , y , z ) [ 0 , 1 ] 3 , and the function g is defined as
g ( x , y , z ) = x 5 y 2 z 2 x 5 20 x 3 + 120 x 120 sin ( x ) ( 1 6 y 3 1 4 y 2 sin ( 4 y ) 1 4 y cos ( 4 y ) + 1 8 sin ( 4 y ) 1 8 sin ( 2 y ) ) z 2 sin ( 3 z ) + 2 z cos ( 3 z ) 2 sin ( 3 z ) + 2 sin ( 2 z ) .
The exact solution to the above equation is as follows:
μ ( x , y , z ) = x 5 y 2 z 2 .
In this example, we implement this calculation process to calculate the maximum absolute error and convergence order of the high-order numerical solution. According to different step sizes h, the maximum absolute error, and the convergence order can be obtained as shown in Table 2.
As can be seen from Table 2, when our step size h is 1 5 , the maximum absolute error of this numerical example is already very small, and the maximum absolute error at this time is about 1.3 × 10 4 . When h = 1 80 , the maximum absolute error is about 3.6 × 10 9 . At the same time, as h gradually decreases, the maximum absolute error becomes smaller and smaller. From the definition of the maximum absolute error, we can also know that the numerical solution in this numerical problem gradually approaches the exact solution when h becomes gradually smaller. From Table 2, we can roughly get that our numerical method is gradually approaching theoretical convergence order.
In the following, we will apply two nonlinear numerical examples to further demonstrate the conclusion of our Theorem 2.
Example 3.
We consider the following 3D nonlinear VIEs:
μ ( x , y , z ) = g ( x , y , z ) + 0 x 0 y 0 z ( x y z + w t + ν 2 ) μ 2 ( w , t , ν ) d ν d t d w ; ( x , y , z ) [ 0 , 1 ] 3 ,
with the function g defined as
g ( x , y , z ) = x 3 y 5 z 3 1 539 x 8 y 12 z 8 1 672 x 8 y 12 z 7 1 693 x 7 y 11 z 9 .
It is easy to show that μ ( x , y , z ) = x 3 y 5 z 3 is the corresponding exact solution of the above equation.
In order to show the theoretical analysis results of Theorem 2 for the 3D nonlinear VIE with a general kernel function, we choose a suitable h to indicate the corresponding maximum absolute error, and the maximum absolute error gradually decreases from 10 4 to 10 8 . Among them, the step size h is the same as the selection of h in Example 1 and Example 2, and it gradually decreases from 1 5 to 1 80 . Through some of numerical calculations, we obtain the relevant data of the 3D nonlinear VIEs as shown in the following Table 3.
From Table 3, we have that the approximate solution of 3D nonlinear VIEs is convergent to the exact solution. Through a large number of tests, we can see the corresponding convergence order gradually approaches 4; that is to say, its optimal convergence order is 4. This is completely consistent with the theoretical analysis.
Next, we show the distribution of function value for the numerical solution and the exact solution. The function value distribution of N = 80 in the following Figure 2, where it is a comparison of the numerical solution and the exact solution corresponding to the 3D-VIEs. Figure 2 shows the three-dimensional surface mesh map of the corresponding μ ( x , y , 1 ) = x 3 y 5 when z = 1 .
In Figure 2, we can hardly see the difference between the exact solution and the numerical solution of this numerical example. Therefore, we can draw the conclusion that the numerical solution of the nonlinear 3D-VIEs is very consistent with the exact solution, which also verifies the correctness of our theoretical analysis.
Example 4.
Let us consider the following 3D nonlinear VIEs with exponential form:
μ ( x , y , z ) = g ( x , y , z ) + 0 x 0 y 0 z e ( x w ) + ( y t ) + ( z ν ) μ 2 ( w , t , ν ) d ν d t d w ; ( x , y , z ) [ 0 , 1 ] 3 ,
with the function g defined as
g ( x , y , z ) = x 3 y 4 z 4 ( x 6 6 x 5 30 x 4 120 x 3 360 x 2 720 x 720 + 720 e x ) ( y 8 8 y 7 56 y 6 336 y 5 1680 y 4 6720 y 3 20160 y 2 40320 y 40320 + 40320 e y ) ( z 8 8 z 7 56 z 6 336 z 5 1680 z 4 6720 z 3 20160 z 2 40320 z 40320 + 40320 e z ) ,
and the exact solution is μ ( x , y , z ) = x 3 y 4 z 4 .
In order to show the constructed numerical scheme’s good convergence, we will verify this example. The error E h and the convergence order of the solution of this 3D-VIE with respect to different h under the maximum norm of error are given in the following Table 4.
It can be seen from Table 4 that the corresponding error is gradually decreasing with respect to h, and the convergence order of numerical scheme is 4.
Summarizing the linear and nonlinear 3D-VIEs described above, we know that the numerical method is convergent with the optimal convergence order of 4. They are completely consistent with the theoretical analysis results of Theorem 2.

7. Conclusions

A higher-order approximate solution for the nonlinear three-dimensional Volterra integral equation is presented in this paper. This numerical scheme is constructed by dividing the region into multiple subregions and applying cubic Lagrangian interpolation on each subregion. The scheme has uniform accuracy with an optimal convergence order of 4. We analyze the error generated by this numerical scheme through detailed steps, and use four numerical examples to verify that the error is within a reasonable range. Then it is concluded that this theoretical analysis is correct. In the future, we hope to construct a high-order approximate solution for a peridynamics plate model based on the idea of [29]. In addition, we intend to construct an efficient high-order numerical solution with uniform accuracy for 3D-VIEs with a generally weak nonlinear singular kernel function based on the ideas of [30,31]. Finally, we will apply a fast algorithm to implement the high-order numerical scheme for large-scale practical engineering problems based on the idea of [32].

Author Contributions

Conceptualization, Z.-Q.W.; methodology, Z.-Q.W. and J.-Y.C.; software, M.-D.L.; validation, M.-D.L.; formal analysis, Z.-Q.W. and J.-Y.C.; investigation, Z.-Q.W. and M.-D.L.; resources, Z.-Q.W. and J.-Y.C.; data curation, M.-D.L.; writing—original draft preparation, Z.-Q.W. and M.-D.L.; writing—review and editing, Z.-Q.W., M.-D.L. and J.-Y.C.; visualization, M.-D.L.; supervision, Z.-Q.W. and J.-Y.C.; project administration, Z.-Q.W. and J.-Y.C.; funding acquisition, Z.-Q.W. and J.-Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Z.-Q.W. of National Natural Science Foundation of China (NSFC) grant number 11961009. This research was funded by J.-Y.C. of National Natural Science Foundation of China (NSFC) grant number 11901135 and Guizhou Provincial Science and Technology Projects grant number [2020]1Y015. This research was funded by Z.-Q.W. and J.-Y.C. of Natural Science Research Project of Department of Education of Guizhou Province (Grant Nos. QJJ2022015 and QJJ2022047).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data were computed using our high-order numerical scheme.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses or interpretation of data, in the writing of the manuscript or in the decision to publish the results.

References

  1. El-Borai, M.; Abdou, M.; Badr, A.; Basseem, M. Singular nonlinear integral equation and its application in viscoelastic nonlinear material. Int. J. Appl. Math. Mech. 2008, 4, 56–77. [Google Scholar]
  2. Brunner, H. Volterra Integral Equations: An Introduction to Theory and Applications; Cambridge University Press: New York, NY, USA, 2017. [Google Scholar]
  3. Souchet, R. An analysis of three-dimensional rigid body collisions with friction by means of a linear integral equation of Volterra. Internat. J. Engrg. Sci. 1999, 37, 365–378. [Google Scholar] [CrossRef]
  4. Katani, R.; Shahmorad, S. Block by block method for the systems of nonlinear Volterra integral equations. Appl. Math. Model. 2009, 34, 400–406. [Google Scholar] [CrossRef]
  5. Tao, 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]
  6. Maleknejad, K.; Aghazadeh, N. Numerical solution of Volterra integral equations of the second kind with convolution kernel by using Taylor-series expansion method. Appl. Math. Comput. 2003, 161, 915–922. [Google Scholar] [CrossRef]
  7. Jaabar, S.; Hussain, A. Solving Volterra integral equation by using a new transformation. J. Interdiscip. Math. 2021, 24, 735–741. [Google Scholar] [CrossRef]
  8. Bazm, S.; Lima, P.; Nemati, S. Analysis of the Euler and trapezoidal discretization methods for the numerical solution of nonlinear functional Volterra integral equations of Urysohn type. J. Comput. Appl. Math. 2021, 398, 113628. [Google Scholar] [CrossRef]
  9. Cao, J.; Xu, C. A high order schema for the numerical solution of the fractional ordinary differential equations. J. Comput. Phys. 2013, 238, 154–168. [Google Scholar] [CrossRef]
  10. Cao, J.; Cai, Z. Numerical analysis of a high-order scheme for nonlinear fractional difffferential equations with uniform accuracy. Numer. Math. Theor. Meth. Appl. 2021, 14, 71–112. [Google Scholar] [CrossRef]
  11. Wu, C.; Wang, Z. The spectral collocation method for solving a fractional integro-differential equation. AIMS Math. 2021, 6, 9577–9587. [Google Scholar] [CrossRef]
  12. Maleknejad, K.; Hashemizadeh, E.; Ezzati, R. A new approach to the numerical solution of Volterra integral equations by using Bernstein’s approximation. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 647–655. [Google Scholar] [CrossRef]
  13. Rahma, A.; Evans, D. The numerical solution of Volterra integral equations on parallel computers. Int. J. Comput. Math. 2007, 27, 103–111. [Google Scholar] [CrossRef]
  14. Maleknejad, K.; Basirat, B.; Hashemizadeh, E. A Bernstein operational matrix approach for solving a system of high order linear Volterra-Fredholm integro-differential equations. Math. Comput. Model. 2012, 55, 1363–1372. [Google Scholar] [CrossRef]
  15. Odibat, Z. Differential transform method for solving Volterra integral equation with separable kernels. Math. Comput. Model. 2008, 48, 1144–1149. [Google Scholar] [CrossRef]
  16. Assari, P.; Dehghan, M. The approximate solution of nonlinear Volterra integral equations of the second kind using radial basis functions. Appl. Numer. Math. 2018, 131, 140–157. [Google Scholar] [CrossRef]
  17. Han, G.; Hayami, K.; Sugihara, K.; Wang, J. Extrapolation method of iterated collocation solution for two-dimensional nonlinear Volterra integral equations. Appl. Math. Comput. 2000, 112, 49–61. [Google Scholar]
  18. McKee, S.; Tang, T.; Diogo, T. An Euler-type method for two-dimensional Volterra integral equations of the first kind. IMA J. Numer. Anal. 2000, 20, 423–440. [Google Scholar] [CrossRef]
  19. Tari, A.; Rahimi, M.; Shahmorad, S.; Talati, F. Solving a class of two-dimensional linear and nonlinear Volterra integral equations by the differential transform method. J. Comput. Appl. Math. 2009, 228, 70–76. [Google Scholar] [CrossRef]
  20. Maleknejad, K.; Rashidinia, J.; Eftekhari, T. Existence, uniqueness, and numerical solutions for two-dimensional nonlinear fractional Volterra and Fredholm integral equations in a Banach space. Comput. Appl. Math. 2020, 39, 1–22. [Google Scholar] [CrossRef]
  21. Assari, P.; Dehghan, M. A meshless local discrete Galerkin (MLDG) scheme for numerically solving two-dimensional nonlinear Volterra integral equations. Appl. Math. Comput. 2019, 350, 249–265. [Google Scholar] [CrossRef]
  22. Wang, Y.; Huang, J.; Wen, X. Two-dimensional Euler polynomials solutions of two-dimensional Volterra integral equations of fractional order. Appl. Numer. Math. 2021, 163, 77–95. [Google Scholar] [CrossRef]
  23. Khosrow, M.; Jalil, R.; Tahereh, E. Numerical solution of three-dimensional Volterra-Fredholm integral equations of the first and second kinds based on Bernstein’s approximation. Appl. Math. Comput. 2018, 339, 272–285. [Google Scholar]
  24. Zaky, M.; Ameen, I. A novel Jacob spectral method for multi-dimensional weakly singular nonlinear Volterra integral equations with nonsmooth solutions. Eng. Comput. 2021, 37, 2623–2631. [Google Scholar] [CrossRef]
  25. Ziqan, A.; Armiti, S.; Suwan, I. Solving three-dimensional Volterra integral equation by the reduced differential transform method. Int. J. Appl. Math. Res. 2016, 5, 103–106. [Google Scholar] [CrossRef] [Green Version]
  26. Bakhshi, M. Three-dimensional differential transform method for solving nonlinear three-dimensional Volterra integral equations. Int. J. Appl. Math. Comput. Sci. 2012, 4, 246–256. [Google Scholar] [CrossRef]
  27. Nawaz, R.; Ahsan, S.; Akbar, M.; Farooq, M.; Sulaiman, M.; Ullah, H.; Islam, S. Semi analytical solutions of second type of three-dimensional Volterra integral equations. Int. J. Appl. Comput. Math. 2020, 6, 3079–3096. [Google Scholar] [CrossRef]
  28. Mirzaee, F.; Hadadiyan, E.; Bimesl, S. Numerical solution for three-dimensional nonlinear mixed Volterra-Fredholm integral equations via three-dimensional block-pulse functions. Appl. Math. Comput. 2014, 237, 168–175. [Google Scholar] [CrossRef]
  29. Wang, Z.; Cui, J. Second-order two-scale method for bending behavior analysis of composite plate with 3-D periodic configuration and its approximation. Sci. China Math. 2014, 57, 1713–1732. [Google Scholar] [CrossRef]
  30. Wang, Z.; Liu, Q.; Cao, J. A higher-order numerical scheme for two-dimensional nonlinear fractional Volterra integral equations with uniform accuracy. Fractal Fract. 2022, 6, 314. [Google Scholar] [CrossRef]
  31. Cao, J.; Zhang, J.; Yang, X. Fully-discrete Spectral-Galerkin scheme with second-order time-accuracy and unconditionally energy stability for the volume-conserved phase-field lipid vesicle model. J. Comput. Appl. Math. 2022, 406, 113988. [Google Scholar] [CrossRef]
  32. Gu, X.; Wu, S. A parallel-in-time iterative algorithm for Volterra partial integro-differential problems with weakly singular kernel. J. Comput. Phys. 2022, 417, 109576. [Google Scholar] [CrossRef]
Figure 1. Comparison of numerical and exact solutions.
Figure 1. Comparison of numerical and exact solutions.
Axioms 11 00476 g001
Figure 2. Comparison of numerical and exact solutions.
Figure 2. Comparison of numerical and exact solutions.
Axioms 11 00476 g002
Table 1. Maximum errors and convergence rate as functions with h.
Table 1. Maximum errors and convergence rate as functions with h.
h Error Rate
1 5 1.629551910207194 × 10 4
1 10 1.301070087844636 × 10 5 3.646704724239447
1 20 9.126369742151752 × 10 7 3.833513767270810
1 40 6.033795529880592 × 10 8 3.918903390618973
1 80 3.877251630868273 × 10 9 3.959959534894697
Table 2. Maximum errors and convergence rate as functions with h.
Table 2. Maximum errors and convergence rate as functions with h.
h Error Rate
1 5 1.306204365916219 × 10 4
1 10 1.122759745353363 × 10 5 3.540259485877905
1 20 8.162257019161956 × 10 7 3.781937294325425
1 40 5.490936638707922 × 10 8 3.893843971083792
1 80 3.558625838451235 × 10 9 3.947660105874005
Table 3. Maximum errors and decay rate as functions with h.
Table 3. Maximum errors and decay rate as functions with h.
h Error Rate
1 5 9.194168478405818 × 10 4
1 10 8.959147119291799 × 10 5 3.359285798830145
1 20 7.183001882626883 × 10 6 3.640702600341250
1 40 5.105040028485774 × 10 7 3.814592715250178
1 80 3.399614900700954 × 10 8 3.908479037192826
Table 4. Maximum errors and decay rate as functions with h.
Table 4. Maximum errors and decay rate as functions with h.
h Error Rate
1 5 2.116544808936638 × 10 4
1 10 1.715380950773060 × 10 5 3.625110122038858
1 20 1.221438151288723 × 10 6 3.811876286531620
1 40 8.137970008981199 × 10 8 3.907768039330447
1 80 5.251625134761184 × 10 9 3.953833119346601
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.-Q.; Long, M.-D.; Cao, J.-Y. A High-Order Approximate Solution for the Nonlinear 3D Volterra Integral Equations with Uniform Accuracy. Axioms 2022, 11, 476. https://doi.org/10.3390/axioms11090476

AMA Style

Wang Z-Q, Long M-D, Cao J-Y. A High-Order Approximate Solution for the Nonlinear 3D Volterra Integral Equations with Uniform Accuracy. Axioms. 2022; 11(9):476. https://doi.org/10.3390/axioms11090476

Chicago/Turabian Style

Wang, Zi-Qiang, Ming-Dan Long, and Jun-Ying Cao. 2022. "A High-Order Approximate Solution for the Nonlinear 3D Volterra Integral Equations with Uniform Accuracy" Axioms 11, no. 9: 476. https://doi.org/10.3390/axioms11090476

APA Style

Wang, Z. -Q., Long, M. -D., & Cao, J. -Y. (2022). A High-Order Approximate Solution for the Nonlinear 3D Volterra Integral Equations with Uniform Accuracy. Axioms, 11(9), 476. https://doi.org/10.3390/axioms11090476

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