Next Article in Journal
Monotone Iterative and Upper–Lower Solution Techniques for Solving the Nonlinear ψ−Caputo Fractional Boundary Value Problem
Previous Article in Journal
A Special Study of the Mixed Weighted Fractional Brownian Motion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Landweber Iterative Regularization Method for Identifying the Initial Value Problem of the Rayleigh–Stokes Equation

School of Science, Lanzhou University of Technology, Lanzhou 730050, China
*
Author to whom correspondence should be addressed.
The main idea of this paper was proposed by Dun-Gang Li, Jun-Liang Fu and Fan Yang prepared the manuscript initially and performed all the steps of the proofs in this research. All authors read and approved the final manuscript.
Fractal Fract. 2021, 5(4), 193; https://doi.org/10.3390/fractalfract5040193
Submission received: 27 August 2021 / Revised: 22 October 2021 / Accepted: 27 October 2021 / Published: 2 November 2021

Abstract

:
In this paper, we study an inverse problem to identify the initial value problem of the homogeneous Rayleigh–Stokes equation for a generalized second-grade fluid with the Riemann–Liouville fractional derivative model. This problem is ill posed; that is, the solution (if it exists) does not depend continuously on the data. We use the Landweber iterative regularization method to solve the inverse problem. Based on a conditional stability result, the convergent error estimates between the exact solution and the regularization solution by using an a priori regularization parameter choice rule and an a posteriori regularization parameter choice rule are given. Some numerical experiments are performed to illustrate the effectiveness and stability of this method.

1. Introduction

Let Ω R d ( d = 1 , 2 , 3 ) be a convex polyhedral domain with its boundary being Ω , and T > 0 is a fixed time. In this paper, we consider the following homogeneous Rayleigh–Stokes problem for a generalized second-grade fluid with a fractional derivative model
t u ( x , t ) ( 1 + γ 0 R t α ) Δ u ( x , t ) = 0 , x Ω , t ( 0 , T ) , u ( x , t ) = 0 , x Ω , t ( 0 , T ) , u ( x , 0 ) = f ( x ) , x Ω , u ( x , T ) = g ( x ) , x Ω ,
where γ > 0 is a fixed constant and 0 R t α is the Riemann–Liouville fractional derivative of the order α ( 0 < α < 1 ) defined by
0 R t α u ( x , t ) = d d t 0 t ω 1 α ( t τ ) u ( x , τ ) d τ , ω α ( t ) = t α 1 Γ ( α ) ,
in which Γ ( · ) is the gamma function.
In problem (1), if f ( x ) is known, we can use some basic methods to solve the direct problem. The inverse problem in this paper is to reconstruct the initial value f ( x ) according to the additional data g ( x ) . However, in practical problems, g ( x ) can only be obtained by measurement, and the measured data will inevitably be disturbed by noise. Therefore, suppose the exact data function g ( x ) and the measured data function g δ ( x ) satisfy
g δ ( · ) g ( · ) δ ,
where · denotes L 2 ( Ω ) norm and δ > 0 is a noise level.
Recently, many scientists have demonstrated that fractional models describe natural phenomena in an accurate, systematic way better than their classic integer-order counterparts with ordinary time derivatives [1,2,3,4,5,6,7]. Recently, the fractional calculus has been employed for the description of many complex biological systems. Although these studies provided better results than the classic integer-order models, a satisfactory precision may not be achieved in the whole of the time duration because of the appearance of a singularity in the definition of traditional fractional derivatives, a fact which makes such operators impractical for the description of nonlocal dynamics. The Rayleigh–Stokes problem has attracted extensive attention in recent years due to its importance in physics. Fractional derivatives are of great value in capturing viscoelastic behavior of flows; see [8,9]. Problem (1) plays an important role in describing the behavior of some non-Newtonian fluids, such as polymer solutions and melts [10]. Regarding the direct problem of Rayleigh–Stokes, the reader is referred to [11,12,13,14,15,16,17,18,19].
However, there are few results on Rayleigh–Stokes inverse problem, and some regularization methods in particular are used to study the Rayleigh–Stokes inverse problem. For identifying the initial value problem, Nguyen et al. [20] investigated a backward problem for the Rayleigh–Stokes problem by using the filter regularization method in Gaussian random noise, with the aim of determining the initial status of some physical fields, such as the temperature for slow diffusion from its present measurement data. Furthermore, based on a priori assumptions, the expectation between the exact solution and the regularization solution under the L 2 and H m norm was established, but the authors did not give the posteriori regularization parameter choice rule. Compared with the priori regularization parameter choice, the posteriori regularization parameter choice rule, which only depends on the measurable data, may be useful in practice. Moreover, the authors did not carry out numerical experiments to show their methods.
In this paper, we mainly use the Landweber iterative regularization method to solve the initial value problem of problem (1). We prove that the initial value problem is ill posed, which means that f ( x ) does not depend on the data g ( x ) continuously. Based on an a priori bound condition, the priori convergent error estimate under an a priori regularization parameter choice rule is given, and the posteriori convergent error estimate under an a posteriori regularization parameter choice rule is obtained. Some numerical examples show the effectiveness of this method.
At present, there are many effective regularization methods for the study of inverse problems, such as the truncation method [21,22], Tikhonov regularization method [23], quasi-boundary value method [24,25,26], quasi-reversibility regularization method [27,28], mollification regularization method [29], Fourier regularization method [30,31,32], and Landweber iterative regularization method [33,34,35], which never shows the saturation phenomenon. These regularization methods have been successfully applied to some inverse problems of mathematical and physical equations, and a lot of research results have been obtained. In this paper, we use the Landweber iterative regularization method, which does not have the saturation phenomenon; that is, it is order optimal for any p > 0 to solve this inverse problem.
The manuscript is organized as follows: The results of the ill-posed analysis and the conditional stability for identifying the initial value problem (1) are given in Section 2. In Section 3, the Landweber iterative regularization method is used to solve the inverse problem, and the priori and posteriori convergent error estimates are obtained. In Section 4, some numerical examples are given to convincingly demonstrate the effectiveness of the Landweber regularization method. A brief conclusion is presented in Section 5.

2. Ill-Posed Analysis and Conditional Stability Results for Problem (1)

In this section, we mainly give the results of the ill-posed analysis and conditional stability for the identification of the initial value problem of Problem (1). Let { λ n } n = 1 and { χ n } n = 1 be the Dirichlet eigenvalues and eigenfunctions of Δ on the domain Ω , respectively, and satisfy
Δ χ n ( x ) = λ n χ n ( x ) , x Ω , χ n ( x ) = 0 , x Ω ,
where 0 < λ 1 λ 2 λ n , lim n λ n = + , and χ n ( x ) H 2 ( Ω ) H 0 1 ( Ω ) is an orthonormal basis L 2 ( Ω ) .
Define
H p ( Ω ) = ϕ L 2 ( Ω ) | n = 1 ( λ n ) p | ( ϕ , χ n ) | 2 < , p > 0 ,
where ( · , · ) is the inner product in L 2 ( Ω ) , and H p ( Ω ) is a Hilbert space with the norm
ϕ H p ( Ω ) : = n = 1 ( λ n ) p | ( ϕ , χ n ) | 2 1 2 .
According to the result of theorem 2.1 in Bazhlekova, Jin, Lazarov, and Zhou’s paper [36], for any f ( x ) L 2 ( Ω ) , there exists a unique solution u ( x , t ) C ( [ 0 , T ] ; L 2 ( Ω ) ) C ( [ 0 , T ] ; H 2 ( Ω ) H 0 1 ( Ω ) ) and the solution for (1) is given by
u ( x , t ) = n = 1 f n u n ( t ) χ n ( x ) ,
where f n = ( f ( x ) , χ n ( x ) ) is the Fourier coefficients.
In Formula (7), the function u n ( t ) satisfies
u n ( t ) = 0 e s t B n ( s ) d s ,
where
B n ( s ) = γ π λ n s α sin α π ( s + λ n γ s α cos α π + λ n ) 2 + ( λ n γ s α sin α π ) 2 .
Letting t = T in (7), we have
g ( x ) = n = 1 f n u n ( T ) χ n ( x ) ,
and further, we have
g n = f n u n ( T ) ,
where g n = ( g ( x ) , χ n ( x ) ) is the Fourier coefficients.
From (9), we obtain the exact solution
f ( x ) = n = 1 g n u n ( T ) χ n ( x ) .
In order to analyze the ill posedness of the inverse problem and give the result of conditional stability, the following lemmas are useful for the whole paper.
Lemma 1
([36]). The functions u n ( t ) , n = 1 , 2 , , have the following properties:
(a)
u n ( 0 ) = 1 , 0 < u n ( t ) 1 , t 0 ;
(b)
u n ( t ) are completely monotone for t 0 ;
(c)
| λ n u n ( t ) | C min { t 1 , t α 1 } , t > 0 ;
where the constant C does not depend on n and t.
Lemma 2
([20]). Assuming that α ( 0 , 1 ) , the following estimate holds for all t [ 0 , T ]
u n ( t ) C ( γ , α , λ 1 ) λ n ,
where
C ( γ , α , λ 1 ) = γ sin α π 0 e s T s α d s γ 2 s 2 α + s 2 λ 1 2 + 1 .
According to Lemma 1, we know
1 u n ( T ) λ n C min { T 1 , T α 1 } .
Due to λ n ( n ) , it can be obtained 1 u n ( T ) . Accordingly, from Formula (10), it can be seen that the small perturbation of g δ ( x ) will cause a great change in f ( x ) ; i.e., it is an ill-posed problem.
Next, we will give the conditional stability result of this inverse problem.
Theorem 1.
When f ( x ) satisfies an a priori bound condition
f ( · ) H p ( Ω ) = n = 1 ( λ n ) p | ( f , χ n ) | 2 1 2 E , p > 0 ,
we have
f ( · ) C 1 E 2 p + 2 g ( · ) p p + 2 ,
where C 1 : = C p p + 2 ( γ , α , λ 1 ) ) is a positive constant.
Proof. 
Due to the Hölder inequality and (10), we have
f ( · ) 2 = n = 1 g n 2 u n 2 ( T ) = n = 1 g n 4 p + 2 u n 2 ( T ) g n 2 p p + 2 n = 1 f n 2 u n p ( T ) 2 p + 2 n = 1 g n 2 p p + 2 = n = 1 f n 2 u n p ( T ) 2 p + 2 g ( · ) 2 p p + 2 .
According to Lemma 2 and (11), we obtain
n = 1 f n 2 u n p ( T ) n = 1 f n 2 λ n p C p ( γ , α , λ 1 ) C p ( γ , α , λ 1 ) E 2 .
From (13) and (14), we have
f ( · ) C 1 E 2 p + 2 g ( x ) p p + 2 ,
where C 1 = C p p + 2 ( γ , α , λ 1 ) .
We completed the proof of Theorem 1. □
In the next section, we mainly introduce the Landweber iterative regularization method to solve this inverse problem. Furthermore, we obtain an a priori and an a posteriori convergence error estimate by using the proving techniques of trigonometric inequality and Hölder inequality [25,26].

3. Landweber Iterative Regularization Method and Convergence Analysis

In this section, we propose the Landweber iterative regularization method to solve the ill-posed problem (1). Under an a priori regularization parameter choice rule and an a posteriori regularization parameter choice rule, we obtain the convergent error estimates between the exact solution and the regularization solution. The inverse problem of identifying the initial value f ( x ) can be converted to solving the following integral equation:
( K f ) ( x ) : = Ω k ( x , ξ ) f ( ξ ) d ξ = g ( x ) ,
where the kernel function
k ( x , ξ ) = n = 1 u n ( T ) χ n ( x ) χ n ( ξ ) .
Due to the kernel function k ( x , ξ ) = k ( ξ , x ) , it can be seen that K is a self-adjoint operator. Further, we have the following theorem:
Theorem 2.
If f ( x ) L 2 ( Ω ) , the integral operator K in (15) is compact from L 2 ( Ω ) to L 2 ( Ω ) and its singular values are u n ( T ) .
Proof. 
From Lemma 1 and (8), we have
g ( · ) H 2 2 = n = 1 f n u n ( T ) χ n ( x ) H 2 2 = n = 1 ( λ n ) 2 | f n u n ( T ) | 2 n = 1 f n 2 C min { T 1 , T α 1 } λ n 2 ( λ n ) 2 = ( C min { T 1 , T α 1 } ) 2 f ( · ) 2 .
From (16), we know if f ( x ) L 2 ( Ω ) , g ( x ) H 2 ( Ω ) . Furthermore, due to H 2 ( Ω ) L 2 ( Ω ) , we can determine that the integral operator K is compact from L 2 ( Ω ) to L 2 ( Ω ) . Let K * be the adjoint of K . Since { χ n } n = 1 is the orthonormal basis in space L 2 ( Ω ) , it is easy to verify
K * K χ n ( ξ ) = ( u n ( T ) ) 2 χ n ( ξ ) .
Hence, the singular values of the compact operator K are u n ( T ) . The proof of Theorem 2 is complete. □
Next, we give the Landweber iterative regularization solution of f ( x ) . It must be borne in mind that f m , δ ( x ) is the Landweber iterative regularization solution of the initial value problem (1). We use the operator equation f = ( I a K * K ) f + a K * g to replace the equation K f = g and obtain the following iterative format:
f 0 , δ ( x ) = 0 , f m , δ ( x ) = ( I a K * K ) f m 1 , δ ( x ) + a K * g δ ( x ) , m = 1 , 2 , 3 , ,
where I is a unit operator; m is the iterative step number, also known as the regularization parameter; and a is called the relaxation factor and satisfies 0 < a < 1 K 2 . Note that the operator R m : L 2 ( Ω ) L 2 ( Ω ) is defined as
R m : = a k = 0 m 1 ( I a K * K ) k K * , m = 1 , 2 , 3 , .
By simple calculation, we have
f m , δ ( x ) = R m g δ ( x ) = a k = 0 m 1 ( I a K * K ) k K * g δ ( x ) .
Since K is a self-adjoint operator, applying the singular values of the operator K and Formula (17), we obtain the Landweber iterative regularization solution of the inverse problem (1) as follows:
f m , δ ( x ) = n = 1 1 ( 1 a u n 2 ( T ) ) m u n ( T ) g n δ χ n ( x ) ,
where g n δ = ( g δ ( x ) , χ n ( x ) ) .
In the following, we give two the convergent error estimates by using an a priori choice rule and an aposteriori choice rule for the regularization parameter.

3.1. The Convergent Error Estimate with an a Priori Parameter Choice Rule

Theorem 3.
Let f m , δ ( x ) given by (18) be the Landweber iterative regularization solution of the exact solution (10). Suppose the priori condition (11) and the noise assumption (3) hold. Choosing the regularization parameter m = b , where
b = E δ 4 p + 2 ,
we have the following convergence error estimate:
f m , δ ( · ) f ( · ) C 2 E 2 p + 2 δ p p + 2 ,
where b denotes the largest integer less than or equal to b and C 2 : = a + p 2 a C 2 ( γ , α , λ 1 ) p 4 is a positive constant.
Proof. 
According to the triangle inequality, we have
f m , δ ( · ) f ( · ) f m , δ ( · ) f m ( · ) + f m ( · ) f ( · ) .
From (3), we have
f m , δ ( · ) f m ( · ) = n = 1 1 ( 1 a u n 2 ( T ) ) m u n ( T ) g n δ χ n ( x ) n = 1 1 ( 1 a u n 2 ( T ) ) m u n ( T ) g n χ n ( x ) = n = 1 1 ( 1 a u n 2 ( T ) ) m u n ( T ) ( g n δ g n ) χ n ( x ) sup n 1 A ( n ) δ ,
where A ( n ) : = 1 ( 1 a u n 2 ( T ) ) m u n ( T ) .
Since u n ( T ) is a singular value of the operator K and 0 < a < 1 K 2 , we obtain 0 < a u n 2 ( T ) < 1 . Due to Bernoulli inequality, we have
1 ( 1 a u n 2 ( T ) ) m 1 ( 1 a u n 2 ( T ) ) m a m u n ( T ) ,
then we obtain
A ( n ) a m .
Further, we can obtain
f m , δ ( · ) f m ( · ) sup n 1 A ( n ) δ a m δ .
From (11), we have
f m ( · ) f ( · ) = n = 1 1 ( 1 a u n 2 ( T ) ) m u n ( T ) g n χ n ( x ) n = 1 1 u n ( T ) g n χ n ( x ) = n = 1 ( 1 a u n 2 ( T ) ) m u n ( T ) g n χ n ( x ) = n = 1 ( 1 a u n 2 ( T ) ) m ( λ n ) p 2 ( λ n ) p 2 f n χ n ( x ) sup n 1 B ( n ) E ,
where B ( n ) : = ( 1 a u n 2 ( T ) ) m ( λ n ) p 2 .
According to Lemma 2, we have
B ( n ) 1 a C 2 ( γ , α , λ 1 ) ( λ n ) 2 m ( λ n ) p 2 .
Let H ( s ) : = 1 a C 2 ( γ , α , λ 1 ) s 2 m s p 2 , s : = λ n . Supposing that s 0 satisfies H ( s 0 ) = 0 , we have
s 0 = a C 2 ( γ , α , λ 1 ) ( 4 m + p ) p 1 2 ,
then we have
H ( s ) H ( s 0 ) = 1 p 4 m + p m a C 2 ( γ , α , λ 1 ) ( 4 m + p ) p p 4 p 2 a C 2 ( γ , α , λ 1 ) p 4 ( m + 1 ) p 4 .
Hence, we have
f m ( · ) f ( · ) B ( n ) E H ( s ) E p 2 a C 2 ( γ , α , λ 1 ) p 4 ( m + 1 ) p 4 E .
Combining (19), (21)–(23), we obtain
f m , δ ( · ) f ( · ) C 2 E 2 p + 2 δ p p + 2 ,
where C 2 = a + p 2 a C 2 ( γ , α , λ 1 ) p 4 .
The proof of Theorem 3 is complete. □

3.2. The Convergent Error Estimate with an a Posteriori Parameter Choice Rule

In this section, we consider an a posteriori regularization parameter choice rule in the Morozov discrepancy principle [37] and obtain the convergent error estimate under an a posteriori regularization parameter choice rule. We assume that τ > 1 is given a fixed constant and stop the algorithm at the first occurrence of m = m ( δ ) N 0 with
K f m , δ ( · ) g δ ( · ) τ δ ,
where g δ τ δ is constant.
Lemma 3.
Let μ ( m ) = K f m , δ ( · ) g δ ( · ) ; then, we have the following conclusions:
(a)
μ ( m ) is a continuous function;
(b)
lim m 0 μ ( m ) = g δ ;
(c)
lim m + μ ( m ) = 0 ;
(d)
μ ( m ) is a strictly decreasing function for any m ( 0 , + ) .
Proof. 
The proof of Lemma 3 can be obtained by Formula (24), so it is omitted here. □
Remark 1.
According to the Lemma 3, we can find that the inequality (24) exists an unique solution.
Lemma 4.
Suppose the priori condition (11) and the noise assumption (3) hold. For fixed τ > 1 , if we choose the regularization parameter by using Morozov’s discrepancy principle (24), then the regularization parameter m = m ( δ ) satisfies
m p + 2 2 a C 2 ( γ , α , λ 1 ) C min { T 1 , T α 1 } τ 1 4 p + 2 E δ 4 p + 2 .
Proof. 
Firstly, we know that | 1 a u n 2 ( T ) | < 1 . From (3), we have
τ δ K f m 1 , δ g δ = n = 1 ( 1 a u n 2 ( T ) ) m 1 g n δ χ n ( x ) n = 1 ( 1 a u n 2 ( T ) ) m 1 ( g n δ g n ) χ n ( x ) + n = 1 ( 1 a u n 2 ( T ) ) m 1 g n χ n ( x ) δ + n = 1 ( 1 a u n 2 ( T ) ) m 1 g n χ n ( x ) .
Further, from (11) we can obtain
n = 1 ( 1 a u n 2 ( T ) ) m 1 g n χ n ( x ) = n = 1 ( 1 a u n 2 ( T ) ) m 1 u n ( T ) ( λ n ) p 2 ( λ n ) p 2 f n χ n ( x ) sup n 1 C ( n ) E ,
where C ( n ) : = ( 1 a u n 2 ( T ) ) m 1 u n ( T ) ( λ n ) p 2 .
According to Lemmas 1 and 2, we have
C ( n ) C min { T 1 , T α 1 } 1 a C 2 ( γ , α , λ 1 ) ( λ n ) 2 m ( λ n ) ( p 2 + 1 ) .
Let I ( s ) : = C min { T 1 , T α 1 } 1 a C 2 ( γ , α , λ 1 ) s 2 m s ( p 2 + 1 ) , s : = λ n . Supposing that s * satisfies I ( s * ) = 0 , we obtain
s * = a C 2 ( γ , α , λ 1 ) ( 4 m + p 2 ) p + 2 1 2 ,
then we have
I ( s * ) = C min { T 1 , T α 1 } 1 p 4 m + p 2 m 1 a C 2 ( γ , α , λ 1 ) ( 4 m + p 2 ) p p + 2 4 C min { T 1 , T α 1 } p 2 a C 2 ( γ , α , λ 1 ) p + 2 4 m p + 2 4 .
Hence, we obtain
n = 1 ( 1 a u n 2 ( T ) ) m 1 g n χ n ( x ) C ( n ) E C min { T 1 , T α 1 } p 2 a C 2 ( γ , α , λ 1 ) p + 2 4 m p + 2 4 E .
Thus, we obtain
m p + 2 2 a C 2 ( γ , α , λ 1 ) C min { T 1 , T α 1 } τ 1 4 p + 2 E δ 4 p + 2 .
The proof of Lemma 4 is complete. □
Theorem 4.
Let f m , δ ( x ) given by (18) be the Landweber iterative regularization solution of the exact solution (10). Suppose the priori condition (11) and the noise assumption (3) hold. If the regularization parameter m = m ( δ ) is chosen by Morozov’s discrepancy principle with stopping rule (24), then we have the following error estimate:
f m , δ ( · ) f ( · ) C 3 E 2 p + 2 δ p p + 2 ,
where C 3 : = p + 2 2 C 2 ( γ , α , λ 1 ) 1 2 C min { T 1 , T α 1 } τ 1 2 p + 2 + C 1 ( τ + 1 ) p p + 2 is a positive constant.
Proof. 
Using the triangle inequality, we have
f m , δ ( · ) f m ( · ) f m , δ ( · ) f m ( · ) + f m ( · ) f ( · ) .
Applying Lemma 4 and (3), we obtain
f m , δ ( · ) f m ( · ) a m δ p + 2 2 C 2 ( γ , α , λ 1 ) 1 2 C min { T 1 , T α 1 } τ 1 2 p + 2 E 2 p + 2 δ p p + 2 .
In addition, combining (3) and (24), we have
K ( f m ( · ) f ( · ) ) = n = 1 ( 1 a u n 2 ( T ) ) m g n χ n ( x ) = n = 1 ( 1 a u n 2 ( T ) ) m ( g n g n δ ) χ n ( x ) + n = 1 ( 1 a u n 2 ( T ) ) m g n δ χ n ( x ) n = 1 ( 1 a u n 2 ( T ) ) m ( g n g n δ ) χ n ( x ) + n = 1 ( 1 a u n 2 ( T ) ) m g n δ χ n ( x ) ( τ + 1 ) δ .
Hence, from (11) we can obtain
f m ( · ) f ( · ) H p ( Ω ) = n = 1 ( λ n ) p 2 ( 1 a u n 2 ( T ) ) m f n χ n ( x ) n = 1 ( λ n ) p 2 f n χ n ( x ) E .
Further, according to Theorem 1 and (30), we have
f m ( · ) f ( · ) C 1 ( τ + 1 ) p p + 2 E 2 p + 2 δ p p + 2 .
Combining (28), (29) and (31), we obtain
f m , δ ( · ) f ( · ) C 3 E 2 p + 2 δ p p + 2 ,
where C 3 = p + 2 2 C 2 ( γ , α , λ 1 ) 1 2 C min { T 1 , T α 1 } τ 1 2 p + 2 + C 1 ( τ + 1 ) p p + 2 .
The proof of Theorem 4 is complete. □

4. Numerical Implementation

In this section, we use Matlab software to give several numerical examples to verify the effectiveness of the Landweber iterative regularization method. The following numerical simulation consists of two parts. First, we use the known function f ( x ) to obtain the additional data value g ( x ) , which is a forward problem. Then, we use the additional data value g ( x ) to solve an inverse problem and obtain the regularization solution; i.e., we use the Landweber iterative algorithm to obtain the regularization solution f m , δ ( x ) . Letting d = 1 , Ω = ( 0 , a 0 ) , we consider a one-dimensional forward problem
t u ( x , t ) ( 1 + γ 0 R t α ) x x u ( x , t ) = 0 , x ( 0 , a 0 ) , t ( 0 , T ) , α ( 0 , 1 ) , u ( 0 , t ) = u ( a 0 , t ) = 0 , t ( 0 , T ) , u ( x , 0 ) = f ( x ) , x ( 0 , a 0 ) , u ( x , T ) = g ( x ) , x ( 0 , a 0 ) ,
where f ( x ) is known. Here, we use the finite difference method to discretize the above problem (32). In addition, two discrete schemes are introduced, i.e., the backward difference (BD) scheme and the implicit numerical optimization scheme (INAS) provided by reference [38]. Both schemes are unconditionally stable.
Define
t k = k Δ t ( k = 0 , 1 , , N ) , x i = i Δ x ( i = 0 , 1 , , M ) ,
where Δ t = 1 N is the step size of temporal direction and Δ x = 1 M is the step size of spatial direction. Then, the approximate value of u at each grid point is recorded as u i k u ( x i , t k ) . First, we give the BD iterative scheme. In the first step, the Riemann–Liouville operator (2) is discretized by the Grünwald–Letnikov formula [39]
0 R t α u ( x , t ) x i , t k = 1 h α j = 0 [ t k / h ] ϖ j ( α ) u ( x i , t k j h ) + O ( h q ) ,
where [ t k / h ] is the integer part of t k / h and ϖ j ( α ) are the coefficients of the generating function ϖ ( z , α ) = j = 0 ϖ j ( α ) z j . In formula (34), if q = 1 and ϖ ( z , α ) = ( 1 z ) α , it is simply called the Grünwald–Letnikov formula [39]. In this case, the coefficients ϖ j ( α ) = ( 1 ) j j α can be calculated by means of the recursive formulae
ϖ 0 ( α ) = 1 , ϖ j ( α ) = ( 1 α + 1 j ) ϖ j 1 ( α ) .
If h = Δ t , the expression (34) is indicated by u ( x , t ) evaluated at the grid points
0 R t α u ( x , t ) x i , t k = 1 ( Δ t ) α j = 0 k ϖ j ( α ) u i k j + O ( Δ t ) .
In the second step, we use backward difference formula to discretize the differential operators
t u ( x , t ) x i , t k = u i k u i k 1 Δ t + O ( Δ t )
and
x x u ( x , t ) x i , t k = u i + 1 k 2 u i k + u i 1 k ( Δ x ) 2 + O ( Δ x ) 2 .
By using Formulas (36)–(38), we obtain the BD iterative form of Problem (32) expressed as
u i k = u i k 1 + μ 1 δ x 2 u i k + μ 2 j = 0 k ϖ j ( α ) δ x 2 u i k j , k = 1 , 2 , , N , i = 1 , 2 , , M ,
where μ 1 : = Δ t / ( Δ x ) 2 , μ 2 : = γ ( Δ t ) 1 α / ( Δ x ) 2 and δ x 2 u i k : = u i + 1 k 2 u i k + u i 1 k .
Secondly, according to the literature [38], we obtain the INAS iteration form of Problem (32) as
u i 1 = u i 0 + ( ν 1 + ν 2 ) δ x 2 u i 1 , u i k + 1 = u i k + ( ν 1 + ν 2 ) δ x 2 u i k + 1 + ν 1 j = 0 k 1 ( ρ j + 1 ρ j ) δ x 2 u i k j , k = 1 , 2 , , N 1 , i = 1 , 2 , , M 1 ,
where ν 1 : = γ ( Δ t ) 1 α / Γ ( 2 α ) ( Δ x ) 2 , ν 2 : = Δ t / ( Δ x ) 2 and b j : = ( j + 1 ) 1 α j 1 α ( j = 0 , 1 , , N 1 ) .
Through the above two discrete formats, we use Matlab software to program and run the program to obtain the data function g. Next, we solve an inverse problem.
In practical application, the data g are obtained by measurements and have certain error. Therefore, in the numerical simulation, we add a random disturbance to the data g. The noisy data g δ are generated by adding random disturbances, i.e.,
g δ = g + ε · r a n d n ( s i z e ( g ) ) ,
where the function r a n d n ( · ) represents the generation of a column of random numbers with a mean value of 0 and a variance of 1, and ε represents the relative error level. The absolute error level δ is expressed as
δ = 1 M + 1 i = 1 M + 1 ( g i g i δ ) 2 .
Finally, the regularization solution f m , δ is obtained by the following formula:
f m , δ ( x ) = a k = 0 m 1 ( I a K * K ) k K * g δ ( x ) ,
where K satisfies K f = g δ .
To see the accuracy of numerical solutions, we compute the relative root mean square errors by
ζ ( f ) = i = 1 n ( f m , δ ( x i ) f ( x i ) ) 2 / i = 1 n f ( x i ) 2 1 / 2 ,
where n is the total number of test points.
The priori regularization parameter is based on the smooth conditions of the exact solution, which is difficult to give in practical problems. The following examples are based on an a posterior regularization parameter choice rule (24) to verify the effectiveness of the Landweber iterative regularization method. By simple calculation, we obtain λ n = ( n π ) 2 and χ n = 2 sin n π x for n = 1 , 2 , in Formula (4). In the numerical calculation of Problem (32), we choose γ = 1 , a 0 = 2 , T = 1 .
Let us take p = 2 , τ = 1.01 in Formulas (11) and (24). Choosing M = 80 , N = 50 , we give the following three examples:
Example 1.
Consider a smooth function f ( x ) = sin ( 2 π x ) , x [ 0 , 2 ] .
Example 2.
Consider a piecewise smooth function
f ( x ) = 0 , 0 x 1 2 , 4 ( x 1 2 ) , 1 2 < x 1 , 4 ( x 3 2 ) , 1 < x 3 2 , 0 , 3 2 < x 2 .
Example 3.
Consider a non-smooth function
f ( x ) = 0 , x [ 0 , 0.4 ] , 1 , x ( 0.4 , 0.8 ] , 0 , x ( 0.8 , 1.2 ] , 1 , x ( 1.2 , 1.6 ] , 0 , x ( 1.6 , 2 ] .
Figure 1, Figure 2 and Figure 3 show a comparison of the exact solution f ( x ) and its approximate solution f m , δ ( x ) between BD and INAS in the iterative form of Example 1 for the relative error levels ε = 0.05 , 0.01 , 0.005 with various values of α = 0.2 , 0.5 , 0.7 , 0.9 . Table 1 shows the relative root mean square error ζ ( f ) differences between the exact solution and the regularization solution of Example 1 for various values of α = 0.2 , 0.5 , 0.7 , 0.9 and ε = 0.05 , 0.01 , 0.005 . Table 2 shows a comparison between the number of iterations ( m ) for the exact solution and the regularization solution of Example 1 for various values of α = 0.2 , 0.5 , 0.7 , 0.9 and ε = 0.05 , 0.01 , 0.005 . Figure 4 shows a comparison between the exact solution and the regularization solution of the two iterative methods for ε = 0.05 , 0.01 , 0.005 under α = 0.6 .
Figure 5, Figure 6 and Figure 7 show a comparison of the exact solution f ( x ) and its approximate solution f m , δ ( x ) between BD and INAS in the iterative form of Example 2 for the relative error levels ε = 0.05 , 0.01 , 0.005 with various values of α = 0.2 , 0.5 , 0.7 , 0.9 . Table 3 shows a comparison of the relative root mean square errors ζ ( f ) for the exact solution and the regularization solution of Example 2 for various values of α = 0.2 , 0.5 , 0.7 , 0.9 and ε = 0.05 , 0.01 , 0.005 . Table 4 shows a comparison of the number of iterations ( m ) for the exact solution and the regularization solution of Example 2 for various values of α = 0.2 , 0.5 , 0.7 , 0.9 and ε = 0.05 , 0.01 , 0.005 . Figure 8 shows a comparison between the exact solution and the regularization solution of the two iterative methods for ε = 0.05 , 0.01 , 0.005 under α = 0.8 .
Figure 9 and Figure 10 show a comparison of the exact solution f ( x ) and its approximate solution f m , δ ( x ) between BD and INAS in the iterative form of Example 3 for the relative error levels ε = 0.05 , 0.01 with various values of α = 0.2 , 0.5 , 0.7 , 0.9 . Table 5 shows a comparison of the relative root mean square errors ζ ( f ) of the exact solution and the regularization solution of Example 3 for various values of α = 0.2 , 0.5 , 0.7 , 0.9 and ε = 0.05 , 0.01 . Table 6 shows a comparison of the number of iterations ( m ) for the exact solution and the regularization solution of Example 3 for various values of α = 0.2 , 0.5 , 0.7 , 0.9 and ε = 0.05 , 0.01 . Figure 11 shows a comparison between the exact solution and the regularization solution of the two iterative methods for ε = 0.05 , 0.01 , 0.005 under α = 0.95 .

5. Conclusions

In this paper, an inverse problem to determine the initial value problem of the homogeneous Rayleigh–Stokes equation for a generalized second-grade fluid with the Riemann–Liouville fractional derivative model is studied. The Landweber iterative regularization method is used to solve the ill-posed problem. Under a conditional stability result, we obtain the convergent error estimates between the exact solution and the regularization solution by using an a priori regularization parameter choice rule and an a posteriori regularization parameter choice rule. Finally, three numerical examples are given to illustrate the effectiveness of this method. In the future, we will consider the other fractional derivative model of the homogeneous Rayleigh–Stokes equation for a generalized second-grade fluid and consider the inverse problem of identifying the unknown source and initial value for 1 < α < 2 .

Author Contributions

The main idea of this paper was proposed by D.-G.L.; J.-L.F. and F.Y. prepared the manuscript initially and performed all the steps of the proofs in this research. D.-G.L.; J.-L.F.; F.Y. and X.-X.L. All authors have read and agreed to the published version of the manuscript.

Funding

The project is supported by the National Natural Science Foundation of China (No.11961044) and the Doctor Fund of Lan Zhou University of Technology, the Natural Science Foundation of Gansu Province (No. 21JR7RA214).

Institutional Review Board Statement

Institutional review board approves this study.

Informed Consent Statement

All study participants all consent prior to study enrollment.

Data Availability Statement

No data is generated during the study.

Acknowledgments

The authors would like to thanks the editor and the referees for their valuable comments and suggestions that improve the quality of our paper.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Baleanu, D.; Jajarmi, A.; Mohammadi, H.; Rezapour, S. A new study on the mathematical modelling of human liver with Caputo-Fabrizio fractional derivative. Chaos Solitons Fract. 2020, 134, 109705. [Google Scholar] [CrossRef]
  2. Rezapour, S.; Mohammadi, H.; Jajarmi, A. A new mathematical model for Zika virus transmission. Adv. Differ. Eq. 2020, 2020, 589. [Google Scholar] [CrossRef]
  3. Baleanu, D.; Mohammadi, H.; Rezapour, S. Analysis of the model of HIV-1 infection of CD4+ T-cell with a new approach of fractional derivative. Adv. Differ. Eq. 2020, 2020, 71. [Google Scholar] [CrossRef] [Green Version]
  4. Tuan, N.H.; Ganji, R.M.; Jafari, H. A numerical study of fractional rheological models and fractional Newell-Whitehead-Segel equation with non-local and non-singular kernel. Chin. J. Phys. 2020, 68, 308–320. [Google Scholar] [CrossRef]
  5. Nikan, O.; Jafari, H.; Golbabai, A. Numerical analysis of the fractional evolution model for heat flow in materials with memory. Alex. Eng. J. 2020, 59, 2627–2637. [Google Scholar] [CrossRef]
  6. Ganji, R.M.; Jafari, H.; Baleanu, D. A new approach for solving multi variable orders differential equations with Mittag-Leffler kernel. Chaos Solitons Fract. 2020, 130, 109405. [Google Scholar] [CrossRef]
  7. Baleanu, D.; Etemad, S.; Rezapour, S. A hybrid Caputo fractional modeling for thermostat with hybrid boundary value conditions. Bound Value Probl. 2020, 2020, 64. [Google Scholar] [CrossRef] [Green Version]
  8. Shen, F.; Tan, W.; Zhao, Y.; Masuoka, T. The Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivative model. Nonlinear Anal. Real World Appl. 2006, 7, 1072–1080. [Google Scholar] [CrossRef]
  9. Fetecau, C.; Jamil, M.; Fetecau, C.; Vieru, D. The Rayleigh–Stokes problem for an edge in a generalized Oldroyd-B fluid. Z. Angew. Math. Phys. 2009, 60, 921–933. [Google Scholar] [CrossRef]
  10. Song, D.Y.; Jiang, T.Q. Study on the constitutive equation with fractional derivative for the viscoelastic fluids-Modified Jeffreys model and its application. Rheol. Acta 1998, 37, 512–517. [Google Scholar] [CrossRef]
  11. Zhao, C.L.; Yang, C. Exact solutions for electro-osmotic flow of viscoelastic fluids in rectangular micro-channels. Appl. Math. Comput. 2009, 211, 502–509. [Google Scholar] [CrossRef]
  12. Khan, M. The Rayleigh–Stokes problem for an edge in a viscoelastic fluid with a fractional derivative model. Nonlinear Anal. Real World Appl. 2009, 10, 3190–3195. [Google Scholar] [CrossRef]
  13. Khan, M.; Anjum, A.; Qi, H.; Fetecau, C. On exact solutions for some oscillating motions of a generalized Oldroyd-B fluid. Z. Angew. Math. Phys. 2010, 61, 133–145. [Google Scholar] [CrossRef]
  14. Chen, C.M.; Liu, F.; Anh, V. Numerical analysis of the Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivatives. Appl. Math. Comput. 2008, 204, 340–351. [Google Scholar] [CrossRef] [Green Version]
  15. Chen, C.M.; Liu, F.; Burrage, K.; Chen, Y. Numerical methods of the variable-order Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivative. IMA J. Appl. Math. 2013, 78, 924–944. [Google Scholar] [CrossRef]
  16. Elyas, S.; Ahmad, J. Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivatives: A stable scheme based on spectral meshless radial point interpolation. Eng. Comput. 2017, 34, 77–90. [Google Scholar]
  17. Dehghan, M.; Abbaszadeh, M. A finite element method for the numerical solution of Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivatives. Eng. Comput. 2017, 33, 1–19. [Google Scholar] [CrossRef]
  18. Nikan, O.; Golbabai, A.; Tenreiro Machado, J.A.; Nikazad, T. Numerical solution of the fractional Rayleigh–Stokes model arising in a heated generalized second-grade fluid. Eng. Comput. 2021. [Google Scholar] [CrossRef]
  19. Chen, C.M.; Liu, F.; Turner, I.; Anh, V. Numerical methods with fourth-order spatial accuracy for variable-order nonlinear Stokes’ first problem for a heated generalized second grade fluid. Comput. Math. Appl. 2011, 62, 971–986. [Google Scholar] [CrossRef] [Green Version]
  20. Hoang, L.N.; Huy, T.N.; Kirane, M.; Xuan, T.D.D. Identifying initial condition of the Rayleigh–Stokes problem with random noise. Math. Methods Appl. Sci. 2019, 42, 1561–1571. [Google Scholar]
  21. Yang, F.; Zhang, P.; Li, X.X. The truncation method for the Cauchy problem of the inhomogeneous Helmholtz equation. Appl. Anal. 2019, 98, 991–1004. [Google Scholar] [CrossRef]
  22. Yang, F.; Fan, P.; Li, X.X. Fourier truncation regularization method for a three-dimensional Cauchy problem of the modified Helmholtz equation with perturbed wave number. Mathematics 2019, 7, 705. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, J.G.; Wei, T.; Zhou, Y.B. Tikhonov regularization method for a backward problem for the time-fractional diffusion equation. Appl. Math. Model. 2013, 37, 8518–8532. [Google Scholar] [CrossRef]
  24. Feng, X.L.; Eldén, L. Solving a Cauchy problem for a 3D elliptic PDE with variable coefficients by a quasi-boundary-value method. Inverse Probl. 2014, 30, 015005. [Google Scholar] [CrossRef]
  25. Yang, F.; Wang, N.; Li, X.X. A quasi-boundary regularization method for identifying the initial value of time-fractional diffusion equation on spherically symmetric domain. J. Inverse Ill-Posed Probl. 2019, 27, 609–621. [Google Scholar] [CrossRef]
  26. Yang, F.; Sun, Y.R.; Li, X.X. The quasi-boundary value method for identifying the initial value of heat equation on a columnar symmetric domain. Numer. Algorithms 2019, 82, 623–639. [Google Scholar] [CrossRef]
  27. Zhang, H.W.; Qin, H.H.; Wei, T. A quasi-reversibility regularization method for the Cauchy problem of the Helmholtz equation. Int. J. Comput. Math. 2011, 88, 839–850. [Google Scholar] [CrossRef]
  28. Yang, F.; Fu, C.L. The quasi-reversibility regularization method for identifying the unknown source for time fractional diffusion equation. Appl. Math. Model. 2015, 39, 1500–1512. [Google Scholar] [CrossRef]
  29. Yang, F.; Fu, C.L.; Li, X.X. A mollification regularization method for unknown source in time-fractional diffusion equation. Int. J. Comput. Math. 2014, 91, 1516–1534. [Google Scholar] [CrossRef]
  30. Xiong, X.T.; Fu, C.L.; Li, H.F. Fourier regularization method of a sideways heat equation for determining surface heat flux. J. Math. Anal. Appl. 2006, 317, 331–348. [Google Scholar] [CrossRef] [Green Version]
  31. Li, X.X.; Lei, J.L.; Yang, F. An a posteriori Fourier regularization method for identifying the unknown source of the space-fractional diffusion equation. J. Inequal. Appl. 2014, 2014, 1–13. [Google Scholar] [CrossRef] [Green Version]
  32. Yang, F.; Fu, C.L.; Li, X.X. The Fourier regularization method for identifying the unknown source for the modified Helmholtz equation. Acta Math. Sin. 2014, 34, 1040–1047. [Google Scholar]
  33. Yang, F.; Zhang, Y.; Li, X.X. Landweber iterative method for identifying the initial value problem of the time-space fractional diffusion-wave equation. Numer. Algorithms 2019, 83, 1509–1530. [Google Scholar] [CrossRef]
  34. Yang, F.; Ren, Y.P.; Li, X.X. Landweber iteration regularization method for identifying unknown source on a columnar symmetric domain. Inverse Probl. Sci. Eng. 2018, 26, 1109–1129. [Google Scholar] [CrossRef]
  35. Yang, F.; Liu, X.; Li, X.X. Landweber iterative regularization method for identifying the unknown source of the modified Helmholtz equation. Bound. Value Probl. 2017, 2017, 1–16. [Google Scholar] [CrossRef]
  36. Bazhlekova, E.; Jin, B.T.; Lazarov, R.; Zhou, Z. An analysis of the Rayleigh–Stokes problem for a generalized second-grade fluid. Numer. Math. 2015, 131, 1–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Engl, H.W.; Hanke, M.; Neubauer, A. Regularization of Inverse Problems; Kluwer Academic Publishers: Boston, MA, USA, 2019. [Google Scholar]
  38. Wu, C. Numerical solution for Stokes’ first problem for a heated generalized second grade fluid with fractional derivative. Appl. Numer. Math. 2009, 59, 2571–2583. [Google Scholar] [CrossRef]
  39. Yuste, S.B. Weighted average finite difference methods for fractional diffusion equations. J. Comput. Phys. 2006, 216, 264–274. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The exact solution and its approximate solution of Example 1 for ε = 0.05 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Figure 1. The exact solution and its approximate solution of Example 1 for ε = 0.05 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Fractalfract 05 00193 g001
Figure 2. The exact solution and its approximate solution of Example 1 for ε = 0.01 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Figure 2. The exact solution and its approximate solution of Example 1 for ε = 0.01 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Fractalfract 05 00193 g002
Figure 3. The exact solution and its approximate solution of Example 1 for ε = 0.005 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Figure 3. The exact solution and its approximate solution of Example 1 for ε = 0.005 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Fractalfract 05 00193 g003
Figure 4. The exact solution and its approximate solution of Example 1 with α = 0.6 for ε = 0.05 , 0.01 , 0.005 . (a) BD: m = 92, 135, 154, ζ ( f ) = 0.0714, 0.0020, 0.0007; (b) INAS: m = 2369, 3181, 3543, ζ ( f ) = 0.0029, 0.0004, 8.9000 × 10 5 .
Figure 4. The exact solution and its approximate solution of Example 1 with α = 0.6 for ε = 0.05 , 0.01 , 0.005 . (a) BD: m = 92, 135, 154, ζ ( f ) = 0.0714, 0.0020, 0.0007; (b) INAS: m = 2369, 3181, 3543, ζ ( f ) = 0.0029, 0.0004, 8.9000 × 10 5 .
Fractalfract 05 00193 g004
Figure 5. The exact solution and its approximate solution of Example 2 for ε = 0.05 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Figure 5. The exact solution and its approximate solution of Example 2 for ε = 0.05 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Fractalfract 05 00193 g005
Figure 6. The exact solution and its approximate solution of Example 2 for ε = 0.01 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Figure 6. The exact solution and its approximate solution of Example 2 for ε = 0.01 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Fractalfract 05 00193 g006
Figure 7. The exact solution and its approximate solution of Example 2 for ε = 0.005 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Figure 7. The exact solution and its approximate solution of Example 2 for ε = 0.005 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Fractalfract 05 00193 g007
Figure 8. The exact solution and its approximate solution of Example 2 with α = 0.8 for ε = 0.05 , 0.01 , 0.005 . (a) BD: m = 16, 26, 30, ζ ( f ) = 0.0694, 0.0027, 7.3236 × 10 4 ; (b) INAS: m = 337, 1043, 1565, ζ ( f ) = 0.5461, 0.2424, 1403.
Figure 8. The exact solution and its approximate solution of Example 2 with α = 0.8 for ε = 0.05 , 0.01 , 0.005 . (a) BD: m = 16, 26, 30, ζ ( f ) = 0.0694, 0.0027, 7.3236 × 10 4 ; (b) INAS: m = 337, 1043, 1565, ζ ( f ) = 0.5461, 0.2424, 1403.
Fractalfract 05 00193 g008
Figure 9. The exact solution and its approximate solution of Example 3 for ε = 0.05 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Figure 9. The exact solution and its approximate solution of Example 3 for ε = 0.05 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Fractalfract 05 00193 g009
Figure 10. The exact solution and its approximate solution of Example 3 for ε = 0.01 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Figure 10. The exact solution and its approximate solution of Example 3 for ε = 0.01 . (a) α = 0.2 , (b) α = 0.5 , (c) α = 0.7 , (d) α = 0.9 .
Fractalfract 05 00193 g010
Figure 11. The exact solution and its approximate solution of Example 3 with α = 0.95 for ε = 0.05 , 0.01 , 0.005 . (a) BD: m = 43,62,71, ζ ( f ) = 0.0443, 0.0031, 5.9184 × 10 4 ; (b) INAS: m = 1841, 20,875, 28,442, ζ ( f ) = 0.3082, 0.2395, 0.2069.
Figure 11. The exact solution and its approximate solution of Example 3 with α = 0.95 for ε = 0.05 , 0.01 , 0.005 . (a) BD: m = 43,62,71, ζ ( f ) = 0.0443, 0.0031, 5.9184 × 10 4 ; (b) INAS: m = 1841, 20,875, 28,442, ζ ( f ) = 0.3082, 0.2395, 0.2069.
Fractalfract 05 00193 g011
Table 1. Comparison between the relative root mean square errors of the exact solution and the regularization solution of Example 1 for various values of α and ε .
Table 1. Comparison between the relative root mean square errors of the exact solution and the regularization solution of Example 1 for various values of α and ε .
α = 0.2 α = 0.5 α = 0.7 α = 0.9
ζ ( f ) ε = 0.05 BD0.04100.03120.02720.0288
INAS0.01410.01900.01930.0151
ε = 0.01 BD0.00990.00630.00650.0059
INAS0.00390.00590.00300.0019
ε = 0.005 BD0.00460.00320.00270.0027
INAS0.00230.00210.00160.0017
Table 2. Comparison of the number of iterations for the exact solution and the regularization solution of Example 1 for various values of α and ε .
Table 2. Comparison of the number of iterations for the exact solution and the regularization solution of Example 1 for various values of α and ε .
α = 0.2 α = 0.5 α = 0.7 α = 0.9
Number of iterations (m) ε = 0.05 BD143623625631
INAS26,36636231529763
ε = 0.01 BD19,9884807339
INAS35,055459221591026
ε = 0.005 BD22,1175348043
INAS37,551536223741070
Table 3. Comparison of the relative root mean square errors for the exact solution and the regularization solution of Example 2 for various values of α and ε .
Table 3. Comparison of the relative root mean square errors for the exact solution and the regularization solution of Example 2 for various values of α and ε .
α = 0.2 α = 0.5 α = 0.7 α = 0.9
ζ ( f ) ε = 0.05 BD0.36720.14410.08450.0734
INAS0.19100.13600.15790.1500
ε = 0.01 BD0.11060.02880.01590.0169
INAS0.09630.08490.08260.0905
ε = 0.005 BD0.07480.01700.00760.0077
INAS0.06410.05760.06380.0690
Table 4. Comparison of the number of iterations for the exact solution and the regularization solution of Example 2 for various values of α and ε .
Table 4. Comparison of the number of iterations for the exact solution and the regularization solution of Example 2 for various values of α and ε .
α = 0.2 α = 0.5 α = 0.7 α = 0.9
Number of iterations (m) ε = 0.05 BD926882114
INAS3822638265172
ε = 0.01 BD44701983623
INAS12,48627381127550
ε = 0.005 BD65022414327
INAS25,34343681538765
Table 5. Comparison of the relative root mean square errors for the exact solution and the regularization solution of Example 3 for various values of α and ε .
Table 5. Comparison of the relative root mean square errors for the exact solution and the regularization solution of Example 3 for various values of α and ε .
α = 0.2 α = 0.5 α = 0.7 α = 0.9
ζ ( f ) ε = 0.05 BD0.18500.05350.03540.0388
INAS0.36930.33160.28840.3108
ε = 0.01 BD0.04260.01000.00710.0066
INAS0.26410.28380.27290.2441
Table 6. Comparison of the number of iterations for the exact solution and the regularization solution of Example 3 for various values of α and ε .
Table 6. Comparison of the number of iterations for the exact solution and the regularization solution of Example 3 for various values of α and ε .
α = 0.2 α = 0.5 α = 0.7 α = 0.9
Number of iterations (m) ε = 0.05 BD14,0902564123
INAS24,384649167142243
ε = 0.01 BD47,5554205632
INAS275,52730,49615,46315,362
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, D.-G.; Fu, J.-L.; Yang, F.; Li, X.-X. Landweber Iterative Regularization Method for Identifying the Initial Value Problem of the Rayleigh–Stokes Equation. Fractal Fract. 2021, 5, 193. https://doi.org/10.3390/fractalfract5040193

AMA Style

Li D-G, Fu J-L, Yang F, Li X-X. Landweber Iterative Regularization Method for Identifying the Initial Value Problem of the Rayleigh–Stokes Equation. Fractal and Fractional. 2021; 5(4):193. https://doi.org/10.3390/fractalfract5040193

Chicago/Turabian Style

Li, Dun-Gang, Jun-Liang Fu, Fan Yang, and Xiao-Xiao Li. 2021. "Landweber Iterative Regularization Method for Identifying the Initial Value Problem of the Rayleigh–Stokes Equation" Fractal and Fractional 5, no. 4: 193. https://doi.org/10.3390/fractalfract5040193

APA Style

Li, D. -G., Fu, J. -L., Yang, F., & Li, X. -X. (2021). Landweber Iterative Regularization Method for Identifying the Initial Value Problem of the Rayleigh–Stokes Equation. Fractal and Fractional, 5(4), 193. https://doi.org/10.3390/fractalfract5040193

Article Metrics

Back to TopTop