![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
A more realistic binge drinking model with time delay is introduced. Time delay is used to represent the time lag of the immunity against drinking in our model. For the model without the time delay, using Routh–Hurwitz criterion, we obtain that the alcohol-free equilibrium is locally asymptotically stable if . We also obtain that the unique alcohol-present equilibrium is locally asymptotically stable if
. For the model with time delay, the local stability of all the equilibria is derived. Furthermore, regardless of the time delay length, using comparison theorem and iteration technique, we obtain that the alcohol-free equilibrium is globally asymptotically stable under certain conditions. Numerical simulations are also conducted to illustrate and extend our analytic results.
1. Introduction
The drinking behaviours of people pose a significant public health concern. US surveys indicate that approximately 90% of college students have consumed alcohol at least once [Citation7], and more than 40% of college students have engaged in binge drinking [Citation9,Citation20]. As we all know, excessive drinking is not only harmful to personal health, but also leads to a range of negative social effects such as violence, antisocial and criminal behaviour. Long-term alcoholism will produce negative change in the brain, such as tolerance and physical dependence. Alcohol damages almost all parts of the body and contributes to a number of human disease including live cirrhosis, pancreatitis, heart disease, sexual dysfunction and eventually death [Citation1]. The growing excessive drinking increases the risk of breast cancer by 13% every year. In addition, there are data show that, on average, heavy drinker die 10-20 years earlier than nondrinkers. 40% and 32% of family divorce are caused by alcohol, in America and China, respectively [Citation6,Citation15,Citation16]. The World Health Organization reports that the harmful use of alcohol causes approximately 3.3 million deaths every year(or 5.9% of all the global deaths), and 5.1% of the global burden of disease is attributable to alcohol consumption. Although there had been many attempts to reduce the problem, alcohol by young people had persisted and in some cases increased over the past several decades [Citation19]. Alcoholism is becoming more and more dangerous and serious as well as a widespread social phenomenon, so the study of alcoholism is an important aspect of social epidemic.
Mathematical models can mimic the process of drinking behaviour. Since the social interaction is considered to be the key factor in spreading the behaviour which can result in adverse health effects, alcoholism can be viewed as a treatable contagious disease. So far, several different mathematical models for drinking have been formulated and studied. Manthey et al. [Citation10] built an epidemiological model to capture the dynamics of campus drinking. Their results showed that the reproductive numbers were not sufficient to predict whether drinking behaviour would persist on campus and that the pattern of recruiting new members played a significant role in the reduction of campus alcohol problems. Xiang et al. [Citation23] dealt with the global property of a drinking model with public health educational campaigns, and concluded that the public health educational campaigns of drinking individuals could slow down the drinking dynamics.
It is well known that taking into account the time delay effect is a typical approach for modelling oscillatory behaviour of population. Huo and Ma [Citation4] developed a nonlinear alcoholism model with awareness programmes and time delay, they incorporated time delay during which a susceptible individual became a heavy drinker. They concluded that the time delay in alcohol consumption habit which developed in susceptible population might result in Hopf bifurcation by increasing the value of time delay. Wang et al. [Citation18] formulated an alcohol quitting model in which they considered the impact of distributed time delay between contact and infection process by characterizing dynamic nature of alcoholism behaviours, and they generalized the infection rate to the general case. Furthermore they considered two control strategies, that is, education and treatment. Other models about drinking, smoking or epidemic can be found in [Citation3,Citation5,Citation8,Citation11–13,Citation17,Citation21,Citation22,Citation24,Citation25].
Motivated by the above works, we set up a more realistic drinking model with two stages and time delay. We divide alcoholics into light problem alcoholics and heavy problem alcoholics, and incorporate time delay which describes the time lag of the ‘immunity’ against drinking in our model. Our goal is to study dynamics of our model. The remainder of this paper is organized as follows. In Section 2, we formulate our mathematical model. The existence of equilibria are obtained in Section 3. In Section 4, the stability of all the equilibria is studied. To support our theoretical results, some numerical simulations are included in Section 5. In Section 6, we give some conclusions and discussions.
2. The model formulation
Research Group of treatment and rehabilitation of the U.S. National Institute on Alcohol Abuse and Alcoholism (NIAAA), said that the daily alcohol consumption for men is not more than four standard volume (standard volume= 10 g of pure alcohol), women are not more than three standard volume. The ceiling of low-risk alcohol consumption per week is 14 standard volume for men, and 7 standard volume for women. If a person whose alcohol consumption is more than daily or weekly drinking ceiling, he/she most likely develop to abuse alcohol or addicted alcohol [Citation14]. So we divide the whole population into four compartments: ,
,
,
.
denotes susceptible individuals who do not drink or only drink moderately. The light problem alcoholics, denoted by
, refers to the drinkers who drink beyond daily or weekly ceiling a little. The heavy problem alcoholics, denoted by
, refers to the drinkers who drink far exceeding daily and weekly ceiling. The recovered individuals, denoted by
, refers to the people who quit problem drinking temporarily. The total number of population at time t is given by
. The transfer diagram is shown in Figure .
Figure leads to the following drinking model with time delay:
(1)
(1)
where Λ is the recruitment rate of susceptible individuals who do not drink or drink only moderately. We assume that if a susceptible individual is infected by alcoholics, he will enter into light problem alcoholics. β is the transmission coefficient of infection for the susceptible individuals from the light problem alcoholics. γ is the transmission coefficient of infection for the susceptible individuals from the heavy problem alcoholics. A light problem alcoholic will become a heavy problem alcoholic at rate p due to excessive drinking. A light problem alcoholic will become a recovered individual at rate m due to treatment. A heavy problem alcoholic will become a recovered individual at rate δ due to treatment. A recuperator individual becomes a susceptible at rate η. We assume that all individuals' death rate is μ.
represents that an individual has survived natural death in a recovery compartment before becoming susceptible again, where
is a constant representing the time lag of the ‘immunity’ against drinking. η is the rate of losing ‘immunity’ against drinking,
implies that the recovered individuals would lose the ‘immunity’,
implies that the recovered individuals are with the permanent ‘immunity’. The parameters are all positive.
Suppose that the initial conditions for system (Equation1(1)
(1) ) take the form
(2)
(2) where
the Banach space of continuous functions mapping the interval
into
where
It is well known that system (Equation1(1)
(1) ) has a unique solution
satisfying the initial conditions (Equation2
(2)
(2) ) by the fundamental theory of functional differential equations [Citation2]. It is easy to show that all solutions of (Equation1
(1)
(1) ) with initial conditions (Equation2
(2)
(2) ) are defined on
and remain positive for all
.
Using the fact that , system (Equation1
(1)
(1) ) is rewritten as
(3)
(3) from Equation (Equation3
(3)
(3) ), we can obtain
So,
. Hence, the feasible region of system (Equation1
(1)
(1) ) is given by the set
which is positively invariant with respect to system (Equation1
(1)
(1) ). We will consider dynamic behaviour of system (Equation1
(1)
(1) ) on the set Ω.
3. Existence of equilibria
According to the condition of system (Equation1(1)
(1) ) has a unique alcohol-present equilibrium, we can define the value of the basic reproductive number
as follows:
(4)
(4) It is easy to know the following theorem.
Theorem 3.1
System (Equation1(1)
(1) ) always has the alcohol-free equilibrium
. Furthermore, system (Equation1
(1)
(1) ) has a unique alcohol-present equilibrium
regardless of the time delay length if
.
Proof.
It is easy to see that the alcohol-free equilibrium always exists. From the third equation of Equation (Equation1
(1)
(1) ), let the right-hand side equal to zero, we have
(5)
(5)
Using the fourth equation of Equations (Equation1(1)
(1) ) and (Equation5
(5)
(5) ),
can be written as
(6)
(6)
Together the first and the second equations of Equation (Equation1(1)
(1) ), using Equations (Equation5
(5)
(5) ) and (Equation6
(6)
(6) ), we obtain
(7)
(7)
Using the second equation of Equations (Equation1(1)
(1) ) and (Equation5
(5)
(5) ),
can be written as
(8)
(8)
If , then
, so
. Hence, if
system (Equation1
(1)
(1) ) has a unique endemic equilibrium
regardless of the time delay length.
4. Stability analysis of equilibria
4.1. Stability analysis of equilibria With ![](//:0)
![](//:0)
Theorem 4.1
The alcohol-free equilibrium of system (Equation1
(1)
(1) ) with
is local asymptotically stable if
.
Proof.
The characteristic matrix of system (Equation1(1)
(1) ) is as follows:
(9)
(9) where
denotes any equilibrium of Equation (Equation1
(1)
(1) ). The associated transcendental characteristic equation of (Equation9
(9)
(9) ) is given by
(10)
(10) It is clear that the associated transcendental characteristic equation of (Equation10
(10)
(10) ) at
becomes
(11)
(11) For
, Equation (Equation11
(11)
(11) ) can be written as
(12)
(12) We can obtain the characteristic roots are
, and roots of the following equation:
(13)
(13) If
, we have
, then,
The real parts of all eigenvalues of
are negative. Thus, the alcohol-free equilibrium
of system (Equation1
(1)
(1) ) with
is local asymptotically stable.
Next, we turn to consider the local stability of the alcohol-present equilibrium with
.
Theorem 4.2
If the alcohol-present equilibrium
of the system (Equation1
(1)
(1) ) with
is local asymptotically stable.
Proof.
The characteristic equation of (Equation10(10)
(10) ) at
is
(14)
(14) where
Since
, then
. It is easy to see
.
Let , we have
By complex calculations, we also can get
,
.
Thus, by using the Routh–Hurwitz criterion, we obtain that all the eigenvalues of Equation (Equation14(14)
(14) ) have negative real parts. Hence, if
, the alcohol-present equilibrium
of the system (Equation1
(1)
(1) ) with
is local asymptotically stable.
4.2. Stability analysis of equilibria with ![](//:0)
![](//:0)
Theorem 4.3
If and
the alcohol-free equilibrium
of system (Equation1
(1)
(1) ) with
is local asymptotically stable.
Proof.
Linearizing the system (Equation1(1)
(1) ) with
around
we can obtain the characteristic roots are
and roots of the following:
(15)
(15) where
If
, it is easy to see
. Then Equation (Equation15
(15)
(15) ) does not have non-negative real solutions. Now we consider whether Equation (Equation15
(15)
(15) ) has purely imaginary solutions or not.
Suppose is a root of Equation (Equation15
(15)
(15) ). Then we have
(16)
(16) Separating the real and imaginary parts, we have the following system:
(17)
(17)
(18)
(18) To eliminate the trigonometric functions we square both sides of each equation above and we add the squared equations (Equation17
(17)
(17) ) and (Equation18
(18)
(18) ) to obtain the following equation:
(19)
(19) i.e.
(20)
(20) Let
, we obtain
(21)
(21) where
If
and
, then
, so Equation (Equation21
(21)
(21) ) does not have positive solutions and purely imaginary solutions. Hence, if
and
, the alcohol-free equilibrium
of system (Equation1
(1)
(1) ) with
is local asymptotically stable.
Next, we will study the global stability of the alcohol-free equilibrium of system (Equation1
(1)
(1) ) with
Theorem 4.4
Let ,
and the (H1), (H2) hold:
then the alcohol-free equilibrium
of system (Equation1
(1)
(1) ) is globally asymptotically stable.
Proof.
Let be any positive solution of system (Equation1
(1)
(1) ) with initial conditions (Equation2
(2)
(2) ). From
Equation (Equation3
(3)
(3) ), we have
(22)
(22)
Hence, for sufficiently small , there is a
such that if
,
For sufficiently small
, it follows from the second and the third equation of system (Equation1
(1)
(1) ), for
,
(23)
(23) Where
. Note that if
and
hold, we have
(24)
(24) Then,
. Hence, for sufficiently small
, there is a
such that if
.
For sufficiently small , by the forth equation of system (Equation1
(1)
(1) ), for
, we have
(25)
(25) Then,
(26)
(26) We derive from the first equation of system (Equation1
(1)
(1) ) that, for
,
(27)
(27) So,
(28)
(28) Together with Equation (Equation22
(22)
(22) ), yields
.
Noting that if , the alcohol-free equilibrium
of system (Equation1
(1)
(1) ) is locally asymptotically stable, we conclude that
is globally asymptotically stable. This completes the proof.
Next, we turn to study of the stability of the alcohol-present equilibrium of system (Equation1
(1)
(1) ) with
.
Theorem 4.5
If , the alcohol-present equilibrium
of the system (Equation1
(1)
(1) ) with
is local asymptotically stable.
Proof.
The characteristic equation of is
(29)
(29) where
If
, it is easy to see
. Then Equation (Equation29
(29)
(29) ) does not have non-negative real solutions. Now we consider whether Equation (Equation29
(29)
(29) ) have purely imaginary solutions or not.
Suppose is a root of Equation (Equation29
(29)
(29) ). Then we have
(30)
(30) Separating the real and imaginary parts, we have the following system:
(31)
(31)
(32)
(32) From Equations (Equation31
(31)
(31) ) and (Equation32
(32)
(32) ), we can get
(33)
(33) i.e.
(34)
(34)
where
Let
again, we obtain
(35)
(35) By complex calculations, we can get
if
Equation (Equation35
(35)
(35) ) does not have positive roots and Equation (Equation35
(35)
(35) ) can not have purely imaginary solutions. Hence, If
, the alcohol-present equilibrium
of the system (Equation1
(1)
(1) ) with
is local asymptotically stable.
5. Numerical simulation
In this section, some numerical results of system (Equation1(1)
(1) ) are presented for supporting the analytic results obtained above. We take
(which are taken from [Citation12], other parameters are estimated as follows
First, we choose p=0.9, and m=0.9, numerical simulation gives
,
, then the alcohol-free equilibrium
is globally asymptotically stable (Figure (a)).
Second, we choose p=0.1, , m=0.1, numerical simulation gives
,
, then the alcohol-present equilibrium
is globally asymptotically stable (Figure (b)).
Third, we choose p=0.9, , m=0.9, numerical simulation gives
, then the alcohol-free equilibrium
is globally asymptotically stable (Figure (a)).
At last, we choose p=0.1, , m=0.1, numerical simulation gives
,
, then the alcohol-present equilibrium
is globally asymptotically stable (Figure (b)).
6. Discussion
We have formulated a more realistic drinking model with two stages and time delay. We obtain the basic reproduction number of the model. For the model without time delay, the alcohol-free equilibrium is locally asymptotically stable if the basic reproduction number
. When
, there exists a unique alcohol-present equilibrium which is locally asymptotically stable under a certain condition. For the model with time delay, the stability of all the equilibria are obtained. Furthermore, regardless of the time delay length, using comparison theorem and iteration technique, we obtain that the alcohol-free equilibrium is globally asymptotically stable. Numerical simulations also confirm the theoretical results.
Direct calculation can show that ,
,
,
, then decreasing contact rate and increasing the treat rate are helpful to reduce the alcohol abuse. By Figure , we can see that
increases with the increase of β or γ, that means higher contact rate can increase the propagation of alcoholism. Figure shows the relation among the basic reproduction number
, the treat coefficient δ and m. From Figure , we find that it is obvious that m works more effective than δ. Therefore, it is important to treat light problem alcoholics in order to control the propagation of alcoholism.
It is easy to see the time delay can influence the dynamic behaviours of the system. From Equation (Equation4(4)
(4) ) we know that
is independent of the time delay. However, the stability of equilibria in the model depends on the time delay. It is obvious that time delay can influence the stability of the alcohol-free equilibrium according to Theorems 4.3 and 4.4.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- M.M. Glavas and J. Weinberg, Stress, alcohol consumption and the hypothalamicpituitary adrenal axis, in Nutrients, Stress and Medical Disorders, Humana Press, New York, 2006, pp. 165–183.
- J.K. Hale, Theory of Functional Differential Equations, Springer, New York, 1997.
- H.F. Huo and L.X. Feng, Global stability for an HIV/AIDS epidemic model with different latent stages and treatment, Appl. Math. Model. 37 (2013), pp. 1480–1489. doi: 10.1016/j.apm.2012.04.013
- H.F. Huo, S.H. Ma and X.Y. Meng, Modelling Alcoholism as a Contagious Disease: A Mathematical Model with Awareness Programs and Time Delay, Discret. Dyn. Nat Soc. 2015 Article ID 260195 (13 pages), 2015.
- H.F. Huo and X.M. Zhang, Complex dynamics in an alcoholism model with the impact of Twitter, Math. Biosci. 281 (2016), pp. 24–35. doi: 10.1016/j.mbs.2016.08.009
- M. Inoue, Association between alcohol consumption and colorectal cancer risk, Cur. Nutrition Rep, 2 (2013), pp. 71–73. doi: 10.1007/s13668-012-0033-z
- L.D. Johnston, P.M. OMalley and J.G. Bachman, National Survey Results on Drug use from the Monitoring the Future Study, 1975–1992, National Institute on Drug Abuse, Rockville, MD, 1993.
- S. Lee, E. Jung and C. Castillo-Chavez, Optimal control intervention strategies in low- and high-risk problem drinking populations, Socio-Econ. Plan. Sci. 44 (2010), pp. 258–265. doi: 10.1016/j.seps.2010.07.006
- P.M. O'Malley and L.D. Johnston, Epidemiology of alcohol and other drug use among American college students, J. Stud. Alcohol 63(14) (2002), pp. 23–39. doi: 10.15288/jsas.2002.s14.23
- J.L. Manthey, A.Y. Aidoob and K.Y. Ward, Campus drinking: an epidemiolodical model, J. Biol. Dyn. 2 (2008), pp. 346–356. doi: 10.1080/17513750801911169
- A. Mubayi, P.E. Greenwood, C. Castillo-Chavez, P.J. Gruenewald and D.M. Gorman, The impact of relative residence times on the distribution of heavy drinkers in highly distinct environments, Socio-Econ. Plan. Sci. 44(1) (2010), pp. 45–56. doi: 10.1016/j.seps.2009.02.002
- G. Mulone and B. Straughan, Modeling binge drinking, Int. J. Biomath. 5(1) (2012), pp. 1250005. (14 pages). doi: 10.1142/S1793524511001453
- S. Mushayabasa and C.P. Bhunu, Modelling the effects of heavy alcohol consumption on the transmission dynamics of gonorrhea, Nonlinear Dyn. 66 (2011), pp. 695–706. doi: 10.1007/s11071-011-9942-4
- National Institute on Alcohol Abuse and Alcoholism, Department of Health and Human Services, United States of America. Available at http://www.niaaa.nih.gov/alcohol-health, last accessed March 11, 2016.
- J. Rehm, The risk associated with alcohol use and alcoholism, Alcohol Res. Health 34 (2011), pp. 135–143.
- J. Ribes, R. Clèries, L. Esteban, V. Moreno, and F.X. Bosch, The influence of alcohol consumption and hepatitis B and Cinfections on the risk of liver cancer in Europe, J. Hepatol. 49(2) (2008), pp. 233–242. doi: 10.1016/j.jhep.2008.04.016
- O. Sharomi and A.B. Gumel, Curtailing smoking dynamics: a mathematical modeling approach, Appl. Math. Comput. 195 (2008), pp. 475–499. doi: 10.1016/j.amc.2007.05.012
- X.Y. Wang, K. Hattaf, H.F. Huo and H. Xiang, Stability analysis of a delayed social epidemics model with general contact rate and its optimal control, J. Indus. Manag. Optim. 12 (2016), pp. 1267–1285. doi: 10.3934/jimo.2016.12.1267
- H. Wechsler, J.E. Lee, T.F. Nelson and M.C. Kuo, Underage college students drinking behavior, access to alcohol, and the influence of deterrence policies, J. Amer. College Health 50(5) (2002), pp. 223–236. doi: 10.1080/07448480209595714
- H. Wechsler, J.E. Lee, M. Kuo and H. Lee, College binge drinking in the 1990s:a continuing problem.Results of the Harvard School of Public Health 1999 College Alcohol Study, J. Amer. College Health 48(5) (2000), pp. 199–210. doi: 10.1080/07448480009599305
- WHO, Global Status Report on Alcohol and Health 2014, WHO, Geneva, Switzerland, 2014.
- H. Xiang, L.X. Feng and H.F. Huo, Stability of the virus dynamics model with Beddington–DeAngelis functional response and delays, Appl. Math. Model. 37 (2013), pp. 5414–5423. doi: 10.1016/j.apm.2012.10.033
- H. Xiang, N.N. Song and H.F. Huo, Modelling effects of public health educational campaigns on drinking dynamics, J. Biol. Dyn. 10 (2016), pp. 164–178. doi: 10.1080/17513758.2015.1115562
- H. Xiang, Y.L. Tang and H.F. Huo, A viral model with intracellular delay and humoral immunity, Bull. Malays. Math. Sci. Soc. (2016) doi: 10.1007/s40840-016-0326-2.
- G. Zaman, Qualitative behavior of giving up smoking models, Bull. Malays. Math. Sci. Soc. 34 (2009), pp. 403–415.