1.3. Basdevant reduction

Writing the nonlinear term as (ωu)\boldsymbol{\nabla}\cdot(\omega\mathbf{u}), we need to apply the inverse fft to three functions, namely on ω^\hat\omega, u^\hat u, v^\hat v, and then two direct fft, namely on ωu\omega u and ωv\omega v, before taking the gradient in spectral space.

Alternatively, we can keep it as ωu\boldsymbol{\nabla}\omega \cdot \mathbf{u} and compute four ifft and one fft.

The best, however, is to perform the Basdevant reduction to reduce to four ffts (two inverse and two direct). This is obtained by writing

(ωu)=(x2y2)(uv)+xy(v2u2). \boldsymbol{\nabla} \cdot (\omega\mathbf{u}) = \left(\partial_x^2 - \partial_y^2\right)(uv) + \partial_{xy}\left(v^2 - u^2\right).

This way, we need to compute the inverse fft of uu and vv, compute the product uvuv and the term v2u2v^2 - u^2 in physical space, and finally compute the FFT of these two functions.



Last modified: March 21, 2022. Built with Franklin.jl, using the Book Template.