A Theoretical Study of Biological Lotka-Volterra Ecological Model Using Comprehensive Thermodynamic Theory of Stability of Irreversible Processes (Cttsip)
Total Page:16
File Type:pdf, Size:1020Kb
International Journal of Knowledge Engineering ISSN: 0976-5816 & E-ISSN: 0976-5824, Volume 3, Issue 1, 2012, pp.-91-94. Available online at http://www.bioinfo.in/contents.php?id=40 A THEORETICAL STUDY OF BIOLOGICAL LOTKA-VOLTERRA ECOLOGICAL MODEL USING COMPREHENSIVE THERMODYNAMIC THEORY OF STABILITY OF IRREVERSIBLE PROCESSES (CTTSIP) RAWAT S.G.1, BHALEKAR A.A.2 AND TANGDE V.M.3* 1Department of Applied Chemistry, Priyadarshini College of Engineering, Nagpur, MS, India. 2Department of Chemistry, Rashtrasant Tukadoji Maharaj Nagpur University, Nagpur, MS, India. 3Department of Applied Chemistry, Smt. Bhagwati Chaturvedi College of Engineering, Nagpur, MS, India. *Corresponding Author: Email- [email protected] Received: February 28, 2012; Accepted: March 06, 2012 Abstract- The dynamical relationship between predator and prey (Lotka-Volterra model or host-parasitoid system)[1] is one of the dominant themes in ecology. It was observed from the population data that interaction between a pair of predator-prey influences the population growth of both the species. This paper presents the study of thermodynamic stability of periodic Lotka-Volterra system against Prey popula- tion perturbation. The thermodynamic stability of representative model of Lotka-Volterra ecosystem has been investigated using proposed thermodynamic Lyapunov function in CTTSIP [2, 3, 4, 5] which follows the steps of Lyapunov's second method (also termed as direct meth- od) of stability of motion[6, 7]. The thermodynamic Lyapunov function, used herein is the excess rate of entropy production in thermodynam- ic perturbation space that conforms well with the dictates of second law of thermodynamics[8, 9]. Moreover, present study reveals the re- gions of stability, asymptotic stability and instability. Keywords- Lotka-Volterra model, perturbation, predator, prey, Lyapunov function, irreversible thermodynamics Citation: Rawat S.G., Bhalekar A.A. and Tangde V.M. (2012) A Theoretical Study of Biological Lotka-Volterra Ecological Model Using Com- prehensive Thermodynamic Theory of Stability of Irreversible Processes (CTTSIP). International Journal of Knowledge Engineering, ISSN: 0976-5816 & E-ISSN: 0976-5824, Volume 3, Issue 1, pp.-91-94. Copyright: Copyright©2012 Rawat S.G., et al. This is an open-access article distributed under the terms of the Creative Commons Attribu- tion License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are cred- ited. Introduction an appropriate thermodynamic space using the entropy source Lyapunov's second method of stability of motion [6, 7, 10, 11, 12] strength function, the latter is a positive definite quantity as per has a generality element, which is very much akin to that of ther- the second law of thermodynamics that appears in Clausius- modynamics. In this direct method there is no need to rigorously Duhem inequality[8, 9, 15], that reads as, solve the differential equation of motion. Indeed, this course is ds adopted in Lyapunov's first method [6, 10, 11, 12], wherein one rs+divJ ss = ³ 0 tackles every case individually. This then allows the investigator to dt (1) map the regions of stability and instability (if any). In essence, Lyapunov's second method involves effecting of a sufficiently r J s small perturbations and then observing the response in the corre- where, is the mass density, s is per unit mass entropy, sponding perturbation space. Previously, the stability aspects of s s some elementary chemical reactions including autocatalytic ones is the entropy flux density, is the entropy source strength [13] and industrial sulphur dioxide oxidation [14] has already been and t is the time. investigated using this generalized framework of CTTSIP[5]. Here 0 in the present case we will be using “Comprehensive thermody- Yti () namic theory of stability of irreversible processes(CTTSIP)” which Let 's be the thermodynamic coordinates of the given uses the fabric of Lyapunov's direct method of stability of motion. Yti () CTTSIP[2, 3, 4, 5] defines a thermodynamic Lyapunov function in process on the unperturbed trajectory and 's be that on International Journal of Knowledge Engineering ISSN: 0976-5816 & E-ISSN: 0976-5824, Volume 3, Issue 1, 2012 Bioinfo Publications 91 A Theoretical Study of Biological Lotka-Volterra Ecological Model Using Comprehensive Thermodynamic Theory of Stability of Irreversible Processes (CTTSIP) i ()t i the perturbed one. The perturbation coordinates are de- punov’s theory of stability of motion20-24 ’s need to be fined as, sufficiently small the higher order terms in the expansion of 0 (t ) ( Y ( t ) Y ( t )); ( i 1,2, , n ) LtS () i i i (2) can be ignored and hence the expressions for it get and hence the equation of unperturbed trajectory is obtained as, restricted to, 0 0. 00 i 0 SS L (3) LS S S > i i iiii ()t (10) Thus i 's determine the thermodynamic perturbation space. 0 The constitutive equations of motion within the perturbation space LSi/ in general are described as , In Eq. (10) each is either positive or negative d L i f( , , , , t ) ( i 1,2, n ). The total time derivative of S then reads as, dt i 111 (4) dL L L d S S S i LS The thermodynamic Lyapunov function, is defined as, dt ti i dt (11) 0 LSSS( ( t ) ( t )) 0 0 (5) Si/0 On the other hand, if and hence 0 S S 0 Where and are the corresponding entropy LSi/0 source strengths. , and as stated above in view of the Obviously, Eq. (26) leads to the following two possibilities, either requirement of Lyapunov’s theory of stability of motion i ’s L L( , , , , t ) 0 need to be sufficiently small the third and higher order terms in the SS1 2 3 ( 6 ) Lt() or expansion of S can be ignored and hence the expressions LSS L(1 , 2 , 3 , , t ) 0 L ( 7 ) for S , get reduced to, 0 2 and on the unperturbed trajectory we have, 1 L () 0 S L L(0,0,0, , t ) 0 S S S i j SS 2 ij, ( 8 ) ij 0 L 2 S 1 L where, , is a differentiable function. In the generalized setup S ˆ 0 ij (12) 2 ij, ij of CTTSIP we expand the entropy source strength, S , and L thermodynamic Lyapunov function, S , in a Taylor expansion in L The total time derivative of S then reads as, i 00 terms of the perturbation coordinates 's as, dL1 22 d S S S i i j i 0 0 dt2 t dt 1 2 ii j j i i j L () 0 sS (13) S S S i i j ii2 i, j i j where we have used the following expression for the local time 0 3 1 LS S derivative of , namely: 6 i j k i,, j k i j k 0 (9 2 LSS1 Thus, even the bilateral perturbations can also be treated by this ij tt2 0 ij, ij (14) Si/0 formalism. However, if and hence Thermodynamic stability is then determined by finding out the sign 0 dL/ dt LSi/0 S , and as per the requirement of Lya- of on the lines of Lyapunov’s direct method of stability International Journal of Knowledge Engineering ISSN: 0976-5816 & E-ISSN: 0976-5824, Volume 3, Issue 1, 2012 Bioinfo Publications 92 Rawat S.G., Bhalekar A.A. and Tangde V.M. of motion[6,7,11,12]. Thus the unperturbed trajectory is said to be dx=- x x stable if , Perturbation in Prey population ( 0 ) The Thermodynamic stability of above model is then analyzed dLS LS .0 using CTTSIP under the perturbation effect in Prey population ( dt (15) The asymptotic stability is obtained if the result is either, ). Let us consider a case when the concentra- dL x L 0, S 0, tion of prey is perturbed by a suffciently small amount say, S dt (16) dx . We then have by the definition of perturbation coordinate, dL S dx=- x x 0 ˜ 0 LS 0, 0, (20) Or dt (17) Now, we investigate the effect of such perturbation on the thermo- dynamic stability of process below. Let us consider the population where, , is a strictly positive number that van- rate equation of prey as ishes only at the origin. The stability under constantly acting small disturbances, as per Malkin’s theorem [6], is obtained if in addition dx =-ax bxy dt LSi/ (21) to Eq. (16) or (17) each is finite. Following the track of CTTSIP let us directly move for deriving a suitable thermodynamic Lyapunov function for the present case, System under investigation as In this study, we investigate the thermodynamic stability of Lotka- Volterra ecosystem model (Prey-Predator interaction). A L()=-a by0 d x In ecosystem, one species do not occur in isolation of another. s T They compete for resources, assists, excludes and kills one an- (22) other. In ecology, predation describes a biological interaction and on solving the steps through-out we finally reach to the total where a predator feeds on its prey. The predator-prey relationship L is important in maintaining the balance between the different spe- time derivative of s as cies. Without prey, there would be no predators and vice-versa. e. d (L ) A g. Cheetahs and gazelles, foxes and rabbits, lion and foxes, foxes s =-2 (a by02 ) d x and goat, birds and insects,goat and grass, etc. The Lotka- dt T Volterra Prey-Predator model involves coupled equations, one (23) which describes how the prey population changes and second Conclusion which describes how the predator population changes, namely: Analytical conclusions regarding the thermodynamic stability of above system against the perturbation in prey population have put forward in following different situations, namely: dx =-ax bxy dt 1.