Initial Value Methods for Parabolic Control Problems - American ...

3 downloads 120 Views 844KB Size Report
boundary value control and where the observation is the final state. The methods ... (1.2) are generated by solving Unear initial value problems of the form (1.1).
MATHEMATICS OF COMPUTATION, VOLUME 34, NUMBER 149 JANUARY 1980, PAGES 115-125

Initial Value Methods for Parabolic Control Problems By Ragnar Winther Abstract.

We study iterative

boundary

value control

based on transforming control) straints).

tial value problems

for parabolic

control

problems

is the final state.

with a Neumann The methods

the original control problem (which may have constraints

into an equivalent The methods

methods

and where the observation

problem of minimizing a strictly convex functional

are semidiscrete

are

on the (no con-

in the sense that we assume that parabolic

ini-

can be solved exactly.

1. Introduction. The purpose of this paper is to study approximations of certain parabolic control problems. In order to describe these problems, let £2 be a bounded domain in Rd with sufficiently smooth boundary, 8Í2, and for a fixed T > 0

let Ô = (0, T) x n and 2 = (0, T) x 3Í2. On the domain £2 let L denote the second order differential operator

Lu = -

i,7=i dxi \

bxjJ

We will assume that ai • and c are sufficiently smooth real-valued functions on

Í2 and that the operator L is strongly elUptic; i.e., there is a constant cx > 0 such

that í,7= i

i'= l

for all x G Í2 and %G Rd. Now let i>! G L2(ü) and /G L2(Q) be given and for any g G ¿2(E), let ug = u (t, x) denote the unique weak solution of the problem

(1.1)

du/dt + Lu = f

on Q,

du/dv=g

on 2,

«(0, •) = vx

on il.

Here d/dv = Ef •_, a¡ n¡(blbxA, where n¡ is the ith component of the outward unit

normal on 8Í2. We shall consider optimal control problems of the form

(1.2)

min{\\ugiT,-)-v2\\l2(n)+a\\g\\2L2

},

g£K

Received November 27, 1978. AMS (MOS) subject classifications

(1970).

Primary 49D05. © 1980 American Mathematical 0025-571 8/80/0000-0006/$03.75

115 License or copyright restrictions may apply to redistribution; see http://www.ams.org/journal-terms-of-use

Society

116

RAGNAR WINTHER

where v2 G Z,2(S2) and a > 0 are given and where A' is a closed convex set of L2(Z).

In Section 2 we shaU consider the case when K = L2ÇZ), while in Section 4 we con-

sider K = {g&L2ÇZ)\g >0 a.e. on Z}. In both cases we shall derive algorithms where approximations of the solution of (1.2) are generated by solving Unear initial value problems of the form (1.1). This will

be done by showing that (1.2) (in both cases mentioned above) is equivalent to minimizing a strictly convex functional over the entire space L2(£l).

The minimum of

this functional can then be found by standard iterative techniques.

In the case K =

Z,2(2) the associated functional is quadratic, and in Section 3 we study how this functional can be minimized by the conjugate gradient method. We shall particularly show certain superlinear convergence estimates for this method by applying some results recently given by the author in [13]. In order to evaluate the functionals mentioned above we have to solve two initial value problems of the form (1.1). Hence, the approximation of (1.2) is reduced

to the well-studied problem of solving problems of the form (1.1) numerically. Such numerical methods will not be considered here. We shall instead assume that problems of the form (1.1) can be solved exactly.

However, we mention that a detailed study of fuUy discrete analogs of the methods considered here was done by the author in [11] and [12] in the case when K = L2(Z), and some numerical examples can be

found in [11]. 2. Preliminaries. We start with some notation. If Hx and H2 axe two Hubert spaces then l(//j, H2) will denote the set of bounded linear maps from Hx to H2. On the spaces L (Í2) and ¿2(9Í2), respectively, we shall use the notation

(, i/O= Jn dx

and

(= jan y\pdo

for the inner products, and the associated norms wiU be denoted by II• II and I-I. If p > 0, let /^(íí) denote the Sobólev space of order p of real-valued functions on £2 (see for example Lions and Magenes [5]). If p < 0, //p(f2) is defined to be the dual of/T~P(Í2) with respect to the inner product of L2(Í2).

Furthermore,

W(0,T) = J/l/GL2(0, T-H'im. JteL2(*°>^/T1^))!

we let

.

where d/dt is taken in the sense of distributions on (0, T) with values in //'(Í2).

We

recaUfrom [5] that W(0, T) C C(0, T; ¿2(S2)). Let also B: HliQ,) x H1^)

—►R denote the biUnear form associated with the

operator L; i.e.,

r

[

"

dip d\b

dx.

The problem (1.1) can now be written in the variational form

(2.1)

(idu/dt, 0 such that \\F(zx) - F(z2)\\ < (5.5)

sup of1 i7" \E(t)(zx - z2)\ \E(t^p\dt il^ll=i "