arxiv: v1 [math.ap] 27 Nov 2018

Transcription

1 arxiv: v1 [math.ap] 27 Nov 2018 Internal observability of the wave equation in tiled domains Anna Chiara Lai Dipartimento di Scienze di Base e Applicate per l Ingegneria, Sapienza Università di Roma via A. Scarpa, 16, Roma, Italy anna.lai@sbai.uniroma1.it November 29, 2018 Abstract We investigate the internal observability of the wave equation with Dirichlet boundary conditions in tilings. The paper includes a general result relating internal observability problems in general domains to their tiles, and a discussion of the case in which the domain is the triangle. Keywords: Internal observability, wave equation, Fourier series, tilings. Mathematics subject classification: 42B05, 52C20. 1 Introduction The aim of the present paper is to investigate internal observability properties of vibrating repetitive structures. Motivated by applications of hexagonal and triangular tilings (and related subtilings) to engineering, the particular case of the half to the equilateral triangle is treated in detail. By a repetitive structure, or tessellation, is meant a structure obtained by the assemblation of identical substructures, or tiles. For instance, two-dimensional lattices and the honeycomb lattice are examples of tessellation of R 2 ; while the regular hexagon (i.e., the tile of the honeycomb lattice) and the rectangle with aspect ratio equal to 3 are bounded domains that can be both tiled with triangles, see Figure 1. The interest in repetitive systems of vibrating membranes is motivated by applications in mechanical, civil and aerospace engineering This research is supported by Sapienza Università di Roma, Dipartimento di Scienze di Base e Applicate per l ingegneria, Assegno di Ricerca n. 7/

2 (a) Honeycomb tiling (b) Triangular tiling (c) Tiling of the regular hexagon (d) Tiling of the rectangle with aspect ratio equal to 3 Figure 1: Some tilings related to the triangle. [5, 24]. Modular structures have indeed the double advantage of a cost-effective manufacturing and construction (due to the repetitivity of the process) as well as a computationally cost-effective design. In particular, structural eigenproblems (e.g., vibrations and buckling) for repetitive structures in general involve a lower number of degrees of freedoms and, consequently, a less computationally demanding numerical solution [26]. Tilings involving regular triangles and hexagons (known as triangular lattice and honeycomb lattice, respectively) find countless applications in engineering, as well [27]. For instance, the use of such structures in architectural engineering is motivated by their mechanical properties, including resistance to external load and energy absorption, see for instance [6, 20] and, for a comprehensive dissertation on the topic, the book [7]. Finally, we mention that honeycomb lattice plays a crucial role in nanosciences and, in particular, in graphene technology [1]. As mentioned above, we are interested in the internal observability of the wave equation, that is the problem of reconstructing initial data from the observation of the evolution of the system in a subregion of the domain. Using folding and tessellation techniques, in the spirit of [23] and [18], we provide a general class of tilings, called admissible tilings, for which some internal observability properties of tiled domains extend to their tiles and under some symmetry assumptions on initial data vice versa. In particular, we show how to bridge the well-established theory concerning rectangular domains [8, 9, 15, 16, 17, 22] to the case of a triangular domain. In the remaining part of this Introduction we discuss in 2

3 K 3 (T 1 ) K 6 (T 1 ) K 2 (T 1 ) K 5 (T 1 ) T 1 K 4 (T 1 ) Figure 2: The tiling of R with T. Note that K 1 is the identity map, hence K 1 (T ) = T. detail this case, while postponing the more technical, general result to Section A case study: observability in a triangular domain We consider the problem u tt u = 0 in R T u = 0 in R T u(t, 0) = u 0, u t (t, 0) = u 1 in T (1) where T is the open triangle with vertices (0, 0), (1/ 3, 0) and (0, 1). Also consider ther rectangle R := (0, 3) (0, 1) and remark that there exists 6 rigid transformations K 1,..., K 6 satisfying the relation cl(r) = 6 K h (cl(t )), where cl(ω) represent the closure of a set Ω see Figure 2. We then say that T tiles 0 R. As it is well known, a complete orthonormal base for L 2 (R) is given by the eigenfunctions of in H 1 0 (R) e k := sin(πk 1 x 1 / 3) sin(πk 2 x 2 ), where k = (k 1, k 2 ), k 1, k 2 N 0 For a precise definition of tilings see Definition 2.1 below, while the explicit definition of the K h s is given in Section 3. 3

4 and the associated eigenvalues are γ k = k k 2 2. In [23], a folding technique (that we recall in detail in Section 3) is used to derive from {e k } an orthogonal base {e k } of L 2 (T ) formed by the eigenfunctions of in H 1 0 (T ). The explicit knowledge of a eigenspace for H 1 0 (T ) allows us to set the problem (1) in the framework of Fourier analysis see [14, 10, 8, 9, 4, 2, 3]. Our goal is to exploit the deep relation between the eigenfunctions for H 1 0 (R) and those of H 1 0 (T ) in order to extend known observability results for R to T. In particular, we are interested in the internal observability of (1), i.e., in the validity of the estimates T u 0 2 L 2 (T ) + u 1 2 H 1 (T ) 0 T 0 u(t, x) 2 dx where T 0 is a subset of T and T is sufficiently large. Here and in the sequel A B means c 1 A B c 2 A with some constants c 1 and c 2 which are independent from A and B. When we need to stress the dependence of these estimates on the couple of constants c = (c 1, c 2 ), we write A c B. Also by writing A c B we mean the inequality ca B while the expression A c B denotes ca B. We have Theorem 1.1. Let u be the solution of u tt u = 0 in R R u = 0 in R R u(t, 0) = u 0, u t (t, 0) = u 1 in R, (2) let R 0 be a subset of R and assume that there exists a constant T 0 0 such that if T > T 0 then there exists a couple of constants c = (c 1, c 2 ) such that u satisfies u 0 2 L 2 (R) + u 1 2 H 1 (R) c T for all (u 0, u 1 ) L 2 (R) H 1 (R). Moreover let T 0 := 6 0 K 1 h (R 0) T. R 0 u(t, x) 2 dxdt (3) Then for each T > T 0 and (u 0, u 1 ) L 2 (T ) H 1 (T ), the solution u of (1) satisfies T u 0 2 L 2 (T ) + u 1 2 H 1 (T ) c u(t, x) 2 dxdt (4) 0 T 0 The result also holds by replacing every occurrence of c with c or c. 4

5 We point out that the time of observability T 0 stated in Theorem 1.1, as well as the couple c of constants in the estimates (3) and (4), are the same for both the domains R and T. Also note that in Section 3 we prove a slightly stronger version of Theorem 1.1, that is Theorem 3.5: its precise statement requires some technicalities that we chose to avoid here, however we may anticipate to the reader that the assumption on initial data (u 0, u 1 ) L 2 (R) H 1 (R) can be weakened by replacing L 2 (R) H 1 (R) with an appropriate subspace. 1.2 Organization of the paper. In Section 2 we consider a generic domain Ω tiling a larger domain Ω : we establish a result, Theorem 2.10, relating the observability properties of wave equation on Ω and on its tile Ω. Section 3 is devoted to the proof of Theorem An observability result on tilings The goal of this section is to state an equivalence between an observability problem on a domain Ω and an observability problem on a larger domain Ω, under the assumption that Ω tiles Ω. We begin with some definitions. Definition 2.1 (Tiling). Let Ω and Ω be two open bounded subsets of R n. We say that Ω tiles Ω if there exists a set {K h } N rigid transformations of Rn such that N cl(ω ) = K h (cl(ω)) and such that K h (Ω) K j (Ω) = for all h j. Definition 2.2 (Foldings and prolongations). Let (Ω, {K h } N ) be a tiling of Ω and δ = (δ 1,..., δ N ) { 1, 1} N. The prolongation with coefficients δ of a function u : Ω R to Ω is the function P δ u : Ω R P δ u(k h x) = δ h u(x) h = 1,..., N. The folding with coefficients δ of a function u : Ω R is the function F δ u : Ω R defined by F δ u(x) = 1 N N 2 δ h u(k h x) h = 1,..., N. When the particular choice of δ is not relevant we omit it in the under scripts and we simply write P and F. Definition 2.3 (Admissible tiling). A tiling (Ω, {K h } N ) of Ω is admissible if there exists δ { 1, 1} N such that F δ ϕ H 1 0 (Ω) ϕ H 1 0 (Ω ). (5) 5

6 T K 2T Figure 3: A non-admissible tiling of R = (0, 1/ 3) (0, 1) with T. Example 2.4. We show in Lemma 3.1 below that the tiling of R with T depicted in Figure 2 is admissible, in particular (5) holds with δ = (1, 1, 1, 1, 1, 1). On the other hand the tiling of R := (0, 1/ 3) (0, 1) given by the transformations K 1 := id and K 2 : (x 1, x 2 ) (x 1, x 2 ) + (1/ 3, 1), see Figure 3, is not admissible. Let indeed v 1 := (1/ 3, 0), v 2 := (0, 1) and x λ := λv 1 + (1 λ)v 2 with λ (0, 1). Then x λ T and K 2 (x λ ) = x 1 λ Therefore it suffices to choose ϕ H 1 0 (R) such that ϕ(x λ ) ±ϕ(x 1 λ ) to obtain F δ ϕ(x λ ) = δ 1 ϕ(x λ ) + δ 2 ϕ(x 1 λ ) 0 for all δ 1, δ 2 { 1, 1}. Consequently F δ ϕ H 1 0 (T ) for all δ { 1, 1} 2. Remark 2.5. We borrowed the notion of prolongation and folding from [23]: while our definition of P δ is exactly as it is given in [23], we introduced a normalizing term 1/N 2 in the definition of F δ in order to enlighten the notations. Note that the following equality holds: for all u : Ω R. F δ (P δ u) = 1 N u (6) 6

7 Also remark that we shall need to prolong and fold also functions u : R Ω R and ū : R Ω R, in this case the definition of P and F naturally extends by applying the transformations K h s to the spatial variables x. For instance if u : R Ω R then its prolongation to R Ω reads P δ u(t, K h x) = δ h u(t, x). We want to establish a relation between solutions of a wave equation with Dirichlet boundary conditions and their prolongation. To this end we introduce the notations P δ L 2 (Ω) := {P δ u u L 2 (Ω)}, and P δ H 1 0 (Ω) := {P δ u u H 1 0 (Ω)} P δ H 1 (Ω) := {P δ u u H 1 (Ω)}. Note that P δ L 2 (Ω) L 2 (Ω ), P δ H0 1 (Ω) H0 1 (Ω ) and P δ H 1 (Ω) H 1 (Ω). All results below hold under the following assumptions on the domains Ω, Ω and on a base {e k } for L 2 (Ω): Assumption 1. (Ω, {K h } N ) is an admissible tiling of Ω. Assumption 2. {e k } is a base of eigenvectors of in H 1 0 (Ω), it is defined on Ω Ω and there exists δ { 1, 1} N such that for each k N. P δ (e k Ω ) = e k Ω Remark 2.6 (Some remarks on Assumption 2). We note that Assumption 2 can be equivalently stated as e k (K h x) = δ h e k (x) for all x Ω, h = 1,..., N, k N. (7) Indeed, by definition of prolongation and noting δh 2 1, we have e k (K h x) = δ 2 he k (K h x) = δ h P δ e k (x) = δ h e k (x). for every x Ω, h = 1,..., N and k N. Also remark that, in view of (6), Assumption 2 also implies F δ e k = 1 N e k. (8) Example 2.7. Let Ω = (0, π) 2 and Ω = (0, 2π) 2. Consider the transformations of R 2 K1 := id, K2 : (x1, x2) ( x1 + 2π, x2), K 3 : (x 1, x 2 ) (x 1, x 2 + 2π), K 4 : (x 1, x 2 ) (x 1, x 2 ) + (2π, 2π). 7

8 Then {Ω, {K h } 4 } is a tiling for Ω. In particular, Assumption 1 is satisfied: indeed setting δ = (1, 1, 1, 1) we have for each ϕ H 1 0 (Ω ) F δ ϕ(x) = 0 x Ω. Also note that the functions e k (x) := sin(k 1 x 1 ) sin(k 2 x 2 ) k = (k 1, k 2 ) N 2 satisfy Assumption 2, indeed they are a base for L 2 (Ω) composed by eigenfunctions of in H 1 0 (Ω) and e k (K h (x)) := δ h e k (x) for all x R 2, h = 1,..., 4 and k N 2. The space P δ L 2 (Ω) in this case coincides with the space of so-called (2, 2)-cyclic functions, i.e., functions in L 2 (Ω ) which are odd with respect to both axes x 1 = π and x 2 = π. We refer to [16] for some results on observability of wave equation with (p, q)-cyclic initial data. Our starting point is to show that, under Assumption 1 and Assumption 2, the base of eigenfunctions {e k } is also a base of eigenfunctions also for an appropriate subspace of L 2 (Ω ), and to compute the associated coefficients. Lemma 2.8. Let Ω, Ω and {e k } satisfy Assumption 1 and Assumption 2. Then {e k } H0 1 (Ω ) and it is also a complete base for P δ L 2 (Ω) formed by eigenfunctions of in P δ H0 1 (Ω ). In particular, for every k N, if u k is the coefficient of u L 2 (Ω) (with respect to e k ) then Nu k is the coefficient of P δ u. Proof. The proof is organized two steps. Claim 1: {e k } is a set of eigenfunctions of in H0 1 (Ω ). Extending a result given in [23], we need to show that, under Assumption 1 and Assumption 2, if e k H0 1 (Ω) is a solution of the boundary value problem e k ϕdx = γ k e k ϕdx ϕ H0 1 (Ω) Ω Ω for some γ k R, then e k is also solution of the boundary value problem on Ω e k ϕdx = Ω γ k e k ϕdx Ω ϕ H0 1 (Ω ). Now, recall from Assumption 1 that if ϕ H 0 1 (Ω ) then F δ ϕ H 0 1 (Ω). Then it follows again from Assumption 1 and from Assumption 2 (in particular by recalling 8

9 that K h s are isometries and (7)) that for all ϕ H0 1 (Ω ) e k (x) ϕ(x)dx = e k (x) ϕ(x)dx Ω N K h(ω) N N = e k (K h x) ϕ(k h x)dx = e k (x) δ h ϕ(k h x)dx Ω = e k (x) F δ ϕ(x)dx = Ω = γ k e k (x)ϕ(x)dx. Ω Ω Ω γ k e k (x)f δ ϕ(x)dx and this completes the proof of Claim 1. Claim 2: completeness of {e k } and computation of coefficients By Assumption 1 and Assumption 2 and by recalling δh 2 = 1 for each h = 1,..., N, we have P δ u(x)e k (x)dx = P δ u(x)p δ e k (x)dx Ω Ω N = P δ u(x)pe k (x)dx = = N N K h (Ω) K h (Ω) Ω δ 2 hu(k h (x))e k (K h (x))dx u(x)e k (x)dx = N Ω u(x)e k (x)dx, where the second to last equality holds because K h s are isometries. Then we may deduce two facts: first if {u k } are the coefficients of u L 2 (Ω) then {Nu k } are coefficients of P δ u. Secondly, {e k } is a complete base for P δ L 2 (Ω), indeed if the coefficients of P δ u are identically null, then also the coefficients of u are identically null: since {e k } is complete for Ω then u 0 and, consequently, P δ u 0, as well. Next result establishes a relation between solutions of wave equations on tiles and their prolongations. Lemma 2.9. Let Ω, Ω and {e k } satisfy Assumption 1 and Assumption 2. Let u be the solution of u tt u = 0 in R Ω u = 0 in R Ω (9) u(t, 0) = u 0, u t (t, 0) = u 1 in Ω 9

10 Then u is well defined in Ω Ω and u = Nu Ω is the solution of u tt u = 0 in R Ω u = 0 in R Ω (10) u(t, 0) = P δ u 0, u(t, 0) = P δ u 1 in Ω Conversely, if ū is the solution of (10) then F δ ū is the solution of (9) and for every h = 1,..., N F δ ū(t, x) = δ h N ū(t, K hx) for each x Ω. (11) Proof. Let {γ k } be the sequence of eigenvalues associated to {e k } and set ω k = γk, for every k N. Expanding u(t, x) with respect to e k we obtain u(t, x) = (a k e iωkt + b k e iωkt )e k (x) k=1 with a k and b k depending only the coefficients c k and d k of u 0 and u 1 with respect to {e k }. In particular a k + b k = c k and a k b k = id k /ω k. We then have that the natural domain of u coincides with the one of {e k } s, hence it is included in Ω Ω. By Lemma 2.8 the coefficients of P δ u 0 and P δ u 1 are Nc k and Nd k, respectively. Then it is immediate to verify that Nu(t, x) = (Na k e iωkt + Nb k e iωkt )e k (x) is the solution of (10). Now, let ū(t, x) = k=1 (ā k e iωkt + b k e iωkt )e k (x) k=1 be the solution of (10), and note that, by the reasoning above, setting a k := 1 N āk and b k := 1 N b k we have that u(t, x) := (a k e iωkt + b k e iωkt )e k (x) = 1 ū(t, x) N k=1 is the solution of (9). Hence to prove that u(t, x) = F δ ū(t, x) it it suffices to note that by Assumption 1 (see in particular (8)) F δ ū(t, x) = (ā k e iωkt + b k e iωkt )F δ e k (x) k=1 = 1 N (ā k e iωkt + b k e iωkt )e k (x) = 1 ū(t, x). N k=1 10

11 Finally, we show (11): for each h = 1,..., N we have ū(t, x) = δ 2 hū(t, x) = δ h = k=1 (ā k e iω kt + b k e iω kt )δ h e k (x) (ā k e iωkt + b k e iωkt )e k (K h x) = δ h ū(t, K h x) k=1 and this concludes the proof. (10). We are now in position to state the main result of this section, that bridges observability of tiles with their prolongations. Theorem Let Ω, Ω and {e k } satisfy Assumption 1 and Assumption 2. Let u be the solution of u tt u = 0 in R Ω u = 0 in R Ω (12) u(t, 0) = u 0, u t (t, 0) = u 1 in Ω with u 0, u 1 L 2 (Ω) H 1 (Ω) and let u be the solution of u tt u = 0 in R Ω u = 0 in R Ω u(t, 0) = P δ u 0, u t (t, 0) = P δ u 1 in Ω. (13) Also let Ω 0 Ω and define Ω 0 := N K 1 h (Ω 0) Ω. Then for every T > 0 and for every couple c = (c 1, c 2 ) of positive constants, the inequalities T u 0 2 L 2 (Ω) + u 1 2 H 1 (Ω) c u(t, x) 2 dxdt. (14) Ω 0 hold if and only if P δ u 0 2 L 2 (Ω ) + P δu 1 2 H 1 (Ω ) c Proof. By Lemma 2.9, u and u satisfy 0 T u(t, x) = δ h N u(t, K hx) for all h = 1,..., N Ω 0 u(t, x) 2 dxdt. (15)

12 Since Ω tiles Ω, then setting Ω h := K 1 h (Ω ) Ω we have Ω 0 = N Ω h and Ω 0 = N K h(ω h ), and that these unions are disjoint. Hence, also recalling δ h 1 and that K h s are isometries, we have N u(t, x) 2 dx = u(t, x) 2 dx I Ω 0 I K h (Ω h ) N = u(t, K h (x)) 2 dx I Ω h N 2 = N 2 δ h N u(t, K h(x)) dx I N = N 2 = N 2 I I Ω h Ω h u(t, x) 2 dx Ω 0 u(t, x) 2 dx Finally, by Lemma 2.8 P δ u 0 2 L 2 (Ω ) = N 2 u 0 2 L 2 (Ω) and P δ u 1 2 H 1 (Ω ) = N 2 u 1 2 H 1 (Ω) and this implies the equivalence between (14) and (15). 3 Proof of Theorem 1.1 The proof of Theorem 1.1 is based on the application of Theorem 2.10 to the particular case Ω = T and Ω = R. We then need to admissibly tile R with T and a base {e k } formed by the eigenfunctions of in H0 1 (T ) satisfying Assumption 2. Such ingredients are provided in [23]: in order to introduce them we need some notations. We consider the Pauli matrix ( ) 1 0 σ z := 0 1 and the rotation matrix ( ) cos α sin α R α := sin α cos α where α := π/3. Now let v 1 := (0, 1/ 3) and v 2 := (0, 1) be two of the three vertices of T and define the transformations from R 2 onto itself K 1 := id; K 4 : x R α (x v 2 ) + 3v 1 K 2 : x R α σ z (x v 2 ) + v 2 ; K 5 : x R α (x v 2 ) + 3v 1 + v 2 (16) K 3 : x R α (x v 2 ) + v 2 ; K 6 : x x + 3v 1 + v 2 12

13 Figure 4: The tiling of R with T, the grey areas correspond to negative δ h s. and note (T, {K h } 6 ) is a tiling for R. Indeed cl(r) = 6 K h cl(t ), (17) and the sets K h T, for h = 1,..., 6, do not overlap see Figure 4 and [23]. We set δ := (1, 1, 1, 1, 1, 1). and, in next result, we prove that T admissibly tiles R. Lemma 3.1. (T, {K h } 6 ) is an admissible tiling of R. Proof. We want to show that if ϕ H 1 0 (R) then F δ ϕ H 1 0 (T ). To this end let v 0 := (0, 0), v 1 := (1/ 3, 0) and v 2 := (0, 1) be the vertices of T and define x λ ij := λv i + (1 λ)v j. so that T = {x λ ij λ [0, 1], 0 i < j 2}. By a direct computation, for all λ [0, 1] K 1 (x λ 01), K 6 (x λ 01) R, and K 2 (x λ 01) = K 4 (x λ 01), K 3 (x λ 02) = K 5 (x λ 02). Since ϕ H 1 0 (R) then F δ ϕ(x λ 01) = 0. Similarly, for all λ [0, 1] K 1 (x λ 02), K 6 (x λ 02) R, K 2 (x λ 02) = K 3 (x λ 02), 13

14 and K 4 (x λ 02) = K 5 (x λ 02) therefore F δ ϕ(x λ 02) = 0 for all λ [0, 1]. Finally for all λ [0, 1] and K 3 (x λ 12), K 4 (x λ 12) R, K 1 (x λ 12) = K 2 (x λ 12), K 5 (x λ 12) = K 6 (x λ 12) therefore we get also in this case F δ ϕ(x λ 12) = 0 for all λ [0, 1] and we may conclude that F δ ϕ H 1 0 (T ). Remark 3.2. Lemma 3.1 was remarked in [23, p.312], but to the best of our knowledge, this is the first time an explicit proof is provided. Now, consider the eigenfunctions of in H 1 0 (R): e k (x 1, x 2 ) := sin(πk 1 x 1 3 ) sin(πk 2 x 2 ), k = (k 1, k 2 ) N 2. We finally define for every k N 2 e k (x) := N 2 F δ e k = 6 δ h e k (K h x). (18) Next result, proved in [23], states that Assumption 2 is satisfied by {e k }. Lemma 3.3. The set of functions {e k } defined in (18) is a complete orthogonal base for T formed by the eigenfunction of in H 1 0 (T ). Furthermore P δ e k (x) = e k (x). Remark 3.4. For each k N 2, the eigenfunctions e k and ē k share the same eigenvalue γ k = π 2 ( k k2), 2 see [23]. Next gives access to classical results on observability of rectangular membranes for the study of triangular domains. Theorem 3.5. Let u be the solution of (1) with u 0, u 1 L 2 (Ω) H 1 (Ω) and let u be the solution of u tt u = 0 on R R u = 0 in R R u(t, 0) = P δ u 0, u t (t, 0) = P δ u 1 in R. 14

15 Also let R 0 R and define T 0 := N K 1 h (T 0) Ω. Then for every T > 0 and for every couple c = (c 1, c 2 ) of positive constants, the inequalities hold if and only if u 0 2 L 2 (T ) + u 1 2 H 1 (T ) c T P δ u 0 2 L 2 (R) + P δu 1 2 H 1 (R) c 0 T 0 T 0 u(t, x) 2 dxdt. (19) R 0 u(t, x) 2 dxdt. (20) Proof. Since T, R and {e k } satisfy Assumption 1 and Assumption 2, then the claim follows by a direct application of Theorem 2.10 with Ω = T and Ω = R. We conclude this section by showing that Theorem 1.1 is a direct consequence of Theorem 3.5: Proof of Theorem 1.1. By Lemma 2.8, if (u 0, u 1 ) L 2 (T ) H 1 (T ) then (P δ u 0, P δ u 1 ) L 2 (R) H 1 (R). The claim hence follows by Theorem Acknowledgements The author is grateful to Vilmos Komornik (Université de Strasbourg) and Paola Loreti (Sapienza Università di Roma) for the fruitful discussions on the topic of internal observability of wave equations and related symmetry properties of the initial data: they provided insight and expertise that greatly assisted the research. References [1] M. J. Allen, V. C. Tung, and R. B. Kaner. Honeycomb carbon: a review of graphene. Chemical reviews 110.1, (2009). [2] C. Baiocchi, V. Komornik, P. Loreti, Ingham type theorems and applications to control theory, Bol. Un. Mat. Ital. B (8) 2 (1999), no. 1, [3] C. Baiocchi, V. Komornik, P. Loreti, InghamBeurling type theorems with weakened gap conditions, Acta Math. Hungar. 97 (2002), 1-2, [4] J.N.J.W.L. Carleson and P. Malliavin, editors, The Collected Works of Arne Beurling, Volume 2, Birkhäuser,

16 [5] D. A. Evensen.Vibration analysis of multi-symmetric structures. AIAA Journal 14.4 (1976): APA [6] S. Frølich, et al. Uncovering Nature s Design Strategies through Parametric Modeling, Multi-Material 3D Printing, and Mechanical Testing. Advanced Engineering Materials 19.6 (2017): e [7] L.J. Gibson, and M. F. Ashby. Cellular solids: structure and properties. Cambridge university press, [8] A. Haraux, Sries lacunaires et contrle semi-interne des vibrations d une plaque rectangulaire, J. Math. Pures Appl. 68 (1989): [9] A. Haraux, On a completion problem in the theory of distributed control of wave equations, Nonlinear partial differential equations and their applications. Collge de France Seminar, Vol. X (Paris, ), , Pitman Res. Notes Math. Ser., 220, Longman Sci. Tech., Harlow, 1991 [10] A.E. Ingham, Some trigonometrical inequalities with applications in the theory of series, Math. Z. 41 (1936): [11] E. G. Karpov, N. G. Stephen, and D. L. Dorofeev. On static analysis of finite repetitive structures by discrete Fourier transform. International Journal of Solids and Structures (2002): [12] V. Komornik, Rapid boundary stabilization of the wave equation, SIAM J. Control Optim., 29 (1991): [13] V. Komornik, P. Loreti, Ingham type theorems for vector-valued functions and observability of coupled linear systems, SIAM J. Control Optim. 37 (1998): [14] V. Komornik, P. Loreti, Fourier Series in Control Theory, Springer-Verlag, New York, [15] V. Komornik, P. Loreti, Multiple-point internal observability of membranes and plates, Appl. Anal. 90 (2011), 10, [16] V. Komornik, P. Loreti, Observability of rectangular membranes and plates on small sets, Evol. Equations and Control Theory 3 (2014), 2, [17] V. Komornik, B. Miara, Cross-like internal observability of rectangular membranes, Evol. Equations and Control Theory 3 (2014), 1, [18] J.R. Kuttler, and V. G. Sigillito, Eigenvalues of the Laplacian in two dimensions, Siam Review 26(1984), 2, [19] P. Loreti, M. Mehrenberger, An Ingham type proof for a two-grid observability theorem, ESAIM Control Optim. Calc. Var. 14 (2008), no. 3,

17 [20] I.G. Masters, and K. E. Evans. Models for the elastic deformation of honeycombs. Composite structures 35.4 (1996): [21] P. Martinez. Boundary stabilization of the wave equation in almost star-shaped domains. SIAM journal on control and optimization 37 (1999) 3, [22] M. Mehrenberger, An Ingham type proof for the boundary observability of a Nd wave equation. C. R. Math. Acad. Sci. Paris 347 (2009), no. 1-2, [23] M. Práger. Eigenvalues and eigenfunctions of the Laplace operator on an equilateral triangle. Applications of mathematics, 43(4), (1998) [24] D. L. Thomas. Dynamics of rotationally periodic structures. International Journal for Numerical Methods in Engineering 14.1 (1979): APA [25] J-C. Sun. On approximation of Laplacian eigenproblem over a regular hexagon with zero boundary conditions. J. Comp. Math., international edition, 22.2 (2004), [26] D. Wang, C. Zhou, and J. Rong. Free and forced vibration of repetitive structures. International Journal of Solids and Structures (2003): APA [27] Q. Zhang, et al. Bioinspired engineering of honeycomb structure Using nature to inspire human innovation. Progress in Materials Science 74 (2015):