1 Introduction
The problem of the dynamics of a liquid film on the surface a rotating horizontal cylinder in a vertical gravity field has originated from the studies of the behavior of
steam condensate in paper machine dryers [48] and has a number of contemporary industrial applications including spin coating, spin casting, and rotational molding technologies [41].
The considerable interest to this kind of systems is also connected with availability of a variety of experimental configurations [30].
A liquid film is called a rimming flow if the fluid is on the inside of a cylinder, and a coating flow if it is outside. Most of the mathematical studies of rimming and coating flows have been done for the 2-dimensional geometry, where the flow is independent of the axial variable.
Works of Pukhnachev and Moffatt in 1970s on the thin film limit for coating flows led to a great amount of research on the lubrication approximations derived from the governing Navier–Stokes equations, see [9, 7, 46, 36, 6, 41] and references therein. The leading part of the resulting evolution equation is in the class of 4th order nonlinear PDEs, where a rigorous theory of weak solutions has been developing intensively starting from the works of Bernis and Friedman (in the context of rimming flows, see the references and discussions in [14, 12]). We would like to note that the lubrication
approximation has been used and studied extensively for a variety of fluid mechanics models [43, 26, 25, 27, 40, 44], and that, without the lubrication approximation, the exact Stokes equations for rimming flow and emergence of shocks were investigated numerically in [8].
The main mathematical interest to the specifics of thin film models for rimming flows is motivated by a variety of nonlinear phenomena and by several open questions [14, 12, 16, 41].
For the cases where the leading order approximation has various additional differential terms (most important, the smoothing surface tension term), steady state solutions were considered in [36, 6] for certain limiting regimes by numerical and approximation methods. Analytic studies of steady states (in particular, nonexistence theorems) go back to Pukhnachev and were continued in [14], where additionally the regularity of generalized weak solutions was investigated, as well as the short-time and long-time existence. The stability of steady states was intensively studied by approximation, numerical, and linearization methods, see [10, 9, 13] and referenced therein.
In the zero surface tension limit, the question on spectral effects behind unusual instability and
ill-posedness effects raised in the paper [10] led to an intensive spectral analysis of nonselfadjoint differential operators appearing after taking of the Fréchet derivative at a steady state,
see [20, 15, 17, 21, 11] and references therein.
It was shown in [16] that in this case the ill-posedness is also a feature of the original nonlinear equation, and so, is not an artefact of the linearization.
This paper investigates analytically the long-time behavior of solutions
to the rimming-flow model
|
|
|
|
(1.1) |
|
|
|
|
(1.2) |
with -periodic in boundary conditions. Here , , and the circle is systematically identified with the interval .
Namely, we study positive steady states, their stability, and periodic in solutions to
(1.1).
Equation (1.1) can be obtained by means of suitable rescaling of the lubrication approximation model of Benilov, Kopteva, and O’Brien [7], which describes the evolution of Newtonian fluid coating the inner surface of a rotating horizontal cylinder in the thin layer limit and was designed to takes into account the effects of viscosity, gravity, surface tension, and hydrostatic pressure.
The function is, roughly speaking, the nondimensional thickness of the thin film of fluid. The term in (1.1)
corresponds to the motion of the fluid due to the rotation of the cylinder.
The term
describes effects of the gravity. The term accounts for the effect of the surface tension. Finally, the term gives a correction corresponding to differences of the hydrostatic pressure at different parts of the fluid.
Observe that (1.1) is a quasi-linear problem of fourth order, which is degenerate parabolic in the sense that it lacks uniform parabolicity if vanishes at some points.
For small positive values of parameters and the equation (1.1)
has a unique steady state close to the constant solution , which the model has in the case (see Theorem 2.1).
The main points of the present paper are:
-
•
We study analytically the stability properties of the stationary solutions to (1.1), see Theorems 2.3–2.4 and Sections
5.2–5.3. These solutions are obtained in Theorem 2.1 by the application of the implicit function theorem around the explicit steady state corresponding to the unperturbed case .
-
•
We prove rigorously that, for suitable choices of small and , there exists a bifurcation of Poincaré–Andronov–Hopf type, i.e., an onset of a family of periodic solutions branching from the equilibrium . The pairs of values of and , where the bifurcations take place, form an analytic curve in the parameter space, see Theorems 2.4 and 3.11.
-
•
In order to obtain the bifurcation result we employ the theory of multi-parameter Hopf-bifurcations, which is extended to infinite-dimensional dynamical systems in Section 3.
-
•
The analysis of positions of eigenvalues for the linearization of (1.1) is performed with the use of the multi-parameter perturbation theory of eigenvalues, see Appendix A.
To finalize the passage to the analytical formulations and proofs, let us describe in more detail the background of the paper [7], which is our starting point, and the transition from the equations of [7] to the equation (1.1).
In [7] the dynamics of a rimming flow was studied for the model
|
|
|
(1.3) |
which is a simplified version of the nondimensional equation derived by Benilov and O’Brien in [9] from the Navier–Stokes system by means of the lubrication approximation. The model of [9] contains more terms that give for certain ranges of parameters additional information about the physics of the problem, specifically, about the fluid’s inertia. In [7], equation (1.3) was studied as a version of the model of [9] with the inertia terms ignored, and the stability and instability of the related steady states was examined by numerical computations of the spectrum for the linearized equation.
The model of [9, 7] describes the evolution of a thin layer of a Newtonian fluid coating the inner wall of a horizontal cylinder of radius . The cylinder rotates counterclockwise at fixed angular velocity around a horizontal axis. The polar coordinates with the origin being located in the center of the cylinder are used to describe spatial positions.
It is assumed that the fluid has constant density , constant kinematic viscosity , and constant surface tension . We denote as the gravity acceleration.
The nondimensionalization
introduces the function
where is the height of the fluid layer and
In the lubrication approximation limit, is much smaller than . So the function , which solves the equation (1.3), is essentially the nondimensional height of the liquid film.
The
first two terms in the (nondimensional) mass flux
in (1.3) describe the leading order effects due to viscosity and gravity. The term with the parameter corresponds to hydrostatic pressure.
The term containing the parameter describes effects caused by surface tension.
Steady states of (1.3) are solutions to the ordinary differential equation
|
|
|
(1.4) |
with periodic boundary conditions, where is the constant of integration and means . Physically, is the nondimensional mass flux.
In this paper our attention is concentrated on bifurcations of the Poincaré–Andronov–Hopf type
(for the the main definitions of the bifurcation theory we refer to [23, 38, 39, 42],
see also Section 3 below).
To avoid additional technical difficulties connected with the degeneration of the leading term of the equation,
we consider a specific scaling, which will modify equation (1.3) to (1.1).
Namely, let us assume that is positive, but small.
We are interested in steady states of (1.3) for which the modified height of the fluid film is of the same order as the small parameter .
In this case, we find that the term ,
reflecting the gravitational effects, is of the same order of magnitude as .
Then, we choose the parameters and in (1.3) such that: (a)
is of the same order as , (b)
the terms and are of the same order of magnitude,
and (c) the terms due to surface tension and viscosity dominate those due to hydrostatic pressure and gravity. By rescaling equations (1.4) and (1.3) via
,
, as well as
,
, ,
we obtain the ordinary differential equation for the rescaled steady states
|
|
|
and, respectively, we obtain the corresponding evolution equation (1.1), which is the main object of the present paper.
We assume to be of order one, while
are assumed to be small.
We are interested in the branching of periodic solutions at bifurcation points
marking the passage from stable to unstable steady-states.
Such bifurcation points are spectrally critical in the sense that a conjugate pair of nonzero eigenvalues of the linearized operator is situated on the pure imaginary line and crosses this line under specific small variations of . This effect is usually associated with the variety of bifurcations of Poincaré–Andronov–Hopf type that distinguish how the eigenvalues cross (tangentially or nontangentially), positions of other eigenvalues, and various types of degeneracy, e.g., connected with higher multiplicities of eigenvalues.
In the context of multi-parameter bifurcations of finite-dimensional dynamical systems and retarded systems of ODEs we refer to [19, 24, 34, 22] and to the detailed descriptions of various types of the Poincaré–Andronov–Hopf bifurcations in the monograph [39].
The case, where the dynamical system is infinite-dimensional, brings a number of additional effects into the bifurcation theory, in particular, because of the necessity to trace the behavior of infinite number of eigenvalues of linearized operators. The study of bifurcations of the Poincaré–Andronov-Hopf type in this context seems to be initiated for the one-parameter case by Crandall and Rabinowitz [18].
Several types of Poincaré–Andronov-Hopf bifurcations for infinite-dimensional dynamical systems were described in monographs [38, 23, 42]. However, the settings of these monographs do not fit to specifics of the rimming-flow equations (1.3) and (1.1)
by two reasons:
-
•
The models of [9, 7], as well as equation (1.1), include several independent parameters, which makes it natural to consider them from the point of view of the multi-parameter bifurcation theory. Multi-parameter Hopf-bifurcations in the infinite-dimensional phase spaces are not fully understood yet (see the discussion in [22]).
-
•
In comparison with the standard spectral diagram for Poincaré–Andronov-Hopf bifurcations, the linearized operators for the dynamical systems (1.3) and (1.1) have always an additional eigenvalue at zero, which corresponds to the mass conservation. This effect can be handled by the restriction of the dynamical system to the invariant hyperplanes
|
|
|
(1.5) |
where is the -inner product, and is the constant function equal on the whole .
This removes the zero eigenvalue, but brings into (1.1) an additional parameter.
In Section 3, we develop an abstract version of the multi-parameter Hopf-bifurcation theory, which deals with infinite–dimensional dynamical systems and is sufficient for the specifics of the rimming-flow model under consideration. We show that the rimming-flow model (1.1) fits to a specific singular, but relatively simple, case when two bifurcation curves cross each other nontangentially (see Sections 3.4 and 5). This provides then the analytic basis for the computation of the tangent vector of physically relevant branch of the Hopf-bifurcation locus in the the first quadrant of the -parameter plane (see Theorem 2.3).
Notation.
By and we denote general Banach spaces over or . By we denote the Banach space of bounded (linear) operators . In particular, . If , we sometimes consider linear operators , as operators in the space with the domain , see [37]. The corresponding notation is .
By the kernel of a linear operator is denoted, and by
its range.
By we denote the identity operator on .
Inner products in real or complex Hilbert spaces are denoted by .
Moreover, we use for an open ball of radius in . In particular, is an open complex disc. A neighborhood of a point is an open set containing . To distinguish several different neighborhoods of , we often index them with the subscript, i.e., , , ….
By the set of natural numbers is denoted, while .
The notation stands for the standard inner product in , while
|
|
|
is an inner product in the complex Hilbert space . By we denote the standard signum function.
2 Main results of the paper
The parameter is fixed through the rest of the paper.
We start from the rigorous description of a 2-dimensional analytic surface of steady states of
the rimming-flow equation
|
|
|
|
(2.1) |
where and .
This equation is parametrized by in a sufficiently small -neighborhood of . As before, we identify the circle with the interval , i.e., solutions are -periodic in . (Recall that, the physically relevant quadrant of the parameter plane is
.)
Steady-state solutions to (2.1) are the time independent solutions to the ordinary differential equation
|
|
|
(2.2) |
where stands for and denotes the -th power of .
Note that in the case , the constant function is an explicit stationary solution to (2.1) satisfying (2.2).
The following functional spaces and norms are used. We identify functions in the (real) Lebesgue spaces
with -periodic functions in .
The real Hilbert space and its complex version are equipped with the inner product
.
The induced norm is denoted by . With this choice of the inner product,
is an orthonormal basis in .
By
we denote the corresponding real Hilbertian Sobolev spaces; in particular, and coincide as linear spaces and have equivalent norms.
We can represent any by its Fourier series
, where the Fourier coefficients are such that
.
Denoting
|
|
|
one sees that is a closed subspace of .
More generally, for ,
the closed subspaces of , are real (respectively, complex) Hilbert spaces.
Theorem 2.1.
There exist constants such that the following statements hold for all in the neighborhood :
-
(i)
the steady state equation (2.2) admits a unique solution satisfying
;
-
(ii)
the -valued map is real analytic in
in ;
-
(iii)
;
-
(iv)
the first terms of the Maclaurin series of the -valued map are given by the asymptotic formula
as .
The proof is given in Section 4 with the use of the implicit function theorem for analytic maps applied around and .
The quantity corresponds to the conserved total mass in the evolution equation (1.3).
Mathematically, the fact that hyperplanes are invariant sets of (2.1) is easily seen in the following way. We write (2.1) as the following dynamical system in the Hilbert space :
|
|
|
(2.3) |
|
, |
|
(2.4) |
where is an -parameter.
One sees that
.
Hence, if for the initial condition
|
|
|
(2.5) |
then a solution to (2.3)-(2.5) satisfies
for all . Here and below we perceive as a real analytic map from an (open) -neighborhood of to and treat solutions to the initial value problem (2.3)-(2.5) as
strict -solutions (we use here the settings of [42]).
It is not difficult to see from the form of the map that, for is a vicinity of , the Fréchet derivatives of with respect to (w.r.t.) at a steady state has an eigenvalue at .
We need to remove this eigenvalue before we pass to the analysis of the stability and Hopf bifurcations. To this end, the dynamical system (2.3) is considered as a family of dynamical systems in the invariant hyperplanes .
For the purposes of the spectral analysis of Fréchet derivatives we also consider the complexified version of (2.3) in the complex Hilbert space .
Let us take an arbitrary fixed in the neighborhood described in Theorem 2.1 and consider the corresponding steady state .
Then with .
Using the representation
we write the restriction of system (2.3) to the invariant hyperplane as a dynamical system
w.r.t. the state in the complex Hilbert space
in the following way
|
|
|
(2.6) |
Then the identically zero function is a steady state of (2.6) that corresponds to the steady state
of the original equation (2.1).
The complexification of the mapping can be defined by the same formula as . However, for this purpose have to be understood as a holomorphic extension of the real analytic function from Theorem 2.1. Such a holomorphic extension to a certain -neighborhood
of with can be obtained by means of the complex Maclaurin series of .
The complexification of the Fréchet derivative is given by the operator
|
|
|
(2.7) |
defined for by , where
|
|
|
(2.8) |
We call the linearized operator (corresponding to the 2-parameter ). To describe some of its spectral properties, the following notation and definitions of spectral theory are used (see [37, 47]).
By we denote a complex open disc of radius with the center at , and by the polydisc that is the direct product of two discs .
For a linear operator in a complex Banach space , we denote by its domain (of definition), by its spectrum, and by
its resolvent set. An eigenvalue of is called isolated if there exists a complex open neighborhood of such that . An eigenvalue is called simple if it has algebraic multiplicity .
Roughly speaking, the open left half-plane
in the complex plane of the spectral parameter is associated with the case of stability (see Theorem 2.4).
The following theorem describes the behavior under variation of near of the pair of eigenvalues of the linearized operator which are closest to the purely imaginary axis .
Theorem 2.3.
For certain small enough , there exist two holomorphic functions
and an increasing real analytic function
,
with the following properties:
-
(i)
For all , the part of the spectrum of the linearized operator
consists of exactly two simple isolated eigenvalues such that
|
, , and ; |
|
moreover, the operator is a generator of quasi-bounded holomorphic semigroup (see Remark 3.3).
-
(ii)
For ,
|
|
|
(2.9) |
|
is equivalent to , |
|
(2.10) |
|
|
|
(2.11) |
-
(iii)
The leading terms of the Maclaurin series of are given by
|
|
|
(2.12) |
This theorem is proved in Section 5.2.
Controlling the rest of eigenvalues for small enough ,
we show that
the graph of the function , ,
separates in the first quadrant the regions of instability and exponential stability for the steady state of the reduced system (2.6).
That is, it separates the regions of instability and exponential stability of the steady state of (2.1) if the dynamics considered in the corresponding hyperplane .
(Note that, because of the invariance of the hyperplanes , the steady states obviously cannot be asymptotically stable in the sense of the original equation (2.1)).
It will be shown in Theorem 2.4 that
the graph of the function is essentially a parametric bifurcation curve
from the point of view of multi-parameter Hopf-bifurcations (see Definition 3.8 and Theorem 3.7 below; for the basics of the theory in finite-dimensional settings see [39]).
The rigorous description of multi-parameter Hopf bifurcations is given in Section 3.2, where we adapt to infinite-dimensional dynamical systems the approach suggested by Diekmann and Korvasová [22] for 2-parameter Hopf bifurcations of delay differential equations.
Namely, if the parameter moves along a -path , , that crosses a Hopf-bifurcation curve nontangentially at a certain point , then 1-parameter Hopf bifurcation takes place at this point w.r.t. the new one-dimensional parameter . That is, a family of periodic nonstationary solutions branches from the curve of steady states at .
These stability and bifurcation results are summarized in the next theorem.
Theorem 2.4.
Assume that and the function are as in Theorem 2.3. Let us denote
and
Then the following statements hold:
-
(i)
If , then is an exponentially stable steady state of the restriction of the dynamical system (2.1) to the invariant hyperplane .
-
(ii)
If , then the steady state of the restriction of (2.1) to the invariant hyperplane is not stable in the sense of Lyapunov.
-
(iii)
Let be a simple -path that crosses
the graph nontangentially at certain point
passing from the region
to the region .
Then, for the reduced system (2.6),
|
a Hopf bifurcation takes place at on the path |
|
(2.13) |
in the sense of Definition 3.6.
In particular, for arbitrary there exists such that, for ,
|
|
|
(2.14) |
This theorem is proved in Section 5.3.
3 Multi-parameter Hopf bifurcations
This section considers multi-parameter Hopf bifurcations in the abstract infinite-dimensional settings and provides the theoretical background for the proofs of the results of Section 2.
Bifurcation theory distinguishes between several types of criticality effects. These include, in particular:
-
(B1)
spectral (or linear) criticality, where the spectrum of an operator linearized around an equilibrium corresponds to the threshold between stability and instability;
-
(B2)
dynamic criticality, i.e., the situation when the actual change in the behavior of the dynamical system happens under a modification of a parameter through its specific critical value.
In the finite-dimensional context of Arnold’s lectures [3], situations of type (B2) fit the notion of a bifurcation. Situations of type (B1) usually play the role of necessary spectral conditions for specific bifurcations of type (B2).
We are interested in Hopf bifurcation
(Andronov–Hopf bifurcation in the terminology of [4, 39])
which has a variety of definitions depending on the dimension of the state space and the dimension of the parameter space
(see, e.g., [3, 4, 39, 23]
for the case and ,
[18, 38, 23, 42] for the case and ,
and the discussion of the case and in [38, 22]).
The corresponding change of the long-time dynamics can be described
as a branching of a family of nontrivial periodic solutions from a curve of equilibria
at the point where a conjugate pair of eigenvalues of a linearized operator crosses the stability threshold .
This phenomenon was observed for -dimensional systems already by Poincaré and studied by Andronov (see [3, 4, 39]).
The one-parameter Hopf bifurcation [33] describes an important particular case of such branching where an additional nondynamical condition is satisfied (Hopf condition); namely, the conjugate pair of eigenvalues crosses with the nonvanishing real part of their derivative. Roughly speaking, the Hopf bifurcation can be seen as a pitchfork bifurcation of a curve of cycles from a curve of equilibria.
In this section, we adapt the definition of multi-parameter Hopf bifurcation suggested for delayed ODEs by Diekmann and Korvasova [22] to the case of a dynamical system in an infinite-dimensional state space in such a way that it fits the specifics of the rimming–flow equation (1.1). This definition, roughly speaking, reduces the multi-parameter Hopf bifurcation to the one-parameter case considering paths , , in the parameter space that cross
-dimensional bifurcation surfaces
in at certain points which are critical in the sense of (B1), i.e., in the sense of the spectrum of the linearized operator (we call such points spectrally critical or -critical).
In these settings, plays the role of the 1-dimensional parameter of the one-parameter Hopf-bifurcation theorems in the style of [18, 38, 42].
Consider an abstract evolution equation
|
|
|
(3.1) |
dependent on an -parameter ,
where the set is open in , and is an open subset of a real Banach space .
It is assumed that a real Banach space and the map satisfy the following hypotheses:
|
|
|
(3.2) |
|
|
|
(3.3) |
The derivative is understood in the sense of the convergence w.r.t. the norm of the Banach space .
The complexification of a real Banach space is defined as a formal sum (in particular, is the complexification of the parameter space).
Then, for , the real part , the imaginary part , and the complex conjugation of any vector are well-defined. If is an inner product space, the real inner product of is extended in a natural way to the complex inner product .
Similarly, a real (linear) operator is extended in a natural way to a complex operator .
Note that, if is an eigenvalue of corresponding to an eigenvector (i.e., and ), then the complex conjugate also belongs to the set of eigenvalues of and
.
Denoting by
the Fréchet derivative of w.r.t. at , one obtains a real analytic map
acting from
to .
For each , we define the operator
|
|
|
(3.4) |
of .
In the sequel, the complexified -derivative is considered also as an operator in the space with domain , which corresponds to the notation
.
3.1 Hypersurfaces of spectrally critical -parameters
We define and study in this subsection the family of spectrally critical -parameters using the theory of multi-parameter perturbations of eigenvalues. The necessary results on multi-parameter perturbations are collected in Appendix A. We prove in this subsection Proposition 3.4, which is crucial for the extension of the multi-parameter Hopf-bifurcation theory to the case of infinite-dimensional state space . Proposition 3.4 shows that the combination of the spectral criticality with the sectoriality properties of is robust under suitable multi-parameter perturbations of . This later allows us to define rigorously Hopf-bifurcation manifolds via Proposition 3.5 and Theorem 3.7.
A point is said to be an equilibrium point if is an equilibrium of
, i.e., if .
Accordingly,
is the set of all equilibrium points of the evolution equation .
Let be an equilibrium point such that , i.e., such that
is boundedly invertible as a densely defined operator in .
Then
the analytic implicit function theorem (see, e.g., [45, Appendix B]) implies that there exists an open neighborhood of such that all equilibria of with
have the form , where is a certain real analytic map satisfying .
So provides a local analytic parametrization of the set of all equilibrium points near . In other words, the graph
of , ,
is a local equilibrium manifold near .
Taking into account the stability of the property (see Remarks A.3–A.4 in Appendix A), one sees that the set of all equilibrium points where is invertible
|
|
|
(3.5) |
is an -dimensional real analytic manifold
in .
Definition 3.1 (cf. [38, 42]).
An equilibrium point of the evolution equation is called
-
(i)
spectrally exponentially stable if the spectrum of the complexified Fréchet derivative of is a subset of the halfplane
for a certain ;
-
(ii)
spectrally unstable if .
Definition 3.2 (cf. [39, 38, 23, 42]).
An equilibrium point is said to be -critical (or spectrally critical in the sense of Poincaré–Andronov) if the following two conditions are satisfied:
-
(i)
lies in the closed left half-plane ;
-
(ii)
the pure imaginary part of the spectrum consists of a pair of complex conjugate simple isolated eigenvalues, which are denoted by in such a way that .
In this case, is called a -critical -parameter (associated with
).
Let be a -critical equilibrium point. Using Definition 3.2 (ii), we define eigenvalues of and infer that . Consequently,
|
|
|
and there exists a local analytic parametrization of all equilibrium points in a certain neighborhood of . The restriction of the -derivative to
the local equilibrium manifold can be written as a real analytic function of , which yields after the complexification
the -valued function depending only on .
It follows from Proposition A.1 and Remark A.4 that
there exist an open
neighborhood of and functions such that
|
|
are holomorphic in , , |
|
(3.6) |
|
|
and are simple isolated eigenvalues of
for all , |
|
(3.7) |
where is an -neighborhood of .
Without loss of generality, the neighborhood can be chosen in such a way that
the Taylor series of
converges in .
Here and below the summations are done w.r.t. the multi-parameter , where
.
In the -neighborhood of , we define the function
|
|
|
(3.8) |
which is real analytic since in .
The condition that generates a quasi-bounded holomorphic semigroup (in the terminology of Kato [37]) is one of the conditions used in the infinite-dimensional stability theory [18, 38, 31, 42] to ensure that, except the two eigenvalues , the spectrum of stays far away from the pure imaginary line for close to .
We use this condition in the multi-parameter settings and adapt for this purpose the perturbation results of [32, 37].
Under -sectorial operators we understand generators of quasi-bounded holomorphic semigroups (essentially, they are sectorial operators in the sense of Henry [31]).
That is, for angles , we consider sectors of the form
|
|
|
where the complex argument of is defined as a unique such that for a certain . Let . Slightly modifying the notation of [37], let us denote by the family of densely defined closed operators in a complex Banach space such that the resolvent set contains the sector and, for any smaller sector
with , the uniform estimate
holds with a constant depending on .
An operator belongs to the family if with a certain and certain . Finally, we introduce the notation
|
|
|
(3.9) |
and we say that an operator is -sectorial if .
The next proposition essentially says that, if is an -sectorial operator in , then the set of -critical -parameters coincides locally with the zero-locus of .
Proposition 3.4.
Let be a -critical point and let . Then:
-
(i)
is -sectorial for in a certain open -neighborhood of .
-
(ii)
There exists an open -neighborhood of such that, for all ,
the equilibrium point is -critical if and only if .
Proof.
Since is -critical and is -sectorial, there exists such that the two simple eigenvalues are the only points of in the half-plane .
Using Remark A.4 and Proposition A.2 one can choose the neighborhoods and with the desired properties. This completes the proof.
∎
Definition 3.3.
Let be an equilibrium point. We say that is an -regular equilibrium point with an -regular -parameter if is a regular point of , i.e., if , where is the gradient at of the function introduced in (3.8).
If , the equilibrium point and the associated -parameter are said to be -singular.
Recall that a hypersurface in an -dimensional manifold is by definition an -dimensional submanifold.
Proposition 3.5.
Let be the set of all -regular -critical equilibrium points
such that is -sectorial.
Then is an analytic -dimensional manifold in .
Proof.
Let be an -regular equilibrium point.
Locally in the neighborhood from Proposition 3.4, operators are -sectorial and the set of
-critical -parameters coincides with .
Applying the implicit function theorem at -regular -parameters , we see that the set
|
|
|
is an analytic hypersurface in the parameter space .
Globally, this implies the statement of the proposition.
∎
3.2 Hopf bifurcations on paths
We need a formalization of the notion of continuous families of cycles (confer Hopf-bifurcation theorems in [38, 42]).
A cycle of the dynamical system is its periodic solution with the corresponding period defined by assuming that periodic solutions produced by phase-shifts are considered to be identical.
It is convenient to consider equilibria as degenerate cycles of period .
A simple path is an injective continuous function defined on a nondegenerate interval , where is a certain Banach space.
We use also a particular type of parametrized closed curves without self-intersections, which we call simple loops. A simple loop is understood in this paper as an injective continuous map ,
where as before is the unit circle .
A simple loop is called an -parametrization of a cycle of period if
for all
with a certain .
An -parametrization of a degenerate cycle
(i.e., an -parametrization of an equilibrium point ) is defined as the constant function , .
Consider a function . One says that
|
|
|
is a continuous family of cycles of passing at through an equilibrium point if the three following conditions are fulfilled:
-
(C0)
is an -parametrization of a cycle of ;
-
(C1)
and
is the degenerate cycle with ;
-
(C2)
and defined by is a homotopy.
Definition 3.4.
Let be a simple path such that, for a certain , the value is a -critical -parameter associated with an equilibrium of .
We say that a Poincaré–Andronov bifurcation takes place at on the path if there exists
a subinterval with having the following properties:
-
(i)
for each , the state is a spectrally exponentially stable equilibrium of ;
-
(ii)
for each , the state is a spectrally unstable equilibrium of ;
-
(iii)
there exist a continuous function and
a continuous family of cycles
passing at through
so that the cycle is nondegenerate for .
A Hopf bifurcation w.r.t. a 1-dimensional parameter is a special type of Poincaré–Andronov bifurcation with an additional condition on the behavior of the perturbed eigenvalues , which can be adapted to bifurcations on paths in the following way.
Definition 3.5 (Hopf transversality condition).
Let be a simple path such that is a -critical point for a certain . Let be the two analytic functions defined in a vicinity of by (3.6)–(3.7) such that and are two simple isolated eigenvalues of with the property that belong to . We say that, on the path , the Hopf transversality condition is satisfied at if
|
|
|
(3.10) |
Under certain assumptions in the 1-parameter case,
the Hopf transversality condition implies the existence of a bifurcation of Poincaré–Andronov–Hopf type, as well as a variety of additional properties of the branching curve of cycles (see the Hopf bifurcation theorems in [38, 23, 42]).
In Definition 3.6 below, we adapt these properties to the -parameter case combining them with the approach of [22].
Slightly modifying the notation of [38], we denote by the Banach space of the functions that are Lipschitz continuous on with values in a Banach space .
The Banach space consists of -functions such that the derivative belongs to (the norms in are introduced in the standard way, see [38]).
A continuous family of cycles passing at through an equilibrium point is said to be a -family of cycles passing through at
if the following conditions are fulfilled additionally to (C0)–(C2):
-
(C3)
there exists such that
is equal to the period of for all ;
-
(C4)
for all ;
-
(C5)
the mapping
defines a -path in the Banach space .
Definition 3.6 (Hopf bifurcations on paths).
Assume that on
a simple -path a Poincaré–Andronov bifurcation takes place at . Assume also that on the Hopf transversality condition (3.10) holds at and that are the eigenvalues of defined as in Definitions 3.2 and 3.5.
Then we say that
a Hopf bifurcation takes place at on if for any interval containing there exist
a -function and
a -family of cycles
of
passing at through such that for small enough the following statements hold true:
-
(iv)
the ‘period function’ defined by (C3) satisfies ;
-
(v)
, , and
for all ;
-
(vi)
for , the existence of a nondegenerate cycle of the dynamical system with
an -parametrization
|
|
|
and with the period
implies
and for a certain (i.e., implies the uniqueness of the branching curve of cycles).
Here is the constant function equal to for all .
Definition 3.7.
If a Hopf bifurcation takes place at on a certain path, we say that is a Hopf–bifurcation point and that is the associated Hopf-bifurcation -parameter.
A more general notion of a Poincaré–Andronov–bifurcation point can be defined similarly. Obviously, every
Hopf–bifurcation point is also a Poincaré–Andronov–bifurcation point. The converse generally is not true.
3.3 Manifolds of Hopf-bifurcation points
The aim of this subsection is to formulate rigorously the notion of Hopf-bifurcation manifold (Hopf-bifurcation curve in the case ) and to prove the corresponding multi-parameter Hopf-bifurcation theorem, which essentially says the set of the -regular -critical points satisfying an additional -sectoriality assumption is a Hopf-bifurcation manifold. A particular feature of this multi-parameter bifurcation theorem is that the Hopf transversality condition is essentially replaced by the -regularity (see Remark 3.6).
A -path is said to be regular if for all .
Definition 3.8.
An -dimensional -manifold in is called a Hopf-bifurcation manifold if for every point there exist its -neighborhood and a -hypersurface in with the following properties:
-
(i)
is the set of all
Hopf-bifurcation points in .
-
(ii)
There exists a -map such that
.
-
(iii)
On every simple regular -path passing through nontangentially to the hypersurface , a Hopf bifurcation takes place at after a possible change of orientation of .
If these conditions are fulfilled, we say that is a local (parametric Hopf-) bifurcation hypersurface that parametrizes the local Hopf-bifurcation manifold
. In the case , we also call a connected -dimensional Hopf-bifurcation manifold a Hopf-bifurcation curve and call a connected local parametric Hopf-bifurcation hypersurface a local parametric Hopf-bifurcation curve.
Theorem 3.7 (multi-parameter Hopf-bifurcation theorem).
Assume (3.2) and (3.3). Let be the set of all -regular -critical equilibrium points
such that is -sectorial. Then:
-
(i)
is a Hopf-bifurcation manifold.
-
(ii)
Assume that a -critical point is -singular (i.e., ). Then is not a Hopf-bifurcation point.
-
(iii)
is the set of all Hopf-bifurcation points such that
is -sectorial.
3.3.1 Proof of Theorem 3.7.
Step 1. Proof of statement (ii). If is an -singular -critical point, then Remark 3.6 implies that is not a Hopf-bifurcation point.
Step 2. Let be an -regular -critical point such that
is -sectorial. Then it follows from Section 3.1 that there exist an -neighborhood of and an analytic map such that and the following properties hold:
-
(a)
the manifold is the set of all equilibrium points in ;
-
(b)
the Taylor series at of the real analytic function converges in ;
-
(c)
the set is a connected analytic hypersurface in and, for every , the operator
is -sectorial;
-
(d)
the -dimensional analytic manifold coincides with
and contains .
Remark 3.6 and statement (ii), which is already proved in Step 1, imply that satisfies condition (ii) of Definition 3.8 and that contains all Hopf bifurcation points in .
Step 3. We now verify property (iii) of Definition 3.8. This property implies that
is exactly the set of all Hopf-bifurcation points in , and so, implies statements (i) and (iii) of Theorem 3.7.
Property (iii) of Definition 3.8 follows from Proposition 3.9 given below, which can be proved essentially by the reduction to the one-parameter arguments of Kielhöfer [38, Section I.8].
Proposition 3.9.
Let be a
-critical point such that is -sectorial.
Let , and let be a simple -path with such that the following transversality condition holds true
|
|
|
(3.12) |
Then a Hopf bifurcation takes place at on the path .
Proof of Proposition 3.9..
Without loss of generality, we can assume that for all and consider the reduced dynamical system
|
|
|
(3.13) |
depending only on the -dimensional parameter . The complexification of the Fréchet derivative is equal to
|
|
|
In particular, is an equilibrium point of (3.13) and .
Since is -critical, consists of two isolated simple eigenvalues of , where is a certain number.
There exists a bounded domain with a positively oriented smooth Jordan curve as its boundary satisfying and .
We denote by the Riesz projection
corresponding to .
The closed subspaces and are the images of the projection and of the complementary projection , respectively. By the spectrum decomposition theorem [37, Theorem III.6.17] each of the subspaces and is an invariant subspace of , and the direct sum decomposition
(not necessarily orthogonal)
takes place. These two invariant subspaces are also invariant for the quasi-bounded holomorphic semigroup (see Remark 3.3). In particular, the restriction of to the invariant subspace is an operator in the subspace .
The following statement holds true:
|
is a self-homeomorphism of . |
|
(3.14) |
(compare [42, Proposition 4.4.8]). Indeed,
since and since is -critical, condition (ii) of Definition 3.2
implies that there exists a constant such that
|
|
|
(3.15) |
Besides,
is a generator of a quasi-bounded holomorphic semigroup in . Combining the spectral mapping theorem for quasi-bounded holomorphic semigroups [29, Corollary IV.3.12] with (3.15), we see that .
This completes the proof of (3.14).
Now, the arguments of the proof of [38, Theorem I.8.2] on Hopf bifurcations w.r.t. a -dimensional parameter can be applied to
(see Remark 3.8).
∎
This completes the proofs of Proposition 3.9 and of Theorem 3.7.
3.4 The two-parameter case and bifurcation curves near -singular points
We have shown that the set of -critical points is described locally by the equation . In the two-parameter case , the branching of the curves of the zero-locus for can be studied with the use of the Newton polygon and Puiseux series, see [49, 38].
Motivated by the rimming-flow equation (2.1), we consider a relatively simple case of an -singular point with a nonzero Hessian determinant of . We describe for this case the parametric Hopf-bifurcation curves contained in the local zero-locus of .
To this end, let be an -singular -critical 2-parameter associated with an equilibrium . We denote by
|
|
|
(3.16) |
the Hessian matrix of at and assume that its determinant is nonzero. In this section,
are the Taylor coefficients of at (see (3.8)) and .
Note that the case is trivial from the point of view of Hopf bifurcations since is an isolated -critical point and preserves its sign in a sufficiently small punctured neighborhood of . Hence, is not a Poincaré–Andronov bifurcation point, and so, is not a Hopf–bifurcation point.
The following lemma, describing the local zero-locus of in the case , can be obtained by a combination of several formulas of [49, Section I.2.7] which consider one-by-one the corresponding cases of the Newton diagram. By we denote the signum function.
Lemma 3.10 ([49]).
Let be a real analytic function with
|
, and
,
|
|
where
is the Hessian matrix of .
Then:
-
(i)
There exists a neighborhood of so that the set
consists of the union of two simple regular real analytic curves having a nontangential intersection at as their only intersection point.
-
(ii)
Regular parametrization of can be chosen such that , , and the tangential vectors at are
|
|
|
(3.17) |
|
|
|
(3.18) |
For the reader’s convenience, an elementary proof of Lemma 3.10 without the use of Newton diagrams is given in Appendix B.
Theorem 3.11.
In the settings of Subsection 3.1, assume that and that the operator is -sectorial at an -singular -critical point with .
Then there exists a neighborhood of such that:
-
(i)
for all and where the curves and their parametrization are as in Lemma 3.10. In particular, (3.17)–(3.18) are satisfied.
-
(ii)
The set is a -dimensional Hopf-bifurcation manifold that contains all Hopf-bifurcation points in and consists of four simple disjoint Hopf-bifurcation curves.
-
(iii)
The corresponding set of Hopf-bifurcation 2-parameters is a union of four simple disjoint regular real analytic curves.
Proof.
It follows from Proposition 3.4 and Lemma 3.10 that
can be chosen such that statement (i) holds.
By Theorem 3.7, all local Hopf-bifurcation 2-parameters near belong to and are characterized in this set by the property that the gradient is nonzero. In particular, is not a Hopf-bifurcation 2-parameter.
Since
as ,
we see from (3.17)–(3.18) that on the curves the following asymptotic formulae hold:
|
|
|
|
|
|
It follows from that for with a certain . This means that, after a possible shrinking of ,
|
|
|
and the set described by this equality is the set of all local Hopf-bifurcation parameters.
This proves statement (iii) and, in turn implies statement (ii).
∎
Appendix B The Puiseux series proof of Lemma 3.10
For the convenience of the reader, a short proof of Lemma 3.10 is given in this section.
It is based on the Weierstrass preparation theorem and [47, Theorem 12.2] on the Puiseux series for roots of a complex polynomial depending on a parameter. In comparison with [49, Section I.2.7] this proof does not require the consideration of several related Newton diagrams.
Without loss of generality, we assume .
In a certain -neighborhood of , we consider the function
with the real Taylor coefficients . That is,
is a complex analytic extension of the real analytic function .
Under the assumptions of the lemma, , , and the Hessian determinant is negative.
In particular, .
Step 1. Let us consider the solutions of in the particular case . In this case, the Weierstrass preparation theorem and [47, Theorem 12.2] imply that small enough satisfying have the form , where two Puiseux series with certain coefficients are convergent in a vicinity of . Moreover, the multi-set of the zeroes of the function in a vicinity of takes the multiplicity of zeros into account and the total multiplicity of such zeroes for small enough is constant, and so equals 2, as in the case .
Step 2. Let us prove statement (i) of the lemma. We consider only the case . Every other case, where , can be reduced to this case by a constant linear transformation of the parameter space . Statement (i) is invariant under such transformations.
It follows from that and for
small enough nonzero and .
Hence, in every -neighborhood of , has a solution in each of the four quadrants . Comparing this with the statement of Step 1, one sees that
the two Puiseux series are real locally near for real and real .
This, in turn, implies that
-
(a)
all the Puiseux coefficients are real and
for all odd ;
-
(b)
can be chosen such that and .
Note that statement (a) follows from the local conservation of the total multiplicity of zeroes of . Indeed, if (a) is not fulfilled, the lowest order nonreal coefficient or the lowest order nonzero coefficient with odd index produces a locally nonreal for small branch of the zero-locus of , which violate the local conservation of total multiplicity.
In other words, are real analytical functions of in an -neighborhood of . In order to prove the nontangentiality in the statement (i) of the lemma, it is enough to show that . One obtains immediately from the assumptions and by plugging into .
Step 3. Let us prove statement (ii) of the lemma.
Consider first the case . Steps 1 and 2 imply that the local zero locus of can be represented as the union of the graphs of two real analytic functions with distinct first coefficients .
Plugging the Taylor series of into ,
one gets and, in turn,
This implies (ii) in the case .
The case can be considered similarly. The case can be reduced to , by a linear transformation. This completes the proof.