The Wave Equation with Dynamic Wentzell Boundary Condition in Polygonal and Polyhedral Domains: Obse...

Tawfik Masrour

  Open Access OPEN ACCESS  Peer Reviewed PEER-REVIEWED

The Wave Equation with Dynamic Wentzell Boundary Condition in Polygonal and Polyhedral Domains: Observation and Exact Controllability

Tawfik Masrour

Research Laboratory (M2.I), Mathematical Modeling for Analysis and Decision Making Research team (M2APD), Moulay Ismail University ENSAM, Meknes, Morocco


We study in this article the boundary observability and the exact controllability for a problem of transmission. The system is governed by the wave equation with Wentzell dynamic artificial condition on the boundary. The geometrical domains considered are polyhedrons or polygons.

At a glance: Figures

Cite this article:

  • Masrour, Tawfik. "The Wave Equation with Dynamic Wentzell Boundary Condition in Polygonal and Polyhedral Domains: Observation and Exact Controllability." International Journal of Partial Differential Equations and Applications 2.1 (2014): 13-22.
  • Masrour, T. (2014). The Wave Equation with Dynamic Wentzell Boundary Condition in Polygonal and Polyhedral Domains: Observation and Exact Controllability. International Journal of Partial Differential Equations and Applications, 2(1), 13-22.
  • Masrour, Tawfik. "The Wave Equation with Dynamic Wentzell Boundary Condition in Polygonal and Polyhedral Domains: Observation and Exact Controllability." International Journal of Partial Differential Equations and Applications 2, no. 1 (2014): 13-22.

Import into BibTeX Import into EndNote Import into RefMan Import into RefWorks

1. Introduction

We consider the problem of observation and control of a system of transmission waves. The system is governed by a three-dimensional (respectively two-dimensional) D’ Alembert wave equation in a bounded domain which comprises corners but which is not fissured. On the boundary we have inside an artificial dynamic Wentzell boundary condition coupled with the internal equation by the normal derivative. Let us start by giving a few motivations of interest in the study of observation, control, stabilization and inverse problems with artificial conditions: The condition at the boundary can be seen like a contribution of energy (kinetic energy and potential energy): through entering flow and which is due to the two-dimensional wave propagation (respectively one-dimensional) on the boundary. The dynamic Wentzell condition can, then, being written in a simple way in the case as of acoustic waves or of vibration longitudinal as follows: .

The coefficients α and β can, in certain cases, be some appropriate operators, they can also be simply constants containing information on the physical characteristics of the problem as well as information about geometrical characteristics of the domain. Indeed the boundary conditions of the type Wentzell involved in many physical problems: electromagnetic, acoustics, structural mechanics or in heat problems. One of the most important reasons to use such conditions is the presence of thin layers. Thin layer, have in general, very small dimensions compared with characteristic values of the main problem. The numerical calculation of systems with thin layers presents many problems at several levels, and the discretization with a size scale of the thin layer becomes difficult and leads to expensive calculations and sometimes very sketchy. A very interesting alternative, then, is to model the problem with a system that replaces the presence of thin layers by a boundary on which we impose an "artificial" conditions. There is a whole family of artificial conditions sometimes coupling the Dirichlet, Neumann and operators boundaries may be the same as the operator governing equations of the system whose dimensions are dominant. We can cite for example the case of anechoic chambers are used to analyze the waves emitted by telecommunication devices: these are rooms whose walls, floors and ceilings are covered with an absorbent so that reflections on the material walls do not interfere with the experiment. The thickness of the absorbent zone is small compared to the characteristic dimensions of the work-piece and can therefore be modeled as a thin layer. We can also cite the example of radar stealth: a special paint covers the plane and it absorbs the waves emitted by the radar, which allows the aircraft to remain invisible. The relationship between the thickness of the paint layer and the characteristic dimensions of the device is also very low.

The Wentzell condition was first introduced by [1]. For the physical meaning of the Wentzell condition and artificial conditions, as well as for the study of some mathematical and numerical aspects of such problems, the interested reader may consult for example: [2-8][2] and references in these Articles, and the list is not exhaustive.

Before stating our main result, it may be interesting to give brief state of the art methods in observation and control to locate this work in context and to make a link with related work. The study of many problems of controllability or stabilization or even inverse problem reduces to finding a priori estimates on the solutions of adjoint problems [9, 10]. These estimates are clearly relevant to the notion of stable observation, and it constitute an important ingredient in the Hilbert Uniqueness Method (H.U.M) of J.L. Lions [11]. In the case of regular open we often have the choice between the use of so-called global or local multiplier techniques (Rellich) or technical finer microlocal analysis: by classical approach consisting in propagation analysis of singularities [12, 13, 14, 15, 16], or by using the microlocal defect measures (or H measures) consisting in analysis propagation analysis of regularity (or compactness) [17, 18]. Both techniques (global or microlocal) give relatively satisfactory results, but with one notable difference. Indeed, global techniques (multipliers) are relatively simple but the results are not always optimal especially for the region of observation or control and the time of observation or action. Classical microlocal techniques are by far more difficult against but give finer and often optimal results. The introduction of microlocal defect measures in some cases considerably simplifies demonstrations while keeping the optimality and the accuracy of the results [19]. In the presence of singularities due to the geometry of the domains the questions of observation, control, Stabilization and inverse problems are much more difficult and there are strict restrictions to extend both technical and results[20, 21]. For example it is not easy to work with the tools of microlocal analysis in the case where open are not polyhedral convex (this is the case discussed in this article), or in the case of open only Lipchitz [22].

For some results on control or stabilization with artificial conditions we can see for example [10,23-28], and others.

Figure 1. Original Geometry with different thin layer
Figure 2. Geometry after transformation with different artificial boundary condition

2. Position of the Problem and Main Result

We note a bounded open domain of and Q the cylinder . We assume that is a polygon (or polyhedron) and we note its boundary.

We assume that there are two open subsets of the boundary and such that: . Where with , designating the union of all portions of the boundary where we impose a Dirichlet condition, and where with , designating the union of all portions of the boundary where we impose a Wentcell condition.

For each portion , , we note its two oriented vertices and the angle between et . We also note for each by the set . Let us then define as the state of the system Q and the trace on ; , satisfying:


We then act on the system through the controls as follows:


where denotes the number of vertices of .

The system data are:


The control problem is then formulated as follows:

Given , can we find for all initial data in a suitable space, a control functions , a geometric region in the boundary, such that the system satisfies:


Theorem 1.

There is a time , a region and a set acting on the system in such that the system is brought back to the idle state (7)-(8).

To prove Theorem 1 we will use Hilbert Uniqueness Method (HUM) [11] and we proceed as follows:

In Section 3 we study the adjoint problem to the system under question, and we establish results regularity. Then, in Section 4 will be given the estimated observability using the technique of Rellich multipliers and regularity results established in section 3. In Section 5 we will implement the method of Hilbert uniqueness to complete the proof of the main result. In Conclusion of this article we give some remarks concerning the controllability of the transmission system in the case of regular domains and some possible prospects.

3. Regularity of the Adjoint Problem

The adjoint system to (1)-(6) is considered and we denote the corresponding solution φ:


Data of homogeneous system are:


The following spaces are then introduced:




Where D(A) is the domain of the operator concerned, we have then the following result:

Theorem 2.

Let in , then the solution of the system (9)-(14) is in space .

The proof of the theorem 2 is divided into several steps.

3.1. Proof of Theorem 2-Step 1: “Static Problem”

We begin by examining the regularity of the following static problem:


3.1.1. Weak Formulation

We endow the space with the norm:


Where denotes the norm induced by on the space of functions vanishing on . Then we have the following result:

Lemma 1.

Given , there exists a unique solution of the following variational formulation:


Proof of lemma 1.

The proof of Lemma 1 is conventional: the bilinear form associated with the formula (22) is coercive and continuous; it then suffices to apply the Lax-Milgram theorem.

3.1.2. Interpretation

The solution u of (22) clearly satisfies the equation (17), and as then .

Now as u is in , this then allows us to give sense of the normal derivative:


Where is the dual space to defined by:


However, for any , we know (cf. [29]) that there is an extension satisfying:


We deduce, then, that the solution u of the variational formulation (22) satisfies the equation of surface waves (19) in the sense of for all .

3.1.3. Regularity Hs

We notes the angles between the edges of the boundary. We have then the following regularity result of the static solution:

Lemma 2.

Let , where , then


Proof of Lemma 2.

Since , the equation (19) then asserts that .

On the other hand we know that , we deduce by interpolation that:

Furthermore, we know that in for all . We therefore deduce that:

We are then led to solve an equation of Laplace with Dirichlet data:


The compatibility conditions being satisfied:

it is, then, sufficient to apply the Grisvard’s result [29] to complete the proof of Lemma 2. For reading paper be self-sufficient we will give an outline of the demonstration:

It is well known that for every closed neighborhood U of the set of vertices such as and , where .

If it locates, then around one of these peaks is considered and the polar coordinates centered at and such that the boundary match and the boundary with .

Singular solutions of the problem then are writing


The conclusion of the lemma is immediate.

Proposition 1.

For any non-cracked domain (i.e. ), the solution u is in fact in space and checks the system (17)-(20) in a classical sense since , and then:


Where .

3.2. Step 2. “End of the Proof of Theorem 2”

The operator A is self-adjoint and positive, we can then introduce a complete system in formed by eigenfunctions of A i.e.:


The homogeneous solution of the wave transmission system (9)-(14) can be written on this basis as follows:


Where the functions satisfy:


The symbol denotes the classic Kronecker symbol and hook denotes the scalar product in .

The operator A is injective and extends in a linear operator of H to H. This provides continuous and linear dependence on the Cauchy data:


For any .

Moreover we know that , the result of Theorem 2 is immediate, indeed, it suffices to consider the particular case .

4. Observability Problem of Homogeneous Transmission

4.1. Identity

Let , and let be a sufficiently smooth vector field, we then have the following equality:


But it is possible to define on Γ and almost everywhere vector fields constants on the portions of the boundary, and such that realizes an orthonormal basis, outside vertices. We, then, have for any sufficiently smooth function u:


Where is the component of .

If we note:


Then is the component of .

Let be the adjoint Operator of and such that:


We, then, consider the surface integral :


However, in , the tangential gradient is identically zero, the second term in the right member of the equation (42) becomes:


At this level we will do a little more attention and distinguish the case of dimension two and the case of the higher dimension .

4.1.1. Case of Dimension Two

However we have: and .

Apply Green's formula on each we obtain the following identity:


where we noted


Regarding the second integral in (43) we have:


We apply again Green's formula and we obtain:


This completes the proof in the two-dimensional case.

4.1.2. Case of Dimension 3

In this case the identity (45) remains valid, provided consider:


Concerning the second integral, we have:


Where we noted the outer normal to the face considered as a variety in .

Now we know that φ is zero on , so the last term in the identity can be rewritten as follows:


Where we extended the normal of into a vector in .

On the other hand we have:


Let us, now, consider and develop the following volume integral:


finally we obtain the identity:


Remark 1.

The identity (53) is, a priori, established only for sufficiently smooth solutions, it then extends weak solutions:

in taking into account the following remark: It is known (cf. (36)) that φ is continuously and linearly dependent on the Cauchy data and , then it follows that there is a sequence of strong solutions in

which converges to in

4.2. Boundary Observability

We will show in this section the inequality so called “Observability Inequality” or “Reverse Inequality” which is the central estimate for the implementation of the Hilbert Uniqueness Method (HUM). To do this we need some definitions.


Selecting a current point x0 of the space and . We note the part of the boundary called the shadow area respect to the source x0 i.e.

Figure 3. Illustration of Lighted Area and Shadow Area of a smooth Domain
Figure 4. Illustration of Lighted Area and Shadow Area of a Polygonal Domain

We note F the space of the initial data such that the solution φ of the system (9)-(14) satisfies:


We will show that defines a norm.

Lemma 3. "Observability"

For any solution φ of the system (9)-(14) associated with the initial data in F, the quantity defines a norm. Moreover, we have the following Observability Inequality:


Where C is a positive constant depending only on the geometry.

Proof of Lemma 3.

We will use the identity (54) with a particular choice of the vector field, ie we will replace by , where , identity then becomes:


On the other hand we know that:


Then we put (59) in (58) and obtain:


If we note

by Y, we then have the following increases:


with , and we have:


Now observe that:


And by choosing we obtain:




The second term on the right in the above inequality is bounded by , which gives the existence of a constant and a time such that:


It is now easy to see that clearly defines a norm finer than the energy norm.

In addition we have the following frame:


Otherwise if we consider the solutions θ of the following system:


where are supposed to be in the space F, and are in the space . Then θ satisfies the following proprieties:


Indeed, it suffices to write the solution , where corresponds to the homogeneous solution of the Cauchy data, and in the second member.

Inclusions above are true for by definition of space F. On the other hand we have:


This then implies that:


Which completes the proof of the announced result.

5. Transposition Method

We begin this section by the second main ingredient for the implementation of Hilbert Uniqueness Method and which is based on the results previously shown:

Lemma 4.

For all in in in in , there is a unique and solutions of:


For all and .

Proof of Lemma 4.

Based on the foregoing, the right side in the identity (80) defines a continuous linear form on , which establishes the existence and uniqueness of the solution u in .

We are now able to prove our main result:

The End of the Proof of Theorem 1.

First we solve the homogeneous problem with Cauchy data in F, there solution of the system (9)-(14). Then we know that there solution of:


For all and θ solution of the system (68)-(73).

We define the operator Λ as follows:


It is a continuous linear operator from F in F which verifies more:


Λ is an isomorphism of F in F.

Given , then there is such as:


We conclude, then, posing:

6. Conclusion

In this article, we established an observation and exact controllability result of a wave boundary problem with Wentzell condition in the presence of geometric singularity of polygonal or polyhedral types. We can notice that the extension of the results in case of smooth domains does not cause any difficulty [9]. However, it would be interesting to examine the problem using microlocal analysis, and eventually find a sharp and sufficient conditions on the control region and time control. Otherwise as we noted in [9] the case where Wentzell waves propagate with a velocity smaller than wave velocity in the interior (i.e. β < 1) is more interesting. Indeed the bicharacteristics solutions of the Hamiltonian flow on the boundary are different from those on the inside, and the habitual geometric control condition is not sufficient in this case to observe and control all the energy of waves of systems. Another interesting problem would be to examine the convergence of Wentzell control system to the Dirichlet control system (respectively Neumann) when the parameter α tends to infinity (respectively to zero) and to extend the results obtained [30] in the case of Robin-Fourier control.


[1]  A. Wentzell, «On boundary conditions for multi-dimensional diffusion processes,» Theory Probab. Appl 4, p. 164-177, 1959.
In article      CrossRef
[2]  G. R. Goldstein, «Derivation and physical interpretation of general boundary conditions,» Advances in Differential Equations, vol. 11, n 14, pp. 457-480, 2006.
In article      
[3]  I. Achdou, O. Pironneau et F. Valentin, «Effective boundary conditions for laminar flows over periodic rough boundaries,» Journal of Computational Physics , vol. 147, n 11, pp. 187-218, 1998.
In article      CrossRef
[4]  W. Arendt, G. Metafune, D. Pallara et S. Romanell, «The Laplacian with Wentzell-Robin boundary conditions on spaces of continuous functions,» Semigroup Forum Springer-Verlag., vol. 67, n 12, pp. 247-261, 2003.
In article      CrossRef
[5]  Y. AMIRAT, G. A. CHECHKIN et R. R. GADYL’SHIN, «Asymptotics of simple eigenvalues and eigenfunctions for the Laplace operator in a domain with an oscillating boundary,» Computational Mathematics and Mathematical Physics, vol. 46, n 11, pp. 97-110, 2006.
In article      
[6]  V. Bonnaillie-Noël, M. Dambrine, F. Hérau et G. Vial, «On generalized Ventcel's type boundary conditions for Laplace operator in a bounded domain,» SIAM Journal on Mathematical Analysis, vol. 42, n 12, pp. 931-945., 2010.
In article      CrossRef
[7]  M. M. A. Khemmoudj, «Exponential decay for the semilinear damped Cauchy–Ventcel problem,» Bol. Soc. Parana. Mat., p. 97-116, 2004.
In article      
[8]  K. Lemrabet, «Problème aux limites de Ventcel dans un domaine non régulier,» C. R. Acad. Sci. Paris Sér I Math. 300 (15), p. 531-534, 1985.
In article      
[9]  T. Masrour, Controlabilité et observation des systèmes distribués, problèmes et méthodes, Paris: Ecole Nationale des Ponts et Chaussées Paris, 1995.
In article      
[10]  E. Zuazua, «Propagation, observation, and control of waves approximated by finite difference methods,» SIAM review, vol. 47, n 12, pp. 197-243, 2005.
In article      CrossRef
[11]  J. Lions, Controlabilité exacte, perturbations et stabilisation de systèmes distribués tome 1, Paris: Masson, 1988, pp. [9] J.L. Lions, Controlabilité exacte, perturbations et stabilisation de systèmes distribués, tome 1, Masson, Paris, 1988.
In article      
[12]  R. B. Melrose et J. Sjöstrand, «Singularities of boundary value problems. I,» Communications on Pure and Applied Mathematics, vol. 31, n 15, pp. 593-617, 1978.
In article      CrossRef
[13]  R. B. Melrose et J. Sjöstrand, «Singularities of boundary value problems. II,» Communications on Pure and Applied Mathematics, vol. 35, n 12, pp. 129-168, 1982.
In article      CrossRef
[14]  C. Bardos, G. Lebeau et J. Rauch, «Sharp and sufficient conditions for the observation, control, and stabilization of waves from the boundary,» SIAM Journal, Control and Optimisation 30, pp. 1024-1065, 1992.
In article      CrossRef
[15]  C. Bardos, T. Masrour et F. Tatout, «Singularités du problème d'élastodynamique,» Comptes rendus de l'Académie des sciences. Série 1, Mathématique, 320 (9), pp. 1157-1160, 1995.
In article      
[16]  C. Bardos, T. Masrour et F. Tatout, «Condition nécessaire et suffisante pour la contrôlabilité exacte et la stabilisation du problème de l'élastodynamique,» Comptes rendus de l'Académie des sciences. Série 1, Mathématique, 320 (10), pp. 1279-12981, 1995.
In article      
[17]  L. Tartar, «H-measures, a new approach for studying homogenisation, oscillations and concentration effects in partial differential equations,» Proc. Roy. Soc. Edinburgh Sect. A, vol. 115, n 13-4, pp. 193-230, 1990.
In article      
[18]  L. Tartar, The general theory of homogenization: a personalized introduction, (Vol. 7). Springer., 2009.
In article      
[19]  C. Bardos et T. Masrour, «Mesures de défaut: observation et contrôle de plaques,» Comptes rendus de l'Académie des sciences. Série 1, Mathématique, 323 (6), pp. 621-626, 1996.
In article      
[20]  P. Grisvard, «Contrôlabilité exacte des solutions de l'équation des ondes en présence de singularités,» Journal de mathématiques pures et appliquées, vol. 68, n 12, pp. 215-259, 1989.
In article      
[21]  S. Nicaise, «About the Lamé system in a polygonal or a polyhedral domain and a coupled problem between the Lamé system and the plate equation. I: Regularity of the solutions,» Annali della Scuola Normale Superiore di Pisa-Classe di Scienze, vol. 19 (3), n 13, 1992.
In article      
[22]  A. Chaïra, «Equation des ondes et régularité sur un ouvert lipschitzien,» Comptes rendus de l'Académie des sciences, Série 1, Mathématique, vol. 316, n 11, pp. 33-36, 1993.
In article      
[23]  M. Cavalcanti, A. Khemmoudj et M. Medjden, «Uniform stabilization of damped Cauchy-Ventcel problem with variable coefficients and dynamic boundary conditions,» J. Math. Anal. Appl. 328, pp. 900-930, 2007.
In article      CrossRef
[24]  M. Cavalcanti et H. Oquendo, «Frictional versus viscoelastic damping in a semilinear wave equation,» SIAM J. Control Optim, vol. 42, n 14, p. 1310-1324, 2003.
In article      CrossRef
[25]  A. Heminna, «Exact controllability of the linear elasticity system with evolutive Ventcel conditions,» Portugaliae Mathematica. Nova Série, vol. 58, n 13, pp. 271-315, 2001.
In article      
[26]  I. Lasiecka et D. Tataru, «Uniform boundary stabilization of semilinear wave equation with nonlinear boundary damping,» Differential Integral Equations , vol. 6, p. 507-533, 1993.
In article      
[27]  I. Lasiecka et R. Triggiani, «Uniform exponential energy decay of wave equations in a bounded region with L^2 (0,∞,L^2 (Γ))-feedback control in the Dirichlet boundary conditions».
In article      
[28]  I. Lasiecka, R. Triggiani et P. Yao, «Inverse observability estimates for second-order hyperbolic equations with variable coefficients,» J. Math. Anal. Appl. 235 (1), p. 13-57, 1999.
In article      CrossRef
[29]  P. Grisvard, Singularities in boundary value problems, Vol. 22. Springer, 1992.
In article      
[30]  T. Masrour, «Convergence des fonctions propres de troisième espèce pour le laplacien,» Comptes rendus de l'Académie des sciences. Série 1, Mathématique,, pp. 309-312, 1995.
In article      
[31]  C. Bardos, T. Masrour et F. Tatout, «Observation and control of elastic waves. In Singularities and Oscillations,» Springer New York, pp. 1-16, 1997.
In article      CrossRef
  • CiteULikeCiteULike
  • Digg ThisDigg
  • MendeleyMendeley
  • RedditReddit
  • Google+Google+
  • StumbleUponStumbleUpon
  • Add to DeliciousDelicious
  • FacebookFacebook
  • TwitterTwitter
  • LinkedInLinkedIn