Startseite Mathematik Lipschitz continuity of the Fréchet gradient in an inverse coefficient problem for a parabolic equation with Dirichlet measured output
Artikel
Lizenziert
Nicht lizenziert Erfordert eine Authentifizierung

Lipschitz continuity of the Fréchet gradient in an inverse coefficient problem for a parabolic equation with Dirichlet measured output

  • Alemdar Hasanov EMAIL logo
Veröffentlicht/Copyright: 25. Januar 2018

Abstract

This paper studies the Lipschitz continuity of the Fréchet gradient of the Tikhonov functional J(k):=(1/2)u(0,;k)-fL2(0,T)2 corresponding to an inverse coefficient problem for the 1D parabolic equation ut=(k(x)ux)x with the Neumann boundary conditions -k(0)ux(0,t)=g(t) and ux(l,t)=0. In addition, compactness and Lipschitz continuity of the input-output operator

Φ[k]:=u(x,t;k)|x=0+,Φ[]:𝒦H1(0,l)H1(0,T),

as well as solvability of the regularized inverse problem and the Lipschitz continuity of the Fréchet gradient of the Tikhonov functional are proved. Furthermore, relationships between the sufficient conditions for the Lipschitz continuity of the Fréchet gradient and the regularity of the weak solution of the direct problem as well as the measured output f(t):=u(0,t;k) are established. One of the derived lemmas also introduces a useful application of the Lipschitz continuity of the Fréchet gradient. This lemma shows that an important advantage of gradient methods comes when dealing with the functionals of class C1,1(𝒦). Specifically, this lemma asserts that if JC1,1(𝒦) and {k(n)}𝒦 is the sequence of iterations obtained by the Landweber iteration algorithm k(n+1)=k(n)+ωnJ(k(n)), then for ωn(0,2/Lg), where Lg>0 is the Lipschitz constant, the sequence {J(k(n))} is monotonically decreasing and limnJ(k(n))=0.

MSC 2010: 35R30; 49J50; 35D30

Dedicated to Prof. Marian Slodička on his 60th birthday


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 the direct problem (1.1).

A.1 Estimates for the weak solution

First of all, we derive a useful trace norm estimate for the weak solution uL2(0,T;H1(0,l)) of the direct problem (1.1) that can be treated as a one-dimensional analogue of the trace norm estimate in [9, Section 1.6, (6.15)].

Lemma A.1.

Let conditions (1.3) hold. Then for the trace norm u(0,)L2(0,T) of the weak solution of the direct problem (1.1) the following estimate holds:

(A.1)0tu2(0,τ)𝑑τ2l0l0tu2(x,τ)𝑑τ𝑑x+l0t0lux2(x,τ)𝑑x𝑑τ.

Proof.

We use the identity

u2(0,t)=(u(x,t)-0xuξ(ξ,t)𝑑ξ)2,t(0,T).

Applying to the right-hand side the inequality (a-b)22(a2+b2), then using the Hölder inequality and integrating on [0,t], t[0,T], we deduce that

0tu2(0,τ)𝑑τ20tu2(x,τ)𝑑τ+2x0t0lux2(x,τ)𝑑x𝑑τ.

Integrating then on [0,l] and dividing both sides by l>0, we arrive at the required estimate (A.1). ∎

Lemma A.2.

Let conditions (1.3) hold. Then for the weak solution uL2(0,T;H1(0,l)) of the direct problem (1.1) with utL2(0,T;H-1(0,l)) the following estimates hold:

(A.2)uC([0,T];L2(0,l))C1gL2(0,T),
(A.3)uxL2(0,T;L2(0,l))C2gL2(0,T),

where

(A.4)C1=(l/(c0))1/2exp(Tc0/l2),C2=(2TC12/l2+l/c02)1/2

and c0>0 is the constant defined in (1.3).

Proof.

Multiply both sides of equation (1.1) by u(x,t), integrate on [0,l] and use the integration by parts formula. Taking into account the boundary conditions in (1.1), integrating then both sides on [0,t], t[0,T], and using the initial condition u(x,0)=0, we obtain the following energy identity:

0lu2(x,t)𝑑x+20t0lk(x)ux2(x,τ)𝑑x𝑑τ=20tg(τ)u(0,τ)𝑑τ

for a.e. t[0,T]. We apply the ε-inequality 2ab(1/ε)a2+εb2 in the right-hand-side integral to get

(A.5)0lu2(x,t)𝑑x+2c00t0lux2(x,τ)𝑑x𝑑τ1ε0tg2(τ)𝑑τ+ε0tu2(0,τ)𝑑τ

for a.e. t[0,T]. To estimate the last integral in (A.5) we use estimate (A.1). Then we conclude that

(A.6)0lu2(x,t)𝑑x+(2c0-lε)0t0lux2(x,τ)𝑑x𝑑τ2εl0l0tu2(x,τ)𝑑τ𝑑x+1ε0tg2(τ)𝑑τ

for a.e. t[0,T] and ε(0,2c0/l). Taking for convenience this parameter as ε=c0/l, we obtain

(A.7)0lu2(x,t)𝑑x+c00t0lux2(x,τ)𝑑x𝑑τ2c0l20l0tu2(x,τ)𝑑τ𝑑x+lc00tg2(τ)𝑑τ.

The first consequence of (A.7) is the inequality

0lu2(x,t)𝑑x2c0l20t0lu2(x,τ)𝑑x𝑑τ𝑑x+lc00tg2(τ)𝑑τ

for a.e. t[0,T]. Applying the Gronwall–Bellman inequality [1], we deduce that

(A.8)0lu2(x,t)𝑑xC120tg2(τ)𝑑τ,t[0,T],

where C1>0 is the constant defined in (A.4). Remember that according to [5, Section 5.9, Theorem 3] the weak solution uL2(0,T;H1(0,l)) with utL2(0,T;H-1(0,l)) is in C([0,T];L2(0,l)). This implies estimate (A.2).

The second consequence of (A.7) is the inequality

c00T0lux2(x,t)𝑑x𝑑t2c0l20T0lu2(x,t)𝑑x𝑑t𝑑x+lc00Tg2(t)𝑑t.

By estimate (A.8), this inequality implies the second required estimate (A.3). ∎

Using estimates (A.3) and (A.8) in (A.1), we can estimate the L2-norm of the output u(0,t).

Corollary A.3.

Let the conditions of Lemma A.2 hold. Then

(A.9)u(0,)L2(0,T)C3gL2(0,T),C3=(2TC12/l+lC22)1/2,

where C1>0 and C2>0 are the constants defined in (A.4).

Corollary A.4.

Let the conditions of Lemma A.2 hold. Then for the weak solution of the adjoint problem (5.3) the following estimates hold:

(A.10)φL2(0,T;L2(0,l))28TC12[C32gL2(0,T)2+fL2(0,T)2],
(A.11)φxL2(0,T;L2(0,l))28C22[C32gL2(0,T)2+fL2(0,T)2],

where C1,C2,C3>0 are the constants defined in (A.4) and (A.9).

Proof.

Using the transformation for the time variable τ=T-t in the adjoint (backward) problem (5.3), we deduce that it becomes the well-posed problem

(A.12){ψτ-(k(x)ψx)x=0,(x,τ)(0,l)×[0,T),ψ(x,0)=0,x(0,l),-k(0)ψx(0,τ)=q^(τ),ψx(l,τ)=0,τ(0,T),

for the function ψ(x,τ):=φ(x,T-t), where q^(τ)=-2[u(0,τ;k)-f(τ)] and q^(0)=0. Applying Lemma A.1 to the well-posed problem (A.12), we conclude that

ψL2(0,T;L2(0,l))24TC12u(0,)-fL2(0,T)2,
ψxL2(0,T;L2(0,l))24C22u(0,)-fL2(0,T)2.

Using here the inequality (a-b)2(a2+b2), applying estimate (A.9) to the norm u(0,)L2(0,T) and then returning to the original function φ(x,t), we arrive at the required estimates (A.10) and (A.11). ∎

Now we derive some estimates for the weak solution v(x,t):=u(x,t;k1)-u(x,t;k2) of problem (2.7).

Lemma A.5.

Let k1,k2L(0,l) and v(x,t):=u(x,t;k1)-u(x,t;k2). Assume that conditions (1.3) hold. Then for the weak solution vL2(0,T;H1(0,l)) of problem (2.7) the following estimates hold:

(A.13)vC([0,T];L2(0,l))C2c0gL2(0,T)k1-k2L(0,l),
(A.14)vxL2(0,T;L2(0,l))C2c0gL2(0,T)k1-k2L(0,l),

with the constants defined in Lemma A.1.

Proof.

As in the proof of Lemma A.2, we multiply both sides of equation (2.7) by v(x,t), integrate on Ωt, t(0,T], and use the integration by parts formula. Then we obtain the following integral identity:

120lv2(x,t)𝑑x+Ωtk1(x)vx2(x,τ)𝑑τ𝑑x-0t(k1(x)vx(x,τ)v(x,τ))x=0x=l𝑑τ
=-Ωtδk(x)u2x(x,τ)vx(x,τ)𝑑τ𝑑x+0t(δk(x)u2x(x,τ)v(x,τ))x=0x=l𝑑τ

for all t[0,T], where δk(x)=k1(x)-k2(x). The terms under the last integrals on the left- and right-hand side are zero at x=l due to the homogeneous boundary condition vx(l,t)=0 in (2.7). Taking into account the Neumann boundary condition -k1(0)vx(0,t)=δk(0)u2x(0,t;k2) in the last left-hand-side integral, we deduce that this term and the term at x=0 under the last right-hand-side integral are mutually exclusive. Then applying the ε-inequality to the first right-hand-side integral, after elementary transformations we obtain the following inequality:

0lv2(x,t)𝑑x+(2c0-ε)Ωtvx2(x,τ)𝑑x𝑑τ1εδkL(0,l)2Ωtu2x2(x,τ)𝑑τ𝑑x

for t[0,T], with ε(0,2c0). Choosing this parameter as ε=c0, we get

vL2(0,l)2+c0vxL2(0,t;L2(0,l))21c0δkL(0,l)2u2x2(x,τ)L2(0,t;L2(0,l))2

for a.e. t[0,T]. Using here estimate (A.3), we obtain from this inequality the desired estimates (A.13) and (A.14). ∎

Finally, we derive a necessary estimate for the gradient norm φxL2(0,T;L2(0,l)) of the weak solution of problem (5.4).

Lemma A.6.

Let the conditions of Lemma A.2 hold. Then for the weak solution

δφ(x,t;k):=φ(x,t;k1)-φ(x,t;k2)

of problem (5.4) the following estimate holds:

(A.15)δφxL2(0,T;L2(0,l))2C62[C42gL2(0,T)2+C52fL2(0,T)2]δkL(0,l)2,

where δk(x)=k1(x)-k2(x),

(A.16){C42=C22(8c02C32+2Tc0+l2)/c0,C52=8lC22/c0,C62=((l+2)T/l2)exp((l+2)Tc0/l2)+1/c0

and c0,C2,C3>0 are the constants defined in Lemma A.1 and Corollary A.3.

Proof.

We employ the same transformation τ=T-t used in Corollary A.4, then use the proof scheme of Lemma A.2 for the well-posed problem

(A.17){δψτ-(k1(x)δψx)x=(δk(x)ψ2x(x,τ))x,(x,τ)ΩT,δψ(x,0)=0,x(0,l),-k1(0)δψx(0,τ)=δk(0)ψ2x(0,τ)-δw(0,τ),δψx(l,τ)=0,τ(0,T),

for the function ψ(x,τ):=φ(x,T-t), where δw(0,τ):=δu(0,T-t). Specifically, multiplying both sides of (A.17) by δψ(x,τ), integrating it over (0,l)×(0,τ) and taking into account the boundary conditions in (A.17), we obtain the following integral identity:

0lδψ2(x,τ)𝑑x+20τ0lk1(x)δψx2(x,ζ)𝑑x𝑑ζ
=-20τ0lδk(x)ψ2x(x,ζ)δψ(x,ζ)𝑑x𝑑ζ+20τδw(0,ζ)δψ(0,ζ)𝑑ζ

for a.e. τ(0,T]. We apply the ε-inequality 2ab(1/ε)a2+εb2 to the right-hand-side integrals with subsequent use of inequality (A.1) in the second right-hand-side integral. After elementary transformations, above integral identity becomes the inequality

0lδψ2(x,τ)𝑑x+(2c0-lε)0τ0lδψx2(x,ζ)𝑑x𝑑ζ
(2+l)εl0τ0lδψ2(x,ζ)𝑑x𝑑ζ+lεδkL(0,l)20τ0lψ2x2(x,ζ)𝑑x𝑑ζ
   +2lε0τ0lδw2(x,ζ)𝑑x𝑑ζ+lε0τ0lδwx2(x,ζ)𝑑x𝑑ζ.

Choosing the arbitrary parameter ε(0,2c0/l) as ε=c0/l, we obtain

0lδψ2(x,τ)𝑑x+c00τ0lδψx2(x,ζ)𝑑x𝑑ζ(l+2)c0l20τ0lδψ2(x,ζ)𝑑x𝑑ζ+lc0δkL(0,l)20τ0lψ2x2(x,ζ)𝑑x𝑑ζ
(A.18)+2c00τ0lδw2(x,ζ)𝑑x𝑑ζ+l2c00τ0lδwx2(x,ζ)𝑑x𝑑ζ

for a.e. τ(0,T]. Evidently, at τ=T the second, third and fourth right-hand-side integrals for the functions ψ2x2(x,ζ), δw2(x,ζ) and δwx2(x,ζ) are equal to the integrals for the functions φ2x2(x,t), δu2(x,t) and δux2(x,t), respectively, and for these integrals we can use estimates (A.11), (A.13) and (A.14). Substituting these estimates into (A.18), we conclude that

0lδψ2(x,τ)𝑑x+c00τ0lδψx2(x,ζ)𝑑x𝑑ζ
(A.19)(l+1)c0l20τ0lδψ2(x,ζ)𝑑x𝑑ζ+[C42gL2(0,T)2+C52fL2(0,T)2]δkL(0,l)2,

where C42,C52>0 are defined in (A.16).

As in the proof of Lemma A.2, we first estimate from inequality (A.19) the norm ψL2(0,T;L2(0,l)), then, substituting it in the inequality

c0ψxL2(0,T;L2(0,l))(l+2)c0lψL2(0,T;L2(0,l))+[C42gL2(0,T)2+C52fL2(0,T)2]δkL(0,l)2,

we estimate the norm ψxL2(0,T;L2(0,l)). Returning then to the function φ(x,t), by the above reason we arrive at the required estimate (A.15). ∎

A.2 Estimates for the regular weak solution

We assume here that the inputs k(x) and g(t) in the direct problem (1.1) satisfy the regularity conditions (2.1). As noted in Section 1, under these conditions there exists a unique regular weak solution with improved regularity defined by (2.2).

Lemma A.7.

Let conditions (2.1) hold. Then for the regular weak solution with improved regularity defined by (2.2) the following estimates hold:

(A.20)utL(0,T;L2(0,l))C1gL2(0,T),
(A.21)uxtL2(0,T;L2(0,l))C2gL2(0,T),

where C1,C2>0 are the constants defined in (A.4).

Proof.

Differentiate formally equation (1.1) with respect to t(0,T), multiply both sides by ut(x,t), integrate over [0,l] and use the integration by parts formula. Taking then into account the boundary conditions, we obtain the following integral identity:

12ddt0lut2(x,t)𝑑x+0lk(x)uxt2(x,t)𝑑x=g(t)ut(0,t).

Integrating this identity on [0,t], t[0,T], and using the limit equation

0lut2(x,0+)𝑑x=limt0+0l((k(x)ux(x,0+))x)2𝑑x=0

with the initial condition u(x,0)=0, we deduce that

0lut2(x,t)𝑑x+20t0lk(x)uxτ2(x,τ)𝑑x𝑑τ=20tg(τ)uτ(0,τ)𝑑τ

for all t[0,T]. Applying now the ε-inequality 2ab(1/ε)a2+εb2 to the last right-hand-side integral, we obtain the inequality

0lut2(x,t)𝑑x+2c00t0luxτ2(x,τ)𝑑x𝑑τε0tuτ2(0,τ)𝑑τ+1ε0t(g(τ))2𝑑τ,t[0,T].

Using the same argument as in the proof of inequality (A.1), we estimate the first right-hand-side integral as follows:

(A.22)0tuτ2(0,τ)𝑑τ2l0l0tuτ2(x,τ)𝑑τ𝑑x+l0t0luxτ2(x,τ)𝑑x𝑑τ

for all t[0,T]. Substituting this in the above inequality, we conclude

0lut2(x,t)𝑑x+(2c0-lε)0t0luxτ2(x,τ)𝑑x𝑑τ2εl0t0luτ2(x,τ)𝑑x𝑑τ+1ε0t(g(τ))2𝑑τ,t[0,T].

This is the same inequality (A.6) with u(x,t) replaced by ut(x,t), and g(t) replaced by g(t). Repeating the same procedure as in the proof of Lemma A.1, we deduce from this inequality the required estimates (A.20) and (A.21). ∎

Using this lemma and inequality (A.22), we can estimate the L2-norm of the partial derivative ut(0,t) of the output.

Corollary A.8.

Let the conditions of Lemma A.6 hold. Then

(A.23)ut(0,)L2(0,T)C3gL2(0,T),

where C3>0 is the constant defined in (A.9).

Corollary A.9.

Let the conditions of Lemma A.6 hold. Assume, in addition, that the measured output f(t) satisfies the regularity condition (5.10), that is, fV(0,T). Then for the regular weak solution with improved regularity of the adjoint problem (5.3) the following estimates hold:

(A.24)φtL2(0,T;L2(0,l))22TC12[C32gL2(0,T)2+fL2(0,T)2],
(A.25)φxtL2(0,T;L2(0,l))22C22[C32gL2(0,T)2+fL2(0,T)2].

Proof.

Using the transformation τ=T-t, we can prove in the same way as in Lemma A.5 that estimates similar to (A.20) and (A.21) hold also for the regular weak solution of the adjoint problem (5.3). That is,

ΩTφt2(x,t)𝑑x𝑑tC120T(ut(0,x;k)-f(t))2𝑑t,
ΩTφxt2(x,t)𝑑x𝑑tC220T(ut(0,x;k)-f(t))2𝑑t,

with the same constants C1,C2>0 defined in (A.4). Using here the inequality (a-b)2a2+b2 and estimate (A.23), we arrive at the required estimates (A.24) and (A.25). ∎

Lemma A.10.

Let conditions (2.1) and (5.9) hold. Then for the regular weak solution with improved regularity of problem (1.1) the following estimate holds:

(A.26)uxxL2(0,T;L2(0,l))2C72g𝒱(0,T)2,

where

C72=2(TC12+C22M1)/c02

and C1,C2,M1>0 are the constants defined in (A.4) and (5.9).

Proof.

Let u(x,t) be the regular weak solution with improved regularity defined by (2.2). Then, integrating the formal identity

(k(x)uxx(x,t))2=(ut(x,t)-k(x)ux(x,t))2

over ΩT and using the inequality (a+b)22(a2+b2), we get

c02ΩTuxx2(x,t)𝑑x𝑑t2ΩTut2(x,t)𝑑x𝑑t+2kL(0,l)2ΩTux2(x,t)𝑑x𝑑t.

Taking into account here estimates (A.3) and (A.20), we obtain the required estimate (A.26). ∎

The following assertion is the analogue of the above lemma, and it can be proved in the very same way using estimates (A.11) and (A.24).

Lemma A.11.

Let the conditions of Lemma A.10 hold. Assume, in addition, that the measured output f(t) satisfies the regularity condition (5.10). Then for the regular weak solution with improved regularity of the adjoint problem (5.3) the following estimate holds:

(A.27)φxxL2(0,T;L2(0,l))24(TC12+4C22M1)[C32g𝒱(0,T)2+f𝒱(0,T)2],

where C1,C2,C3,M1>0 are the constants defined by (A.4), (A.9) and (5.9).

Finally, we derive some estimates for the regular weak solutions with improved regularity of problems (2.7) and (5.4).

Lemma A.12.

Let conditions (2.1) hold and let k1,k2K. Denote by v(x,t):=u(x,t;k1)-u(x,t;k2) the regular weak solution of problem (2.7) with improved regularity, defined by (2.2). Then the following estimates hold:

(A.28)vtL(0,T;L2(0,l))C2c0gL2(0,T)k1-k2C[0,l],
(A.29)vxtL2(0,T;L2(0,l))C2c0gL2(0,T)k1-k2C[0,l],

where C2>0 and c0>0 are the constants defined in Lemma A.1.

Proof.

Differentiate formally equation (2.7) with respect to t(0,T), multiply both sides by vt(x,t), integrate over [0,l] and use the integration by parts formula. Then we have

12ddt0lvt2(x,t)𝑑x+0lk1(x)vxt2(x,t)𝑑x=-0lδk(x)u2xt(x,t)vxt(x,t)𝑑x+(δk(x)u2xtvt)x=0x=l+(k1(x)vxtvt)x=0x=l,

where δk(x)=k1(x)-k2(x). Due to the boundary conditions -k1(0)vx(0,t)=δk(0)u2x(0,t), vx(l,t)=0 in (2.7), the last two terms on the right-hand side drop out. Integrating above integral identity over (0,t) and then applying the ε-inequality, after elementary transformations we obtain

0lvt2(x,t)𝑑x+(2c0-ε)0t0lvxτ2(x,τ)𝑑x𝑑τ1εδkC[0,l]0t0lu2xt2(x,τ)𝑑x𝑑τ,

where 2c0-ε>0. Choosing here ε=c0>0 and using estimate (A.21), we arrive at the estimate

(A.30)vtL2(0,l)2+c0vxτL2(0,t;L2(0,l))2C22c0δkC[0,l]2gL2(0,T)2,t(0,T].

This estimate implies the desired results (A.28) and (A.29). ∎

Lemma A.13.

Let conditions (2.1) and (5.9) hold. Then for the regular weak solution with improved regularity of problem (2.7) the following estimate holds:

(A.31)vxxL2(0,T;L2(0,l))C8g𝒱(0,T)k1-k2H1(0,l),

where

C82=4c02[(Tc0+M1)C22/c02+(1+l)C72]

and C7>0 is the constant defined in Lemma A.7.

Proof.

As in the proof of Lemma A.10, we use the integral inequality

c02ΩTvxx2(x,t)dxdt4[ΩTvt2(x,t)dxdt+δkC[0,l]2ΩTu2xx2(x,t)dxdt
(A.32)+δkL2(0,l)2supx(0,l)0Tu2x2(x,t)dt+k1L(0,l)2ΩTvx2(x,t)dt].

To estimate the third right-hand-side integral, we integrate over (0,T) the identity

u2x2(x,t)(xlu2ξξ(ξ,t)𝑑ξ)2,(x,t)ΩT,

which holds due to the boundary condition ux(l,t)=0. Applying then to the right-hand side the Hölder inequality, we get

0Tu2x2(x,t)𝑑tl0T0lv2xx2(x,t)𝑑x𝑑t

for all x(0,l). Taking into account estimate (A.26), we obtain

(A.33)supx(0,l)0Tu2x2(x,t)𝑑tlC72gL2(0,T)2.

To estimate the first right-hand-side integral of (A.32) we use (A.30). We have

(A.34)ΩTvt2(x,t)𝑑x𝑑tTC22c0δkC[0,l]2gL2(0,T)2.

Using in (A.32) estimates (A.33) and (A.34), with estimates (A.26) and (A.14) for the second and fourth right-hand-side integrals, we arrive at estimate (A.31). ∎

Now we prove that similar results to those given in Lemma A.12 and Lemma A.13 also hold for the regular weak solution with improved regularity of problem (5.4), whose Neumann data contains also the increment of the output δ(0,t;k).

Lemma A.14.

Let conditions (2.1) hold. Assume, in addition, that the measured output f(t) satisfies the regularity condition (5.10). Then for the regular weak solution with improved regularity of problem (5.4) the following estimates hold:

(A.35)δφtL2(0,T;L2(0,l))2[C92gL2(0,T)2+C102fL2(0,T)2]k1-k2C[0,l]2,
(A.36)δφxtL2(0,T;L2(0,l))[C112gL2(0,T)2+C122fL2(0,T)2]k1-k2C[0,l]2,

where

(A.37){Ce2=exp(Tβ32),β32=(l+2)c0/l2,β42=(2lc02C32+2Tc0+l2)C22/c03,β52=2lC22/c0,C92=TCe2β42,C102=TCe2β52,C112=(Tβ32Ce2+1)β42/c0,C122=(Tβ32+1)β52/c0,

and the constants C2,C3>0 are defined in Lemma A.1 and Corollary A.3.

Proof.

As in the proof of Lemma A.5, we first employ the transformation τ=T-t to transform the backward problem (5.4) to the well-posed problem. Then we differentiate formally the transformed equation δψτ-(k1(x)δψx)x=(δk(x)ψ2x(x,τ))x with respect to τ(0,T), multiply both sides by δψτ(x,τ), integrate over [0,l] and use the integration by parts formula. Using the proof scheme of Lemma A.5, we obtain the integral inequality

0lδψt2(x,τ)𝑑x+c00τ0lδψxζ2(x,ζ)𝑑x𝑑ζ(l+2)c0l20τ0lδψζ2(x,ζ)𝑑x𝑑ζ+lc0δkC[0,l]0τ0lψ2xζ2(x,ζ)𝑑x𝑑ζ
+2c00τ0lδwζ2(x,ζ)𝑑x𝑑ζ+l2c00τ0lδwxζ2(x,ζ)𝑑x𝑑ζ,

which is similar to the integral inequality (A.18). Use now estimates (A.25), (A.28) and (A.29) for the second, third and fourth right-hand-side integrals for the functions ψ2xζ2(x,ζ), δwζ2(x,ζ) and δwxζ2(x,ζ), which are equal to the integrals for the functions φ2xζ2(x,t), δut2(x,t) and δuxt2(x,t), respectively, at τ=T:

0lδψt2(x,τ)𝑑x+c00τ0lδψxζ2(x,ζ)𝑑x𝑑ζβ320τ0lδψζ2(x,ζ)𝑑x𝑑ζ+[β42gL2(0,T)2+β52fL2(0,T)2]δkC[0,l]2,

where the constants are defined in (A.37).

Estimates (A.35) and (A.36) are obtained from this integral inequality by the same argument as in the proof of Lemma A.2. ∎

Lemma A.15.

Let conditions (2.1) hold. Assume, in addition, that the measured output f(t) satisfies the regularity condition (5.10). Then for the regular weak solution with improved regularity δφ(x,t) of problem (5.4) the following estimate holds:

(A.38)δφxxL2(0,T;L2(0,l))[C132g𝒱(0,T)2+C142f𝒱(0,T)2]k1-k2H1(0,l)2,

where

(A.39){β62=4(l+1)(TC12+4C22+k1L(0,l)2),C132=4(C92+β62C32+C42C62)/c02,C142=4(C102+β62+C52C62)/c02,

and the constants C1,C2,C3,C4,C5,C6,C9,C10>0 are defined in Lemma A.1, Corollary A.3, Lemma A.6 and Lemma A.14.

Proof.

As in the proof of Lemma A.7, we use the integral inequality

c02ΩTδφxx2(x,t)dxdt4[ΩTδφt2(x,t)dxdt+δkC[0,l]2ΩTφ2xx2(x,t)dxdt
(A.40)+δkL2(0,l)2supx(0,l)0Tφ2x2(x,t)dt+k1L(0,l)2ΩTδφx2(x,t)dt]

for the regular weak solution δφ(x,t) of problem (5.4).

In the same way as estimate (A.33), we can estimate the third right-hand-side integral by using estimate (A.27). We have

sup(0,l)0Tφ2x2(x,t)𝑑t4l(TC12+4C22k1L(0,l)2)[C32g𝒱(0,T)2+f𝒱(0,T)2].

Using this estimate with estimates (A.35), (A.27) and (A.15) for the first, second and fourth, respectively, right-hand-side integrals of (A.40), we obtain

c02δφxxL2(0,T;L2(0,l))24{[C92gL2(0,T)2+C102fL2(0,T)2]δkC[0,l]2
+4(TC12+4C22k1L(0,l)2)[C32g𝒱(0,T)2+f𝒱(0,T)2]δkC[0,l]2
+4l(TC12+4C22k1L(0,l)2)[C32g𝒱(0,T)2+f𝒱(0,T)2]δkL2(0,l)2
C62[C42gL2(0,T)2+C52fL2(0,T)2]k1L(0,l)2δkL(0,l)2}.

This leads to the desired estimate (A.38) with the constants C13,C14>0 defined by (A.39).∎

Acknowledgements

The author would like to thank Professors Vladimir G. Romanov, Marian Slodička and Cristiana Sebu for their valuable comments and suggestions.

References

[1] R. Bellman, Asymptotic series for the solutions of linear differential-difference equations, Rend. Circ. Mat. Palermo (2) 7 (1958), 261–269. 10.1007/BF02849324Suche in Google Scholar

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

[3] P. DuChateau, Monotonicity and invertibility of coefficient-to-data mappings for parabolic inverse problems, SIAM J. Math. Anal. 26 (1995), no. 6, 1473–1487. 10.1137/S0036141093259257Suche in Google Scholar

[4] H. W. Engl, M. Hanke and A. Neubauer, Regularization of Inverse Problems, Math. Appl. 375, Kluwer Academic Publishers, Dordrecht, 1996. 10.1007/978-94-009-1740-8Suche in Google Scholar

[5] L. C. Evans, Partial differential equations, 2nd ed., Grad. Stud. Math. 19, American Mathematical Society, Providence, 2010. Suche in Google Scholar

[6] A. Hasanov, P. DuChateau and B. Pektaş, 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), no. 5, 435–463. 10.1515/156939406778247615Suche in Google Scholar

[7] A. Hasanov Hasanoğlu and V. G. Romanov, Introduction to Inverse Problems for Differential Equations, Springer, Cham, 2017. 10.1007/978-3-319-62797-7Suche in Google Scholar

[8] S. I. Kabanikhin, A. Hasanov and A. V. Penenko, The gradient-based method for solving the inverse coefficient heat-conduction problem, Sib. Zh. Vychisl. Mat. 11 (2008), no. 1, 41–54. 10.1134/S1995423908010047Suche in Google Scholar

[9] O. A. Ladyzhenskaya, The Boundary Value Problems of Mathematical Physics, Appl. Math. Sci. 49, Springer, New York, 1985. 10.1007/978-1-4757-4317-3Suche in Google Scholar

[10] L. A. Liusternik and V. J. Sobolev, Elements of Functional Analysis, John Wiley & Sons, New York, 1961. Suche in Google Scholar

[11] E. Zeidler, Nonlinear Functional Analysis and its Applications. II/A: Linear Monotone Operators, Springer, New York, 1990. 10.1007/978-1-4612-0981-2Suche in Google Scholar

Received: 2017-11-18
Accepted: 2017-12-15
Published Online: 2018-1-25
Published in Print: 2018-6-1

© 2018 Walter de Gruyter GmbH, Berlin/Boston

Heruntergeladen am 10.12.2025 von https://www.degruyterbrill.com/document/doi/10.1515/jiip-2017-0106/pdf
Button zum nach oben scrollen