DECAY RATES FOR ELASTIC-THERMOELASTIC STAR-SHAPED NETWORKS

. This work discusses the asymptotic behaviour of a transmission problem on star-shaped networks of interconnected elastic and thermoelastic rods. Elastic rods are undamped, of conservative nature, while the thermoelastic ones are damped by thermal eﬀects. We analyse the overall decay rate depending of the number of purely elastic components entering on the system and the irrationality properties of its lengths. First, a suﬃcient and necessary condition for the strong stability of the thermoelastic-elastic network is given. Then, the uniform exponential decay rate is proved by frequency domain analysis techniques when only one purely elastic undamped rod is present. When the network involves more than one purely elastic undamped rod the lack of exponential decay is proved and nearly sharp polynomial decay rates are deduced under suitable irrationality condi- tions on the lengths of the rods, based on Diophantine approximation arguments. More general slow decay rates are also derived. Finally, we present some numerical simulations supporting the analytical results.

A closely related interesting issue is the asymptotic behaviour of a system consisting of two different materials joined together at the interface, one being purely elastic and the other thermoelastic (see Figure 1). In these systems, other than the intrinsic coupling effects of thermal and elastic components, typical in thermoelasticity, the thermoelastic and purely elastic rods are also coupled through the interface. Exponential and polynomial decay properties for this kind of systems were proved by Rivera et al. in [12], [28], [31] and Messaoudi et al. in [30], using the energy multiplier method. Han and Xu in [15] got a sharp polynomial decay rate for a thermoelastic transmission problem with joint mass, based on a detailed spectral analysis and resolvent operator estimates. We also refer to [29] for the analysis of the large time behaviour of transmission problems of multi-dimensional thermoelasticity.
In this work, we consider similar transmission problems in multi-connected networks. More precisely, we are interested in the large time behaviour of star-shaped networks constituted by coupled thermoelastic and purely elastic rods (see Figure  2).
The natural energy of this system is as follows Hence, the energy of system (5) is decreasing. Moreover, from (6), it is easy to see that the dissipation mechanism only acts in the thermoelastic rods. This motivates the problem of whether or not the dissipation is strong enough to make the total energy of the network decay to zero, and with which rate.
In this paper, the large time behaviour of system (5) is mainly discussed based on frequency domain analysis ( [5], [6], [13], [19], [25], [26] and [34]). In [37], Shel showed the exponential stability of networks of thermoelastic and elastic materials for some special cases by similar methods. However, it was assumed that there was no heat exchange between the thermoelastic rods connected at the common nodes. In this paper, we allow for the heat exchange between thermoelastic rods at common nodes, which is a natural assumption.
By estimating the resolvent operator along the imaginary axis and employing multiplier techniques, we get a necessary and sufficient condition for system (5) to decay uniformly exponentially, namely that there is no more than one purely elastic rod entering in the network. If this condition fails, the system lacks exponential decay and we further show that the decay rate of the networks can not be faster than t −1 . Moreover, for a very special case that there are two purely elastic rods involved in the network, the optimal polynomial decay result is obtained and for the general case that there are more purely elastic rods involved, a nearly optimal polynomial decay result is also derived under certain Diophantine approximation conditions, that refer to the lengths of the purely elastic rods involved in the network.
To discuss the sharpness of slow decay rates it is useful to get explicit information on the spectrum of the system, and compare its real and imaginary parts (see [6] and [41]). However, spectra of PDE networks are often difficult to calculate. Thus, we prove the optimality by estimating the norm of the resolvent operator along the imaginary axis (see [1]). But resolvent estimates are hard to be achieved due to the thermoelastic coupling. Thus, we employ diagonalisation argument to deal with the resolvent problem. This allows building explicit approximations of solutions ensuring that the polynomial decay rates we get are nearly optimal.
The rest of the paper is organised as follows. In section 2, the main result of this paper is given. Section 3 is devoted to show the well-posedness and strong asymptotic stability of the system (5). In section 4, we prove the exponential and nearly optimal slow decay rates, under different conditions. Section 5 is devoted to discuss some more general slow decay rates for system (5). Especially, for the special case N − N 1 = 2, the condition to achieve optimal decay of the network is obtained. In section 6, the numerical simulations of the dynamical behaviour of system (5) are presented.
The results in this paper contribute to the understanding of the decay properties of wave and thermoelastic wave networks, a topic in which important issues are still to be understood.
There has been an extensive literature on other closely related issues such as the large time behaviour and controllability properties of elastic networks with node and boundary feedback controls. We refer, among others (the present list of references is by no means complete), to Lagnese et al. [21] for the modelling and control of elastic networks; Ammari et al. [2], [3] and [4] and Nicaise et al. [32] for stabilisation problems on networks of wave and Euler-Bernoulli beams with star-shaped and tree-shaped configurations; Dáger and Zuazua [8], [9] and [10] for boundary controllability of wave networks; Xu et al. [16], [17] and [40] for the stabilisation and spectral properties of the wave networks.
2. Main results. This section is devoted to state the main result of this paper.
As we will see later, system (5) can be rewritten as an abstract Cauchy problem in an appropriate Hilbert space H: where U (t) = (u, u t , θ) T and U (0) = (u (0) , u (1) , θ (0) ) T ∈ H are given. The problem of whether the energy of solutions tends to zero as time goes to infinity or not has a simple answer: Theorem 2.1. Operator A generates a C 0 semigroup of contractions on H. Moreover, the energy of the system (5) decays to zero as t → ∞ if and only if one of the following two conditions is fulfilled, We now obtain explicit decay rates for network (5).
It is well known that if all the components in the network are thermoelastic, that is N = N 1 , the energy of the system decays exponentially to zero.
In fact, in Propositions 1 and 2 in section 4, we obtain the following necessary and sufficient condition for the exponential decay.
Theorem 2.2. The energy of system (5) decays to zero exponentially if and only if N − N 1 ≤ 1, that is, if no more than one purely elastic undamped rod is involved in the network.
Accordingly, when N − N 1 > 1, one can only expect slow decay rate for network (5). In order to address this issue we need the following definition from [36] and [11]. [11]) Real numbers 1 , 2 , · · · , m are said to verify the conditions (S), if 1 , 2 , · · · , m are linearly independent over the field Q of rational numbers; and the ratios i / j are algebraic numbers for i, j = 1, 2, · · · , m.
The notation (S) for this condition was introduced in [11] to refer to the fundamental contribution by Schmidt [36], that defined this class of irrational numbers for the simultaneous approximation by rational ones.
It should be noted that the condition (S) denotes a narrow class of irrationals, since the set of algebraic numbers is countable and has Lebesgue measure zero.
By a detailed frequency domain analysis, we have the following explicit polynomial decay rate for system (5).
Remark 2. Using the method of proof of Theorem 2.4, we can obtain more general slow decay rates (polynomial, logarithmic or arbitrarily slow decay), which will be presented in section 5. Especially, for a very special case that there are two purely elastic rods involved in the network (N − N 1 = 2), we show that the network can achieve the optimal decay rate t −1 if the mutual ratio of the lengths of the purely elastic rods belongs to the set of irrational numbers having a continuous fraction expansion [a 0 , a 1 , ..., a n , ...] with bounded (a n ). (8) is well-known for wave-like equations with velocity damping in the case of one single string with damping on an internal point, which is equivalent to a simple star-like network constituted only by two strings ( [20], [11]).

Remark 3. The optimality result in
Here we show the same lower bound on the decay rate for the more general system involving thermoelastic rods. However, for general cases that more than two purely elastic rods entering in the networks, it is still an open problem to find the condition to guarantee achieving a sharp polynomial decay rate.
3. Well-posedness and strong stability. This section is devoted to show the well-posedness of network (5) by the semigroup theory and prove Theorem 2.1, with a necessary and sufficient condition for strong stability of this system.
Let us first introduce an appropriate Hilbert space setting for the well-posedness of the system.
(H, · H ) is a Hilbert space. Then, define the system operator A in H as follows: Thus, system (5) can be rewritten as the evolution equation (7) in H.
It is easy to check that A is dissipative in H. Moreover, A is injective and surjective and hence 0 ∈ ρ(A). Then, by Lummer-Phillips theorem (see [33]), A generates a C 0 semigroup of contractions S(t) on H. Now, we focus on proving the strong stability property of the system in Theorem 2.1. The proof by contradiction is mainly used here. Proof of Theorem 2.1.
Remark 4. For the proof of "Sufficiency", we also can use the unique continuation property for wave networks in Dager and Zuazua [11] (See Corollary 5.28, p.135). Indeed, (10), in the absence of thermal components, corresponds to the eigenproblem associated with the pure wave system and, according to the results in [11], its unique solution is the trivial one, which contradicts that (u T is an eigenvector. Then the desired result follows. 4. Decay rates. This section is devoted to achieve explicit decay rates of the total energy of solutions of system (5). The exponential decay and slow decay rates are deduced under different assumptions of the various components of the network. 4.1. Exponential decay rate: Case N − N 1 = 1. In this subsection, we analyse the decay rate of network (5) when N − N 1 = 1, namely, there is only one purely elastic rod involved in the network. To do this, let us introduce the following lemma as described in [13], [19], [25] and [34]. and We are then in conditions to prove the following proposition, which is one of the main statements in Theorem 2.1: Proposition 1. When N − N 1 = 1, the energy of system (5) decays exponentially to zero.
Proof. It is sufficient to show that the conditions in Lemma 4.1 are fulfilled. It should be noted that although the idea of this proof is similar to the one in [37], some different multipliers are employed to get certain estimates so as to deal with the transmission conditions in the present paper.
By the argument of the proof of Theorem 2.1, it is easy to see that, in the present case, there is no eigenvalue of A on the imaginary axis, that is, the condition (11) in Lemma 4.1 holds. Now we proceed to prove condition (12) and to deduce the exponential decay rate.
Taking the inner product of (16) Note that Thus, by (17), (24) and (26), we get Φ n → 0, in H, which is in contradiction with Φ n = 1. The desired result follows.

4.2.
Lack of exponential decay: Case N − N 1 > 1. In this subsection, we shall show that if more than one purely elastic undamped rod is involved in the network, the exponential decay rate does not hold. Proposition 2. If N − N 1 > 1 the network (5) lacks the property of exponential decay.
Proof. Note that from Lemma 4.1, it is sufficient to show that the norm of the resolvent operator of system (5) along the imaginary axis is necessarily unbounded when N − N 1 > 1.
We rewrite (35) in the vector form In order to diagonalize the matrix A j , we employ the transformation Y j := P j Z j , j = 1, 2, · · · , N 1 , where We then have In this way we obtain the following system By the boundary and transmission conditions in (5), together with (34), we get the following estimate (the technical details are given as in Appendix): Choosing σ n N1+1 = nπ, n → +∞, then sin σ n N1+1 = 0. Since the irrational numbers always can be approximated by rational ones, we can find a subsequence σ n k , such that sin σ n k j → 0, j = N 1 + 2, N 1 + 3, · · · , N.

4.3.
Lower bounds on the polynomial decay rate. This subsection is devoted to get the lower bounds on the polynomial decay rate as stated in Theorem 2.4.
We need the following result from [6] (see also [26]). Proposition 3. The decay rate of the energy of system (5) can be, at most, polynomial of order t −1 .

ZHONG-JIE HAN AND ENRIQUE ZUAZUA
Proof. From Lemma 4.2, it is sufficient to show that there exists at least one sequence (σ n , F n ) such that σ n → +∞, n → ∞ and where F n ∈ H and F n H is bounded, C being some positive constant. For simplicity, we consider the case N − N 1 = 2, that is, there are two purely elastic rods involved in the network. Based on the proof of Proposition 2, we choose the sequence (σ n k , F n k ) satisfying (38) in the following three steps.

4.4.
Explicit polynomial decay rate: Case N − N 1 > 1 (Proof of Theorem 2.4). In subsection 4.2 we proved that if there are more than one purely elastic rods involved in the network (5), the system can not achieve the exponential decay. We have also shown that the decay rate can be, at most, polynomial of order t −1 .
The mutual-irrationality of the radii of j , j = N 1 + 1, N 1 + 2, · · · , N in Theorem 2.1 is a sufficient and necessary condition for strong stability. So, for the rest of this section, we always assume this condition to be fulfilled. Now, we prove effective and explicit polynomial decay results when N − N 1 > 1, and j , j = N 1 + 1, N 1 + 2, · · · , N satisfy the conditions (S) in Definition 2.3.

Proof of Theorem 2.4.
Similarly to the proof of Theorem 2.2, we argue by contradiction, on the basis of Lemma 4.2.
Let Φ n be a sequence such that Φ n = (U n , V n , Θ n ) ∈ D(A), with Φ n H = 1, where U n = (u n j ) N j=1 , V n = (v n j ) N j=1 , Θ n = (θ n j ) N1 j=1 , and a sequence σ n ∈ R with σ n → ∞ such that where is given as in Theorem 2.4. Hence, and (46) Thus, by Poincaré inequality, σ (1+ ) n θ n j → 0, in L 2 (0, j ), j = 1, 2, · · · , N 1 . Similar to the proof of Theorem 2.2, we can get that (48) Thus, by the transmission conditions at the common node, we get Then we will prove that Φ n → 0 in H, which leads to a contradiction. To do this, we discuss it by the following two cases: Case 1). Assume that sin σ n j = 0, j = N 1 + 1, N 1 + 2, · · · , N . By (42), (45), we can easily get that u n j,xx + σ 2 n u j = − g n j +iσnf n j σ 2(1+ ) n .
Note that there exists at most one j n 0 , N 1 + 1 ≤ j 0 ≤ N satisfying sin σ n j0 = 0 due to the irrationality of the mutual ratii between j , N 1 + 1 ≤ j ≤ N . In this case, where γ n is given as in (53). Thus, by the same arguments as in Step 1), we can get u n j,x (0) → 0, N 1 + 1 ≤ j ≤ N, j = j n 0 , n → ∞, which together with the transmission condition (50), implies

Then by (57)-(59) in
Step 1), we get u n j,x , v n j → 0, in L 2 (0, j ), j = N 1 + 1, N 1 + 2, · · · , N. Hence, the same contradiction holds as in Step 1). Therefore, by Lemma 4.2, we get the polynomial decay rate of system (5), that is , ∀t ≥ 0. Note that in Proposition 3, we have proved that the polynomial decay order of the energy of system (5) is at most t −1 . Hence, t − 1 1+ is the nearly sharp. The proof of Theorem 2.4 is complete.

5.
More general slow decay rates. In the last section, we have shown that, when N − N 1 > 1, namely, when the network (5) involves more than one purely elastic rod, the system can not achieve the exponential decay rate. Then we further derived the nearly optimal polynomial decay rate, if the conditions (S) in Definition 2.3 are fulfilled.
In fact, the discussion of the proof of Theorem 2.4 (see subsection 4.3), shows that the slow decay rate of the system is determined by γ n given as (53), which depends on the property of the lengths of the purely elastic rods entering in the network. Hence, other more general slow decay rates can be obtained when different conditions on j , j = N 1 + 1, N 1 + 2, · · · , N are imposed.
To do this, let us introduce the following definition on the irrational sets (see p. 209 in [11]), which is deduced from [7] (see Theorem I, p. 120). 2. Set F: the set of all real numbers ρ such that ρ / ∈ Q and so that its expansion as a continued fraction [0, a 1 , a 2 , · · · , a n , · · · ] is such that (a n ) is bounded. In particular F is contained in the sets B for every > 0. (see p. 209 in [11] for more details). It contains all quadratic algebraic irrational numbers and also some transcendental numbers. This set F has Lebesgue measure zero and is not denumerable.
We have the following result: Corollary 1. For any (u (0) , u (1) , θ (0) ) ∈ D(A), there always exists a constant C > 0 such that the energy of network (5) satisfies where s is given as follows: Since the proof is similar to the one of Theorem 2.4, we only give a sketch of it.
If the decay rate is not fulfilled, there exists a sequence Φ n = (U n , V n , , and a sequence σ n ∈ R with σ n → ∞ such that where s is given as in Corollary 1. By Diophantine approximation (see [11]), different estimates can be gotten for γ n , when the irrational numbers i / j belong to the sets B ε and F (as in as in Definition (5.1), respectively. By Corollary A. 10 in [11] ( see also [39]), we get that (2) if i j ∈ F, i, j = N 1 + 1, N 1 + 2, · · · , N, i = j, then γ n ≥ c/σ N −N1−1 n , n ≥ 1. Thus, proceeding as in the discussion of the proof of Theorem 2.4, finally we can get Φ n → 0, in H, n → ∞. This contradicts Φ n H = 1. Hence, the desired result follows. Remark 6. By Corollary 1, it is easy to see that when N −N 1 = 2, that is there are two purely elastic rods in the network, the system can achieve optimal polynomial decay rate t −1 , if N 1 +1 N ∈ F. More generally, by the similar proof as the one for Theorem 2.4, together with the so called M log −Theorem in [5], we obtain the general slow decay rate of system (5) as follows. If where M −1 log (·) is the inverse function of M log . Especially, if M (s) ≤ ce as , a, c > 0, then the network (5) achieves logarithmic decay rate.
Proof. We can still use the proof of Theorem 2.4 to derive lim sup which together with M log −Theorem in Batty and Duyckaerts [5] yields (61).

Remark 7.
From Corollary 2, we see that in order to obtain an explicit decay rate, it is very important to estimate the lower bound of γ n . Some techniques such as Diophantine approximations can be used to estimate it. However, it is still open to get a sharp estimate for it.
6. Numerical simulations. This section is devoted to present some numerical simulations on the dynamical behaviour of system (5) to support the results obtained above. The backward Euler method in time ( time step: dt = 0.01 ) and the Chebyshev spectral method in space (spatial grid size S = 40) were employed, in a MatLab environment (see [38]).
In these figures, we can see that the behaviour of each u j , j = 1, 2, 3 and θ j , j = 1, 2 are convergent to zero very fast. This shows numerically that the energy of system (5) decays to zero exponentially, which is consistent with the result on exponential decay rate in previous sections. Case D. N 1 = 1 (One thermoelastic rod and two purely elastic rods) In this case, the network consists of one thermoelastic rod and two purely elastic ones. The lengths of rods are still given as 1 = 1, 2 = 2, 3 = 1. Simulations are plotted in Figure D-1, 2, 3, 4. We observe the lack of decay as predicted by the theory. Case E. N 1 = 1, 2 = √ 2 (One thermoelastic rod and two purely elastic ones) In this case, the network is still constituted by one thermoelastic rod and two purely elastic ones. This time the length 2 = √ 2, which leads to 2 3 = √ 2 / ∈ Q. In Figure E-1, 2, 3, 4, we observe a very slow decay rate. It implies that the energy of this system decays to zero but lacks exponential growth rate.
Moreover, we presented Figure F-1, 2 to compare the decay rate of the energy for each case. Figure F-1 shows the dynamical behaviours of the logarithmic scale of the energy for Case A, B and C. Figure F-2 shows the dynamical behaviours of the energies for Case D and E. From these two figures, we can clearly see the behaviours of energies for each case respect time, which are consistent with our theoretical results obtained in this paper.
Appendix A. Proof of (37). This appendix is devoted to show how to get (37).