# Multiple-scale analysis

(Redirected from Method of multiple scales)

In mathematics and physics, multiple-scale analysis (also called the method of multiple scales) comprises techniques used to construct uniformly valid approximations to the solutions of perturbation problems, both for small as well as large values of the independent variables. This is done by introducing fast-scale and slow-scale variables for an independent variable, and subsequently treating these variables, fast and slow, as if they are independent. In the solution process of the perturbation problem thereafter, the resulting additional freedom – introduced by the new independent variables – is used to remove (unwanted) secular terms. The latter puts constraints on the approximate solution, which are called solvability conditions. Mathematics research from about the 1980s proposes that coordinate transforms and invariant manifolds provide much sounder support for multiscale modelling (for example, see center manifold and slow manifold).

## Example: undamped Duffing equation

### Differential equation and energy conservation

As an example for the method of multiple-scale analysis, consider the undamped and unforced Duffing equation:[1]

$\frac{d^2 y}{d t^2} + y + \varepsilon y^3 = 0,$   $y(0)=1, \qquad \frac{dy}{dt}(0)=0,$

which is a second-order ordinary differential equation describing a nonlinear oscillator. A solution y(t) is sought for small values of the (positive) nonlinearity parameter 0 < ε ≪ 1. The undamped Duffing equation is known to be a Hamiltonian system:

$\frac{dp}{dt}=-\frac{\partial H}{\partial q}, \qquad \frac{dq}{dt}=+\frac{\partial H}{\partial p}, \quad \text{ with } \quad H = \tfrac12 p^2 + \tfrac12 q^2 + \tfrac14 \varepsilon q^4,$

with q = y(t) and p = dy/dt. Consequently, the Hamiltonian H(pq) is a conserved quantity, a constant, equal to H = ½ + ¼ ε for the given initial conditions. This implies that both y and dy/dt have to be bounded:

$\left| y(t) \right| \le \sqrt{1 + \tfrac12 \varepsilon} \quad \text{ and } \quad \left| \frac{dy}{dt} \right| \le \sqrt{1 + \tfrac12 \varepsilon} \qquad \text{ for all } t.$

### Straightforward perturbation-series solution

A regular perturbation-series approach to the problem gives the result:

$y(t) = \cos(t) + \varepsilon \left[ \tfrac{1}{32} \cos(3t) - \tfrac{1}{32} \cos(t) - \underbrace{\tfrac38\, t\, \sin(t)}_\text{secular} \right] + \mathcal{O}(\varepsilon^2).$

The last term between the square braces is secular: it grows without bound for large |t|, making the perturbation solution valid for only small values of the time t.

### Method of multiple scales

To construct a globally valid solution, the method of multiple-scale analysis is used. Introduce the slow scale t1:

$t_1 = \varepsilon t\,$

and assume the solution y(t) is a perturbation-series solution dependent both on t and t1, treated as:

$y(t) = Y_0(t,t_1) + \varepsilon Y_1(t,t_1) + \cdots.$

So:

\begin{align} \frac{dy}{dt} &= \left( \frac{\partial Y_0}{\partial t} + \frac{dt_1}{dt} \frac{\partial Y_0}{\partial t_1} \right) + \varepsilon \left( \frac{\partial Y_1}{\partial t} + \frac{dt_1}{dt} \frac{\partial Y_1}{\partial t_1} \right) + \cdots \\ &= \frac{\partial Y_0}{\partial t} + \varepsilon \left( \frac{\partial Y_0}{\partial t_1} + \frac{\partial Y_1}{\partial t} \right) + \mathcal{O}(\varepsilon^2), \end{align}

using dt1/dt = ε. Similarly:

$\frac{d^2 y}{d t^2} = \frac{\partial^2 Y_0}{\partial t^2} + \varepsilon \left( 2 \frac{\partial^2 Y_0}{\partial t\, \partial t_1} + \frac{\partial^2 Y_1}{\partial t^2} \right) + \mathcal{O}(\varepsilon^2).$

Then the zeroth- and first-order problems of the multiple-scales perturbation series for the Duffing equation become:

\begin{align} \frac{\partial^2 Y_0}{\partial t^2} + Y_0 &= 0, \\ \frac{\partial^2 Y_1}{\partial t^2} + Y_1 &= - Y_0^3 - 2\, \frac{\partial^2 Y_0}{\partial t\, \partial t_1}. \end{align}

### Solution

The zeroth-order problem has the general solution:

$Y_0(t,t_1) = A(t_1)\, e^{+it} + A^\ast(t_1)\, e^{-it},$

with A(t1) a complex-valued amplitude to the zeroth-order solution Y0(tt1) and i2 = −1. Now, in the first-order problem the forcing in the right hand side of the differential equation is

$\left[ -3\, A^2\, A^\ast - 2\, i\, \frac{dA}{dt_1} \right]\, e^{+it} - A^3\, e^{+3it} + c.c.$

where c.c. denotes the complex conjugate of the preceding terms. The occurrence of secular terms can be prevented by imposing on the – yet unknown – amplitude A(t1) the solvability condition

$-3\, A^2\, A^\ast - 2\, i\, \frac{dA}{dt_1} = 0.$

The solution to the solvability condition, also satisfying the initial conditions y(0) = 1 and dy/dt(0) = 0, is:

$A = \tfrac12\, \exp \left(\tfrac38\, i \, t_1 \right).$

As a result, the approximate solution by the multiple-scales analysis is

$y(t) = \cos \left[ \left( 1 + \tfrac38\, \varepsilon \right) t \right] + \mathcal{O}(\varepsilon),$

using t1 = εt and valid for εt = O(1). This agrees with the nonlinear frequency changes found by employing the Lindstedt–Poincaré method.

Higher-order solutions – using the method of multiple scales – require the introduction of additional slow scales, i.e.: t2 = ε2 t, t3 = ε3 t, etc. However, this introduces possible ambiguities in the perturbation series solution, which require a careful treatment (see Kevorkian & Cole 1996; Bender & Orszag 1999).[2]

Alternatively, modern sound approaches derive these sorts of models using coordinate transforms[3] as also described next.

### Coordinate transform to amplitude/phase variables

We seek a solution $y\approx r\cos\theta$ in new coordinates $(r,\theta)$ where the amplitude $r(t)$ varies slowly and the phase $\theta(t)$ varies at an almost constant rate, namely $d\theta/dt\approx 1$. Straightforward algebra finds the coordinate transform

$y=r\cos\theta +\frac1{32}\varepsilon r^3\cos3\theta +\frac1{1024}\varepsilon^2r^5(-21\cos3\theta+\cos5\theta)+\mathcal O(\varepsilon^3)$

transforms Duffing's equation into the pair that the radius is constant $dr/dt=0$ and the phase evolves according to

$\frac{d\theta}{dt}=1 +\frac38\varepsilon r^2 -\frac{15}{256}\varepsilon^2r^4 +\mathcal O(\varepsilon^3).$

That is, Duffing's oscillations are constant amplitude but a different frequencies depending upon the amplitude.[4]

More difficult examples are better treated using a time dependent coordinate transform involving complex exponentials (as also invoked in the previous multiple time scale approach). A web service will perform the analysis for a wide range of examples.[5]