This paper investigates an observability estimate for the parabolic equations with inverse square potential in a C2 bounded domain Ω⊂Rd, which contains 0. The observation region is a product set of a subset E⊂(0,T] with positive measure and a non-empty open subset ω⊂Ω with 0∉ω. We build up this estimate by a delicate result in measure theory in [
Citation: Guojie Zheng, Baolin Ma. Observability estimate for the parabolic equations with inverse square potential[J]. AIMS Mathematics, 2021, 6(12): 13525-13532. doi: 10.3934/math.2021785
[1] | Varadharaj Dinakar, Natesan Barani Balan, Krishnan Balachandran . Identification of Source Terms in a Coupled Age-structured Population Model with Discontinuous Diffusion Coefficients. AIMS Mathematics, 2017, 2(1): 81-95. doi: 10.3934/Math.2017.1.81 |
[2] | Yilihamujiang Yimamu, Zui-Cha Deng, Liu Yang . An inverse volatility problem in a degenerate parabolic equation in a bounded domain. AIMS Mathematics, 2022, 7(10): 19237-19266. doi: 10.3934/math.20221056 |
[3] | Yashar Mehraliyev, Seriye Allahverdiyeva, Aysel Ramazanova . On one coefficient inverse boundary value problem for a linear pseudoparabolic equation of the fourth order. AIMS Mathematics, 2023, 8(2): 2622-2633. doi: 10.3934/math.2023136 |
[4] | Zuliang Lu, Fei Cai, Ruixiang Xu, Chunjuan Hou, Xiankui Wu, Yin Yang . A posteriori error estimates of hp spectral element method for parabolic optimal control problems. AIMS Mathematics, 2022, 7(4): 5220-5240. doi: 10.3934/math.2022291 |
[5] | Hao Zhang, Zejian Cui, Xiangyu Xiao . Decay estimate and blow-up for a fourth order parabolic equation modeling epitaxial thin film growth. AIMS Mathematics, 2023, 8(5): 11297-11311. doi: 10.3934/math.2023572 |
[6] | Dengming Liu, Luo Yang . Extinction behavior for a parabolic $ p $-Laplacian equation with gradient source and singular potential. AIMS Mathematics, 2022, 7(1): 915-924. doi: 10.3934/math.2022054 |
[7] | Jin Zhong, Yilin Zhang . Dual group inverses of dual matrices and their applications in solving systems of linear dual equations. AIMS Mathematics, 2022, 7(5): 7606-7624. doi: 10.3934/math.2022427 |
[8] | Hassen Aydi, Bessem Samet, Manuel De la Sen . On Hermite-Hadamard-type inequalities for second order differential inequalities with inverse-square potential. AIMS Mathematics, 2024, 9(7): 17955-17970. doi: 10.3934/math.2024874 |
[9] | Nataliya Gol'dman . On mathematical models with unknown nonlinear convection coefficients in one-phase heat transform processes. AIMS Mathematics, 2019, 4(2): 327-342. doi: 10.3934/math.2019.2.327 |
[10] | Zui-Cha Deng, Fan-Li Liu, Liu Yang . Numerical simulations for initial value inversion problem in a two-dimensional degenerate parabolic equation. AIMS Mathematics, 2021, 6(4): 3080-3104. doi: 10.3934/math.2021187 |
This paper investigates an observability estimate for the parabolic equations with inverse square potential in a C2 bounded domain Ω⊂Rd, which contains 0. The observation region is a product set of a subset E⊂(0,T] with positive measure and a non-empty open subset ω⊂Ω with 0∉ω. We build up this estimate by a delicate result in measure theory in [
Let Ω be a C2 bounded domain in Rd (d≥3) such that 0∈Ω. Let ω⊂Ω be a non-empty open subset, and 0∉ω. Let T>0, E⊂(0,T] be a Lebesgue measurable subset with positive measure, and denote by χE the characteristic function of E.
In this paper, we study the observability estimate for a parabolic equation with a singular potential term, which is described by
{∂ty(x,t)−△y(x,t)−μy(x,t)|x|2=0, in Ω×(0,T],y(x,t)=0, on ∂Ω×(0,T],y(x,0)=ϕ(x), in Ω, | (1.1) |
where the initial datum ϕ(x)∈L2(Ω), and the parameter μ satisfies
μ<μ∗:=(d−2)24. | (1.2) |
It should be emphasized that condition (1.2) is crucial for Eq (1.1). In [2], it is proved that if the initial datum ϕ(x)∈L2(Ω) is non-negative, then Eq (1.1) admits a unique global solution under condition (1.2), and when μ>μ∗, the local solution may not exist. In recent years, the study on Eq (1.1) has been a hot topic in the theory of PDEs due to its important applications in various branches of applied science and engineering. Extensive related references can be found in [4,10,11] and the rich works cited therein. In particular, the well-posedness of Eq (1.1) is further discussed from the point of view of semigroup theory in [11]. More specifically, it shows that when μ<μ∗, for each ϕ∈L2(Ω), there exists a unique
y∈C([0,+∞);L2(Ω))∩L2(0,+∞;H10(Ω)), |
which is a weak solution of the Eq (1.1).
In this paper, we mainly study the observability estimate for Eq (1.1). The main difficulty in this problem lies in the singularity of Eq (1.1) caused by the potential term. It turns out that Carleman inequalities, frequency functions, and the spectral inequality are the main frameworks in the literature on the observability estimate for the parabolic equations in recent decades. The application of Carleman estimates was described systematically in the works [3,14]. About frequency functions, we would like to mention references [1,7]. In this article, we apply the spectral inequality for the Schrödinger operator A=−△−μ|x|2 to build this estimate.
We denote ‖⋅‖ and ⟨⋅,⋅⟩ to be the usual norm and the inner product in L2(Ω), respectively. The main result of the paper is presented as follows:
Theorem 1.1. Let E⊆(0,T] be a measurable set with a positive measure, and let ω be a non-empty open subset of Ω with 0∉ω. Suppose that
μ≤{72.33, if d=3,(d−1)(d−3)4, if d≥4. | (1.3) |
Then, there exists a constant C=C(Ω,T,ω,E) such that for each datum ϕ∈L2(Ω), the solution to Eq (1.1) satisfied
‖y(x,T)‖2L2(Ω)≤C(Ω,ω,T,E)∫E∫ω|y(x,t)|2dxdt. | (1.4) |
Remark 1.1. Several notes on Theorem 1.1 are given in order.
(a) Inequality (1.4) is called the observability estimate. Observability estimate (1.4) in Theorem 1.1 allows estimating the total energy of the solutions to Eq (1.1) at time T in terms of the partial energy localized in the observation region ω×E.
(b) With the aid of observability estimate (1.4), we can study the bang-bang property for the time optimal control problem governed by Eq (1.1) with internal control.
(c) From the perspective of control theory, this inequality is equivalent to the null controllability for its adjoint equation. (See [5].) We mention [3,12,13] in this direction.
(d) In [16], a similar problem was investigated for Eq (1.1) when Ω is a convex domain of Rd and 0∈ω, where ω is a non-empty open subset of Ω. It obtains in [16] that there exist two positive numbers α=α(Ω,ω)∈(0,1), C=C(Ω,ω) such that for each T>0, the solution of (1.1) satisfies
∫Ω|y(x,T)|2dx≤CeCT(∫Ω|ϕ(x)|2dx)1−α(∫ω|y(x,T)|2dx)α, | (1.5) |
where ϕ(x)∈L2(Ω) is the initial datum. Clearly, (1.5) is a Hölder type interpolation inequality for Eq (1.1), from which and the method in Theorem 4 of [8], we can also get (1.4) in this case.
We organize the paper as follows. In Section 2, we give some preliminary results. Section 3 is devoted to the proof of Theorem 1.1.
We first introduce certain notations. Let {λi}∞i=1, with 0<λ1<λ2≤⋯≤λi→+∞ as i→∞, be the set of eigenvalues of the Schrödinger operator A:=−△−μ|x|2 with homogeneous Dirichlet boundary conditions, and let {ei(x)}∞i=1 be the set of corresponding eigenfunctions with ‖ei(x)‖L2(Ω)=1, (i=1,2,3,⋯) which constitutes an orthonormal basis of L2(Ω) (e.g., see Theorem 3.1 in [12]). Now, for each r>0, we define two subspaces of L2(Ω) as follows
Xr:=span{ei(x)}λi≤r and X⊥r:=span{ei(x)}λi>r. |
Clearly, we have the following orthogonal direct sum decomposition
L2(Ω)=Xr⊕X⊥r. |
Let S(t) be the compact semigroup generated by the operator −A on L2(Ω). Then, the unique solution to Eq (1.1) corresponding to the initial datum ϕ(x)∈L2(Ω) can be written as (e.g., see section 4.1 in [9])
y(x,t)=S(t)ϕ(x). |
In addition,
S(t)ϕ(x)=∞∑n=1e−λnt⟨ϕ,en⟩en(x). | (2.1) |
It follows from (2.1) that
‖S(t)ϕ(x)‖L2(Ω)≤e−λ1t‖ϕ(x)‖L2(Ω). | (2.2) |
Remark 2.1. In fact, it follows from (2.2) that the solution to Eq (1.1) possesses the energy decay property.
As a direct consequence of (2.1) and (2.2), we have
Lemma 2.1. For each positive number r>0 and each η∈X⊥r, the following conclusion is true:
‖S(t)η‖L2(Ω)≤e−rt‖η‖L2(Ω), for t≥0. | (2.3) |
Next, the following result is brought from [6] (see Theorem 1.2 in [6]).
Lemma 2.2. Let ω be a non-empty open subset of Ω with 0∉ω. Suppose (1.3) holds. Then, there exist two positive constants C1>1 and C2>0, which only depend on Ω,ω, such that
∑λi≤r|ai|2≤C1eC2√r∫ω|∑λi≤raiei(x)|2dx | (2.4) |
for each positive number r>0 and each choice of the coefficients {ai}λi≤r with ai∈R.
Remark 2.2. (i) Inequality (2.4) is called the spectral inequality for the Schrödinger operator A=−△−μ|x|2.
(ii) By Lemma 2.2, we have that for each positive number r>0 and each ψ∈Xr,
‖ψ‖L2(Ω)≤C1eC2√r‖ψ‖L2(ω). | (2.5) |
Throughout the rest of this paper, the following notation will be used. For each measurable set A⊆R, we denote by |A| its Lebesgue measure in R. The following lemma is quoted from [7,15] (e.g., see Proposition 2.1 in [7]).
Lemma 2.3. Let E⊆[0,T] be a measurable set with a positive measure, and let l be a density point of E. Then, for each z>1, there exists l1∈(l,T) such that the sequence {li}∞i=1 given by
li+1:=l+1zi(l1−l), | (2.6) |
satisfies
|E∩(li+1,li)|≥12(li−li+1). | (2.7) |
Proof. The proof will be organized in three steps.
Step 1. We first construct a subset sequence of (0,T].
Let l be a density point for E. We arbitrarily fix a positive number z>1. By Lemma 2.3, there exists l1∈(l,T] and a sequence {li}∞i=1 satisfying (2.6) and (2.7). We now define a subset sequence of (0,T] as follows:
Ei:=E∩(li+1+li−li+14,li). | (3.1) |
It is obvious that
|Ei|=|E∩[(li+1,li)∖(li+1,li+1+li−li+14)]|. |
Then, it follow from (2.7) that
|Ei|≥li−li+14. | (3.2) |
It follows from (2.6) that
li−li+1=1zi(z−1)(l1−1),i=1,2,3,⋯. | (3.3) |
Thus,
li−li+1li+1−li+2=z,i=1,2,3,⋯. | (3.4) |
Let b>1 be a positive number satisfying
(bz)i(z−1)(l1−l)>4C2bi+12+4ln(8C1z),i=2,3,4⋯, | (3.5) |
where C1 and C2 are the positive numbers given in Lemma 2.2. Taking ri=bi (i=2,3,4⋯), and using (3.3) and (3.5), we obtain that
rili−li+14−C2√ri+1>ln(8C1z),i=2,3,4,⋯. |
It, together with (3.4), indicates
li+1−li+24C1eC2√ri+1>2(li−li+1)e−rili−li+14,i=2,3,4,⋯. | (3.6) |
Step 2. For each i∈{2,3,4,⋯}, we estimate the term li−li+14‖S(li)ϕ‖L2(Ω), where ϕ∈L2(Ω).
From (3.2) and (2.2), it follows that for each φ∈Xri (i∈{2,3,4⋯}),
li−li+14‖S(li)φ‖L2(Ω)≤∫lili+1+li−li+14χEi(t)‖S(t)φ‖L2(Ω)dt. | (3.7) |
For each initial datum ϕ∈L2(Ω), there exists a unique decomposition
ϕ=ϕ1+ϕ2, | (3.8) |
where ϕ1∈Xri and ϕ2∈X⊥ri. Taking φ=ϕ1 in (3.7), and then using (2.5), we obtain
li−li+14‖S(li)ϕ1‖L2(Ω)≤∫lili+1+li−li+14χEi(t)‖S(t)ϕ1‖L2(Ω)dt≤C1eC2√ri∫lili+1+li−li+14χEi(t)‖S(t)ϕ1‖L2(ω)dt. |
By the triangle inequality and (3.1), we deduce
li−li+14‖S(li)ϕ1‖L2(Ω)≤C1eC2√ri∫lili+1+li−li+14χE(t)(‖S(t)ϕ‖L2(ω)+‖S(t)ϕ2‖L2(Ω))dt. |
Together with (2.2), we have
li−li+14‖S(li)ϕ1‖L2(Ω)≤C1eC2√ri∫lili+1+li−li+14χE(t)‖S(t)ϕ‖L2(ω)dt+C1eC2√ri(li−li+1)‖S(li+1+li−li+14)ϕ2‖L2(Ω). |
This, along with Lemma 2.1, indicates that
li−li+14‖S(li)ϕ1‖L2(Ω)≤C1eC2√ri∫lili+1χE(t)‖S(t)ϕ‖L2(ω)dt+C1eC2√ri(li−li+1)e−rili−li+14‖S(li+1)ϕ2‖L2(Ω). | (3.9) |
By Lemma 2.2 and ϕ2∈X⊥ri, we have
li−li+14‖S(li)ϕ2‖L2(Ω)≤li−li+14e−ri(li−li+1)‖S(li+1)ϕ2‖L2(Ω). | (3.10) |
Applying (3.8), and the triangle inequality, we obtain
li−li+14‖S(li)ϕ‖L2(Ω)≤li−li+14‖S(li)ϕ1‖L2(Ω)+li−li+14‖S(li)ϕ2‖L2(Ω). |
Then, it, together with (3.9), and (3.10), indicates
li−li+14‖S(li)ϕ‖L2(Ω)≤C1eC2√ri∫lili+1χE(t)‖S(t)ϕ‖L2(ω)dt+C1eC2√ri(li−li+1)e−rili−li+14‖S(li+1)ϕ2‖L2(Ω)+li−li+14e−ri(li−li+1)‖S(li+1)ϕ2‖L2(Ω), |
Thus,
li−li+14‖S(li)ϕ‖L2(Ω)≤C1eC2√ri∫lili+1χE(t)‖S(t)ϕ‖L2(ω)dt+(li−li+1)e−rili−li+14(C1eC2√ri+1)‖S(li+1)ϕ‖L2(Ω), |
which leads to
li−li+14C1eC2√ri‖S(li)ϕ‖L2(Ω)−C1eC2√ri+1C1eC2√ri(li−li+1)e−rili−li+14‖S(li+1)ϕ‖L2(Ω)≤∫lili+1χE(t)‖S(t)ϕ‖L2(ω)dt. | (3.11) |
Step 3. We prove (1.4).
Summing (3.11) from 2 to ∞, it follows that
l2−l34C1eC2√r2‖S(l2)ϕ‖L2(Ω)+∞∑i=2ki‖S(li+1)ϕ‖L2(Ω)≤∫T0χE(t)‖S(t)ϕ‖L2(ω)dt, | (3.12) |
where
ki=li+1−li+24C1eC2√ri+1−C1eC2√ri+1C1eC2√ri(li−li+1)e−rili−li+14(i=2,3,4,⋯). |
By (3.6), we obtain
ki>0, for i=2,3,4,⋯. |
This, together with (3.12) and (2.2), indicates that
‖S(T)ϕ‖L2(Ω)≤‖S(l2)ϕ‖L2(Ω)≤4C1eC2√r2l2−l3∫E‖S(t)ϕ‖L2(ω)dt. |
That is
‖y(x,T)‖L2(Ω)≤4C1eC2√r2l2−l3∫E‖y(x,t)‖L2(ω)dt, |
from which and Cauchy's inequality, it follows (1.4). This completes the proof.
Using the Lebeau-Robbiano strategy, we obtain the observation estimate (1.4) for Eq (1.1) in this paper. From (1.4), we have that if the solution y=0 over ω×E, then y(⋅,T)=0 over Ω. It, together with (2.1), shows that y≡0 over Ω×[0,T]. From the application point of view, we can recover the initial state and the evolution history, according to the observation of y in ω×E.
The authors would like to express their sincere thanks to the referees for their valuable suggestions. This work was partially supported by the Natural Science Foundation of Henan Province (No. 202300410248), and by the Natural Science Foundation of Henan Normal University (No. 2019PL15).
The authors declare that they have no competing interests in this paper.
[1] |
J. Apraiz, L. Escauriaza, G. Wang, C. Zhang, Observability inequalities and measurable sets, J. Eur. Math. Soc., 16 (2014), 2433–2475. doi: 10.4171/JEMS/490
![]() |
[2] |
P. Baras, J. A. Goldstein, The heat equation with a singular potential, Trans. Amer. Math. Soc., 284 (1984), 121–139. doi: 10.1090/S0002-9947-1984-0742415-3
![]() |
[3] |
S. Ervedoza, Control and stabilization properties for a singular heat equation with an inverse square potential, Commun. Part. Diff. Eq., 33 (2008), 1996–2019. doi: 10.1080/03605300802402633
![]() |
[4] | J. A. Goldstein, Q. S. Zhang, Linear parabolic equations with strong singular potentials, Trans. Amer. Math. Soc., 355 (2003), 197–211. |
[5] |
J. L. Lions, Exact controllability, stabilization and perturbations for distributed system, SIAM Rev., 30 (1988), 1–68. doi: 10.1137/1030001
![]() |
[6] |
K. D. Phung, Carleman commutator approach in logarithmic convexity for parabolic equations, Math. Control Relat. F., 8 (2018), 899–933. doi: 10.3934/mcrf.2018040
![]() |
[7] |
K. D. Phung, G. Wang, An observability estimate for the parabolic equations from a measurable set in time and its applications, J. Eur. Math. Soc., 15 (2013), 681–703. doi: 10.4171/JEMS/371
![]() |
[8] |
K. D. Phung, L. Wang, C. Zhang, Bang-bang property for time optimal control of semilinear heat equation, Ann. I. H. Poincare-An., 31 (2014), 477–499. doi: 10.1016/j.anihpc.2013.04.005
![]() |
[9] | A. Pazy, Semigroups of linear operators and applications to partial differential equations, Springer-Verlag, New York, 1983. |
[10] |
I. Peral, J. L. Vazquez, On the stability or instability of the singular solution of the semilinear heat equation with exponential reaction term, Arch. Ration. Mech. An., 129 (1995), 201–224. doi: 10.1007/BF00383673
![]() |
[11] |
J. L. Vazquez, E. Zuazua, The Hardy inequality and the asymptotic behaviour of the heat equation with an inverse-square potential, J. Funct. Anal., 173 (2000), 103–153. doi: 10.1006/jfan.1999.3556
![]() |
[12] |
J. Vancostenoble, E. Zuazua, Null controllability for the heat equation with singular inverse square potentials, J. Funct. Anal., 254 (2008), 1864–1902. doi: 10.1016/j.jfa.2007.12.015
![]() |
[13] |
G. Wang, L∞-null controllability for the heat equation and its consequences for the time optimal control problem, SIAM J. Control Optim., 47 (2008), 1701–1720. doi: 10.1137/060678191
![]() |
[14] |
G. Wang, L. Wang, The Carleman inequality and its application to periodic optimal control governed by semilinear parabolic differential equations, J. Optimz. Theory App., 118 (2003), 429–461. doi: 10.1023/A:1025459624398
![]() |
[15] |
C. Zhang, An observability estimate for the heat equation from a product of two measurable sets, J. Math. Anal. Appl., 396 (2012), 7–12. doi: 10.1016/j.jmaa.2012.05.082
![]() |
[16] | G. Zheng, K. Li, Y. Zhang, Quantitative unique continuation for the heat equations with inverse square potential, J. Inequal. Appl., 1 (2018), 1–17. |
1. | Mahmoud Baroun, Hind El Baggari, Ilham Ouled Driss, Said Boulite, Impulse controllability for the heat equation with inverse square potential and dynamic boundary conditions, 2023, 40, 0265-0754, 353, 10.1093/imamci/dnad012 |