A higher-order nonlinear Boussinesq system with a time-dependent boundary delay is considered. Sufficient conditions are presented to ensure the well-posedness of the problem by utilizing Kato’s variable norm technique and the Fixed Point Theorem. More significantly, the energy decay for the linearized problem is demonstrated using the energy method.
The Boussinesq system comprises a set of nonlinear partial differential equations that model wave dynamics in fluids with small amplitude and long wavelengths. Originally formulated by the French mathematician Joseph Boussinesq in the 19th century to describe shallow water waves (Boussinesq, 1871), since then, the system has been recognized as a model for various physical phenomena, including ocean currents, atmospheric circulation, and heat transfer in fluids. Consequently, the Boussinesq system remains an essential tool in numerous fluid dynamics, with broad applications in fields such as meteorology, oceanography, and engineering.
In more recent studies, Bona et al. (2002, 2004) introduced a four-parameter family of Boussinesq systems to describe the motion of small-amplitude long waves on the surface of an ideal fluid under gravity, particularly in scenarios where the motion is predominantly two-dimensional. In particular, Bona et al. (2002, 2004) investigated the following system:
In this context, represents the elevation of the fluid surface from its equilibrium position, while denotes the horizontal velocity of the flow at a height , where is the undisturbed depth of the fluid and is a constant within the interval . The variables and correspond to space and time, respectively, and the physical parameters must satisfy the following relationships:
Stabilization results for the higher-order system (1.1) on the periodic domain were established in Bautista and Pazoto (2020b) under the conditions , with general damping applied to each equation. Furthermore, the local exact controllability of system (1.1) was investigated in Bautista et al. (2021), where the control is localized within the interior of the domain and influences only one equation.
Negative controllability results are explored in Bautista and Pazoto (2020a) and Sierra Fonseca and Pazoto (2022) when the third- and fifth-order Korteweg–de Vries (KdV) terms are removed from the system mentioned above, that is, (1.1) with . In this case, the system consists of two coupled Benjamin–Bona–Mahony-type equations. The authors demonstrated that the linear model is approximately controllable but not spectrally controllable. This implies that although any state can be brought arbitrarily close to another, no finite linear combination of eigenfunctions, other than zero, can be driven to zero.
Let us now consider and make a scaling argument to obtain the fifth-order Boussinesq system
In the above system, we note that . Thus, we consider the following case:
It is important to highlight that, to the best of our knowledge, there are no existing results that combine a damping mechanism with a boundary time-varying delay to achieve stabilization of the higher-order Boussinesq system associated with (1.2). This gap forms the main motivation for the present work. It is also worth emphasizing the practical prevalence of time delays in control systems, a phenomenon that is virtually unavoidable due to factors such as the lag between sensors, actuators, and data processing. In light of this, significant efforts have been made to mitigate or eliminate the effects of constant time delays. For example, stability results under smallness conditions on the domain length and the initial data have been established in Baudouin et al. (2019) and Valein (2022) for the KdV equation and in Capistrano-Filho et al. (2023) for the Kawahara equation. In the case of time-dependent delays, similar results have been obtained in Parada et al. (2023) for the KdV equation and Capistrano-Filho et al. (2024) for the coupled KdV–KdV system.
Notations and Main Results
This article is concerned with the following system:
where the parameters verify (1.3). Moreover, we assume that there exist positive constants , , and such that the time-dependent delay function satisfies the following standard conditions:
Finally, the feedback gains and must obey the following constraint:
The condition (1.6) ensures the dissipation of system (1.4). It is worth mentioning that a similar condition is used for other types of delayed dispersive systems (see, for instance, Capistrano-Filho et al., 2024; Parada et al., 2023). Furthermore, recalling that and , one can readily check that (1.6) is fulfilled if, for instance, the feedback gains and satisfy
Next, let and the state space
equipped with the inner product
for any . Moreover, we shall consider the space
whose norm is
To present our first result, let us introduce the following space:
The first result of this manuscript ensures the local well-posedness of system (1.4).
Let and suppose the parameters verify (1.3). Then, there exists such that, for every satisfying
the system (1.4) admits a unique solution . Moreover
for some positive constant, .
Our second result is closely related to the total energy associated with system (1.4) that is defined in by
Indeed, the second result of the article guarantees that the energy associated with the following system:
decays exponentially, even in the presence of delay, and provides an estimate of the decay rate. The result is expressed as follows:
Let the parameters verify (1.3) and . Suppose also that the time-dependent delay function satisfies (1.5). Then, there exist two positive constants
and
such that the energy given by (1.8) associated with system (1.9) satisfies
Here and are two positive constants small enough to be well-chosen.
The main contribution of this work is to establish the local well-posedness of system (1.2) and prove the exponential stability of system (1.9). These results extend and refine those obtained in Capistrano-Filho et al. (2024) and Parada et al. (2023) in several significant directions. More specifically, unlike Capistrano-Filho et al. (2024), where the nonlinear coupling appears only in one equation through the term , our model introduces an additional nonlinear coupling term of the form . Furthermore, while Capistrano-Filho et al. (2024) feature a single uncoupled nonlinear term, our system includes four additional uncoupled nonlinear terms of higher order, making the analysis more intricate and requiring careful handling of the extra terms in the computations. In particular, in our case, the higher-order spatial derivatives (of order three and five) appear with positive signs, leading to conflicts between the and -norm terms during integration by parts. Compared to Parada et al. (2023), the situation is even more complex: the problem in Parada et al. (2023) involves a single uncoupled equation with only one nonlinearity and a highest-order derivative of three. Finally, and importantly, in contrast to Capistrano-Filho et al. (2024) and Parada et al. (2023), we employ the transposition method (Capistrano-Filho et al., 2019) to address the well-posedness of our system.
Outline
The structure of the paper is as follows. In Section 2, we establish the well-posedness of the nonlinear problem (1.4), namely, we show Theorem 1.1 starting with an analysis of the linear system (1.9) using the variable norm technique of Kato, followed by the application of the Fixed Point Theorem to prove well-posedness of the full nonlinear problem. Section 3 focuses on the stability result presented in Theorem 1.2, along with a discussion of the optimal decay rate. Finally, we conclude the paper with further remarks in Section 4.
Well-Posedness Results
In the sequel, we will assume and in (1.4) as well as in (1.9). We will first examine the well-posedness of the linear system (1.9) and subsequently analyze the properties of the nonlinear problem (1.4) in suitable spaces.
Linear Problem
Consider the following linear Cauchy problem:
where is densely defined, and is independent of time , that is, for all The next theorem ensures the existence and uniqueness of the Cauchy problem (2.1).
Now, we pick up and consider the time-dependent operator
given by
with a domain defined by
This allows us to write problem (1.9) in the abstract form (2.1) by using (2.2)–(2.4). Additionally, it is noteworthy that is independent of time since .
Now, taking the triplet , with , for some fixed and , we can state and prove the well-posedness result of (2.1) related to .
Let the parameters verify (1.3). Assume that and are real constants such that (1.6) holds. Taking , there exists a unique solution to (2.1) whose operator is defined by (2.3)–(2.4). Moreover, if , then
The result will be proved in a standard way (see, for instance, Nicaise et al., 2009). First, it is not difficult to see that is a dense subset of and , for all . Thus, the requirement (1) of Theorem 2.1 is fulfilled.
Concerning the condition (2) of Theorem 2.1, let us note that simple integrations by parts together with the boundary conditions yield
For all , the operator is maximal, or equivalently, we have that is surjective, for some .
In fact, let us fix . Given , we seek a solution of the equation , that is,
One can readily verify that is given by
Thereby, , in which
and
Combining the latter with (2.5), it follows that and are solutions of the system
and satisfy the boundary conditions
Now, let be a function such that and . Next, we define a function and let . This, together with (2.6), implies that and satisfy
as well as the boundary conditions
Let us mention that for the sake of simplicity, we still use after translation. Then, we can verify that (see, for instance, Capistrano-Filho et al., 2024). Thus, thanks to (1.6), we deduce that . Consequently, showing the Claim 1 is equivalent to proving that is surjective, where is given by
with a dense domain
Now, observe that adjoint of , denoted by , is defined by
with
Since
and
we can claim that the operators and are dissipative. Therefore, the desired result follows from the Lummer–Phillips theorem (see, e.g., Pazy, 1983). This shows the Claim 1. Consequently, generates a strongly semigroup on and is a stable family of generators in , whose stability constant is independent of . Thus, the condition (2) of Theorem 2.1 is satisfied.
Lastly, since for all , we reach that
is bounded on for all and
Moreover, the coefficient of is bounded on , and the regularity (3) of Theorem 2.1 is satisfied.
To sum up, we verified the assumptions of Theorem 2.1 and hence for each , the Cauchy problem
has a unique solution and . Thus, the solution of (2.1) is explicitly given by .
We also have the following result.
Let the parameters verify (1.3). Suppose and are real constants such that (1.6) holds. Then, for any mild solution of (2.1), the energy defined by (1.8) is nonincreasing and
We now proceed to prove the Kato smoothing property, along with several a priori estimates. These results are crucial for establishing the well-posedness of the system (1.4). In the following, represents the two-parameter semigroup of contractions associated with the operator . We are now prepared to state the following result:
Let the parameters verify (1.3) and and are real constant such that (1.6) holds. Then, the following estimate holds:
Furthermore, for every initial condition , we have that
On the other hand, for the initial datum, we have the following estimates:
and
Finally, for , the Kato smoothing effect is verified
and the map
is well-defined and continuous.
Using (2.7) and the fact that is a symmetric negative definite matrix, we deduce the existence of a positive constant , such that
Now, we show the inequality (2.12) provided that . Initially, multiplying the first equation of (1.9) by and the second one by . Next, adding the results, then integrating by parts over and invoking (2.8) and (2.9), we obtain
for some positive constant . Since , from Poincaré inequality, there exists , such that
In this subsection, we show the well-posedness of the nonlinear problem (1.4) by using the approach of Capistrano-Filho et al. (2019), where the solutions are obtained via the transposition method and the existence and uniqueness by using the Riesz-representation theorem.
To prove the well-posedness result for system (1.4), we consider the nonhomogeneous system
where the parameters verify (1.3). Remember the definition of given by (1.7), and also consider the following set:
We define a solution by transposition as follows, see Lions and Magenes (1968, 1986), to justify the choice of formula (2.17) below:
(Solution by transposition)
Let , and
A solution of problem (2.16) is a function such that, for all and the following identity holds
where the pair is the solution of
Thanks to Capistrano-Filho et al. (2019, Corollary 2.5 and Proposition 2.6), the following well-posedness result for system (2.18) is established:
For all system (2.18) admits a unique solution which satisfies
Additionally, we have that
The following result gives us the existence and uniqueness of the solution for system (2.16).
Let and . There exists a unique solution of system (2.16). Moreover, there exists a positive constant , such that
for all .
Let us define as the linear functional given by the right-hand side of (2.17), that is
We infer from (2.19), (2.20), and the Cauchy–Schwarz inequality that
and we obtain that Thus, from the Riesz-representation theorem, there exists one and only one such that
and we obtain the uniqueness of the solution to problem (2.16). Now, in order to prove the estimate (2.21), we define the map as
Now, we pass to show the well-posedness of the nonhomogeneous feedback linear system associated with (2.16)
Let . Then, for every in and in , there exists a unique solution of system (2.16) such that with , where and belong to . Moreover, for some positive constant , we have
for all
Note that if from the trace theorems, it follows that
We claim that: there exists a positive constant such that
Indeed, note that
By using the conditions in (1.5), we deduce the existence of some positive constant such that
giving the claim.
Now, let to be determined later. For each consider the map
where is the solution of (2.16) with By Lemma 2.7 and (2.23), the linear operator is well defined. Furthermore, there exists a positive constant such that
Hence, is a contraction, and by Banach Fixed Point Theorem, we obtain a unique such that and
Since is independent of the standard continuation extension argument yields that the solution belongs to and the proof ends.
The first main result of the article ensures the existence of local solutions to (1.4) and is proved below.
Proof of Theorem 1.1.
Let and where will be determined later. We know from Capistrano-Filho et al. (2019) that for there exists a positive constant such that the following inequalities hold true:
and
Thus, the nonlinearities
belong to and
Taking this into consideration, we define the following map:
we obtain that Finally, following the same argument as done in Lemma 2.8, we can conclude that is a contraction in , then, the Banach Fixed Point Theorem guarantees the existence of a unique such that and
achieving the proof.
In (1.5), the time-dependent delay is assumed to be positive for all . This requirement is relaxed in Nicaise et al. (2011) since is allowed to degenerate. Notwithstanding, the problem in Nicaise et al. (2011) is linear and hence simpler than ours. The key idea of the proof in Nicaise et al. (2011) is to consider a new delay defined by
where , for some . Therefore, satisfy (1.5) and hence problem (2.1) has a unique solution . The whole task is to tend to under more regularity on the solution. We have tried to adopt this approach, but we faced difficulties because of the nonlinearities in our problem.
Long-Time Behavior of Solutions
In this section, we are in a position to prove the second main result of our work. First, we demonstrate that the energy associated with (1.9) is exponentially stable. Moreover, we establish that the solutions decay at an optimal rate.
Proof of Theorem 1.2
Recall that Theorem 2.2 (see also Proposition 2.3) guarantees the a priori estimate for the linear system (2.1) whose operator is defined by (2.3)–(2.4). Therefore, the solutions of system (1.9) are globally well-posed. Whereupon, we can treat the exponential stability for this system.
To proceed, consider the following Lyapunov functional
where will be chosen later. Here, is the total energy given by (1.8), while
and
Observe that,
by assuming and .
On the other hand, using system (1.9) and the boundary conditions, we get that
In addition, from (2.2) and by integration by parts, we deduce that
The objective is to show that . To do so, let us deal with each term in (3.4), for .
Estimate of : Since the matrix (see (1.6)) is negative definite, it follows from the continuity of the trace and determinant functions that one can choose sufficiently small so that the new matrix is also negative definite. Thus,
Estimate of : Observe that using Poincaré inequality, we get
for and fulfilling (1.10) and (1.11), respectively. This achieves the proof of the theorem. □
Decay Rate: An Optimal Result
We can optimize the value of in Theorem 1.2 to obtain the best decay rate for the linear system (1.9) in the following way:
If the constant given in Theorem 1.2 is chosen as follows:
then has the largest possible value.
Define the functions and
by
and
respectively. On the other hand, let us consider . Thus, we have the following claims.
The function (respectively, ) is increasing (respectively, decreasing) in the interval
A simple computation shows that
and hence for
Furthermore, one can rewrite as follows:
and thus
This ascertains the Claim 2.
There exists only one point , satisfying (3.5) such that .
Indeed, since
and
the existence of this point is a direct consequence of the mean value theorem, applied to the function . The uniqueness follows from the fact that the function is increasing in this interval, and Claim 3 holds.
Lastly, thanks to the Claims 2 and 3, the maximum value of the function is obtained when satisfies (3.5), where , and the proof of Proposition 3.1 is achieved.
Conclusion
This paper establishes the existence and uniqueness of a solution for a higher-order nonlinear Boussinesq system in a bounded domain, even when a time-dependent delay is present in one of the boundary conditions. Additionally, we prove that solutions to the linearized problem are exponentially stable, both results being obtained under certain conditions related to the system’s parameters and the delay. These findings extend the results of the second and third authors in Capistrano-Filho et al. (2024) for a higher-order dispersive system. Further comments on our results are provided below.
It is worth mentioning that the solutions of system (1.4) obtained in Theorem 1.1 are local. Proving the global existence of solutions remains a challenge due to the absence of an a priori estimate. Specifically, it is challenging to tackle this problem within the energy space for the nonlinear system that includes a delay term.
Observe that the restriction in Theorem 1.2 arises from the Kato smoothing effect, which does not occur in the lower-order Boussinesq system (see, e.g., Capistrano-Filho et al., 2024). This difference is because in system (1.4), we have spatial derivatives of order three and five, both with positive signs. Thus, after performing some integration by parts, the left-hand side of (2.14) contains the -norm with a negative sign and the -norm with a positive sign. To recover the -norm, the Poincaré inequality must be applied, which imposes this restriction on the size of .
A version of the higher-order Boussinesq system was proposed by Olver (1984, equations (4.7) and (4.8), p. 283) and is given by:
Through scaling, we arrive at the following system:
where and . System (4.1) was studied in Capistrano-Filho et al. (2019). Using the same boundary conditions as in problem (1.4), we believe that similar results proved in our work can be obtained for system (4.1) without the restriction over since the sign of the third derivatives in (4.1) is negative instead of positive unlike our case (see system (1.4)).
It is important to point out that system (1.4) is shown to be locally well-posed, and hence we are unable to establish any exponential stability for the nonlinear problem. One interesting research avenue is to show the stability for the nonlinear problem.
Footnotes
Acknowledgments
The authors are grateful to the associate editor and the referees for the careful reading of this paper and their valuable suggestions and comments.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Roberto de A Capistrano-Filho was partially supported by CAPES/COFECUB grant number 88887.879175/2023-00, CNPq grant numbers 421573/2023-6 and 307808/2021-1, and Propesqi - QUALIS A (UFPE). George J Bautista was supported by Universidad Tecnológica de los Andes, Abancay-Peru. Oscar Sierra Fonseca was supported by FAPERJ (Rio de Janeiro, Brazil) grant number SEI 260003/000175/2024 and by Escola de Matemática Aplicada, Fundação Getúlio Vargas (Rio de Janeiro - Brazil).
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Data Availability Statement
It does not apply to this article as no new data were created or analyzed in this study.
Notes
References
1.
BaudouinL.CrépeauE.ValeinJ. (2019). Two approaches for the stabilization of nonlinear KdV equation with boundary time-delay feedback. IEEE Transactions on Automatic Control, 64(4), 1403–1414.
2.
BautistaG. J.MicuS.PazotoA. F. (2021). On the controllability of a model system for long waves in nonlinear dispersive media. Nonlinearity, 34, 989–1013.
3.
BautistaG. J.PazotoA. F. (2020a). On the controllability of a Boussinesq system for two-way propagation of dispersive waves. Journal of Evolution Equations, 20, 607–630.
4.
BautistaG. J.PazotoA. F. (2020b). Decay of solutions for a dissipative higher-order Boussinesq system on a periodic domain. Communications on Pure and Applied Analysis, 19, 747–769.
5.
BonaJ. L.ChenM.SautJ.-C. (2002). Boussinesq equations and other systems for small-amplitude long waves in nonlinear dispersive media. I. Derivation and linear theory. Journal of Nonlinear Science, 12, 283–318.
6.
BonaJ. L.ChenM.SautJ.-C. (2004). Boussinesq equations and other systems for small-amplitude long waves in nonlinear dispersive media. II. The nonlinear theory. Nonlinearity, 17, 925–952.
7.
BoussinesqJ. V. (1871). Théorie générale des mouvements qui sont propagés dans un canal rectangulaire horizontal. Comptes Rendus de l’Acadèmie des Sciences, Paris, 72, 755–759.
8.
Capistrano-FilhoR. A.ChentoufB.Gonzalez MartinezV. H.MuñozJ. R. (2024). On the boundary stabilization of the KdV–KdV system with time-dependent delay. Nonlinear Analysis: Real World Applications, 79, 104122.
9.
Capistrano-FilhoR. A.ChentoufB.SousaL. S.Gonzalez MartinezV. H. (2023). Two stability results for the Kawahara equation with a time-delayed boundary control. Zeitschrift für Angewandte Mathematik und Physik, 74, 16.
10.
Capistrano-FilhoR. A.GallegoF. A. (2018). Asymptotic behavior of Boussinesq system of KdV–KdV type. Journal of Differential Equations, 6, 2341–2374.
11.
Capistrano-FilhoR. A.GallegoF. A.PazotoA. F. (2019). On the well-posedness and large-time behavior of higher order Boussinesq system. Nonlinearity, 5, 1852–1881.
12.
KatoT. (1970). Linear evolution equations of “hyperbolic” type. Journal of the Faculty of Science, University of Tokyo. Section I: Mathematics, Astronomy, Physics, Chemistry, 17, 241–258.
13.
LionsJ.-L.MagenesE. (1968). Problémes aux limites non Homogénes et applications, vol. 1. Dunod.
14.
LionsJ.-L.MagenesE. (1968). Problémes aux limites non Homogénes et applications, vol. 2. Dunod.
15.
NicaiseS.PignottiC. (2006). Stability and instability results of the wave equation with a delay term in the boundary or internal feedbacks. SIAM Journal on Control and Optimization, 5, 1561–1585.
16.
NicaiseS.PignottiC.ValeinJ. (2011). Exponential stability of the wave equations with boundary time-varying delays. Discrete and Continuous Dynamical Systems — Series S, 4, 693–722.
17.
NicaiseS.ValeinJ.FridmanE. (2009). Stability of the heat and of the wave equations with boundary time-varying delays. Discrete Continuous Dynamical Systems-S, 3, 559–581.
18.
OlverP. J. (1984). Hamiltonian and non-Hamiltonian models for water waves. Lecture Notes in Physics, 195, 273–96.
19.
ParadaH.TimimounC.ValeinJ. (2023). Stability results for the KdV equation with time-varying delay. Systems & Control Letters, 177, 105547.
20.
PazyA. (1983). Semigroups of linear operators and applications to partial differential equations, volume 44 of applied mathematical sciences. Springer-Verlag.
21.
Sierra FonsecaO.PazotoA. F. (2022). On the lack of controllability of a higher-order regularized long-wave system. Communications in Mathematical Sciences, 20, 1243–1278.
22.
ValeinJ. (2022). On the asymptotic stability of the Korteweg–de Vries equation with time-delayed internal feedback. Mathematical Control & Related Fields, 12(3), 667–694.
23.
XuG. Q.YungS. P.LiL. K. (2006). Stabilization of wave systems with input delay in the boundary control. ESAIM: Control, Optimisation and Calculus of Variations , 12, 770–785.