Innovation diffusion model in patch environment. (Q1855916)

A mathematical model is proposed to describe the dynamics of users of one product in two different patches. Advertisement force, contact rate between users and non-users, population dispersal rate and returning rate from user class to non-user class are chosen as key parameters. For the product with a long life-span, it is assumed that the returning rate is proportional to the number of users. It is shown that this model has a unique positive equilibrium which is a globally stable. For the product with a short life-span, time delays are introduced to represent the duration of the product in two different markets and the stability of a positive equilibrium is analyzed in one reasonable case. Periodic advertisements are also incorporated and the existence and uniqueness of positive periodic solutions are investigated. The model designs in assuming, that two patches are isolated, and the dynamics of both classes are governed by \[ \begin{aligned} \frac{dN_1} {dt}&=-\gamma_1N_1- \lambda_1A_1N_1-\delta_1N_1+ \nu_1A_1+\beta_1,\\ \frac{dA_1} {dt}&=\gamma_1N_1+ \lambda_1A_1N_1-\delta_1A_1- \nu_1A_1,\end{aligned} \] where \(N_1\) is the number of non-users in the first patch, \(A_1\) is the number of users in the first patch, \(\beta_1\) is the birth rate of the population in the first patch, \(\delta_1\) is the death rate of the population in the first patch, \(\gamma_1\) represents the intensity of advertisement for promoting product-users in the first patch, and \(\lambda_1\) is the valid contact rate of users of the product with non-users in the first patch. Similarly, it is supposed that the dynamics of non-user-class and user-class in the second patch are governed by \[ \begin{aligned} \frac{dN_2 } {dt}&=-\gamma_2 N_2 - \lambda_2 A_2 N_2 -\delta_2 N_2 + \nu_2 A_2 +\beta_2 , \\\frac{dA_2 } {dt}&=\gamma_2 N_2 + \lambda_2 A_2 N_2 -\delta_2 A_2 - \nu_2 A_2 ,\end{aligned} \] where \(N_2 \) is the number of non-users in the second patch, \(A_2 \) is the number of users in the second patch, \(\beta_2 \) is the birth rate of the population in the second patch, \(\delta_2 \) is the death rate of the population in the second patch, \(\gamma_2 \) represents the intensity of advertisement for promoting product-users in the second patch, and \(\lambda_2 \) is the valid contact rate of users of the product with non-users in the second patch. But in case, when the two patches are connected, \(\theta_1\) is the probability that an individual migrates from the first patch to the second patch, \(\theta_2 \) the probability that an individual migrates from the second patch to the first patch. Further, it is assumed that \(k_2 \) is the fraction of new immigrants from the second patch who remain in user-class and that \(k_1\) is the fraction of new immigrants from the first patch who remain in user-class. Under these assumptions, the model is: \[ \begin{aligned} \frac{dN_1} {dt}&=-\gamma_1N_1- \lambda_1A_1N_1-\delta_1N_1+ \nu_1A_1+\beta_1+ \theta_2 (N_2 +(1-k_2 )A_2 )- \theta_1N_1,\\ \frac{dA_1} {dt}&=\gamma_1N_1+ \lambda_1A_1N_1-\delta_1A_1- \nu_1A_1-\theta_1A_1+ k_2 \theta_2 A_2 ,\\ \frac{dN_2 } {dt}&=-\gamma_2 N_2 - \lambda_2 A_2 N_2 -\delta_2 N_2 + \nu_2 A_2 +\beta_2 + \theta_1(N_1+(1-k_1)A_1)- \theta_2 N_2 ,\\ \frac{dA_2 } {dt}&=\gamma_2 N_2 + \lambda_2 A_2 N_2 -\delta_1A_2 - \nu_2 A_2 -\theta_2 N_2 + k_1\theta_1A_1.\end{aligned} \] In order to incorporate duration of the product into the model, described above was modified. The assumption that the returning rate is proportional to the number of users seems reasonable when the length of the life-span of the product is large. Otherwise, the effect of the length will be significant. Because of this, time delays due to the life-span of the product was introduced into the model and it was obtained: \[ \begin{multlined}\frac{dN_1} {dt}=-\gamma_1N_1(t)- \lambda_1A_1(t) N_1(t) -\delta_1N_1(t) +e^{-\delta_1 \tau_1} (\gamma_1N_1(t-\tau_1)\\+\lambda_1A_1(t-\tau_1) N_1(t-\tau_1)) +\beta_1+ \theta_2 (N_2 +(1-k_2 )A_2 )- \theta_1N_1,\end{multlined} \] \[ \begin{multlined}\frac{dA_1} {dt}=\gamma_1N_1(t) +\lambda_1A_1(t) N_1(t) -\delta_1A_1(t) - e^{-\delta_1 \tau_1} (\gamma_1N_1(t-\tau_1)\\ +\lambda_1A_1(t-\tau_1) N_1(t-\tau_1))-\theta_1A_1(t) +k_2 \theta_2 A_2 (t) ,\end{multlined} \] \[ \begin{multlined} \frac{dN_2 } {dt}=-\gamma_2 N_2 (t) - \lambda_2 A_2 (t) N_2 (t) -\delta_2 N_2 (t) + e^{-\delta_2 \tau_2 } (\gamma_2 N_2 (t-\tau_2 )\\+\lambda_2 A_2 (t-\tau_2 ) N_2 (t-\tau_2 ))+\beta_2 +\theta_1(N_1+(1-k_1)A_1)- \theta_2 N_2 ,\end{multlined} \] \[ \begin{multlined}\frac{dA_2 } {dt}=\gamma_2 N_2 (t) + \lambda_2 A_2 (t) N_2 (t) -\delta_2 A_2 (t) - e^{-\delta_2 \tau_2 } (\gamma_2 N_2 (t-\tau_2 )\\ +\lambda_2 A_2 (t-\tau_2 ) N_2 (t-\tau_2 ))-\theta_2 A_2 (t) +k_1\theta_1A_1(t).\end{multlined} \] Here \(\tau_1\) represents the life-span of the product in the first patch and \(\tau_2 \) represents the life-span of the product in the second patch.