On drug resistance and metronomic chemotherapy: A mathematical modeling and optimal control approach

1.   Introduction

Synchronization in complex networks has been a focus of interest for researchers from different disciplines[1,2,4,8,15]. In this paper, we investigate synchronous phenomena in an ensemble of Kuramoto-like oscillators which is regarded as a model for power grid. In [9], a mathematical model for power grid is given by

where $P_{source}^{i}>0$ is power source (the energy feeding rate) of the $i$th node, $I\ge 0$ is the moment of inertia, $K_{D}>0$ is the ratio between the dissipated energy by the turbine and the square of the angular velocity, and $a_{il}>0$ is the maximum transmitted power between the $i$th and $l$th nodes. The phase angle $\theta_{i}$ of $i$th node is given by $\theta_i(t) = \Omega t+\tilde\theta_i(t)$ with a standard frequency $\Omega$ (50 or 60 Hz) and small deviations $\tilde\theta_i(t)$. Power plants in a big connected network should be synchronized to the same frequency. If loads are too strong and unevenly distributed or if some major fault or a lightening occurs, an oscillator (power plant) may lose synchronization. In that situation the synchronization landscape may change drastically and a blackout may occur. The transient stability of power grids can be regarded as a synchronization problem for nonstationary generator rotor angles aiming to restore synchronism subject to local excitations. Therefore, the region of attraction of synchronized states is a central problem for the transient stability.

By denoting $\omega_{i} = \frac{P_{source}^{i}}{K_{D}}$ and assuming $\frac{a_{il}}{K_{D}} = \frac{K}{N} \,(\forall i,l\in\{1,2,\dots,N\})$ and $I = 0$ in (1), Choi et. al derived a Kuramoto-like model for the power grid as follows [5]:

Here, the setting $\dot{\theta}_{i}>0$ was made in accordance to the observation in system (1) that the grid is operating at a frequency close to the standard frequency with small deviations. In order to ensure that the right-hand side of (2) is positive, it is reasonable to assume that $\omega_{i}>K,\; i = 1,2,\dots,N$. This model can be used to understand the emergence of synchronization in networks of oscillators. As far as the authors know, there are few analytical results. In [5], the authors considered this model with identical natural frequencies and prove that complete phase synchronization occurs if the initial phases are distributed inside an arc with geodesic length less than $\pi/2$.

If $(\dot{\theta}_{i})^2$ in (2) is replaced by $\dot{\theta}_{i}$, the model is the famous Kuramoto model [13]; for its synchronization analysis we refer to [3,6,10], etc. Sakaguchi and Kuramoto [16] proposed a variant of the Kuramoto model to describe richer dynamical phenomena by introducing a phase shift (frustration) in the coupling function, i.e., by replacing $\sin(\theta_{l}-\theta_{i})$ with $\sin(\theta_{l}-\theta_{i}+\alpha)$. They inferred the need of $\alpha$ from the empirical fact that a pair of oscillators coupled strongly begin to oscillate with a common frequency deviating from the simple average of their natural frequencies. Some physicists also realized via many experiments that the emergence of the phase shift term requires larger coupling strength $K$ and longer relaxation time to exhibit mutual synchronization compared to the zero phase shift case. This is why we call the phase shift a frustration. The effect of frustration has been intensively studied, for example, [11,12,14]. In power grid model using Kuramoto oscillators, people use the phase shift to depict the energy loss due to the transfer conductance [7]. In this paper, we will incorporate the phase shift term in (2) and study the Kuramoto-like model with frustration $\alpha\in(-\frac{\pi}{4},\frac{\pi}{4})$:

We will find a trapping region such that any nonstationary state located in this region will evolve to a synchronous state.

The contributions of this paper are twofold: First, for identical oscillators without frustration, we show that the initial phase configurations located in the half circle will converge to complete phase and frequency synchronization. This extends the analytical results in [5] in which the initial phase configuration for synchronization needs to be confined in a quarter of circle. Second, we consider the nonidentical oscillators with frustration and present a framework leading to the boundness of the phase diameter and complete frequency synchronization. To the best of our knowledge, this is the first result for the synchronization of (3) with nonidentical oscillators and frustration.

The rest of this paper is organized as follows. In Section 2, we recall the definitions for synchronization and summarize our main results. In Section 3, we give synchronization analysis and prove the main results. Finally, Section 4 is devoted to a concluding summary.

Notations. We use the following simplified notations throughout this paper:

2.   Preliminaries

In this paper, we consider the system

Next we introduce the concepts of complete synchronization and conclude this introductory section with the main result of this paper.

Definition 2.1. Let $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a dynamical solution to the system (4). We say

1. it exhibits asymptotically complete phase synchronization if

2. it exhibits asymptotically complete frequency synchronization if

For identical oscillators without frustration, we have the following result.

Theorem 2.2. Let $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a solution to the coupled system (4) with $\alpha = 0,\omega_{i} = \omega_0,i = 1,2,\dots,N$ and $\omega_0>K$. If the initial configuration

then there exits $\lambda_1,\lambda_2>0$ and $t_{0}>0$ such that

and

Next we introduce the main result for nonidentical oscillators with frustration. For $\bar{\omega}< \frac{1}{2\sin^2|\alpha|}$, we set

For suitable parameters, we denote by $D_{1}^{\infty}$ and $D_{*}^{\infty}$ the two angles as follows:

Theorem 2.3. Let $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a solution to the coupled system (4) with $K>K_{c}$ and $[\underline{\omega},\bar{\omega}]\subset \left[K+1,\frac{1}{2\sin^2|\alpha|}\right)$. If

then for any small $\varepsilon>0$ with $D_{1}^{\infty}+\varepsilon<\frac{\pi}{2}$, there exists $\lambda_{3}>0$ and $T>0$ such that

Remark 1. If the parametric conditions in Theorem 2.3 are fulfilled, the reference angles $D_{1}^{\infty}$ and $D_{*}^{\infty}$ are well-defined. Indeed, because $K>K_{c}$ and $\bar{\omega}<\frac{1}{2\sin^2 |\alpha|}$, we have

This implies

Then, by $\underline{\omega}\ge K+1$ and $K\le \bar{\omega}$ we obtain

Remark 2. In order to make $1<K+1<\frac{1}{2\sin^2|\alpha|}$, it is necessary to assume $\alpha\in \left(-\frac{\pi}{4},\frac{\pi}{4}\right)$. This is the reason for the setting $\alpha\in \left(-\frac{\pi}{4},\frac{\pi}{4}\right)$.

3.   Synchronization analysis

3.1.   Synchronization estimates: Identical oscillators without frustration

In this subsection we consider the system (4) with identical natural frequencies and zero frustration:

To obtain the complete synchronization, we need to derive a trapping region. We start with two elementary estimates for the transient frequencies.

Lemma 3.1. Suppose $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a solution to the coupled system (8), then for any $i,j\in\{1,2,\dots,N\},$ we have

Proof. It is immediately obtained by (8).

Lemma 3.2. Suppose $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a solution to the coupled system (8) and $\Omega>K$, then we have

Proof. It follows from (8) and $\dot{\theta}_{i}>0$ that we have

Next we give an estimate for trapping region and prove Theorem 2.2. For this aim, we will use the time derivative of $D(\theta(t))$ and $D(\nu(t))$. Note that $D(\theta(t))$ is Lipschitz continuous and differentiable except at times of collision between the extremal phases and their neighboring phases. Therefore, for the collision time $t$ at which $D(\theta(t))$ is not differentiable, we can use the so-called Dini derivative to replace the classic derivative. In this manner, we can proceed the analysis with differential inequality for $D(\theta(t))$. For $D(\nu(t))$ we can proceed in the similar way. In the following context, we will always use the notation of classic derivative for $D(\theta(t))$ and $D(\nu(t))$.

Lemma 3.3. Let $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a solution to the coupled system (8) and $\Omega>K$. If the initial configuration $\theta(0)\in \mathcal{A}$, then $\theta(t)\in \mathcal{A},\; t\ge 0$.

Proof. For any $\theta^0\in \mathcal A$, there exists $D^\infty\in (0,\pi)$ such that $D(\theta^0)<D^{\infty}$. Let

Since $D(\theta^0)<D^{\infty}$, and $D(\theta(t))$ is a continuous function of $t$, there exists $\eta>0$ such that

Therefore, the set $\mathcal T$ is not empty. Let $T^*: = \sup \mathcal T.$ We claim that

Suppose to the contrary that $T^*<\infty$. Then from the continuity of $D(\theta(t))$, we have

We use Lemma 3.1 and Lemma 3.2 to obtain

Here we used the relations

and

Therefore, we have

which implies that

This is contradictory to $D(\theta(T^*)) = D^{\infty}.$ The claim (9) is proved, which yields the desired result.

Now we can give a proof for Theorem 2.2.

Proof of Theorem 2.2.. According to Lemma 3.3, we substitute $T^* = \infty$ into (10), then (5) is proved with $\lambda_1 = \frac{K\cos\frac{D^{\infty}}{2}}{\pi\sqrt{\omega_0+K}}$.

On the other hand, by (5) there exist $t_{0}$ and $\delta (0<\delta<\frac{\pi}{2})$ such that $D(\theta(t))\le \delta$ for $t\ge t_{0}$. Now we differentiate (8) to find

Using Lemma 3.2, we now consider the temporal evolution of $D(\nu(t))$:

This implies that

and proves (6) with $\lambda_2 = \frac{K\cos\delta}{2\sqrt{\omega_0+K}}$.

Remark 3. Theorem 2.2 shows, as long as the initial phases are confined inside an arc with geodesic length strictly less than $\pi$, complete phase synchronization occurs exponentially fast. This extends the main result in [5] where the initial phases for synchronization need to be confined in an arc with geodesic length less than $\pi/2$.

3.2.   Synchronization estimates: Nonidentical oscillators with frustration

In this subsection, we prove the main result for nonidentical oscillators with frustration.

Lemma 3.4. Let $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a solution to the coupled system (4), then for any $i,j\in\{1,2,\dots,N\},$ we have

Proof. By (4) and for any $i,j\in\{1,2,\dots,N\}$

the result is immediately obtained.

Lemma 3.5. Let $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a solution to the coupled system (4) and $\underline{\omega}\ge K+1$, then we have

Proof. From (4), we have

and also because $\dot{\theta}_{i}>0$. The following lemma gives a trapping region for nonidentical oscillators.

Lemma 3.6. Let $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a solution to the coupled system (4) with $K>K_{c}$ and $[\underline{\omega},\bar{\omega}]\subset [K+1,\frac{1}{2\sin^2|\alpha|})$. If the initial configuration $\theta^0\in \mathcal{B}$, then $\theta(t)\in \mathcal{B}$ for $t\ge 0.$

Proof. We define the set $\mathcal T$ and its supremum:

Since $\theta^0\in \mathcal{B}$, and note that $D(\theta(t))$ is a continuous function of $t$, we see that the set $\mathcal T$ is not empty and $T^*$ is well-defined. We now claim that

Suppose to the contrary that $T^*<\infty$. Then from the continuity of $D(\theta(t))$, we have

We use Lemma 3.4 to obtain

For $I$, we have

We now consider two cases according to the sign of $\alpha$.

(1) $\alpha\in[0,\frac{\pi}{4}).$ In this case, we have

(2) $\alpha\in(-\frac{\pi}{4},0).$ In this case, we have

Here we used the relations

and

Since $D(\theta(t))+|\alpha|<D_{*}^{\infty}<\pi$ for $t\in[0,T^*),$ we obtain

By (12) and Lemma 3.5 we have

Then we obtain

i.e.,

Here we used the definition of $\sin D_{*}^{\infty}$. By Gronwall's inequality, we obtain

Thus

Let $t\to T^{*-}$ and we have

which is contradictory to $D(\theta(T^*)) = D_{*}^{\infty}-|\alpha|.$ Therefore, we have

That is,

Lemma 3.7. Let $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a solution to the coupled system (4) with $K>K_{c}$ and $[\underline{\omega},\bar{\omega}]\subset [K+1,\frac{1}{2\sin^2|\alpha|})$. If the initial configuration $\theta(0)\in \mathcal{B}$, then

Proof. It follows from (11) and Lemma 3.5, Lemma 3.6 and that we have

The proof is completed.

Lemma 3.8. Let $\theta: = (\theta_{1},\theta_{2},\dots,\theta_{N})$ be a solution to the coupled system (4) with $K>K_{c}$ and $[\underline{\omega},\bar{\omega}]\subset [K+1,\frac{1}{2\sin^2|\alpha|})$. If the initial configuration $\theta(0)\in \mathcal{B}$, then for any small $\varepsilon>0$ with $D_{1}^{\infty}+\varepsilon<\frac{\pi}{2}$, there exists some time $T>0$ such that

Proof. Consider the ordinary differential equation:

It is easy to find that $y_{*} = D_{1}^{\infty}$ is a locally stable equilibrium of (13), while $y_{**} = D_{*}^{\infty}$ is unstable. Therefore, for any initial data $y^0$ with $0<y^0<D_{*}^{\infty}$, the trajectory $y(t)$ monotonically approaches $y_{*}$. Then, for any $\varepsilon>0$ with $D_{1}^{\infty}+\varepsilon<\frac{\pi}{2}$, there exists some time $T>0$ such that

In particular, $y(t) <y_*+\varepsilon$ for $t\ge T.$ By Lemma 3.7 and the comparison principle, we have

which is the desired result.

Remark 4. Since

we have $D_{1}^{\infty}>|\alpha|$.

Proof of Theorem 2.3. It follows from Lemma 3.8 that for any small $\varepsilon>0$, there exists some time $T>0$ such that

We differentiate the equation (4) to find

We now consider the temporal evolution of $D(\nu(t))$:

where we used

Thus we obtain

and proves (7) with $\lambda_{3} = \frac{K\cos(D_{1}^{\infty}+\varepsilon)}{2\sqrt{\bar{\omega}+K}}$.

4.   Conclusions

In this paper, we presented synchronization estimates for the Kuramoto-like model. We show that for identical oscillators with zero frustration, complete phase synchronization occurs exponentially fast if the initial phases are confined inside an arc with geodesic length strictly less than $\pi$. For nonidentical oscillators with frustration, we present a framework to guarantee the emergence of frequency synchronization.

Acknowledgments

We would like to thank the anonymous referee for his/her comments which helped us to improve this paper.