SOLVING A CROP PROBLEM BY AN OPTIMAL CONTROL METHOD

Transcription

1 SOLVING A CROP PROBLEM BY AN OPTIMAL CONTROL METHOD Maïtine Bergounioux, Hem Raj Joshi, Suzanne Lenhart To cite this version: Maïtine Bergounioux, Hem Raj Joshi, Suzanne Lenhart. SOLVING A CROP PROBLEM BY AN OPTIMAL CONTROL METHOD. Natural Resource Modeling, Rocky Mountain Mathematics Consortium, 2005, 8, number 3 pp <hal > HAL Id: hal Submitted on 3 Mar 2006 HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

2 Solving a Crop Problem by an Optimal Control Method Hem Raj Joshi Suzanne Lenhart Maïtine Bergounioux (Draft) Last Update : May 0 th, 2004 ccsd , version - 3 Mar 2006 Résumé A system of ordinary differential equations coupled with a parabolic partial differential equation is studied in order to understand an interaction between two crops and a pathogen. Two different types of crops are planted in same field in some pattern so that the spread of pathogen can be controlled. The pathogen prefers to eat one crop. The other crop, which is not preferred by pathogen, is introduced to control the spread of pathogen in the farming land. The optimal initial planting pattern is sought to maximize plant yields while minimizing the variation in the planting pattern. The optimal pattern is characterized by a variation inequality involving the solutions of the optimality system. Numerical examples are given to illustrate the results. Key Words : Crop Problem, Optimal Control, Pathogen Spread AMS Classification : 35K55, 49K20, 92D25 University of Tennessee, Department of Mathematics, Knoxville, TN Départment de Mathématiques, Université d Orléans, Laboratoire MAPMO, UFR Sciences, BP 7659, Orléans, France.

3 Crop Problem 2 Introduction We consider a system with two ordinary differential equations (ODEs) and one parabolic partial differential equation (PDE) modeling the planting of two different types of crops (u and v) in same field in some pattern so that the spread of pathogen (w) can be controlled. The pathogen prefers to eat crop u. The crop v, which is not preferred by pathogen, is introduced to control the spread of pathogen in the farming land (Ω = (0,)). The state system is the following : du dt = r u(k u) k 3 uw in Q = Ω (0,T) dv dt = r 2 v(k 2 v) k 4 vw w t = d w xx + α k 3 uw + α 2 k 4 vw µw with initial and boundary conditions : u(x,0) = u 0 (x), v(x,0) = v 0 (x) = a u 0 (x) w(x,0) = w 0 (x) for x Ω w(0,t) = 0 = w(,t) The coefficients and terms can be interpreted as : u(x,t) v(x,t) w(x,t) r,r 2 d k,k 2 α,α 2,k 3,k 4 = crop preferred by pathogen (first state variable). = crop not preferred by pathogen (second state variable). = pathogen (third state variable). = growth rates. = diffusion coefficient. = the carrying capacities. = interaction coefficients. µ = pathogen natural death rate. The control set is U = {u 0 (x) H 0(0,) 0 u 0 (x) a} on Ω (0,T). We seek to maximize the objective functional over u 0 U : [ J(u 0 ) = (A u + A 2 v)(x,t) ( (u 2 0 ) 2 + B u B 2 (a u 0 ) 2)] dx. (.) 0 The positive constants A and A 2 represent the relative importance of the terms u and v respectively and B and B 2 are multipliers of the cost of implementing the control. Minimizing the u 0 term represents low variation in u 0. A planting pattern with high

4 Crop Problem 3 variation would be unrealistic to implement. The goal would be to maximize plant yields (subject to relative importance of the two crops) while minimizing the variation in the initial planting pattern. Intercropping for weed and pest management has been considered in a variety of contexts [5, 8, 9], but mostly in the setting of systems of ODEs differential equations or difference equations. Here we have a combination of a parabolic PDE for diffusion of the pathogen with ODEs for the crops. Due to inclusion of the u variation term in the objective functional, our characterization of the optimal planting pattern is a variational inequality [7], instead of simply an algebraic expression in terms of the state and adjoint variables. Such a variational inequality is somewhat novel in control of PDE problems and requires an unusual numerical algorithm. In section 2, we discuss the existence of an optimal control, i.e. the optimal planting pattern. The optimality system, which characterizes the optimal control, is derived in section 3. The optimality system involves the state system and the adjoint system together with the characterization of the optimal control given by a variational inequality. Section 4 treats the uniqueness of the optimal control by obtaining the uniqueness of the solutions of the optimality system for T sufficiently small. Finally in section 5, we discuss our numerical algorithm and illustrate numerical examples for our problem. 2 Existence of an Optimal Control The following assumptions are made throughout this paper.. α 2 k 4 < µ 2. r,r 2,d,k 3,k 4,α,α 2,µ are positive constants and α k 3 > α 2 k < a. Assumption means that the u 0 would cause the pathogen to decay when it only eats the second crop. Assumption 2 means that the consumption of crop u contributes more to the growth of the pathogen than would comsumption of crop v. The underlying solution space for system (.) is V = (L 2 (Q)) 2 L 2 (0,T;H 0 (Ω)). Definition 2. We say a triple of functions u,v,w V, with w t L 2 (0,T;H (Ω)), is a weak solution of (.) with given boundary and initial conditions provided u(x,t) = u 0 (x) + t 0 [r u(k u) k 3 uw] (x,s)ds v(x,t) = v 0 (x) + t 0 [r 2v(k 2 v) k 4 vw] (x,s)ds T 0 < w t,φ > dt + d Q u φdxdt

5 Crop Problem 4 = Q (α k 3 uw + α 2 k 4 vw µw) φdxdt for all φ L 2 (0,T;H 0 (Ω)), where the <, > inner product is the duality between H (Ω) and H0 (Ω), and u(x,0) = u 0 (x), v(x,0) = v 0 (x) w(x,0) = w 0 (x), w(0,t) = 0 = w(,t) (2.) Remark : Since u,v C(0,T;L 2 (Ω)) from Evans [], the initial conditions (2.) make sense. We will prove an existence and uniqueness result for the state system (.). This result will be established in Theorem 2.. Theorem 2. Given u 0 U, there exists a unique solution (u,v,w) in V solving (.) and (2.). Proof : We consider the following problems du dt = r k U U(x,0) = u 0 (x) dv dt = r 2 k 2 V V (x,0) = v 0 (x) W t = d W xx + α k 3 UW + α 2 k 4 V W µw W(x,0) = w 0 (x) W(0,t) = 0 = W(,t) for x Ω for x Ω for x Ω on Ω (0,T). The functions U,V,W are supersolutions of system equations in (.), which are L bounded in Q. To obtain the existence, we will construct three sequences by means of iteration, using the above supersolutions. Define : u = U,v = V,w 2 = W,u 0 = 0,v 0 = 0,w = 0, where the superscripts denote the iteration step. For i = 2,3,, we define u i,v i and w i+ as the solution of the following problems respectively : u i t + Rui = f(u i 2,w i ) in Q u i (x,0) = u 0 (x) for x Ω vt i + Rv i = g(v i 2,w i ) in Q v i (x,0) = v 0 (x) for x Ω wt i d wxx i + Rwi = h(u i,v i,w i 2 ) in Q w i (x,0) = w 0 (x) for x Ω w i (0,t) = 0 = w i (,t) on Ω (0,T).

6 Crop Problem 5 where f(u i 2,w i ) = Ru i 2 + r u i 2 ( k u i 2) k 3 u i 2 w i g(v i 2,w i ) = Rv i 2 + r 2 v i 2 ( k 2 v i 2) k 4 v i 2 w i h(u i,v i,w i 2 ) = Rw i 2 + α k 3 u i w i 2 + α 2 k 4 v i w i 2 µw i 2. R is a constant satisfying sup((k 3 + k 4 + µ) W + 2r U + 2r 2 V ). Q Since (2.2),(2.2) and (2.2) are linear problems, the solutions u i,v i,w i, for i =,2,, exist. Claim : For i =,2,, 0 u i U, 0 v i V, 0 w i W in Q. Claim can be proved by induction using the maximum principle. Note that for u i 0,v i 0,w i 0, i =,2,, Claim 2 : f is increasing in u i 2, and decreasing in w i. g is increasing in v i 2, and decreasing in w i. h is increasing in u i,v i and decreasing in w i 2. There exists u,v,w,u,v,w in V such that the following monotone pointwise convergence holds. Proof of Claim 2 : u 2i ր u, u 2i+ ց u in Q v 2i ր v, v 2i+ ց v in Q w 2i ց w, w 2i+ ր w in Q. To prove the convergence we use the induction method. Since w = 0, we have (w 5 w 3 ) t d (w 5 w 3 ) xx + R(w 5 w 3 ) = w 3 (R + α k 3 u 4 + α 2 k 4 v 4 µ) 0 and thus 0 w w 3 w 5. Similarly we have u 3 u,v 3 v, w 2 w 4,u 0 u 2 and v 0 v 2. Also u 2 u 3,v 2 v 3 and w 4 w 3 since (u 2 u 3 ) t + R(u 2 u 3 ) = u (R + r (k u ) k 3 w 3 ) 0 (v 2 v 3 ) t + R(v 2 v 3 ) = v (R + r 2 (k 2 v ) k 4 w 3 ) 0 (w 4 w 3 ) t d (w 4 w 3 ) xx + R(w 4 w 3 ) = w 2 (R + α k 3 u 3 + α 2 k 4 v 3 µ) 0. Fix i, assume that for all j i such that u 2j ր v 2j ր u 2j+ ց v 2j+ ց w 2j ց w 2j+ ր.

7 Crop Problem 6 Now we compare w 2i+ and w 2i. From (2.2), we have (w 2i+2 w 2i ) t d (w 2i+2 w 2i ) xx + R(w 2i+2 w 2i ) = h(u 2i+,v 2i+,w 2i ) h(u 2i,v 2i,w 2i 2 ) 0. Then w 2i+2 w 2i. Other cases can be proved using similar arguments. Hence by boundedness of u i,v i,w i and monotone properties of f,g,h, we get the pointwise convergence. Claim 3 : The subsequences of {u i }, {v i }, {w i } satisfy : u 2j u, u 2j+ u, v 2j v, v 2j+ v, w 2j w, w 2j+ w weakly in V. Proof of Claim 3 : Since RHS of (2.2), (2.2) and (2.2) are bounded in L (Q), then the sequences {( u 2i,v 2i,w 2i)}, {( u 2i+,v 2i+,w 2i+)} are uniformly bounded in V. Hence using the weak compactness of the sequences in V, we get u 2j u, u 2j+ u, v 2j v, v 2j+ v, w 2j w, w 2j+ w weakly in V. Since we have pointwise convergence on each sequence by claim 2, this weak convergence is also on the whole (even or odd) sequence (not just on a subsequence). Claim 4 The subsequences of {u i }, {v i }, {w i } satisfy u 2j u, u 2j+ u, v 2j v, v 2j+ v, w 2j w, w 2j+ w strongly in L 2 (Q). Proof of Claim 4 : Using the V boundedness on u 2i,u 2i+,v 2i,v 2i+,w 2i,w 2i+ and the system in (2.2)- (2.2), we obtain that w 2i t, w2i+ t are bounded in L 2 (0,T;H (Ω)). Hence, using weak compactness again, we have w 2i t w t, w 2i+ t w t weakly in L 2 (0,T;H (Ω)). Using a compactness result by Simon [0] for w and, we have w 2i w, w 2i+ w strongly in L 2 (Q).

8 Crop Problem 7 Since equations (2.2), (2.2) are ODEs in time with parameters x and bounded RHS, the sequence is uniformly Lipschitz in t for each x, we can pass to limit in u,v by using weak solution definition u 2i u, u 2i+ u, v 2i v, v 2i+ v pointwise uniformly in t, for each x Now we prove the uniqueness, i. e., u = u,v = v and w = w. Passing to the limit in the u 2i,u 2i+,v 2i,v 2i+,w 2i, and w 2i+ PDE s, we obtain : u t = r k u r u 2 k 3 uw in Q u(x,0) = u 0 (x) for x Ω u t = r k u r u 2 k 3 uw in Q u(x,0) = u 0 (x) for x Ω v t = r 2 k 2 v r 2 v 2 k 4 vw in Q v(x,0) = v 0 (x) for x Ω v t = r 2 k 2 v r 2 v 2 k 4 vw in Q v(x,0) = v 0 (x) for x Ω w t d w xx = α k 3 uw + α 2 k 3 vw µw in Q w(0,t) = 0 = w(,t) on Ω (0,T) w(x,0) = w 0 (x) for x Ω (2.2) (2.3) (2.4) (2.5) (2.6) w t d w xx = α k 3 uw + α 2 k 3 vw µw in Q w(0,t) = 0 = w(,t) w(x,0) = w 0 (x) on Ω (0,T) for x Ω (2.7) Let u = e λt f, u = e λt f, v = e λt g, v = e λt g, w = e λt h and w = e λt h, where λ > 0 is to be chosen. To illustrate the transformed system, we write equation (2.2) : f t + λf = r k f r e λt f 2 k 3 e λt fh in Q (2.8) We consider the weak formulation of the f f, g g and h h problems, and after adding both weak formulations, we obtain on Q = Ω (0,t ) : Q {(f f) t (f f) + λ(f f) 2 + (g g) t (g g) + λ(g g) 2 }dxdt + Q {(h h) t (h h) + d (h h) 2 + λ(h h) 2 }dxdt = Q {r k (f f) 2 r e λt (f 2 f 2 )(f f) k 3 e λt (f f)(fh f h)}dxdt + Q {r 2 k 2 (g g) 2 r 2 e λt (g 2 g 2 )(g g) k 4 e λt (g g)(gh g h)}dxdt + Q {α k 3 (fh fh)(h h) + α 2 k 4 (gh gh)(h h) µ(h h) 2 }dxdt.

9 Crop Problem 8 We obtain 2 Ω {[f f]2 (x,t) + [g g] 2 (x,t) + [h h] 2 (x,t)}dx + Q {[d ( (h h) ) 2 ]dxdt + ( λ C C 2 e λt ) Q [(f f) 2 + (g g) 2 + (h h) 2 ]dxdt 0 where C,C 2 depend on the coefficients and the solution bounds. If we choose λ > C + C 2 and t such that t < ( ) λ λ ln C, then inequality (2.9) holds if and only if f = f, g = g, w = w a.e. in Q. Similarly the proof can be completed for time intervals ([t,2t ],[2t,3t ]...). Therefore, u = u, v = v and v = v a.e. in Q, and u,v,w solve the state system (.). Hence the solution to the state system (.) exists and the uniqueness of u,v,w as solutions of (.) follows similarly as in the above argument. C 2 Theorem 2.2 There exists an optimal control in U that maximizes the functional J(u 0 ). Proof : sup{j(u 0 ) u 0 U} < since the state variables and controls are uniformly bounded. Thus there exists a maximizing sequence u 0 U such that lim n J(un 0) = sup {J(u 0 ) u 0 U}. By the existence and uniqueness of solutions to the state system (.), we define u n = u(u n 0), v n = v(u n 0), w n = w(u n 0) for each n. On a subsequence, as n, u n 0 u 0 in L2 (Q) and (u n 0 ) (u 0 ) weakly in L 2 (Q) Passing to the limit in the u n, v n, w n system and using the convergences as in Theorem 2., we have that (u,v,w) is weak solution of (.) associated with u 0. Since the payoff functional is upper semi-continuous with respect to the weak convergence, we have J(u 0 ) sup {J(un 0 ) un 0 U}. (2.9) Therefore u 0 is an optimal control that maximizes the payoff functional. 3 Derivation of the Optimality System We now derive the optimality system which consists of the state system coupled with the adjoint system. In order to obtain the necessary conditions for the optimality system

10 Crop Problem 9 we differentiate the objective functional with respect to the control. As our objective functional also depends on state variables, we differentiate the state variables with respect to the control u 0 Theorem 3. The mapping u 0 U (u,v,w) V is differentiable in the following sense : u(u 0 + l) u(u 0 ) v(u 0 + l) v(u 0 ) w(u 0 + l) w(u 0 ) ψ weakly in L 2 (Q) ψ 2 weakly in L 2 (Q) ψ 3 weakly in L 2 (0,T;H0 (Ω)) as 0 for any u 0 U and l L (Q) s.t. (u o + l) U for small. Also ψ,ψ 2,ψ 3 (depending on u,v,w,u 0,l) satisfy the following system : (ψ ) t = r (k 2u)ψ k 3 (uψ 3 + wψ ) in Q (ψ 2 ) t = r 2 (k 2 2v)ψ 2 k 4 (vψ 3 + wψ 2 ) in Q (ψ 3 ) t = d (ψ 3 ) xx + α k 3 (uψ 3 + wψ ) + α 2 k 4 (vψ 3 + wψ 2 ) µψ 3 in Q ψ (x,0) = l, ψ 2 (x,0) = l, ψ 3 (x,0) = 0 for x Ω ψ 3 (0,t) = 0 = ψ 3 (,t) on Ω (0,T) Proof : Define u = u(u 0 + l), v = v(u 0 + l), u = u(u 0 ), v = v(u 0 ) and w = w(u 0 ). We do a change of variables : u = e λt f, u = e λt f, v = e λt g, v = e λt g, w = e λt h, w = e λt h, where λ > 0 is to be chosen below. On the set Q = Ω (0,t ) for 0 < t T, we illustrate the h equation : [ (h h Q = [ Q +α 2 k 4 e λt (g ( h h ) ( h ) h t + λ( ) h h 2 ( + d h ) h 2] dxdt x α k 3 e (f λt ( ) h h 2 ( ) + h f f (h ) ) h ) 2 ( ) + h g g (h ) ) h µ ( ) h h 2 ] dxdt Continuing to estimate using L bounds on the coefficients and f,g,h,f,g and h, we have 2 Ω {t } [ (f f ) 2 + ( g g ) ] 2 ( + h ) h 2 + ( λ ( C + C 2 e λt )) Q [ (f f C Ω l2 dx. ) 2 + ( g g dx + ( Q d h ) h x 2 dxdt ) ] 2 ( + h ) h 2 dxdt

11 Crop Problem 0 For λ > C + C 2 and t small such that t < λ lnλ C, we conclude C 2 f f 2 L 2 (Q ) + g g 2 L 2 (Q ) + h h 2 C L 2 (0,t ;H0 (Ω) Ω l 2 dx. Similarly this estimate can be carried out on intervals [t,2t ],[2t,3t ], and the estimate finally holds on [0,T]. This estimate justifies the convergence of f,g and h quotients, and hence u u v v w w ψ weakly in L 2 (Q) ψ 2 weakly in L 2 (Q) ψ 3 weakly in L 2 (0,T;H0 (Ω)) Similarly we obtain and ( w ) w t (ψ 3) t weakly in L 2 (0,T;H (Ω)) w w ψ 3, strongly in L 2 (Q) u u, v v are uniformly Lipschitz in t for each x. These convergences also give u u, v v, w w strongly in L 2 (Q). To see the system satisfied by ψ,ψ 2,ψ 3, consider terms from the system satisfied by u u, v v and w w ; for example ) r ((u ) 2 u 2 = r (u u)(u + u) 2r uψ as 0, k 3 (u w uw) = k 3 (u (w w) + w (u u)) k 3 (uψ 3 + wψ ), since u u, w w as 0. The above estimates justify passing the limit in the system satisfied by u u, v v and w w, and we conclude that ψ,ψ 2,ψ 3 solves (3.). To derive the optimality system and to characterize the pairs of optimal controls, we need adjoints and adjoints of the operators associated with ψ,ψ 2,ψ 3 system as ψ 0 L ψ 2 = 0, 0 ψ 3

12 Crop Problem where L ψ ψ 2 ψ 3 = L ψ L 2 ψ 2 L 3 ψ 3 L ψ L 2 ψ 2 L 3 ψ 3 = M (ψ ) t (ψ 2 ) t ψ ψ 2 ψ 3 (ψ 3 ) t d (ψ 3 ) xx r (k 2u) k 3 w 0 k 3 u M = 0 r 2 (k 2 2v) k 4 w k 4 v α k 3 w α 2 k 4 w α k 3 u + α 2 k 4 v µ p 0 We define the adjoint PDE system as L q = 0 r 0 where p L L q = p p L 2 q Mτ q r L 3 r r with L p L 2 q L 3 r = p t q t r t d r xx. (3.) and M τ is the transpose of matrix M. Note that A,A 2 from the objective functional will occur as final time values for p,q. To clarify the characterization of our optimal control, we make the following definition involving a variational inequality with upper and lower obstacles [3] : Definition 3. u 0 U is a weak solution of the following bilateral variational inequality min{max( (p(x,0) q(x,0) + (u 0 ) xx B u 0 + B 2 (a u 0 )),u 0 a),u 0 0} = 0 if for all v 0 U, [(u 0 ) x (v 0 u 0 ) x + {q(x,0) p(x,0) + B u 0 B 2 (a u 0 )} (v 0 u 0 )] dx 0. Ω Theorem 3.2 Existence of Weak Solution Given an optimal control u 0 and corresponding solution (u,v,w) = (u(u 0 ),v(u 0 ),w(u 0 ))

13 Crop Problem 2 there exists a weak solution (p,q,r) (L 2 (Q)) 2 L 2 (0,T;H0 (Ω)) satisfying the adjoint system : and transversality conditions L p = (r (k 2u) k 3 w)p + α k 3 wr L 2 q = (r 2(k 2 2v) k 4 w)q + α 2 k 4 wr L 3 r = k 3up k 4 vq + (α k 3 u + α 2 k 4 v µ)r p(x,t) = A, q(x,t) = A 2, r(x,t) = 0 where x Ω r(0,t) = 0 = r(,t) where t [0,T]. And furthermore u 0 (x) must satisfy the following variational inequality in the weak sense. min{max( (p(x,0) q(x,0) + (u 0 ) xx B u 0 + B 2 (a u 0 )),u 0 a),u 0 0} = 0 Proof : Let u 0 (x) be an optimal control (which exists by Theorem 2.2) and (u,v,w) be its corresponding state solution. Let (u 0 (x) + l) U for > 0, and u,v,w be the corresponding weak solution of state system (.). Since the adjoint equations are linear, there exists a weak solution p,q,r satisfying ( ). We compute the directional derivative of the objective functional J(u 0 ) with respect to u 0 in the direction l at u 0. Since J(u 0 ) is the maximum value, we have 0 lim J(u 0(x)+l) J(u 0 ) 0 + { [ ( = lim u ) ( u A + v )] v A2 (x,t)dx [ (u 0 B 0 +l) 2 u 2 0 (a (u + B 0 +l)) 2 (a u 0 ) 2 2 ((u 0 +l) t) 2 (u 0 ) 2 x ] } dx = 0 [(A ψ + A 2 ψ 2 )(x,t) (u 0 ) x l x B u 0 l + B 2 (a u 0 )l] dx = 0 [(pψ + qψ 2 )(x,t) (u 0 ) x l x B u 0 l + B 2 (a u 0 )l] dx = 0 (pψ (x,0) + qψ 2 (x,0) (u 0 ) x l x B u 0 l + B 2 (a u 0 )l)dx + T 0 0 [p(ψ ) t + (r (k 2u) k 3 w)p + α k 3 wr]ψ dxdt + T 0 0 [q(ψ 2) t + (r 2 (k 2 2v) k 2 w)q + α 2 k 4 wr] ψ 2 dxdt + T 0 0 [r(ψ 3) t + d (ψ 3 ) x r x k 3 up k 4 vq + (α k 3 u + α 2 k 4 v µ)r] ψ 3 dxdt = 0 [(pψ + qψ 2 )(x,0) (u 0 ) x l x B u 0 l + B 2 (a u 0 )l] dx + Q (p,q,r) L = 0 [lp(x,0) lq(x,0) (u 0) x l x B u 0 l + B 2 (a u 0 )l] dx dx ψ ψ 2 ψ 3 dxdt

14 Crop Problem 3 The above inequality can be re-written in a formal way as (after integrating by parts the (u 0 ) x l x term) : 0 l [p(x,0) q(x,0) + (u 0 ) xx B u 0 + B 2 (a u 0 )] dx 0 On the set where u 0 = 0, we choose variation l with support on this set and l 0, which implies (p(x,0) q(x,0) + (u 0 ) xx B u 0 + B 2 (a u 0 )) 0. 0 < u 0 < a, l arb sign, (p(x,0) q(x,0) + (u 0 ) xx B u 0 + B 2 (a u 0 )) = 0 u 0 = a, l 0, (p(x,0) q(x,0) + (u 0 ) xx B u 0 + B 2 (a u 0 )) 0 This can be written in the compact form as min{max ( (p(x,0) q(x,0) + (u 0 ) xx B u 0 + B 2 (a u 0 )),u 0 a),u 0 0} = 0. Our optimality system (OS) is : L u = r u(k u) k 3 uw L 2 v = r 2 v(k 2 v) k 4 vw L 3 w = α k 3 uw + α 2 k 4 vw µw L p = (r (k 2u) k 3 w)p + α k 3 wr L 2 q = (r 2(k 2 2v) k 4 w)q + α 2 k 4 wr L 3r = k 3 up k 4 vq + (α k 3 u + α 2 k 4 v µ)r (3.2) min{max ( (p(x,0) q(x,0) + (u 0 ) xx B u 0 + B 2 (a u 0 )),u 0 a),u 0 0} = 0 u(x,0) = u 0 (x), v(x,0) = v 0 (x), w(x,0) = w 0 (x) for x Ω w(0,t) = 0 = w(,t), r(0,t) = 0 = r(,t) on (0,T) p(x,t) = A, q(x,t) = A 2, r(x,t) = 0 for x Ω r(0,t) = 0 = r(,t) where t [0,T]. The weak solution of the optimality system exists by Theorem 2.2 and 3.2. For small T, we now prove the uniqueness of weak solutions of optimality system, which gives the characterization (3.2) of the unique optimal control in terms of the solutions of the optimality system. Note that such a small T restriction is common in optimal control problems involving parabolic PDEs ; see the uniqueness results in optimal control of the PDE/ODE systems in [2, 6]. 4 Uniqueness of the Optimality System Theorem 4. For T sufficiently small and B + B 2 sufficiently large, weak solutions of the optimality system are unique.

15 Crop Problem 4 Proof : Suppose u,v,w,p,q,r,u 0 and u,v,w,p,q,u 0 are two solutions of the optimal system (3.2). We change the variables for λ > 0 to be chosen such that u = e λt u u = e λt u 2 v = e λt v v = e λt v 2 w = e λt w w = e λt w 2 p = e λt p p = e λt p 2 v = e λt q q = e λt q 2 r = e λt r r = e λt r 2. The variational inequality for u 0 becomes [ { } ] (u 0 ) x (v 0 u 0 ) x + e λt (q p )(x,0) + B u 0 B 2 (a u 0 ) (v 0 u 0 ) dx 0. Ω Substituting v 0 = u 0 in the u 0 variational inequality, and adding the resulting inequalities, we obtain Ω (u 0 u 0 ) 2 x + (B + B 2 )(u 0 u 0 ) 2 Ω e λt (u 0 u 0 )((q q 2 ) (p p 2 )) (x,0)dx. Using the weak form of the state and adjoint systems and the above inequality, gives Ω [ (p p 2 ) 2 + (q q 2 ) 2 + (r r 2 ) 2] (x,0)dx + [ Ω (u u 2 ) 2 + (v v 2 ) 2 + (w w 2 ) 2] (x,t)dx +λ [ Q (u u 2 ) 2 + (v v 2 ) 2 + (w w 2 ) 2 + (p p 2 ) 2 + (q q 2 ) 2 + (r r 2 ) 2] dxdt + ] Ω [((u 0 u 0 ) x ) 2 + (B + B 2 )(u 0 u 0 ) 2 dx + ( Q d (w w 2 ) x 2 + (r r 2 ) x 2) dxdt + Ω 2 ( C + C 2 e 2λT)[ Q { (u u 2 ) 2 + (v v 2 ) 2 + (w w 2 ) 2 +(p p 2 ) 2 + (q q 2 ) 2 + (r r 2 ) 2} dxdt ] [ (p p 2 ) 2 + (q q 2 ) 2] (x,0)dx + 2 Ω (u 0 u 0 ) 2 dx + 2 Ω (u 0 u 0 ) 2 dx, where the last term comes from the initial conditions on u, v, u, v. If we take λ large and then T small, we have λ C C 2 e 2λT > 0. If we also assume B + B 2 > 5, then we obtain the uniqueness. 2

16 Crop Problem 5 5 Numerical Realization We perform the numerical realization of the problem using the optimality system (3.2) that we recall thereafter : du = r u(k u) k 3 uw, u(0) = u o dt dv = r 2 v(k 2 v) k 4 uw, v(0) = a u o dt w 2 w = d t x 2 + (α k 3 u + α 2 k 4 v µ)w, w x=0,x= = 0, w(0) = w o dp dt dq dt r t = (r (k 2u) k 3 w) p + α k 3 w r, p(t) = A = (r 2 (k 2 2v) k 4 w) q + α 2 k 4 w r, q(t) = A 2 = d 2 r x 2 + (α k 3 u + α 2 k 4 v µ)r k 3 up k 4 v q, r x=0,x= = 0, min 0 u o a 2 0 r(t) = 0. (5.) (5.2) [ ] 2 duo + (B + B 2 )u 2 o 2 (ab 2 + p(0) q(0)) u o dx. (5.3) dx In the sequel, we set B = B + B 2 and f = (B 2 + p(0) q(0)). These equations are coupled and we are going to solve this system via a relaxation method that can be roughly described a follows :. Initialization step Choose u o 2. Iteration n : u n is known. Relaxation forward-backward method (a) Solve the forward system (5.) : we get (u n,v n,w n ). (b) Solve the backward system (5.2) with (u n,v n,w n ) : we get (p n,q n,r n ). (c) Solve the Variational Inequality (5.3) with f n = B 2 + p n (0) q n (0) : we get u n+. 3. Check a stopping criterion and set n = n + if necessary. The discretization of the system is done via finite difference methods with respect to the time variable and the space variable. The implicit Euler scheme is used to solve the forward system ODE s and the space-discretized part of PDE s and the Crank-Nicholson scheme is used for the backward (linear) system. This choice has been made using many numerical tests : though the Crank-Nicholson scheme is inconditionnally stable we could not avoid scattering for some examples that were particularly ill-conditionned. The state system (5.) is a non-linear system, we use Newton s method to solve it. The initial point

17 Crop Problem 6 is chosen as the previous iterate so that the convergence is quite fast. Note that the ratio between the time discretization step and the space discretization step has to be small enough (CFL condition) to avoid numerical unstability. To solve the Variational Inequality (5.3), we discretize the energy functional with the trapezoidal integration rule (for example) and use the classical projected gradient method to solve it. Indeed the functional is quadratic and constraints are bound constraints. The discretized problem turns out to be min 2 X MX F X, 0 X i a,i =,,N + where X = (u o (x i )) i=,n is the discretized control function, X denotes the transpose of X, F is the (space) discretized vector for f = B 2 + p(0) q(0) and M = A + B Id. Also A is the discretized D- Laplacian matrix : A = 0 2. ( x) We have performed numerical tests with the following parameters : r r 2 d k k 2 α α 2 k 3 k 4 µ a A A 2 B B 2 T We have set w o (x) = 0x( x), N = 200, K = 0 and the initial guess for u o is equal to 0.5 ; the tolerance has been set to 0 4 and the parameter of projected gradient method is ρ = 0.5. We report hereafter the results. Global iteration n # of projected gradient iterations u n o uo n.852e e e e-05 So the global iterations number is 3 and the value of the cost functional is J =

18 Crop Problem 7 Références [] Ames, W. F., Numerical Methods for Partial Differential Equations, Barnes and Noble, 969. [2] S. Chawla and S. Lenhart, Application of Optimal Control Theory to in situ Bioremediation, IMA Preprint # 49, Computational and Applied Math Journal 4(2000), [3] M. Chipot,Variational Inequalities and Flow in Porous Media, Springer-Verlag, New York, 984. [4] Evans, L.C., Partial Differential Equations, American Mathematical Society, Providence, 998. [5] Federer,W. T., Statistical Design and Analysis for Intercropping Experiments Volume I : Two Crops, Springer-Verlag, New York, 993. [6] N. Handagama and S. Lenhart, Optimal Control of a PDE/ ODE System Modeling a GasPhase Bioreactor, referred Proceedings of the Conference on Mathematical Models in Medical and Health Sciences, Vanderbilt University, 998, [7] D. Kinderlehrer and G. Stampacchia, An Introduction to Variational Inequalities and their Applications, Academic Press, New York, 980. [8] Liebman, M. and E. Dyck, Crop Rotation and Intercropping Strategies for Weed Management, Ecological Applications 3 (993), [9] Plant, R. E. and M. Mangel, Modeling and Simulation in Agricultural Pest Management, SIAM Review 29 (987), [0] Simon, J., Compact sets in the space L 2 (0,T;B), Annali di Matematica Pura et Applicata CXLVI (987),