Next Article in Journal
Characterization Results of Solution Sets Associated with Multiple-Objective Fractional Optimal Control Problems
Next Article in Special Issue
Tomographic Reconstruction: General Approach to Fast Back-Projection Algorithms
Previous Article in Journal
Fixed Point Results in Soft Fuzzy Metric Spaces
Previous Article in Special Issue
Combination of Multigrid with Constraint Data for Inverse Problem of Nonlinear Diffusion Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Reconstruction of a Space-Dependent Reaction Coefficient and Initial Condition for a Multidimensional Wave Equation with Interior Degeneracy

1
Department of Mathematics, Faculty of Sciences, University of Nouakchott Al Aasriya, Nouakchott Bp 6093, Mauritania
2
Department of Mathematics and Statistics, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
3
Department of Applied Mathematics, National Research Centre, Cairo 12622, Egypt
4
Beheer en Algemene Directie, Ghent University Hospital, Corneel Heymanslaan 10, B-9000 Ghent, Belgium
5
Research Department, Ghent University, Sint-Pietersnieuwstraat 25, B-9000 Ghent, Belgium
6
Department of Computational Mathematics and Computer Science, Institute of Natural Sciences and Mathematics, Ural Federal University, 19 Mira St., Yekaterinburg 620002, Russia
7
Department of Mathematics, Faculty of Science, Benha University, Benha 13511, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(14), 3186; https://doi.org/10.3390/math11143186
Submission received: 22 April 2023 / Revised: 8 July 2023 / Accepted: 10 July 2023 / Published: 20 July 2023
(This article belongs to the Special Issue Inverse Problems and Imaging: Theory and Applications)

Abstract

:
A simultaneous reconstruction of the initial condition and the space-dependent reaction coefficient in a multidimensional hyperbolic partial differential equation with interior degeneracy is of concern. A temporal integral observation is utilized to achieve that purpose. The well-posedness, existence, and uniqueness of the inverse problem under consideration are discussed. The inverse problem can be reformulated as a least squares minimization and the Fréchet gradients are determined, using the adjoint and sensitivity problems. Finally, an iterative construction procedure is developed by employing the conjugate gradient algorithm while invoking the discrepancy principle as a stopping criterion. Some numerical experiments are given to ensure the performance of the reconstruction scheme in one and two dimensions.

1. Introduction

Numerous significant applications in cosmology [1], data science [2], remote sensing [3], medicine [4], and geophysics [5] are modeled in inverse problems involving the determination of unknown coefficients of partial differential equations based on limited information about the system over a finite period. Degenerate wave models are gaining more focus these days in many physical applications [6,7,8]. A survey of the numerical techniques that have been applied to direct and inverse problems with integer or fractional order derivatives indicates a lot of focus in recent days [9,10,11]. A reconstruction of missing sole terms in different styles of time-dependent fractional diffusion problems has been seen in [12,13].
The problem under consideration invokes the one proposed in [14] in which the reconstruction is only targeted to the potential term. We address here the reconstruction of two factors (initial condition and potential) in a multidimensional wave problem with interior degeneracy.
t t φ · ( a φ ) + b φ = Q in I , φ = 0 in ( R d χ ) × ( 0 , S ] , φ ( · , 0 ) = u 0 ( x ) in χ , φ t ( · , 0 ) = v 0 ( x ) in χ ,
where I : = χ × ( 0 , S ] . Endowed with the final observation data φ ( · , S ) in the domain χ , where S > 0 and χ R d ( d 1 ) , Q L 2 ( χ × ( 0 , S ) ) , b L 2 ( χ ) i s a C 1 ( χ ) a function that degenerates into a point x 0 inside the spatial domain χ with { 1 a ( x ) 0 , x χ } . The Hilbert space is defined as
H a 1 ( χ ) : = φ W 0 1 , 1 ( χ ) : a φ L 2 ( χ ) ,
with the inner product
φ , v H a 1 : = χ a φ · v d x + χ φ v d x .
An important application of the inverse problem (1) is to distinguish between various types of seismic events, such as implosion, explosion, or earthquake, which generate waves that propagate through the Earth and can be recorded using seismometers. In [15], a seismic source modeled as a point moment tensor forcing in the elastic wave equation for the displacement was estimated by minimizing the gap between the time-dependent measured/recorded and computed waveforms (see [16]). The weak formulation of (1) is:
χ t t φ v d x + χ a ( x ) φ · v d x + χ b φ v d x = χ Q v d x , v H 0 1 ( χ ) .
In the literature, the inverse problem of determining the coefficient b ( x ) over a large scale from time t = S and its time-integratation for temperature, among other data, has been studied and used to prove its existence and uniqueness when the source term Q, the Dirichlet boundary conditions, and the initial condition u 0 are known [17,18,19,20].
In addition, several numerical methods have been proposed. Examples include the standard regularization method of the Tikhonov type [21], the method of Armijo combined with the finite element method [22], the (NAG E04FCF) combined with the method of finite differences (FDM) [23], and the conjugate gradient method (CGM) [24] to reconstruct the coefficient b ( x ) numerically from the additional measurements.
Extensive research has been conducted on the inverse problem of determining the initial condition from time-integral temperature measurements and a final instant when the reaction potential and the source term are known (see, for example, [25,26]).
In this paper, we study the generalization of the conditions used in [27], which amounts to knowing the direct solution at certain instants of the time domain and at its end. This is difficult in practice to implement with precision at recordings with an average time, which reduces the possibility of significant measurement errors in the direct solution φ . More specifically, using the average weighted integral observations given in (3) and (4), we study the inverse problem to determine the pair ( b ( x ) and u 0 ( x ) ) in (1). Let ρ 1 ( t ) , ρ 2 ( t ) be two approximations to the delta function at t = S , such that ρ 1 ( t ) , ρ 2 ( t ) C 1 ( 0 , T ) , which are given, and ψ 1 ( x ) and ψ 2 ( x ) are the average weighted integral observations that are also given. Let
0 S ρ 1 ( t ) φ ( x , t ) d t = ψ 1 ( x ) , x χ
0 S ρ 2 ( t ) φ ( x , t ) d t = ψ 2 ( x ) , x χ .
Taking into account the generalizations of the final observations φ ( t = S ) , the integral observations (3) and (4) can be thought of in this way. However, note that the selection of the weighted functions (3) and (4) has a crucial role in obtaining the relevant data to recover the two unknown values ( b ( x ) , u 0 ) ; for more information, refer to [28]. On the other hand, in the literature, the determination of the reaction potential and the source term have been studied from the final observation of time in [27] and the measurement of the integral observation in time in [29] for non-degenerate parabolic problems. In the inverse problems (1), (3), and (4), these approaches can also be used to simultaneously calculate the response coefficient b ( x ) and the initial condition u 0 ( x ) .
The manuscript is arranged to have the uniqueness of the inverse problem in Section 2. The stability and the regularity results, the proof of the Fréchet differentiability of the objective functional, and the conjugate gradient and sensitivity problems are presented in Section 3. Using the conjugate gradient method (CGM) regularized by the discordance principle [30], the inverse problem is solved numerically in a stable way. Section 4 is devoted to the numerical simulations and their good agreement with the theoretical analysis. The well-posedness of the direct problem of (1) is discussed in the following theorem.
Theorem 1. 
Assume that v 0 L 2 ( χ ) , 0 b C 1 and u 0 H a 1 ( χ ) , Q L 2 ( 0 , S , χ ) . The problem (1) has the following unique weak solution
φ X 0 = L 2 0 , S ; H a 1 ( χ ) L 0 , S ; L 2 ( χ ) , t φ L 2 ( 0 , S ; L 2 ( χ ) ) ,
with
sup t 0 , S φ ( t ) L 2 ( χ ) 2 + t φ L 2 ( 0 , S ; L 2 ( χ ) ) 2 + a ( x ) φ L 2 ( 0 , S ; L 2 ( χ ) ) 2 C u 0 H a 1 ( χ ) 2 + v 0 L 2 ( χ ) 2 + Q L 2 ( 0 , S ; L 2 ( χ ) ) 2 .
The constant C depends on χ and S.
Proof. 
The proof of the following theorem relies on the outlines used to prove the existence and uniqueness theorem for the degenerate linear viscose-elastic problem presented in [31].    □

2. Well-Posedness of the Inverse Problem

Firstly, we introduce the following admissible set:
E = { b L ( χ ) : 0 b 1 b ( x ) b 2 , a . e . x χ }
A reformulation of the inverse problems (1), (3), and (4) as a nonlinear non-classical parabolic problem is targeted to simplify the proof of existence and uniqueness. By multiplying the first equation in (1) by ρ 1 ( t ) and ρ 2 ( t ) , respectively, integrating the resulting relations with respect to t from 0 to S, and using (3) and (4), we have
ρ i φ ( x , T ) + ρ i u 0 ( x ) + 0 T ρ i φ ( x , t ) d t = . a . ψ i b ( x ) φ ( x , t ) + 0 T ρ i Q ( x , t ) d t
Using (7), we have
q ( x ) = A 2 ( x ) ( a 1 ( x ) + φ ¯ 1 ( x ) ) A 1 ( x ) ( a 2 ( x ) + φ ¯ 2 ( x ) )
u 0 ( x ) = B 2 ( x ) ( a 1 ( x ) + φ ¯ 1 ( x ) ) B 1 ( x ) ( a 2 ( x ) + φ ¯ 2 ( x ) )
where
a i ( x ) = . a ( x ) ψ i + 0 S Q ( x , t ) ρ i ( t ) d t A i ( x ) = ρ i ( 0 ) ψ 2 ρ 1 ( t ) ψ 1 ρ 2 ( t ) B i ( x ) = ψ i ( x ) ψ 2 ρ 1 ( t ) ψ 1 ρ 2 ( t ) φ ¯ i = ρ i φ ( x , T ) 0 T ρ i φ ( x , t ) d t
We introduce the following assumptions:
(a)
ψ 1 , ψ 2 H 2 ( χ ) L ( χ ) and ρ 1 , ρ 2 C 2 [ 0 , T ] ;
(b)
ψ 2 ρ 1 ( t ) ψ 1 ρ 2 ( t ) 0 ;
(c)
A 1 a 1 A 1 a 2 K 1 and B 2 a 1 B 1 a 2 K 2 a.e. in χ ¯ , for some positive constants K 1 and K 1 .
Inserting (8) and (9) into (1), we obtain
t t φ · ( a φ ) + A 2 ( x ) ( a 1 ( x ) + φ ¯ 1 ( x ) ) A 1 ( x ) ( a 2 ( x ) + φ ¯ 2 ( x ) ) φ = 0 in I , φ = 0 in ( R d χ ) × ( 0 , S ] , φ t ( · , 0 ) = 0 in χ , φ ( · , 0 ) = B 2 ( x ) ( a 1 ( x ) + φ ¯ 1 ( x ) ) B 1 ( x ) ( a 2 ( x ) + φ ¯ 2 ( x ) ) in χ ,
Thus, the solution to inverse problems (1), (3), and (4) is equivalent to obtaining the solution φ ( x , t ) to the nonlinear parabolic problem (10). Utilize the technique in [19], and consider the following two auxiliary hyperbolic problems:
t t Φ · ( a Φ ) + ( A 2 ( x ) a 1 ( x ) A 1 ( x ) a 2 ( x ) Φ = Q in I , Φ = 0 in ( R d χ ) × ( 0 , S ] , Φ t ( · , 0 ) = 0 in χ , Φ ( · , 0 ) = 0 in χ ,
and
t t Ψ · ( a Ψ ) + A 2 a 1 A 1 a 2 + c ( A 2 ( Ψ ¯ 1 + Φ ¯ 1 ) A 1 ( Ψ ¯ 2 + Φ ¯ 2 ) ) Ψ c A 2 ( Ψ ¯ 1 + Φ ¯ 1 ) A 1 ( Ψ ¯ 2 + Φ ¯ 2 ) Φ Φ = 0 in I , Ψ = 0 in ( R d χ ) × ( 0 , S ] , Ψ t ( · , 0 ) = 0 in χ , Ψ ( · , 0 ) = B 2 a 1 B 1 a 2 + c ( B 2 ( Ψ ¯ 1 + Φ ¯ 1 ) B 1 ( Ψ ¯ 2 + Φ ¯ 2 ) ) Ψ c B 2 ( Ψ ¯ 1 + Φ ¯ 1 ) B 1 ( Ψ ¯ 2 + Φ ¯ 2 ) Φ in χ ,
where c is a Lipschitz continuous function on R defined by
c ( x ) = x i f x K 0 K 0 i f x > K 0 K 0 i f x < K 0
and
Φ ¯ i = ρ i ( t ) Φ ( x , T ) 0 T ρ i Φ ( x , t ) d t , Ψ ¯ i = ρ i ( t ) Ψ ( x , T ) 0 T ρ i Ψ ( x , t ) d t
Therefore Φ , and Ψ , satisfy
Φ ( x , t ) L 2 a . e . ( x , t ) I , Ψ ( x , t ) L 3 a . e . ( x , t ) I ,
with L 1 > 0 .
To prove the existence and uniqueness of the solution to the inverse problems (1), (3), and (4). We set
K 3 : = ρ 1 ( t ) + ρ 2 ( t ) + 0 S ρ 1 d t + 0 S ρ 2 d t
K 4 = max x χ ¯ { A 1 , A 2 , B 1 , B 2 }
L 3 : = 2 K 3 K 4 ( L 1 + L 2 ) K 1 2 K 3 K 4 ( L 1 + L 2 )
Theorem 2. 
Let a , v 0 , ψ 1 , ψ 2 L ( χ ) , and ρ 1 , ρ 2 L ( 0 , S ) , Q L ( I ) , and, for the functional spaces inverse problem solution, we put Y 0 = L 2 0 , S ; H a 1 ( χ ) L 0 , S ; L 2 ( χ ) H 2 , 1 ( I ) .
Suppose that assumptions (a)–(c) are satisfied. Assume that there exists a number K 0 ( 0 , K 1 ) satisfying
2 K 3 K 4 ( L 1 + L 2 ) K 0 , L 3 < 1 .
Then, there exists at most one solution ( φ ( x , t ) , b ( x ) , u 0 ( x ) ) X 0 × L ( χ ) × L ( χ ) and b ( x ) > 0 a.e. x χ ¯ to the inverse problems (1), (3), and (4).
Proof. 
Using (13), we have
Φ ¯ 1 ρ 1 ( t ) + 0 T ρ 1 d t Φ L ( χ ) K 3 L 2
Likewise, we have Ψ ¯ 2 K 3 L 2 ,    Ψ ¯ 2 K 3 L 1 , which imply that
Φ ¯ i + Ψ ¯ i + K 3 ( L 1 + L 2 ) , f o r i = 1 , 2 .
A 2 ( Ψ ¯ 1 + Φ ¯ 1 ) A 1 ( Ψ ¯ 2 + Φ ¯ 2 ) K 3 ( L 1 + L 2 ) ( | A 1 | + | A 2 | ) 2 K 3 K 4 ( N 1 + N 2 ) ,
B 2 ( Ψ ¯ 1 + Φ ¯ 1 ) B 1 ( Ψ ¯ 2 + Φ ¯ 2 ) K 3 ( L 1 + L 2 ) ( | B 1 | + | B 2 | ) 2 K 3 K 4 ( N 1 + N 2 ) ,
Using inequality (17) and the definition of function c ( . ) , we obtain
c A 2 ( Ψ ¯ 1 + Φ ¯ 1 ) A 1 ( Ψ ¯ 2 + Φ ¯ 2 ) = A 2 ( Ψ ¯ 1 + Φ ¯ 1 ) A 1 ( Ψ ¯ 2 + Φ ¯ 2 ) ,
c B 2 ( Ψ ¯ 1 + Φ ¯ 1 ) B 1 ( Ψ ¯ 2 + Φ ¯ 2 ) = B 2 ( Ψ ¯ 1 + Φ ¯ 1 ) B 1 ( Ψ ¯ 2 + Φ ¯ 2 ) .
Hence, problem (12) becomes
t t Ψ · ( a Ψ ) + A 2 a 1 A 1 a 2 + A 2 ( Ψ ¯ 1 + Φ ¯ 1 ) A 1 ( Ψ ¯ 2 + Φ ¯ 2 ) Ψ A 2 ( Ψ ¯ 1 + Φ ¯ 1 ) A 1 ( Ψ ¯ 2 + Φ ¯ 2 ) Φ = 0 in I , Ψ = 0 in ( R d χ ) × ( 0 , S ] , Ψ t ( · , 0 ) = 0 in χ , Ψ ( · , 0 ) = B 2 a 1 B 1 a 2 + ( B 2 ( Ψ ¯ 1 + Φ ¯ 1 ) B 1 ( Ψ ¯ 2 + Φ ¯ 2 ) ) Ψ B 2 ( Ψ ¯ 1 + Φ ¯ 1 ) B 1 ( Ψ ¯ 2 + Φ ¯ 2 ) Φ in χ ,
By taking φ ( x , t ) = Ψ ( x , t ) + Φ ( x , t ) , it is easy to obtain that φ satisfies the estimate
φ ( x , t ) L 1 + L 2 a . e . ( x , t ) I ¯
Let w ( x , t ) and v ( x , t ) be two solutions to problem (10), and let w ( x , t ) = w ( x , t ) v ( x , t ) . Then, φ 3 ( x , t ) , satisfies the following problem:
t t z · ( a z ) + A 2 a 1 A 1 a 2 + A 2 ( w ¯ 1 A 1 ( w ¯ 2 z ( A 2 ( z ¯ 1 A 2 ( z ¯ 2 ) v = 0 in I , z = 0 in ( R d χ ) × ( 0 , S ] , z t ( · , 0 ) = 0 in χ , z ( · , 0 ) = B 2 ( z ¯ 1 B 1 ( z ¯ 2 in χ ,
with
w ¯ i = ρ i ( t ) w ( x , T ) 0 T ρ i w ( x , t ) , v ¯ i = ρ i ( t ) v ( x , T ) 0 T ρ i v ( x , t ) .
Using (17) and (18), we obtain
A 2 a 1 A 1 a 2 + A 2 w ¯ 1 A 1 w ¯ 2 K 1 max x χ ¯ A 2 w ¯ 1 A 1 w ¯ 2 K 1 2 K 3 K 4 ( L 1 + L 2 ) > K 0 2 K 3 K 4 ( L 1 + L 2 ) 0 ,
and, as z ¯ i K 3 z ¯ i L 2 ( I ) for i = 1 , 2 , then
max | A 2 z ¯ 1 A 1 z ¯ 2 | , | B 2 z ¯ 1 B 1 z ¯ 2 | 2 K 3 K 4 z L ( I )
z L ( I ) 2 K 3 K 4 z L ( I ) v L ( I ) K 1 2 K 3 K 4 ( L 1 + L 2 ) L 3 z L ( I ) .
Using the fact that L 3 < 1 , we obtain that z L ( I ) = 0 . This implies the uniqueness of the solution to problem (10). This means that the potential b ( x ) given by (8) and the initial condition u 0 ( x ) given by (9) only satisfy the inverse problems (1), (3), and (4). The proof is completed.    □
In this paper, we generalize the constant-time observations to mean-time records, which can be difficult to achieve in practice, and these records smooth out potentially large measurement errors in the direct solution φ . More precisely, we are looking for the triple ( φ ( x , t ) , b ( x ) , u 0 ( x ) ) of problem (1) with average weighted integral observations. Let φ ( x , t ; b , u 0 ) be the solution to the direct problem (1). Now, let us reformulate our inverse problem as an optimization problem. In reality, the average weighted integral observations ψ 1 , ψ 2 L ( χ ) , given in (3) and (4), respectively, may contain noise. Due to the poor posing of the inverse problem, which causes tiny inaccuracies in the input data (3) and (4) to lead to substantial errors in the output coefficients b ( x ) and u 0 ( x ) . This creates the main challenge numerically in the reconstruction of the solution. Numerically, we are looking for an approximation of the answer using measurements with noise. Let us consider ψ 1 δ , ψ 2 δ L ( χ ) , satisfying
ψ i δ ψ L 2 ( χ ) 2 δ , for i = 1 , 2 .
Hence, we will simultaneously reconstruct the coefficient b ( x ) and the condition u 0 ( x ) under the following noisy data: ( ψ 1 δ ( x ) , ( ψ 2 δ ( x ) )
0 S ρ 1 ( t ) φ ( x , t ) d t = ψ 1 δ ( x ) , x χ ,
0 S ρ 2 ( t ) φ ( x , t ) d t = ψ 2 δ ( x ) , x χ .
We minimize the objective functional M [ b , u 0 ] : E × L 2 ( χ ) R defined in (25) to obtain the solution to the inverse problem
M [ b , u 0 ] = 1 2 0 S ρ 1 ( t ) φ ( · , b , u 0 ) d t ψ 1 δ L 2 ( χ ) 2 + 1 2 0 S ρ 2 ( t ) φ ( · , b , u 0 ) d t ψ 2 δ L 2 ( χ ) 2 .
In this paper, our technique does not depend on the “regularize then discretize approach”; more precisely, we adopt the approach used in [32].
Theorem 3. 
The optimization problem in (25) admits at least one solution.
Proof. 
We know that inf E × L 2 ( χ ) M [ b , u 0 ] = : M 0 0 . This gives the existence of a minimizing sequence { ( b n , u 0 n ) : n N } E × L 2 such that
lim n M [ b n , u 0 n ] = M 0 .
Hence, the subsequence { ( b n , u 0 n ) : n N } is uniformly bounded in L ( χ ) × L 2 ( χ ) , which proves the existence of the subsequence ( b n , u 0 n ) which converges weakly to ( b * , u 0 * ) in L ( χ ) × L 2 ( χ ) . Using estimate (5), we find that the sequence { φ n = φ ( b n , u 0 n ) : n N } is uniformly bounded in H 0 1 ( I ) , So we can extract a subsequence { φ n : n N } , which converges weakly to φ * H 0 1 ( I ) in H 0 1 ( I ) .
Let θ H 0 1 ( I ) such that θ ( . , T ) = 0 . The variational formulation of problem (1) for ( φ n , b n , u 0 n ) gives
I φ n t t θ + I a φ n · θ + I b n φ n θ = I Q θ + I u 0 n t θ ( x , 0 ) .
We can write I b n φ n θ = I b * φ n θ + I ( b n b * ) φ n θ . Using estimate (5) and the fact that b n converges weakly to b * in L ( χ ) , we find that I ( b n b * ) φ n θ 0 . Consequently, by letting n tend to infinity in (26), we obtain
I φ * t t θ + I a φ * · θ + I b * φ * θ = I Q θ + I φ 0 * t θ ( x , 0 ) ,
and, since H 0 1 injects into L 2 ( I ) , and the solution to the direct problem is unique, we deduce that φ * = φ ( b * , u 0 * ) . Now, by the weak lower semi-continuity of the norm, we have
M [ b * , u 0 * ] = 1 2 0 S ρ 1 ( t ) φ ( · , t ) d t ψ 1 δ L 2 ( χ ) 2 + 1 2 0 S ρ 2 ( t ) φ ( · , t ) d t ψ 2 δ L 2 ( χ ) 2 1 2 lim n 0 S ρ 1 ( t ) φ n ( · , t ) d t ψ 1 δ L 2 ( χ ) 2 + 1 2 lim n 0 S ρ 2 ( t ) φ n ( · , t ) d t ψ 2 δ L L 2 ( χ ) 2 lim inf n M [ b n , u 0 n ] = inf E × L 2 ( χ ) M [ b , u 0 ] .
This shows that ( b * , u 0 * ) is a minimizer of (25) on the set E × L 2 ( χ ) .    □
To calculate the gradient of the functional M , we use the conjugate gradient method introduced in Section 3. First, we show that this functional M is differentiable as stated in the following lemma:
Lemma 1. 
Let the coefficient b and the initial condition be perturbed by a small δ b and δ u 0 , with ( b + δ b , u 0 + δ u 0 E × L 2 ( χ ) , and let φ be the weak solution to (1) corresponding to ( b , u 0 ) . The function E × L 2 ( χ ) H 0 1 ( I ) × H 0 1 ( I ) continues, i.e.,
φ ( b + δ b , u 0 ) φ ( b , u 0 ) L 2 ( 0 , T , H 0 1 ( χ ) ) c δ b L ( χ ) ,
φ ( b , u 0 + δ u 0 ) φ ( b , u 0 ) L 2 ( 0 , T , H 0 1 ( χ ) ) c δ u 0 L 2 ( χ ) .
Proof. 
The proof is very simple. Just use estimate (5) and replace φ with φ ( b + δ b , u 0 ) and φ ( b , u 0 + δ u 0 ) in (1).    □
We rely on the following lemma to show the differentiability of the map u φ ( b , u 0 ) .
Lemma 2. 
Assume that ( b , u 0 ) E × L 2 ( χ ) . The map ( b , u 0 ) u ( b , u 0 ) is Fréchet differentiable; more precisely, there are two operators O b , O u 0 : E × L 2 H 0 1 ( I ) × H 0 1 ( I ) , such that
lim δ b L ( χ ) 0 φ ( b + δ b , u 0 ) φ ( b , u 0 ) O b δ b L 2 ( 0 , T , H 0 1 ( χ ) ) δ b L ( χ ) = 0 ,
lim δ u 0 L 2 ( χ ) 0 φ ( b , u 0 + δ u 0 ) φ ( b , u 0 ) O u 0 δ u 0 L 2 ( 0 , T , H 0 1 ( χ ) ) δ b L 2 ( χ ) = 0 .
Proof. 
Let b , δ b E , and φ ( b , u 0 ) be the solution to problem (1). We can verify that φ satisfies the following problem:
t t φ b · ( a φ b ) + b φ b + δ b φ ( b , u 0 ) = 0 in I , φ b = 0 in ( R d χ ) × ( 0 , S ] , φ b ( · , 0 ) = 0 in χ , φ t ( · , 0 ) = 0 in χ ,
Then, problem (31) has a unique solution φ b which depends linearly on δ b . Let z b : = φ ( b + δ b , u 0 ) φ ( b , u 0 ) O b δ b = δ φ b φ b . So, δ φ b satisfies the following problem
t t δ φ b · ( a δ φ b ) + b δ φ b + δ b ( δ φ b + φ ( b , u 0 ) ) = 0 in I , δ φ b = 0 in ( R d χ ) × ( 0 , S ] , δ φ b ( · , 0 ) = 0 in χ , δ ( φ b ) t ( · , 0 ) = 0 in χ .
Using (31) directly shows that z verifies    □
t t z b · ( a z b ) + b z b + δ b δ φ b = 0 in I , z b = 0 in ( R d χ ) × ( 0 , S ] , z b ( · , 0 ) = 0 in χ , ( z b ) t ( · , 0 ) = 0 in χ .
Now, the application of estimate (5) to (33) gives
z b L 2 ( 0 , S , H a 1 ( χ ) ) c b , u 0 δ b δ φ b L 2 ( I ) c b , u 0 δ b L ( I ) δ φ b L 2 ( I ) c b , u 0 δ b L ( I ) δ φ b L 2 ( 0 , S , H a 1 ( χ ) ) .
From (27), we conclude that
φ ( b + δ b , u 0 ) φ ( b , u 0 ) O b δ b L 2 ( 0 , S , H a 1 ( χ ) ) c b , u 0 δ b L ( χ ) 2 .
Hence, the map ( b , u 0 ) φ ( b , u 0 ) is Fréchet differentiable. This shows relation (29). Now, let δ φ u 0 = φ ( b , u 0 + δ u 0 ) φ ( b , u 0 ) be the solution to the problem
t t φ u 0 · ( a φ u 0 ) + b φ u 0 + b δ φ b = 0 in I , φ u 0 = 0 in ( R d χ ) × ( 0 , S ] , δ φ u 0 ( · , 0 ) = 0 in χ , δ t φ u 0 ( · , 0 ) = 0 in χ .
In the same way, we can directly prove (30) for the coefficient b.
To calculate the gradient of M [ b , u 0 , we need the following lemma, which gives the existence and uniqueness of the following adjoint problem of (1)
t t p · ( a p ) + b p = k = 1 2 ρ k ( t ) 0 S ρ k φ ( x , τ ) d τ ψ k δ in I , p = 0 in ( R d χ ) × ( 0 , S ] , p ( · , 0 ) = 0 in χ , t p ( · , 0 ) = 0 in χ .
Lemma 3. 
Assume that ρ 1 , ρ 2 L ( χ ) and b L ( χ ) , ψ 1 , ψ 2 L ( 0 , S ) . Based on these assumptions, problem (35) admits a unique solution. Moreover, there is c > 0 such that
p L 2 ( I ) c k = 1 2 ρ k L ( 0 , S ) ψ k δ L ( χ ) + ρ k L ( 0 , S ) φ L ( I ) .
Proof. 
To show the existence and uniqueness of p, we use the same approach for problem (1). Multiplying both sides of the first equation of (35) by p t and integrating it over χ , we obtain
χ t t p t p χ · ( a ( p ) t p + χ b p = χ k = 1 2 ρ k ( t ) 0 S ρ k φ ( x , τ ) d τ ψ k δ t p .
Now, we integrate for t [ 0 , S ] , and we use the fact that p ( x , t ) = 0 t t p ( x , τ ) d τ + p ( x , 0 ) to obtain
p L 2 ( I ) 2 + I a ( φ ) 2 + I b p 2 c k = 1 2 I p ρ k ( t ) 0 S ρ k φ ( x , τ ) d τ ψ k δ .
since
I p ρ k ( t ) 0 S ρ k φ ( x , τ ) d τ ψ k δ c ψ k δ L ( χ ) ρ k L ( 0 , S ) + ρ k L ( 0 , S ) 2 φ L ( I ) p L 2 ( I ) ,
dividing by p L 2 ( I ) , we find that
p L 2 ( I ) c k = 1 2 ρ k L ( 0 , S ) ψ k δ L ( χ ) + ρ k L ( 0 , S ) φ L ( I ) .
this proves Lemma 3.    □
Theorem 4. 
The functional M [ b , u 0 ] is Fréchet differentiable, and its gradients are given respectively by
M b [ b , u 0 ] = 0 S φ ( x , t ) p ( x , t ) d t , x χ ,
M u 0 [ b , u 0 ] = p ( x , 0 ) , x χ ,
such that p is the solution to the adjoint problem (35).
Proof. 
Let δ b be a small perturbation of b such that b + δ b E , denoted by δ u b = φ ( b + δ b , u 0 ) φ ( b , u 0 ) . Then, δ b is the solution to the following problem
t t δ φ b · ( a δ φ b ) + b δ b δ φ b + δ b φ ( b + δ b , u 0 ) = 0 in I , δ φ b = 0 in ( R d χ ) × ( 0 , S ] , δ φ b ( · , 0 ) = 0 in χ , t δ φ b ( · , 0 ) = 0 in χ ,
Then, we have
M [ b + δ b , u 0 ) M [ b , u 0 ) = I δ u b ρ 1 ( t ) 0 S ρ 1 φ ( x , τ ) d τ ψ 1 δ + ρ 2 ( t ) 0 S ρ 2 φ ( x , τ ) d τ ψ 2 δ + 1 2 k = 1 2 0 S ρ k ( t ) δ φ b L 2 ( χ ) 2 = t t p · ( a p ) + b p + 1 2 k = 1 2 0 S ρ k ( t ) δ φ b ( . , τ ) d τ L 2 ( χ ) 2 .
Using problems (35) and (39) and integrating over I, we obtain
I δ u b t t p · ( a p ) + b p = I p t t δ φ b · ( a δ φ b ) + b δ φ b = I p δ b φ ( b + δ b , u 0 ) = I p δ b δ φ b I p δ b φ ( b , u 0 ) .
Hence,
M [ b + δ b , u 0 ) M [ b , u 0 ) = I p δ b δ φ b I p δ b φ ( b , u 0 ) + 1 2 k = 1 2 0 S ρ k ( t ) δ φ b ( . , τ ) d τ L 2 ( χ ) 2 .
According to an estimate similar to that of (5) for δ φ b , we have
0 S ρ k ( t ) δ φ b ( . , τ ) d τ L 2 ( χ ) 2 c ρ k ( t ) L ( 0 , S ) δ b L ( χ ) ,
and
I p δ b δ φ b c δ b L ( χ ) 2 p L 2 ( I ) .
Therefore,
M [ b + δ b , u 0 ) M [ b , u 0 ] = χ 0 S p φ d t d x .
Hence,
M b [ b , u 0 ] = 0 S φ ( x , t ) p ( x , t ) d t , x χ .
Using the same approach used to calculate the gradient with respect to b, we calculate the gradient with respect to the initial condition u 0 , such that the direct solution δ φ u 0 = φ ( b , u 0 + δ u 0 φ ( b , u 0 ) is the solution to the following problem
t t δ φ u 0 · ( a δ φ u 0 ) + b δ φ u 0 = 0 in I , δ φ u 0 = 0 in ( R d χ ) × ( 0 , S ] , δ φ u 0 ( · , 0 ) = δ u 0 ( x ) in χ , t δ φ u 0 ( · , 0 ) = 0 in χ ,
so we find (38).    □

3. Conjugate Gradient Method

To reconstruct simultaneously the coefficient b and the initial condition u 0 , we apply a process based on the conjugate gradient method (CGM), which consists of minimizing the functional in Equation (25). This process is given by:
b i + 1 = b i + β b i γ b i , u 0 i + 1 = u 0 i + β u 0 i γ u 0 i i = 0 , 1 , ,
where b 0 ( x ) , u 0 0 ( x ) are the initial guesses for b ( x ) , u 0 ( x ) where i is the number of iterations, and β i denotes the step size. Moreover, we define the search directions γ b i and γ u 0 i by
γ b i = M b 0 f o r i = 0 , M b i + ρ b i M b i 1 f o r i 1 , γ u 0 i = M u 0 0 f o r i = 0 , M u 0 i + ρ u 0 i M u 0 i 1 f o r i 1 ,
Additionally, we can use Fletcher–Reeves formula in [33] to obtain the step sizes β b i , β u 0 i
ρ b i = M b i L 2 ( χ ) M b i 1 L 2 ( χ ) , ρ u 0 i = M u 0 i L 2 ( χ ) M u 0 i 1 L 2 ( χ ) , f o r i = 1 , 2 , .
In the same way, we can find β b i and β u 0 i
M ( b i + 1 , u 0 i + 1 ) = min β b , β u 0 0 M ( b i + β b γ b i , u 0 i + β u 0 γ u 0 i ) .
According to the arguments used in [34], the fact that b ( b , u 0 ) , u 0 ( b , u 0 ) is Frécher differentiable, and the fact that the functional M ( b i + 1 , u 0 i + 1 ) is a monotone decreasing convergent sequence, we obtain the following result:
Theorem 5. 
The CGM (41)–(44) either terminates at a stationary point or converges in the following senses:
lim inf i b ( b , u 0 ) L 2 ( χ ) = lim inf i u 0 ( b , u 0 ) L 2 ( χ ) .
Since the errors can be noticed at the average weighted integral observations (3) and (4). The iteration process given by (42) cannot be preformed by the conjugate gradient method, and so the method is not well-posed because we do not have the regularization term. However, the method can become well-posed if we apply a divergence criterion so that the iteration procedure is stopped. This criterion is given by:
M ( b , u 0 ) 1 2 ψ δ ψ L 2 ( χ ) 2 δ 2 .
Thus, the iterations of this algorithm, based on CGM for the numerical reconstruction of the coefficient b ( x ) and the initial condition u 0 ( x ) , are as given by Algorithm 1.
Algorithm 1 CG algorithm for the minimizer of (25)
1.
Set i = 0 and initiate b 0 and u 0 0 for the coefficient b and initial condition u 0 .
2.
Determine numerically, using the finite difference method, the solution to the direct problem (1) φ ( x , t , b i , u 0 i ) and the objective functional (25).
3.
Determine numerically, using the finite difference method, the solution to the adjoint problem (35) p ( x , t , b i , u 0 i ) and the gradient of the objective functional (37) and (38). Calculate the coefficients ( ρ b i , ρ u 0 i ) and ( γ b i , γ u 0 i ) given in Equations (43) and (42), respectively.
4.
Determine numerically, using the finite difference method, the solution to the sensitivity problems (39) and (39) δ φ b ( x , t , b i , u 0 i ) and δ φ u 0 ( x , t , b i , u 0 i ) by using δ φ b i = γ b i and δ φ u 0 i = γ u 0 i , where the step sizeS β b i and β u 0 i are given in (44)
5.
Update b i , u 0 i using (41).
6.
If the condition (46) is satisfied, go to Step 7. Else, set i = i + 1 and return to Step 2.
7.
End.

4. Numerical Experiments

Here, we apply the CG method [35] in one and two dimensions in order to identify simultaneously the coefficient b and the initial condition u 0 in 1. We discretize this problem using the finite element method in space and the Crank–Nicholson scheme in time direction. For the two-dimensional problem, we apply the alternating direction implicit (ADI) method as described in [28]. We create the noisy data, and a random perturbation is added, i.e.,
ψ 1 δ ( x ) = ψ 1 + η × r a n d o m ( 1 ) , ψ 2 δ ( x ) = ψ 2 + η × r a n d o m ( 1 ) ,
where η = p / 100 × max x χ { ψ 1 , ψ 2 } , and p represents the percentage of noise.
We calculate the approximate L 2 error by the following formula to demonstrate the precision of the numerical solution
E b i = b i b ex L 2 ( χ ) ,
E u 0 = u 0 i u 0 ex L 2 ( χ ) ,
where ( b i , u 0 i ) are the initial guesses reconstructed at the kth iteration, and ( b ex , u 0 ex ) are the exact values. The residual R i at the ith iteration is given by
R 1 i = u ( b i ) ψ 1 δ L 2 ( χ ) , R 2 i = u ( u 0 i ) ψ 2 δ L 2 ( χ ) .

4.1. One-Dimensional Problem

We fix χ = [ 0 , 1 ] without a loss of generality. We also fix δ t = 1 100 and δ x = 1 100 , respectively, and let a ( x ) = ( x 0.5 ) 2 . We choose ρ 1 ( t ) and ρ 2 ( t ) in (3) and (4) as
ρ i ( t ) = 1 q π exp ( t t i ) 2 10 6 t [ 0 , S ] , i = 1 , 2 ,
where q is a small positive constant, and t 1 t 2 [ 0 , S ] . It is clear that ρ i ( t ) δ ( t t i ) for small values of q, where δ ( . ) is the Dirac delta function. Then, according to the properties of the Dirac delta function, Equations (3) and (4)would become
ψ i = 0 S ρ i ( t ) φ ( x , t ) 0 S δ ( t t i ) φ ( x , t ) = φ ( x , t ) .
For all the three numerical examples presented later, we choose the weight functions as
ρ 1 ( t ) = 1000 π exp ( t 0.3 ) 2 10 6 , ρ 2 ( t ) = 1000 π exp ( t 1 ) 2 10 6 .
Example 1. 
To validate this choice, we apply our proposed algorithm for the reconstruction of the coefficient and the initial condition defined on χ by
b ex ( x ) = 1 + x 2 , u 0 ex ( x ) = x ( 1 x ) .
We take the initial guesses b 0 ( x ) = 1 , u 0 0 ( x ) = 0 .
Figure 1 is devoted to showing the variation of the functional according to the number of iterations i for the simultaneous determination of the coefficient b and the initial condition u 0 , in the cases of no noise ( p = 0 ) and with noise ( p = 1 ) , ( p = 5 ) . It can be observed that if p = 0 , then the function M ( b , u 0 ) quickly converges to a very small value because it is a decreasing function according to the values of i. We see that the number of iterations for the algorithm to stop is i * = 20 in the case of no noise and i * = 3 in the case of noise. These numbers of iterations are obtained via the discrepancy in (46). Errors (47) and (48) associated with b ( x ) and u 0 ( x ) have been found for E b = 0.00008 / 0.0284 / 0.0274 and E u 0 = 0.000025 / 0.00322 / 0.00223 with the noise levels p = 0 / 1 / 5 , respectively. Therefore, we can conclude that our numerical process is reasonably accurate in determining the coefficient b ( x ) and the initial condition u 0 ( x ) . Similarly, the norms of the gradients of the functional are obtained for M b L 2 ( χ ) = 2 × 10 5 / 0.0009 / 0.0014 and M u 0 3 L ( 2 ( χ ) = 1.64 × 10 5 / 0.00018 / 0.00015 with p = 0 / 1 / 5 .
In Figure 2, we illustrate the comparison between the functions of the recovered coefficient b and the initial condition u 0 and their exact values with the noise levels p = 0 / 1 / 5 .
Figure 2 illustrates the numerical reconstruction of the coefficient b and the initial condition u 0 with the number of final iterations i * given in Figure 1 with different values for the noise levels, p = 0 / 1 / 5 . From Figure 1, we can notice that the exact and the numerical solutions are almost equal; that is, the proposed algorithm to determine the two coefficients converges to the required solutions.
In the following example, we will apply our method to reconstruct the coefficient b and the initial condition defined by
Example 2. 
We consider
b ex ( x ) = x 3 + sin ( π x ) , u 0 ex ( x ) = sin ( π x ) .
Similarly, for Example 1, using the discrepancy in (46), the stopping iterations are, i * = 20 / 3 / 3 for the noise levels p = 0 , 1 , and 5. The errors (47) and (48) associated with b ( x ) and u 0 ( x ) have been found for E b = 0.000025 / 0.0081 / 0.0092 and E u 0 = 0.00018 / 0.00254 / 0.00353 with noise level p = 0 , 1 and 5, respectively, Under these stopping iterations, Figure 3 illustrates the exact and numerical solutions for Example 2 for different values of the noise levels ( p = 0 , 1 , and 5).

4.2. Two-Dimensional Problem

The space–time region ¯ χ × [ 0 , S ] : = [ 0 , 1 ] 2 × [ 0 , 1 ] is partitioned into 40 × 40 × 80 equidistant meshes. We testify the numerical performance of the Algorithm 1 for the reconstruction of the coefficient b and the initial condition b with the following degeneracy
a ( x ) : = a ( x 1 , x 2 ) = ( x 1 0.3 ) 2 + ( x 2 0.3 ) 2 .
We take the initial guesses b 0 ( x 1 , x 2 ) = 1 , u 0 0 ( x 1 , x 2 ) = 0 .
Example 3. 
Suppose that the exact coefficient and initial condition for the degenerate wave problem (1) are given by:
b ex ( x 1 , x 2 ) = cos ( 2 π x 1 ) cos ( 3 x 2 ) , u 0 ex ( x 1 , x 2 ) = ( 1 x 2 ) sin ( π x 1 ) .
For this example, the number of the final iterations of CGM is i * = 60 / 3 / 3 , for the noise levels p = 0 , 1 , and 5 (Figure 4). Errors (47) and (48), associated with b ( x ) and u 0 ( x ) , have been found for E b = 0.01 , 0.05 , and 0.21 . Additionally, E u 0 = 0.051 , 0.154 , and 0.254 .
The exact coefficient, initial condition functions, and recovered solution are shown in Figure 5 and Figure 6. The absolute errors between the exact coefficient and initial condition functions and their numerical reconstruction are shown in Figure 7 and Figure 8. We can notice that the recovered terms are very close to the exact solutions; this shows the effectiveness of our proposed method.
The following Table 1 and Table 2 are devoted to the different values of E b i and E u 0 i for p = 0 , p = 1 , and p = 5 .
Here, we summarize the results of the numerical experiments. For one-dimensional problems, as considered in Example 1, the resulting figures can be described as follows. Figure 1a represents the graph of the functional M without a noise level (0 percent) and with the iteration number i * = 20 , while Figure 1b represents the graph of the functional M with a noise level (1 percent) with the number of iterations i * = 3 , Figure 1c represents the graph of the functional M with a noise level (5 percent) with the number of iterations i * = 3 .  Figure 2a represents the graph of the reaction coefficient b without a noise level (0 percent) and with a noise level (1 and 5 percent). Additionally, Figure 2b represents the graph of the initial condition u 0 without a noise level (0 percent) and with a noise level (1 and 5 percent). For Example 2, Figure 3a represents the graph of the reaction coefficient b without a noise level (0 percent) and with a noise level (1 and 5 percent). Moreover, Figure 3b represents the graph of the initial condition u 0 without a noise level (0 percent) and with a noise level (1 and 5 percent).
For two-dimensional problems, as given in Example 3, Figure 4a represents the graph of the functional M without a noise level (0 percent) and with the iteration number i * = 60 . Figure 4b represents the graph of the functional J with a noise level (1 percent) and with the number of iterations i * = 3 . Figure 4c represents the graph of the functional J with a noise level (5 percent) and with the number of iterations i * = 3 .  Figure 5 represents the exact solution to the recovered reaction coefficient. Figure 9a, Figure 9b, and Figure 9c represent the graphs of the reaction coefficient b without a noise level (0 percent) and with noise levels (1 and 5 percent), respectively. Figure 7a, Figure 7b, and Figure 7c represent the graphs of the absolute error between the exact and numerical reaction coefficient b without a noise level (0 percent) and with noise levels (1 and 5 percent), respectively. Figure 6 gives the exact solution to the initial condition. Figure 10a, Figure 10b, and Figure 10c represent the graph of the reaction coefficient b without a noise level (0 percent) and with noise levels (1 and 5 percent), respectively. Figure 8a, Figure 8b, and Figure 8c represent the graphs of the absolute error between the exact and numerical reaction coefficient b without a noise level (0 percent) and with noise levels (1 and 5 percent), respectively.

5. Conclusions

In this article, we were concerned with the simultaneous determination of the reaction coefficient, which depends on space, and the initial condition of a hyperbolic problem with degeneracy within the spatial domain from temporal integral observations. The existence, uniqueness, and stability of the inverse problem are examined. The CG conjugate gradient method has been proposed with adjoint and sensitivity problems for simultaneously reconstruction of the two unknown functions by minimizing the objective least squares functional. In the simulation part, we consider some numerical experiments in which we reconstruct numerically the initial condition and potential for wave problems in one and two dimensions.

Author Contributions

Conceptualization, H.O.S., M.A.Z., R.H.D.S., and A.S.H.; Formal analysis, H.O.S., M.A.Z., R.H.D.S., and A.S.H.; Funding acquisition, R.H.D.S. and A.S.H.; Methodology, H.O.S., M.A.Z., R.H.D.S., and A.S.H.; Supervision, R.H.D.S. and A.S.H.; Visualization, H.O.S. and M.A.Z.; Writing—original draft, H.O.S., M.A.Z., and A.S.H.; Writing—review and editing, R.H.D.S. All authors have read and agreed to the published version of the manuscript.

Funding

Ahmed Hendy wishes to acknowledge the RSF, Russia grant, project 22-21-00075.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the respected reviewers for their valuable comments which significantly improved the paper.

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. Frisch, U.; Matarrese, S.; Mohayaee, R.; Sobolevski, A. A reconstruction of the initial conditions of the universe by optimal mass transportation. Nature 2002, 417, 260–262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kalnay, E. Atmospheric Modeling, Data Assimilation and Predictability; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  3. Rodgers, C.D. Inverse Methods for Atmospheric Sounding: Theory and Practice; World Scientific: Singapore, 2000; Volume 2. [Google Scholar]
  4. Kuchment, P.; Kunyansky, L. Mathematics of thermoacoustic tomography. Eur. J. Appl. Math. 2008, 19, 191–224. [Google Scholar] [CrossRef] [Green Version]
  5. Gouveia, W.P.; Scales, J.A. Bayesian seismic waveform inversion: Parameter estimation and uncertainty analysis. J. Geophys. Res. Solid Earth 1998, 103, 2759–2779. [Google Scholar] [CrossRef]
  6. Beauchard, K.; Zuazua, E. Some controllability results for the 2D Kolmogorov equation. In Annales de l’IHP Analyse Non Linéaire; Elsevier: Paris, France, 2009; Volume 26, pp. 1793–1815. [Google Scholar]
  7. Buchot, J.M.; Raymond, J.P. A linearized model for boundary layer equations. In Proceedings of the Optimal Control of Complex Structures: International Conference in Oberwolfach, 4–10 June 2000; Springer: Berlin/Heidelberg, Germany, 2002; pp. 31–42. [Google Scholar]
  8. Shimakura, N. Partial Differential Operators of Elliptic Type; Translated from Japanese by N. Shimakura; American Mathematical Society: Providence, RI, USA, 1992. [Google Scholar]
  9. Zaky, M.A.; Taha, T.R.; Suragan, D.; Hendy, A.S. An L1 type difference/Galerkin spectral scheme for variable-order time-fractional nonlinear diffusion–reaction equations with fixed delay. J. Comput. Appl. Math. 2023, 420, 114832. [Google Scholar] [CrossRef]
  10. Hendy, A.S.; De Staelen, R.H.; Aldraiweesh, A.A.; Zaky, M.A. High order approximation scheme for a fractional order coupled system describing the dynamics of rotating two-component Bose-Einstein condensates. AIMS Math. 2023, 8, 22766–22788. [Google Scholar] [CrossRef]
  11. Ameen, I.G.; Zaky, M.A.; Doha, E.H. Singularity preserving spectral collocation method for nonlinear systems of fractional differential equations with the right-sided Caputo fractional derivative. J. Comput. Appl. Math. 2021, 392, 113468. [Google Scholar] [CrossRef]
  12. Hendy, A.S.; Van Bockstal, K. On a reconstruction of a solely time-dependent source in a time-fractional diffusion equation with non-smooth solutions. J. Sci. Comput. 2022, 90, 41. [Google Scholar] [CrossRef]
  13. Hendy, A.S.; Van Bockstal, K. A solely time-dependent source reconstruction in a multiterm time-fractional order diffusion equation with non-smooth solutions. Numer. Algorithms 2022, 90, 809–832. [Google Scholar] [CrossRef]
  14. Sidi, H.O.; Zaky, M.A.; Qiu, W.; Hendy, A.S. Identification of an unknown spatial source function in a multidimensional hyperbolic partial differential equation with interior degeneracy. Appl. Numer. Math. 2023, 192, 1–18. [Google Scholar] [CrossRef]
  15. Petersson, N.A.; Sjögreen, B. Super-grid modeling of the elastic wave equation in semi-bounded domains. Commun. Comput. Phys. 2014, 16, 913–955. [Google Scholar] [CrossRef] [Green Version]
  16. Tarantola, A. Inversion of seismic reflection data in the acoustic approximation. Geophysics 1984, 49, 1259–1266. [Google Scholar] [CrossRef]
  17. Isakov, V. Inverse parabolic problems with the final overdetermination. Commun. Pure Appl. Math. 1991, 44, 185–209. [Google Scholar] [CrossRef]
  18. Prilepko, A.I.; Kostin, A.B. On certain inverse problems for parabolic equations with final and integral observation. Matematicheskii Sbornik 1992, 183, 49–68. [Google Scholar] [CrossRef]
  19. Kozhanov, A.I. A Nonlinear Loaded Parabolic Equation and a Related Inverse Problem. Math. Notes 2004, 76, 784–795. [Google Scholar] [CrossRef]
  20. Kamynin, V.; Kostin, A. Two inverse problems of finding a coefficient in a parabolic equation. Differ. Equ. 2010, 46, 375–386. [Google Scholar] [CrossRef]
  21. Deng, Z.C.; Yang, L.; Yu, J.N. Identifying the radiative coefficient of heat conduction equations from discrete measurement data. Appl. Math. Lett. 2009, 22, 495–500. [Google Scholar] [CrossRef] [Green Version]
  22. Chen, Q.; Liu, J. Solving an inverse parabolic problem by optimization from final measurement data. J. Comput. Appl. Math. 2006, 193, 183–203. [Google Scholar] [CrossRef] [Green Version]
  23. Trucu, D.; Ingham, D.B.; Lesnic, D. Space-dependent perfusion coefficient identification in the transient bio-heat equation. J. Eng. Math. 2010, 67, 307–315. [Google Scholar] [CrossRef]
  24. Cao, K.; Lesnic, D. Reconstruction of the space-dependent perfusion coefficient from final time or time-average temperature measurements. J. Comput. Appl. Math. 2018, 337, 150–165. [Google Scholar] [CrossRef]
  25. Erdem, A.; Lesnic, D.; Hasanov, A. Identification of a spacewise dependent heat source. Appl. Math. Model. 2013, 37, 10231–10244. [Google Scholar] [CrossRef]
  26. Johansson, B.T.; Lesnic, D. A procedure for determining a spacewise dependent heat source and the initial temperature. Appl. Anal. 2008, 87, 265–276. [Google Scholar] [CrossRef]
  27. Yamamoto, M.; Zou, J. Simultaneous reconstruction of the initial temperature and heat radiative coefficient. Inverse Probl. 2001, 17, 1181. [Google Scholar] [CrossRef]
  28. Cao, K.; Lesnic, D.; Liu, J. Simultaneous reconstruction of space-dependent heat transfer coefficients and initial temperature. J. Comput. Appl. Math. 2020, 375, 112800. [Google Scholar] [CrossRef]
  29. Cao, K.; Lesnic, D. Simultaneous identification and reconstruction of the space-dependent reaction coefficient and source term. J. Inverse Ill-Posed Probl. 2021, 29, 867–894. [Google Scholar] [CrossRef]
  30. Alifanov, O.M. Inverse Heat Transfer Problems; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  31. Atifi, K.; Essoufi, E.H.; Khouiti, B. New approach to identify the initial condition in degenerate hyperbolic equation. Inverse Probl. Sci. Eng. 2019, 27, 484–512. [Google Scholar] [CrossRef]
  32. Engl, H.W.; Groetsch, C.W. Projection-regularization methods for linear operator equations of the first kind. In Special Program on Inverse Problems; Australian National University, Mathematical Sciences Institute: Canberra, Australia, 1988; Volume 17, pp. 17–32. [Google Scholar]
  33. Fletcher, R.; Reeves, C.M. Function minimization by conjugate gradients. Comput. J. 1964, 7, 149–154. [Google Scholar] [CrossRef] [Green Version]
  34. Dai, Y.H.; Yuan, Y.x. Convergence properties of the Fletcher-Reeves method. Ima J. Numer. Anal. 1996, 16, 155–164. [Google Scholar] [CrossRef]
  35. Nazareth, J.L. Conjugate gradient method. Wiley Interdiscip. Rev. Comput. Stat. 2009, 1, 348–353. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Graph of the functional M ( b i , u 0 i ) to Example 1 with (a) p = 0 ; (b) p = 1 ; (c) p = 5 .
Figure 1. Graph of the functional M ( b i , u 0 i ) to Example 1 with (a) p = 0 ; (b) p = 1 ; (c) p = 5 .
Mathematics 11 03186 g001aMathematics 11 03186 g001b
Figure 2. Reconstruction results for Example 1 with p = 0 , 1 , and 5. (a) coefficient b; (b) the initial condition u 0 .
Figure 2. Reconstruction results for Example 1 with p = 0 , 1 , and 5. (a) coefficient b; (b) the initial condition u 0 .
Mathematics 11 03186 g002
Figure 3. Reconstruction results for Example 2 with p = 0 , 1 , and 5. (a) coefficient b; (b) the initial condition u 0 .
Figure 3. Reconstruction results for Example 2 with p = 0 , 1 , and 5. (a) coefficient b; (b) the initial condition u 0 .
Mathematics 11 03186 g003aMathematics 11 03186 g003b
Figure 4. Graph of the functional M ( b i , u 0 i ) for Example 3 with (a) p = 0 ; (b) p = 1 ; (c ) p = 5 .
Figure 4. Graph of the functional M ( b i , u 0 i ) for Example 3 with (a) p = 0 ; (b) p = 1 ; (c ) p = 5 .
Mathematics 11 03186 g004aMathematics 11 03186 g004b
Figure 5. True coefficient function (53) for Example 3.
Figure 5. True coefficient function (53) for Example 3.
Mathematics 11 03186 g005
Figure 6. True initial condition function (53) for Example 3.
Figure 6. True initial condition function (53) for Example 3.
Mathematics 11 03186 g006
Figure 7. Absolute error between exact and numerical coefficient b with different noise levels (a) p = 0 ; (b) p = 1 ; (c) p = 5 for Example 3.
Figure 7. Absolute error between exact and numerical coefficient b with different noise levels (a) p = 0 ; (b) p = 1 ; (c) p = 5 for Example 3.
Mathematics 11 03186 g007
Figure 8. Absolute error between exact and numerical initial condition u 0 with different noise levels (a) p = 0 ; (b) p = 1 ; and (c) p = 5 for Example 3.
Figure 8. Absolute error between exact and numerical initial condition u 0 with different noise levels (a) p = 0 ; (b) p = 1 ; and (c) p = 5 for Example 3.
Mathematics 11 03186 g008
Figure 9. Numerical coefficient function reconstruction results for Example 3 with different noise levels (a) p = 0 ; (b) p = 1 ; (c) p = 5 .
Figure 9. Numerical coefficient function reconstruction results for Example 3 with different noise levels (a) p = 0 ; (b) p = 1 ; (c) p = 5 .
Mathematics 11 03186 g009aMathematics 11 03186 g009b
Figure 10. Numerical initial condition function reconstruction results for Example 3 with different noise levels (a) p = 0 ; (b) p = 1 ; (c) p = 5 .
Figure 10. Numerical initial condition function reconstruction results for Example 3 with different noise levels (a) p = 0 ; (b) p = 1 ; (c) p = 5 .
Mathematics 11 03186 g010aMathematics 11 03186 g010b
Table 1. Numerical results (accuracy error E b i ).
Table 1. Numerical results (accuracy error E b i ).
p = 0 p = 1 p = 5
Example 10.000080.02840.0377
Example 30.0000250.003020005040
Example 30.010.050.21
Table 2. Numerical results (accuracy error E u 0 i ).
Table 2. Numerical results (accuracy error E u 0 i ).
p = 0 p = 1 p = 5
Example 10.000250.003220.00524
Example 30.002210.001020.00417
Example 30.0510.02540.0354
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sidi, H.O.; Zaky, M.A.; De Staelen, R.H.; Hendy, A.S. Numerical Reconstruction of a Space-Dependent Reaction Coefficient and Initial Condition for a Multidimensional Wave Equation with Interior Degeneracy. Mathematics 2023, 11, 3186. https://doi.org/10.3390/math11143186

AMA Style

Sidi HO, Zaky MA, De Staelen RH, Hendy AS. Numerical Reconstruction of a Space-Dependent Reaction Coefficient and Initial Condition for a Multidimensional Wave Equation with Interior Degeneracy. Mathematics. 2023; 11(14):3186. https://doi.org/10.3390/math11143186

Chicago/Turabian Style

Sidi, Hamed Ould, Mahmoud A. Zaky, Rob H. De Staelen, and Ahmed S. Hendy. 2023. "Numerical Reconstruction of a Space-Dependent Reaction Coefficient and Initial Condition for a Multidimensional Wave Equation with Interior Degeneracy" Mathematics 11, no. 14: 3186. https://doi.org/10.3390/math11143186

APA Style

Sidi, H. O., Zaky, M. A., De Staelen, R. H., & Hendy, A. S. (2023). Numerical Reconstruction of a Space-Dependent Reaction Coefficient and Initial Condition for a Multidimensional Wave Equation with Interior Degeneracy. Mathematics, 11(14), 3186. https://doi.org/10.3390/math11143186

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