Share This Article:

Mathematical Analysis of an Optimal Control Problem of Surface Water Pollution

Full-Text HTML XML Download Download as PDF (Size:471KB) PP. 164-179
DOI: 10.4236/am.2017.82014    827 Downloads   1,035 Views  

ABSTRACT

We present in this paper a new technique based on Gelfand’s triplet [1] and include differential theory to make a theoretical analysis of an optimal control problem with constraints governed by coupled partial differential equations. This technique allowed us to give some theoretical results of existence and uniqueness of the solution of constraints and characterize the optimal control.

1. Introduction

The objective of this paper is to make a theoretical analysis of an optimal control problem with constraints governed by coupled partial differential equations.

The interest of this work is two-fold. Indeed firstly the resolution of optimal control problem with partial differential equations constraints is a challenge for the current research in the field both theoretically and numerically. Secondly, the mathematical model studied is of utmost importance in practice since the physical phenomenon is modeled through this paper concerning the pollution of surface water. The studied problem is very complex since the objective function is not explicitly defined in term of control directly. Hence there is a need to solve equations that govern constraints that are also complex as we couple the model which describes the pollutants dissolution with Navier-Stokes equations.

To overcome these difficulties, we built a new technique based on the Gelfand’s triplet and theory of differential solutions to transform the partial differential equations into ordinary differential equations [1] [2] .

The outline of this paper is as follows: in this second section we present the mathematical model; the third section is devoted to the main results of existence uniqueness and also estimation of the system; in the fourth section we conduct the characterization of the optimal control and we end this work with a conclusion and perspectives.

2. Mathematical Model

Let Ω be a bounded domain of 3 , with boundary Γ sufficiently regular. We denote by x = ( x , y , z ) an element of Ω .

We consider the following model:

min J ( v ) = 1 2 0 T Ω | D c ( v ) z d | 2 d x d t + N 2 0 T Ω | v | 2 d x d t , (1)

subject to

{ c t + i = 1 3 ( u i c ) x i d Δ c = g + v in Ω × ] 0 , T [ ( 2. a ) c = 0 on Γ × ] 0 , T [ ( 2. b ) c ( , 0 ) = c 0 , in Ω ( 2. c ) u t + ( u ) u ν Δ u + p = f in Ω × ] 0 , T [ ( 2. d ) div u = 0 in Ω × ] 0 , T [ ( 2. e ) u = 0 , on Γ × ] 0 , T [ ( 2. f ) u ( , 0 ) = u 0 in Ω ( 2. g ) (2)

where:

・ zd is an observation function; D is a given linear continuous operator of observation; N is a given positif real;

・ c denotes the concentration of pollutants, g is a source term, d is a diffusion coefficient, v represented the control which allows to act on the system;

u = ( u 1 , u 2 , u 3 ) and p = p ( x , t ) denote respectively the velocity and the pressure of water, f external forces, ν the kinematic viscosity;

・ Equations (2.a)-(2.c) modelling the transport and dissolution of pollutants;

・ Systems (2.d)-(2.g) is Navier-Stokes equations.

We rewrite problems (1)-(2) in the following compact form:

{ Find u U a d such that : J ( u ) = min v U a d J ( v ) (3)

where Uad denotes set of admissible controls. In this form it is easier to establish results of existence and uniqueness concerning the control when knowing some properties of J and Uad [3] [4] [5] .

3. Main Results

To make the theoretical study of system (2), we define the following Gelfand triplet:

V H V , where the injections are continuous, dense and compact.

We choose in this paper V = H 0 1 ( Ω ) , H = L 2 ( Ω ) and V = ( H 0 1 ( Ω ) ) = H 1 ( Ω ) .

We assume that there exist two positive constants A and A* such that the following inequality hold:

1 A V H A V . (4)

We also recall the following result concerning the Navier-Stokes equations whose proof can be found in [6] .

Theorem 1. Let u 0 belonging to L 2 ( Ω ) a divergence field. Then there exists a unique solution u of the Navier-Stokes equations associated to the initial data u 0 , such that

u L ( [ 0 , T ] ; L 2 ( Ω ) ) L 2 ( [ 0 , T ] ; H 0 1 ( Ω ) )

and

t 0, 1 2 u ( t ) L 2 ( Ω ) 2 + 0 T u ( θ ) L 2 ( Ω ) 2 d θ 1 2 u 0 L 2 ( Ω ) 2 (5)

3.1. Main Result 1

This first result concerning existence and uniqueness of solution of the model that governs the dissolution of pollutants:

{ c t + i = 1 3 ( u i c ) x i d Δ c = g + v on Ω × ] 0 , T [ c = 0 in Γ × ] 0 , T [ c ( , 0 ) = 0 , on Ω (6)

For this we introduce Gelfand’s triplet for the following operators:

Δ : = ( 2 x 2 0 0 0 2 y 2 0 0 0 2 z 2 ) : V V , (7)

and

U : = ( u 1 x 0 0 0 u 2 y 0 0 0 u 3 z ) : V H . (8)

We consider a function C : [ 0, T ] V defined by:

C ( t ) : = ( k c ( t , ) k c ( t , ) k c ( t , ) ) where k = 1 3 (9)

and

G = ( k ( g + v ) k ( g + v ) k ( g + v ) ) = G 0 + V (10)

With these notations Equation (6) becomes:

{ C ˙ ( t ) + U C ( t ) + d Δ C ( t ) = G ( t ) a . e t [ 0 , T ] C ( 0 ) = C 0 on H , (11)

where C ˙ denotes the derivative in the sense of the distributions of C and C 0 = ( k c 0 k c 0 k c 0 ) .

We define the space of solutions by:

W ( 0 , T ) = { w L 2 ( 0 , T ; H 0 1 ( Ω ) ) , w ˙ L 2 ( 0 , T ; H 1 ( Ω ) ) } , (12)

where L 2 ( 0 , T ; H 0 1 ( Ω ) ) and L 2 ( 0, T ; H 1 ( Ω ) ) are equipped with the respective norms:

f L 2 ( 0, T ; H 0 1 ( Ω ) ) = ( 0 T f ( t ) H 0 1 ( Ω ) 2 d t ) 1 2 (13)

and

f L 2 ( 0, T ; H 1 ( Ω ) ) = { 0 T ( sup ϕ H 0 1 ( Ω ) ( f ( t ) , ϕ ( t ) ) ϕ H 0 1 ( Ω ) ) 2 d t } 1 2 ϕ ( t ) 0. (14)

Theorem 2. If U L ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) , C 0 H and G L 2 ( 0, T ; H ) then there exists a unique state:

C W ( 0 , T ) L ( 0 , T ; H ) , (15)

solution of the problem (11).

Proof. To proof uniqueness we first establish the following intermediate result:

U L 2 ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) and C L 2 ( 0, T ; V ) then for all t [ 0, T ] :

( U ( t ) C ( t ) , C ( t ) ) V , V = 0. (16)

Indeed,

( U ( t ) C ( t ) , C ( t ) ) V , V = Ω U ( t ) C ( t ) C ( t ) = Ω i = 1 3 U i ( t ) C ( t ) x i C ( t ) = Γ i = 1 3 U i C ( t ) C ( t ) Ω i = 1 3 ( U i ( t ) C ( t ) x i + C ( t ) U i ( t ) x i ) C ( t ) = Ω U ( t ) . C ( t ) C ( t ) = ( U ( t ) . C ( t ) , C ( t ) ) V , V (17)

2 ( U ( t ) . C ( t ) C ( t ) ) V , V = 0 ( U ( t ) . C ( t ) , C ( t ) ) V , V = 0.

Assuming existence of two solutions C 1 and C 2 of (11), let denote by C = C 1 C 2 then C satisfies the following equation:

{ C ˙ ( t ) + d Δ C 1 ( t ) = U 2 C 2 ( t ) U 1 C 1 ( t ) a . e t [ 0 , T ] C ( 0 ) = 0 on H . (18)

by multiplying the first equation of (18) by C and integrating over [ 0, t ] , we obtain:

0 t ( C ˙ ( θ ) , C ( θ ) ) d θ + d 0 t ( Δ C ( θ ) , C ( θ ) ) d θ = 0 t ( U 2 C 2 ( t ) U 1 C 1 ( t ) , C ( θ ) ) d θ (19)

1 2 | C ( t ) | 2 1 2 | C ( 0 ) | 2 = d 0 t ( Δ C 1 ( θ ) Δ C 2 ( θ ) , C ( θ ) ) d θ 0 t ( U 2 . C 2 ( t ) , C 2 ( θ ) ) d θ + 0 t ( U 1 . C 1 ( t ) , C 1 ( θ ) ) d θ

using the monotony of Δ ( t , ) , we deduce that, ( Δ C 1 ( θ ) Δ C 2 ( θ ) , C ( θ ) ) 0 and using relation (16) we deduce that:

1 2 | C ( t ) | 2 1 2 | C ( 0 ) | 2 0 (20)

| C ( t ) | 2 | C ( 0 ) | 2 , (21)

since C ( 0 ) = 0 , we have | C ( t ) | = 0 C = 0 implying that C 1 = C 2 .

To prove existence of the solution we first prove that the problem approached of the problem model has a unique solution and that this approximate solution converges to the exact solution. Let consider the following approached problem ( P m ) of Equation (11) given by:

{ C ˙ m ( t ) + U C m ( t ) + d Δ C m ( t ) = G m ( t ) a . e t [ 0 , T ] C m ( 0 ) = C 0 m (22)

where G m ( : ) = i = 1 m ( G ( ) , e j ) e j and C 0 m = i = 1 m ( C 0 , e j ) e j .

We assume that this approximated problem is posed on the subspace V m = V e c t { e 1 , e 2 , , e m } of the space V which is separable, then there exists a sequence ( e m ) m such that:

m , The first terms e 1 , e 2 , , e n are linearly independent;

・ The set of linear combinations of elements of this sequence is dense in V (and therefore in H and V ).

Proposition 3. Under the assumptions of theorem (2), m 1 , (22) admits a unique solution C m which belongs to V m .

Proof. The proof is based on the Carathéodory’s theorem which we apply to the following function:

F m : [ 0 , T ] × V m m defined by

F m ( t , X ) = i = 1 m [ ( d Δ C m ( t ) U X ( t ) , e j ) + G m ( t ) ] . (23)

F m is measurable as Δ ( , v ) and U ( , v ) are measurable respectively in V and H , this function is also continuous in W [7] [8] ;

・ the condition of minoration on any compact of V m is assured by the fact that Δ ( t , ) is bounded a . e t [ 0, T ] in V and U ( t , ) is bounded in H a . e t [ 0, T ] .

Then, for every X 0 Ω , there exists at least one solution to the system:

{ X ˙ = F m ( t , X ) X ( 0 ) = X 0 (24)

which is therefore the solution of (22).

In order to demonstrate the convergence of the approximated solution to the exact solution, we establish the following priori estimation results.

Proposition 4. If C is a solution of (11) then:

C L 2 ( 0, T ; V ) 1 2 d | C 0 | + 1 d G L 2 ( 0, T ; V ) (25)

C L ( 0, T ; H ) | C 0 | + 2 d G L 2 ( 0, T ; V ) (26)

C ˙ L 2 ( 0 , T ; V ) [ d 2 | C 0 | + G L 2 ( 0 , T ; V ) ] [ A + 2 d U L 2 ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) ] + G L 2 ( 0, T ; V ) (27)

Proof. To prove relation (25) we form the scalar product of the first equation of Equation (11) with C and integrate over [ 0, T ] , we obtain:

( C ˙ ( t ) , C ( t ) ) V , V + ( U C ( t ) , C ( t ) ) V , V + d ( Δ C ( t ) , C ( t ) ) V , V = ( G ( t ) , C ( t ) ) V , V (28)

( C ˙ ( t ) , C ( t ) ) V , V + d ( Δ C ( t ) , C ( t ) ) V , V = ( G ( t ) , C ( t ) ) V , V .

Integrate over [ 0, T ]

0 T ( C ˙ ( t ) , C ( t ) ) V , V d t + d 0 T ( Δ C ( t ) , C ( t ) ) V , V d t = 0 T ( G ( t ) , C ( t ) ) V , V d t .

then

0 T ( C ˙ ( t ) , C ( t ) ) V , V d t = ( C ( T ) , C ( T ) ) L 2 ( Ω ) ( C ( 0 ) , C ( 0 ) ) L 2 ( Ω ) 0 T ( C ˙ ( t ) , C ( t ) ) V , V d t

0 T ( C ˙ ( t ) , C ( t ) ) V , V d t = 1 2 | C ( T ) | 2 1 2 | C ( 0 ) | 2 .

We have: Δ C V α 1 ( t ) + α 2 ( t ) C ( t ) V C ( t ) V :

0 T ( Δ C ( t ) , C ( t ) ) V , V d t 0 T Δ C V 2 d t 0 T C ( t ) V 2 d t = C L 2 ( 0 , T ; V ) 2

and we deduce the following successive inequalities

0 T ( G ( t ) + V ( t ) , C ( t ) ) V , V d t 0 T G ( t ) V C ( t ) V d t G L 2 ( 0 , T ; V ) C L 2 ( 0 , T ; V ) .

1 2 | C ( T ) | 2 1 2 | C ( 0 ) | 2 + d C L 2 ( 0, T ; V ) 2 G L 2 ( 0, T ; V ) C L 2 ( 0, T ; V ) .

d C L 2 ( 0 , T ; V ) 2 G L 2 ( 0 , T ; V ) C L 2 ( 0 , T ; V ) 1 2 | C ( 0 ) | 2 1 2 | C ( T ) | 2 1 2 | C ( 0 ) | 2

[ d C L 2 ( 0 , T ; V ) 1 2 d G L 2 ( 0 , T ; V ) ] 2 1 2 | C ( 0 ) | 2 + 1 4 d G L 2 ( 0 , T ; V ) [ 1 2 | C ( 0 ) | + 1 2 d G L 2 ( 0 , T ; V ) ] 2

d C L 2 ( 0 , T ; V ) 1 2 d G + V L 2 ( 0 , T ; V ) [ 1 2 | C ( 0 ) | + 1 2 d G L 2 ( 0 , T ; V ) ]

d C L 2 ( 0 , T ; V ) [ 1 2 | C ( 0 ) | + 1 d G L 2 ( 0 , T ; V ) ]

C L 2 ( 0 , T ; V ) 1 2 d | C 0 | + 1 d G L 2 ( 0 , T ; V ) .

To establish relation (26) we integrate the scalar product over [ 0, t ] [ 0, T ] :

0 t ( C ˙ ( s ) , C ( s ) ) V , V d s + d 0 t ( Δ C ( s ) , C ( s ) ) V , V d s = 0 t ( G ( s ) , C ( s ) ) V , V d s . (29)

1 2 | C ( t ) | 2 1 2 | C ( 0 ) | 2 + d C L 2 ( 0 , T ; V ) 2 G L 2 ( 0 , T ; V ) C L 2 ( 0 , T ; V ) .

1 2 | C ( t ) | 2 1 2 | C ( 0 ) | 2 d C L 2 ( 0 , T ; V ) 2 + G L 2 ( 0 , T ; V ) C L 2 ( 0 , T ; V ) 1 2 | C ( 0 ) | 2 + G L 2 ( 0 , T ; V ) C L 2 ( 0 , T ; V ) ,

Using (25) we have:

1 2 | C ( t ) | 2 1 2 | C ( 0 ) | 2 + G + V L 2 ( 0 , T ; V ) [ 1 2 d | C 0 | + 1 d G L 2 ( 0 , T ; V ) ] 1 2 | C ( 0 ) | 2 + 1 2 d | C 0 | G L 2 ( 0 , T ; V ) + 1 d G L 2 ( 0 , T ; V ) 2

| C ( t ) | 2 | C ( 0 ) | 2 + 2 d | C 0 | G L 2 ( 0 , T ; V ) [ 2 d G L 2 ( 0 , T ; V ) ] 2 [ | C ( 0 ) | + 2 d G L 2 ( 0 , T ; V ) ] 2

| C ( t ) | 2 | C ( 0 ) | + 2 d G L 2 ( 0 , T ; V ) ,

then we have

C L ( 0 , T ; H ) | C 0 | + 2 d G L 2 ( 0 , T ; V )

To proof relation (27) we use the fact that the Laplacian is a bounded linear operator and the following lemma:

Lemma 3.1. w [ H 0 1 ( Ω ) ] 3 and for all z H 0 1 ( Ω ) there exist two constants β and γ such as:

| w z | γ β w [ H 0 1 ( Ω ) ] 3 z H 0 1 ( Ω ) . (30)

Proof. For all w [ H 0 1 ( Ω ) ] 3 and for all z H 0 1 ( Ω ) , we have:

| w z | 2 = Ω | i = 1 3 w i z x i | 2 d x Ω | i = 1 3 w i | 2 | i = 1 3 z x i | 2 w [ L 2 ( Ω ) ] 3 2 z H 1 ( Ω ) 2 . (31)

According to Poincaré inequality there exist γ > 0 such as w [ L 2 ( Ω ) ] 3 γ w [ L 2 ( Ω ) ] 3 :

| w z | 2 γ 2 | w | [ L 2 ( Ω ) ] 3 2 z H 1 ( Ω ) 2

| w z | 2 γ 2 w [ H 0 1 ( Ω ) ] 3 2 z H 1 ( Ω ) 2

Norms H 0 1 and H 1 are equivalent in H 1 ( Ω ) , so that there exist a positive constant β such as:

| w z | 2 γ 2 β 2 w [ H 0 1 ( Ω ) ] 3 2 z H 0 1 ( Ω ) 2 ,

then we have

| w z | γ β w [ H 0 1 ( Ω ) ] 3 z H 1 ( Ω )

Lemma 3.2. If U L 2 ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) and C L 2 ( 0, T ; V ) then U C L 2 ( 0, T ; H ) . Furthermore

U C L 2 ( 0 , T ; H ) γ β U L 2 ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) C L 2 ( 0 , T ; V ) (32)

Proof. For all U L 2 ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) and for all C L 2 ( 0, T ; V ) , we have:

U C L 2 ( 0 , T ; H ) 2 = 0 T | U C | 2 d t . (33)

According to Lemma (3.1) we have:

U C L 2 ( 0 , T ; H ) 2 γ 2 β 2 0 T U ( t ) [ H 0 1 ( Ω ) ] 3 2 d t 0 T C ( t ) H 0 1 ( Ω ) 2 d t , U C L 2 ( 0 , T ; H ) γ β U [ H 0 1 ( Ω ) ] 3 2 d t C H 0 1 ( Ω ) , (34)

Consider ϕ L 2 ( 0, T ; V ) and we form the scalar product of the first equation of (11) with ϕ :

( C ˙ ( t ) , ϕ ( t ) ) V , V + ( U C ( t ) , ϕ ( t ) ) V , V + d ( Δ C ( t ) , ϕ ( t ) ) V , V = ( G ( t ) , ϕ ( t ) ) V , V

( C ˙ ( t ) , ϕ ( t ) ) V , V = ( U C ( t ) , ϕ ( t ) ) V , V d ( Δ C ( t ) , ϕ ( t ) ) V , V + ( G ( t ) , ϕ ( t ) ) V , V (35)

| 0 T ( C ˙ ( t ) , ϕ ( t ) ) V , V | d t d | 0 T ( Δ C ( t ) , ϕ ( t ) ) V , V d t | + | 0 T ( U . C ( t ) , ϕ ( t ) ) V , V d t | + | 0 T ( G ( t ) , ϕ ( t ) ) V , V | d t

| 0 T ( Δ C ( t ) , ϕ ( t ) ) V , V | d t 0 T | Δ C ( t ) | ϕ ( t ) Δ C V A C L 2 ( 0 , T ; V )

| 0 T ( U C ( t ) , ϕ ( t ) ) V , V d t | γ β U L 2 ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) C L 2 ( 0 , T ; V ) .

| 0 T ( U C ( t ) , ϕ ( t ) ) V , V d t | U L 2 ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) C L ( 0 , T ; H ) .

Finally we get

| 0 T ( C ˙ ( t ) , φ ( t ) ) V , V | d t [ d A C L 2 ( 0 , T ; V ) + U L 2 ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) C L ( 0 , T ; H ) + G L 2 ( 0 , T ; V ) ] ,

Using estimates (25) and (26), we obtain:

C ˙ L 2 ( 0 , T ; V ) A d 2 d | C ( 0 ) | + A d d G L 2 ( 0 , T ; V ) + U L 2 ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) + 2 d U L 2 ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) G L 2 ( 0 , T ; V ) + G L 2 ( 0 , T ; V ) ,

To achieve the proof of Theorem 2 it must be shown that the approximate solution ( C m ) converges to the exact solution in W ( 0, T ) . By construction, the result of initial states ( C 0 m ) m converges to C 0 in H .

Using the relation (25), we deduce that the terms of the sequence ( C m ) is bounded in L 2 ( 0, T ; V ) which is a reflexive Banach space. We can extract a subsequence that is weakly convergent in σ ( L 2 ( 0 , T ; V ) , L 2 ( 0 , T ; V ) ) .

Let denote by C ϕ ( n ) this subsequence where ϕ : is a function strictly increasing such that:

C ϕ ( n ) C 1 weakly in L 2 ( 0 , T ; V ) (36)

From relation (26) we deduce that the sequence ( C m ) are bounded in L ( 0, T ; H ) which is the dual of L 1 ( 0, T ; H ) . Then there exists a sub-sequence which converges weakly towards C 2 such as:

C ϕ ( n ) C 2 weakly in L ( 0, T ; V ) ; (37)

Using relation (27) we obtain the weak convergence C ˙ ϕ ( n ) to α L 2 ( 0, T ; V ) i.e:

C ˙ ϕ ( n ) α weakly in L 2 ( 0, T ; V ) (38)

According to [9] [10] , we deduce that C = C 1 = C 2 and C ˙ = α .

On the other hand let j fixed, m denoting by ( C m ) the unique solution of ( P m ) then when m + we obtain:

( U C m ( t ) d Δ C m ( t ) G m ( t ) , e j ) ( U C ( t ) d Δ C ( t ) G ( t ) , e j )

and

( C ˙ m ( t ) , e j ) ( C ˙ ( t ) , e j ) . (39)

Since the family { e j , j } is dense in V thus:

C ˙ ( t ) + U C ( t ) + d Δ C ( t ) = G ( t ) p . p t [ 0 , T ] .

Also to [10] , we can extract a subsequence noted again Cn such as p . p t [ 0, T ] , ( C n ( t ) ) converges to C ( t ) . Considering the continuous representative C , we obtain strong convergence ( C n ( 0 ) ) towards C ( 0 ) .

3.2. Regularity of the Solution

Let C ¯ be the unique solution of (11) and C ¯ + q a perturbed one of (11). We suppose that C 0 and λ H , F , G L 2 ( 0, T , V ) and U a , U b L ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) , then q satisfies the following equation:

{ q ˙ ( t ) + U a q ( t ) + U b ( C ¯ ( t ) + q ( t ) ) + d Δ q ( t ) = F ( t ) a . e t [ 0 , T ] q ( 0 ) = λ (40)

Effect the scalar product of (40) by q ( t ) and integrate over [ 0, t ] [ 0, T ] , we obtain:

1 2 | q ( t ) | 2 1 2 | λ | 2 + 0 t ( U a q ( s ) , q ( s ) ) d s + 0 t ( U b . ( C ¯ ( s ) + q ( s ) ) , q ( s ) ) d s + d 0 t ( Δ q ( t ) , q ( t ) ) = 0 t ( F ( s ) , q ( s ) ) d s . (41)

Δ being monotonous then:

1 2 | q ( t ) | 2 1 2 | λ | 2 + 0 t ( U b ( C ¯ ( t ) + q ( t ) , q ( t ) ) ) d s = 0 t ( F ( t ) , q ( t ) ) d s .

By Lemma 3.2, we have:

| q ( t ) | 2 | λ | 2 + 2 β γ 0 t U b [ H 0 1 ( Ω ) ] 3 C ¯ ( s ) + q ( s ) H 0 1 ( Ω ) q ( s ) H 0 1 ( Ω ) d s + 2 0 t F ( s ) H 1 ( Ω ) q ( s ) H 0 1 ( Ω )

Let:

η ( s ) = 2 β γ U b [ H 0 1 ( Ω ) ] 3 C ¯ ( s ) + q ( s ) H 0 1 ( Ω ) + F ( s ) H 1 ( Ω )

then

| q ( t ) | 2 | λ | 2 + 0 t η ( s ) q ( s ) H 0 1 ( Ω ) .

According to Willet-Wong’s inequality we deduce:

| q ( t ) | | λ | + 1 2 0 t η ( s ) d s ,

Using relation (25), one can have a priori estimation of q solution of (40). There exist a constant M λ , W , F as if | λ | , F L 2 ( 0 , T , H 1 ( Ω ) ) are increased then:

| q ( t ) | | λ | + β γ 0 t M λ , W , H U b [ H 0 1 ( Ω ) ] 3 d s + 0 t F ( s ) H 1 ( Ω ) d s

| q ( t ) | | λ | + β γ M λ , W , H U b L ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) + F ( s ) L 2 ( 0 , T ; H 1 ( Ω ) ) .

Let ρ = max { 1 , β γ M λ , W , H } , then:

| q ( t ) | ρ | λ | + U b L ( 0 , T ; [ H 0 1 ( Ω ) ] 3 ) + F ( s ) L 2 ( 0 , T ; H 1 ( Ω ) ) (42)

3.3. Main Result 2

The space of admissible controls considered in this paper is U a d = L 2 ( 0 , T ; H ) .

Theorem 5 The optimal control problem (3) has a unique solution W U a d And we have the following inequality: V U a d :

W V L 2 ( 0 , T ; H ) 4 N [ J ( V ) J ( W ) ] . (43)

Proof. to proof this result, one can proof that the functional J is inferior semi- continuous, strongly convex and differentiable.

Let us proof that J is inferior semi-continuous. We consider a sequence ( C n , V n ) in L 2 ( 0 , T ; H 0 1 ( Ω ) ) × L 2 ( 0 , T ; H ) which converges strongly to ( C , V ) , we can extract a subsequence also denoted by ( C n , V n ) such that:

( C n ( t ) , V n ( t ) ) ( C ( t ) , V ( t ) ) a . e in H 0 1 ( Ω ) × L 2 ( Ω ) . (44)

Let denote by L ( C ( t ) , V ( t ) ) = X ( C ( t ) ) + Y ( V ( t ) ) with:

X : H 0 1 ( Ω ) [ 0 , + [ C ( t ) X ( C ( t ) ) = Ω | D C ( v ) z d | 2 d x Y : L 2 ( Ω ) [ 0 , + [ V ( t ) Y ( V ( t ) ) = Ω | V | 2 d x

X and Y are inferior semi-continuous then:

L ( C ( t ) , V ( t ) ) lim inf n L ( C n ( t ) , V n ( t ) ) (45)

Apply Fatou’s lemma to the sequence L ( C n ( t ) , V n ( t ) ) we obtain:

0 T lim inf n L ( C n ( t ) , V n ( t ) ) d t lim inf n 0 T L ( C n ( t ) , V n ( t ) ) d t ,

according to (45) we have:

1 2 0 T Ω | D C ( v ) z d | 2 d x d t + N 2 0 T Ω | V | 2 d x , d t 1 2 0 T lim inf n Ω | D C n ( v ) z d | 2 d x d t + N 2 0 T lim inf n Ω | V n | 2 d x d t

J ( V ) lim inf n J ( V n ) .

So that J is inferior semi-continuous.

To proof the strong convexity, we choose W and V L 2 ( 0 , T ; H ) , then:

J ( W + V 2 ) = N 2 W + V 2 L 2 ( 0 , T ; H ) 2 + 1 2 0 T Ω | 1 2 ( D C ( V ) z d ) + 1 2 ( D C ( W ) z d ) | 2 d x d t . (46)

Applying the first unequally Clarkson to the first term in right member, we obtain:

J ( W + V 2 ) N 2 [ 1 2 W L 2 ( 0 , T ; H ) 2 + 1 2 V L 2 ( 0 , T ; H ) 2 W V 2 L 2 ( 0 , T ; H ) 2 ] + 1 2 0 T Ω | 1 2 ( D C ( v ) z d ) + 1 2 ( D C ( w ) z d ) | 2 d x d t

By applying the equality of the parallelogram to the second term the right member, we obtain:

J ( W + V 2 ) N 4 W L 2 ( 0 , T ; H ) 2 + N 4 V L 2 ( 0 , T ; H ) 2 N 8 W V L 2 ( 0 , T ; H ) 2 + 1 4 0 T Ω | D C ( W ) z d | 2 d x d t + 1 4 0 T Ω | D C ( V ) z d | 2 d x d t N 8 D C ( W ) D C ( V ) L 2 ( 0 , T ; H ) 2 .

J ( W + V 2 ) N 4 W L 2 ( 0 , T ; H ) 2 + N 4 V L 2 ( 0 , T ; H ) 2 N 8 W V L 2 ( 0 , T ; H ) 2 + 1 4 0 T Ω | D C ( W ) z d | 2 d x d t + 1 4 0 T Ω | D C ( V ) z d | 2 d x d t

J ( W + V 2 ) J ( W ) + J ( V ) 2 N 8 W V L 2 ( 0 , T ; H ) 2 . (47)

And the differentiability of J is proved by calculating the following limit:

lim θ 0 J ( V + θ W ) J ( V ) θ . (48)

Let us denoted by C the trajectory corresponding to V ( t ) + θ W ( t ) then C satisfies the following equation:

{ C ˜ ( t ) + U ˜ C ˜ ( t ) + d Δ C ˜ ( t ) = G 0 ( t ) + V ( t ) + θ W ( t ) a . e t [ 0 , T ] C ˜ ( 0 ) = C ˜ 0 in H . (49)

J ( V + θ W ) J ( V ) = 1 2 0 T Ω | D C ˜ z d | 2 d x d t + N 2 0 T | V + θ W | 2 d x d t 1 2 0 T Ω | D C z d | 2 d x d t N 2 0 T | V | 2 d x d t

J ( V + θ W ) J ( V ) = 1 2 0 T Ω [ | D C ˜ z d | 2 | D C z d | 2 ] d x d t + N 2 0 T [ | V + θ W | 2 | V | 2 ] d x d t = 1 2 0 T [ ( D C ˜ z d , D C ˜ D C + D C z d ) ( D C z d , D C z d ) ] d t + N 2 0 T [ ( V + θ W , V + θ W ) ( V , V ) ] d t = 1 2 0 T [ ( D C ˜ z d , D ( C ˜ C ) + D C z d ) ( D C z d , D C z d ) ] d t + N 2 0 T [ 2 θ ( V , W ) + θ 2 ( W , W ) ] d t = 1 2 0 T [ ( D C ˜ z d , D ( C ˜ C ) ) + ( D ( C ˜ C ) , D C z d ) ] d t + N 2 0 T [ 2 θ ( V , W ) + θ 2 ( W , W ) ] d t

Let us designed by

C ^ = lim θ 0 C ˜ C θ ,

then

lim θ 0 J ( V + θ W ) J ( V ) θ = 1 2 0 T [ ( D C ˜ z d , D C ^ ) + ( D C ^ , D C z d ) ] d t + N 0 T ( V , W ) d t

J ( V ) W = 1 2 0 T [ ( D C ˜ z d , D C ^ ) + ( D C ^ , D C z d ) ] d t + N 0 T ( V , W ) d t

Note that if θ 0 then C ˜ C

J ( V ) W = 0 T Ω ( D C z d ) D C ^ d x d t + N 0 T Ω V W d x d t (50)

According to these properties of J we conclude that problem (3) admits an unique minimum W on U a d .

From relation (47) we have the inequality:

W V L 2 ( 0 , T ; H ) 4 N [ J ( V ) J ( W ) ] .

As U a d = L 2 ( 0, T ; H ) then

J ( V ) W = 0 T Ω ( C z d ) C ^ d x d t + N 0 T Ω V W d x d t (51)

search C ^ , i.e. the tangent model: In fact ( 49 ) ( 11 ) θ , we obtain:

C ˜ C ˙ θ + U ˜ C ˜ U . C θ d Δ C ˜ C θ = w , (52)

by tender θ 0 , we obtain

{ C ^ + U C ^ d Δ C ^ = W C ^ ( , 0 ) ) = 0 , (53)

since U ˜ U when θ 0 .

3.4. Main Result 3

Theorem 6. If D is a surjective operator of the space L 2 ( 0, T ; V ) in L 2 ( 0, T ; V ) or D is the identity application of L 2 ( 0, T ; V ) then the control V is given by:

V = N 1 P (54)

where P is the solution of the adjoint equation.

First Case: We assume that D is an injective operator define from the space in L 2 ( 0, T ; V ) L 2 ( 0, T ; H ) and z d L 2 ( 0 , T ; V ) , then the adjoint equation is given by:

{ P ˙ U P d Δ P = C z d P ( T ) = 0 (55)

Seconde Case: D is the identity operator define from L 2 ( 0, T ; V ) L 2 ( 0, T ; V ) when z d L 2 ( 0, T ; V ) then adjoint equation is given by:

{ P ˙ U P d Δ P = C z d P ( T ) = 0 (56)

when z d L 2 ( 0, T ; H ) the adjoint equation becomes:

{ P ˙ U P d Δ P = ( Δ x + I ) ( C z d ) P ( T ) = 0 (57)

Proof. The derivative given by (51) is unusable since for each test W we must solving the tangent model (53). The introduction of the adjoint state allows us to obtain a explicit expression of J [11] [12] .

We first perform the scalar product of the first equation of (53) with P and then integrate by parts we obtain:

0 T ( C ^ ( t ) , P ) d t = C ^ ( T ) P ( T ) C ( 0 ) P ( 0 ) 0 T ( P ( t ) , C ^ ( t ) ) d t (58)

0 T ( U C ^ ( t ) , P ) d t = 0 T Ω ( U C ^ ) P d t = 0 T Γ U C ^ P d t 0 T Ω ( U P ) C ^ d x d t = 0 T ( U P , C ^ ) d x d t

0 T ( Δ C ^ , P ) d t = 0 T Ω Δ C ^ P d x d t = 0 T Γ C ^ P d x d t 0 T Ω C ^ P d x d t = 0 T Γ C ^ P d x d t + 0 T Ω Δ P C ^ d x d t .

with additional conditions on P given by:

P ( x , t = T ) = 0 , P ( x , t ) = 0 x Γ

we obtain:

0 T Ω ( P ˙ U P d Δ P ) C ^ d x d t = 0 T Ω W P d x d t

By identifying this result with the first term of J ( V ) W , one obtains the adjoint equation according to the nature of the operator D given by (55)-(57)

To characterize the control we consider the two cases given by relations (56), and (57). Then if the adjoint equation is given by the proof is the same for both cases D . The adjoint equation is given by:

{ P ˙ U P d Δ P = ( C z d ) P ( T ) = 0

In this case relation (51) becomes:

J ( V ) W = 0 T Ω ( C z d ) C ^ d t + N 0 T Ω V W d x d t (59)

We form the scalar product between the first equation of the adjoint equation with C ^ and we integrate over [ 0, T ] , we obtain:

0 T ( P ˙ ( t ) , C ^ ) d t = C ( T ) P ( T ) C ( 0 ) P ( 0 ) 0 T ( C ^ ( t ) , P ( t ) ) (60)

0 T ( U P ( t ) , C ^ ) d t = 0 T Ω ( U P ) C ^ d t = 0 T Γ U P C ^ d t 0 T Ω ( U C ^ ) P d x d t = 0 T ( U C ^ , P ) d x d t

0 T ( Δ P , C ^ ) d t = 0 T Ω Δ P C ^ d x d t = 0 T Γ P C ^ d x d t 0 T Ω P C ^ d x d t = 0 T Γ P C ^ d x d t + 0 T Ω Δ C ^ P d x d t .

We have therefore:

0 T ( P ˙ U P d Δ P , C ^ ) d t = 0 T ( P , C ^ ( t ) + U C ^ + Δ C ^ ) d t = 0 T ( P , W ) d t = 0 T Ω P W d x d t .

So that thus one can write J as follows (51) give:

J ( V ) W = 0 T Ω ( P ˙ U P d Δ P ) C ^ d x d t + N 0 T Ω V W d x d t

J ( V ) W = 0 T Ω P W d x d t + N 0 T Ω V W d x d t

J ( V ) W = 0 T Ω ( P + N V ) W d x d t = 0 (61)

His equation had to be satisfied whatever the disturbance on W , especially if W 0 . Necessarily, we have:

P + N V = 0. (62)

4. Conclusions

The aim of this paper is a mathematical analysis of an optimal control problem of surface water pollution by using a triplet of evolution adapted to order of derivations. This technique allowed us to propose theorems not only of existence and uniqueness of the control but also of the solution of the equation that governs the constraints which is nonlinear. Then we solved the control problem directly and gave an optimal control characterization by using the adjoint equation.

It is essential to remember that this study only allowed us to characterize the optimal control without giving an analytical expression of it. In a future work we will proceed to a numerical approximation of this model in order to propose solution to the decision makers.

Cite this paper

Moustapha, D. , Haoua, H. and Bisso, S. (2017) Mathematical Analysis of an Optimal Control Problem of Surface Water Pollution. Applied Mathematics, 8, 164-179. doi: 10.4236/am.2017.82014.

References

[1] Gelfand, I.M. and Ya, N. (1964) Vilenkin; Generalized Functions. Applications of Harmonic Analysis. Vol. 4, Academic Press, New York.
[2] Pyatnitsky, E.S. and Rapoport, L.B. (1996) Criteria of Asymptotic Stability of Differential Inclusions and Periodic Motions of Time-Varying Nonlinear Control Systems. IEEE Transactions on Circuits and Systems, 43, 219-229. https://doi.org/10.1109/81.486446
[3] Ciarlet, P.G. (1982) Introduction à l’analyse numérique matricielle et à l’optimisation. Masson, Paris.
[4] Culioli, J.C. (1994) Introduction à l’optimisation. Editions Ellipses, Paris.
[5] Faurre, P. (1988) Analyse numérique. Notes d’optimisation, Ecole Polytechnique. Ellipse, Edition Marketing, Paris.
[6] Cannone, M. (2000) Nombres de Reynolds, stabilité et Navier-Stokes. Vol. 52, Banach Center Publications, Warszawa.
[7] Brezis, H. (1968) Equations et inéquations non linéaires dans les espaces vectoriels en dualité. Annales de l’Institut Fourier (Grenoble), 18, 115-175.
[8] Bonnans, J.F. (1998) Second Order Analysis for Control Constrained Optimal Control Problems of Semilinear Elliptic Systems. Applied Mathematics & Optimization, 38, 303-325.
https://doi.org/10.1007/s002459900093
[9] Aubin, J.P. (1963) Un théorème de compacité. Comptes Rendus de l’Académie des Sciences, Paris, 256.
[10] Clérin, J.-M. (2009) Problèmes de controle optimal du type bilinéaire gouvernés par des équations aux dérivées partielles d’évolution. Thèse, LANLG, Université d’Avignon.
[11] Lions, J.L. (1968) Controle optimal de systèmes gouvernés par des équations aux dérivées partielles. Dunod, Paris.
[12] Lions, J.-L. (1969) Quelques méthodes de résolution des problèmes aux limites non linéaires. Gauthier-Villars, Paris.

  
comments powered by Disqus

Copyright © 2017 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.