Next Article in Journal
Special Issue on Algorithms for PID Controllers 2021
Next Article in Special Issue
Local Convergence Analysis of a One Parameter Family of Simultaneous Methods with Applications to Real-World Problems
Previous Article in Journal
Sybil in the Haystack: A Comprehensive Review of Blockchain Consensus Mechanisms in Search of Strong Sybil Attack Resistance
Previous Article in Special Issue
Some Inexact Version of the Iterative Exponential Method for Solving Nonsmooth Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solving of the Inverse Boundary Value Problem for the Heat Conduction Equation in Two Intervals of Time

by
Bashar Talib Al-Nuaimi
1,
H.K. Al-Mahdawi
2,
Zainalabideen Albadran
3,
Hussein Alkattan
4,
Mostafa Abotaleb
4,* and
El-Sayed M. El-kenawy
5
1
Computer Science Department, University of Diyala, Diyala 32001, Iraq
2
Electronic Computer Centre, University of Diyala, Diyala 32001, Iraq
3
Computer Technical Engineering Department, University of Alkafeel, Najaf 60061, Iraq
4
Department of System Programming, South Ural State University, 454080 Chelyabinsk, Russia
5
Department of Communications and Electronics, Delta Higher Institute of Engineering and Technology, Mansoura 35111, Egypt
*
Author to whom correspondence should be addressed.
Algorithms 2023, 16(1), 33; https://doi.org/10.3390/a16010033
Submission received: 23 November 2022 / Revised: 23 December 2022 / Accepted: 30 December 2022 / Published: 6 January 2023

Abstract

:
The boundary value problem, BVP, for the PDE heat equation is studied and explained in this article. The problem declaration comprises two intervals; the (0, T) is the first interval and labels the heating of the inside burning chamber, and the second (T, ∞) interval defines the normal cooling of the chamber wall when the chamber temperature concurs with the ambient temperature. It is necessary to prove the boundary function of this problem has its place in the space H 1 0 , in order to successfully apply the Fourier transform method. The applicability of the Fourier transform for time to this problem is verified. The method of projection regularization is used to solve the inverse boundary value problem for the heat equation and to obtain an evaluation for the error between the approximate and the real solution. These results are new and of practical interest as shown in the numerical case study.

1. Introduction

A heat speared problem handles the estimate of unidentified numbers appearing in the mathematics of physical in thermal knowledges, by means of the dimensions or measurement of the temperature, radiation intensities, heat flux, etc.
The inverse problem for the heat PDE system can be solved by many methods; for example, the method of Tikhonov [1], the method of Lavrentiev [2], Ivanov [3], and many others. The inverse problems in the heat PDE system can be grouped as two types depending on the e unknown function or vector for the initial part or the boundary part conditions, and many studies of these problems are considered in many works [4,5,6,7,8,9,10,11,12]. Various methods for solving this type of inverse problem have been proposed in many works [13,14,15,16,17]. In the article [13], the BVP for the PDE heat equation in a hollow cylinder was solved by using the Fourier projection method. Papers [14,16] studied the multigrid method with the iterative method to find the solution for the inverse problem, IP, in the heat PDE system. In [15,17], the iterative methods with necessary analyses were studied for solving the inverse linear operator equation and the case study in this paper was the inverse heat PDE system problem.
The successfully accomplished approaches for resolving the IPs are dependent, to a large degree, on the deep insight into the mathematical problems related to the algorithms and statements and the definition of the specific difficulties in their solving [18,19,20,21,22,23].
The goal of this article is to provide the approximation solution for the BVP in the PDE for the heat equation system with the mixed interval for time. Hence, the result of this problem (BVP) is not contingent continuously on the known data in the field, which means the solution is not stable; therefore, this problem is known as an ill-posed inverse problem. The proving of the boundary function of this problem belonged to the class H 1 0 , necessary for applying the projection regularization method by using the Fourier transform. For solving the ill-posed problems, a central role is played by the error estimations between the approximation and real solutions. We obtain the estimate solution by applying the projection regularization method with the Fourier transform, making these results new and interesting.

2. Materials and Methods Direct Formulation of the Problem on Interval 0 , T

We considered the case of the heat equation on a segment with inhomogeneous boundary data.
u 1 x , t t = 2 u 1 x , t x 2 ,   0 < x < 1 ,   t 0 , T ,
u 1 x , 0 = 0 ,   0 x 1 ,
u 1 0 , t = 0 ,   0 t T ,
u 1 1 , t = u t ,   0 t T ,
Assume the q t function is defined as the following
q 1 t C 3 0 , T , q 1 0 = q 1 0 = q 1 0 = q 1 T = q 1 T = 0 ,
by using Duhamel’s principle method ([24], p. 109)
u 1 ( x , t ) = 0 t w t ( x , t τ ) q 1 ( τ ) d τ ,
integration by parts for the right part for (6) once, we obtain
u 1 ( x , t ) = 0 t w ( x , t τ ) q 1 ( τ ) d τ + q 1 ( 0 ) w ( x , t ) .
Now, we can decide to obtain the solution for w ( x , t ) as the following
w ( x , t ) = x + 2 n = 1 1 n n π e ( n π ) 2 t sin ( n π x ) ,
by substituting (8) in (7), to obtain a solution to a non-stationary problem, from (5) q 1 0 = 0 ,
u 1 ( x , t ) = x + n = 1 1 n C n ( t ) sin ( n π x ) ,
C n ( t ) = 2 n π 0 t e ( n π ) 2 ( t τ ) q 1 ( τ ) d τ ,
where x 0 , 1 , and t 0 , T .
Lemma 1.
Let q 1 ( t ) satisfy condition (5). Then, there exists a solution u 1 ( x , t ) for problem (1)–(5) such that u 1 ( x , t ) satisfies the Equation (1) on the set 0 , 1 × 0 , T , initial condition (2), boundary conditions (3), (4) and u 1 x , t C 0 , 1 × 0 , T C 2 , 1 0 , 1 × 0 , T .
Proof. 
By integrating the right side of the Formula (10) in parts twice, we obtain
C n ( t ) = 2 n π 5 0 t e ( n π ) 2 ( t τ ) q 1 ( τ ) d τ + 2 n π 3 q 1 ( t ) 2 n π 5 q 1 ( t ) ,
since C n ( t ) sin ( n π x ) C n ( t ) for any n and from the Cauchy–Bunyakovsky inequality
0 t e ( n π ) 2 ( t τ ) q 1 ( τ ) d τ q 1 ( τ ) L 2 ( 0 , T ] 1 2 n π ,
by means of (5), (4) and (12) for any t 0 , T and for any n we obtain
C n ( t ) 2 q 1 ( t ) L 2 n π 6 + 2 n π 3 q 1 ( t ) 2 n π 5 q 1 ( t ) .
Using Equations (11)–(13) and convergence of the series n 1 1 n 3 , n 1 1 n 5 , n 1 1 n 6 with the Weierstrass criterion follows the unchanging convergence of the above series on 0 , 1 × 0 , T .
Since the functions e ( n π ) 2 ( τ ) q 1 ( τ ) L 2 0 , T , obtaining
0 t e ( n π ) 2 ( t τ ) q 1 ( τ ) d τ = e ( n π ) 2 t 0 t e ( n π ) 2 τ q 1 ( τ ) d τ C 0 , T .
Thus, with q 1 ( t ) , q 1 ( t ) and q 1 ( t ) C 0 , T in addition to Equations (11) and (14), we take C n ( t ) C 0 , T . From this condition and the convergence of (19) in domain 0 , 1 × 0 , T , we have u 1 ( x , t ) C ( 0 , 1 × 0 , T ) . Differentiating a C n ( t ) sin ( n π x ) with x and by using (13), we obtain
C n ( t ) sin n π x x 2 q 1 ( t ) L 2 n π 5 + 2 n π 3 q 1 ( t ) 2 n π 4 q 1 ( t ) .
From the above relation, we obtain the convergence of the n 0 C n ( t ) sin n π x x in 0 , 1 × 0 , T , from (8) we have u 1 ( x , t ) x = n 0 C n ( t ) sin n π x x in 0 , 1 × 0 , T and u 1 ( x , t ) x ( 0 , 1 × 0 , T ) .
Now, let us examine the function 2 u 1 ( x , t ) x 2 .
Differentiating the function C n ( t ) sin ( n π x ) by x twice and using (11), we obtain C n ( t ) sin ( n π x ) x x = 2 sin ( n π x ) n π 3 0 t e ( n π ) 2 ( t τ ) q 1 ( τ ) d τ + 2 sin ( n π x ) n π q 1 ( t ) 2 sin ( n π x ) n π 3 q 1 ( t ) , since the number series n 1 1 n 3 , converge according to the Weierstrass criterion, the functional series 2 n = 0 sin ( n π x ) n π q 1 ( t ) , converge absolutely and uniformly on 0 , 1 × 0 , T .
Then, we need to check the convergence for n = 0 sin ( n π x ) n π to any ε > 0 in this series, related to the Dirichlet criterion, the convergence is consistently on ε , 1 ε .
Meanwhile, any ε > 0 series n 0 C n ( t ) sin ( n π x ) x x converges on ε , 1 ε × 0 , T and the parts of this series are nonstop, we obtain
2 u 1 ( x , t ) x 2 = n 0 C n ( t ) sin n π x x x , 2 u 1 ( x , t ) x 2 C ( 0 , 1 ) × 0 , T .
The lemma is proofed. □
Now, let us examine the function u 1 ( x , T ) .
Lemma 2.
Function u 1 ( x , T ) , defined by formulas (9) and (11), belongs to space H 4 0 , 1 .
Proof. 
From (5), (9) and (11) it follows that
u 1 ( x , T ) = x + n = 1 1 n C n ( T ) sin ( n π x ) ,
where
C n ( T ) = 2 n π 5 0 T e ( n π ) 2 ( T τ ) q 1 ( τ ) d τ + 2 n π 3 q 1 ( T ) 2 n π 5 q 1 ( T ) .
Since the conditions
u 1 ( x , T ) x = n = 0 n π C n ( T ) cos ( n π x ) , 2 u 1 ( x , T ) x 2 = n = 0 n π 2 C n ( T ) sin ( n π x ) , 3 u 1 ( x , T ) x 3 = n = 0 n π 3 C n ( T ) cos ( n π x ) , 4 u 1 ( x , t ) x 4 = n = 0 n π 4 C n ( T ) sin ( n π x ) ,
are right, then, form (15) and (16) by means of the Weierstrass criterion which leads to the convergence of the series, therefore
u 1 ( x , T ) , u 1 ( x , T ) x , 2 u 1 ( x , T ) x 2 , 3 u 1 ( x , T ) x 3 C [ 0 , 1 ] .
We will show that 4 u 1 ( x , T ) x 4 L 2 [ 0 , 1 ] . From (16) and (17), we obtain
n π 4 C n ( T ) = 2 n π 0 T e ( n π ) 2 ( T τ ) q 1 ( τ ) d τ + 2 n π q 1 ( T ) 2 n π q 1 ( T ) .
From (12), it follows that
2 n π 0 T e ( n π ) 2 ( T τ ) q 1 ( τ ) d τ q 1 ( t ) L 2 [ 0 , T ] 2 ( n π ) 2 ,
First series
n = 0 2 n π 0 T e ( n π ) 2 ( T τ ) q 1 ( τ ) d τ .
Second series
2 q 1 ( T ) n = 0 n π < .
Third series
2 q 1 ( T ) n = 0 1 n π ,  
absolutely converges on 0 , 1 × 0 , T then 4 u 1 ( x , T ) x 4 L 2 [ 0 , 1 ] . □

3. Expansion of the Direct Problem (1)–(5) on T ,

Let us study the following PDE system in the interval T , .
u 2 x , t t = 2 u 2 x , t x 2 ,   0 x 1 ,   t T , ,
u 2 x , T = u 1 x , T = f ( x ) ,   0 x 1 ,
u 2 0 , t = 0 ,   T t < ,
u 2 x 1 , t + κ u 2 1 , t = 0 ,   T t < , κ > 0 .
Assume that
f ( x ) W 2 2 0 , 1 ; f ( 0 ) = 0 ,   f ( 1 ) + κ f ( 1 ) = 0 ,
We obtain the following solution by applying the separation of variables as a way for solving problem (18)–(21)
u 2 ( x , t ) = n = 1 A n e λ n 2 ( t T ) sin λ n x , 0 x 1 , t T ,
where cot λ n = λ n κ ,
A n = 4 2 λ n sin 2 λ n 0 1 f ( x ) sin λ n x d x .
By integrating the right side of (24) twice, we obtain
A n = 4 λ n 2 λ n sin 2 λ n 0 1 f ( x ) sin λ n x d x ,
From (22) and (25), we define a number c 1 > 0 such that for any n
A n c 1 λ n 2 .
From (23) and (26), any t T + 1
u 2 ( x , t ) c 1 n = 1 λ n 2 e λ n 2 ( t T ) ,
u 2 ( x , t ) x c 1 n = 1 λ n 1 e λ n 2 ( t T ) ,
then,
2 u 2 ( x , t ) x 2 c 1 n = 1 e λ n 2 ( t T ) .
Since
e λ n 2 ( t 1 ) = e λ n 2 e λ n 2 ( t T 1 ) ,
Let us consider there exists the numbers c 2 and c 3 such that for any n
c 2 ( n + 1 ) λ n c 3 ( n + 1 ) ,
and, it follows from (30) and (31), that
e λ n 2 e c 1 2 n ,
then, it follows from (18), (27)–(32) that there is c 4 > 0 known as a number such that for any t T + 2
sup 0 x 1 u 2 ( x , t ) , u 2 ( x , t ) x , u 2 ( x , t ) t , 2 u 2 ( x , t ) x 2 c 4 e ( t T 1 ) .
Now, let us examine the behavior u 2 ( 1 , t ) t ,
Lemma 3.
Let A n be defined by the formula (24). Then
A n = 4 2 λ n sin 2 λ n 1 λ n 4 0 1 f ( 4 ) ( x ) sin λ n x d x f ( 1 ) sin λ n + 1 λ n 3 f ( 1 ) cos λ n ,
where f ( 4 ) ( x ) is the fourth derivative with respect to x for function f ( x ) .
Proof. 
A n defined by the Equation (24), and integrating 0 1 f ( x ) sin λ n x d x in parts twice, we obtain
0 1 f ( x ) sin λ n x d x = sin λ n λ n 2 f ( 1 ) + cos λ n λ n f ( 1 ) 1 λ n f ( 0 ) 1 λ n 2 0 1 f ( x ) sin λ n x d x ,
from (3) and (19)
f ( 0 ) = u 1 ( 0 , T ) = 0 .
Since
sin λ n λ n 2 f ( 1 ) + cos λ n λ n f ( 1 ) = sin λ n λ n 2 u 2 x 1 , t + cos λ n λ n u 2 1 , t = cos λ n λ n 2 u 2 x 1 , t + λ n u 2 1 , t = λ n κ cos λ n κ λ n u 2 x 1 , t + κ u 2 1 , t = 0 ,
as a result, we obtain
0 1 f ( x ) sin λ n x d x = 1 λ n 2 0 1 f ( x ) sin λ n x d x ,
Integrating the right part of the previous equation twice in parts, it leads to
1 λ n 2 0 1 f ( x ) sin λ n x d x = 1 λ n 3 cos λ n f ( 1 ) f ( 0 ) 1 λ n 4 sin λ n f ( 1 ) 0 1 f ( 4 ) ( x ) sin λ n x d x ,
The lemma is proofed. □
From Lemmas 2 and 3, the series n = 0 λ n 2 A n < ; hence, from (23), we obtain
u 2 ( 1 , t ) t = n = 0 λ n 2 A n e λ n 2 ( t T ) sin λ n x
u 2 ( 1 , t ) t C T , , and   there   exists ,   lim t T u 2 ( 1 , t ) t < .
Denote u 2 ( 1 , t ) = q 2 ( t ) , from (34) and (35), it follows that q 2 ( t ) H 1 T , .
Lemma 4.
Let the function u 2 ( 1 , t ) t be defined by Equation (34). Then, d 1 > 0 such that for any t T + 1
u 2 ( 1 , t ) t d 1 e ( t T 1 ) + A 0 e λ 0 2 ( t T ) .
Proof. 
From (34) and (35), it follows that
u 2 ( 1 , t ) t d 1 n = 0 λ n 2 e λ n 2 ( t T ) ,
where d 1 some number.
Let us assume that t T + 1 and n > 0
λ n 2 e λ n 2 e λ n 2 ( t T 1 ) λ n 2 e λ n 2 e ( t T 1 ) .
From d 1 n λ n d 2 n it follows that, for n > 0 and numbers d 1 , d 2 > 0
λ n 2 e λ n 2 d 2 n 2 e d 1 2 n ,
from (36), it follows that n = 0 n 2 e d 1 2 n < . Hence there is a number d 3 for any t T + 2
u 2 ( 1 , t ) t d 3 e ( t T 1 ) + A 0 e λ 0 2 ( t T ) ,
from (35) and Lemma 4, it follows that
u 2 ( 1 , t ) t C T , L 1 T , L 2 T , .
Now, let us introduce the notation
h ( t ) = q 1 ( t ) ,   t 0 , T q 2 ( t ) ,   t T , ,   u ( x , t ) = u 1 ( x , t ) ,   t 0 , T , u 2 ( x , t ) ,   t T , ,   x 0 , 1 ,
From (33) and (37), it follows that, for any ε > 0 there is χ ε ( t ) which is defined as a function such that, for any t 0
sup 0 x 1 u ( x , t ) , u ( x , t ) x , 2 u ( x , t ) x 2 χ ε ( t ) ,
where
χ ε ( t ) = d 6 ( ε ) , 0 t T + 2 , d 5 ( ε ) n = 1 λ n 2 e λ n 2 ,   t T + 2 .
Since χ ε ( t ) L 1 0 , , then the Fourier transform for t can be used for the combined direct problem (1)–(5) and (18)–(21).
The lemma is proofed. □
From Lemma 1 and Equation (38), we obtain the following theorem.
Theorem 1.
Let Φ ( t ) C 0 , and Φ ( t ) is limited over this line. Then, the following relations are true
0 u ( x , t ) x Φ ( t ) d t = x 0 u ( x , t ) Φ ( t ) d t , 0 2 u ( x , t ) x 2 Φ ( t ) d t = 2 x 2 0 u ( x , t ) Φ ( t ) d t
Lemma 5.
Let u ( x , t ) be a solution of the combined problem (1)–(5) and (18)–(21). Then, the following relations are true
lim x 0 0 u ( x , t ) h ( t ) d t = lim x 1 0 u ( x , t ) u ( 1 , t ) d t = lim x 1 0 u x ( x , t ) u x ( 1 , t ) d t = 0
.
Proof. 
It follows from Lemma 1 and (35) that, for any t > 0
lim x 0 u ( x , t ) = h ( t ) ,   lim x 1 u ( x , t ) = u ( 1 , t )   and   lim x 1 u x ( x , t ) = u x ( 1 , t ) .
Let the number d 8 be defined by the formula
d 8 = max u ( x , t ) + u x ( x , t ) : 0 t T + 2 .
Then, let us denote by s ( t ) the function defined by the formula
s ( t ) = d 8 ; 0 t T + 2 , c 4 e ( t T 1 ) ;   t > T + 2 .
Since 0 s ( t ) d t < and for any t > 0
u ( x , t ) s ( t ) , u x ( x , t ) s ( t ) ,
then, given (39); by the Lebesgue theorem on the passage to the limit under the integral sign, the assertion of the lemma is proved. □

4. Solution of the Inverse BVPs (1)–(5) and (18)–(21)

Let us assume that the function h ( t ) in the combined problem (1)–(5) and (18)–(21) is unknown, and, instead, the function is given as g ( t ) = u ( x 0 , t ) , where x 0 0 , 1 , t 0 .
Let us adopt that, for g ( t ) = g 0 ( t ) , there is a function h 0 ( t ) H 1 0 , such that, when it is substituted into the boundary of (1)–(5) and (18)–(21), we obtain a real solution u 0 ( x , t ) which is defined as the following
u 0 ( x 0 , t ) = g 0 ( t ) .
Function g 0 ( t ) unknown, and, instead, we have g δ ( t ) and δ > 0 such that
g δ ( t ) g 0 ( t ) L 2 0 , δ .
It is necessary to use the given data g δ ( t ) and δ inverse BVP (1)–(5) and (18)–(21) in order to find an approximate solution h δ ( t ) and obtain an error estimate h δ ( t ) h 0 ( t ) L 2 0 , .

5. Solution of the Inverse BVP (1)–(5) and (18)–(21) by the Projection Regularization Method

Let H ¯ = L 2 ; + + i L 2 ; + be the interval on the area of complex numbers, and the set of correction class M r H ¯ demarcated by the following
M r = h ( t ) : h ( t ) H ¯ ,   0 h ( t ) 2 d t + h ( t ) 2 d t r 2 ,
r known positive number.
In order to resolve the problem (1)–(5) and (18)–(21), we present F , as the operator which is mapping from H ¯ to H ¯ and we named as the operator via the Fourier transform
h ^ ( τ ) = F h ( t ) = 1 2 π 0 + h ( t ) e i τ t d t ,   τ ,   h ( t ) H ¯ L ¯ 1 , .
There L ¯ 1 , —interval on the of complex numbers set.
Denote by F ¯ operator continuation F in H ¯ . Following from Plancherel’s theorem, the operator F ¯ has isometric mapping H ¯ into H ¯ .
Let h ^ ( τ ) H ¯ . Then, we have
F ¯ 1 h ^ ( τ ) = 1 2 π lim N N + N h ^ ( τ ) e i t τ d τ ,   < t < + ,
where the way to the limit has the sense of the convergence of root-mean-square.
Using transform F , (1)–(5) and (18)–(21) come down to the following problem
2 u ^ ( x , τ ) x 2 = i τ u ^ ( x , τ ) ,   x ( 0 , 1 ) ,   < τ < + ,
u ^ ( 0 , τ ) = 0 ,   u ^ ( x 0 , τ ) = g ^ ( τ ) , < τ < + ,
where u ^ ( x , τ ) = F u ( x , t ) ,   g ^ ( τ ) = F g ( t ) .
Solutions (45) and (46) are of the form
u ^ ( x , τ ) = D 1 ( τ ) e μ 0 τ x + D 2 ( τ ) e μ 0 τ x ,   τ > 0 ,
where μ 0 = 1 2 ( 1 + i ) ,   D 1 ( τ ) and D 2 ( τ ) are functions that satisfy (40) and (46).
With u ^ ( 1 , τ ) = h ^ ( τ ) we obtain
h ^ ( τ ) = sh   μ 0 τ sh   μ 0 τ x 0 g ^ ( τ ) ,   τ > 0 , 0 , τ 0 ,
Therefore, the problem (45) and (46) reduces to the equation
A h ^ ( τ ) = g ^ ( τ ) ,   0 < τ < .
Let g ^ 0 ( τ ) = F g 0 ( t ) , g ^ δ ( τ ) = F g δ ( t ) and, from the Formula (41), it follows that
g ^ δ ( τ ) g ^ 0 ( τ ) H ¯ 2 δ .
Let M ^ r denote a set of H ¯ such that M ^ r F M r and
M ^ r = h ^ ( τ ) : h ^ ( τ ) H ¯ ,   0 1 + τ 2 h ^ ( τ ) d τ 2 r 2 .
Since h 0 ( t ) M r , then h ^ 0 ( τ ) M ^ r .
In order to find the approximation solution for (49)–(51) we use the regularizing family of operators R α , which are defined by
h ^ δ α ( τ ) = R α g ^ ( τ ) = sh   μ 0 τ sh   μ 0 τ x 0 g ^ δ ( τ ) ,   0 τ α , 0 , τ > α ,
For selecting a regularization parameter α ^ = α ^ ( g ^ δ , δ ) in Equation (52) from the initial data ( g ^ δ , δ ) , use the equation A h ^ δ α ( τ ) g ^ δ ( τ ) 2 = 16 δ 2 .
Let us describe an estimated solution for (49) by the formulation of h ^ δ ( τ ) = h ^ δ α ^ ( g δ , δ ) ( τ ) .
This follows from the theorem formulated in the article [25] [c. 284], that
h ^ δ ( τ ) h ^ 0 ( τ ) 7 ω ( δ , r ) ,
where ω ( δ , r ) = h ^ ( τ ) : h ^ ( τ ) M ^ r , A h ^ ( τ ) δ .
Let us describe R α ^ ( g ^ δ , δ ) : 0 < δ δ 0 as the operator for use in the regularization method in order to obtain the approximate solution for the problem. (49) in M ^ r . Now, let us introduce δ 0 , δ 0 as the quantitative characteristic of the accuracy of this method on the set M ^ r .
Δ δ R α ^ ( f ^ δ , δ ) = sup f ^ 0 , f ^ δ R α ^ ( g ^ δ , δ ) g ^ δ ( τ ) h ^ 0 ( τ ) : h ^ 0 ( τ ) M ^ r , g ^ δ ( τ ) H , g ^ δ ( τ ) A h ^ 0 ( τ ) δ .
From the theorem proved in [23], it follows that the following estimate holds
Δ δ R α ^ ( g ^ δ , δ ) ω ( δ , r ) .
Let
Φ 2 ( α ) = sup α h ^ 0 ( τ ) 2 d τ : h ^ 0 ( τ ) M ^ r .
From (51) and (55), we obtain Φ 2 ( α ) = 2 r 2 1 + α 2 for h ^ 0 ( τ ) M ^ r .
Lemma 6.
Let α 0 = 1 2 x 0 2 ln 2 2 . Then, for α α 0 the ratio is true
1 4 e ( 1 x 0 ) α 2 R α 4 e ( 1 x 0 ) α 2 .
Lemma 6 tracks from the explanation of the operator norm. According to [26], lemma 2, to compute the modulus of continuity, ω ( δ , r ) we need to solve
r α G ( α ) = δ .
Solving α ¯ ( δ ) is replaced into the function G ( α ) parameter determined by
G ˜ ( β ) = 1 1 + β 2 ,   α = e x 0 1 β 2 .
From (56) and (57), it follows that
ω ( δ , r ) = r G ( α ¯ ( δ ) ) .
Therefore, from (53), (57) and (58), we obtain the estimate
h ^ δ ( τ ) h ^ 0 ( τ ) 7 r G ( α ¯ ( δ ) ) .
In order to simplify the assessment (59), consider the equations
e ( x 0 1 ) α 2 = r δ ,   e 2 ( x 0 1 ) α 2 = r δ .
Let α ¯ 1 ( δ ) and α ¯ 2 ( δ ) , respectively, be solutions of the Equation (60).
Then, from (56), (60), we find that, for sufficiently small δ , defined α ¯ 2 ( δ ) , the following relations are valid
α ¯ 2 ( δ ) α ¯ ( δ ) α ¯ 1 ( δ ) ,
where α ¯ 1 ( δ ) = 2 x 0 1 2 ln 2 r δ , α ¯ 2 ( δ ) = 2 2 x 0 1 2 ln 2 r δ and, from the resulting inequality, we have
α ¯ ( δ ) ln 2 δ   a t   δ 0 .
From the theorem proved in [26], it follows that
G ( α ¯ 2 ( δ ) ) G ( α ¯ ( δ ) ) G ( α ¯ 1 ( δ ) ) ,
where
G ( α ¯ 1 ( δ ) ) = 4 1 + 4 x 0 1 4 ln 4 δ ,   G ( α ¯ 2 ( δ ) ) = 1 4 1 + 1 x 0 1 4 ln 4 δ ,
from (54) we find that this is an exact ordinal estimate,
sup h ^ δ α ¯ ( δ ) ( σ ) h ^ 0 ( σ ) : h ^ 0 ( σ ) M ^ r , g ^ δ ( σ ) g ^ 0 ( σ ) δ r 4 1 + 1 4 ( x 0 1 ) 4 ln 4 r δ .
From lemma 5, (53) and (63) we obtain
Theorem 2.
For method R α ^ ( g ^ δ , δ ) : 0 < δ δ 0 we have an exact estimate of the order error
r 4 1 + 1 4 ( x 0 1 ) 4 ln 4 δ Δ δ R α ^ ( f δ , δ ) 74 r 1 + 1 4 ( x 0 1 ) 4 ln 4 δ .
Applying к q ^ δ ( τ ) transformation
h δ ( t ) = Re F ¯ 1 h ^ δ ( τ ) ,   t 0 , 0 , t < 0 ,
where F ¯ 1 is the inverse Fourier transform operator, we obtain an estimated solution for the problem (1)–(5) and (18)–(21).
Thus, for an approximate solution h δ ( t ) for problem (1)–(5) and (18)–(21), we have a precise error estimation by
h δ ( t ) h 0 ( t ) 28 r 1 + 4 x 0 1 4 ln 4 δ .

6. Case Study

Consider the function h 0 ( t ) = t sin ( π t ) ,   t 0 , T , 0 , t T , , suppose T = 1 , x 0 = 0.5   and   N = 100 .
From the solution of the direct problem (1)–(5) and (18)–(21), we find u 0 ( x 0 , t ) = g 0 ( t ) . We set a partition of the time interval 0 , T with the number of nodes N such that
0 = t 0 < t 1 < t 2 < < t n = T ,   t k = k N ,   k = 0 , N ¯ .
This simulates the one-dimensional nature of the heat equation using the Fast Fourier Transform, FFT, as shown in Figure 1. In this example, the PDE system is linear, and it is possible to advance the system directly in the frequency domain.
From Figure 1 we find g δ ( t ) , introducing an error level δ 1 = 0.05 and δ 2 = 0.02 in g 0 ( t ) by the following
g δ k = g 0 ( t k ) + δ j T . γ , j = 1 , 2 ,
where the error level can compute by
g δ g 0 .
Figure 2 and Figure 3 show the visualization of the function as a solution for the inverse problem with δ = 0.02 and δ = 0.05 , respectively. The real solution is shown by a dotted line and the approximate solution is shown by a line.

7. Conclusions

In this work, the inverse BVP of a thermal conductivity equation in two different intervals of time was solved. The heating process for an object was definitively separated into two intervals: the first one, by the heating of the boundary part or place in domain, and the second one, by the free cooling of the object. It has been verified that the boundary function or condition fits the space H 1 0 , . This means we can use the projection regularization method for solving this problem by using the Fourier transform for time. The error estimate was obtained for the solution.

Author Contributions

Conceptualization, M.A., H.K.A.-M. and B.T.A.-N.; software, M.A.; H.K.A.-M. and H.A.; validation, Z.A. and H.A.; formal analysis, E.-S.M.E.-k.; investigation, M.A. and E.-S.M.E.-k.; writing—original draft preparation, H.K.A.-M., M.A. and B.T.A.-N.; writing—review and editing, M.A. and H.A.; visualization, H.A.; supervision, M.A. and B.T.A.-N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tikhonov, A.N. On the Regularization of Ill-Posed Problems. Proc. USSR Acad. Sci. 1963, 153, 49–52. [Google Scholar]
  2. Lavrent’ev, M.M. On Certain Ill-Posed Problems of Mathematical Physics; Sobolev Division of the USSR Academy Science: Novosibirsk, Russia, 1962. [Google Scholar]
  3. Ivanov, V.K.; Vasin, V.V.; Tanana, V.P. Theory of Linear Ill-Posed Problem and Application; Nauok: Moscow, Russia, 1978. [Google Scholar]
  4. John, F. Numerical solution of the equation of heat conduction for preceding times. Ann. di Mat. pura ed Appl. 1955, 40, 129–142. [Google Scholar] [CrossRef]
  5. Cheniguel, A. Numerical method for the heat equation with Dirichlet and Neumann conditions. Proc. Int. MultiConference Eng. Comput. Sci. 2014, 1, 12–14. [Google Scholar]
  6. Alifanov, O.M.; Artioukhine, E.A.; Rumyantsev, S.V. Extreme Methods for Solving Ill-Posed Problems with Applications to Inverse Heat Transfer Problems; Begell House: New York, NY, USA, 1995; ISBN 156700038X. [Google Scholar]
  7. Bergman, T.L.; Incropera, F.P.; DeWitt, D.P.; Lavine, A.S. Fundamentals of Heat and Mass Transfer; John Wiley & Sons: Hoboken, NJ, USA, 2011; ISBN 0470501979. [Google Scholar]
  8. Dmitriev, V.I.; Stolyarov, L.V. Numerical method for the inverse boundary-value problem of the heat equation. Comput. Math. Model. 2017, 28, 141–147. [Google Scholar] [CrossRef]
  9. Sidikova, A.I. Approximate Solution of Inverse Boundary Value Problem for Heat Conduction Equation. In Proceedings of the 2018 International Russian Automation Conference (RusAutoCon), Sochi, Russia, 9–16 September 2018; pp. 1–7. [Google Scholar]
  10. Sidikova, A.I. A Study of an Inverse Boundary Value Problem for the Heat Conduction Equation. Numer. Anal. Appl. 2019, 12, 70–86. [Google Scholar] [CrossRef]
  11. Tanana, V.P.; Markov, B.A. The control problem for the heat equation in the case of a composite material. J. Phys. Conf. Ser. 2021, 1715, 12049. [Google Scholar] [CrossRef]
  12. Makhtoumi, M. Numerical solutions of heat diffusion equation over one dimensional rod region. arXiv 2018, arXiv:1807.09588. [Google Scholar]
  13. Sidikova, A.I.; Al-Mahdawi, H.K. The solution of inverse boundary problem for the heat exchange for the hollow cylinder. AIP Conf. Proc. 2022, 2398, 60050. [Google Scholar]
  14. Al-Mahdawi, H.K.I.; Abotaleb, M.; Alkattan, H.; Tareq, A.-M.Z.; Badr, A.; Kadi, A. Multigrid Method for Solving Inverse Problems for Heat Equation. Mathematics 2022, 10, 2802. [Google Scholar] [CrossRef]
  15. Al-Mahdawi, H.K.; Sidikova, A.I. Iterated Lavrent’ev regularization with the finite-dimensional approximation for inverse problem. AIP Conf. Proc. 2022, 2398, 60080. [Google Scholar]
  16. Al-Mahdawi, H.K.; Sidikova, A.I.; Alkattan, H.; Abotaleb, M.; Kadi, A.; El-kenawy, E.-S.M. Parallel Multigrid Method for Solving Inverse Problems. MethodsX 2022, 9, 101887. [Google Scholar] [CrossRef] [PubMed]
  17. Al-Mahdawi, H.K.I.; Alkattan, H.; Abotaleb, M.; Kadi, A.; El-kenawy, E.-S.M. Updating the Landweber Iteration Method for Solving Inverse Problems. Mathematics 2022, 10, 2798. [Google Scholar] [CrossRef]
  18. Glasko, V.B.; Kulik, N.I.; Shklyarov, I.N.; Tikhonov, A.N. An inverse problem of heat conductivity. Zhurnal Vychislitel’noi Mat. I Mat. Fiz. 1979, 19, 768–774. [Google Scholar]
  19. Belonosov, A.S.; Shishlenin, M.A. Continuation problem for the parabolic equation with the data on the part of the boundary. Siber. Electron. Math. Rep. 2014, 11, 22–34. [Google Scholar]
  20. Kabanikhin, S.I.; Hasanov, A.; Penenko, A.V. A gradient descent method for solving an inverse coefficient heat conduction problem. Numer. Anal. Appl. 2008, 1, 34–45. [Google Scholar] [CrossRef]
  21. Yagola, A.G.; Stepanova, I.E.; Van, Y.; Titarenko, V.N. Obratnye zadachi i metody ikh resheniya. Prilozheniya k geofizike. In Inverse Problems and Methods for their Solution: Applications to Geophysics; BKL Publishers: Moscow, Russia, 2014. [Google Scholar]
  22. Kabanikhin, S.I.; Krivorot’ko, O.I.; Shishlenin, M.A. A numerical method for solving an inverse thermoacoustic problem. Numer. Anal. Appl. 2013, 6, 34–39. [Google Scholar] [CrossRef]
  23. Tanana, V.P. On the order-optimality of the projection regularization method in solving inverse problems. Sib. Zhurnal Ind. Mat. 2004, 7, 117–132. [Google Scholar]
  24. Farlow, S.J. Partial Differential Equations for Scientists and Engineers; Courier Corporation: North Chelmsford, MA, USA, 1993; ISBN 048667620X. [Google Scholar]
  25. Tanana, V.P.; Bredikhina, A.B.; Kamaltdinova, T.S. On an error estimate for an approximate solution for an inverse problem in the class of piecewise smooth functions. Tr. Inst. Mat. I Mekhaniki UrO RAN 2012, 18, 281–288. [Google Scholar]
  26. Tanana, V.P.; Rudakova, T.N. The optimum of the M. M. Lavrent’ev method. J. Inverse Ill-Posed Probl. 2011, 18, 935–944. [Google Scholar] [CrossRef]
Figure 1. Evolution of the direct heat Equations (45) and (46).
Figure 1. Evolution of the direct heat Equations (45) and (46).
Algorithms 16 00033 g001
Figure 2. Visualization of the answer with δ = 0.02 .
Figure 2. Visualization of the answer with δ = 0.02 .
Algorithms 16 00033 g002
Figure 3. Visualization of the answer with δ = 0.05 .
Figure 3. Visualization of the answer with δ = 0.05 .
Algorithms 16 00033 g003
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

Al-Nuaimi, B.T.; Al-Mahdawi, H.K.; Albadran, Z.; Alkattan, H.; Abotaleb, M.; El-kenawy, E.-S.M. Solving of the Inverse Boundary Value Problem for the Heat Conduction Equation in Two Intervals of Time. Algorithms 2023, 16, 33. https://doi.org/10.3390/a16010033

AMA Style

Al-Nuaimi BT, Al-Mahdawi HK, Albadran Z, Alkattan H, Abotaleb M, El-kenawy E-SM. Solving of the Inverse Boundary Value Problem for the Heat Conduction Equation in Two Intervals of Time. Algorithms. 2023; 16(1):33. https://doi.org/10.3390/a16010033

Chicago/Turabian Style

Al-Nuaimi, Bashar Talib, H.K. Al-Mahdawi, Zainalabideen Albadran, Hussein Alkattan, Mostafa Abotaleb, and El-Sayed M. El-kenawy. 2023. "Solving of the Inverse Boundary Value Problem for the Heat Conduction Equation in Two Intervals of Time" Algorithms 16, no. 1: 33. https://doi.org/10.3390/a16010033

APA Style

Al-Nuaimi, B. T., Al-Mahdawi, H. K., Albadran, Z., Alkattan, H., Abotaleb, M., & El-kenawy, E. -S. M. (2023). Solving of the Inverse Boundary Value Problem for the Heat Conduction Equation in Two Intervals of Time. Algorithms, 16(1), 33. https://doi.org/10.3390/a16010033

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