Asymptotic solution of primary resonance equation
in bifurcation layer
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 investigate the hard mode of stimulating of two phase oscillations for
the equation:
| ei y¢+
|y|2y
= exp |
æ
ç
è |
|
it2
2e |
|
ö
÷
ø |
, 0 < e << 1 |
|
(1) |
The simplest kind of the asymptotic solutions for this equation is the
solutions oscillating with the frequency of the external force. An equation
for the amplitude is:
| ei U¢+
|U|2 U
- tU = 1, where U = yexp |
æ
è |
-it2/(2e) |
ö
ø |
. |
|
(2) |
We investigate a saddle-center bifurcation for the slowly varying equilibriums
of this equation and construct a matching asymptotic solution unformly
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 sepapraprix 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 Painleve-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
If we suppose that the derivative in the equation (2)
is bounded then we obtain a nonlinear algebraic equation for the main term
of asymptotics U(t):
| | |
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* = 3(1/2)2/3 there
is one simple root and one multiple root U* = -(1/2)1/3.
At t < t* the equation (3)
has the alone root.
One can see these roots are slowly varying equilibrium positions for the
equation (2). The question is:
what happened with the asymptotic solution of the equation
(2) when two roots of the equation (3)
coalesces?
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
| y(t,e)
= exp |
æ
ç
è |
|
it2
2e |
|
ö
÷
ø |
|
¥
å
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 (2)
give the picture:
Asypmptotic analysis
 |
In the domain (t-t*)e-4/5
>> 1
the asymptotics has the form:
y(t,e) = exp([(it2)/(2e)])å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:
| y = |
æ
ç
è |
U* +e2/5 |
¥
å
n = 0 |
e2n/5 |
æ
è |
|
n
a
|
(t)+ ie1/5 |
n
b
|
(t) |
ö
ø |
|
ö
÷
ø |
exp |
æ
ç
è |
|
it2
2e |
|
ö
÷
ø |
, |
|
(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 least 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 :
| y(t,e)
= |
æ
ç
è |
U*+ |
0
w
|
(q0)+e4/5 |
¥
å
n = 1 |
e(n-1)/5 |
n
w
|
(q0) |
ö
÷
ø |
exp |
æ
ç
è |
|
it2
2e |
|
ö
÷
ø |
. |
|
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 lternating asymptotics. Let us call
them by ïntermediate" 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:
y(t,e) = (U*+e1/3ån
= 0¥ei/30(Ak+ie1/6Bk))exp([(it2)/(2e)]).
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 term is 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 |
y(t,e) = (U*+e1/3ån
= 0¥en/6(Ak
+ ie1/6Bk))exp([(it2)/(2e)]).
| 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 :
| y = |
æ
ç
è |
U*+ |
0
W
|
(qk)+e4/5 |
¥
å
n = 1 |
e(n-1)/30 |
n
W
|
(qk) |
ö
÷
ø |
exp |
æ
ç
è |
|
it2
2e |
|
ö
÷
ø |
. |
|
| 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 becomes two-phase. The amplitude of
the stimulated oscillations in the solution of (1) oscillates
fast. The form of the solution is:
y = (U(t1,t,e)
+ eU(t1,t,e)
+ e2 U(t1,t,e))exp([(it2)/(e)]),
where t1
is a new fast variable t1 = S(t)/e+f(t).
The main term of the symptotics 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. |
|
-
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