A system ofN harmonic oscillators having a quadratic nonlinear coupling was numerically analysed. Solutions of the nonlinear problem are decomposed in terms of the normal modes of the linear system. A number of initial conditions were considered. We computed the energy in a mode as a function of time and the corresponding average energy. We found that, when initially all the energy of the mode system was in the lowest-frequency mode and the oscillators were at rest, the final time-averaged energy distribution showed that, as we increased the number of oscillators, the first fewlow odd modes received more energy on the average than neighboring even modes. When the energy was initially equipartitioned among the modes and the oscillators were initially at rest, we found that the time-averaged energy distribution of approximately the firstN/4 modes fell off monotonically with increasing frequency, with the other 3N/4 modes keeping essentially all of their energy. Finally, when the energy was initially distributed in a «Planck-like» distribution, the final time-averaged energy distribution was also a «Planck-like» distribution which was quite similar to the initial distribution. In fact, there was very little energy sharing among the modes in this case. When the space-time variation of the actual string was analysed, we found that the «Planck-like» distribution with the oscillators initially at rest was, at least to a good approximation, a «normal mode» of the nonlinear system.
Read full abstract