  We present the applications of nonlinear local harmonic analysis methods to
the modelling of beam-beam interaction. Our approach is based on methods
provided the possibility to work with dynamical beam localization in phase
space. The consideration of Fokker-Planck or Vlasov-Maxwell models is based on
a number of anzatzes, which reduce initial problems to a number of dynamical
systems (with constraints) and on variational-wavelet approach to
polynomial/rational approximations for reduced nonlinear dynamics. We calculate
contribution to full dynamics (partition function) from all underlying
subscales via nonlinear eigenmodes decomposition.
