In this paper we consider a linear KdV equation with a transport term posed on a finite interval with the boundary conditions considered by Colin and Ghidaglia. The main results concern the behavior of the cost of null controllability with respect to the dispersion coefficient when the control acts on the left endpoint. In particular, for any final time we prove that this cost grows exponentially as the dispersion coefficient vanishes and the transport coefficient is negative.
Let , and . We consider the following controlled Korteweg–de Vries (KdV) equation posed on a finite domain:
where is the dispersion coefficient, is the transport coefficient, is the initial condition and stands for the control. Here and all along the paper we use the notation, for a given function ,
The boundary conditions in (1.1) were proposed by Colin and Ghidaglia in [4] (see also [3]) as a model for propagation of surface water waves in the situation where a wave maker is putting energy on a finite-length channel from the left extremity and the right one is free.
Most results for (1.1) are related to its well-posedness of its nonlinear version (see [13] and the references therein). However, let us observe that if , this system is not dissipative and therefore the existence of solutions needs to be justified. This will be done later in Section 2.
As for the controllability problem, in [2] the authors considered controls in all the boundary conditions and all the possible combinations. For controllability results of the KdV equation in a finite interval with another boundary conditions, see [1,15,16] and the references therein.
In this paper, we are mainly interested in the study of how the size of the control behaves as a function of the dispersion coefficient ε. To this end, we define the quantity
which stands for the cost of null controllability of (1.1). Notice that is the best constant such that, for all and such that , the estimate
is satisfied.
Our first result states an improvement of the cost of the control with respect to the one in [12]. From this work, one can deduce that there exists such that
The first main result of this paper is the following theorem.
Letand. Then, for any, there exists a controlsuch that the associated solution of (1.1) satisfies. Furthermore, we have the estimate if, andif, whereis a constant independent of T, M and ε anddepends at most polynomially on,, T and.
We remark that (1.4)–(1.5) say that the cost of the control is at most of order , whereas in (1.3) is of order . This difference becomes of great importance when studying its behavior in the limit .
Before stating the second and third main results of this paper, let us consider the transport equation
Since (1.6) is controllable if and only if (see, for instance, [5, Theorem 2.6, p. 29]), with a control if and a control if . Furthermore, the cost of null controllability is equal to zero. Indeed, the solution of (1.6) can be brought to zero at time T just by taking when , and by taking when . Thus, we should expect that the cost would decrease to zero in this case as , or at least if the final time T is large enough. On the other hand, if it is expected that the cost of the control would explode as ε tends to zero.
In [10], the authors consider this problem for the classical boundary conditions
and in [9] with controls in all the boundary terms. We refer also to [6] and [11] for the case of vanishing viscosity in one and arbitrary space dimension, respectively.
In these works, the strategy relies on the combination of a suitable Carleman inequality, which gives an observability constant that explodes with ε, with an exponential dissipation estimate for the adjoint equation such that for T large enough counteracts the previous constant. It has been pointed out in [9] and [10] that such a result can only be expected for (1.1) when due to the asymmetric effect of the dispersion term.
We will prove that the cost of null controllability explodes as , even if T is large. Actually, we will prove it for the quantity
where the space is defined as follows: for any , let
endowed with the norm
The following theorem states a lower bound for this new cost of null controllability.
Letand. Then,whereand,andare given by (4.2) and (4.25), respectively.
From the expression of in (4.2), one can see that it depends polynomially on ε and .
From Theorem 1.2, we can deduce the following result that shows that the cost defined in (1.7) blows up as ε vanishes.
Letand. Then, there exists such that for allwe have
Notice that this implies that the cost defined in (1.2) goes to infinity as ε goes to zero since . What is interesting about this result is that it differs from the one obtained in [10] (and [6] for that matter). Notice that for and in (1.1), it seems very difficult (if not impossible) to prove convergence in some sense to a solution of the transport equation (1.6) since we do not know the value at the right-end of the interval . Actually, this convergence question has not been addressed even for the classical boundary conditions (see [10]). Nevertheless, one could have expected to obtain an appropriate dissipation estimate as in [9], but Corollary 1.4 shows that this is not the case.
Finally, when we are able to obtain an explosion result when T is smaller than .
Letand. Then, there exist a constant (independent of ε),and initial conditionssuch that, ifis a control such that the solution y of (1.1) satisfies, then, for every ,In particular,
The rest of the paper is organized as follows. In Section 2 we prove the existence of solutions of equation (1.1). In Section 3, we prove Theorem 1.1. We prove an observability inequality for the adjoint system (3.2) and then applying the Hilbert Uniqueness Method (HUM). The observability inequality is proved by means of a suitable Carleman estimate. In Section 4, we prove Theorem 1.2. In Section 5, we prove Theorem 1.5. Finally, we give the proof of the Carleman estimate in Appendix A.
Existence of solutions when
In this section we will prove the existence (and uniqueness) of a solution of (1.1) when , and . Let us first remark that if y solves (1.1), then
would solve
with
Then, if we prove the existence (and uniqueness) of solution z of (2.1), we would have proved the existence (and uniqueness) of a solution of (1.1) by simply defining
Letand. Then, there exists a unique solution z of (2.1) which belongs to.
We will use the contracting map fixed-point theorem. We consider, for any , the following problem
Notice that the operator defined by
is dissipative with
It is easy to check that the adjoint of A is also dissipative, thus the existence of a unique z solution of (2.4) is ensured.
Now, we multiply the equation in (2.4) by and integrate by parts in space. We obtain
Using Young’s and Poincaré’s inequalities, we get
where depends only on L. We integrate between 0 and t, and take the supremum:
Using Young’s inequality once more and (2.2) we obtain
where C depends only on L.
We consider the norm
Let us now prove the existence of such that the map
is a contraction. In view of (2.5), we have
Let now, for :
with
Choosing β such that
we have for any :
Then, for any :
Therefore, is a contraction mapping on and admits a unique fixed point in which is easy to check that is the solution of (2.1) in . Since β does not depend on v nor , we can repeat this argument in the time intervals , for any , with initial condition , where is the solution of (2.1) in . Thus, and is a solution of (2.1). □
From (2.1)–(2.3), one sees that the solution of (1.1) belongs to . In the sequel, we will also need to establish the existence of solutions of (1.1) when .
Let,and. Then, there exists a unique solution of (1.1).
Let . Then (see for instance [7, Theorem 1, p. 283]), there exist such that
Furthermore,
Then,
where p and q are the solutions of
and
respectively.
From Lemma 2.1 and the computations above, one deduces the desired result. □
Proof of Theorem 1.1
In this section, we will prove Theorem 1.1 by applying the Hilbert Uniqueness Method (HUM) (see, for instance, [14]), that is, we prove the following observability inequality:
Here, φ is the solution of the adjoint equation
with and is a constant independent of φ. Indeed, (3.1) is equivalent to prove that for every , there exists a control such that satisfying
where y is the solution of (1.1). Thus, once an inequality like (3.1) is established, the proof of Theorem 1.1 is finished since
The observability inequality (3.1) is proved by means of Carleman and energy estimates, which are the goals of the following sections.
Furthermore, we notice that from (3.5) and the boundary conditions on in (3.2), for every we have the following initial value ordinary differential equation
We can actually find an explicit formula for φ in terms of ϕ, but we need to distinguish the cases and :
Case
It is not difficult to show that the solution is given by
Thus,
Taking the -norm in this inequality, we get
In this case, φ is given by
The same computations as for the case show that
and
Carleman estimates
To establish the Carleman estimate, we introduce some weight functions. Let
where . Notice that there exist positive constants and that do not depend on T such that
and
for every . This kind of weights has been introduced in [8] and widely used in the literature (in particular, in [9,10,12]).
We now are in position to present our Carleman inequality whose proof is given at the end of the paper (Appendix A).
Let,and. There exists a positive constant C independent of T, ε and M such that, for any solution ϕ of (3.4), we have for all.
Furthermore, we can deduce from Proposition 3.1 and (3.6)–(3.11) a Carleman estimate for the solutions of (3.2).
Let,and. There exists a positive constant C independent of T, ε and M such that, for any solution φ of (3.2), we have for allanddepends at most polynomially on, ε,and.
The lack of homogeneous Dirichlet condition on plays an important role in the choice of the power m of the weight function to prove (3.15). Indeed, a similar inequality was proved in [12] where was needed to estimate a trace term on . From this inequality one can deduce that the cost of the null controllability is bounded by .
Here, by means of the change of variable (3.5), which satisfies , we manage to take the optimal power as in [9,10].
It would be interesting to know if a Carleman estimate can be obtained for the solutions of (3.2) for without using this change of variables.
Dissipation estimates
To prove (3.1), we will combine (3.14) with a dissipation estimate. For , it is easy to check that
for every . For we can prove the following result.
Let. Then, for every pair such thatand every solution ϕ of (3.4) the following inequality is satisfiedwhereis a constant independent of M, T and ε, and
We proceed in two steps. First, we multiply (3.4) by and integrate in . We obtain after some integration by parts
Next, we take the derivative with respect to x of (3.4), multiply by and proceed as before. Straightforward computations lead to
Adding (3.18) and (3.19) we obtain, after neglecting the positive terms on the left-hand side,
Now, notice that
and therefore we obtain
Estimate (3.17) can be easily deduced by integrating in . □
Observability inequality
In this section we combine the Carleman inequality (3.14) and the dissipation estimates obtained in the previous paragraph to finish the proof of (3.1) and therefore the proof of Theorem 1.1. Let us separate the cases and .
In this section we prove Theorem 1.2. For ease of comprehension, we have divided the proof in two parts presented in the following paragraphs.
Previous estimates
Here we will establish several estimates on y (solution of (1.1)) in terms of a constant depending only on T, L, ε and M times the norm of v and the norm of . The principal result of this part is the following lemma.
For any, any , any and any, the solution y of (1.1) satisfieswherewithandgiven by Lemma4.2.
The proof of Lemma 4.1 relies on the following estimates for Eq. (1.1):
There existssuch that for any, any , any and any, the solution y of (1.1) satisfieswhere
Furthermore, if , then where
To start the proof, we perform the change of variables
Then, satisfies system (1.1) for
with
For y we will prove the following estimates:
and, if ,
from where inequalities (4.3) and (4.4) are easily deduced by going back to the original variable. Here and in the following, C denotes a positive constant independent of all parameters.
For the sake of clearness we will denote y, t and x instead of , and until the end of the proof of Lemma 4.2.
The proof of (4.5) and (4.6) is divided in several steps:
A previous computation.
Let us first perform a previous computation which will be useful in the proof. Let be a positive increasing function such that . After some integrations by parts we obtain, for all ,
Estimate of .
In this paragraph we lift the boundary and initial conditions. Let
Then, z satisfies
where
We multiply the equation in (4.9) by and integrate in space. After integration by parts we have
To gain more derivatives of z and estimate the boundary term on , we use again (4.7) with and (). Notice that from (4.9) and , we have that . Thus:
Again from (4.9) and , we have the relation . Using it in this last inequality, we obtain
Using (4.20) to estimate the last term and applying Gronwall’s Lemma, we get
Finally, using the equation in (1.1), we obtain from (4.20) and (4.24):
From here, we obtain (4.6) as in the end of Step 2. □
As in the proof of Corollary 2.2, we write , for some , and , where p and q are the solutions of (2.9) and (2.10), respectively. From Lemma 4.2, using (4.3) for p and (4.4) for q, we have from (2.8) that
for every such that (2.6) is satisfied. Therefore, from (2.6)–(2.7) we obtain (4.1). □
Finally, let us state another technical result whose proof is given in Appendix B.
Let. For any , there exists such thatwhere
Auxiliary problem and conclusion
Now, we introduce the following auxiliary control problem:
For this problem, we define the cost associated to the null controllability as follows:
We find now a lower bound of .
Letand. Then
We remark that is the smallest constant such that
for any , where φ is the solution of
Now, let
It is straightforward to check that ψ satisfies the partial differential equation and the boundary conditions in (4.30). Furthermore, since belongs to , we have
On the other hand,
Consequently, from (4.29) we deduce that (4.28) holds. □
We argue by contradiction, i.e., we suppose that for any , there exists such that the solution of (1.1) satisfies and
Let where . We choose to be the extension of given by Lemma 4.3. In particular, we have
Observe that y solves (4.26) with and . Let us take . Then, using Lemma 4.1 we have
From (4.31) and (4.32):
which implies that
This and (4.28) show that , which is a contradiction. □
Proof of Theorem 1.5
The proof of Theorem 1.5 relies on finding a particular solution of (3.2) such that decays exponentially as and behaves like a constant in ε. The proof we perform here is inspired by [10, Theorem 1.4]. The main difference with respect to [10] is that the boundary condition on is not homogeneous.
First, notice that we have the dissipation estimate:
Now, we choose such that
and a nonnegative function such that
Let be the solution of (5.1) associated to as initial condition. We will prove that
and
where depends on at most polynomially.
Let us explain how (5.5) and (5.6) allow us to conclude Theorem 1.5. Let be a control which drives the solution y of (1.1) from to 0 (we know such a v exists by Theorem 1.1 and [12]). We multiply (1.1) by and integrate by parts to get
Setting and using (5.4), (5.5) and (5.6) in this last inequality we obtain (1.8).
Let us define in Q. It follows from (5.3) and (5.4) that for all .
We multiply (5.1) by θ and after integration by parts we obtain
On the other hand, from the definition of θ, (5.3) and (5.4), it is easy to see that
and
Using these elements in (5.7), together with Young’s inequality, we obtain (5.5) for ε small enough depending on T and R. □
The objective is to prove that
from where (5.6) will readily follow. This is done by proving the estimate
and then applying an internal regularity result proved in [10] to conclude.
We consider a cut-off function such that
We set and multiply (5.1) by , where is to be chosen later on. We perform several integrations by parts, but observe that from (5.3) and (5.10), we have that for all , so there are no boundary terms. We get, after neglecting the positive terms,
where is a polynomial function of degree 2 in r and we have used (5.10) to restrict the limits in the integral. Multiplying by and using that β is decreasing, we have
By (5.2), we obtain
Integrating in , we get:
where we have used the fact that for all .
Now, notice that for all thanks to (5.3), so we have
and thus
Again from (5.3), we obtain
We finish the proof of (5.9) by choosing such that it minimises the expression inside the exponential, that is,
To prove (5.8), we will use the following lemma, which corresponds to Proposition 3.3 in [10].
Letand. Consider a solution w of for someand. Then,, with the estimate for some constantdepending at most polynomially inand.
Let and apply Lemma 5.1 with and instead of and , respectively. Notice that with this setting we have and . Thus, from (5.12) we have
Now, we estimate the term in the right-hand side in a slightly larger interval. To do this, we multiply (5.1) by and integrate in . We obtain
Since in , we get by integrating between 0 and T
Combining this with (5.13) and (5.9), we obtain (5.8). □
Great part of this work has been accomplished while the first author was a PhD student at Laboratoire Jacques-Louis Lions, Université Pierre et Marie Curie, Paris, France. He would like to thank them for their support and hospitality during those years. This work has been partially supported by Fondecyt 3150089 (N. Carreño).
Proof of Proposition 3.1
We now follow the steps of [9] and [10]. Let . Using Eq. (3.4), we get
where we have denoted
and
Notice that we have the following boundary values for ψ:
Taking the -norm we have
In the following, our efforts will be devoted to computing the double product in the previous equation. Let us denote by the jth term of .
Computing.
For the first term, we integrate by parts twice in space:
Using the properties (3.12), (3.13), (A.2) and (A.3), we obtain
We integrate by parts again in space in the second term, and using the boundary values for and :
Here, we have used the boundary values (A.2), (A.3), the properties in (3.12) and Young’s inequality.
For the third term, we proceed in a similar manner:
Putting together these computations, we obtain
Computing.
For the first term, we integrate by parts in time:
The second terms gives:
In the third term, we integrate by parts first in space and then in time. We obtain
where we have used the boundary value (A.2).
Putting together this inequalities, we have
Computing.
We integrate by parts in space and the properties (3.12)–(3.13) to treat the first term:
The second term simply gives:
As for the third term, we integrate by parts in space once and use the boundary value (A.2). We obtain:
Putting together these estimates, we get:
Computing.
For the first term, we integrate by parts in space:
The second term gives directly:
The third and final term gives, after integration by parts:
Putting together these expressions, we obtain:
The entire product.
Adding inequalities (A.5)–(A.8), we find four positive terms, namely:
In the following, we explain how to estimate the nonpositive integrals coming from the addition of (A.5)–(A.8) in terms of .
Let us start with the terms concerning in Q. We can easily check that they can all be bounded by
which by taking can be absorbed by .
The integral of , can be bounded by
Taking , this term can be absorbed by .
Furthermore, taking shows that all the integrals concerning can be estimated by
Finally, let us treat the terms containing in Q. They can be estimated by
Now, similarly as the previous steps, integration by parts in space shows that:
Notice that here we have used (A.2) and Young’s inequality. By taking , the first two integrals can be absorbed by and , respectively, and the last one can be estimated by (A.9).
Finally, all these estimations give
for every .
Coming back to (A.4), and together with the fact that
we obtain
for every .
Coming back to the original variable.
Let us now go back to the original variable ϕ. First, we point out that the same computations made in (A.10), show that
as long as . This means that we can add this term to the left-hand side of (A.11), and together with , we have directly from (A.11) that
for every .
Now, from , we find that
and taking the -norm, we see that we can add
to the left-hand side of (A.12). Similarly, from
we can add
to the left-hand side of (A.12) if . Finally, using (A.1), we can add to the left-hand side of (A.12) the respective boundary integrals and obtain
for every . The proof of Proposition 3.1 is complete.
Proof of Lemma 4.3
We define by:
where is given by
It is easy to check that and
Let us estimate the term concerning p. Notice that
Taking the square and using Cauchy–Schwarz inequality we get
where we have used that in the last inequality.
Combining this with (B.1) and using Young’s inequality, we obtain the desired result.
References
1.
E.Cerpa and E.Crépeau, Boundary controllability for the nonlinear Korteweg–de Vries equation on any critical domain, Ann. Inst. H. Poincaré Anal. Non Linéaire26(2) (2009), 457–475.
2.
E.Cerpa, I.Rivas and B.-Y.Zhang, Boundary controllability of the Korteweg–de Vries equation on a bounded domain, SIAM J. Control Optim.51(4) (2013), 2976–3010.
3.
T.Colin and J.-M.Ghidaglia, Un problème mixte pour l’équation de Korteweg–de Vries sur un intervalle borné, C. R. Acad. Sci. Paris Sér. I Math.324(5) (1997), 599–603.
4.
T.Colin and J.-M.Ghidaglia, An initial-boundary value problem for the Korteweg–de Vries equation posed on a finite interval, Adv. Differential Equations6(12) (2001), 1463–1492.
5.
J.-M.Coron, Control and Nonlinearity, Mathematical Surveys and Monographs, Vol. 136, American Mathematical Society, Providence, RI, 2007.
6.
J.-M.Coron and S.Guerrero, Singular optimal control: A linear 1-D parabolic–hyperbolic example, Asymptot. Anal.44(3–4) (2005), 237–257.
7.
L.Evans, Partial Differential Equations, Graduate Studies in Mathematics, Vol. 19, American Mathematical Society, Providence, RI, 1998.
8.
A.V.Fursikov and O.Yu.Imanuvilov, Controllability of Evolution Equations, Lecture Notes Series, Vol. 34, Global Analysis Research Center, Research Institute of Mathematics, Seoul National Univ., Seoul, 1996.
9.
O.Glass and S.Guerrero, Some exact controllability results for the linear KdV equation and uniform controllability in the zero-dispersion limit, Asymptot. Anal.60(1-2) (2008), 61–100.
10.
O.Glass and S.Guerrero, Uniform controllability of a transport equation in zero diffusion–dispersion limit, Math. Models Methods Appl. Sci.19(9) (2009), 1567–1601.
11.
S.Guerrero and G.Lebeau, Singular optimal control for a transport–diffusion equation, Comm. Partial Differential Equations32(10–12) (2007), 1813–1836.
12.
J.-P.Guilleron, Null controllability of a linear KdV equation on an interval with special boundary conditions, Math. Control Signals Syst.26(3) (2014), 375–401.
13.
E.Kramer, I.Rivas and B.-Y.Zhang, Well-posedness of a class of non-homogeneous boundary value problems of the Korteweg–de Vries equation on a finite domain, ESAIM Control Optim. Calc. Var.19(2) (2013), 358–384.
14.
J.-L.Lions, Exact controllability, stabilization and perturbations for distributed systems, SIAM Rev.30(1) (1988), 1–68.
15.
L.Rosier, Exact boundary controllability for the Korteweg–de Vries equation on a bounded domain, ESAIM Control Optim. Calc. Var.2 (1997), 33–55.
16.
L.Rosier, Control of the surface of a fluid by a wavemaker, ESAIM Control Optim. Calc. Var.10(3) (2004), 346–380.