Asymptotic description of separatrix crossing
near a saddle-center point
http://arxiv.org/abs/math/0105011
O.M.Kiselev, S.G.Glebov,
Institute of Mathematics
of Ufa Sci. Centre of RAS,
Ufa State Petroleum University
We consider a formal asymptotic solution of the main resonance
equation:
|
ieU¢+(|U|2-t)U = 1, 0 < e << 1. |
| (1) |
This equation defines a behaviour of an nonlinear system with
cubic resonance and therefore it is one of major equations
for nonlinear studies. The equation has been written in the
form (1) which means the solution depend on
fast (t) and slow (et) typical scales.
We investigate a saddle-center bifurcation for
the slowly varying equilibriums of this equation and
construct a matching asymptotic solution uniformly as
e®0 before, inside and after the bifurcation layer.
This problem may be considered as a separatrix crossing in
confluent point. The passing through a separatrix of the
second order equations in a general position was considered by
A.V. Timofeev in 1979, A.I. Neishtadt in 1986.
The separatrix crossing for the second order equations in the
confluent point was considered in a preliminary fashion by
R.Haberman in 1979 and D.C.Diminnie
and R.Haberman in 2000 in more detail. The asymptotic
solution crossing the separatrix in the confluent point was
constructed by O.M.Kiselev in 1999,2001 for the
Painlevé-2 equation.
Using our approach one can construct the uniform
asymptotic solution crossing the separatrix in the confluent
point for second order equation in general case.
Algebraic analysis
Let us seek an asymptotic solution of (1) in the
form of the formal asymptotic series:
|
U(t;e) = |
¥ å
n = 0
|
en |
n
U
|
(t). |
| (2) |
After a natural supposition about boundedness of derivative in
the equation (1) we obtain the nonlinear
equation for the main term of the asymptotic expansion:
|
| |
0
U
|
|2 |
0
U
|
- t |
0
U
|
= 1. |
| (3) |
The number of the roots of this algebraic equation depends on
a parameter t. There exist a value of the parameter t
equals to t* = 3(1/2)2/3 so that the equation (3)
has three real roots at t > t*. At t = t* there is one
simple root and one multiple root U* = -(1/2)1/3. At
t < t* the equation (3) has the alone root.
What happened with the asymptotic
solution of the equation (1) when two roots of
the equation (3) coalesce?
Qualitative analysis
To obtain an answer on the precision question one can
consider an autonomous equation with a "frozen"
coefficient T:
|
|
 |
iV'+(|V|^2-T)V=1.
This equation has three equilibrium positions T > t*. There
are U1 < U2 < U3, where U1 is a saddle, U2 and U3
are centers.
|
|
|
| At T = t* the saddle-node
bifurcation takes place and there exist center U3 and
confluent saddle-center point U*. |
|
|
|
 | When T < t* there exists along center U3.
|
Statement of the problem
We will construct the formal asymptotic solution of the
equation (1) in the interval t Î [t*-C,t*+C]
where C = const > 0 uniform on e . We suppose that the
solution in the domain t > t* has the form
|
U(t,e) = |
¥ å
n = 0
|
en |
n
U
|
(t), where |
0
U
|
(t) = U2(t) |
|
.
The qualitative analysis shows that this
asymptotic solution oscillates when t < t*. Our problem is
to study the transition layer between the nonoscillating
asymptotics when t > t* and the oscillating asymptotics when
t < t*.
Numeric evaluations
The numeric evaluations for the special solution
of the equation (1) is given the picture:
Asymptotic analysis
|
|
 |
In the domain
(t-t*)e-4/5 >> 1 the asymptotics has the form:
U(t,e) = ån = 0¥enU(t). Here
U(t) = U2(t) and corrections
U(t) are algebraic functions of
t.
|
|
In the domain |t-t*| << 1 the asymptotics is
defined by four various expansions of different types. First
of them is: |
 |
|
U(t,e) = U* +e2/5 |
¥ å
n = 0
|
e2n/5 |
æ è
|
|
n
a
|
(t)+ ie1/5 |
n
b
|
(t) |
ö ø
|
, |
| (4) |
where t = (t-t*)e-4/5. The leader term
a(t) is a special solution of
the Painlevé-1 equation :
with the given
asymptotics as t®¥:
|
|
0
a
|
(t) = |
å
n ³ 0
|
ant-[((5n-1))/ 2], where a0 = |
1
Ö3
|
, a1 = |
1
24
|
. |
|
In the domain t > -¥ this solution has poles on the
real axis of t. Denote the largest of them by t0.
The asymptotics (4) is valid as
(t-t0)e-1/5 >> 1.
In the neighborhood of t = t0 the coefficients of the
asymptotic expansion depend on one more fast time scale
q = (t-t0)e-1/5. Denote by
|
q0 = q+ |
¥ å
n = 1
|
en/5 |
n
q
|
0
|
, |
|
where q0 = const. Then in the domain
-e-1/5 << q0 << e-1/10 the formal asymptotic
solution has the form :
|
U(t,e) = U*+ |
0
w
|
(q0)+e4/5 |
¥ å
n = 1
|
e(n-1)/5 |
n
w
|
(q0). |
|
The main term of asymptotics w(q0) is
the separatrix solution of the autonomous
equation :
|
i |
0
w
|
¢ +U* |
æ ç
è
|
2| |
0
w
|
|2+ |
0
w
|
2 |
ö ÷
ø
|
+U*2 |
æ ç
è
|
|
0
w
|
* - |
0
w
|
|
ö ÷
ø
|
+| |
0
w
|
|2 |
0
w
|
= 0, |
| (5) |
 |
namely:
w(q0) = [(-2)/((q0-iU*)2)].
In the domain -q0 >> 1 the asymptotic
solution is defined by a sequence of two alternating
asymptotics. Let us call them by "intermediate"
and "separatrix" asymptotics. To obtain the
intermediate asymptotics let us introduce one more slow
variable:
|
|
Tk = qk-1e1/6, k = 1,2,.... |
|
An asymptotic solution in the intermediate domain for not too
large values k << e-1/7 has the form:
U(t,e) = U*+e1/3ån = 0¥ei/30(Ak+ie1/6 Bk). The
leader term satisfies to the equation :
and can be expressed by the Weierstrass
Ã-function:
Ak = -2Ã(Tk;0,g3(k)), g3(k) = [1/56](g3(k-1)+p/2).
 |
Here
g3(0) = [(a4)/56], a4 is the coefficient as
(t-t0)4 in the Laurent expansion of
a(t). The intermediate expansion with the
leader termis valid in the domain between two poles Tk = 0 and Tk = Wk of
theÃ-function:
-e-1/6Tk >> 1, e-2/15(Tk+Wk) >> 1.
|
At the large values of k the intermediate
asymptotics are constructed in the form
U(t,e) = U*+e1/3ån = 0¥en/6(Ak + ie1/6 Bk).
|
The main term satisfies: |
0
A
|
k
|
¢¢+3 |
0
A
|
k
|
2 = lk, where |
|
|
lk(e) = e1/6 |
æ ç
è
|
|
k å
j = 1
|
Wj + |
¥ å
n = 1
|
e(n-1)/30 |
k å
j = 1
|
|
n
x
|
j
|
+ |
ö ÷
ø
|
. |
|
The main term of the asymptotics is:
|
|
0
A
|
k
|
(Tk) = -2Ã(Tk,lk/2,g3(k,e)), |
|
where
g3(k,e) = g3(k)+ån = 1¥en/30g3(k).
The intermediate
expansion with the leader term is valid in the
domain between the poles of the Weierstrass function as
|
-e-1/6Tk >> 1, e-2/15(Tk+Wk) >> 1. |
|
The separatrix expansions are valid in a small neighborhood
of the Weierstrass function poles. Denote:
qk = (Tk+Wk-1/4ån = 1¥en/30xk+)e-1/6, k = 1,2,....
When |qk|e1/6 << 1 the formal
asymptotic solution of equation (1) has the form
:
|
U(t,e) = U*+ |
0
W
|
(qk)+e4/5 |
¥ å
n = 1
|
e(n-1)/30 |
n
W
|
(qk). |
|
|
The leader term of the asymptotics
W(qk) is a separatrix
solution of the autonomous equation (5):
W(qk) = [(-2)/((qk-iU*)2)].
The sequence of the alternating intermediate expansions and
separatrix asymptotics is valid as e-1/6(t*-t) << 1.
|
 |
In the domain (t*-t)e-2/3 >> 1 the
asymptotic solution breaks into oscillation. The
amplitude of the stimulated oscillations in the solution of
(1) oscillates fast. The form of the solution is:
U(t,e) = U(t1,t,e) + eU(t1,t,e) + e2U(t1,t,e),
|
where t1 is
a new fast variable t1 = S(t)/e+f(t). The
main term of the asymptotics U lies on the
curve G(t):
1/2|y|4 -t |y|2 - (y + [`y]) = E(t), and satisfies to the Cauchy problem for the equation
|
iS¢¶t1U +(|U|2 - t)U = 1,
with an initial condition U|t1 = 0 = u0,
such, that
Á(u0) = 0, Â(u0) = miny Î G(t)(Â(y)). The function S(t) is a solution for
the Cauchy problem
|
iS¢ |
ó õ
|
G(t)
|
|
dy
|
|
| __________________ Ö3y3+(2E+t2)y2+2ty+1
|
|
|
= T,S|t = 0 = 0, |
|
where T = const > 0. The function E(t) is the solution of
the transcendental equation :
The phase shift f is defined by initial problems for the
equation:
|
|
f¢
¶E S
|
¶E I = f1 = const, f(t*) = f0. |
|
References
-
Kuzmak G E, Asymptotic Solutions of Nonlinear Differential
Second-Order Equations with Variable Coefficients, Appl.
Math. and Mech., 1959, V.23, N 3, 515-526.
-
Timofeev A.V.. To a question about constancy of an adiabatic
invariant as changing of the type of the motion. ZhETP, 1978,
v.75, n4, pp.1303-1308.
-
Neishtadt A.I., About change of adiabatic invariant as passing
through a separatrix. Fizika plazmy, 1986, v. 12, n8, pp.
992-1001.
-
Haberman R.. Slowly varying jump and transmision phenomena
associated with aljebraic bifurcation problems. SIAM J. Appl.
Math., v.37, n1, pp. 69-106.
-
Bourland F.J., Haberman R., The modulated phase shit of strongly
nonlinear, slowly varying and weakly damped oscillators. SIAM J.
Appl. Math., 1988, v.48, n3, pp.737-748.
-
Diminnie D.C., Haberman R.. Slow passage through a saddle-center
bifurcation. J. Nonlinear Sci., v.10,
pp.197-221(2000).
-
Kiselev O.M. Asymptotic approach for the rigid condition of
appearance of the oscillations in the solution of the
Painlevé-2 equation.
http://xxx.lanl.gov/solv-int/9902007
-
Kiselev O.M., Hard loss of stability in Painlevé-2 equation.
J. of Nonlin. Math. Phys., 2001, v.8, n1, pp.65-95.
-
Kiselev O.M., Glebov S.G., Asymptotic description of nonlinear
resonance.
http://arxiv.org/abs/math/0105011