Home Simultaneously identifying the thermal conductivity and radiative coefficient in heat equation from Dirichlet and Neumann boundary measured outputs
Article
Licensed
Unlicensed Requires Authentication

Simultaneously identifying the thermal conductivity and radiative coefficient in heat equation from Dirichlet and Neumann boundary measured outputs

  • Alemdar Hasanov EMAIL logo
Published/Copyright: September 5, 2020

Abstract

This paper deals with an inverse coefficient problem of simultaneously identifying the thermal conductivity k(x) and radiative coefficient q(x) in the 1D heat equation ut=(k(x)ux)x-q(x)u from the most available Dirichlet and Neumann boundary measured outputs. The Neumann-to-Dirichlet and Neumann-to-Neumann operators Φ[k,q](t):=u(,t;k,q), Ψ[k,q](t):=-k(0)ux(0,t;k,q) are introduced, and main properties of these operators are derived. Then the Tikhonov functional

J(k,q)=12Φ[k,q]-νL2(0,T)2+12Ψ[k,q]-φL2(0,T)2

of two functions k(x) and q(x) is introduced, and an existence of a quasi-solution of the inverse coefficient problem is proved. An explicit formula for the Fréchet gradient of the Tikhonov functional is derived through the weak solutions of two appropriate adjoint problems.

MSC 2010: 65N21; 80A23

Funding statement: This research has been supported by The Scientific and Technological Research Council (TUBITAK) of Turkey through the Incentive Program for International Scientific Publications (UBYT).

Appendix A Some necessary estimates for the solutions of the direct and adjoint problems

We derive here some necessary estimates for the weak and regular weak solutions of direct problem (1.1).

A.1 Estimates for the weak and regular weak solution of direct problem (1.1)

It is known that, under conditions (2.1), there exists a unique weak solution of problem (1.1) defined as uL2(0,T;V(0,)), utL2(0,T;L2(0,)), uttL2(0,T;H-1(0,)), where V(0,):={vH1(0,):v(0)=0} (see [3]).

Theorem A.1

Let conditions (2.1) hold. Then, for the weak solution of direct problem (1.1), the following estimates hold:

(A.1)uL2(0,T;L2(0,))C1fL2(0,),
(A.2)uxL2(0,T;L2(0,))C2fL2(0,),
where C12=[exp(C0T)-1]k0C0, C22=exp(C0T)k0, C0=2k0.

Proof

Multiply both sides of (1.1) by u(x,t), integrate on Ωt:=(0,)×(0,t) and use the integration by parts formula. Taking into account the initial and boundary conditions, we obtain the energy identity

(A.3)0u2(x,t)dx+2Ωtk(x)ux2(x,τ)dxdτ+2Ωtq(x)u2(x,τ)dxdτ=20tf(τ)u(,τ)dτ

for a.e. t[0,T]. Apply the 𝜀-inequality to the last integral, and use the inequality

0tu2(,τ)dτ2Ωtu2(x,τ)dxdτ+Ωtux2(x,τ)dxdτ,t[0,T],

to get

20tf(τ)u(,τ)2εΩtu2(x,τ)dxdτ+εΩtux2(x,τ)dxdτ+1ε0tf2(τ)dτ.

Use this inequality in (A.3), and choose ε>0 as ε=k0, where k0>0 is the constant defined in (2.1). After elementary transformations, we obtain the integral inequality

(A.4)0u2dx+k0Ωtux2dxdτ+2Ωtq(x)u2dxdτ2k02Ωtu2dxdτ+k00Tf2(t)dt

for a.e. t[0,T]. Estimates (A.1) and (A.2) follow from this inequality. ∎

The following results are proved in the same way as Theorem A.1

Theorem A.2

Let conditions (2.1) hold. Assume in addition that fH1(0,T). Then, for the regular weak solution uL2(0,T;H2(0,)), with utL2(0,T;V(0,)), uttL2(0,T;L2(0,)), of direct problem (1.1), the following estimates hold:

(A.5)utL2(0,T;L2(0,))C1fL2(0,T),
(A.6)uxtL2(0,T;L2(0,))C2fL2(0,T),
where C1,C2>0 are the same constants defined in Theorem A.1.

Theorem A.3

Let conditions (2.1) hold. Assume in addition that fH2(0,T). Then, for the regular weak solution with improved regularity uL2(0,T;H4(0,)), with utL2(0,T;H2(0,)), uttL2(0,T;V(0,)), utttL2(0,T;L2(0,)), of direct problem (1.1), the following estimates hold with the same constants C1,C2>0 defined in Theorem A.1:

(A.7)uttL2(0,T;L2(0,))C1f′′L2(0,T),
uxttL2(0,T;L2(0,))C2f′′L2(0,T).

A.2 Estimates for the weak solution of adjoint problem (4.5)

In this subsection, we derive necessary estimates for the weak solution of adjoint problem (4.5) with the input p(t) defined in (4.8), i.e. for the solution of the problem

(A.8){ϕt=-(k(x)ϕx)x+q(x)ϕ,(x,t)(0,)×[0,T),ϕ(x,T)=0,x(0,l),ϕ(0,t)=0,k()ϕx(,t)=u(,t)-ν(t),t(0,T),

Evidently, estimates (A.1) and (A.2) can also be used for the solution of problem (A.8) with f(t) replaced by p(t)=u(,t)-ν(t). In view of the inequalities

pL2(0,T)u(,)L2(0,T)+νL2(0,T)andu(,)L2(0,T)uxL2(0,T;L2(0,)),

from estimates (A.1) and (A.2), we deduce that

(A.9)ϕL2(0,T;L2(0,))C1[C2fL2(0,T)+νL2(0,T)],
(A.10)ϕxL2(0,T;L2(0,))C2[C2fL2(0,T)+νL2(0,T)].

Theorem A.4

Let conditions (2.1) hold. Assume that νL2(0,T). Then, for the weak solution of adjoint problem (A.8), estimates (A.9) and (A.10) hold.

A.3 Estimates for the weak solution of adjoint problem (4.8)

Here we derive necessary estimates for the weak solution of the second adjoint problem (4.7) with the Dirichlet input μ(t) defined in (4.10), i.e. for the solution of the problem

(A.11){ψt=-(k(x)ψx)x+q(x)ψ,(x,t)(0,)×[0,T),ψ(x,T)=0,x(0,l),ψ(0,t)=k(0)ux(0,t)+φ(t),ψx(,t)=0,t(0,T),
Theorem A.5

Let conditions (2.1) hold. Assume that fH2(0,T) and, in addition, the Neumann measured output satisfies the condition φH1(0,T). Then, for the weak solution of adjoint problem (A.11), the following estimates hold:

(A.12)ψL2(0,T;L2(0,))C^3(2C3fH2(0,)2+φH1(0,)2),
(A.13)ψxL2(0,T;L2(0,))C^2(2C3fH2(0,)2+φH1(0,)2),
where

(A.14){C^12=2CT2[exp(T)-1],C^22=2k0exp(T),CT2=max{2T;1+T},C^32=2max{C^12;},

and C3>0 is the constant defined in (2.4).

Proof

We use the transformation

(A.15)v(x,t)=ψ(x,t)-μ(t),(x,t)ΩT¯,

where μ(t)=k(0)ux(0,t)+φ(t). Then the function v(x,t) solves the problem

(A.16){vt=-(k(x)vx)x+q(x)v-μ(t),(x,t)(0,)×[0,T),v(x,T)=-μ(T),x(0,l),v(0,t)=0,vx(,t)=0,t(0,T).

In the same way as in Theorem A.3, we can show that the energy identity corresponding to the transformed problem (A.16) is as follows:

0v2(x,t)dx+2tT0k(x)vx2(x,τ)dxdτ+2tT0q(x)u2(x,τ)dxdτ=2tT0μ(τ)v(x,τ)dxdτ+μ2(T),t[0,T].

We employ the inequality 2aba2+b2 for the first right-hand side integral and then use the inequality

[0T(μ(t))2dt+μ2(T)]CT2μH1(0,T)2,

with the constant CT>0 defined in (A.14) to deduce the main integral inequality

0v2dx+2k0tT0vx2dxdτ+2tT0q(x)u2dxdτtT0v2dxdτ+CT2μH1(0,T)2,t[0,T].

This inequality leads to the estimates

vL2(0,T;L2(0,))C^1μH1(0,T),vxL2(0,T;L2(0,))C^2μH1(0,T),

with the constants C^1,C^2>0 defined in (A.14).

In view of transformation (A.15), these estimates imply

ψL2(0,T;L2(0,))C^3μH1(0,T),ψxL2(0,T;L2(0,))C^2μH1(0,T),

where C^3>0 is defined in (A.14). Taking into account that the input μ(t) in the above estimates is defined as μ(0,t)=k(0)ux(0,t)+φ(t) and using estimates (2.4) and (2.5), we arrive at the required estimates (A.12) and (A.13). ∎

Acknowledgements

The author would like to thank Dinh Nho Hào for useful discussion.

References

[1] J. R. Cannon and P. DuChateau, An inverse problem for a nonlinear diffusion equation, SIAM J. Appl. Math. 39 (1980), 272–289. 10.1137/0139024Search in Google Scholar

[2] M. Choulli and M. Yamamoto, Generic well-posedness of an inverse parabolic problem – the Hölder space approach, Inverse Problems 12 (1996), no. 3, 195–205. 10.1088/0266-5611/12/3/002Search in Google Scholar

[3] L. C. Evans, Partial Differential Equations, AMC, Rhode Island, 2002. Search in Google Scholar

[4] D. N. Hào, Methods for Inverse Heat Conduction Problems, Peter Lang, Frankfurt/Main, 1998. 10.1080/174159798088027675Search in Google Scholar

[5] A. Hasanov, P. DuChateau and B. Pektas, An adjoint problem approach and coarse-fine mesh method for identification of the diffusion coefficient in a linear parabolic equation, J. Inverse Ill-Posed Probl. 14 (2006), 435–463. 10.1515/156939406778247615Search in Google Scholar

[6] A. Hasanov and Romanov, Introduction to Inverse Problems for Differential Equations, Springer, Heidelberg, 2017. 10.1007/978-3-319-62797-7Search in Google Scholar

[7] S. I. Kabanikhin, A. Hasanov and A. V. Penenko, The gradient-based method for solving the inverse coefficient heat-conduction problem, Siberian J. Num. Math. 11 (2008), 41–54. Search in Google Scholar

[8] R. Kohn and M. Vogelius, Determining conductivity by boundary measurements, Comm. Pure Appl. Math. 37 (1984), no. 3, 289–298. 10.1002/cpa.3160370302Search in Google Scholar

[9] E. Zeidler, Applied Functional Analysis, Main Principles and Their Applications, Springer, New York, 1995. Search in Google Scholar

Received: 2020-04-22
Accepted: 2020-08-23
Published Online: 2020-09-05
Published in Print: 2021-02-01

© 2020 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 28.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/jiip-2020-0047/html
Scroll to top button