Delay-induced blow-up in a planar oscillation model

In this paper we study a system of delay differential equations from the viewpoint of a finite time blow-up of the solution. We prove that the system admits blow-up solutions, no matter how small the length of the delay is. In the non-delay system every solution approaches to a stable unit circle in the plane, thus time delay induces blow-up of solutions, which we call “delay-induced blow-up” phenomenon. Furthermore, it is shown that the system has a family of infinitely many periodic solutions, while the non-delay system has only one stable limit cycle. The system studied in this paper is an example that arbitrary small delay can be responsible for a drastic change of the dynamics. We show numerical examples to illustrate our theoretical results.


Introduction
In various disciplines of the science, mathematical modelling offers a description of phenomena. In some phenomena, the history of the state, not only the current state, affects the change of the state, thus it is reasonable to consider the effect of time delay [7]. Up to now the theory of delay differential equations have been intensively developed [9,17].
In this paper we study a system of delay differential equations from the viewpoint of the blow-up solutions. Here we use the terminology "blow-up" as a finite time blow-up of the solutions, i.e., the solution diverges (in a suitable topology) in finite time. The blow-up phenomenon has been widely investigated in partial differential equations, see [2,11,18,19] and references therein. There are also extensive studies in Volterra integral equation, as an alternative formulation of partial differential equation of parabolic type with a point source term [3,18,21]. Compactification of the phase space is a method to study the blow-up solutions of ordinary differential equations [6,14]. Numerical analysis has been an unavoidable tool for understanding the blow-up phenomenon. Numerical method to compute the blow-up solutions for polynomial systems of ordinary differential equations has been proposed for ordinary differential equations with application to partial differential equations, see [11] and references therein. Recently, numerical validation for the existence of the blow-up solutions is proposed [20].
The blow-up phenomenon in delay differential equations has not been much studied, as far as we know, except for a few studies [1,8,10,13,22,23]. Perhaps the reason is that many examples of delay differential equations, which appear in population biology, control theory, etc, negative feedback condition is usually imposed which excludes blow-up solutions. One also sees that the following delay differential equation does not have a blow-up solution (at least if the initial function is continuous) for > 0 and is reduced to a famous example of the ordinary differential equation having a blow-up solution when = 0 . Thus one may speculate that time delay inhibits the blow-up phenomenon in general (which is certainly not).
In the paper [8] the authors study the existence of the blow-up solutions for a class of delay differential equations. The authors are interested if adding (and multiplying) a delay term to an ordinary differential equation affects the solution behavior. Using the comparison principle, the authors obtain conditions that the delay term does not change the qualitative properties concerning the global existence and blow-up of the solution. In [23] the authors study the blow-up phenomenon of differential equations with piecewise constant arguments, in comparison with the corresponding ordinary differential equations. The blow-up phenomenon is studied in Volterra integro-differential equations. See [1,10,13,22] as applications to a parabolic type partial differential equations to study of Volterra integro-differential equations.
It seems that the blow-up phenomenon that stems from the time delay has not been reported, to the best of our knowledge. Our motivations are to demonstrate whether the time delay itself induces blow-up of the solution, and to understand the mechanism of blow-up of solutions. In this paper, we propose an example model that time delay drastically changes the solution behavior and induces the blow-up solution together with a family of infinitely many periodic solutions, where most of solutions are shown to be unstable. In our equation, a blow-up solution exists no matter how small the length of the delay is, while non-delay equation does not have a blow-up solution. We thus call these phenomena "delay-induced blow-up".

3
Delay-induced blow-up in a planar oscillation model This paper is organized as follows. In Sect. 2, we introduce a planar system of delay differential equations which we study in this paper. In the absence of time delay, the system becomes a planar system of ordinary differential equations. It can be seen that every solution except the trivial solution approaches to the limit cycle, thus no solution blows up. Concerning the system of delay differential equations, we present our two main theorems that show blow-up of solutions is possible due to time delay and that the system admits infinitely many periodic solutions. In this section we use a special initial function in order to show blow-up of solutions for any > 0 . We demonstrate numerical examples of blow-up solutions and global solutions. See Remark 3.1, Fig. 1 and also Sect. 6. In Sect. 3, by a careful estimation of the solution we show that there is a blow-up solution for the system of delay differential equations and provide a proof of Theorem 2.1. In Sect. 4, we study existence of a periodic solution with constant radius and constant angular velocity in the plane. It is shown that there exist infinitely many periodic solutions which appear due to time delay. In Sect. 5, we study a characteristic equation which characterizes stability of the periodic solutions. It is shown that most of periodic solutions are not stable except for the only one periodic solution that is a continuation of the periodic solution of the non-delay model. In Sect. 6, we provide numerical examples which illustrate our theoretical results. In Sect. 7, we discuss our results.

A planar system of delay differential equations and main results
Let ≥ 0 be a parameter for time delay. In this paper we consider the following planar system of delay differential equations For the special case = 0 , the system (2.1) is reduced to the following system of ordinary differential equations with a = 1 : Here a ∈ ℝ is a parameter. This system is a famous model for Hopf bifurcation at a = 0 . (See [12], for instance.) From the elementary calculation, one can see that every solution of (2.2) except the trivial solution (x(t), y(t)) ≡ (0, 0) tends to a periodic solution of minimal period 2 and satisfies i.e., the limit cycle of (2.2) is the unit circle.
For the system (2.1) we prove that the system admits blow-up solutions due to the presence of time delay. We prove the following theorem in Sect. 3.

Theorem 2.1
For any > 0 , there exist blow-up solutions for the planar system of delay differential Eq. (2.1).
Then in Sects. 4 and 5, we further investigate the system (2.1) and show the existence of infinitely many periodic solutions. We also study the stability of the periodic solutions. The following theorem is proved in Sects. 4 and 5.

Theorem 2.2
For any > 0 , there exist infinitely many unstable periodic solutions of (2.1) with constant radius and angular velocity. Moreover, there exists a positive * such that the system (2.1) admits only one asymptotically stable periodic solution with constant radius and angular velocity for 0 < < * .

Delay-induced blow-up
We consider a polar coordinate system. Let By the change of the variables, from (2.1), we obtain the following polar coordinate system Now we prove that there exists a solution such that r blows up in a finite time with → 4 as r → ∞ . Therefore, we can conclude that x and y blow up in a finite time.

3
Delay-induced blow-up in a planar oscillation model For (3.1), we consider the following initial condition where R > 0 and s(t) is a continuous function such that  1 2 , from the initial condition (3.2), we have Then, for t ∈ 0, 1 2 , the solution of the system of delay differential equations (3.1) with the initial condition (3.2) is given by the following system of ordinary differential equations with the initial condition We are going to prove that the solution of the system of ordinary differential equations (3.3) with the initial condition (3.4) blows up in t ≤ 1 2 . Figure 2 shows each step for blow-up of solutions. The process for blow-up of solutions is divided into the following 3 steps: Step 1 The angle of the solutions is monotonically increasing from − ∕2 to − ∕4 in t ∈ [0, T 1 ] for some T 1 > 0 . In this region, the radius of the solutions may decrease, and thus we establish a decay estimate of the solutions in Lemma 3.1.
Step 2 Both radius and angle of the solutions monotonically increase in t ∈ [T 1 , for some T 2 > T 1 . The angle varies from − ∕4 to 0 and the radius grows up beyond a threshold for an emergence of a nullcline of the angle. (See Lemmas 3.2, 3.3.) Step 3 The final stage for blow-up of solutions. The angle monotonically reaches to the nullcline of the angle which appears in Step 2. Then we show in Lemma 3.4 that for large R there exists a blow-up time For the exposition, we define In the proof below we need to estimate the solution r. For the estimation we use the following equation which is obtained from (3.6).

Step 1: The decay estimate of the radius
The solution of the system of ordinary differential equation (3.6) with the initial condition (3.7) exists for sufficiently small t. We show that there exists T 1 > 0 such that the solution exists for t ∈ 0, T 1 satisfying that monotonically increases from ∕4 to ∕2 for 0 ≤ t ≤ T 1 . We first obtain an a priori estimate for r for ∈ 1 4 , 1 2 . From (3.8) it follows that One also obtains the following estimation Integrating the inequalities (3.10) and (3.11), we obtain the following estimation Using the a priori bound (3.12), from the equation (3.6b), we see that is an increasing function and that � (t) ≥ 1 , provided ∈ 1 4 , 1 2 . Therefore, there exists T 1 > 0 such that monotonically increases from ∕4 to ∕2 for t ∈ 0, T 1 . The inequality (3.9) holds from the estimation (3.12). Since, from the inequality (3.9), we have R exp(− ∕4) ≤ r for t ∈ 0, T 1 , the following estimation

Step 2: Emergence of the nullcline of the angle
Next we have the following estimation.
One also has that r monotonically increases for t ∈ T 1 , T 2 and Moreover, there exists R 2 > 0 such that Proof First, we derive an a priori estimate for r, provided ∈ 1 2 , 3 4 . Let ∈ 1 2 , 3 4 . We have r � (t) ≥ r from (3.6a), thus r monotonically increases. Then one has Since, from (3.8), it holds that . Hence, there exists T 2 > T 1 such that monotonically increases for t ∈ T 1 , T 2 and (T 2 ) = 3 4 . Therefore, the lower estimate (3.15) is valid for t ∈ [T 1 , T 2 ] and by virtue of (3.17) the boundedness of r(t) also holds for t ∈ [T 1 , T 2 ] . That is, the inequality (3.14) holds for t ∈ [T 1 , T 2 ]. Finally, from the inequality (3.14), the following estimation holds: We are ready to show that the solution of (3.6) blows up. Note that the solution exists and that and r increase and thus Fig. 3, we plot the graph of , ′ . We show that, for sufficiently large R, a nullcline for exists in , 3 2 , where the right-hand side of the -equation (3.6b) becomes 0. We see in Lemma 3.4 that has a upper bound for suitable large R and then r(t) blows up in a finite time. We now let R be a sufficiently large number such that R 2 exp(− 4 ) > √ 2 . Then, for r > R exp(− 4 ) , there is a -nullcline in , 5 4 that is given as It is easy to obtain the following elementary lemma.

3
Hence, as long as the solution exists, for t > T 2 , r(t) ≥ r(T 2 ) ≥ R exp − 4 increases and increases and tends to * (r(t)) . Observe that stays in the interval 3 4 , 5 4 , thus the sign of cos in the right hand side of (3.6a) is fixed and is positive. From Lemma 3.2, we have

Step 3: Blow-up of solutions
Finally we show the blow-up of solutions.

Lemma 3.4
There exists T 3 (> T 2 ) such that and r monotonically increase for t ∈ T 2 , T 3 and Furthermore, there exists R 3 > 0 such that Delay-induced blow-up in a planar oscillation model Proof We consider the system (3.6) for t > T 2 . Let R be a sufficiently large number such that R > max R 1 , R 2 and that R 2 exp(− 4 ) > from Lemmas 3.1 and 3.2, we have (T 2 ) = 3 4 and T 2 < 3 8 . We derive an a priori estimate for , provided r < ∞ . Since we have R 2 exp(− 4 ) > √ 2 , one sees that the equation (3.6b) has an equilibrium, * (r) ∈ ( , 5 4 ) from Lemma 3.3. Suppose that there exists t * > T 2 such that � (t * ) = 0 while r(t * ) < ∞ . Note that (t * ) ∈ ( , 5 4 ) and by (3.8) we have dr d | t=t * = ∞ , that is, the solution orbit crosses the curve of -nullcline {( , r)|1 + Rr sin = 0, < < 5 4 } vertically. Then, the solution enters the region where � (t) < 0 , in what follows (t) < 5 4 . Note that the intersection of the solution orbit and -nullcline may occur at most once because of the shape of the nullcline and the fact of vertical crossing. See Fig. 4.

Existence of periodic solutions
We consider a periodic solution with a constant radius, i.e. r(t) ≡ > 0 , for the system (3.1). Then, � (t) is also a constant from Eq. (3.1). Thus, we treat a periodic solution of the form where ∈ ℝ is an angular velocity. From (3.1), the periodic solution satisfies Remark that for the special case = 0 , from (4.2), we obtain for periodic solution (r(t), (t)) = (1, t) corresponding to the system (3.1) with = 0. From (4.2) one has cos > 0 . Then we have the following equations In Fig. 5a we plot the functions y = − 1 and y = tan . For ∈ (− 2 + 2n , 2 + 2n )(n ∈ ℤ), intersections of the two functions correspond to the roots of (4.4a) satisfying cos > 0 . First we study the roots of (4.4a) for 0 ≤ < 2 . The implicit function (4.4a) for 0 ≤ < 2 defines a function which attains a unique maximum at = * where � ( * ) = 0 . One can compute that from which we can numerically compute * ≈ 2.22913 ⋯ . Let * = ( * ) . Then we numerically obtain
For < * there is at least 2 periodic solutions where (r(t), (t)) = j , j t , j ∈ {0, 1}, For > 0 , one can see that the Eq. (4.2) has infinitely many roots. It is elementary to prove the following result, thus we omit the proof. See also Fig. 5a. Therefore, we obtain roots j for j ∈ ℤ if ≤ * and for j ∈ ℤ ⧵ {0, 1} if > * of the Eq. (4.4a). In Fig. 5b we plot the branches for as a function of . Once of the periodic solution (4.1) is given, is determined from (4.4b). The radius can be determined as j given as, similar to (4.6), This implies that the periodic solutions of the form (4.1) has larger radius than the unit circle which is the trajectory of the periodic solution for = 0 . We also note that there is a root = 1 if = 2n , n ∈ ℕ . In Fig. 6 we plot the branches for the radius as a function of .
Summarizing the above findings, we obtain the following result for the existence of the periodic solution (4.1) for the system (3.1). The result is not intuitive and not expected that the delay induces many periodic solutions from the non-delay system (2.1).

Theorem 4.1
For each > 0 , the Eq. (4.2) has infinitely many roots, which are countable. Thus the system (3.1) has infinitely many periodic solutions, which are countable, of the form (4.1).

Analysis of the characteristic equation for the stability of the periodic solutions
To analyze stability of the periodic solution obtained in Section 3, we study a system of a delay differential equation and an integral equation, employing the principle of linearized stability for the coupled systems of renewal equations and delay differential equations [4].

3
Delay-induced blow-up in a planar oscillation model Let v(t) = � (t) . From the system (3.1), we get the following system of a delay differential equation and a renewal equation The system (5.1) has equilibria corresponding to the periodic solutions of the system (3.1), given as in Sect. 3.
In the Appendix A, by linearization of the system (5.1) at the equilibrium (5.2), we obtain the following characteristic equation with = j , j ∈ ℕ . Note that we have a family of the characteristic equations (5.3) which are indexed by the equilibrium about which we linearize.
Define ∶= . Then, from (5.3), we obtain the following equation where f ∶ ℂ → ℂ defined as We study (5.4) to analyze the distribution of the roots of (5.3) with respect to the imaginary axis. First let us study the existence of real roots for the characteristic equation. For the uniqueness of the root, we study the equation 0 = f ( ) . Here

3
Delay-induced blow-up in a planar oscillation model Therefore, one sees that g 1 is a downward-convex function (attaining minimum at = log ) and g 2 is upward-convex function. Hence, the intersection of the functions g 1 and g 2 except 0 is unique. Thus we obtain the conclusion. ◻ For each j ∈ ℕ we obtain the estimation of 4 j with respect to 1 as follows.

Lemma 5.2
The following statements are true.
-It holds 4 0 < 1 and 4 1 > 1 for < * and 4 0 = 4 1 = 1 for = * , and -it holds 4 Proof For any j ∈ ℤ , from the condition (4.4), Thus holds. One can see that for < * and for = * . Thus we obtain the first statement. For j ∈ ℤ ⧵ {0, 1} it is clear that Thus we obtain the conclusion. ◻ From Lemmas 5.1 and 5.2 and the principle of linearized stability [4], we obtain the following result concerning instability of the periodic solutions.

Theorem 5.1
The periodic solution of the form (4.1) for j = 1 is unstable for < * . The periodic solution of the form (4.1) for j ∈ ℤ ⧵ {0, 1} is unstable.
Let us consider stability of the periodic solution of the form (4.1) for j = 0 . For j = 0 the characteristic equation (5.4) may have an imaginary root with positive real part. Thus so far we cannot determine stability of the periodic solution from Lemmas 5.1 and 5.2. Below we exclude this possibility to conclude that the periodic solution of the form (4.1) for j = 0 is asymptotically stable for < * .
First let us show that there is a compact region in the complex plane for the existence of a root with positive real part. By the straightforward estimation |e − | < 1 and | | 1 − e −2 | | ≤ 1 + | | e −2 | | < 2 using Re > 0 , we obtain the estimation (5.6). From the inequality (5.6), one gets from which we obtain the conclusion. ◻ Substituting = + i , ( , ) ∈ ℝ 2 into (5.4), we obtain the following two equations.
From (5.7b) it is immediate to obtain the following lemma. We omit the proof.  Note that we have 0 ≥ 1 from (4.6). Since, for any a ∈ ℝ , it holds that 1 − 2 a + 2 2 a 2 > 0, we see that   which is a contradiction to (5.9b). Therefore, we obtain the conclusion. ◻ Therefore, from the application of the Rouche's theorem (see e.g. Lemma 2.8 in Chapter XI of [5]), we obtain the following conclusion.

Numerical simulations
In this section we demonstrate numerical solutions of the delay differential equation (2.1) with the special initial conditions First we fix = 0.392 . For < * ≈ 0.398284 ⋯ , we show that the periodic solution of the form (4.1) with j = 0 is asymptotically stable (Theorem 5.2). From (4.4), we can compute the radius of the stable periodic solution as r = 0 ≈ 1.18547 ⋯ . In Fig. 7, a trajectory of the solution for = −36 is plotted in (x, y) plane, which shows that the periodic solution with the radius r = 0 ≈ 1.18547 ⋯ attracts the solution. Observe that there is an unstable periodic solution in the vicinity of the stable periodic solution (in Fig. 7, the trajectory of the asymptotically stable periodic solution and of the unstable periodic solution are illustrated as the dashed orange circle and as the dashed blue circle, respectively). In Fig. 8, a trajectory of the solution with the initial condition = −37 is plotted. In this case, the solution winds around the unstable periodic solution (dashed blue circle) and then leaves and goes far away.
In Fig 9, we plot log r(t) for several initial conditions . The numerical experiment suggests that the solution exists globally for −36 < ≤ 0 and blows up for < −37 : in short, the solution blows up for large | | . We can numerically observe many blow-up solutions, which blow up even after t = ∕2 . In this paper, we prove the existence of blow-up solutions, which blow up in the time interval (0, ∕2) . The numerical simulation suggest that many solutions blow up in finite times.
We also numerically compute the solution for > * . Figure 10

Discussion
In this paper we show an example of a blow-up phenomenon in a planar system of delay differential equations. In our system, the delay completely changes the system: many blow-up solutions and periodic solutions suddenly appear, no matter how small the length of the delay is. In neutral delay differential equations, it is known that arbitrary small delay can destabilize the system (see Chapter 1.7 in [9]). Here we find that arbitrary small delay can induce blow-up solutions in delay differential equations. Numerical simulations suggest that many solutions either tend to the stable periodic solution or blow up in a finite time. It is not obvious if more complicated solution behavior exists. The transient behavior observed in Fig. 10 is interesting, as it looks that the solution tries to find a stable periodic solution which does not exist in this parameter setting. Many results concerning the existence of the blow-up solutions are available in Volterra integral equations (see [2,3,15,16] and references therein). When delay Fig. 9 We plot the growth of log r(t) for several . Here is fixed ( = 0.392) . The numerical simulation suggests that larger | | causes the solution to blow up faster Fig. 10 Transient behavior of the solution for = 0.3985 > * . The solution stays around the origin for a while, then it blows up in a finite time differential equations can be formulated as Volterra integral equations, we may apply the blow-up results for Volterra integral equations to delay differential equations, see e.g. [1,21] for a relation between Volterra type integral equations and integro-differential equations. However, we cannot apply the results for Volterra integral equations to a system of integral equations which is rewritten formally from our target system of delay differential equations. Since delay differential equations form an infinite dimensional dynamical system, it is also not straightforward to apply the results established in ordinary differential equations.
In the context of mathematical modelling, our example suggests that arbitrary small delay can be responsible for a drastic change of the dynamics, thus one should be careful when ignoring small delay. Other blow-up mechanisms in delay differential equations will be explored in our future work.

A The characteristic equation
Here, by linearization of the system (5.1) about the equilibrium (5.2), we derive the characteristic equation (5.3), which characterizes stability of the periodic solutions. Fixing j ∈ ℤ , we let Then, applying the Taylor expansion, we get Ignoring the higher order terms, we obtain the following linearized equation Substituing the exponential solution (s(t), u(t)) = e t , where ∈ ℂ and ∈ ℂ 2 , we get the following characteristic equation Using (4.2a), we have cos ∫  1 + e − + − ∫ 0 e − s ds + 4 j 1 + e − = 0.