wavebnd2 details:

Keywords: analytic solution

Global classification: linear-quadratic, convex

Functional: convex quadratic

Geometry: easy, fixed

Design: coupled via boundary values 0th order

Differential operator:

Design constraints:

State constraints:

Mixed constraints:



Submitted on 2013-03-01 by Martin Gugat. Published on 2013-03-11


wavebnd2 description:


Introduction

Even for smooth data, the optimal state in boundary control problems for the wave equation is in general discontinuous. The continuity of the optimal state can be formulated as an additional requirement. In Gugat [2006], a boundary control problem for the 1D wave equation is considered, with a constraint that the state is at rest at terminal time T. The continuity of the state is achieved by imposing conditions which require the compatibility of the boundary controls with the initial data. Moreover, the objective function (the sum of the L2 norms of the time derivatives of the controls) requires the continuity in time of the controls.

This problem and analytical solution appear as [Gugat2006, Example 13.2].

Variables & Notation

Unknowns

u0 H1(0,T) control at the boundary x = 0 uL H1(0,T) control at the boundary x = L y H1((0,L) × (0,T))state variable

Given Data

The given data is chosen in a way which admits an analytic solution.

L = 1 length of spatial domain (0,L) T = 2 final time y0(x) = 1initial elongation at time zero y1(x) = 0 initial velocity at time zero

Problem Description

Minimize0T u 0(s)2 + u L(s)2ds s.t. ytt(x,t) = yxx(x,t)in (0,L) × (0,T) y(x,0) = y0(x) in (0,L) yt(x,0) = y1(x) in (0,L) y(0,t) = u0(t) in (0,T) y(L,t) = uL(t) in (0,T) y(x,T) = 0 in (0,L) yt(x,T) = 0 in (0,L) and y(0,0) = u0(0) y(L,0) = uL(0) y(0,T) = u0(T) y(L,T) = uL(T).

The state is forced to rest at terminal time T. The final block of constraints represents the compatibility conditions which lead to the continuity of the state, together with the continuity of the controls induced by the objective.

Supplementary Material

The optimal controls are given by

u0(t) = uL(t) = 1 + t 2,t [0,T].

Figure 0.1 shows two views of the optimal state.


PIC PIC

Figure 0.1: Optimal state viewed from different angles as given in [Gugat2006, Figure 1(a) and 1(b)]

The state can be obtained from [Gugat2006, Section 11.2.3] as

y(x,t) = 1 2α(x + t) + β(x t)

with the functions

α(s) = 1 s [0,1] s 2s [1,2] 0 s [2,3] β(s) = 0 s [2,1] s 1s [1,0] 1 s [0,1].

References

   M. Gugat. Optimal boundary control of string to rest in finite time with continuous state. Zeitschrift für Angewandte Mathematik und Mechanik, 86:134–150, 2006. doi: 10.1002/zamm.200410236.