Home Mathematics On finding a cavity in a thermoelastic body using a single displacement measurement over a finite time interval on the surface of the body
Article
Licensed
Unlicensed Requires Authentication

On finding a cavity in a thermoelastic body using a single displacement measurement over a finite time interval on the surface of the body

  • Masaru Ikehata ORCID logo EMAIL logo
Published/Copyright: January 30, 2018

Abstract

A mathematical formulation of an estimation problem of a cavity inside a three-dimensional thermoelastic body by using time domain data is considered. The governing equation of the problem is given by a system of equations in the linear theory of thermoelasticity which is a coupled system of the elastic wave and heat equations. A new version of the enclosure method in the time domain which is originally developed for the classical wave equation is established. For a comparison, the results in the decoupled case are also given.

Funding statement: The author was partially supported by Grant-in-Aid for Scientific Research (C) (no. 17K05331) of Japan Society for the Promotion of Science. Some part of this work was initiated during the author’s visit at ICUB (Research Institute of the University of Bucharest), Bucharest, Romania in November 2016, with support by ICUB’s visiting professors fellowship program.

A Appendix

A.1 Explicit computation of integrals of (2.38) and (2.41)

The following formulae have been derived in [22, (A.5), (A.6)]. Note that the first formula of (A.1) below is also a result of the mean value theorem for the modified Helmholtz equation [2].

Proposition A.1 ([22]).

Let xR3B¯. Then we have

(A.1){I0(x;τ)=4πφ0(τη)τ3e-τ|x-p||x-p|,I1(x;τ)=4πφ1(τη)τ4e-τ|x-p||x-p|,

where

(A.2){φ0(s)=scoshs-sinhs,φ1(s)=(s2+2)coshs-2ssinhs-2.

In this subsection, we add the following formulae.

Proposition A.2.

We have

(A.3){𝑰0(x;τ)=πe-τ|x-p||x-p|2Kτ0(|x-p|,τη)x-p|x-p|,𝑰1(x;τ)=4πτ2e-τ|x-p||x-p|2Kτ1(|x-p|,τη)x-p|x-p|,I2(x;τ)=4πφ2(τη)τ5e-τ|x-p||x-p|,

where

Kτ0(ξ,s)=2τ{(1-1τ)ξ2-2τ2ξ-2τ3}(coshs-1)+2τ3(1-1τ){(s2+2)coshs-2ssinhs-2}
(A.4)+4τ3(ssinhs-coshs),
Kτ1(ξ,s)={(1-1τ)ξ2-2τ2(ξ+1τ)}(scoshs-sinhs)+1τ2(1-1τ){s2(s+6)coshs-3(s2+2)sinhs}
(A.5)+2τ2(ξ+1τ){(s2+1)sinhs-2scoshs}

and

(A.6)φ2(s)=s2(s+6)coshs-3(s2+2)sinhs.

Proof.

It suffices to consider the case when p=0. Let 𝑨x be the orthogonal matrix such that (𝑨x)Tx=|x|𝒆3. The change of variables y=rω(0<r<η,ωS2) and a rotation give us

𝑰0(x;τ)=0ηr2𝑑rS2e-τ|x-rω||x-rω|ω𝑑ω
=0ηr2𝑑rS2e-τ||x|𝒆3-rω|||x|𝒆3-rω|𝑨xω𝑑ω

=0ηr2𝑑r02π𝑑θ0πsinφdφe-τ|x|2-2r|x|cosφ+r2|x|2-2r|x|cosφ+r2𝑨x(sinφcosθsinφsinθcosφ)
=2π0ηr2𝑑r0πsinφdφe-τ|x|2-2r|x|cosφ+r2|x|2-2r|x|cosφ+r2cosφ𝑨x𝒆3
(A.7)=2π0ηU(|x|,r)r2𝑑rx|x|,

where

U(ξ,r)=0πe-τξ2-2rξcosφ+r2ξ2-2rξcosφ+r2sinφcosφdφ,ξ>η, 0<r<η.

Fix ξ]η,[ and r]0,η[. The change of variable

s=ξ2-2rξcosφ+r2,φ]0,π[,

gives

sinφcosφdφ=ξ2+r2-s22r2ξ2sds.

Thus we have

U(ξ,r)=12r2ξ2ξ-rξ+re-τs(ξ2+r2-s2)𝑑s=12r2ξ2{(ξ2+r2)(-e-τ(ξ+r)+e-τ(ξ-r))-ξ-rξ+re-τss2𝑑s}.

Using the formula

ξ-rξ+re-τss2𝑑s=-1τ{(ξ+r)2+2τ(ξ+r)+2τ2}e-τ(ξ+r)+1τ{(ξ-r)2+2τ(ξ-r)+2τ2}e-τ(ξ-r),

we obtain

U(ξ,r)=e-τ,ξ2ξ2r2(A(ξ,r)eτr-B(ξ,r)e-τr)

and

A(ξ,r)=(ξ2+r2)-1τ(ξ-r)2-2τ2(ξ-r)-2τ3,
B(ξ,r)=(ξ2+r2)-1τ(ξ+r)2-2τ2(ξ+r)-2τ3.

From this and (A.7) we obtain

(A.8)𝑰0(x;τ)=πe-τξ0η(A(ξ,r)eτr-B(ξ,r)e-τr)dr|ξ=|x|x|x|3.

We have

0ηe±τrr𝑑r=±1τ(η1τ)e±τη+1τ2,
0ηe±τrr2𝑑r=±1τ(η22ητ+2τ2)e±τη2τ3,
0ηe±τrr3𝑑r=±1τe±τηη3-3τ2(η22ητ+2τ2)e±τη+6τ4,

and thus

(A.9){0η(eτr-e-τr)r𝑑r=2τ2(τηcoshτη-sinhτη),0η(eτr+e-τr)r𝑑r=2τ2(τηsinhτη-coshτη),0η(eτr+e-τr)r2𝑑r=2τ3{(η2τ2+1)sinhτη-2ητcoshτη},0η(eτr-e-τr)r2𝑑r=2τ3{(η2τ2+2)coshτη-2ητsinhτη-2},0η(eτr-e-τr)r3𝑑r=2τ4{η2τ2(ητ+6)coshτη-3(η2τ2+2)sinhτη}.

Now writing

(A.10){A(ξ,r)=(1-1τ)ξ2-2τ2ξ-2τ3+(1-1τ)r2+2τ(ξ+1τ)r,B(ξ,r)=(1-1τ)ξ2-2τ2ξ-2τ3+(1-1τ)r2-2τ(ξ+1τ)r,

and using (A.9), we obtain

0η(A(ξ,r)eτr-B(ξ,r)e-τr)𝑑r
={(1-1τ)ξ2-2τ2ξ-2τ3}0η(eτr-e-τr)𝑑r+(1-1τ)0η(eτr-e-τr)r2𝑑r+2τ0η(eτr+e-τr)r𝑑r
=2τ{(1-1τ)ξ2-2τ2ξ-2τ3}(coshτη-1)+2τ3(1-1τ){(η2τ2+2)coshτη-2ητsinhτη-2}
   +4τ3(τηsinhτη-coshτη).

A combination of this and (A.8) gives the desired formula for 𝑰0(x;τ).

Using the same changes of variable as used in the computation of 𝑰0(x;τ), we have

(A.11)𝑰1(x;τ)=πe-τξ0η(A(ξ,r)eτr-B(ξ,r)e-τr)rdr|ξ=|x|x|x|3.

Then from (A.9) and (A.10) we obtain

0η(A(ξ,r)eτr-B(ξ,r)e-τr)r𝑑r=2τ2{(1-1τ)ξ2-2τ2ξ-2τ3}(τηcoshτη-sinhτη)
+2τ4(1-1τ){η2τ2(ητ+6)coshτη-3(η2τ2+2)sinhτη}
+4τ4(ξ+1τ){(η2τ2+1)sinhτη-2ητcoshτη}.

A combination of this and (A.11) gives the dersired formula for 𝑰1(x;τ).

Similarily to the above, one obtains the expression

I2(x;τ)=2πτe-τ|x||x|0η(eτr-e-τr)r3𝑑r.

Now (A.9) gives the desired formula for I2(x;τ). ∎

A.2 Lower estimates for volume integrals

Recall that B is an open ball centered at p with radius η and satisfies B¯D¯=; we have dist(D,B)=dD(p)-η. In the following lemma, it is assumed that D is C2.

Lemma A.3.

There exist positive constants C, τ0 and κ such that

(A.12)τκe2τρ/μdist(D,B)De-2τρ/μ(|x-p|-η)|x-p|x-p|×𝒂|2𝑑xC

for all ττ0.

Proof.

Choose a point qD such that |q-p|=dD(p). Since D is C2, one can find an open ball B with radius δ and centered at q-δ𝝂q such that BD and BD={q}. Then dist(B,B)=dist(D,B). Thus, it suffices to prove (A.12) in the case when D=B.

Set d=dD(p) and -𝝂q=𝝎0. Let B′′ be the open ball with radius d+δ centered at p. First we give a parametrization of the domain B′′B. Given s]0,δ[, we find the set of all unit vectors 𝝎 satisfying p+(d+s)𝝎B. Since the center of B has the expression p+(d+δ)𝝎0, this condition is equivalent to the equation

|(d+s)𝝎-(d+δ)𝝎0|<δ.

This is equivalent to the condition

𝝎𝝎0>(d+δ)2+(d+s)2-δ22(d+s)(d+δ).

Thus we have

B′′B=0<s<δ{p+(d+s)ωωS(s)},

where

S(s)={𝝎S2|𝝎𝝎0>(d+δ)2+(d+s)2-δ22(d+s)(d+δ)}.

Choose two linearly independent vectors 𝒃 and 𝒄 in such a way that 𝒃𝒄=0 and 𝒃×𝒄=𝝎0. We denote by θ(s)]0,π2[ the unique solution of

(A.13)cosθ=(d+δ)2+(d+s)2-δ22(d+s)(d+δ).

Given s]0,δ[, r]0,(d+s)sinθ(s)[ and γ[0,2π[, set

𝚼(s,r,γ)=p+(d+s)cosθ(s)𝝎0+r(cosγ𝒃+sinγ𝒄)+h𝝎0,

where h is an unknown parameter to be determined by the equation

|𝚼(r,s,γ)-p|=d+s.

By solving this, we obtain

h=-(d+s)cosθ(s)+(d+s)2-r2.

Thus we have

𝚼(s,r,γ)=p+(d+s)2-r2𝝎0+r(cosγ𝒃+sinγ𝒄).

Note that |𝚼(s,r,γ)-p|=d+s and the unit vector (𝚼(s,r,γ)-p)/(d+s) belongs to S(s) for each fixed s. It is easy to check also that the map

G(s,r,γ)𝚼(s,r,γ)B′′B

is bijective, where

G={(s,r,γ)(s,r)G,γ[0,2π[}

and

G={(s,r)0<s<δ, 0<r<(d+s)sinθ(s)}.

A simple computation gives

det𝚼(s,r,γ)=r(d+s)(d+s)2-r2.

Here we have

(A.14)(𝚼(s,r,γ)-p)×𝒂=(d+s)2-r2𝝎0×𝒂+r(cosγ𝒃×𝒂+sinγ𝒄×𝒂)𝑨(s,r,γ)𝒂.

Using the change of variables formula and B′′BB, we obtain

e2τρ/μ(d-η)Be-2τρ/μ(|x-p|-η)|x-p|x-p|×𝒂|2𝑑x
(A.15)0δ𝑑s0(d+s)sinθ(s)𝑑r02π𝑑γe-2τρ/μsd+s|𝑨(s,r,γ)𝒂|2r(d+s)2-r2.

It is easy to see that from (A.14) one gets

02π|𝑨(s,r,γ)𝒂|2𝑑γ=2π{(d+s)2-r2}|𝝎0×𝒂|2+πr2(|𝒃×𝒂|2+|𝒄×𝒂|2).

Also,

0(d+s)sinθ(s)r(d+s)2-r2𝑑r=(d+s)33(1-cos3θ(s)),
0(d+s)sinθ(s)r3dr(d+s)2-r2=(d+s)3{(1-cosθ(s))+1-cos3θ(s)3}.

Thus,

0(d+s)sinθ(s)rdr(d+s)2-r202π|𝑨(s,r,γ)𝒂|2𝑑γ
=2π3(d+s)3(1-cos3θ(s))|𝝎0×𝒂|2+π(d+s)3{(1-cosθ(s))-1-cos3θ(s)3}(|𝒃×𝒂|2+|𝒄×𝒂|2).

From (A.13) we have

(d+s)(1-cosθ(s))=s(2δ-s)2(d+δ)δs2(d+δ).

This gives

(d+s)3(1-cos3θ(s))δ(d+s)2s2(d+δ)(1+cosθ(s)+cos2θ(s))δd(d+s)s2(d+δ).

Also,

(d+s)3{(1-cosθ(s))-1-cos3θ(s)3}δ(d+s)2s2(d+δ)(1-1+cosθ(s)+cos2θ(s)3)
=δ(d+s)2s6(d+δ)(1-cosθ(s))(2+cosθ(s))
δ2(d+s)s212(d+δ)2(2+cosθ(s))
δ2(d+s)s26(d+δ)2.

Thus we have

0s𝑑se-2τρ/μsd+s0(d+s)sinθ(s)rdr(d+s)2-r202π|𝑨(s,r,γ)𝒂|2𝑑γ
(A.16)π3δdd+δ0δse-2τρ/μs𝑑s|𝝎0×𝒂|2+π6(δd+δ)20δs2e-2τρ/μs𝑑s(|𝒃×𝒂|2+|𝒄×𝒂|2).

Here we have, for j=1,2,

0δsje-2τρ/μs𝑑s=1(2τρ/μ)j+1+O(τ-1e-2τρ/μδ).

Therefore, from these and (A.15) and (A.16) one can conclude that if 𝝎0×𝒂𝟎, then one can choose κ=2 in (A.12), and if 𝝎0×𝒂=𝟎, then 𝒂=±𝝎0 and thus |𝒃×𝒂|2+|𝒄×𝒂|2>0, and one can choose κ=3 in (A.12). ∎

Lemma A.4.

There exist positive constants C and τ0 such that

τ2e2τdist(D,B)De-2τ(|x-p|-η)𝑑xC

for all ττ0.

Proof.

From the proof of Lemma A.3 we have

e2τ(d-η)De-2τ(|x-p|-η)𝑑xe2τ(d-η)Be-2τ(|x-p|-η)𝑑x
2π0δ𝑑s0(d+s)sinθ(s)𝑑re-2τsr(d+s)(d+s)2-r2
=2π0δ(d+s)e-2τs𝑑s0(d+s)sinθ(s)r(d+s)2-r2𝑑r
=π0δ(d+s)2(1-cosθ(s))e-2τs𝑑s
π0δ(d+s)δs2(d+δ)e-2τs𝑑s
πδd2(d+δ)0δse-2τs𝑑s.

Since

0δse-2τs𝑑s=14τ202τδξe-ξ𝑑ξ,

we obtain the desired conclusion.∎

Acknowledgements

The author thanks Hiromichi Itou for pointing out the reference [4] and for useful discussions.

References

[1] D. E. Carlson, Linear thermoelasticity, Mechanics of Solids II, Springer, Berlin (1984), 297–345. 10.1007/978-3-642-69567-4_2Search in Google Scholar

[2] R. Courant and D. Hilbert, Methoden der Mathematischen Physik. II, Springer, Berlink, 1937. 10.1007/978-3-642-47434-7Search in Google Scholar

[3] C. M. Dafermos, On the existence and the asymptotic stability of solutions to the equations of linear thermoelasticity, Arch. Ration. Mech. Anal. 29 (1968), 241–271. 10.1007/BF00276727Search in Google Scholar

[4] W. Dan, On the Neumann problem of linear hyperbolic parabolic coupled systems, Tsukuba J. Math. 18 (1994), no. 2, 371–410. 10.21099/tkbjm/1496162608Search in Google Scholar

[5] R. Dautray and J.-L. Lions, Mathematical Analysis and Numerical Methods for Science and Technology. Vol. 5: Evolution Problems I, Springer, Berlin, 1992. Search in Google Scholar

[6] G. Duvaut and J.-L. Lions, Inequalities in Mechanics and Physics, Grundlehren Math. Wiss. 219, Springer, Berlin, 1976. 10.1007/978-3-642-66165-5Search in Google Scholar

[7] M. E. Gurtin, The linear theory of elasticity, Mechanics of Solids II, Springer, Berlin (1984), 1–295. 10.1007/978-3-662-39776-3_1Search in Google Scholar

[8] M. Ikawa, Mixed problems for hyperbolic equations of second order, J. Math. Soc. Japan 20 (1968), 580–608. 10.2969/jmsj/02040580Search in Google Scholar

[9] M. Ikehata, Enclosing a polygonal cavity in a two-dimensional bounded domain from Cauchy data, Inverse Problems 15 (1999), no. 5, 1231–1241. 10.1088/0266-5611/15/5/308Search in Google Scholar

[10] M. Ikehata, Reconstruction of the support function for inclusion from boundary measurements, J. Inverse Ill-Posed Probl. 8 (2000), no. 4, 367–378. 10.1515/jiip.2000.8.4.367Search in Google Scholar

[11] M. Ikehata, Extracting discontinuity in a heat conductive body. One-space dimensional case, Appl. Anal. 86 (2007), no. 8, 963–1005. 10.1080/00036810701460834Search in Google Scholar

[12] M. Ikehata, The enclosure method for inverse obstacle scattering problems with dynamical data over a finite time interval, Inverse Problems 26 (2010), no. 5, Article ID 055010. 10.1088/0266-5611/26/5/055010Search in Google Scholar

[13] M. Ikehata, The framework of the enclosure method with dynamical data and its applications, Inverse Problems 27 (2011), no. 6, Article ID 065005. 10.1088/0266-5611/27/6/065005Search in Google Scholar

[14] M. Ikehata, An inverse acoustic scattering problem inside a cavity with dynamical back-scattering data, Inverse Problems 28 (2012), no. 9, Article ID 095016. 10.1088/0266-5611/28/9/095016Search in Google Scholar

[15] M. Ikehata, The enclosure method for inverse obstacle scattering problems with dynamical data over a finite time interval: II. Obstacles with a dissipative boundary or finite refractive index and back-scattering data, Inverse Problems 28 (2012), no. 4, Article ID 045010. 10.1088/0266-5611/28/4/045010Search in Google Scholar

[16] M. Ikehata, The enclosure method for inverse obstacle scattering problems with dynamical data over a finite time interval: III. Sound-soft obstacle and bistatic data, Inverse Problems 29 (2013), no. 8, Article ID 085013. 10.1088/0266-5611/29/8/085013Search in Google Scholar

[17] M. Ikehata, Extracting the geometry of an obstacle and a zeroth-order coefficient of a boundary condition via the enclosure method using a single reflected wave over a finite time interval, Inverse Problems 30 (2014), no. 4, Article ID 045011. 10.1088/0266-5611/30/4/045011Search in Google Scholar

[18] M. Ikehata, On finding an obstacle embedded in the rough background medium via the enclosure method in the time domain, Inverse Problems 31 (2015), no. 8, Article ID 085011. 10.1088/0266-5611/31/8/085011Search in Google Scholar

[19] M. Ikehata, The enclosure method for inverse obstacle scattering using a single electromagnetic wave in time domain, Inverse Probl. Imaging 10 (2016), no. 1, 131–163. 10.3934/ipi.2016.10.131Search in Google Scholar

[20] M. Ikehata, A remark on finding the coefficient of the dissipative boundary condition via the enclosure method in the time domain, Math. Methods Appl. Sci. 40 (2017), no. 4, 915–927. 10.1002/mma.4021Search in Google Scholar

[21] M. Ikehata, On finding an obstacle with the Leontovich boundary condition via the time domain enclosure method, Inverse Probl. Imaging 11 (2017), no. 1, 99–123. 10.3934/ipi.2017006Search in Google Scholar

[22] M. Ikehata, The enclosure method for inverse obstacle scattering over a finite time interval: IV. Extraction from a single point on the graph of the response operator, J. Inverse Ill-Posed Probl. 25 (2017), no. 6, 747–761. 10.1515/jiip-2016-0023Search in Google Scholar

[23] M. Ikehata and H. Itou, On reconstruction of a cavity in a linearized viscoelastic body from infinitely many transient boundary data, Inverse Problems 28 (2012), no. 12, Article ID 125003. 10.1088/0266-5611/28/12/125003Search in Google Scholar

[24] M. Ikehata and M. Kawashita, The enclosure method for the heat equation, Inverse Problems 25 (2009), no. 7, Article ID 075005. 10.1088/0266-5611/25/7/075005Search in Google Scholar

[25] M. Ikehata and M. Kawashita, On the reconstruction of inclusions in a heat conductive body from dynamical boundary data over a finite time interval, Inverse Problems 26 (2010), no. 9, Article ID 095004. 10.1088/0266-5611/26/9/095004Search in Google Scholar

[26] M. Ikehata and M. Kawashita, An inverse problem for a three-dimensional heat equation in thermal imaging and the enclosure method, Inverse Probl. Imaging 8 (2014), no. 4, 1073–1116. 10.3934/ipi.2014.8.1073Search in Google Scholar

[27] Y. Shibata, On the Neumann problem for some linear hyperbolic systems of second order, Tsukuba J. Math. 12 (1988), no. 1, 149–209. 10.21099/tkbjm/1496160641Search in Google Scholar

[28] K. Yosida, Functional Analysis, 3rd ed., Grundlehren Math. Wiss. 123, Springer, New York, 1971. 10.1007/978-3-662-00781-5Search in Google Scholar

Received: 2017-6-25
Accepted: 2018-1-6
Published Online: 2018-1-30
Published in Print: 2018-6-1

© 2018 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 10.12.2025 from https://www.degruyterbrill.com/document/doi/10.1515/jiip-2017-0066/html
Scroll to top button