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|2y = exp æ 
ç 
è 
it
2e
ö 
÷ 
ø 
,   0 < e <<
 
(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 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): 

|

 
|2

 
- t

 
= 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 æ 
ç 
è 
it
2e
ö 
÷ 
ø 
¥ 
å 
n = 0 
en

 
(t),   where 

 
(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 >>
 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 æ 
è 

a 
 
(t)+ ie1/5
b 
 
(t) ö 
ø 
ö 
÷ 
ø 
exp æ 
ç 
è 
it
2e
ö 
÷ 
ø 
,
 
(4)
 
where t = (t-t*)e-4/5. The leader term a(t) is a special solution of the Painlevé-1 equation :
 

a 
 
¢¢ -3
a 
 
2+t = 0, 
 
with the given asymptotics as -¥:

a 
 
(t) =  å 
n ³ 0 
ant-[((5n-1))/ 2],   where    a0
Ö3
a1
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
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*+

 
(q0)+e4/5 ¥ 
å 
n = 1 
e(n-1)/5

 
(q0) ö 
÷ 
ø 
exp æ 
ç 
è 
it
2e
ö 
÷ 
ø 
.
 
The main term of asymptotics w(q0) is the separatrix solution of the autonomous equation : 
i

 
¢ +U* æ 
ç 
è 
2|

 
|2+

 
2 ö 
÷ 
ø 
+U*2  æ 
ç 
è 


 
*

 
ö 
÷ 
ø 
+|

 
|2 

 
= 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 :



 
k  ¢¢+3

 
k  2 = 0
 
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: 

 
k  ¢¢+3

 
k  2 = lk,   where
 
 
lk(e) = e1/6 æ 
ç 
è 
k 
å 
j = 1 
Wj ¥ 
å 
n = 1 
e(n-1)/30 k 
å 
j = 1 


 
j  + ö 
÷ 
ø 
 
The main term of the asymptotics is: 


 
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*+

 
(qk)+e4/5 ¥ 
å 
n = 1 
e(n-1)/30

 
(qk) ö 
÷ 
ø 
exp æ 
ç 
è 
it
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 :
i ó 
õ 
G(t)  u*du = p
 
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