arxiv: v1 [math.ap] 20 Dec 2018

Size: px
Start display at page:

Download "arxiv: v1 [math.ap] 20 Dec 2018"

Transcription

1 arxiv: v1 [math.ap] Dec 18 Radial solutions of scaling invariant nonlinear elliptic equations with mixed reaction terms arie-françoise Bidaut-Véron, arta Garcia-Huidobro Laurent Véron Abstract We study global properties of positive radial solutions of u = u p + u p in R where p > 1 and is a real number. We prove the existence or the non-existence of ground states and of solutions with singularity at. 1 athematics Subject Classification. 35J6, 35B8, 684. Key words. elliptic equations; ground states; linearization; homoclinic and heteroclinic orbits; Floquet integral. Contents 1 Introduction General properties of the system 8.1 Reduction to autonomous equation and system The standard reduction The geometry of the vector field H Graphic representation of the vector field H Other reduction Regular solutions and ground states Explicit singular solutions Upper estimate of the regular solutions Linearization of the system 1.17) near equilibria Linearization at, ) Linearization at a fixed point P := X, Y ) Energy and Lyapunov functionals for system 1.17) Comparison results Laboratoire de athématiques et Physique Théorique, Université de Tours, 37 Tours, France. veronmf@univ-tours.fr Departamento de atematicas, Pontifica Universidad Catolica de Chile Casilla 37, Correo, Santiago de Chile. mgarcia@mat.puc.cl Laboratoire de athématiques et Physique Théorique, Université de Tours, 37 Tours, France. veronl@univ-tours.fr 1

2 Quasilinear elliptic equations with source reaction 3 Study of ground states of 1.17) when > Behaviour near equilibrium Existence or nonexistence of ground states Study of ground states of 1.17) when < Ground states and large solutions when < The case <, 3 and p The case <, 3 and < p < The case < and 1 < p < The case <, = and p > The case <, = 1 and p > Appendix: on-existence of cycle in the case = Introduction The aim of this article is to study local and global properties of positive radial solutions of the equation u = u u + u p, 1.1) in R or R \ {} where p > 1 and is a real parameter. This is a particular case of the following class of equations u = u u + u q, 1.) where q > 1 which has been the subject or many works in the radial case when <, where a basic observation is that the two terms u u and u q are in competition. The first work in that case is due to Chipot and Weissler [14] who, in particular, solved completely the case = 1, then Serrin and Zou [19] performed a very detailed analysis. uch less is known in the case >. Under the scaling transformation T k defined for k > by 1.) becomes Therefore, if q p, 1.) can be reduced to u k := T k [u]x) = k ukx), 1.3) u k = u k u k + k p q) u k q, 1.4) u = u u ± u p. 1.5) oreover, when q < p, the limit equation of 1.4) when k is the Lane-Emden equation the exponent p is dominant. The other scaling transformation u = u u, 1.6) v k := S k [u]x) = k q q 1 ukx), 1.7)

3 Quasilinear elliptic equations with source reaction 3 transforms 1.) into v k = k q p q) q 1 v k v k + v k q, 1.8) and if q > p, the limit equation of 1.8) when k is the Riccati equation v = v q, 1.9) the exponent q is dominant. In [14] and [19] most of the study deals with the case q p. In the critical case i.e. when q = p p + 1, 1.1) then not only the sign of but also its value plays a fundamental role, with a delicate interaction with the exponent p. otice that an equivalent form of 1.1) is v = λ v v ± v q 1.11) with λ >. In the critical case first studies in the case < are due to Chipot and Weissler [14] for = 1. The case was left open by Serrin and Zou [19] and the first partial results are due to Fila and Quittner [16] and Voirol [, 3]. The case > was not considered. The equation 1.1) is the stationary part of the associated parabolic equation t u u u u u q =. 1.1) which is studied in [14] and [15], one of the aim being to find conditions for the blow-up of positive solutions. A general survey with several open problems can be found in []. In the non radial case an important contribution dealing with a priori estimates of local positive solutions of 1.) and existence or non-existence of entire positive solution in R is due to the authors [7]. In this paper we complete the results of [7] in presenting a quite exhaustive study of the radial solutions of 1.1) for any real number. The radial solutions of 1.1) are functions r ur) defined in, ) where they satisfy u rr 1 u r = u u + u r p. 1.13) r Because of the invariance of 1.13) under the transformation T k there exists an autonomous variant of 1.1) obtained by putting ur) = r xt) with t = ln r, 1.14) then with K = x tt + Lx t K x + x x + )p and L = p K x x t = 1.15) )p + ) = K. 1.16)

4 Quasilinear elliptic equations with source reaction 4 Putting yt) = r u r r), then xt), yt)) satisfies the system x t = x y y t = Ky + x x + y p. 1.17) We are mainly interested in the trajectories of the system which remain in the first quarter Q={x, y) R : x >, y > }. Among these trajectories, the ones corresponding to ground states, i.e. positive C solutions u of 1.13) are defined on [, ). They verify u r ) = and actually they are C on, )). Using the invariance of the equation under T k all the ground states can be derived by scaling from a unique one which satisfies u) = 1. Since it is easy to prove that such a solution u is decreasing, in the variables x, y), a ground state is a trajectory of 1.17) in Q, defined on R and satisfying denoted by T reg. lim t yt) xt) =. The corresponding trajectory is Contrarily to the case of the Lane-Emden equation 1.6), there exists no natural Lyapunov function when. This makes the study much more delicate and it is based upon a phase plane analysis. The solutions of 1.13) invariant under T k for any k > correspond to constant solutions of 1.15) and have the form where X is a positive root of X + Ur) = Xr for all r >, 1.18) ) p X + K =. 1.19) This equation plays a fundamental role in the description of the set of solutions of 1.13). The following constant, defined for = 1, and p > 1 or 3 and 1 < p has an important role in the description of the set roots of 1.19), ) p )p µ := µ ) = p + 1). 1.) p This set is described in the following proposition. Proposition 1 1- If, equation 1.19) admits a positive root, necessarily unique, if and only if 3 and p >. - If < and p, equation 1.19) admits a unique positive root X. 3- If < and either = 1, and p > 1 or 3 and 1 < p, there exists no positive root of 1.19) if µ < <, a unique positive root if = µ < and two positive roots X 1, < X, if < µ. We also set Y = X and P = X, Y ) resp. Y j, = X j, and P j, = X j,, Y j, ), for j=1,) and define the corresponding singular solutions U r) = X r resp. U j, r) = X j, r ).

5 Quasilinear elliptic equations with source reaction 5 Recall briefly the description of the positive solutions of the Lane-Emden equation 1.6), i.e. = : there exists radial ground states if and only if 3 and p + +. If p = these ground states are explicit and they satisfy lim r r ur) = c >. There exist infinitely many singular solutions u ondulating around U. ote that a ground state corresponds to a homoclinic orbit at for system 1.17) and these singular solutions are cycles surrounding P. We recall that an orbit of 1.17) which connects two different equilibria resp. the same equilibrium) when t R is called heteroclinic resp. homoclinic). Without the radiality assumptions, and using a delicate combination of refined Bernstein techniques and Keller-sserman estimate we have obtained in [7, Th. C, Th. D] a series of general a priori estimates for any positive solution of 1.1), in an arbitrary domain of R in the case 1, p > 1 and q = > := no ground state. ) p ) 4p and >. In particular we proved there that if p > 1 and ) p, or if, p < +3 1 and > equation 1.1) admits In the sequel we describe the ground states and the singular global solutions of 1.13) in R \ {}. Concerning the ground states, we discuss according to the sign of and the value of p. The next value of appears when we linearize the system 1.17) at the equilibrium P, =, p) = p + 1) )p ) 4p) p )) + 4) ) Then is positive resp. negative) if p > + + resp. p < and we set µ = ). It is easy to see that if = then the characteristic values of the linearized operator at P are purely imaginary. otice that is positive resp. negative) according p > + resp. p < + ). Theorem A Let 1, p > 1 and >. 1- If 1 < p + if 3, and any p > 1 if = 1,, then equation 1.13) admits no ground state. - If 3 and p > +, there exist constants min, max verifying < < min max, such that - if < < min there exist ground states u satisfying ur) U r) when r. - if = min or = max there exists a ground state u minimal at infinity, that is satisfying lim r r ur) = c >. - for > max there exists no radial ground state. The values of min and max appear as transition values for which the ground state still exists but it is smaller than the others at infinity; it is of order r instead of r. They are not explicit but they can be estimated in function of and p. It is a numerical evidence that min = max in the phase plane analysis of system 1.17) and we conjecture that this is true. When = min or max, the system 1.17) admits homoclinic trajectories. We prove that the system 1.17) admits a Hopf bifurcation when =. When p > + we also prove the existence of different types of positive singular solutions

6 Quasilinear elliptic equations with source reaction 6 Theorem A Let If < p + for any > there exist a unique always up to scaling transformation) positive singular solutions u of 1.13) satisfying lim r ur) = X and lim r ur) = c > r r - If p > +, then i) If > max, there exists a unique singular solution u of 1.13) with the same behaviour as in 1. ii) If < < min there exist positive singular solutions u ondulating around U on R. In terms of the system 1.17) the 1) and -i) correspond to the existence of a heteroclinic orbit in Q connecting P to, ) and ii) to the existence of a cycle in Q surrounding P. When is negative, the precise description of the trajectories of 1.17) depends also on the value of p with respect to +. It is proved in [7, Th. B, E] that for 3 and 1 < p < there exists ɛ > such that if ɛ equation 1.1) admits no positive solution in R. The same conclusion holds if 3, 1 < p or = and p > 1) and > µ. We first consider the case p for which there exists a unique explicit singular solution U, and the results present some similarity with the ones of Theorem A. Theorem B Let 3, p and <. Then 1- If p +, then equation 1.13) admits ground states u. oreover they satisfy ur) U r) as r. - If p < +, there exist numbers µ min and µ max verifying < µ < µ min µ max < µ 1), such that i) for < µ max there exist ground states u such that ur) U r) when r. ii) for = µ min or for = µ max there exist ground states minimal at infinity in the sense that ur) cr when r, c >. iii) for µ min < < there exists no radial ground state. Here also the value of µ min, µ max are not explicit and we conjecture that they coincide. The next result presents some similarity with Theorem A. Theorem B Let 3 and < p < +. i) If < < there exists a unique up to scaling) positive singular solution u of 1.13), such that ur) U r) when r and ur) cr when r for some c >. ii) If µ min < < there exist positive singular solutions u ondulating around U on [, ) and singular solution ondulating around U in a neighbourhood of and satisfying ur) cr for some c > when r. In terms of the system 1.17), i) corresponds to a heteroclinic orbit connecting P and, ), while ii) to the existence of a periodic solution in Q around P, and the existence of a solution in Q converging to, ) at and having a limit cycle at t = which is a periodic orbit around P.

7 Quasilinear elliptic equations with source reaction 7 The situation is more complicated when 1 < p < and < µ because there exist two explicit singular solutions U 1, and U, which coincide when = µ. Theorem C Let <, 3 and 1 < p < constants µ min and µ max verifying such that, or = and p > 1. Then there exist two µ µ < µ min µ max < µ 1), 1- If < µ max then equation 1.13) admits ground states u either ondulating around U, or such that ur) U, r) as r. - If = µ min or = µ max there exists a ground state u such that ur) U 1, r) as r. 3- If µ min < < there exists no radial ground state. Here again µ min and µ max appear as transition values for which the ground state still exists but it is smaller than the others at infinity: it behaves like U 1 instead of U. The proof of this theorem is very elaborate in particular in the case =. In the case = 1 the result is already proved in [14]. The nonexistence of ground state, not necessarily radial for > µ is proved in [1] and independently in [7] with a different method. In the radial case it was obtained much before in the case = 1 in [14] and then by Fila and Quittner [16] who raised the question whether the condition µ min < < is optimal for the non-existence of radial ground state. This question received a negative answer in the work of Voirol [] who extended the domain of non-existence to µ ɛ < µ. The next result is the counterpart of Theorem C when dealing with singular solutions. Theorem C Let <, 3 and 1 < p <. i) If < µ there exist positive singular solutions u such that ur) U 1, r) as r and ur) cr with c > when r. ii) If < µ there exists a unique up to scaling positive singular solution u, such that ur) = U, r) as r and ur) = U 1, r) as r. Furthermore ur) > U 1, r) for all r >. iii) If ˆµ min < < µ there exist positive singular solutions u ondulating around U, at and such that ur) U 1, r) as r, and positive singular solutions u ondulating around U, on R. iv) If = µ min or = ˆµ max there exists a positive singular solutions u different from U 1, such that ur) U 1, r) when r and r. v) If µ min < < ˆµ max there exists a positive singular solution u such that lim r ur) = r c > and either ondulating around U, or such that ur) U, r) when r. vi) If 3 and = µ, there exist positive singular solutions u satisfying lim r ur) = r c > and ur) U µ r) as r. In terms of the system 1.17) i) corresponds to a heteroclinic orbit connecting P 1, to, ); ii) to a heteroclinic orbit connecting P, to P 1, ; iii) to a trajectory having a periodic orbit around P, for limit cycle at and converging to P, at and to a periodic orbit around

8 Quasilinear elliptic equations with source reaction 8 P, ; iv) corresponds to homoclinic orbit at P 1, ; v) corresponds to a trajectory connecting, ) at and either converging to P, at or having a periodic orbit around P, for limit set at ; vi) corresponds to a heteroclinic orbit connecting from, ) to P µ. Acknowledgements This article has been prepared with the support of the collaboration programs ECOS C14E8 and FODECYT grant for the three authors. General properties of the system.1 Reduction to autonomous equation and system.1.1 The standard reduction We recall that if u is a C 3 function defined on some interval I [, ) verifying 1.13) and if then x satisfies the autonomous equation ur) = r xt) with t = ln r, x tt + Lx t K x + x x + p x x t =,.1) on lni) where K and L are defined in 1.16). Setting u r = r yt), then xt), yt)) satisfies { xt = H 1 x, y).) y t = H x, y), where H 1 x, y) = x y H x, y) = Ky + x x + y p. and we denote by H the vector field of R with components H 1 and H..3).1. The geometry of the vector field H Let us denote by Q := {x, y) : x >, y > } the first quadrant. The vector field is inward in resp. outward of) Q on the axis {x, y) : x >, y = } resp. {x, y) : x =, y > }). We set { L := x, y) Q : y = x } { ) and C := x, y) Q : x = Ky y p 1 } p.4) ) and ψy) = Ky y p 1 p. Then x t = on L and y t = on C. The curves L and C have [ zero, one or two intersections in Q according the value of K and. If, K >, then ) 1 ) C, p p K ] [ ) p) p, ] ) K. The points, ), P and, ) K ) belong to C. The function ψ is increasing on, K q ) ) ) and decreasing on K q, ) K ).

9 Quasilinear elliptic equations with source reaction 9 If < and K >, ψ is concave and increasing on, ) with ψy) = ) 1 p y 1 + o1)) as y. If < and K <, ψ is still concave and increasing on ) K, ) with the same asymptotic as above. We quote below the possible connected components of Q \ L C). { } } x A = x, y) : y < {x, y) : Ky + x p + y p < = {x, y) : x t <, y t < }. { } } x B = x, y) : y > {x, y) : Ky + x p + y p < = {x, y) : x t >, y t < }. { } } x C = x, y) : y > {x, y) : Ky + x p + y p > = {x, y) : x t >, y t > }. { } } } x D = x, y) : y < {x, y) : Ky + x p + y p > {x, y) : x > X }, = {x, y) : x t <, y t > } {x, y) : x > X,. { } } } x E = x, y) : y < {x, y) : Ky + x p + y p > {x, y) : x < X } 1, = {x, y) : x t <, y t > } {x, y) : x < X 1,. These connected components are A, B, C, D if K, < or K, >. A, C, D if K < and µ < <. A, C, D, E if K < and = µ. A, B, C, D, E if K < and < µ..1.3 Graphic representation of the vector field H We present below some graphics of the vector field H associated to system 1.17). 3 y C D) x t < y t > L A) x t < y t < P C) x t > y t > B) x t > y t < o x Figure 1: >, K >.

10 C) x t > y t > Quasilinear elliptic equations with source reaction 1 o x y C D) x t < y t > A) x t < y t < L P x t > y t < C) x t > y t > B) o x Figure : <, K. y C D) x t < y t > A) x t < y t < L C) x t > y t > o x y A) x t < y t < D) Figure 3: µ < <, K <. C x t < y t > L P x t > y t < C) x t > y t > B) o x

11 x t > y t < B) C) x t > y t > x E) t < y t > P 1, Quasilinear elliptic equations with source reaction 11 o x y C D) x t < y t > A) x t < y t < L E) x t < y t > P µ C) x t > y t > o 1 x Figure 4: = µ, K <. y A) x t < y t < C D) x t < y t > L P, x t > y t < B) C) x t > y t > x t < y t > E) P 1, o x.1.4 Other reduction y Figure 5: < µ, K <. A) x t < y t < The following change of unknowns, already used in [9] when =, C D) x t < y t > L σt) = yt) xt) = r u rr) ur) E) x t < y t > and zt) P µ = x xt) = r u ur),.5) yt) u r r) C) x t > y t > o 1 x

12 Quasilinear elliptic equations with source reaction 1 and valid if u r, transforms 1.17) into a Kolmogorov system with vector field V = V 1, V ) σ t = σ σ + + z + σ 1 ) σz = V 1 σ; z) z t = z pσ z σ 1 ).6) σz = V σ; z). Since σ and z are in factor the two axis {σ = } and {z = } are trajectories, actually not admissible for.) in view of.5). The system is singular on these two axis however it can be desingularized in setting σ = σ k+1 and z = z k+1 for some integer k > p + 1, which transforms.6) into a new nonsingular Kolmogorov system, σ t = 1 k + 1 σ z t = 1 k + 1 z σ k z k+1 + p σ k+1 z k+1 σ σ) k+1 z 1 ) k+1 σ σ) k+1 z 1 ) k+1. Therefore no other trajectory can intersect them in finite time and the quadrant Q := {σ, z) : σ >, z > } is invariant. Furthermore σz = r u. It is noticeable that if = the initial system is quadratic and regular. The system.6) will be used in the most delicate cases. It corresponds to the differentiation of the initial equation.1).. Regular solutions and ground states Definition.1 A regular solution of 1.13) is a C solution defined on some maximal interval [, r ) satisfying u) = u > and u r ) =. A ground state is a nonnegative C solution defined on [, ). The existence and uniqueness of a regular solution is standard by the Cauchy-Lipschitz integral method. If u is a C solution it satisfies u r < on, r ). Indeed r 1 u r r) is decreasing near, hence u r < on some maximal interval, r 1 ), r ) and u r r 1 ) = if r 1 < r. If ur 1 ) = then u by uniqueness. If ur 1 ) > then u rr < near r 1 which would imply that ur) < ur 1 ) for r 1 ɛ r < r 1 which contradict the negativity of u r on, r 1 ). Hence ur 1 ) < which implies that u r r) < on the maximal interval, r ) where u >. Thus, if u is a ground state u r < on, ). Hence the trajectory of a ground state expressed in the system 1.17) lies in Q and expressed in the system.6) it lies in the quadrant Q. It is easy to check that the regular solution u := u u, such that u) = u satisfies ) ur) = u 1 u r p + 1) u r ) 1 + o1)) as r..8) 4p + ) + )p + ) q Under the scaling transformation T k, u can be transformed into the regular solution 1.13) u := u 1 satisfying u) = 1. If one considers the system 1.17) the transformation T k becomes the time shift which transforms t xt), yt)) into t xt + ln k), yt + ln k)), and the.7)

13 Quasilinear elliptic equations with source reaction 13 trajectory xt), yt)) of 1.17) corresponding to a ground state is therefore uniquely determined and denoted by T reg and satisfies yt) lim xt), yt)) =, ), lim t t xt) = yt) and lim t x p t) = 1..9) Hence in the system.5) there holds on the corresponding trajectory lim σt) =, lim zt) =..1) t t.3 Explicit singular solutions Explicit self-similar solutions of 1.13), necessarily under the form u = Ar, play a fundamental role in the study, whenever they exist. The following result covers Proposition 1. Proposition. 1- Let. Then there exists a unique self-similar solution of 1.13) if and only if 3 and p >. We denote it by U r) = X r, where X > depends also on, and p. To this solution is associated the equilibrium X, X ) of the system 1.17). Furthermore the mappings X is continuous and decreasing on [, ), X increasing and there holds i) X = ii) ) 1 K, ) K 1 1 ) p ) K p ) + X is ).11) K. - Let <. If 3 and p there exists a unique self-similar solution of 1.13) U r). The mapping X is continuous and decreasing, X holds { ) 1 ) } K max, p) X is decreasing and there ) 1 ) ) K + p)..1) 3- Let <. If = 1, and p > 1, or 3 and 1 < p <, there exists no self-similar solution of 1.13) if µ < < where µ = µ ) > is defined in 1.). If = µ there exists a unique self-similar solution U µ r). If < µ < there exist two-self-similar solutions U 1, r) = X 1, r and U, r) = X, r with X 1, < X,. Furthermore the mappings X 1, and X, are continuous, respectively increasing and decreasing on

14 Quasilinear elliptic equations with source reaction 14, µ ), while X1, and X, Furthermore there holds for large enough; are respectively decreasing and increasing. i) X j, µ = ii) iii) ) p K p ) K < X1, < ) ) p) ) 1, ) K K 1 1 p ) 1 K < X, < )K ) ), ) ) p)..13) Proof. The function U = X r is a self-similar solution of 1.13) if and only if X is a positive root of ) p f x) := x + x K =..14) Equivalently P = X, Y ) = Y, Y ) is a fixed point of system 1.17), where Y is the positive root of ) p f y) = y + y K =..15) The use of the variable y is a little easer than x. Since X is explicit if =, we shall study the cases. 1- Case >. If K >, equivalently p >, and, f is an increasing function tending to at and negative at y =. Hence Y is the unique positive root of.15). If K <, f is positive on [, ), hence no such solution exists. Since ) p Y + Y K =, by the implicit function theorem, Y is C 1. For > >, f y) > f y) for all y >. Hence Y is decreasing on [, ). Actually the expression of f shows more, namely that Y and ) p ) K + Y is increasing on [, ). Furthermore Y < K = Y < > K = Y > ) K ) K, 1 ) p ) 1 K p ), +

15 Quasilinear elliptic equations with source reaction 15 and we get.11). - Case <. Clearly f has a minimum at y = y, where y, = p + 1 ) 1 ) p and f y, ) = p ) p K..16) p + 1 -i) If 3 and p, then f ) <, f is decreasing on, y, ) and increasing on y,, ), hence f y) = has a unique positive root Y > y,. Since for < < and y >, f y) > f y), the mapping Y is continuous and decreasing and Y is increasing. Then the left-hand side of.1) is clear. ext we put A p = ) p, B = K, ξ = X, a = B A and η = ξ A. Then φη) = η η a =. Since φ1 + a) = 1 + a) 1 a a 1 + a) 1 a a p a 1 + a) 1 + a) + a p )) we derive η 1 + a, which implies the right-hand side of.1). -ii) If = 1, or 3 and 1 < p <, then K <. Hence, if f y, ) >, or equivalently µ <, the equation f y) = has no positive root, if = µ, it has a double root Y µ where Y µ = ) p K p ) 1 = µ Y µ = p + 1 p K,.17) and if < µ the equation.15) has two positive roots < Y 1, < Y, and since f does not vanish at Y j,, they are C 1 functions of, ), respectively increasing and decreasing. Since, K <, we obtain from f Y ) =, and ) p Y Y 1, > ) K,, < Y, = Y, < For a sharper estimate, we have for large enough, Ỹ = ) K 1 K )K ) ) p) ) ) = f Ỹ ) <.

16 Quasilinear elliptic equations with source reaction 16 Hence ) ) K K < Y1, < 1 K )K ) ). Similarly ) ) p) 1 K 1 p ) < Y, < ) ) p). The estimates.13) follow..3.1 Upper estimate of the regular solutions We first recall the following estimate in the case, consequence of the fact that the positive solutions of 1.1) are superharmonic and proved in a more general setting in [7, Prop..1 ]. Proposition.3 1- There exists no positive solution of 1.13) in R, ), R if and either = 1, and p > 1 or 3 and 1 < p. In particular there exists no ground state. - If 3, p >, and u is a positive solution of 1.13) in R, ), R. Then there exists ρ R such that and ur) min u r r) min { ) 1, { ) ) 1 p ) ), p ) ) Furthermore, if R =, inequalities.18), and.19) hold with ρ =. ) } r, for all r > ρ,.18) ) } r for all r > ρ..19) The next estimate is verified by any ground state, independently of the sign of. Proposition.4 Let p > 1 and 1. Then the ground state u of 1.13) with u) = 1 satisfies } ur) min {1, c,p, r and u r r) c,p, r r >..) Proof. The trajectory T reg starts from, ) and enters the region C. If it stays in C, then x is increasing on R. Since y < x, if xt) tends to some finite limit as t, it implies that the limit set of the trajectory exists. It cannot be a cycle since x is monotone, thus it is one of the equilibrium of the system. Hence xt) X j, for j = 1 or if K < and < µ or xt) X if K and either > or µ. This implies that.) holds. If xt) tends to, so y does because yt) x p t) K xt) xt) ) p as t.

17 Quasilinear elliptic equations with source reaction 17 Therefore y t Cy p for some C > which would imply that t y 1 p t)+c)t is increasing, which is impossible. ext we suppose that the trajectory leaves C by crossing the line L. Since it cannot enters B through C in the case K <, µ ), it leaves C by intersecting L and we denote by t 1 the first time where T reg intersects L. Then x t t 1 ) = and x tt t 1 ) = y t t 1 ) <. Therefore t 1 is a local maximum. ow the trajectory cannot cross again the half line { x, y) : x = xt 1 ), y > xt 1) for any t t 1. } because on it there holds x t = xt 1) y <. Hence xt) xt 1) In the same way, either y is increasing on R and since x t = x y and x is bounded, y cannot tend to infinity when t, thus yt) y >, or y is not monotone and T reg crosses C at a first value t, necessarily larger than t 1 and where x t t ) <. Then y tt t ) = p xt ) x t t ) < and t is a local maximum of y. Therefore t t xt), yt)) remains in the subset of Q bordered by the portion of trajectory of T reg for t t and {x, y) : < x xt ), y = yt )}. This implies that yt) yt ) for all t R. oticing that ur) 1 since u is decreasing, we get the conclusion. Remark. The above method does not give an explicit estimate of the upper bounds of x and y and such a bound can be estimated in some cases. If it follows from Proposition.3 that for any p > 1 there holds ) 1 ur) r, thus this estimate is independent of. If µ ) < <, it is easy to check that the function v = ln u satisfies v ae )v.1) with a = 1 µ )) >. We derive that the function wr) = r 1 v r r)) is increasing, with limit l, ]. Hence This implies e 1 p)vr) e 1 p)v) + a)r 4 } ur) min {1, c a 1 r..).4 Linearization of the system 1.17) near equilibria.4.1 Linearization at, ) The linearized system at, ) is x t = x y y t = Ky,.3)

18 Quasilinear elliptic equations with source reaction 18 which admits the eigenvalues ote that λ λ 1 =. λ 1 = K λ =..4) a) Assume that 3 and p >, equivalently K >, Then, ) is a saddle point. There exists a unique unstable trajectory T unst such that yt) lim xt), yt)) =, ) and lim =,.5) t t xt) and more precisely, from.8), xt) = e t 1 et p + 1) 4p+ ) e t 1 + o1)) p + 1) + p) p + 1) ) yt) = e pt 1 p + 1) + p + 1) + p e )t 1 + o1)) as t..6) From Definition. and the lines which follow, the unstable trajectory T unst coincides with the regular trajectory T reg. This is included in the region C for t < T for some < T. There exists also a unique stable trajectory T st such that yt) lim xt), yt)) =, ) and lim t t xt) = + K =..7) Since >, T st belongs to the region A for t > T 1 for some T 1 <. Remark. If T reg Q satisfies lim xt), yt)) =, ), then T reg = T st and the corresponding t solution is a ground state. The same conclusion holds, if T st Q satisfies xt), yt)) = lim t, ). Such a solution is called a homoclinic orbit at, ). Because of the uniqueness of the stable and unstable trajectories of a saddle point, it is unique in the class of solutions satisfying.7). Equivalently the class of ground states u of 1.13) satisfying ur) cr for some c > is then a one parameter family characterized by u) = u. b) Assume that 3 and 1 < p <. Then K < and < λ 1 < λ. Hence, ) is a source and all the trajectories of 1.17) in some neighbourhood of, ) converge to, ) when t. Among those trajectories there exists one fast trajectory which satisfies yt) lim t xt) =. It is actually the regular trajectory T reg. There exist also infinitely many slow trajectories which satisfy yt) lim t xt) =. c) If p =, then K = and λ 1 = < λ =. We still find the regular trajectory T reg associated to λ and the corresponding eigenvector 1, ). By the central manifold theorem

19 Quasilinear elliptic equations with source reaction 19 corresponding to λ 1 there exists an invariant curve passing through, ) with slope. Using the matched asymptotic expansion method, one finds that if < there exists a solution x of.1) such that xt) C 1 1 t 1 when t, i.e. ur) C 1 1 r ln r)) 1 when r, and if > there exists a solution x of.1) such that xt) C 1 1 t) 1 when t, equivalently ur) C 1 1 r ln 1 r ))1 when r. d) If =,, ) is a source with λ 1 = λ =, with corresponding eigenspace 1, ). The linearized problem is equivalent to equation x tt 4 q 1 x 4 t + q 1) x = with general solutions xt) = ae t q 1 + bte t q 1 for some real parameters a, b. Hence there exists infinitely trajectories of.1) tending to when t and they are tangent to 1, ) at, ). The regular trajectory T reg corresponds to b = and the other trajectories correspond to singular solutions u of 1.1). They satisfy ur) b ln r as r and there holds u = u p + u p πbδ, in the sense of distributions in B ɛ for some ɛ >. ext we give a general result in case the system admits only one equilibrium in Q. and R. If u is a regular solution the following tetra- Lemma.5 Let 3, p > chotomy occurs: i) either lim r r ur) = c for some c >, ii) or ur) U r) as r, iii) or ur) has an ω-limit cycle surrounding P, iv) or ur) changes sign for some r >. Proof. By assumption P is the unique equilibrium. The trajectory T reg starts from, ) and remains in the region C where x t, y t > for t t. If t =, u is a ground state, hence it is bounded from Proposition.4. Its ω-limit set is non-empty. Because x and y are monotone, it converges when t to some point which is necessarily P. If t <, then at t = t the trajectory leaves C through L since it cannot enter in B, and it enters the region D where x t <, y t >. oreover x t t ) = and xt ) > X. Then three possibilities occur: α) either xt) X monotonicaly when t ; thus the trajectory converges to P. β) either xt) monotonicaly when t. Since, ) is a saddle point, then T reg =T st. This implies that T reg is a homoclinic trajectory at, ). γ) or there exists t 1 > t such that x t t 1 ) =. Then xt 1 ) < X. Hence T reg enters the region B and by continuity there exists t < t such that xt ) = xt 1 ) and yt ) < yt 1 ). Therefore the bounded region of R bordered by the segment I = {x, y) : x = xt 1 ), yt ) < y < yt 1 )} and the portion of T reg defined by {xt), yt)) T reg : t t t 1 } is positively invariant notice that x t > on I) and it contains P and no other equilibrium. Therefore either the trajectory converges to P or it admits an ω-limit cycle which is a closed orbit surrounding P.

20 Quasilinear elliptic equations with source reaction.4. Linearization at a fixed point P := X, Y ) Suppose that P or P j, ) exists. Then setting x, y) = X, Y )+x, y), the linearized system at this point is x t = x y y t = px x + p Y K ) y,.8) The eigenvalues of its matrix are the roots of ) p T λ) = λ p + 1 Y L λ + K p p + 1 Y,.9) using equation.15). If is such that p p = L, then K = K L = K + =, and we denote by such a value of which is characterized by Y Y p p + 1 Y = L..3) Since Y is a positive root of.15), we get ) p Y = A := K p + 1 p L = )) + 4, p) hence, well defined for, is given by = p + 1)L p 1 A ) p ) 1 p + 1) )p ) =..31) 4p) p )) + 4) 1 That is 1.1). Remark. We see that > resp. < ) if and only if 3 and p > + resp. 1 < p < + if 3 and any p > 1 if = ). If = and 3, the eigenvalues of the linearized system are purely imaginary. If p = +, then =. It is known that in that case the point P = ), ) ) is a center for the system 1.17) associated to u = u + and that there exist infinitely many cycles turning around P with equation 1 y + x xy = E..3) It can be verified that E <, in particular using the function F defined in.35). For and 1 < p <, there always holds µ and more precisely Lemma.6 If 3 and 1 < p <, then < µ. If = and p > 1, then = µ.

21 Quasilinear elliptic equations with source reaction 1 Proof. From Proposition. and identities.17) and.9), there holds, µ < µ µ Y µ < µy K < L. But L K = series of equivalences. and the conclusion follows. If = we just replace < by = in the above Lemma.7 Assume 3. If p > + and = > then P is a weak sink and a Andronov-Hopf bifurcation point. If p < + and = < then P is a weak source. Proof. We recall that a weak sink is an asymptotically stable equilibrium which attracts the nearby points as t at a rate slower than the usual exponential rate. A weak source is a weak sink of the system obtained by changing t into t see [17, Chap. 9]). We write x = x X, ȳ = y Y and obtain the new nonlinear system where with c 1 = x t = x ȳ ȳ t = px x + p Y h x, ȳ) = c 1 x + c ȳ + c 3 x 3 + c 4 ȳ , p) X p, c p) = p + 1) Y ȳ Kȳ + h x, ȳ),.33) p)p ) c 3 = X p 3 6, c p) 4 = 3p + 1) 3 Y Setting α = 4 +, we have px ) = α, since X = Y and satisfies.3), thus c 1 = α )L )L, c =, c 3 = Y p + 1)Y 3p + 1) Y..34) In order to compute the Lyapunov coefficients we transform the system by setting The new system is γ =, s = γt, γw = x ȳ. x s = w ) w s = x 1 h x, x γ γw. By integrating the first line and using the expansion of h, we obtain [ w s = x 1 ) x ) x 3 γ c 1 x + c γw + c 3 x 3 + c 4 γw +... p+3.

22 Quasilinear elliptic equations with source reaction This can be written in the following way w s = x + ν, x + ν 1,1 xw + ν, w + ν 3, x 3 + ν,1 x w + ν 1, xw + ν,3 w By [17, Th. 9..3], the Lyapunov coefficient is given by Λ = ν, + 3ν,3 + ν 1,1 ν, + ν, ) which yields by computation γλ = α 3c 4 4c ) c ) + 1) α + c 1 ) = α p + 1) Y L. If L > resp. L < ) P is a weak sink resp. a weak source)..4.3 Energy and Lyapunov functionals for system 1.17) If xt) is a solution of.1) we set Then, if xt), yt)) Q, we have x F t t) = Lx t x t ) ) x x = L y + y q Hence, if L >, F is monotone in the region F t) := x t + x p + 1 Kx,.35) ) p xt = Lx t + y q ) x t ) y. { x, y) Q : x y > }, located under L. Remark. a) This function was introduced classically in the case =, leading in particular to.3) when p = +. b) With this function we can deduce an upper estimate of regular solution, completing, Proposition.3 and Proposition.4, namely if > and p + any ground state satisfies ) 1 p + 1 ur) K r for all r >..36) ext we construct a Lyapunov functional adapting the method initiated by [] and already used in [3] and [4]. Lemma.8 We define on R J x, y) = Kx x p + 1 ) p 3 p + 1) x 1 3p + 1 ) x y,.37) and, if xt), yt)) is a solution of.38), Vt) = J xt), yt)). If xt), yt)) = x, y) Q we have ) [ ) ) p )] x x x V t t) = y L y p y ) x ) p ) ).38) x = L y p y x y.

23 Quasilinear elliptic equations with source reaction 3 Consequently, the function t Vt) is decreasing if > and 1 < p +, and increasing if < and p +. Proof. We recall the ansatz introduced in [] for finding a Lyapunov function for a system of the form { xt = fx, y).39) y t = gx, y). If fx, y) = y = hx), then consider the function Lx, y) = y hx) fx, t)dt x gt, ht))dt..4) In the case of system 1.17), hx) = x, and we find Lx, y) = J x, y). Then.38) and the conclusion follow. Remark. If L > we set H = x, y) Q : ) p x y p x y = L ) 1, or x L = y = q q 1. Then R := θ 1 V t t) = xt), yt)) H { x, y) Q : oreover H is starshapped with respect to and we set θh = x, y) Q : } x = y. ) p x y p x y L ) 1, or x L = y q q 1. If > we have V t t) if xt), yt)) R and V t t) if xt), yt)) R c Q. If <, the signs of V t t) in the same regions are opposite. to In [7] we also used a function introduced in [19] for equation 1.). When q = p it reduces p + 1 Zr) = r a u r + u + a uu ) r r + u u r p,.41) where a = ) 1) p+3. Since r = e t, we find that in Q, Zr) = e )L p+3 t The function satisfies the relation ) p + 1 y + x axy + xy p..4) Z r p p + 1 u r Z = r a 1 U,.43)

24 Quasilinear elliptic equations with source reaction 4 where U = 1)) uu r p + 3) r = 1)) p + 3) xy pp + 3) L + L + ote that U has a constant sign in Q if L <. p + 1) ru r pp + 3) y p + 1) ) ) r p+3.44).4.4 Comparison results Lemma.9 Let 1, p > 1 and, R such that <. Then, as long as they lie below the line L, i.e. x t >, the regular trajectories T reg and T reg associated to and respectively do not intersect. Furthermore T reg is below T reg. Proof. We use the expansion.6) and deduce, for i = or, y i t) = e t 1 x i t) + i e t + p + o1)) )1 ).45) as t. Hence T reg is below T reg for t t for some t R. Suppose that the trajectories intersect for a first time at some point x, y ) below L. Since the system is autonomous, there will exist two solutions of the systems relative to and satisfying at the same time t, x t ) = x t ) = x and y t ) = y t ) = y. From.), x t t ) = x t ) > and p y t t ) = y t t ) + )y < y t t ). Hence the intersection of the two trajectories is transverse and the slope of T reg at this intersection point is smaller that the one of T reg which is impossible. The second assertion follows immediately. 3 Study of ground states of 1.17) when > When and either 3 and 1 < p, or = 1, and p > 1, there exists no ground state by Proposition.3; in this section we assume 3 and p > 3.1 Behaviour near equilibrium. Since and p > there exists a unique equilibrium P Q. The sign of the real part of the characteristic roots of the linearization of the system.) at P depends on the position of with respect to defined in.31). Lemma 3.1 Assume > and 3. 1) If < p +, then P is a source. ) If p > + ; then P is a sink when < and it is a source when >. 3) If p > +, then P a weak sink. Furthermore, if < < ɛ, for ɛ small enough, there exists a periodic trajectory in Q surrounding P. t

25 Quasilinear elliptic equations with source reaction 5 Proof. Step 1: Assume < p +. The linearized system is given in.8). Because >, the product of the characteristic roots given by equation.9) is positive since it is given by K p p + 1 Y = p + 1 K + p ) p Y p + 1. The sum or the real part) of the characteristic roots is equal to positive, as p + implies L <. Hence P is a source. Step : Assume p > + p p + 1 Y L which is, hence L >. As in Step 1, the product of the characteristic roots is positive. By Proposition.-1), Y is an increasing function of, then the sum or p the real part of) of the characteristic roots, given by p + 1 Y L is also increasing and vanishes if =. It is negative if < < and positive if >. It implies assertion. Step 3: If =, and p > +, then P is a weak sink by Lemma.7. The appearance of the limit cycle, which is the called the Andronov-Hopf bifurcation, occurs for > when is small enough see [17, Chap. 9]). Remark. The product of the characteristic roots is also expressed by K p p + 1 Y = p + 1 Y + ) Y. Hence it is positive for any, p > 1. ext we give some sufficient conditions for nonexistence of a periodic solution or a homoclinic orbit at, ). Lemma 3. Assume 3. If > and 1 < p + the system 1.17) admits no closed orbit in Q. If < and p > +, the system 1.17) admits no cycle in Q surrounding P. Proof. If γ is a non-trivial closed orbit it corresponds either to a T -periodic solution or a solution such that xt), yt)) = lim xt), yt)). lim t The function V defined in.37) is monotone and it satisfies either V) = VT ) or lim Vt) = t lim Vt). Hence it is constant and by.38) it implies yt) = t xt) for all t, a contradiction. ext we suppose p > +, < and that there exists a T -periodic solution xt), yt) with trajectory γ Q surrounding P, hence P belongs to the bounded connected component Γ of R \ {γ} bordered by γ. Since P is a sink or a weak sink by Lemma 3.1, there exists a neighbourhood O of P such that all the trajectories issued from O converge to P as t. Hence any trajectory issued from O, necessarily contained in Γ, has an α-limit set in Γ which is either a stationary point different from P, which is excluded, or a limit cycle t

26 Quasilinear elliptic equations with source reaction 6 {xt), yt)} t [,τ) := γ Γ τ is its period). This limit cycle is not stable, hence, by Floquet s theory τ H1 xt), yt)) + H ) xt), yt)) τ ) p dt = x y p + 1 y t) L dt 3.1) We perform the change of unknowns x = x X, ȳ = y Y used in Lemma.7 which leads to the system.33). The explicit value of the remaining term is h x, ȳ) = x + X ) p X p px x + ȳ + Y ) p Y = Φ x) + Ψȳ), p p Y where Φ and Ψ are defined accordingly. It is positive by convexity because >. Since from.33), τ τ xt)dt ȳt)dt =, and τ = px xt)dt + p τ τ τ p + 1 Y ȳt)dt K ȳt)dt + h x, ȳ)dt we derive p) p p Y + p ) τ p + 1 Y K ȳt)dt <. Using the equation.15) satisfied by Y it yields ) K ) τ p + 1 Y ȳt)dt <. By.11)-ii), K τ τ p + 1 Y >, hence ȳt)dt <, therefore τ 1 yt)dt < Y and by concavity, 1 τ yt)) 1 τ ) dt yt)dt < Y τ τ. 3.3) Combining 3.1) and 3.3), we obtain Since Y which contradicts 3.4). ȳ 3.) L < p p + 1 Y. 3.4) is decreasing by Proposition.-1-ii) we have for = p p + 1 Y L p p + 1 Y L, Remark. Up to changing the sense of variation of Vt), the proof of the first assertion shows that there exists no closed orbit in Q if < and p +. However the proof of the second assertion is not valid when <. The nonexistence of any periodic solution can also be proved when the equilibrium is a node i.e. the two characteristic values are real with the same sign).

27 Quasilinear elliptic equations with source reaction 7 Lemma Let 3 and p > +. There exists a unique and explicit > such that for any, P is a repelling node, degenerate if =. If < <, P is a repelling spiraling point. If 3 1 and < <, P is an attracting spiraling point. - If 11 and + < p < there exists a unique and explicit < 1 < such that if P is an attracting node if < < 1 and an attracting spiraling point if 1 < <. 3- If > or if < < 1, there exist no periodic trajectory in Q around P, neither no homoclinic trajectory in Q at, ) surrounding P. Proof. The characteristic values are real if and only if ) p D = p + 1 Y L 8 p =. K p + 1 Y L p ) p + 1 Y ) p p + 1 Y L + 1 ) 3.5) Hence 3.5) is satisfied if one of the following conditions holds: i) ii) Y p + 1 ) L 1, p Y p + 1 ) 3.6) L + 1. p It is easy to check that for p > +, one has < L + 1 < qk. Since >, by Proposition.-1), the mapping Y [, ). Therefore there exists a unique such that Using.15) we find = p + 1 p is continuous, increasing and from [, ) into Y = p + 1 ) L ) p Concerning the lower bound in i), there holds ) p K p + 1 p K + ) 1 p ) ) L 1 > p > and ) Then we can define 1 by 1 Y 1 = p + 1 ) L ) p

28 Quasilinear elliptic equations with source reaction 8 which leads to ote that p Hence 1Y 1 1 = p + 1 p ) p = L 1 < L = p K p 1 p K + + ) 1 p Y ) ) L + 1 = p Y. 1 < <. 3.1) ext we prove 3) by adapting an argument introduced in [1] for quadratic systems. We return to system.33) that we write under the form { xt = a x ȳ 3.13) ȳ i = c x + dȳ + Φ x) + Ψȳ) which defines a, c, d with Φ and Ψ given by 3.13), and the trinomial.9) for characteristic values endows the form T λ) = λ d + a)λ + ad + c. In the range of values of, the discriminant D = d + a) 4ad + c) = d a) 4c is positive. We consider the intersection of a straight line l passing through P with equation ȳ = A x with a trajectory xt), ȳt)). Then { U = ȳt A x t = A + d a)a + c) x + Φ x) + Ψȳ) = A + d a)a + c) x + Φ x) + ΨA x). We can choose A such that A + d a)a + c = since D >. Since Φ and Ψ achieves positive values, we derive from the expression of h that U > for x. This proves that any closed orbit around P or passing by P can intersect l only one time which is a contradiction. 3. Existence or nonexistence of ground states Proposition 3.4 Let 3, < p + and >. Then there exists no radial ground state, and there exists a singular solutions ur) such that lim r ur) = X r lim r r ur) = c for some c >. Proof. Assume that T reg remains in Q, by Lemma.5 we have three possibilities: α) either T reg converges to P when t, which is impossible since P is a source, β) or T reg has a limit cycle at, and this is impossible by Lemma 3., γ) or T reg converges to, ) when t, hence it is a homoclinic orbit. The function V defined in.37) is decreasing by Lemma.8. Since V ) = V ) =, it is identically and and

29 Quasilinear elliptic equations with source reaction 9 so is V t. This implies that xt) yt) = for all t which is a contradiction. Hence T reg does not remain in Q. We denote by Ξ the connected region of Q bordered by the semi-axis {x, y) : x =, y > } and T reg. Since H is outward on the semi-axis, Ξ is negatively invariant. By Section 3.4.1, yt), ) is a saddle point, hence the stable trajectory T st = {xt), yt)} satisfies lim t xt) = which implies ur) cr for some c >. Its α-limit set αt st ) cannot be a limit cycle as we have seen it above. If it contains, ) it implies again that V, which is monotone, is equal to, hence V t and xt) yt), which is impossible. Hence αt st) contains P. Since P is a source it implies that T st converges to P when t. Proposition 3.5 Let 3 and p > +. Then for any < there exists a ground state u which satisfies u U at. Proof. If < < resp. = ), P is a sink resp. a weak sink). Suppose first that the trajectory T reg does not stay in Q, then it leaves Q at some point, y s ) with y s >. As a consequence, the stable trajectory T st at, ) remains in the negatively invariant region Ξ defined in the proof of Proposition 3.4. Since it cannot converge to P when t it admits a limit cycle surrounding P which contradicts Lemma 3.. Therefore T reg Q and, again using Lemma 3., either it converges to P when t and the proof is complete, or to, ) and T reg = T st is a homoclinic trajectory. The trace of the linearized system.3) at, ) is equal to K = L <. Therefore, from [17, Th ] the connection is attracting and the trajectories inside the bounded region T bordered by the homoclinic trajectory T reg spiral towards it when t. Hence any such trajectory inside T either has a limit cycle when t which is impossible by Lemma 3. or converges to P which is also impossible. Consequently there exists no homoclinic trajectory at, ) which ends the proof. ext we study the case where is large enough. We have already proved in [7] that for any p > 1, there exists =, p) > see introduction) such that if > there exists no ground state, radial or non-radial. In the radial case we have a more precise result. Proposition 3.6 Let 3 and p > +. Then for any there exists no ground state, but there exist singular solutions u which satisfy lim r ur) = X and lim r ur) = c >. r r Proof. Since >, P is a source by Lemma 3.1, and there exist no periodic orbit neither a homoclinic trajectory at, ) by Lemma 3.3. Thus T reg leaves Q through the semiaxis {, y) : y > } by Lemma.5. As in the proof of Proposition 3.5 the stable trajectory T st at, ) remains in the negatively invariant region Ξ already defined. Then it converges necessarily to P when t. ext we study the case < <. Theorem 3.7 Let 3, p > +. There exist two positive real numbers min and max such that < min max < such that, 1- For < < min there exist ground states ondulating around U when r and positive singular solutions ondulating around U on [, ).

30 Quasilinear elliptic equations with source reaction 3 - For = min and for = max there exist ground states u such that lim r r ur) = c >. 3- For max < < there exists no ground state and there exist singular solutions such that lim r ur) = X or turning around P when r, and lim r ur) = c >. r r Proof. Since is subject to vary, we put it in exponent in the different specific trajectories of the system. For < <, P is a source and the trajectories converging toward this point when t are spiralling. We have three possibilities: i) either T reg leaves Q at some point, y s ), y s >, ii) or T reg has a ω-limit set which is a periodic orbit surrounding P, iii) or T reg converges to, ) when t. If i) holds, then T st remains in the region Ξ := Ξ of Q bordered by T reg and the semi-axis {, y) : y > }. Then, either it converges to P when t, or it admits an unstable from outside) α-limit cycle. We denote by x st, y st ) the first backward intersection of T st with the straight line L. Furthermore T reg intersects L at some point necessarily unique) x reg, y reg). Because of the relative position of T st and T reg, there holds g) = x reg x st >. If ii) holds, then we claim that T st leaves Q at some point x s, ) with x s >. Indeed, if T st Q, it is bounded by Proposition.3. Hence the α-limit set is non-empty. It cannot be, ) since T st T reg and T st cannot converge to P or have a limit cycle around P since it would imply that T st T reg { }. Hence T st intersects L at some point x st, y st ) and T reg intersect the first time L and x reg, y reg). Because T reg lies in the region of Q bordered by T st and the semi-axis {x, ) : x > }, there holds If iii) holds then g) =. g) = x reg x st <. The function g defined for any, ) is continuous. We know that there exists ɛ > such that for any, + ɛ), T reg has a ω-limit cycle surrounding P, hence g) <. If =, T reg leaves Q at some point, y s ) with y s > from Proposition 3.6 and the intersection of T reg with the semi-axis {, y) : y > } is transverse at, y s ). By continuity with respect to the parameter, for any ɛ, ], for ɛ > small enough, T reg leaves Q at some point of the semi-axis. So we are in situation i) and g) >. Then, by continuity of g there exist min and max such that < min max < such that g min ) = g max ) =. If = min or = max, T st = T reg and the trajectory T reg is a homoclinic at, ). For < < min we are in situation ii), which proves 1. For > max we are in situation i) and either T st converges to P or has a α-limit cycle surrounding P when t. Theorems A and A follow from the previous results.

Travelling waves. Chapter 8. 1 Introduction

Travelling waves. Chapter 8. 1 Introduction Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below. 54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and

More information

Half of Final Exam Name: Practice Problems October 28, 2014

Half of Final Exam Name: Practice Problems October 28, 2014 Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A

More information

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12 Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider

More information

A plane autonomous system is a pair of simultaneous first-order differential equations,

A plane autonomous system is a pair of simultaneous first-order differential equations, Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

More information

Nonlinear Autonomous Systems of Differential

Nonlinear Autonomous Systems of Differential Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such

More information

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential

More information

Stability of Dynamical systems

Stability of Dynamical systems Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u

More information

ẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.

ẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D. 4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R

More information

Nonlinear dynamics & chaos BECS

Nonlinear dynamics & chaos BECS Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes

More information

MATH 215/255 Solutions to Additional Practice Problems April dy dt

MATH 215/255 Solutions to Additional Practice Problems April dy dt . For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the

More information

2.10 Saddles, Nodes, Foci and Centers

2.10 Saddles, Nodes, Foci and Centers 2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one

More information

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre

More information

Problem Set Number 5, j/2.036j MIT (Fall 2014)

Problem Set Number 5, j/2.036j MIT (Fall 2014) Problem Set Number 5, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept.,Cambridge, MA 02139) Due Fri., October 24, 2014. October 17, 2014 1 Large µ limit for Liénard system #03 Statement:

More information

STABILITY OF PLANAR NONLINEAR SWITCHED SYSTEMS

STABILITY OF PLANAR NONLINEAR SWITCHED SYSTEMS LABORATOIRE INORMATIQUE, SINAUX ET SYSTÈMES DE SOPHIA ANTIPOLIS UMR 6070 STABILITY O PLANAR NONLINEAR SWITCHED SYSTEMS Ugo Boscain, régoire Charlot Projet TOpModel Rapport de recherche ISRN I3S/RR 2004-07

More information

7 Two-dimensional bifurcations

7 Two-dimensional bifurcations 7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed

More information

Convergence Rate of Nonlinear Switched Systems

Convergence Rate of Nonlinear Switched Systems Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the

More information

7 Planar systems of linear ODE

7 Planar systems of linear ODE 7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

More information

MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/793: Analysis and Control of Nonlinear Systems MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear

More information

Problem set 7 Math 207A, Fall 2011 Solutions

Problem set 7 Math 207A, Fall 2011 Solutions Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase

More information

1 Lyapunov theory of stability

1 Lyapunov theory of stability M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability

More information

Lyapunov Stability Theory

Lyapunov Stability Theory Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous

More information

Now I switch to nonlinear systems. In this chapter the main object of study will be

Now I switch to nonlinear systems. In this chapter the main object of study will be Chapter 4 Stability 4.1 Autonomous systems Now I switch to nonlinear systems. In this chapter the main object of study will be ẋ = f(x), x(t) X R k, f : X R k, (4.1) where f is supposed to be locally Lipschitz

More information

On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I

On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I IOP PUBLISHING Nonlinearity 2 (28) 923 972 NONLINEARITY doi:.88/95-775/2/5/3 On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I S V Gonchenko, L P Shilnikov and D V Turaev

More information

THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL

THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd

More information

Asymptotic Convergence of the Steepest Descent Method for the Exponential Penalty in Linear Programming

Asymptotic Convergence of the Steepest Descent Method for the Exponential Penalty in Linear Programming Journal of Convex Analysis Volume 2 (1995), No.1/2, 145 152 Asymptotic Convergence of the Steepest Descent Method for the Exponential Penalty in Linear Programming R. Cominetti 1 Universidad de Chile,

More information

Part II. Dynamical Systems. Year

Part II. Dynamical Systems. Year Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2

More information

Section 9.3 Phase Plane Portraits (for Planar Systems)

Section 9.3 Phase Plane Portraits (for Planar Systems) Section 9.3 Phase Plane Portraits (for Planar Systems) Key Terms: Equilibrium point of planer system yꞌ = Ay o Equilibrium solution Exponential solutions o Half-line solutions Unstable solution Stable

More information

ENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)

ENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2) ENGI 940 4.06 - Linear Approximation () Page 4. 4.06 Linear Approximation to a System of Non-Linear ODEs () From sections 4.0 and 4.0, the non-linear system dx dy = x = P( x, y), = y = Q( x, y) () with

More information

The Liapunov Method for Determining Stability (DRAFT)

The Liapunov Method for Determining Stability (DRAFT) 44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level

More information

Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains

Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains J. Földes Department of Mathematics, Univerité Libre de Bruxelles 1050 Brussels, Belgium P. Poláčik School

More information

LMI Methods in Optimal and Robust Control

LMI Methods in Optimal and Robust Control LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear

More information

Autonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous.

Autonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous. Autonomous equations Autonomous systems Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous. i f i(x 1, x 2,..., x n ) for i 1,..., n As you

More information

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.

More information

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t

More information

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability

More information

Calculus and Differential Equations II

Calculus and Differential Equations II MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,

More information

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence

More information

Mathematical Economics. Lecture Notes (in extracts)

Mathematical Economics. Lecture Notes (in extracts) Prof. Dr. Frank Werner Faculty of Mathematics Institute of Mathematical Optimization (IMO) http://math.uni-magdeburg.de/ werner/math-ec-new.html Mathematical Economics Lecture Notes (in extracts) Winter

More information

(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);

(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively); STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend

More information

Discretized Fast-Slow Systems near Pitchfork Singularities

Discretized Fast-Slow Systems near Pitchfork Singularities Discretized Fast-Slow Systems near Pitchfork Singularities Luca Arcidiacono, Maximilian Engel, Christian Kuehn arxiv:90.065v [math.ds] 8 Feb 09 February 9, 09 Abstract Motivated by the normal form of a

More information

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1 Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 14, 2019, at 08 30 12 30 Johanneberg Kristian

More information

Nonlinear Control Lecture 5: Stability Analysis II

Nonlinear Control Lecture 5: Stability Analysis II Nonlinear Control Lecture 5: Stability Analysis II Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 5 1/41

More information

Mathematical Modeling I

Mathematical Modeling I Mathematical Modeling I Dr. Zachariah Sinkala Department of Mathematical Sciences Middle Tennessee State University Murfreesboro Tennessee 37132, USA November 5, 2011 1d systems To understand more complex

More information

Stability Implications of Bendixson s Criterion

Stability Implications of Bendixson s Criterion Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 1998 Stability Implications of Bendixson s Criterion C. Connell McCluskey Wilfrid Laurier University,

More information

Prof. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait

Prof. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use

More information

Homoclinic and Heteroclinic Connections. for a Semilinear Parabolic Equation. Marek Fila (Comenius University)

Homoclinic and Heteroclinic Connections. for a Semilinear Parabolic Equation. Marek Fila (Comenius University) Homoclinic and Heteroclinic Connections for a Semilinear Parabolic Equation Marek Fila (Comenius University) with Eiji Yanagida (Tohoku University) Consider the Fujita equation (E) u t =Δu + u p 1 u in

More information

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,

More information

arxiv: v2 [math.ap] 28 Nov 2016

arxiv: v2 [math.ap] 28 Nov 2016 ONE-DIMENSIONAL SAIONARY MEAN-FIELD GAMES WIH LOCAL COUPLING DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES arxiv:1611.8161v [math.ap] 8 Nov 16 Abstract. A standard assumption in mean-field game

More information

2D-Volterra-Lotka Modeling For 2 Species

2D-Volterra-Lotka Modeling For 2 Species Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose

More information

EN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015

EN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions

More information

Chapter III. Stability of Linear Systems

Chapter III. Stability of Linear Systems 1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,

More information

MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation

MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation 1 Bifurcations in Multiple Dimensions When we were considering one-dimensional systems, we saw that subtle changes in parameter

More information

Centre for Mathematics and Its Applications The Australian National University Canberra, ACT 0200 Australia. 1. Introduction

Centre for Mathematics and Its Applications The Australian National University Canberra, ACT 0200 Australia. 1. Introduction ON LOCALLY CONVEX HYPERSURFACES WITH BOUNDARY Neil S. Trudinger Xu-Jia Wang Centre for Mathematics and Its Applications The Australian National University Canberra, ACT 0200 Australia Abstract. In this

More information

Nonlinear Control Lecture 2:Phase Plane Analysis

Nonlinear Control Lecture 2:Phase Plane Analysis Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53

More information

There is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)

There is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth) 82 Introduction Liapunov Functions Besides the Liapunov spectral theorem, there is another basic method of proving stability that is a generalization of the energy method we have seen in the introductory

More information

Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R

Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear

More information

Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling

Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling 1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples

More information

Global unbounded solutions of the Fujita equation in the intermediate range

Global unbounded solutions of the Fujita equation in the intermediate range Global unbounded solutions of the Fujita equation in the intermediate range Peter Poláčik School of Mathematics, University of Minnesota, Minneapolis, MN 55455, USA Eiji Yanagida Department of Mathematics,

More information

The Higgins-Selkov oscillator

The Higgins-Selkov oscillator The Higgins-Selkov oscillator May 14, 2014 Here I analyse the long-time behaviour of the Higgins-Selkov oscillator. The system is ẋ = k 0 k 1 xy 2, (1 ẏ = k 1 xy 2 k 2 y. (2 The unknowns x and y, being

More information

ODE Final exam - Solutions

ODE Final exam - Solutions ODE Final exam - Solutions May 3, 018 1 Computational questions (30 For all the following ODE s with given initial condition, find the expression of the solution as a function of the time variable t You

More information

Introduction to the Phase Plane

Introduction to the Phase Plane Introduction to the Phase Plane June, 06 The Phase Line A single first order differential equation of the form = f(y) () makes no mention of t in the function f. Such a differential equation is called

More information

2 A Model, Harmonic Map, Problem

2 A Model, Harmonic Map, Problem ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or

More information

Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form

Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form d x x 1 = A, where A = dt y y a11 a 12 a 21 a 22 Here the entries of the coefficient matrix

More information

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix

More information

Periodic Sinks and Observable Chaos

Periodic Sinks and Observable Chaos Periodic Sinks and Observable Chaos Systems of Study: Let M = S 1 R. T a,b,l : M M is a three-parameter family of maps defined by where θ S 1, r R. θ 1 = a+θ +Lsin2πθ +r r 1 = br +blsin2πθ Outline of Contents:

More information

Threshold behavior and non-quasiconvergent solutions with localized initial data for bistable reaction-diffusion equations

Threshold behavior and non-quasiconvergent solutions with localized initial data for bistable reaction-diffusion equations Threshold behavior and non-quasiconvergent solutions with localized initial data for bistable reaction-diffusion equations P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455

More information

Hybrid Control and Switched Systems. Lecture #11 Stability of switched system: Arbitrary switching

Hybrid Control and Switched Systems. Lecture #11 Stability of switched system: Arbitrary switching Hybrid Control and Switched Systems Lecture #11 Stability of switched system: Arbitrary switching João P. Hespanha University of California at Santa Barbara Stability under arbitrary switching Instability

More information

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 08, 2018, at 08 30 12 30 Johanneberg Kristian

More information

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient

More information

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form

More information

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point

More information

A CONVEX-CONCAVE ELLIPTIC PROBLEM WITH A PARAMETER ON THE BOUNDARY CONDITION

A CONVEX-CONCAVE ELLIPTIC PROBLEM WITH A PARAMETER ON THE BOUNDARY CONDITION A CONVEX-CONCAVE ELLIPTIC PROBLEM WITH A PARAMETER ON THE BOUNDARY CONDITION JORGE GARCÍA-MELIÁN, JULIO D. ROSSI AND JOSÉ C. SABINA DE LIS Abstract. In this paper we study existence and multiplicity of

More information

THE CRITICAL WAVE SPEED FOR THE FISHER-KOLMOGOROV-PETROWSKII-PISCOUNOV EQUATION WITH CUT-OFF

THE CRITICAL WAVE SPEED FOR THE FISHER-KOLMOGOROV-PETROWSKII-PISCOUNOV EQUATION WITH CUT-OFF THE CRITICAL WAVE SPEED FOR THE FISHER-KOLMOGOROV-PETROWSKII-PISCOUNOV EQUATION WITH CUT-OFF FREDDY DUMORTIER, NIKOLA POPOVIĆ, AND TASSO J. KAPER Abstract. The Fisher-Kolmogorov-Petrowskii-Piscounov (FKPP)

More information

Course 212: Academic Year Section 1: Metric Spaces

Course 212: Academic Year Section 1: Metric Spaces Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........

More information

One Dimensional Dynamical Systems

One Dimensional Dynamical Systems 16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar

More information

Canards at Folded Nodes

Canards at Folded Nodes Canards at Folded Nodes John Guckenheimer and Radu Haiduc Mathematics Department, Ithaca, NY 14853 For Yulij Il yashenko with admiration and affection on the occasion of his 60th birthday March 18, 2003

More information

Practice Problems for Final Exam

Practice Problems for Final Exam Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,

More information

Differential equations, comprehensive exam topics and sample questions

Differential equations, comprehensive exam topics and sample questions Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions

More information

1 2 predators competing for 1 prey

1 2 predators competing for 1 prey 1 2 predators competing for 1 prey I consider here the equations for two predator species competing for 1 prey species The equations of the system are H (t) = rh(1 H K ) a1hp1 1+a a 2HP 2 1T 1H 1 + a 2

More information

Stability of Stochastic Differential Equations

Stability of Stochastic Differential Equations Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010

More information

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables. Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These

More information

Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)

Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x

More information

REUNotes08-ODEs May 30, Chapter Two. Differential Equations

REUNotes08-ODEs May 30, Chapter Two. Differential Equations Chapter Two Differential Equations 4 CHAPTER 2 2.1 LINEAR EQUATIONS We consider ordinary differential equations of the following form dx = f(t, x). dt where x is either a real variable or a vector of several

More information

ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang

ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 11, Number 1, July 004 pp. 189 04 ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS Tian Ma Department of

More information

Example of a Blue Sky Catastrophe

Example of a Blue Sky Catastrophe PUB:[SXG.TEMP]TRANS2913EL.PS 16-OCT-2001 11:08:53.21 SXG Page: 99 (1) Amer. Math. Soc. Transl. (2) Vol. 200, 2000 Example of a Blue Sky Catastrophe Nikolaĭ Gavrilov and Andrey Shilnikov To the memory of

More information

Department of Mathematics IIT Guwahati

Department of Mathematics IIT Guwahati Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,

More information

Perturbation Theory of Dynamical Systems

Perturbation Theory of Dynamical Systems Perturbation Theory of Dynamical Systems arxiv:math/0111178v1 [math.ho] 15 Nov 2001 Nils Berglund Department of Mathematics ETH Zürich 8092 Zürich Switzerland Lecture Notes Summer Semester 2001 Version:

More information

8.1 Bifurcations of Equilibria

8.1 Bifurcations of Equilibria 1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations

More information

In essence, Dynamical Systems is a science which studies differential equations. A differential equation here is the equation

In essence, Dynamical Systems is a science which studies differential equations. A differential equation here is the equation Lecture I In essence, Dynamical Systems is a science which studies differential equations. A differential equation here is the equation ẋ(t) = f(x(t), t) where f is a given function, and x(t) is an unknown

More information

Lecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL)

Lecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL) Lecture 5 Numerical continuation of connecting orbits of iterated maps and ODEs Yu.A. Kuznetsov (Utrecht University, NL) May 26, 2009 1 Contents 1. Point-to-point connections. 2. Continuation of homoclinic

More information

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University. Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ

More information

Math 273 (51) - Final

Math 273 (51) - Final Name: Id #: Math 273 (5) - Final Autumn Quarter 26 Thursday, December 8, 26-6: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all

More information

BIFURCATION FOR NON LINEAR ORDINARY DIFFERENTIAL EQUATIONS WITH SINGULAR PERTURBATION

BIFURCATION FOR NON LINEAR ORDINARY DIFFERENTIAL EQUATIONS WITH SINGULAR PERTURBATION Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 275, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu BIFURCATION FOR NON LINEAR ORDINARY DIFFERENTIAL

More information

SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction

SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.

More information