## 2.7.1 Fourier Analysis of PDEs

Measurable Outcome 2.2, Measurable Outcome 2.11

We will develop Fourier analysis in one dimension. The basic ideas extend easily to multiple dimensions. We will consider the convection-diffusion equation,

We will assume that the velocity, \(u\), and the viscosity, \(\mu\) are constant.

The solution is assumed to be periodic over a length \(L\). Thus,

where \(m\) is any integer.

A Fourier series with periodicity over length \(L\) is given by,

\(k_ m\) is generally called the wavenumber, though \(m\) is the number of waves occurring over the length \(L\). We note that \(\hat{U}_ m(t)\) is the amplitude of the \(m\)-th wavenumber and it is generally complex (since we have used complex exponentials). Substituting the Fourier series into the convection-diffusion equation gives,

Noting that \(i^2 = -1\) and collecting terms gives,

The next step is to utilize the orthogonality of the different Fourier modes over the length \(L\), specifically,

By multiplying Equation (2.129) by \(e^{-i k_ n x}\) and integrating from \(0\) to \(L\), the orthogonality condition gives,

Thus, the evolution of the amplitude for an individual wavenumber is independent of the other wavenumbers. The solution to Equation (2.131),

The convection term, which results in the complex time dependent behavior, \(e^{-i u k_ n t}\), only oscillates and does not change the magnitude of \(\hat{U}_ n\). The diffusion term causes the magnitude to decrease as long as \(\mu > 0\). But, if the diffusion coefficient were negative, then the magnitude would increase unbounded with time. Thus, in the case of the convection-diffusion PDE, as long as \(\mu \geq 0\), this solution is stable.

## Exercise

For the convection-diffusion equation analyzed above, how does the rate of damping compare for a wavelength \(L\) mode to a wavelength \(L/10\) mode?