Martin Gugat
email: martin.gugat@fau.de

Meizhi Qian

Jan Sokolowski
email: jan.sokolowski@univ-lorraine.fr

Network design and control: Shape and topology optimization for the turnpike property for wave equation

Keywords: Turnpike property, wave equation, optimal control, optimality system, spectral method, network optimum design

published or not. The documents may come

Introduction

In this paper, the optimum design problem for networks is considered. The dynamics on the networks are governed by hyperbolic PDEs with boundary or distributed control. The optimal control is uniquely determined by the optimality system for the distributed parameter control problem. The structure of the optimization model is divided into two levels. The lower level is an optimal control problem which delivers an optimal control and the optimal value of the cost. The higher level is a design problem that uses the optimal cost for the shape functional. In other words, the shape functional which is the objective functional on the higher level problem is defined using the optimal value of the lower level problem.

In general, the existence of an optimal shape is a difficult challenge for geometrical domains. However, for one dimensional state equations under reasonable assumptions, the existence of an optimal shape can be shown. The dynamics on the networks are covered by the framework presented in [START_REF] Nicaise | Problèmes de Cauchy posés en norme uniforme sur les espaces ramifiés élémentaires[END_REF]. Thus, we consider bilevel optimization on graphs. The mathematical question that we consider in this paper is the possibility to approximate the dynamic optimal control for evolutionary state equations by the static optimal control obtained for the steady state equation. The approximation error decreases with growing time horizon T for the dynamic problem with fixed initial conditions. This is a significant simplification of the solution strategy when it can be applied. The proposed approximation is useful for large time horizons. The novelty of this paper is the exploitation of the turnpike structure in the framework of optimal shape design.

We consider two optimal control problems. The first problem (OCE) is governed by an evolution equation. In order to define the second problem (OCS), the evolution state equation is reduced to the steady state equation. The optimal controls are given by the appropriate optimality systems, see [START_REF] Lions | Contrôle Optimal de Systèmes Gouvernés Par des équations aux Dérivées partielles[END_REF]. In order to justify the approximation of (OCE) by (OCS), we study the turnpike property of the couple (OCE)-(OCS). For a recent survey on the turnpike property in optimal control, see [START_REF] Grüne | Dissipativity and optimal control: examining the turnpike phenomenon[END_REF] and the references therein. Several forms of the turnpike phenomenon have been studied in detail, for example, the exponential turnpike property (see [START_REF] Trélat | Steady-state and periodic exponential turnpike property for optimal control problems in hilbert spaces[END_REF]) and the interval turnpike property, see [START_REF] Faulwasser | The interval turnpike property for adjoints[END_REF]. Numerical issues and the turnpike phenomenon in optimal shape design with parabolic PDEs has been studied in [START_REF] Lance | Numerical issues and turnpike phenomenon in optimal shape design. In: Optimization and Control for Partial Differential Equations-uncertainty Quantification, Open and Closed-loop Control, and Shape Optimization[END_REF].

Note that that in our paper, we study the application of the turnpike property to shape optimization. This is a novel approach. The turnpike property for control problems in variable domain setting has been considered in [START_REF] Lance | Shape turnpike for linear parabolic PDE models[END_REF]. In this paper, the optimal control determined from the optimality system for (OCS) is used to construct an approach for the determination of the optimal shape that is computationally tractable. Simpler optimum design problems with linear models are considered, e.g., in [START_REF] Klarbring | Topology optimization of flow networks[END_REF]. Our approach is more general.

We specify our approach for the networks, but it could also be used for standard control problems if the existence of an optimal shape is known. In the case of networks the existence of an optimal shape Ω is simpler compared to the case of PDEs defined in geometrical domains in R n for n = 2, 3, see [START_REF] Bucur | Continuity of the elasticity system solutions with respect to the geometrical domain variations[END_REF] for the elasticity system.

The outline of the paper is the following. We introduce two optimal control problems (OCE) and (OCS) with the controls u = u(t, x) and v = v(x) and the optimal controls û and v, respectively. The geometrical domain of the network is designated by Ω. Thus û = û(T, Ω; t, x) and v = v(Ω; x) for the spatial variable x ∈ Ω and the time horizon T for (OCE). The cost functional for (OCE) takes the form

u -→ 1 2 T 0 Ω (y(u) -y d ) 2 + γ (y t (u)) 2 dx dt + 1 2 T 0 Γ (u -u d ) 2 dΓ(x)dt (1)
where y(u) denotes the state that is generated by the control u. For systems that are governed by the wave equation, a regularization parameter γ > 0 should be chosen, since we need the velocity part in the cost. In the parabolic case, we do not need such a term hence γ = 0 can be chosen.

Our aim is to study the shape optimization problem for the modified shape functional

Ω -→ 1 2 T 0 Ω (y(v) -y d ) 2 + γ (y t (v)) 2 dx dt + 1 2 T 0 Γ (v -u d ) 2 dΓ(x)dt. (2)
In this approach, we replace the optimal control for (OCE) by the optimal control for (OCS) and as a result we get the shape functional [START_REF] Lions | Contrôle Optimal de Systèmes Gouvernés Par des équations aux Dérivées partielles[END_REF] for the optimum design of the network. Such a procedure will be justified by the turnpike property of the couple (OCE) and (OCS).

We want to mention some recent contributions about shape optimization problems from the literature. While the time-dependent shape optimization problem for the heat equation has already been studied in the literature (see e.g. [START_REF] Brügger | On the numerical solution of a timedependent shape optimization problem for the heat equation[END_REF]), fewer results are known for the hyperbolic case (see e.g. [START_REF] Cagnol | Hidden shape derivative in the wave equation[END_REF]). In [START_REF] Schlegel | Shape optimization for the mitigation of coastal erosion via porous shallow water equations[END_REF] the problem of mitigation of coastal erosion is studied where the water flow is modeled by the porous shallow water equations.

Optimal control problems on graphs

Let us consider the finite connected graph G = (V, E) with V = {P j |j ∈ J } the set of vertices and E = {E i |i ∈ I} the set of edges. We consider the evolution equations on the edges of G. The dynamics of the edges are coupled by suitable node conditions for example Kirchhoff-type. Therefore we assume that the edges correspond to curves Ω i in R 1 , R 2 , or R 3 .

Define Ω = ∪ i∈I Ω i . We are interested in control problems defined on G. The controls are denoted by v := v(x) with x ∈ Ω, and by u T := u T (x, t) with x ∈ Ω and t ∈ (0, T ). There are two control problems, the first is the static control problem, and the second is the evolution control problem. The static state equation gives the state z := z(v; x) = z(x) that is determined by the state equation for

z ∈ H : a(z, ϕ) = (L(v), ϕ) + (f, ϕ) ∀ϕ ∈ H , (3)
where x → z(v; x) lives in the Hilbert space H. In the applications usually H is a Sobolev space. We assume that the bilinear form a(•, •) is symmetric and coercive with a bounded linear operator L from the space of controls to H and f ∈ H.

As an example, consider the evolution equation the governs the state denoted by y := y(u; x, t) = y(x, t) in variational form for t → y(u; x, t):

∂ 2 y ∂ 2 t
(t), φ + a(y(t), φ) = (L(u)(t), φ) + (F (t), φ)

for all φ ∈ H and a.e. for t ∈ (0, T ), along with the initial conditions y(0, x) = y 0 (x) and y t (0, x) = y 1 (x). Here H is the Hilbert space that contains the functions with H 1 -regularity on the edges that are compatible with the node conditions that are prescribed on the vertices. Moreover, F (t) = f for all t ∈ (0, T ).

In both cases, the controls are introduced as the distributed controls on edges and the boundary controls on vertices. Remark 1 In order to solve the state equation (4) we will use the spectral method under some assumptions. There are two different types of networks under consideration. The first graph is a tree and the shape optimization with respect to the length of edges is analyzed. In such a case the first eigenvalue of the elliptic-steady state operator is supposed to be strictly positive uniformly on the set of admissible shapes. The set of admissible shapes for the tree is compact in the finite dimensional space R N . The second graph includes cycles. Remark 2 For the purposes of topology optimization a small cycle replaces a central node P 0 in the tree, the central node is selected by using the topological derivative of the cost, see Figures 11 and12.

Let us consider for the static problem the Laplacian -∂ xx on the edges of the graph. Under appropriate geometric assumptions the eigenvalues supported by a small cycle of size ε → zero are given e.g., by µ n = nΠ ε The fundamental question is if the spectral proof of turnpike property holds in such a case. We need some regularity of solutions to the wave equation uniform with respect to the size of the cycle. The required regularity assumption on the initial conditions of the wave equation assures the appropriate regularity of solutions with respect to time variable t. Such an assumption implies that the constant displacement on the cycle is not allowed, however, the constant velocity is allowed, see (53) for details.

Indeed, let us consider the constant displacement and denote by E(ε) the energy associated with the initial conditions in function of the size of the cycle. Such an energy associated with the constant initial state y(0) = 1 at the cycle can be evaluated using (53), where the coefficients of the expansion over the eigenfunctions read

a ε k = ε 0 ψ k (x)dx = ε 0 sin kΠx ε dx = - ε kΠ cos kΠx ε ε 0
Therefore, the energy in (53) is unbounded, thus the constant initial state is not allowed for the cycle. Similarly, the positive result is shown for the constant velocity. Remark 3 The singular perturbations of the geometry for networks are of different nature for the static state equations and for the wave equations. Therefore, we have constructive results in the case of elliptic problems which are combined with the turnpike property in order to design the topology of the network. This is a new result to our best knowledge.

We need as well the cost functional, the simplest possibility is the quadratic cost with the appropriate choice of norms in Hilbert function spaces. For the static problem, it is

I(v) = 1 2 ∥z -z d ∥ 2 L 2 (Ω) + 1 2 ∥v -v d ∥ 2 L 2 (Γ) . (5)
For the evolution problem, it is

J T (u) = 1 2 T 0 ∥y -y d ∥ 2 L 2 (Ω) + γ ∥∂ t (y -y d )∥ 2 L 2 (Ω) dt + 1 2 T 0 ∥u -u d ∥ 2 L 2 (Γ) dt (6)
where γ ∈ (0, ∞). For the sake of simplicity, we consider the control problems without constraints.

The optimality conditions are necessary and sufficient and can be obtained by using the Lagrangian formalism. To this end, we need the adjoint states, which simplify the formulas for the gradients of the cost with respect to controls. We ask two questions now: the first is the shape and topology optimization for the graph with respect to the control problem, this results in the bi-level optimization problem for the graph. The second question is the turnpike property for the two-level optimization problem for the graph. To our best knowledge, such a problem has not been considered in the literature.

Optimality system for the static problem

We denote by (v, ẑ, p) the unique optimal solution of static control problem. Let us note that the uniqueness of the optimal control û follows by convexity of variational problem under consideration with quadratic cost functional and the linear state equation. The optimal solution is given by the optimality system which depends on the shape, or design Ω of the network as an infinite dimensional factor or parameter to be selected at the upper level of optimization over the class of admissible shapes. The class of admissible shapes is denoted by S ad , and the continuous variation of the shape is denoted by Ω τ ∈ S ad with the real parameter τ ∈ (-δ, δ) for some δ > 0 for shape variations. We use as well the singular perturbations of the shape denoted by Ω ϵ ∈ S ad for the topology variations with ϵ → 0. Remark 4 In the case of a network, there are at least two possibilities of shape variations. The first is the change of the lengths of edges, it corresponds to the boundary variations in the classical shape optimization. The second, which corresponds to the topology variations in the shape optimization means the presence of a small cycle within the network with the small size of the cycle ϵ → 0.

The optimality system for the static control problem is equivalent to the vanishing of the gradient for the cost, hence

min v {J(v)} = J(v) iff the following optimality system is verified ẑ ∈ H : a(ẑ, φ) = (L(v), φ) + (f, φ) ∀φ ∈ H , (7)
p ∈ H :

a(p, ϕ) = (z d -ẑ, ϕ) ∀ϕ ∈ H , (8)
(v d -v, v) U = (L ′ (v) • v, p) ∀v ∈ U . (9)
In the linear case we have

L ′ (v) • v = L(v).
The optimality system is derived using the Lagrangian formalism,

L(v, z, ϕ) = 1 2 ∥z -z d ∥ 2 L 2 (Ω) + 1 2 ∥v -v d ∥ 2 L 2 (Γ) + a(z, ϕ) -(L(v), ϕ) -(f, ϕ)
Then the adjoint state p ∈ H is introduced

∂ z L(v, z, p)(ϕ) = (z -z d , ϕ) Ω + a(p, ϕ) = 0 ∀ϕ ∈ H
and the gradient of the cost is obtained

dI(v; η) = ∂ v L(v, z, p)(η) = (v -v d , η) Γ -(p, η) Γ
which leads to the optimality condition v -v d = p a.e. on Γ.

In the case of distributed control, we have Γ = Ω. Proposition There exists the unique solution (v, ẑ, p) to the optimality system (7)- [START_REF] Bucur | Continuity of the elasticity system solutions with respect to the geometrical domain variations[END_REF]. The optimal value of the cost J(v) := J(v(Ω)) is defined as a shape functional over the set S ad . Therefore, we consider the optimum design of the network inf

Ω∈S ad J(v(Ω)) . (10)
The analysis of such a variational problem requires:

1. The existence of solutions; 2. The necessary optimality conditions; 3. Finally, numerical methods for solution.

In particular, we perform the shape calculus and determine the shape gradient of the optimal control cost Ω -→ J(v(Ω))

as well as the topological derivative obtained at ϵ := 0 + for the mapping

ϵ -→ J(v(Ω ϵ )).
Remark 5 It is useful for applications to introduce a random right-hand side f := f (ω; x) to the state equation (3).

Optimality system for the evolution problem

We derive the optimality system for the optimal control problem with the evolution state equation. To this end, we introduce the Lagrangian

L(u, y, φ) = 1 2 T 0 ∥y -y d ∥ 2 L 2 (Ω) dt + 1 2 T 0 ∥u -u d ∥ 2 L 2 (Γ) dt + γ 2 T 0 ∥∂ t (y -y d )∥ 2 L 2 (Ω) dt + T 0 ∂ 2 y ∂ 2 t (t), φ L 2 (Ω) dt + T 0 a(y(t), φ)dt - T 0 (L(u)(t), φ) L 2 (Γ) dt - T 0 (F (t), φ) L 2 (Ω) dt.
We recall the integration by parts formula

T 0 ∂ 2 y ∂ 2 t (t), φ L 2 (Ω) dt := T 0 (y tt , φ) Ω dt = T 0 (y, φ tt ) Ω dt + (y 1 , φ(0)) Ω -(y t (T ), φ(T )) Ω -(y 0 , φ t (0)) Ω + (y(T ), φ t (T )) Ω .
The adjoint state p = p(u; x, t) is obtained by differentiation of Lagrangian with respect to the state, thus with

Q(T ) = (0, T ) × Ω (p tt , φ) Q(T ) + T 0 a(p, φ) dt = (y d -y, φ) Q(T ) + γ (∂ t (y d -y), φ t ) Q(T ) ∀φ ∈ H(Q(T )) p(T ) = p t (T ) = 0 .
The gradient of the cost with respect to the control is obtained by differentiation of the Lagrangian with respect to the control

dJ du (u), η = T 0 (u -u d , η) L 2 (Γ) dt - T 0 (L(η)(t), p(t)) L 2 (Γ) dt = T 0 (u -u d , η) Γ dt - T 0 (L(η)(t), p(t)) Γ dt.
The following lemma contains the necessary optimality conditions for the dynamic problem where J T (u) as defined in ( 6) is minimized. Lemma 1 The optimality system for the optimal control of evolution control problem is verified for a.e. t ∈ (0, T ):

(ŷ tt , φ) Q(T ) + T 0 a(ŷ(t), φ) dt = T 0 (L(û)(t), φ) Γ dt + (F (t), φ) Q(T ) (11)
∀φ ∈ H(Q(T ))

ŷ(0) = y 0 , ŷt (0) = y 1 (12) (p tt , φ) H(Q(T )) + T 0 a(p, φ) dt = (y d -ŷ + γ ŷtt , φ) Q(T ) ( 13
)
∀φ ∈ H(Q(T )) p(T ) = pt (T ) = 0 ( 14
) (û -u d , η) Γ -(L(η)(t), p(t)) Γ = 0 ∀η ∈ L 2 (Σ) (15)
The optimality system ( 11)-( 15) admits a unique solution (û, ŷ, p).

5 The difference of the static and the dynamic optimality systems for distributed control

In this section, we study a system that is satisfied by the ordered pair that has the difference between the optimal dynamic state for the time horizon T and the optimal static state as the first component and the difference between the optimal dynamic adjoint state for the time horizon T and the optimal static adjoint state as the second component. The question of long time versus steady state optimal control has already been studied in [START_REF] Porretta | Long time versus steady state optimal control[END_REF] where the focus is on the turnpike property of the state and the control without the adjoint state. The turnpike phenomenon for optimal boundary control problems with first order hyperbolic systems is considered in [START_REF] Gugat | On the turnpike phenomenon for optimal boundary control problems with hyperbolic systems[END_REF]. We assume that Γ = Ω and the operators L(v) and L(u) are identity operators in L 2 (Ω) and F (t) = f , u d (t) = v d , y d = z d for t ∈ [0, T ] almost everywhere. The optimality system for the static problem reads

ẑσ ∈ H : a(ẑ σ , φ) = (v σ , φ) L 2 (Ω) + (f, φ) L 2 (Ω) ∀φ ∈ H, (16)
pσ

∈ H : a(p σ , ϕ) = (z d -ẑσ , ϕ) L 2 (Ω) ∀ϕ ∈ H, (17)
vσ -v d = pσ a.e. in Ω.

The optimality system for the evolution problem implies for t ∈ [0, T ] almost everywhere (with

F (t) = f ) (ŷ T tt (t), φ) L 2 (Ω) + a(ŷ T (t), φ) = (û T (t), φ) L 2 (Ω) + (F (t), φ) L 2 (Ω) ∀φ ∈ H (19) ŷT (0) = y 0 , ŷT t (0) = y 1 (20) (p T tt (t), φ) L 2 (Ω) + a(p T (t), φ) = (y d -ŷT (t) + γ ŷtt , φ) L 2 (Ω) ∀φ ∈ H (21)
pT (T ) = pT t (T ) = 0 [START_REF] Gugat | Optimal Boundary Control and Boundary Stabilization of Hyperbolic Systems[END_REF] ûT (t) -u d = pT (t) a.e. in Ω × (0, T ). [START_REF] Gugat | Proceedings of the 27th International Conference on Methods and Models in Automation and Robotics (MMAR)[END_REF] Define the differences

ω T = ŷT -ẑσ , µ T = pT -pσ , ν T = ûT -vσ .
Then for all φ ∈ H we have the initial condition

ω T (0) = y 0 -ẑσ , ω T t (0) = y 1 , (24)
the terminal conditions µ T (T ) = -p σ , µ T t (T ) = 0, (25) the dynamics

(ω T tt (t), φ) L 2 (Ω) + a(ω T (t), φ) = (ν T (t), φ) L 2 (Ω) = (µ T (t), φ) L 2 (Ω) (26)
and with the assumption that y d = z d and

u d = v d (µ T tt (t), φ) L 2 (Ω) + a(µ T (t), φ) = -(ω T (t), φ) L 2 (Ω) + γ (ω T tt (t), φ) L 2 (Ω) . (27)
Note that for the difference system, the existence of a solution is implied by the construction as the difference between two systems, for which solutions exist.

Note that for the energy

E(t) = (ω T t (t), ω T t (t)) L 2 (Ω) + a(ω T (t), ω T (t)) 2 + (µ T t (t), µ T t (t)) L 2 (Ω) + a(µ T (t), µ T (t)) 2
due to [START_REF] Plotnikov | Geometric aspects of shape optimization[END_REF] and [START_REF] Chen | Animal shapes, modal analysis, and visualization of motion (I): Horse and camel[END_REF] for γ = 0 we have

E ′′ (t) ≥ (ω T (t), ω T (t)) L 2 (Ω) + (µ T (t), µ T (t)) L 2 (Ω) ≥ 0.
Thus E is convex on [0, T ]. Now we perform a spectral analysis to show the exponential turnpike property.

Assume that there exists a complete orthonormal sequence (ψ k ) ∞ k=1 of eigenfunctions with a(ψ k , φ) = λ k (ψ k , φ) L 2 (Ω) for all k ∈ {0, 1, 2, 3, ..} where

λ k ≥ γ > 0 ( 28
)
is a real number. Remark 6 In the case of optimal design, the bilinear form depends on Ω. In the case of a graph, this means the dependence of the lengths of the edges. In this case, the eigenvalues and eigenfunctions depend on these parameters that we denote by ℓ.

Therefore a meaningful analysis has to take into account the sensitivity with respect to ℓ.

Our assumption on the feasible designs is that the smallest eigenvalues is greater than or equal to the given strictly positive lower-bound γ > 0 uniformly on the set of admissible designs. In this way, we ensure that the turnpike property is valid for the bilevel optimization problem that we consider in this paper.

It is well known that the smallest eigenvalue that can be characterized as the Rayleigh quotient depends smoothly on the parameters, see [START_REF] Henrot | Extremum Problems for Eigenvalues of Elliptic Operators[END_REF]. In our analysis, the particular structure of the spectrum is not relevant.

Then we can use the representations

ω T = ∞ k=0 a k (t) ψ k (x), µ T = ∞ k=0 b k (t) ψ k (x) (29)
to show that ω T and µ T have the turnpike property.

It is clear that the functions a k and b k depend on T as a parameter, so a more precise notation would be a k,T (t) and b k,T (t). However, in order to make the text more concise we continue with the shorter notation a k and b k .

From ( 26) and ( 27) we obtain

a ′′ k = -λ k a k + b k and b ′′ k = -λ k b k -a k + γ a ′′ k .
Thus we have the sequence of differential equations

a (4) k + (2 λ k -γ) a (2) k + (λ 2 k + 1) a k = 0, b (4)
k + (2 λ k -γ) b (2) k + (λ 2 k + 1) b k = 0. (30)
Define the characteristic polynomial

p k (z) = z 4 + (2 λ k -γ) z 2 + (λ 2 k + 1).
Since p k is a polynomial in z 2 , for the roots z (l)

k of p k we obtain (z (l) k ) 2 = γ 2 -λ k ± i 2 4 + 4λ k γ -γ 2 for l ∈ {1, 2, 3, 4} which implies (z (l) k ) 2 = 1 + λ 2 k .
Thus there are two pairs of complex conjugate roots and we have the representation

p k (z) = (z -z (1) k )(z -z (1) k )(z + z (1) k )(z + z (1) k ). We have |z (l) k | 4 = λ 2 k + 1, Re((z (l) ) 2 k ) = γ 2 -λ k and |Im((z (l) ) 2 k )| = √ 4+4 λj γ-γ 2 2
. Moreover, we have

|Re(z (l) k )| = γ 4 + 1 2(λ k + 1 + λ 2 k )
and |Re(z

(l) k )| ≥ √ γ 2 . ( 31
)
The initial condition [START_REF] Soko Lowski | Modelling of topological derivatives for contact problems[END_REF] yields the values for a k (0) and a ′ k (0). The terminal condition (25) yields the value for b k (T ) and b ′ k (T ) = 0. Note that [START_REF] Novotny | An Introduction to the Topological Derivative Method[END_REF] implies that the value of b k (T ) is also independent of T , just like b ′ k (T ). For the sake of conciseness, in the sequel we use the notation

z k = z (1) k .

Representation of the solution

The solution of the optimality system means that for the coefficients b k,T of µ T as defined in (29) we solve a boundary value problem with the ODE of order four (30) i.e.

           b (4) k,T + (2 λ k -γ) b (2) k,T + (λ 2 k + 1) b k,T = 0 b k,T (T ) = β k b ′ k,T (T ) = 0 -b ′′ k,T (0) + (γ -λ k ) b k,T (0) = (1 + γ λ k ) a k,T (0) -b ′′′ k,T (0) + (γ -λ k ) b ′ k,T (0) = (1 + γ λ k ) a ′ k,T (0). ( 32
)
where the value of β k is determined by the terminal condition µ T (T ) = -p σ in [START_REF] Novotny | An Introduction to the Topological Derivative Method[END_REF].

We represent the solution in the form

b k,T (t) = F k,T (t) β k + (1 + γ λ k )[G k,T (t) a k,T (0) + H k,T (t) a ′ k,T (0)]. ( 33
)
The following Lemma contains explicit representations of F k,T , G k,T and H k,T . In the representation, the numbers d(k, T ) appear as multipliers, therefore it is important that for T sufficiently large we have d(k, T ) ̸ = 0. Since in the study of the turnpike phenomenon we are interested in large time horizons, the assumption that the time horizon T is large is not restrictive for us. We introduce the notation

Ξ k := γ -λ k -z 2 k = γ 2 ∓ i 2 4 + 4 λ k γ -γ 2 . ( 34
)
Lemma 2 For k ∈ {0, 1, 2, ...} and T sufficiently large define the numbers

d(k, T ) = -2Re Ξ 2 k + 2|Ξ k | 2 Re cosh(z k T ) cosh(z k T ) -z k z k sinh(z k T ) sinh(z k T ) . (35)
Then we have

d(k, T ) F k,T (t) = 2Re -Ξ k 2 + |Ξ k | 2 cosh(z k T ) cosh(z k T ) cosh(z k (t -T )) + 2Re |Ξ k | 2 cosh(z k T ) sinh(z k T ) sinh(z k (t -T )) -2|Ξ k | 2 Re z k z k sinh(z k T ) [sinh(z k T ) cosh(z k (t -T )) + cosh(z k T ) sinh(z k (t -T ))] .
Furthermore, it holds

d(k, T ) G k,T (t) = -2Re (Ξ k [cosh(z k T ) cosh(z k (t -T )) + sinh(z k T ) sinh(z k (t -T ))]) + 2Re Ξ k cosh(z k T ) cosh(z k (t -T )) + z k z k sinh(z k T ) sinh(z k (t -T )) (36)
and

d(k, T ) H k,T (t) = 2Re Ξ k z k sinh(z k T ) cosh(z k (t -T )) + Ξ k z k cosh(z k T ) sinh(z k (t -T )) -2Re Ξ k z k [sinh(z k T ) cosh(z k (t -T )) + cosh(z k T ) sinh(z k (t -T ))] .
Proof: Due to the properties of the hyperbolic functions that can be found for example in [START_REF] Bronstein | Handbook of Mathematics[END_REF][START_REF] Bronstein | Taschenbuch der Mathematik. Übers. aus dem Russischen. Mit 1 CD-ROM[END_REF] we have

d(k, T ) = -2 Re(Ξ 2 k ) (37) +|Ξ k | 2     1 + 1 -γ 2λ k 1 + 1 λ 2 k   cosh(2 Re(z k ) T ) +   1 - 1 -γ 2λ k 1 + 1 λ 2 k   cos(2 Im(z k ) T )   .
This implies lim T →∞ |d(k, T )| = ∞, hence for all sufficiently large T we have d(k, T ) ̸ = 0. In the remaining part of the proof, we assume that d(k, T ) ̸ = 0. Moreover, we see that d(k, T ) grows exponentially fast with T with the growth rate

2Re(z k ) ≥ √ γ (38)
due to (31).

To be more precise, note that Re(

Ξ 2 k ) = γ 2 2 -1 + 5 2 γ λ k + λ 2 k grows with k with the order λ 2 k and |Ξ k | 2 = 1 + γ λ k
grows with the order λ k . At this point, elementary computations show that F k,T (t), G k,T (t) and H k,T (t) satisfy the differential equation and the boundary conditions where (β k , a k,T (0), a ′ k,T (0)) have the values (1, 0, 0), (0, 1, 0), (0, 0, 1) respectively. Now we present details of the verification of these basic functions.

Verification of the basis function F k,T (t)

We start with the function F k,T (t) from (33) that is multiplied with β k , i.e. d(k, T ) F k,T (t). We have

d(k, T ) F k,T (t) = 2Re -Ξ k 2 + |Ξ k | 2 cosh(z k T ) cosh(z k T ) cosh(z k (t -T )) + 2Re |Ξ k | 2 cosh(z k T ) sinh(z k T ) sinh(z k (t -T )) -2Re z k z k |Ξ k | 2 sinh(z k T ) sinh(z k T ) cosh(z k (t -T )) -2Re z k z k |Ξ k | 2 cosh(z k T ) sinh(z k T ) sinh(z k (t -T )) .
Hence for t = T we have

d(k, T ) F k,T (T ) = 2Re -Ξ k 2 + |Ξ k | 2 cosh(z k T ) cosh(z k T ) = -2Re Ξ 2 k + 2|Ξ k | 2 Re cosh(z k T ) cosh(z k T )) - z k z k sinh(z k T ) sinh(z k T )) = d(k, T ). Since d(k, T ) ̸ = 0, this implies F k,T (T ) = 1.
For the derivative we obtain

d(k, T ) F ′ k,T (t) = 2Re z k -Ξ k 2 + |Ξ k | 2 [cosh(z k T ) cosh(z k T )] sinh(z k (t -T )) + 2Re z k |Ξ k | 2 cosh(z k T ) sinh(z k T ) cosh(z k (t -T )) -2Re z k |Ξ k | 2 sinh(z k T ) sinh(z k T ) sinh(z k (t -T )) -2Re z k |Ξ k | 2 cosh(z k T ) sinh(z k T ) cosh(z k (t -T )) .
Thus we have

d(k, T ) F ′ k,T (T ) = 2Re z k |Ξ k | 2 cosh(z k T ) sinh(z k T ) -2Re z k |Ξ k | 2 cosh(z k T ) sinh(z k T ) = 0.
Therefore we have F ′ k,T (T ) = 0. For t = 0 we have

d(k, T ) F k,T (0) = 2Re -Ξ k 2 + |Ξ k | 2 [cosh(z k T ) cosh(z k T )] cosh(z k T ) -2Re |Ξ k | 2 cosh(z k T ) sinh(z k T ) sinh(z k T ) = 2Re -Ξ k 2 cosh(z k T ) + |Ξ k | 2 cosh(z k T ) cosh(z k T ) 2 -sinh 2 (z k T ) = 2Re -Ξ k 2 cosh(z k T ) + |Ξ k | 2 cosh(z k T ) .
For the second derivative, we have

d(k, T ) F ′′ k,T (t) = 2Re z 2 k -Ξ k 2 + |Ξ k | 2 [cosh(z k T ) cosh(z k T )] cosh(z k (t -T )) + 2Re z 2 k |Ξ k | 2 cosh(z k T ) sinh(z k T ) sinh(z k (t -T )) -2Re z k z k |Ξ k | 2 sinh(z k T ) sinh(z k T ) cosh(z k (t -T )) -2Re z k z k |Ξ k | 2 cosh(z k T ) sinh(z k T ) sinh(z k (t -T )) . For t = 0 this yields d(k, T ) F ′′ k,T (0) = 2Re z 2 k -Ξ k 2 + |Ξ k | 2 cosh(z k T ) cosh(z k T ) cosh(z k T ) -2Re z 2 k |Ξ k | 2 cosh(z k T ) sinh(z k T ) sinh(z k T ) = 2Re z 2 k -Ξ k 2 cosh(z k T ) + |Ξ k | 2 cosh(z k T ) .
Hence we have

d(k, T ) -F ′′ k,T (0) + (γ -λ k ) F k,T (0) = 2Re Ξ k -Ξ k 2 cosh(z k T ) + |Ξ k | 2 cosh(z k T ) = 2|Ξ k | 2 Re -Ξ k cosh(z k T ) + Ξ k cosh(z k T ) = 0. Thus we have -F ′′ k,T (0) + (γ -λ k ) F k,T (0) = 0. To proceed, let us observe that d(k, T ) F ′ k,T (0) = 2Re z k Ξ k 2 -|Ξ k | 2 cosh(z k T ) cosh(z k T ) sinh(z k T ) + 2Re z k |Ξ k | 2 cosh(z k T ) sinh(z k T ) cosh(z k T ) + 2Re z k |Ξ k | 2 sinh(z k T ) sinh(z k T ) sinh(z k T ) -2Re z k |Ξ k | 2 cosh(z k T ) sinh(z k T ) cosh(z k T ) = 2Re z k Ξ k 2 sinh(z k T ) -2Re z k |Ξ k | 2 sinh(z k T ) .
For the third derivative, we have

d(k, T ) F ′′′ k,T (t) = 2Re z 3 k -Ξ k 2 + |Ξ k | 2 cosh(z k T ) cosh(z k T ) sinh(z k (t -T )) + 2Re z 3 k |Ξ k | 2 cosh(z k T ) sinh(z k T ) cosh(z k (t -T )) -2Re z 2 k z k |Ξ k | 2 sinh(z k T ) sinh(z k T ) sinh(z k (t -T )) -2Re z 2 k z k |Ξ k | 2 cosh(z k T ) sinh(z k T ) cosh(z k (t -T )) .
This yields

d(k, T ) F ′′′ k,T (0) = -2Re z 3 k -Ξ k 2 + |Ξ k | 2 cosh(z k T ) cosh(z k T ) sinh(z k T ) + 2Re z 3 k |Ξ k | 2 cosh(z k T ) sinh(z k T ) cosh(z k T ) -2Re z 2 k z k |Ξ k | 2 sinh(z k T ) = 2Re z 3 k Ξ k 2 sinh(z k T ) -2Re z 2 k z k |Ξ k | 2 sinh(z k T ) .
Hence we have

d(k, T ) -F ′′′ k,T (0) + (γ -λ k ) F ′ k,T (0) = -2Re z 3 k Ξ k 2 sinh(z k T ) + 2Re z 2 k z k |Ξ k | 2 sinh(z k T ) + (γ -λ k ) 2Re z k Ξ k 2 sinh(z k T ) -(γ -λ k ) 2Re z k |Ξ k | 2 sinh(z k T ) = 2Re z k Ξ k Ξ k 2 sinh(z k T ) -2Re Ξ k z k |Ξ k | 2 sinh(z k T ) = 0.
Thus we have shown that

F k,T (T ) = 1, F ′ k,T (T ) = 0, -F ′′ k,T (0) + (γ -λ k ) F k,T (0) = 0, -F ′′′ k,T (0) + (γ -λ k ) F ′ k,T (0) = 0.
Verification of the basis function G k,T (t).

We continue with G k,T (t) from (33) that is multiplied with

(1 + γ λ k ) a k,T (0) , i.e. d(k, T ) G k,T (t) = -2 Re (Ξ k cosh(z k T ) [1 + tanh(z k T ) tanh(z k (t -T ))] cosh(z k (t -T ))) + 2 Re Ξ k cosh(z k T ) 1 + z k z k tanh(z k T ) tanh(z k (t -T )) cosh(z k (t -T )) .
Due to the definition of z k , the function G k,T (t) satisfies the ODE (30). Moreover, for t = T we have

d(k, T ) G k,T (T ) = -2Re (Ξ k cosh(z k T )) + 2Re Ξ k cosh(z k T ) = 0.
In addition, for the derivative we obtain

d(k, T ) G ′ k,T (t) = -2 Re (Ξ k z k [cosh(z k T ) sinh(z k (t -T )) + sinh(z k T ) cosh(z k (t -T ))]) + 2 Re Ξ k [z k cosh(z k T ) sinh(z k (t -T )) + z k sinh(z k T ) cosh(z k (t -T ))] .
Hence for t = T we have the derivative

d(k, T ) G ′ k,T (T ) = -2Re (Ξ k z k sinh(z k T )) + 2Re Ξ k z k sinh(z k T ) = 0.
For the second derivative, we obtain

d(k, T ) G ′′ k,T (t) = -2 Re Ξ k z 2 k [cosh(z k T ) cosh(z k (t -T )) + sinh(z k T ) sinh(z k (t -T ))] + 2 Re Ξ k z k [z k cosh(z k T ) cosh(z k (t -T )) + z k sinh(z k T ) sinh(z k (t -T ))] . For t = 0 this implies d(k, T ) G ′′ k,T (0) = -2 Re Ξ k z 2 k cosh 2 (z k T ) -sinh 2 (z k T ) + 2 Re Ξ k z k [z k cosh(z k T ) cosh(z k T ) -z k sinh(z k T ) sinh(z k T )] = -2 Re Ξ k z 2 k + 2 Re Ξ k z k [z k cosh(z k T ) cosh(z k T ) -z k sinh(z k T ) sinh(z k T )] .
For t = 0, equation (36) implies

d(k, T ) G k,T (0) = -2 Re (Ξ k ) + 2 Re Ξ k cosh(z k T ) cosh(z k T ) - z k z k sinh(z k T ) sinh(z k T ) .
Hence in view of (34) we obtain

d(k, T ) -G ′′ k,T (0) + (γ -λ k ) G k,T (0) = 2Re Ξ k z 2 k -2Re Ξ k z k [z k cosh(z k T ) cosh(z k T ) -z k sinh(z k T ) sinh(z k T )] -(γ -λ k )2Re (Ξ k ) + (γ -λ k )2Re Ξ k cosh(z k T ) cosh(z k T ) - z k z k sinh(z k T ) sinh(z k T ) = -2Re Ξ 2 k + 2Re Ξ k Ξ k cosh(z k T ) cosh(z k T ) - z k z k Ξ k Ξ k sinh(z k T ) sinh(z k T ) = -2Re Ξ 2 k + 2|Ξ k | 2 Re cosh(z k T ) cosh(z k T ) - z k z k sinh(z k T ) sinh(z k T ) = d(k, T ).
due to the definition (35). Since d(k, T ) ̸ = 0, this yields

-G ′′ k,T (0) + (γ -λ k ) G k,T (0) = 1.
For t = 0, the derivative satisfies the equation

d(k, T ) G ′ k,T (0) = -2 Re (Ξ k z k [-cosh(z k T ) sinh(z k T ) + sinh(z k T ) cosh(z k T )]) + 2 Re Ξ k [-z k cosh(z k T ) sinh(z k T ) + z k sinh(z k T ) cosh(z k T )] = 2 Re Ξ k [-z k cosh(z k T ) sinh(z k T ) + z k sinh(z k T ) cosh(z k T )] .
For the third derivative, we obtain

d(k, T ) G ′′′ k,T (t) = -2 Re Ξ k z 3 k [cosh(z k T ) sinh(z k (t -T )) + sinh(z k T ) cosh(z k (t -T ))] + 2 Re Ξ k z 2 k [z k cosh(z k T ) sinh(z k (t -T )) + z k sinh(z k T ) cosh(z k (t -T ))] .
Hence for t = 0 we have

d(k, T ) G ′′′ k,T (0) = -2 Re Ξ k z 3 k [-cosh(z k T ) sinh(z k T ) + sinh(z k T ) cosh(z k T )] + 2 Re Ξ k z 2 k [-z k cosh(z k T ) sinh(z k T ) + z k sinh(z k T ) cosh(z k T )] = 2 Re Ξ k z 2 k [-z k cosh(z k T ) sinh(z k T ) + z k sinh(z k T ) cosh(z k T )] .
In order to verify the fourth boundary condition in (32), in view of (34) we obtain the equation

d(k, T ) -G ′′′ k,T (0) + (γ -λ k ) G ′ k,T (0) = -2 Re Ξ k z 2 k [-z k cosh(z k T ) sinh(z k T ) + z k sinh(z k T ) cosh(z k T )] + (γ -λ k )2 Re Ξ k [-z k cosh(z k T ) sinh(z k T ) + z k sinh(z k T ) cosh(z k T )] = -2 Re Ξ k Ξ k [z k cosh(z k T ) sinh(z k T ) -z k sinh(z k T ) cosh(z k T )] = -2 |Ξ k | 2 Re ([z k cosh(z k T ) sinh(z k T ) -z k sinh(z k T ) cosh(z k T )]) = 0.
Thus we have shown that

G k,T (T ) = 0, G ′ k,T (T ) = 0, -G ′′ k,T (0) + (γ -λ k ) G k,T (0) = 1, -G ′′′ k,T (0) + (γ -λ k ) G ′ k,T (0) = 0.
Verification of the basis function H k,T (t).

Now we consider H k,T (t) from (33) that is multiplied with

(1 + γ λ k ) a ′ k,T (0) , i.e. d(k, T ) H k,T (t) = 2Re Ξ k z k sinh(z k T ) cosh(z k (t -T )) + Ξ k z k cosh(z k T ) sinh(z k (t -T )) -2Re Ξ k z k [sinh(z k T ) cosh(z k (t -T )) + cosh(z k T ) sinh(z k (t -T ))] .
Due to the definition of z k , the function H k,T (t) satisfies the ODE (30). Moreover, for t = T we have

d(k, T ) H k,T (T ) = 2Re Ξ k z k sinh(z k T ) -2Re Ξ k z k sinh(z k T ) = 0.
For the derivative we have

d(k, T ) H ′ k,T (t) = 2Re Ξ k z k z k sinh(z k T ) sinh(z k (t -T )) + Ξ k cosh(z k T ) cosh(z k (t -T )) -2Re (Ξ k [sinh(z k T ) sinh(z k (t -T )) + cosh(z k T ) cosh(z k (t -T ))]) .
Hence for t = T we obtain H ′ k,T (T ) = 0. For the second derivative, we have

d(k, T ) H ′′ k,T (t) = 2Re Ξ k z 2 k z k sinh(z k T ) cosh(z k (t -T )) + Ξ k z k cosh(z k T ) sinh(z k (t -T )) -2Re (Ξ k z k [sinh(z k T ) cosh(z k (t -T )) + cosh(z k T ) sinh(z k (t -T ))]) .
For t = 0 this implies

d(k, T ) H ′′ k,T (0) = 2Re Ξ k z 2 k z k sinh(z k T ) cosh(z k T ) -Ξ k z k cosh(z k T ) sinh(z k T ) .
We have

d(k, T ) H k,T (0) = 2Re Ξ k z k sinh(z k T ) cosh(z k T ) - Ξ k z k cosh(z k T ) sinh(z k T ) .
Hence in view of (34) we obtain

-d(k, T ) H ′′ k,T (0) + (γ -λ k ) det(M (T )) H k,T (0) = 2Re -Ξ k z 2 k z k sinh(z k T ) cosh(z k T ) + Ξ k z k cosh(z k T ) sinh(z k T ) +(γ -λ k )2Re Ξ k z k sinh(z k T ) cosh(z k T ) - Ξ k z k cosh(z k T ) sinh(z k T ) = 2Re Ξ k Ξ k 1 z k sinh(z k T ) cosh(z k T ) -Ξ k Ξ k z k cosh(z k T ) sinh(z k T ) = 0.
Now we verify that the fourth boundary condition in (32) is satisfied. We have

d(k, T ) H ′ k,T (0) = 2Re -Ξ k z k z k sinh(z k T ) sinh(z k T ) + Ξ k cosh(z k T ) cosh(z k T ) -2Re (Ξ k ) .
For the third derivative we have

d(k, T ) H ′′′ k,T (t) = 2Re Ξ k z 3 k z k sinh(z k T ) sinh(z k (t -T )) + Ξ k z 2 k cosh(z k T ) cosh(z k (t -T )) -2Re Ξ k z 2 k [sinh(z k T ) sinh(z k (t -T )) + cosh(z k T ) cosh(z k (t -T ))] . For t = 0 this yields d(k, T ) H ′′′ k,T (0) = 2 Re -Ξ k z 3 k z k sinh(z k T ) sinh(z k T ) + Ξ k z 2 k cosh(z k T ) cosh(z k T ) -2 Re Ξ k z 2 k [-sinh(z k T ) sinh(z k T ) + cosh(z k T ) cosh(z k T )] = 2 Re -Ξ k z 3 k z k sinh(z k T ) sinh(z k T ) + Ξ k z 2 k cosh(z k T ) cosh(z k T ) -2 Re Ξ k z 2 k .
Using again the definition (34) of Ξ k it follows that

-d(k, T ) H ′′′ k,T (0) + (γ -λ k ) d(k, T ) H ′ k,T (0) = 2Re Ξ k z 3 k z k sinh(z k T ) sinh(z k T ) -Ξ k z 2 k cosh(z k T ) cosh(z k T ) + 2Re Ξ k z 2 k +(γ-λ k )2Re -Ξ k z k z k sinh(z k T ) sinh(z k T ) + Ξ k cosh(z k T ) cosh(z k T ) -(γ-λ k )2Re (Ξ k ) = -2Re Ξ k 2 + 2Re -Ξ k Ξ k z k z k sinh(z k T ) sinh(z k T ) + Ξ k Ξ k cosh(z k T ) cosh(z k T ) = -2Re Ξ k 2 +|Ξ k | 2 2Re cosh(z k T ) cosh(z k T ) - z k z k sinh(z k T ) sinh(z k T ) = d(k, T )
by the definition (35) of d(k, T ). Thus we have shown that

H k,T (T ) = 0, H ′ k,T (T ) = 0, -H ′′ k,T (0) + (γ -λ k ) H k,T (0) = 0, -H ′′′ k,T (0) + (γ -λ k ) H ′ k,T (0) = 1.
The proof of Lemma 2 is complete.

In the next section, it is shown that each of the functions F k,T , G k,T , H k,T in the representation (33) of the solution b k,T that are provided in Lemma 2 satisfies an exponential turnpike inequality on the interval [0, T ] in the sense that for I ∈ {F, G, H} for all t ∈ [0, T ] the following inequality holds:

|I k,T (t)| ≤ C @ exp - √ γ 2 t + exp - √ γ 2 (T -t) .
Here C @ is a constant that is independent of T and k.

The turnpike inequality that we have obtained is used for the applications in shape optimization as described in the subsequent sections. An important point is that the inequality is independent of the properties of the sequence of eigenvalues λ k as long as [START_REF] Royden | Real Analysis[END_REF] holds. Therefore the inequality is valid over a compact set of perturbations of the state equation (see Theorem 1 and also Remark 7). [START_REF] Lance | Numerical issues and turnpike phenomenon in optimal shape design. In: Optimization and Control for Partial Differential Equations-uncertainty Quantification, Open and Closed-loop Control, and Shape Optimization[END_REF] The turnpike property by the spectral method for trees.

In this section, we continue our analysis of the structure of the optimal solutions, in particular for the adjoint states. We have in mind the shape optimization problems for the trees. Therefore, we restrict our analysis to the networks in the form of trees. The case of small cycles for the purposes of topology optimization is considered separately. In the latter case, the spectrum is of a specific structure with the branch supported exclusively on the cycle. We show that the difference between the static optimal adjoint state and the dynamic optimal adjoint state satisfies an exponential turnpike inequality as explained above. Our method is to show that all three basis functions in the representation of the difference between the static optimal adjoint state and the dynamic optimal adjoint state that we have obtained from the respective optimality systems satisfy such an exponential turnpike inequality. More precisely, we show that the basis function with a non-zero value at the time t = 0 (namely G k,T and H k,T ) decay exponentially fast with t. The basis function with a non-zero value at the time t = T (namely F k,T ) decay exponentially fast as a function of T -t, that is backward in time.

The turnpike inequalities for the basis functions

Exponential turnpike inequality for the basis function F k,T (t).

In order to verify the turnpike property we use the following representation where the hyperbolic tangent appears:

d(k, T ) F k,T (t) = -2Re Ξ k 2 cosh(z k (T -t)) + |Ξ k | 2 2Re (cosh(z k T ) cosh(z k T )[1 -tanh(z k T ) tanh(z k (T -t))]] cosh(z k (T -t))) + |Ξ k | 2 2Re cosh(z k T ) cosh(z k T ) z k z k [tanh(z k (T -t)) -tanh(z k T )] tanh(z k T ) cosh(z k (T -t)) .
For our analysis we define the auxiliary functions

S 1 k,T (t) = -2Re Ξ k 2 cosh(z k (T -t)) d(k, T ) , S 2 k,T (t) = |Ξ k | 2 2Re (cosh(z k T ) cosh(z k T )[1 -tanh(z k T ) tanh(z k (T -t))]] cosh(z k (T -t))) d(k, T ) , and
S 3 k,T (t) = |Ξ k | 2 2Re cosh(z k T ) cosh(z k T ) z k z k [tanh(z k (T -t)) -tanh(z k T )] tanh(z k T ) cosh(z k (T -t)) d(k, T ) . F k,T (t) = S 1 k,T (t) + S 2 k,T (t) + S 3 k,T (t). ( 39
)
The numbers d(k, T ) can be represented as

d(k, T ) = -2Re Ξ 2 k + |Ξ k | 2 1 - γ 2 -λ k √ 1+λ 2 k cosh(2Re (z k ) T ) + 1 + γ 2 -λ k √ 1+λ 2 k cos(2Im (z k ) T ) .
Hence we obtain the upper bound

|S 1 k,T (t)| ≤ 2| cosh(z k (T -t))| 1 - γ 2 -λ k √ 1+λ 2 k cosh(2Re (z k ) T ) -1 + γ 2 -λ k √ 1+λ 2 k cos(2Im (z k ) T ) -2Re Ξ 2 k |Ξ 2 k |
.

This implies the inequality

|S 1 k,T (t)| ≤ M 1 exp(-|Re(z k )|(T -t)) (40)
with a constant M 1 ≥ 1 that is independent of k and T . We have

1 -tanh(z k T ) tanh(z k (T -t)) = 2[exp(2 z k T ) + exp(2 z k (T -t))] (exp(2 z k T ) + 1) (exp(2 z k (T -t)) + 1) = 2 (1 + exp(-2 z k T )) (exp(2 z k (T -t)) + 1) + 2 (1 + exp(-2 z k (T -t))) (exp(2 z k T ) + 1)
.

Hence for T → ∞ this term converges to zero exponentially fast. More precisely, due to (31) for all k ∈ {0, 1, 2, ..} for T sufficiently large we have the inequality

|1 -tanh(z k T ) tanh(z k (T -t))| ≤ 16 exp(-2 |Re(z k )| (T -t)). (41)
Thus we have

|S 2 k,T (t)| ≤ 32|Re (cosh(z k T ) cosh(z k T ) cosh(z k (T -t))|) exp(-2 |Re(z k )| (T -t))| 1 - γ 2 -λ k √ 1+λ 2 k cosh(2Re (z k ) T ) -1 + γ 2 -λ k √ 1+λ 2 k cos(2Im (z k ) T ) -2Re Ξ 2 k |Ξ 2 k | .
This implies the inequality

|S 2 k,T (t)| ≤ M 2 exp(-|Re(z k )|(T -t)) (42)
with a constant M 2 ≥ 1 that is independent of k and T .

In addition, we have

tanh(z k (T -t)) -tanh(z k T ) = 2[exp(2 z k T ) -exp(2 z k (T -t))] (exp(2 z k T ) + 1) (exp(2 z k (T -t)) + 1)
so as above for T sufficiently large we obtain the inequality

| tanh(z k (T -t)) -tanh(z k T )| ≤ 16 exp(-2 |Re(z k )| (T -t)). (43)
This yields the bound

|S 3 k,T (t)| ≤ 32|Re (cosh(z k T ) cosh(z k T ) cosh(z k (T -t)) tanh(z k T )) | exp(-2 |Re(z k )| (T -t)) 1 - γ 2 -λ k √ 1+λ 2 k cosh(2Re (z k ) T ) -1 + γ 2 -λ k √ 1+λ 2 k cos(2Im (z k ) T ) -2Re Ξ 2 k |Ξ 2 k |
.

Thus we obtain the inequality

|S 3 k,T (t)| ≤ M 3 exp(-|Re(z k )|(T -t))
with a constant M 3 ≥ 1 that is independent of k and T . Together with (40) and (42) due to (39) and (31) this implies that for F k,T (t) we have

|F k,T (t)| ≤ (M 1 + M 2 + M 3 ) exp - √ γ 2 (T -t) . ( 44
)
Exponential turnpike inequality for the basis function G k,T (t). Now we derive the turnpike inequality for G k,T (t). Again we use a representation with the hyperbolic tangent. We have

d(k, T ) G k,T (t) = -2Re (Ξ k cosh(z k T ) [1 + tanh(z k T ) tanh(z k (t -T ))] cosh(z k (t -T ))) + 2Re Ξ k cosh(z k T ) 1 + z k z k tanh(z k T ) tanh(z k (t -T )) cosh(z k (t -T )) .
We have lim

T →∞ 1 + z k z k tanh(z k T ) tanh(z k (t -T )) = 1 - z k z k .
Similarly, we have lim

T →∞ [1 + tanh(z k T ) tanh(z k (t -T ))] = 0
and for the absolute value we have the bound (41).

For our analysis, we define the auxiliary functions

S1 k,T (t) = -2Re (Ξ k cosh(z k T ) [1 + tanh(z k T ) tanh(z k (t -T ))] cosh(z k (t -T ))) d(k, T ) , S2 k,T (t) = 2Re Ξ k cosh(z k T ) 1 + z k z k tanh(z k T ) tanh(z k (t -T )) cosh(z k (t -T )) d(k, T ) .
Then we have

G k,T (t) = S1 k,T (t) + S2 k,T (t).
(45) Due to (41) we obtain the inequality

| S1 k,T (t)| ≤ 32 Re Ξ k |Ξ 2 k | cosh(z k T ) cosh(z k (t -T )) exp(-2 |Re(z k )| (T -t)) 1 - γ 2 -λ k √ 1+λ 2 k cosh(2Re (z k ) T ) -1 + γ 2 -λ k √ 1+λ 2 k cos(2Im (z k ) T ) -2Re Ξ 2 k |Ξ 2 k | .
This implies the inequality

| S1 k,T (t)| ≤ 1 |Ξ k | M1 exp(-|Re(z k )|(T -t)) (46)
with a constant M1 > 0 that is independent of k and T . Moreover, we have

| S2 k,T (t)| ≤ 32 Re Ξ k |Ξ 2 k | cosh(z k T ) cosh(z k (t -T )) 1 - γ 2 -λ k √ 1+λ 2 k cosh(2Re (z k ) T ) -1 + γ 2 -λ k √ 1+λ 2 k cos(2Im (z k ) T ) -2Re Ξ 2 k |Ξ 2 k | .
This implies the inequality

| S2 k,T (t)| ≤ 1 |Ξ k | M2 [exp(-|Re(z k )| t) + exp(-|Re(z k )|(T -t))] (47)
with a constant M2 > 0 that is independent of k and T . With ( 46) and (47) due to (45) and (31) this implies that for G k,T (t) we have

|G k,T (t)| ≤ 1 |Ξ k | ( M1 + M2 ) exp(- √ γ 2 t) + exp(- √ γ 2 (T -t)) . ( 48
)
Exponential turnpike inequality for the basis function H k,T (t).

Finally we show the exponential turnpike inequality for the third basis function.

For our analysis we define the auxiliary functions

Ŝ1 k,T (t) = 2Re cosh(z k T ) Ξ k z k tanh(z k T ) + Ξ k z k tanh(z k (t -T )) cosh(z k (t -T )) d(k, T ) and Ŝ2 k,T (t) = 2Re cosh(z k T ) Ξ k z k [tanh(z k T ) + tanh(z k (t -T ))] cosh(z k (t -T )) d(k, T ) .
Then we have

H k,T (t) = Ŝ1 k,T (t) -Ŝ2 k,T (t). ( 49
)
For the factor with the hyperbolic tangent that appears in Ŝ1 k,T (t) for T sufficiently large we have the bound

Ξ k z k tanh(z k T ) + Ξ k z k tanh(z k (t -T )) ≤ Ξ k z k [| tanh(z k T )| + | tanh(z k (t -T ))|] ≤ 4 γ + 1 λ k .
This yields the inequality

| Ŝ1 k,T (t)| ≤ 8 |Ξ k | 2 Re cosh(z k T ) γ + 1 λ k cosh(z k (t -T )) 1 - γ 2 -λ k √ 1+λ 2 k cosh(2Re (z k ) T ) -1 + γ 2 -λ k √ 1+λ 2 k cos(2Im (z k ) T ) -2Re Ξ 2 k |Ξ 2 k | .
This implies the inequality

| Ŝ1 k,T (t)| ≤ 1 |Ξ k | 2 M1 [exp(-|Re(z k )| t) + exp(-|Re(z k )|(T -t))] (50)
with a constant M1 > 0 that is independent of k and T . For the factor with the hyperbolic tangent that appears in Ŝ2 k,T (t) for T sufficiently large we have the bound

|tanh(z k T ) + tanh(z k (t -T ))| ≤ 4.

This yields the inequality

| Ŝ2 k,T (t)| ≤ 8 |Ξ k | 2 Re cosh(z k T ) γ + 1 λ k cosh(z k (t -T )) 1 - γ 2 -λ k √ 1+λ 2 k cosh(2Re (z k ) T ) -1 + γ 2 -λ k √ 1+λ 2 k cos(2Im (z k ) T ) -2Re Ξ 2 k |Ξ 2 k | .
This implies the inequality

| Ŝ2 k,T (t)| ≤ 1 |Ξ k | 2 M2 [exp(-|Re(z k )| t) + exp(-|Re(z k )|(T -t))] (51)
with a constant M2 > 0 that is independent of k and T . With (50) and (51) due to (49) and (31) this implies that for H k,T (t) we have

|H k,T (t)| ≤ 1 |Ξ k | 2 ( M1 + M2 ) exp(- √ γ 2 t) + exp(- √ γ 2 (T -t)) . (52)

Turnpike Theorem between the dynamic and static Optimality Systems

We have unique solutions for the optimality systems governed by static and dynamic state equations. Now, we compare the elements of optimality systems and obtain the turnpike inequalities for the Optimality Systems. Indeed, the turnpike inequalities that we have derived lead to the useful for applications Turnpike Theorem which means that :

The difference ν T of the optimal dynamic control and the optimal static control and the corresponding differences ω T for the state and µ T the adjoint state admits an exponential turnpike property.

To this end, we need some additional regularity assumptions. Theorem 1 Assume that (28) holds and that the initial state satisfies the regularity condition

∞ k=0 λ k |a k (0)| 2 + |a ′ k (0)| 2 < ∞, (53)
that is the initial state belongs to the energy space of the elliptic problem defined by the bilinear form a(•, •). If Ω = Γ then there exists a constant D = D(y 0 , y 1 , p σ ) that is independent of T and t such that for all t ∈ [0, T ]

∥ω T (t)∥ 2 L 2 (Ω) + ∥ν T (t)∥ 2 L 2 (Ω) + ∥µ T (t)∥ 2 L 2 (Ω) ≤ D e -√ γ t + e -√ γ(T-t) . (54)
Moreover, the constant D depends on Ω only as a function of the energy norm for the initial state that is determined by Ω as in (53). Remark 7 In the case of boundary control problems on networks, the Turnpike Property for optimality systems can be shown. The difference between the boundary control and the distributed control is the appearance of a linear operator in the optimality conditions. In the optimality conditions of the boundary control problem, a linear operator maps the adjoint state into the optimal control. For the optimal cost the exponential turnpike inequality (54) implies the so-called integral turnpike property (see e.g. [START_REF] Gugat | On the turnpike phenomenon for optimal boundary control problems with hyperbolic systems[END_REF])

sup T >0 T 0 ∥ω T (t)∥ 2 L 2 (Ω) + ∥ν T (t)∥ 2 L 2 (Ω) + ∥µ T (t)∥ 2 L 2 (Ω) dt < ∞.
This implies in turn

lim T →∞ 1 T J T (û T ) = I(v), (55)
see for example [START_REF] Gugat | Turnpike properties for partially uncontrollable systems[END_REF].

In shape sensitivity analysis [START_REF] Sokolowski | Zolésio: Introduction to Shape Optimization[END_REF] we avoid, if possible, the dependence of the shape gradient of the cost with respect to given data including e.g., the given initial conditions. Thus, we define the initial conditions, say y 0 (x) and y 1 (x), x ∈ Ω τ , for variable domains τ → Ω τ used for the derivation of shape gradient of the cost. If Ω τ is e.g., the cross with variable length of the edges, we select the elements y 0 , y 1 as the restrictions to Ω τ of functions Y 0 (x), Y 1 (x) defined on the cross with the edges of maximal lengths.

The shape derivatives of the restrictions are zero. Therefore, there is no contribution from the initial conditions to the shape gradient of the cost. Remark 8 We use the material derivative method [START_REF] Sokolowski | Zolésio: Introduction to Shape Optimization[END_REF] for the purposes of shape sensitivity analysis for networks. The general rule for the data of initial-boundary value problems is to select the given functions by the restriction to actual domain of some functions defined everywhere. The shape gradients of such initial data are simply zero and the material derivatives are given by the gradients, thus some regularity is required. In this way the shape gradient of the cost is independent of the initial data.

In our case, this selection can be used for the initial conditions of the displacement and the velocity. The material derivatives of such initial conditions Y = Y (x) take the form Y ′ (x)V(0, x), where V(τ, x) is the velocity field of the material derivative method.

In another words, the initial conditions are selected in such a way that there is no contribution of the initial conditions to the shape gradient of the cost function.

The following corollary states a shape-turnpike result. It is a relation between the optimal values of a dynamic optimal shape problem for large time horizons and the optimal values for the static optimal shape problem. For the proof we suggest to proceed by contradiction. Assumption 1 Let us consider a tree G = {E, V } with the set of edges E i = [0, L i ], i = 1, . . . , N , and denote by ℓ = col (L 1 , . . . , L N ). The set of admissible trees Ω(ℓ) ∈ A is defined by the conditions

M min i ≤ L i ≤ M max i , where 0 < M 0 ≤ M min i < M max i < M 1 < ∞.
The set of admissible trees is convex and compact, therefore, for the minimizing sequence ℓ n of optimization problem there is a subsequence, still denoted by the same symbol such that we have ℓ n → ℓ ∞ in R N , and in addition Ω(ℓ ∞ ) ∈ A. Corollary 1 Let Assumption 1 hold. Let a sequence of shape parameters (ℓ n ) n and a bounded sequence of controls (û n ) n with ûn (t) ∈ L 2 (Ω(ℓ n )) for all n ∈ {1, 2, 3, ...} be given. Let ŷn denote the generated state and pn the corresponding adjoint state. Assume that for all n ∈ {1, 2, 3, ...} we have

∥ω T n (t)∥ 2 L 2 (Ω(ℓn)) +∥ν T n (t)∥ 2 L 2 (Ω(ℓn)) +∥µ T n (t)∥ 2 L 2 (Ω(ℓn)) ≤ D e -√ γ t + e -√ γ(T-t) (56)
where ω T n = ŷT n -ẑσ n , µ T n = pT n -pσ n , ν T n = ûT n -vσ n and (v n , ẑn , pn ) is optimal for Ω(ℓ n ) Assume that lim n→∞ ℓ opt (T n ) = l and that ûn converges weakly to u opt . Assume that the optimal shape problem with (OCE) has a solution (ℓ opt , u opt ). Then u opt is a solution of (OCE).

Let vopt denote the solution of (OCS) for ℓ opt . Then we have

∥ω T opt (t)∥ 2 L 2 (Ω(ℓopt)) +∥ν T opt (t)∥ 2 L 2 (Ω(ℓopt)) +∥µ T opt (t)∥ 2 L 2 (Ω(ℓopt)) ≤ D e -√ γ t + e -√ γ(T-t) .
(57) Assume that for a subsequence we have lim Tn→∞ ℓ opt (T n ) = l. Then we have lim

n→∞ 1 T n J Tn (û Tn opt , ℓ opt (T n )) = I(v, l). ( 58
)
and v is the optimal control for the network defined in Ω( l), the limit shape parameter l is optimal for the static problem ℓ → (OCS)(ℓ), the optimal shape reads Ω( l).

In general, an optimal shape is not unique but it does exist for the compact set of admissible shapes. Remark 9 Under Assumption 1 the set of admissible shapes is convex and compact in R N . Denote by l an optimal shape for the static problem, note that the optimal shape is not unique, and let the admissible sequence of shapes ℓ n be convergent to an optimal shape for n → ∞,

ℓ n → l then I(v, l) = lim n→∞ 1 T n J Tn (û Tn opt , ℓ n ) ≥ lim n→∞ 1 T n J Tn (û Tn opt , ℓ opt (T n )) = I(v, l) (59) therefore I(v, l) = I(v, l).
Proof of the corollary. The constant control with the value v can be considered as an element of L 2 (Ω(ℓ opt (T n ))). We have the inequality

∥û Tn opt (ℓ opt (T n )) -u d ∥ 2 L 2 (Ω(ℓopt(Tn))) -∥v -u d ∥ 2 L 2 (Ω(ℓopt(Tn))) ≤ ∥û Tn opt (ℓ opt (T n ))-v∥ L 2 (Ω(ℓopt(Tn))) ∥û Tn opt (ℓ opt (T n )) -u d ∥ L 2 (Ω(ℓopt(Tn))) + ∥v -u d ∥ L 2 (Ω(ℓopt(Tn)))
. Integration from 0 to T n and division by T n yields Tn

0 ∥û Tn opt (ℓ opt (T n )) -u d ∥ 2 L 2 (Ω(ℓopt(Tn))) dt T n -∥v -u d ∥ 2 L 2 (Ω(ℓopt(Tn))) ≤ Tn 0 ∥û Tn opt (ℓopt(Tn))-v∥ L 2 (Ω(ℓ opt (Tn ))) ∥û Tn opt (ℓopt(Tn))-u d ∥ L 2 (Ω(ℓ opt (Tn ))) +∥v-u d ∥ L 2 (Ω(ℓ opt (Tn ))) dt Tn ≤ Tn 0 ∥û Tn opt (ℓ opt (T n )) -v∥ 2 L 2 (Ω(ℓopt(Tn))) dt T n + 2∥v -u d ∥ L 2 (Ω(ℓopt(Tn)))
Tn

0 ∥û Tn opt (ℓ opt (T n )) -v∥ L 2 (Ω(ℓopt(Tn))) dt T n .
For the state we obtain a similar inequality, namely

Tn 0 ∥ŷ Tn opt (ℓ opt (T n )) -y d ∥ 2 L 2 (Ω(ℓopt(Tn))) dt T n -∥ẑ -y d ∥ 2 L 2 (Ω(ℓopt(Tn))) ≤ Tn 0 ∥ŷ Tn opt (ℓ opt (T n )) -ẑ∥ 2 L 2 (Ω(ℓopt(Tn))) dt T n + 2∥ẑ -y d ∥ L 2 (Ω(ℓopt(Tn)))
Tn

0 ∥ŷ Tn opt (ℓ opt (T n )) -ẑ∥ L 2 (Ω(ℓopt(Tn))) dt T n .
For the optimal cost the exponential turnpike inequality (54) implies the integral turnpike property (see e.g. [START_REF] Gugat | On the turnpike phenomenon for optimal boundary control problems with hyperbolic systems[END_REF])

sup T >0 sup ℓ T 0 ∥ω T (t)∥ 2 L 2 (Ω) + ∥ν T (t)∥ 2 L 2 (Ω) + ∥µ T (t)∥ 2 L 2 (Ω) dt < ∞.
Moreover, we have sup

T >0 sup ℓ T 0 ∥ω T (t)∥ L 2 (Ω) + ∥ν T (t)∥ L 2 (Ω) + ∥µ T (t)∥ L 2 (Ω) dt < ∞.
Thus adding up the inequalities for the control and the state and taking the limit for T n → ∞ yields (58).

Remark 10 For tree-shaped graphs often the system is exactly controllable in some finite time t min . Recently it was shown that if there is also control action at the interior nodes of the graph, exact controllability is also possible for a graph with cycles, see [START_REF] Avdonin | Controllability for the wave equation on graph with cycle and delta-prime vertex conditions[END_REF].

In this case we can choose for all n ∈ {1, 2, 3, ...} a control function u

(n) init (t) ∈ L 2 0, t min ; L 2 (Ω(ℓ opt (T n ))
) that steers the system to the constant state y(t min , •) = ŷ, y t (t min , •) = 0 and satisfies

max n ∥u (n) init ∥ L 2 (0, tmin;L 2 (Ω(ℓopt(Tn)))) < ∞. ( 60
)
Such a control can be determined using the classical method of moments as described for example in [START_REF] Rolewicz | On controllability of systems of strings[END_REF].

We define the control ṽ(n

) (t) = u (n) init (t), t ∈ (0, t min ), v, t ≥ t min . Then ṽ(n) is feasible for (OCE)(T n , Ω(ℓ opt (T n ))
) and thus we have

J Tn (û Tn opt , ℓ opt (T n )) ≤ J Tn (ṽ (n) , ℓ opt (T n )). ( 61
)
Moreover, (58) implies that for all ε > 0 if T n > 0 is sufficiently large we have

lim Tn→∞ 1 T n J Tn (ṽ (n) , ℓ opt (T n )) ≤ 1 T n J Tn (û Tn opt , ℓ opt (T n )) + ε. ( 62
)
This can be seen as follows. Since lim n→∞ T n = ∞ due to (63) the contribution of the integral on the time interval (0, t min ) vanishes in the limit, that is we have

lim n→∞ tmin 0 ∥y(ṽ (n) ) -y d ∥ 2 L 2 (Ω) + γ ∥∂ t (y(ṽ (n) ) -y d )∥ 2 L 2 (Ω) + ∥ṽ (n) -u d ∥ 2 L 2 (Γ) dt T n = 0
and since lim Tn→∞ ℓ opt (T n ) = l we have

lim n→∞ J Tn (ṽ (n) , ℓ opt (T n )) T n = lim n→∞ T n -t min T n I(v, ℓ opt (T n )) = I(v, l).
Assuming the exact controllability of the system in the finite time t min allows to choose for all n ∈ {1, 2, 3, ...}, a control function u

(n) init (t) ∈ L 2 0, t min ; L 2 (Ω(ℓ opt (T n ))
) that steers the system to the constant state y(t min , •) = ŷ, y t (t min , •) = 0 and satisfies

max n ∥u (n) init ∥ L 2 (0, tmin;L 2 (Ω(ℓopt(Tn)))) < ∞. (63)
Such a control can be determined using the classical method of moments as described for example in [START_REF] Rolewicz | On controllability of systems of strings[END_REF].

We define the control ṽ(n

) (t) = u (n) init (t), t ∈ (0, t min ], v, t ≥ t min . For all n ∈ {1, 2, 3, ...} this implies the inequality 1 T n J Tn (û Tn opt , ℓ opt (T n )) ≤ 1 T n J Tn (ṽ (n) , ℓ opt (T n )).
This yields (61). Assume that the optimal shape problem with (OCE) has a solution (ℓ opt , u opt ). Then u opt is a solution of (OCE). Let vopt denote the solution of (OCS) for ℓ opt . Then Theorem 1 implies that (57) holds.

Proof of Theorem 1.

Define M F = M 1 + M 2 + M 3 , M G = M1 + M2 , M H = M1 + M2 .
Due to (33), ( 44), ( 48) and ( 52) for all k ∈ {0, 1, 2, 3, ...} we have

|b k,T (t)| ≤ M F e - √ γ 2 (T -t) |b k,T (T )| (64)
+[ 1 + γ λ k M G |a k (0)| + M H |a ′ k (0)|] exp(- √ γ 2 t) + exp(- √ γ 2 (T -t)) .
For the square, this yields the bound

|b k,T (t)| 2 ≤ 3M 2 F e -√ γ(T -t) |b k,T (T )| 2 (65)
+[6(1 + γ λ k ) M 2 G |a k (0)| 2 + 6M 2 H |a ′ k (0)| 2 ] [exp(- √ γ t) + exp(- √ γ(T -t))] .
Due to Parseval's equation this implies for all t ∈ [0, T ]

∥µ T (t)∥ 2 L 2 (Ω) = ∞ k=0 |b k (t)| 2 ≤ ∞ k=0 3 M 2 F e -√ γ(T -t) |b k (T )| 2 +[6(1 + γ λ k ) M 2 G |a k (0)| 2 + 6M 2 H |a ′ k (0)| 2 ] [exp(- √ γ t) + exp(- √ γ(T -t))] .
Since the initial state (y 0 , y 1 ) satisfies

∞ k=0 λ k |a k (0)| 2 + |a ′ k (0)| 2 < ∞,
this yields an exponential turnpike property for the adjoint state. To be precise, we have for all t ∈ [0, T ]

∥µ T (t)∥ 2 L 2 (Ω) ≤ C(y 0 , y 1 , p σ ) [exp(- √ γ t) + exp(- √ γ(T -t))] (66)
with a real number C(y 0 , y 1 , p σ ) that is independent of T and k.

For the optimal controls, for Γ = Ω we have ν T = µ T , hence we have a similar inequality as (66) for ∥ν T (t)∥ 2 L 2 (Ω) . Note that in proof of (44), ( 48) and (52) in the estimates we did not take advantage of the real part that appears in the representations. In fact, we have always used upper bound for the modulus of the complex number whose real part appears in the expressions. In the representation of the second order derivatives, the only change is that the factor z 2 k appears in the complex variable representations compared to the original expression. Therefore our proof also yields the inequality

|b ′′ k,T (t)| 2 ≤ 3|z k | 2 M 2 F e -√ γ(T -t) |b k,T (T )| 2 (67) +|z k | 2 [6(1 + γ λ k ) M 2 G |a k (0)| 2 + 6M 2 H |a ′ k (0)| 2 ] [exp(- √ γ t) + exp(- √ γ(T -t))] .
We have

a k = -b ′′ k + (γ -λ k ) b k 1 + γ λ k . (68)
Define

M α := sup k |z k | 2 1 + γ λ k = sup k 1 + λ 2 k 1 + γ λ k < ∞, M β := sup k λ k -γ 1 + γ λ k < ∞ and M κ := M F |b k,T (T )| + 1 + γ λ k M G |a k (0)| + M H |a ′ k (0)|.
Due to (68) we also have the following exponential inequality for the coefficients in the expansion of the optimal state:

|a k (t)| ≤ (M α + M β ) M κ e - √ γ 2 t + e - √ γ 2 (T -t) .
This yields the turnpike inequality for the state ω T :

∥ω T (t)∥ 2 L 2 (Ω) ≤ C(y 0 , y 1 , p σ ) [exp(- √ γ t) + exp(- √ γ(T -t))] (69)
with a real number C(y 0 , y 1 , p σ ) that is independent of T . Note that in proof of (44), ( 48) and (52) in the estimates we did not take advantage of the structure of the spectrum. In fact, the constants in the turnpike inequalities (44) and M G = M1 + M2 in (48) and M H = M1 + M2 in (52) are independent of k. Remark 11 From mathematical point of view, it is of interest to analyze the nucleation of a cycle at the internal node of network. This means that the internal node is replaced by a small cycle, see Figures 11,12. The question which internal node is selected can be solved in the steady state case, the topological derivative of the cost is introduced to this end. The domain decomposition method is applied in order to derive the form of the topological derivative. The associated Steklov-Poincaré operator is represented by a matrix ε → Λ(ε) which is semidefinite positive and differentiable at ε = 0+, the derivative is denoted Λ ′ (0) for the sake of simplicity. In the case of wave equation such a result is not known, however the numerical experiments show that dependence is regular. Therefore, we assume that for a given initial conditions the solution of wave equation enjoys the properties of the steady state boundary value problem for the nucleation of the small cycle. Let us note that for ε > 0 the wave equation is well defined and the regularity of initial conditions required for the turnpike property are already given. Here we assume that the limit of the cost for ε → 0+ is well defined and the value of the cost is continuous ε = 0. This assumption does not imply the differentiability of the cost at ε → 0+. Assumption 2 Let us consider the wave equation on network and let J (Ω) be the cost for given initial conditions y 0 , y 1 and given time horizon T > 0. At the internal node P 0 of the network a cycle of size ε > 0 is introduced, which leads to the cost J (Ω ε ). The topological variation of the network is admissible provided we have

lim ε→0+ J (Ω ε ) = J (Ω).

A simple example

Let real numbers L > 0, T > 0, c > 0 and γ > 0 be given. Let y 0 ∈ H 1 (0, L) with y 0 (0) = 0 and

y 1 ∈ L 2 (0, L), z ∈ H 1 (0, L) with ζ = z ′ (L) be given. Consider the problem min T 0 L 0 |y(t, x) -z(x)| 2 + γ |y t (t, x)| 2 dx + |u(t) -ζ| 2 dt subject to            y(0, x) = y 0 (x), x ∈ (0, L), y t (0, x) = y 1 (x), x ∈ (0, L), y tt (t, x) = c 2 y xx (t, x), (t, x) ∈ (0, T ) × (0, L), y(t, 0) = 0, t ∈ (0, T ), y x (t, L) = u(t), t ∈ (0, T ). ( 70
)
The solution to the initial boundary value problem for a control u ∈ L 2 (0, T ) is stated in [START_REF] Gugat | Optimal Boundary Control and Boundary Stabilization of Hyperbolic Systems[END_REF], Theorem 2.3, p. 17 in the form

y(t, x) = ∞ n=0 α n (t) φ n (x) (71)
with the eigenfunctions

φ n (x) = √ 2 √ L sin π 2 + nπ x L , n ∈ {0, 1, 2, ...}.
For the eigenvalues we have

λ n = 1 L 2 π 2 + nπ 2 and the minimal eigenvalue is λ 0 = π 2 4 1 L 2 . For n ∈ {0, 1, 2, ...}, let α 0 n = L 0 y 0 (x) φ n (x) dx, α 1 n = L 0 y 1 (x) φ n (x) dx. (72)
We have

α n (t) = α 0 n cos ( π 2 + nπ) t t 0 + α 1 n t 0 π 2 + nπ sin ( π 2 + nπ) t t 0 (73) +(-1) n c 2 √ 2 √ L t 0 π 2 + nπ t 0 u(s) sin ( π 2 + nπ) t -s t 0 ds,
where t 0 = L c . Since Parseval's identity states that almost everywhere on [0, T ], we have

L 0 y(t, x) 2 dx = ∞ n=0 |α n (t)| 2 and L 0 y t (t, x) 2 dx = ∞ n=0 |α ′ n (t)| 2 .
We discuss the regularization for the third term of α n (t) in appendix A. We can represent the objective functional in the form

T 0 L 0 |y(t, x) -z(x)| 2 + γ |y t (t, x)| 2 dx + |u(t) -ζ| 2 dt = T 0 ∞ n=0 |α n (t)| 2 + γ|α ′ n (t)| 2 + |u(t) -ζ| 2 dt + T 0 L 0 z 2 (x) + ∞ n=0 2α n (t)φ n (x)z(x) dxdt.
If the control space L 2 (0, T ) is replaced by a finite dimensional space of piecewise constant control functions, this yields a finite dimensional quadratic optimization problem. Since in this case, the necessary optimality conditions are a finite dimensional system of linear equations, this can be used to obtain numerical approximations of the optimal control.

Numerical solutions of three examples

In this section, we discuss the numerical solution. We consider the following problem data: c := 1; γ := 0.1; T := 1, 10, 100; L := 1;

We choose y 1 (x) := 0 (x ∈ (0, 1)).

1. y 0 (x) := x, z = 0, ζ = 0 (example 1); 2. y 0 (x) := π -1 sin(πx), z = 0, ζ = 0 (example 2); 3. y 0 (x) := π -1 sin(πx), z = x, ζ = 1 (example 3);

The coefficients α 0 n and α 1 n in (72) are obtained as the Fourier coefficients of the chosen functions y 0 (x), y 1 (x).

In order to solve the optimal boundary control problem numerically, a finitedimensional approximation is used on two sides of (70) simultaneously:

First, the series expansion in the objective functional has to be cut after N terms which leads us to the consideration of the problem (OPT)(γ, N ). Second, we compute approximations for the optimal controls u ∈ L 2 (0, T ) in the space of piecewise constant functions. Let a grid 0 = t 0 < t 1 < t 2 < ... < t M = T be given. For i ∈ {1, . . . , M } let v j (t) := 1 if t ∈ [t j-1 , t j ), 0 elsewhere, and, define the finite dimensional space X M (T ) by X M (T ) := span{v j (•) : j = 1, . . . , M }.

For any u ∈ X M (T ) we use the representation

u(t) = M j=1 u(t j-1 )v j (t), t ∈ [0, T ),
where v j (t) stands for the characteristic function χ [tj-1,tj ) (t) of the interval [t j-1 , t j ) and the approximation of control is defined by the vector U = col(U 1 , . . . , U M ) ∈ R M . Hence, we are finally led to solve the problem

(D opt )(γ, N, M ) min u∈X M (T ) M j=1 (t j -t j-1 ) (u(t j-1 ) 2 + 2u(t j-1 )ζ + ζ 2 ) + M j=1 N n=0 tj tj-1 |α n (t)| 2 + γ|α ′ n (t)| 2 dt + M j=1 N n=0 tj tj-1 L 0 2α n (t)φ n (x)z(x) + z 2 (x) dxdt (74)
with α j (t) as defined in (73). Problem (D opt )(γ, N, M ) can be equivalently formulated as a quadratic programming problem in R M min

U ∈R M U ⊤ Q U + q ⊤ U + W
where the matrix Q(N ) M ×M , depending on the fixed number N , and the vector q ∈ R M are to be assembled for fixed N from the cost as stated in (74). The assemblage is described below. We use the notation U j = u(t j-1 ), (j ∈ {1, . . . , M }) and take into account a constant term W that is independent of U . We have

U ⊤ QU + q ⊤ U + W = M j=1 (t j -t j-1 ) (u(t j-1 ) 2 + 2u(t j-1 )ζ + ζ 2 ) + M j=1 N n=0 tj tj-1 |α n (t)| 2 + γ|α ′ n (t)| 2 dt + M j=1 N n=0 tj tj-1 L 0 2α n (t)φ n (x)z(x) + z 2 (x) dxdt.
For the convergence of the approximation, it is important to increase both N and M simultaneously. Otherwise, if only M → ∞ convergence to the optimal control in general does not occur (due to a possible spillover effect). For t ∈ (t j-1 , t j ) we have

α n (t) = α 0 n cos ( π 2 + nπ)t + (-1) n √ 2 π 2 + nπ t 0 u(s) sin ( π 2 + nπ)(t -s) ds = α 0 n cos ( π 2 + nπ)t + (-1) n √ 2 π 2 + nπ j k=1 U k t 0 v k (s) sin ( π 2 + nπ)(t -s) ds = α 0 n cos ( π 2 + nπ)t + (-1) n √ 2 π 2 + nπ V ⊤ j U,
where

V j = col(V j,1 , V j,2 , • • • , V j,j , 0, • • • , 0), V j,k =      t k t k-1 sin ( π 2 + nπ)(t -s) ds if k < j, t t k-1 sin ( π 2 + nπ)(t -s) ds if k = j, 0 if k > j.
and

|α n (t)| 2 = (α 0 n ) 2 cos 2 ( π 2 + nπ)t + 2 ( π 2 + nπ) 2 U ⊤ V j V ⊤ j U + (-1) n 2 √ 2α 0 n π 2 + nπ cos ( π 2 + nπ)t V ⊤ j U. (75)
This implies

M j=1 N n=0 tj tj-1 |α n (t)| 2 = M j=1 N n=0 (α 0 n ) 2 tj tj-1 cos 2 ( π 2 + nπ)t dt + M j=1 N n=0 2 ( π 2 + nπ) 2 U ⊤ tj tj-1 V j V ⊤ j dt U + M j=1 N n=0 (-1) n 2 √ 2α 0 n π 2 + nπ tj tj-1 cos ( π 2 + nπ)t V ⊤ j dt U = W 1 + U ⊤ Q 1 U + q ⊤ 1 U. (76)
Moreover, for the derivatives α ′ n (t), we have for t ∈ (t j-1 , t j )

α ′ n (t) = -( π 2 + nπ)α 0 n sin ( π 2 + nπ)t + (-1) n √ 2 j k=1 U k t 0 v k (s) cos ( π 2 + nπ)(t -s) ds = -( π 2 + nπ)α 0 n sin ( π 2 + nπ)t + (-1) n √ 2 Ṽ ⊤ j U, ( 77
) where Ṽj = col( Ṽj,1 , Ṽj,2 , • • • , Ṽj,j , 0, • • • , 0), Ṽj,k =      t k t k-1 cos ( π 2 + nπ)(t -s) ds (k < j), t t k-1 cos ( π 2 + nπ)(t -s) ds (k = j), 0 (k > j).
This yields

|α ′ n (t)| 2 = (α 0 n 2 ( π 2 + nπ) 2 sin 2 ( π 2 + nπ)t + 2U ⊤ Ṽj Ṽ ⊤ j U + (-1) n+1 2 √ 2α 0 n ( π 2 + nπ) sin ( π 2 + nπ)t Ṽ ⊤ j U.
Hence we obtain

γ M j=1 N n=0 tj tj-1 |α ′ n (t)| 2 = γ M j=1 N n=0 (α 0 n ) 2 ( π 2 + nπ) 2 tj tj-1 sin 2 ( π 2 + nπ)t dt + γ M j=1 N n=0 2U ⊤ tj tj-1 Ṽj Ṽ ⊤ j dt U + γ M j=1 N n=0 (-1) n+1 2 √ 2α 0 n ( π 2 + nπ) tj tj-1 sin ( π 2 + nπ)t Ṽ ⊤ j dt U = γ(W 2 + U ⊤ Q 2 U + q ⊤ 2 U ).
And we have

M j=1 N n=0 tj tj-1 α n (t)φ n (x)z(x) = M j=1 N n=0 tj tj-1 α 0 n cos ( π 2 + nπ)t dt L 0 z(x)φ(x)dx + M j=1 N n=0 (-1) n √ 2 π 2 + nπ tj tj-1 V ⊤ j dt U = W 3 + q ⊤ 3 U.
(78) For the objective function, this implies

M j=1 (t j -t j-1 ) (u(t j-1 ) 2 + 2u(t j-1 )ζ + ζ 2 ) := U ⊤ Q 3 U + 2q ⊤ τ U + W τ , M j=1 N n=0 tj tj-1 |α n (t)| 2 + γ|α ′ n (t)| 2 dt := U + q ⊤ 1 U + W 1 + γ(U ⊤ Q 2 U + q ⊤ 2 U + W 2 ), M j=1 N n=0 tj tj-1 L 0 2α n (t)φ n (x)z(x) + z 2 (x) dxdt = 2q ⊤ 3 U + W 3 .
Thus

U ⊤ QU +q ⊤ U +W = U ⊤ (Q 1 +γQ 2 +Q 3 )U +(q 1 +γq 2 +q 3 +q τ ) ⊤ U +W 1 +γW 2 +W 3 +W τ .
Here we use the notation

Q 1 = M j=1 N n=0 2 ( π 2 + nπ) 2 tj tj-1 V j V ⊤ j dt, Q 2 = M j=1 N n=0 2 tj tj-1 Ṽj Ṽ ⊤ j dt, Q 3 = diag(t 1 -t 0 , • • • , t k -t k-1 , • • • , t M -t M -1 ), q 1 = M j=1 N n=0 (-1) n 2 √ 2α 0 n π 2 + nπ tj tj-1 cos ( π 2 + nπ)t V j dt, q 2 = M j=1 N n=0 (-1) n+1 2 √ 2α 0 n ( π 2 + nπ) tj tj-1 sin ( π 2 + nπ)t Ṽj dt, q 3 = M j=1 N n=0 (-1) n √ 2 π 2 + nπ tj tj-1 V ⊤ j dt, q τ = (t 1 -t 0 , t 2 -t 1 , • • • , t M -t M -1 ) ⊤ , W 1 = M j=1 N n=0 (α 0 n ) 2 tj tj-1 cos 2 ( π 2 + nπ)t dt, W 2 = M j=1 N n=0 (α 0 n ) 2 tj tj-1 sin 2 ( π 2 + nπ)t dt, W 3 = M j=1 N n=0 tj tj-1 α 0 n cos ( π 2 + nπ)t dt L 0 z(x)φ(x)dx.
We employ Matlab for the computational analysis of all examples. Figure 1 and Figure 2 illustrate the optimal control and state, respectively, for varying values of T in example 1. Notably, as T increases significantly, the control variable u converges to ζ. Moreover, with the increment in T , there is a discernible trend towards stabilization in the norms of both u and y, as evidenced in Figures 3. Fig. 4-6 present the results obtained from example 2, whereas Fig. 7 to 9 depicts the outcomes of example 3. Additionally, Fig. 10 presents the quotient of the optimal value of (D opt )(γ, N, M ), which adheres to the turnpike property, thereby providing valuable insights into the behavior of the system concerning varying T values. The topology of network for the purposes of optimal control problems is selected in the framework of the topological derivative method for static problems, we refer to [START_REF] Gugat | Proceedings of the 27th International Conference on Methods and Models in Automation and Robotics (MMAR)[END_REF] for an elementary example. It turns out, that that the topological derivatives for a class of cost functions can be determined by using the domain decomposition technique for the state equation [START_REF] Soko Lowski | Modelling of topological derivatives for contact problems[END_REF]. We describe in details the topological derivative method and present numerical results for examples. Let us consider the network static problem. We define the topological derivatives for optimal cost of network control problems with respect to nucleation of a small cycle. We present also numerical examples. The simplest example of a network is the three-star graph with one central vertex P 0 and three boundary vertices P 1 , P 2 , P 3 , thus V = {P 0 , P 1 , P 2 , P 3 } (See Fig. 11). There are three edges E = {E 1 , E 2 , E 3 }.

For the steady state problem, singular domain perturbations of the shape are considered. The topological derivatives of the shape functional are defined. The shape and topology optimization is performed. The network is singularly perturbed by a small cycle of the size ε → 0 (See Fig. 12). In such a case the domain decomposition technique is used and the Steklov-Poincaré operator is introduced. The topological We introduce multiple perturbations of network represented in Figure 13. The Steklov-Poincaré operator Λ ε replaces the subgraph G ε in the state equation of the network. In this way the topological derivative of the cost for optimal control problem on perturbed network is obtained for the nucleation of multiple cycle in the three-star graph.

x = 0 x = l 3 x = l 2 x = l 1

Shape and topology optimization on networks

We recall briefly the shape and topological derivatives of a given cost for the network. We restrict ourselves to static problems however the shape and topology optimization can also be performed for dynamic optimal control problems on networks. For the optimal control problems with the Turnpike Property the analysis of static problem is useful for the solution of dynamic problem. In particular, the topology of the network is designed using the static problem.

E ε,1 P 1 E ε,6 Q 5 E ε,5 P 5 E ε,7 P 6 E ε,4 E ε,8 P 4 Q 4 E ε,9 E ε,2 P 2 E ε,3 P 3 Q 1 Q 2 Q 3 G ε
The shape Ω of the network for fixed topology is governed by the finite dimensional vector ℓ, which contains the lengths of edges, ℓ = col (L 1 , . . . , L N ) where N = #E = {E i |i ∈ I}. Therefore, Ω := Ω(ℓ), and the cost ℓ → I(ℓ) := J (Ω(ℓ)) is defined by the optimal cost of control problem J (Ω) := J(û(Ω)) for evolution problem. For the steady state problem the optimal cost of control problem is denoted by ℓ → J(v(Ω(ℓ))), where v(Ω(ℓ)) is steady state optimal control in the domain defined by the shape Ω(ℓ). We consider the Neumann control in the numerical examples presented for the wave equation.

We use the standard technique for shape [START_REF] Sokolowski | Zolésio: Introduction to Shape Optimization[END_REF] and topology optimization [START_REF] Novotny | An Introduction to the Topological Derivative Method[END_REF]. Namely, the material derivatives are employed in the shape sensitivity analysis in the framework of speed method [START_REF] Sokolowski | Zolésio: Introduction to Shape Optimization[END_REF]. The parameter τ → 0 is used for the boundary variations of τ → J (Ω τ ), at τ = 0, the material derivatives of the cost and the state with respect to the parameter τ are denoted with the dot

J (Ω) = lim τ →0 1 τ (J (Ω τ ) -J (Ω))
which leads to the expansion

J (Ω τ ) = J (Ω) + τ J (Ω) + o(τ ).
In addition, the topology optimization is considered in the case of nucleation of the small cycle of the size ε → 0. The topological derivative of the cost is defined as follows

T (Ω) = lim ε→0 1 ε (J (Ω ε ) -J (Ω))
which leads to the expansion

J (Ω ε ) = J (Ω) + εT (Ω) + o(ε).
The topological derivative T (Ω) for nucleation of small cycles of the size ε → 0 is evaluated for two examples of static problems, and applied to the topological optimization for dynamic problem in the first example.

Examples of topological derivatives for networks

Two examples are presented of singular network perturbation by nucleation of a small cycle. In the first example, the topological derivative is evaluated for the optimal control problem with the static state equation, and then for the dynamic state equation, the optimal size of the cycle is determined. In the second example, the topological derivative is evaluated for multiple singular perturbations of the three-star network.

Example for one cycle: In this example the optimal control u is computed for the geometry depicted in (79)

E ε,1 P 1 P 5 E ε,5
E ε,4

P 6 P 4 E ε,6 E ε,2
E ε,3 By the Lagrange method, the optimality system is given by

P 2 P 3 Q 1 Q 2 Q 3 G ε
   3 i=1
Li-εmax 0

y i ϕ i + a(Ω 0 ; p, ϕ) -ϕ(L i -ε max ) ⊤ Λ ε p(L i -ε max ) = 3 i=1
Li-εmax 0 z i ϕ i , a(Ω 0 ; y, ϕ) -p 1 (0)ϕ 1 (0) -y(L i -ε max ) ⊤ Λ ε ϕ(L i -ε max ) = -ζϕ 1 (0), where ϕ ∈ H = ϕ i , ϕ ′ i ∈ L 2 (0, L i ) , ϕ 2 (0) = ϕ 3 (0) = 0, continuity at interior vertices. } and

Λ ε = 1 2ε -3   2 -1 -1 -1 2 -1 -1 -1 2   .
(80) Fig. 15 shows the shape functional with respect to ε. The derivative of J with respect to ε is consistently negative as ε approaches 0. This is the information that allows for the topology variations by nucleation of a small cycle at an interior vertex of the graph. Another observation emerges: as ε approaches ε 0 , the cost functional (J) converges to 0, indicating an optimal length for the introduced cycle. This implies that at ε = ε 0 , the network experiences an optimal configuration, emphasizing the critical nature of this parameter in shaping the network. Now, considering the dynamic example in Figure 14, where the damping parameter is set to α = 0, the cost functional is defined as: The state equation is given by:

J(u) =
       (y i ) tt -(y i ) xx + αy i = 0, t ∈ [0, T ], x ∈ [0, L i ], i = 1, 2, • • • , 6
, y i (0, x) = y 0 i (x), (y i ) t (0, x) = y 1 i (x), i = 1, 2, • • • , 6, (y 1 ) ′ (t, 0) = u(t), y 2 (t, 0) = 0, y 3 (t, 0) = 0, Continuity and Kirchhoff conditions.

Conclusions

In this paper, we exploit the properties of shape and topology optimization problems in one space dimension. The optimal control problems are considered in static and dynamic cases. For the tree network, we show the Turnpike Property for the wave equation and consider the geometric shape optimization. The new case is the tree with a small cycle. In such a case the topology of network is determined by using the topological derivatives obtained in static case. Numerically, the case of the cycle also does not pose a problem. This is new: Also here shape optimization works! This also holds for a graph with multiple cycles as long as in the shape optimization no topology change occurs, that is no cycle vanishes.

The Turnpike property also holds if a cycle disappears. We have shown this under the assumption that the initial state is not supported on the cycle. We expect that this condition is not sharp.

The convergence of gradient flow for shape optimization is also relevant to networks. The modeling of shapes of animals is another possibility for the spectral methods. We refer to [START_REF] Plotnikov | Geometric aspects of shape optimization[END_REF], and [START_REF] Chen | Animal shapes, modal analysis, and visualization of motion (I): Horse and camel[END_REF] for the related results.

Further research will concern the corresponding shape optimization problems in the general three dimensional case.

Acknowledgments

2 for n = 1 , 2 ,

 212 . . . with the eigenfunctions nontrivial only on the edges of the cycle ψ n (x) = sin nΠx ε . The singular behaviour of the branch {µ n } n of the spectrum means that for fixed n, µ n → ∞ with ε → 0+.

Fig. 1 :Fig. 2 :Fig. 3 :Fig. 4 :Fig. 5 :Fig. 6 :Fig. 7 :Fig. 8 :Fig. 9 :Fig. 10 :

 12345678910 Fig. 1: Optimal control for different values of T in example 1

Fig. 11 :Fig. 12 :

 1112 Fig. 11: The three-star graph

Fig. 13 :

 13 Fig. 13: Multiple perturbations of the three-star graph for domain decomposition technique.

Figure 14 .

 14 The variables z and ζ satisfy   -z ′′ i = 0, x ∈ [0, L i ], i = 1, • • • , 6, z ′ 1 (0) = ζ, z 2 (L 2 ) = z 3 (L 3 ) = 0,Continuity and Kirchhoff Condition.

Fig. 14 :

 14 Fig. 14: Domain decomposition for tripod directed network with an elementary small cycle.

Fig. 15 :

 15 Fig. 15: The shape functional for ε ∈ [0, 1] (Example for one cycle)

  t) -ζ) 2 dt.

.

  ε . The Steklov-Poincaré operator for the small, double cycle is:For numerical results, refer to Fig.17

  . The topological derivative at ε = 0 + is negative. And the optimal size of the cycle is ε = ε 0 = 0.5.

Fig. 17 :

 17 Fig. 17: The shape functional for the size of cycle ε ∈ [0, 1] (Example for multiple cycle).

  . This work was supported by Deutsche Forschungsgemeinschaft (DFG) in the Collaborative Research Centre CRC/Transregio 154, Mathematical Modelling, Simulation and Optimization Using the Example of Gas Networks, Projects C03 and C05, Projektnr. 239904186.Thus H(t) is differentiable almost everywhere and thus absolutely continuous on [0, T ] see[START_REF] Royden | Real Analysis[END_REF], Theorem 14. Almost everywhere on [0, T ] we have the derivativeH ′ (t) =

	A Differentiability properties
	Define the function
													H(t) =		0	t	u(s) sin (	π 2	+ nπ)	t -s t 0	ds
	and for h > 0
		Diff h H(t) =	H(t + h) -H(t) h
						=	1 h	t	t+h	u(s) sin (	π 2	+ nπ)	t + h -s t 0	ds
							+	1 h		0	t	u(s) sin (	π 2	+ nπ)	t + h -s t 0	-sin (	π 2	+ nπ)	t -s t 0	ds.
	Hence we have
	Diff h H(t) -	0	t	u(s)	π 2 + nπ t 0	cos (	π 2	+ nπ)	t -s t 0	ds
	≤	t	t+h	|u(s)|		sin ( π 2 + nπ) t+h-s t0 h	ds
		+	0	t	u(s)	 		sin ( π 2 + nπ) t+h-s t0	-sin ( π 2 + nπ) t-s t0 h	-	π 2 + nπ t 0	cos (	π 2	+ nπ)	t -s t 0	  ds
	≤	t	t+h	|u(s)|		π 2 + nπ t 0	ds
		+	0	t	u(s)	1 h			t	t+h	π 2 + nπ t 0	cos (	π 2	+ nπ)	τ -s t 0	-cos (	π 2	+ nπ)	t -s t 0	dτ ds
	≤	π 2 + nπ t 0	t	t+h	|u(s)| ds + h	( π 2 + nπ) 2 0 t 2	0	t	|u(s)| ds.
														π 2 + nπ t 0	0	t	u(s) cos (	π 2	+ nπ)	t -s t 0	ds.

The adjoint equation is:

Continuity and Kirchhoff conditions, (82) where χ(x) is the characteristic function of [0, L i -ε max ], i.e., χ(x) = 1 on [0, L i -ε max ] and χ(L i -x) = 0.

The optimal control u(t) is obtained from the optimality system as

This formulation provides insights into the dynamical optimization of the system. In Fig. 16, we illustrate the temporal evolution of the shape functional across varying ε ∈ [0, 1]. Subfigures (a), (b), and (c) correspond to T = 1, T = 10, and T = 100, respectively, offering insights into the functional's behavior over different time scales. It is evident from the results that the minimization of the cost functional consistently occurs when ε = ε 0 , emphasizing the pivotal role of this parameter in optimizing the system's response. Example for multiple cycles: In this example, the geometry is depicted in Figure 13. Here,

). The parameters, cost functional, and optimality system remain consistent with the Example for one cycle (static), with the sole distinction being the form of This implies in particular that α ′ n (t) ∈ L 2 (0, T ). Theorem 3.2 in [START_REF] Gugat | Optimal Boundary Control and Boundary Stabilization of Hyperbolic Systems[END_REF] contains a traveling wave representation of the state.