arxiv: v1 [math.ap] 8 Nov 2014

Transcription

1 arxiv: v1 [math.ap] 8 Nov 014 ALL INVARIANT REGIONS AND GLOBAL SOLUTIONS FOR m-component REACTION-DIFFUSION SYSTEMS WITH A TRIDIAGONAL SYMMETRIC TOEPLITZ MATRIX OF DIFFUSION COEFFICIENTS SALEM ABDELMALEK Abstract. The purpose of this paper is the construction of invariant regions in which we estabish the goba existence of soutions for m-component reaction-diffusion systems with a tridiagona symmetric toepitz matrix of diffusion coefficients and with nonhomogeneous boundary conditions. The proposed technique is based on invariant regions and Lyapunov functiona methods. The noninear reaction term has been supposed to be of poynomia growth. 1. Introduction In recent years, the goba existence for soutions for noninear paraboic systems of partia differentia equations have received considerabe attention. On of the most promising wors that can be found in the itarature is that of Jeff Morgan [0], where a the components satisfy the same boundary conditions (Neumann or Dirichet), and the reaction terms are poynomiay bounded and satisfy m inequaities. Hois ater competed the wor of Morgan and estabished goba existence in the presence of mixed boundary conditions subject to certain structure requirements of the system. In 007, Abdemae and Kouachi [] show that soutions of the m-component reaction diffusion systems with a diagona diffusion matrix exist gobay (for m ). The resuts obtained in this wor represent the proof of the goba existence of soutions with Neumann, Dirichet, nonhomogeneous Robin and a mixture of Dirichet with nonhomogeneous Robin conditions. The reaction terms are again assumed to be of poynomia growth and satisfy a singe inequaity. The diffusion matrix is a simpe symmetric tridiagona one. A aong the paper, we wi use the foowing notations and assumptions: we denote by m the number of equations of the system (i.e. an m-component system): u 1 t a u 1 b u = f 1 (U), u t b u 1 a u b u +1 = f (U); =,,, u m b u a u m = f m (U), t in {t > 0}, (1.1) 000 Mathematics Subject Cassification. Primary 35K45, 35K57. Key words and phrases. Reaction-Diffusion Systems, Invariant Regions, Goba Existence. 1

2 SALEM ABDELMALEK with the boundary conditions: and the initia data: where: αu +(1 α) η u = β, = 1,,m, on {t > 0}, (1.) u (0,x) = u 0 (x), = 1,,m, on, (1.3) (i) for nonhomogeneous Robin boundary conditions, we use 0 < α < 1, β R, = 1,,m, or (ii) for homogeneous Neumann boundary conditions, we use α = β = 0, = 1,,m, or (iii) for homogeneous Dirichet boundary conditions, we use 1 α = β = 0, = 1,,m. Here is an open bounded domain of cass C 1 in R N with boundary, η denotes the outward norma derivative on, and U = (u ) m. The constants a and b are supposed to be positive non nu and satisfying the condition: π bcos < a. (1.4) The initia data are assumed to be in the regions: { { m } (u 0 Σ L,Z = 1,,u 0 m) R m such that =1 u0 ( )π sin 0, L m =1 u0 zπ sin 0,z Z, (1.5) with { m =1 β sin ( )π 0, L m =1 β sin ( )π 0,z Z, where { L Z = L Z = {1,,,m}. Hence, we can see that there are m regions. The subsequent wor is simiar for a of these regions as wi be shown at the end of the paper. Let us now examine the first region and then comment on the remaining cases. The chosen region is the case where L = {1,,,m} and Z = : we have { } (u 0 Σ L, = 1 m),,u0 R m such that u 0 ( )π sin 0, L, with =1 =1 β sin ( )π 0, L. (1.6) The aim is now to study the goba existence of sutions for the reaction-diffusion system in (1.1) in this region. In order to achieve this aim, we need to diagonaize the diffusion matrix, see formua (4.1). First, et us define the reaction diffusion functions as: F (w 1,w,,w m ) = =1 f (U)sin ( )π, (1.7)

3 ALL INVARIANT REGIONS AND GLOBAL SOLUTIONS 3 where the variabe w is given by w = u sin ( )π. (1.8) =1 The defined function must satisfy the foowing three conditions: (A1) The functions F are continuousy differentiabe on R m + for a = 1,,m, satisfying F (w 1,,w 1,0,w +1,,w m ) 0, for a u 0; = 1,,m. (A) The functions F are of poynomia growth (see Hois and Morgan [9]), which means that for a = 1,,m with integer N 1: F (W) C 1 (1+ ) N w on (0,+ ) m. (1.9) (A3) The foowing inequaity: ) D F (W)+F m (W) C (1+ w ; (1.10) hods for a w 0; = 1,,m and a constants D D ; = 1,,m where D ; = 1,,m are positive constants sufficienty arge. Note that C 1 andc are positive and uniformy bounded functions defined on R m +.. Preiminary Observations and Notations The usua norms in spaces L p (), L () and C() are denoted respectivey by: u p p = 1 u(x) p dx; u = esssup u(x), (.1) x u C() = max u(x). x It is we-nown that to prove the goba existence of soutions to a reactiondiffusion system (see Henry [7]), it suffices to derive a uniform estimate of the associated reaction term on [0;T max ) in the space L p () for some p > n/. Our aim is to construct poynomia Lyapunov functionas aowing us to obtain L p bounds on the components, which eads to goba existence. Since the reaction terms are continuousy differentiabe on R m +, it foows that for any initia data in C(), it is easy to chec directy their Lipschitz continuity on bounded subsets of the domain of a fractiona power of the operator λ λ 0 O = (.) 0 0 λ m Under these assumptions, the foowing oca existence resut is we nown (see Friedman [4] and Pazy []). Remar 1. Assumption (A1) contains smoothness and quasipositivity conditions that guarantee oca existence and nonnegativity of soutions as ong as they exist, via the maximum principe (see Smoer [4]). Assumption (A3) is the usua

4 4 SALEM ABDELMALEK poynomia growth condition necessary to obtain uniform bounds from p dependent L P estimates. (see Abdemae and Kouachi [], and Hois and Morgan [11]). 3. Some Properties of diffusion matrix Lemma 1. Considering the proposed reaction-diffusion system in (1.1), the resuting diffusion matrix can be given by: a b b a b A = 0 b a..... (3.1) b b a b b a This matrix is said to be positive definite if the condition in (1.4) is satisfied. Proof. The proof of this emma can be found in [13]. Note that if the matrix is positive definite, it foows that deta > 0. Lemma. The eigenvaues (λ < λ 1 ; =,,m) of A are positive and are given by ( ) π λ = a+bcos, (3.) with the corresponding eigenvectors being v = ( t, sin π π mπ,sin,,sin ) for = 1,,m. Hence, we concude that A is diagonaizabe. In the remainder of this wor we require an ascending order of the eigenvaues. In order to simpify the indices in the formuas to come we define ( ( )π λ = λ = a+bcos thus ( λ < λ +1 ; =,,m ). ) ; =,,m, (3.3) Proof. Reca that the diffusion matrix is positive definite, hence it the eigenvaues arenecessariypositive. For an eigenpair(λ,x), the components in (A λi)x = 0 are bx 1 +(a λ)x +bx +1 = 0, = 1,,m, with x 0 = x = 0, or equivaenty, ( ) a λ x + + x +1 +x = 0, = 0,,. b These are second-order homogeneous difference equations, and soving them is simiar to soving anaogous differentia equations. The technique is to see soutions of the form x = ξr for constants ξ and r. This produces the quadratic equation ( ) a λ r + r+1 = 0, b withrootsr 1 andr.itcanbearguedthatthegenerasoutionofx + + ( ) a λ x+1 + x = 0 is x = { αr 1 +βr, if r 1 r, αρ +βρ, if r 1 = r = ρ, b

5 ALL INVARIANT REGIONS AND GLOBAL SOLUTIONS 5 where α and β are arbitrary constants. For the eigenvaue probem at hand, r 1 and r must be distinct -otherwise x = αρ + βρ, and x 0 = x = 0 impies that each x = 0, which is impossibe because X is an eigenvector. Hence, x = αr1 +βr, and x 0 = x = 0 yieds { 0 = α+β 0 = αr 1 +βr ( r1 r ) = β α = 1 r 1 r = e iπ, therefore, r 1 = r e iπ for some 1 m. This couped with: ( ) { a λ r r1 r + r+1 = (r r 1 )(r r ) = 1 b r 1 +r = ( ) a λ, b eads to r 1 = e iπ, r = e iπ, and ( ) λ = a+b e iπ +e iπ The eigenvaues of A can, therefore, be given by ( ) π λ = a+bcos ; ( ) π = a+bcos. for = 1,,m. Since these λ s are a distinct (cosθ is a stricty decreasing function of θ on (0,π), and b 0), A is necessariy diagonaizabe. Finay, the th component of any eigenvector associated with λ satisfies x = αr 1 +βr with α+β = 0, thus x = α ( ) e iπ e iπ ( ) = iαsin π. Setting α = 1 i yieds a particuar eignvector associated with λ given by ( ( ) ( ) ( )) t 1π π mπ v = sin,sin,,sin. Because the λ s are distinct, {v 1,v,,v m }, is a compeat ineary independent set, so (v 1 v v m ) diagonaizes A. Now, et us prove that λ < λ 1 ; =,,m. We have > 1. Dividing by () and mutipying by π, we have π > ( 1)π, The cosine function cosθ is stricty decreasing in θ on (0,π), thus we have ( ) ( ) π ( 1)π cos < cos. Finay, mutipying both sides of the inequaity by b and adding a yieds ( ) ( ) π ( 1)π λ = a+bcos < a+bcos = λ 1.

6 6 SALEM ABDELMALEK 4. Main Resuts Proposition 1. The eigenvectors of the diffusion matrix associated with the eigenvaues λ are defined as v = (v 1,v,,v m ) t. Mutipying each of the m equations in (1.1) by the corresponding eement of the th eigenvector and adding the equations together yieds: w t λ w = F (w 1,w,,w m ), (4.1) and αw +(1 α) η w = ρ on {t > 0}. (4.) The reaction term F, the components w, and the ascending order eigenvaues λ have been defined perviousy in this paper. Note that condition (1.4) guarantees the paraboicity of the proposed reactiondiffusion system in (1.1)-(1.3), which impies that the system described by (4.1)- (4.) is equivaent to it in the region: Σ L, = { (u 0 1,,u 0 m) R m such that with { ρ 0 = =1 { w 0 = =1 β sin ( )π u 0 sin ( )π 0, L. This impies that the components w are necessariy positive. 0, L, Proposition. The system (4.1)-(4.) admits a unique cassica soution (w 1 ;w ;..,w m ) on (0,T max ). If T max < then im w (t,.) =, (4.3) tրt max where T max ( w 0 1, w 0,, w 0 m ) denotes the eventua bow-up time. The main resut of the paper reads as foows. Theorem 1. Suppose that the functions F ; = 1,,m are of poynomia growth and satisfy condition (1.10) for some positive constants D ; = 1,,m sufficienty arge. Let (w 1 (t,.),w (t,.),,w m (t,.)) be a soution of (4.1)-(4.) and L(t) = H pm (w 1 (t,x),w (t,x),,w m (t,x))dx, (4.4) where H pm (w 1,,w m ) = p m p =0 p Cp p m with p m a positive integer and C p Aso suppose that the foowing condition is satisfied where K p κ = K r = K r 1 r 1 Kr 1 Cp p1 θ p 1 1 θp () () wp1 1 wp p1 wm pm p, p κ! p!(p κ p )!. > 0; =,,m, (4.5) [ H r 1 ] ;r = 3,,,

7 and ALL INVARIANT REGIONS AND GLOBAL SOLUTIONS 7 ( ) H r = det (a,κ ),r+1 1,κ κ 1,..r K = λ 1 1 λ =1 θ(p +1) = θ(p +) }{{} positive vaue =r (r ) =1 (det[]) r = 3,, 1, 1 = θ(p +) +(p +1) = θ(p +) [ ] 1 =1 θ A 1, H = λ 1 λ λ θ (p1+1) 1 [ θ1 A ] A 1 A 1. }{{} positive vaue ( ) det (a,κ ),r+1 is denoted determinant of r square symmetric matrix obtained from (a,κ ) 1,κ m 1,κ κ 1,..r by removing the (r+1) th,(r+) th,, th rows and the r th,(r+1) th,,( 1) th coumns. where det[1],,det[m] are the minors of the matrix (a,κ ) 1,κ m. The eements of the matrix are: a κ = λ + λ κ θ p 1 1 θp ( 1) ( 1) θ(p +1) θ (p (κ 1)+1) κ 1 θ (pκ+) κ θ (p ()+) (). (4.6) where λ in (3.)- (3.3). Note that A κ = λ + λ κ λ λκ for a,κ = 1,,m. and θ ; = 1,,() are positive constants. It foows from these conditions that the functiona L is uniformy bounded on the interva [0,T ],T < T max. Coroary 1. Under the assumptions of theorem 1, a soutions of (4.1)-(4.) with positive initia data in L () are in L (0,T ;L p ()) for some p 1. Proposition 3. Under the assumptions of theorem 1 and given that the condition (1.4) is satisfied, a soutions of (4.1)-(4.) with positive initia data in L () are goba for some p > Nn. 5. Proofs For the proof of theorem 1, we first need to define some preparatory Lemmas. Lemma 3. With H pm being the homogeneous poynomia defined by (4.4), differentiating in w 1 yieds p w1 H pm = p m p =0 p.. C p p Cp1 p θ (p1+1) 1 θ (p ()+1) () 1 wp p1 w p3 p 3 w m (p) p. (5.1) Simiary for =,,, we have p w H pm = p m p =0 p C p p Cp1 p θ p 1 1 θp ( 1) 1 θ (p +1) θ (p ()+1) () 1 wp p1 w p3 p 3 w (p) p m. (5.)

8 8 SALEM ABDELMALEK Finay, differentiating in w m yieds p wm H pm = p m p =0 p C p p Cp p 3 Cp p1 θ p 1 1 θp θp () () 1 wp p1 w p3 p 3 w m (p) p. (5.3) Lemma 4. The second partia derivative of H pm in w 1 is given by Simiary, we obtain for a =,,, w 1 H n = p m (p m 1) θ (p1+) p m p =0 1 θ (p ()+) () w H n = p m (p m 1) w w κ H n = p m (p m 1) p 3 p p =0 C p p m Cp1 p 1 wp p1 w (pm ) p m. (5.4) p m p =0 p C p p m Cp1 θ p 1 1 θp θp 1 1 θ(p +) θ (p ()+) () p 1 wp p1 w (pm ) p m. (5.5) θ p 1 1 θp 1 1 θ(p +1) p m p =0 p C p p m Cp1 p θ (pκ 1+1) κ 1 θ κ (pκ+) θ (p ()+) () 1 wp p1 w (pm ) p m (5.6) for a 1 < κ m. Finay, the second derivative in w m is given by w m H n = p m (p m 1) p m p =0 p C p p m Cp1 p θ p 1 1 θp () () 1 wp p1 w (pm ) p m. (5.7) Lemma 5 (see Abdemae and Kouachi []). Let A be the m-square symetric matrix defined by A = (a κ ) 1,κ m then we the foowing property arises: Km m = det[m] =m (det[]) (m ), m > =1 (5.8) K = det[] where K m = K 1 1 K 1 H m = det 1,κ m K m = a 11 a mm (a 1m ), H m = a 11 a m a 1 a 1m. m ( Hm 1 ) ; = 3,,m, ( ) (a,κ ) m,+1 κ,.. = ( ) =1 (det[]) = 3,,,

9 ALL INVARIANT REGIONS AND GLOBAL SOLUTIONS 9 Proof of Theorem 1. TheaimofthisworistoprovethatL(t)isuniformybounded on the interva [0,T ],T < T max. Let us start by differentiating L with respect to t: L (t) = t H pm dx and and where = = = w w H pm t dx w H pm (λ w +F )dx λ w H pm w dx+ = I +J, I = J = w H pm F dx λ w H pm w dx, (5.9) w H pm F dx. (5.10) Using Green s formua, we can divide I into two parts I 1 and I where I 1 = λ w H pm η w dx, (5.11) [( (λ ) ) ] +λ κ I = wκw H pm T Tdx (5.1) 1,κ m forp 1 = 0,,p,p = 0,,p 3 p = 0,,p m andt = ( w 1, w,, w m ) t. Appying emmas 3 and 4 yieds ( λ + λ κ wκw H pm )1,κ m = p m (p m 1) pm ( p C p p m Cp1 p (a κ ) 1,κ m ) 1 w(pm ) p m, p =0 (5.13) where (a κ ) 1,κ m is the matrix defined in formua (4.6). Now the proof of positivity for I simpifies to proving that there exists a positive constant C 4 independent of t [0,T max ) such that and that I 1 C 4 for a t [0,T max ), (5.14) I 0, (5.15) for severa boundary conditions. First, et us prove the formua in (5.14): (i) If = 1,m : 0 < α < 1, then using the boundary conditions (1.) we get I 1 = λ w H pm (γ σ w )dx,

10 10 SALEM ABDELMALEK whereσ = α andγ = β,for = 1,m. SinceH (W) = m λ w H pm (γ σ w ) = 1 α 1 α P n 1 (W) Q n (W), where P n 1 and Q n are poynomias with positive coefficients and respective degrees n 1 and n, and since the soution is positive, it foows that im sup w + H (W) =, (5.16) which proves that H is uniformy bounded on R m + and consequenty (5.14). (ii) If for a = 1,m : α = 0, then I 1 = 0 on [0,T max ). (iii) The case of homogeneous Dirichet conditions is trivia, since in this case the positivity of the soution on [0,T max ) impies η w 0, = 1,m on [0,T max ). Consequenty, one gets obtains the same resut in (5.14) with C 4 = 0. Hence, the proof of (5.14) is compete. Now, we move to the proof of (5.15): Reca the matrix (a κ ) 1,κ m which was defined in formua (4.6). The quadraticforms(withrespectto w, = 1,,m)associatedwiththematrix(a κ ) 1,κ m, with p 1 = 0,,p, p = 0,..,p 3 p = 0,,p m, is positive definite since its minors det[1], det[], det[m] are a positive. Let us examine these minors and prove their positivity by induction: The first minor det[1] = λ 1 θ (p1+) 1 θ (p+)..θ (p ()+) () > 0 is trivia for p 1 = 0,,p, p = 0,,p 3 p = 0,,p m. For the second minor det[], according to emma 5, we get: det[] = K = λ 1λ θ (p1+1) 1 = θ(p +) [ θ 1 A ] 1. Using (4.5) for = we get det[] > 0. Simiary, for the third minor det[3], and again using emma 5, we have: K 3 3 = det[3]det[1]. Since det1 > 0, we concude that sign(k 3 3) = sign(det[3]). Again, using (4.5) for = 3 yieds det[3] > 0. To concude the proof, et us suppose det[] > 0 for = 1,,, 1 and show that det[] is necessariy positive. We have det[] > 0, = 1,,( 1) = ( ) =1 (det[]) > 0. (5.17) From emma 5, we obtain K = ( ) det[]= =1 (det[]), and from (5.17), we get: sign(k ) = sign(det[]). Since K > 0 according to (4.5) then det[] > 0 and the proof of (5.15) is concuded. It foows from (5.14) and (5.15) that I is bounded. Now, et us divert our attention to proving that J in(5.10) is bounded. Substituting the expressions of the partia derivatives given by 3 in the second integra of (5.10)

11 ALL INVARIANT REGIONS AND GLOBAL SOLUTIONS 11 yieds J = ( p m p p =0 θ(p +1) = p F 1 + p m p p =0 θ(p +1) θp = ( p m p θ (p +1) θ p C p p Cp1 κ 1 =1 θp =κ θ(p +1) κ= p F 1 + p =0 κ= C p p Cp1 κ 1 =1 θp =κ θ(p +1) p F 1 + κ= C p p 1 wp p1 wm p p F κ + θp F m ) dx p 1 wp p1 wm p p θp p Cp1 =κ θ (p +1) θ p Hence, using condition (1.10), we deduce that p p J C 5 Cp p1 C p p =0 F κ +F m θp p 1 wp p1 wm p p ) F κ +F m θp p wp1 1 wp p1 wm p p dx. ( 1+ dx ) w dx. To prove that the functiona L is uniformy bounded on the interva [0,T ], et us first write ( ) p p Cp p1 C p 1+ w p =0 = R pm (W)+S p(w), p wp1 1 wp p1 wm p p where R pm (W) and S p(w)are two homogeneous poynomias of degrees p m and p m 1, respectivey. Since a of the poynomias H pm and R pm are of degree p m, there exists a positive constant C 6 such that R pm (W)dx C 6 H pm (W)dx. (5.18) Appying Höder s inequaity to the integra S p (W)dx, one obtains ( S p(w)dx (meas) 1 pm (S p(w)) pm Since for a w 1,w,,w 0 and w m > 0, (S p(w)) pm p H pm (W) )p pm p dx. = (S p (x 1,x,,x,1)) pm p, H pm (x 1,x,,x,1)

12 1 SALEM ABDELMALEK where {1,,,} : x = u u +1 and im x + (S p(x 1,x,,x,1)) pm p H pm (x 1,x,,x,1) one asserts that there exists a positive constant C 7 such that (S p(w)) pm p H pm (W) < +, C 7, for a w 1,w,,w m 0. (5.19) Hence, the functiona L satisfies the differentia inequaity which for Z = L 1 pm can be written as L (t) C 6 L(t)+C 8 L p pm (t), p m Z C 6 Z +C 8. (5.0) A simpe integration gives the uniform bound of the functiona L on the interva [0,T ]. This ends the proof of the theorem. Proof of Coroary 1. The proof of this coroary is an immediate consequence of 1 and the inequaity ( m p w (t,x)) dx C 9 L(t) on [0,T ]. (5.1) for some p 1. Proof of Proposition 3. From coroary 1, there exists a positive constant C 10 such that ( m p w (t,x)+1) dx C 10 on [0,T max ). (5.) From (1.9), we have {1,,,m} : ( m p F (W) p N C 11 (W) W (t,x)) on [0,T max ). (5.3) Since w 1,w,,w m are in L (0,T ;L p ()) and p N > n, then as discussed in the preiminary observations section the soution is goba. 6. Fina Remars Reca that the eigenvectors of the diffusion matrix associated with the eigenvaue λ is defined as v = (v 1,v,,v m ) t. It is important to note that if v is an eigenvectorthen so is ( 1)v. In the region considered in previous sections, we ony used the positive v. The remainder of the m regions can be formed using negative versions of the eigenvectors. In each region, the reaction-diffusion system with a diagonaized diffusion matrix is formed by mutipying each of the m equations in (1.1) by the corresponding eement of either v or ( 1)v and then adding the m equations together. The equations mutipied by eements of v form a set L,

13 ALL INVARIANT REGIONS AND GLOBAL SOLUTIONS 13 whereas the equations mutipied by eements of ( 1)v form a set Z. Hence, we can write define the region in the form: { { (u 0 Σ L,Z = 1,u 0,,u 0 m) R m w such that 0 = m =1 u0 v 0, L wz 0 = ( 1) m =1 u0 v z 0,z Z,, with { ρ 0 = m =1 β v ( ) 0, L ρ 0 = ( 1) m =1 β v ( z) 0,z Z. Using emma yieds { { (u 0 Σ L,Z = 1,u0 m),,u0 R m such that with { and w 0 = m =1 u0 ( )π sin w 0 z = ( 1) m =1 u0 ρ 0 = m =1 β sin ( )π 0, L ρ 0 z = ( 1) m =1 β sin ( z)π 0,z Z, { L Z = L Z = {1,,,m}. References 0, L ( z)π sin 0,z Z, [1] S. Abdemae, Invariant Regions and Goba Existence of Soutions for Reaction-Diffusion Systems with a Tridiagona Matrix of Diffusion Coefficients and Nonhomogeneous Boundary Conditions, Journa of Appied Mathematics, Vo )007), [] S. Abdemae and S. Kouachi, A simpe proof of Syvester s (determinants) identity, Appied Mathematica Sciences, (008) vo., ,. [3] S. Abdemae and S. Kouachi, Proof of existence of goba soutions for m-component reaction-diffusion systems with mixed boundary conditions via the Lyapunov functiona method, Journa of Physics A, (007) vo. 40, [4] A. Friedman, Partia Differentia Equations of Paraboic Type. Prentice Ha Engewood Chiffs. N. J [5] W. B. Fitzgibbon, S. Hois and J. Morgan, Stabiity and Lyapunov functions for reactiondiffusion systems., SIAM J. Math. Ana. (1997) Vo. 8, [6] A. Haraux and A. Youana, On a Resut of K. Masuda Concerning Reaction-Diffusion Equations. Tôhou. Math. J. 40 (1988), [7] D. Henry, Geometric Theory of Semi-inear Paraboic Equations. Lecture Notes in Mathematics 840, Springer-Verag, New-Yor, [8] S. L. Hois, On the Question of Goba Existence for Reaction-Diffusion Systems with Mixed Boundary Conditions. Quartery of Appied Mathematics LI, number, June 1993, pp [9] S. L. Hois and J. J. Morgany, On the Bow-up of Soutions to Some Semiinear and Quasiinear Reaction-diffusion Systems, Rocy Moutain J. Math. vo 14. (1994),, [10] S. L. Hois, R. H. Martin and M. Pierre, Goba Existence and Boundedness in Reaction Diffusion Systems. SIAM. J. Math. Ana, Vo. 18, [11] S. L. Hois and J. J. Morgany, Party dissipative reaction-di usion systems and a mode of phosphorus diffusion in siicon, Noninear Ana. T. M. A., 19 (199), [1] S.L. Hois and J. Morgan, Interior estimates for a cass of reaction-diffusion systems from L 1 a priori estimates, JDE, (199), vo 9, [13] C. R. Johnson, M. Neumann & M. J. Tsatsomeros, Conditions for the positivity of determinants, Linear and Mutiinear Agebra, (1996), Vo 40, [14] M. Kirane and S. Kouachi, Goba soutions to a system of strongy couped reaction-diffusion equations. Noninear Anaysis Theory, Methods and Appications. Vo 6, number 8 (1996). [15] S. Kouachi, Existence of goba soutions to reaction-diffusion systems via a Lyapunov functiona. Eectron. J. Diff. Eqns Vo. 001(001), 1-10.

14 14 SALEM ABDELMALEK [16] S. Kouachi, Existence of goba soutions to reaction-diffusion systems with nonhomogeneous boundary conditions via a Lyapunov functiona. Eectron. J. Diff. Eqns Vo. 00(00), No. 88, pp [17] S. Kouachi and A. Youana, Goba existence for a cass of reaction-diffusion systems. Buetin of the Poish Academy of Sciences, Vo. 49, Number 3, (001). [18] S. Maham and J. Xin, Goba soutions to a reactive Boussinesq system with front data on an in nite domain. Comm. Math. Phys. 193 (1998), no., [19] R. H. Martin and M. Pierre, Noninear reaction-diffusion systemsn, Noninear equations in the appied sciences, Math. Sci. Engrg., 185, Academic Press, Boston, MA, , 199. [0] J. Morgan, Goba Existence for Semiinear Paraboic Systems, SIAM J. Math. Ana. 0, (1989). [1] J. J. Morgan and S. L. Hois, The Existence of Periodic Soutions to Reaction-Diffusion Systems with Periodic Data, SIAM J. Math. Ana. (1995) Vo. 6, [] A. Pazy, Semigroups of Linear Operators and Appications to Partia Differentia Equations, Appied Math. Sciences 44, Springer-Verag, New Yor (1983). [3] F. Rothe, Goba Soutions of Reaction-Diffusion Systems, Lecture Notes in Math. 107, Springer-Verag, Berin (1984). [4] J. Smoer, Shoc Waves and Reaction-Diffusion Equations, Springer-Verag, New Yor (1983). S. ABDELMALEK: Department of Mathematics, Coege of Sciences, Yanbu, Taibah University, Saudi Arabia. E-mai address: sam@gmai.com