In this study we use a fast Fourier spectral technique to simulate the Navier-Stokes equations with no-slip boundary conditions. This is enforced by an immersed boundary technique called volumepenalization. The approach has been justified by analytical proofs of the convergence with respect to the penalization parameter. However, the solution of the penalized Navier-Stokes equations is not smooth on the surface of the penalized volume. Therefore, it is not a priori known whether it is possible to actually perform accurate fast Fourier spectral computations. Convergence checks are reported using a recently revived, and unexpectedly difficult dipole-wall collision as a test case. It is found that Gibbs oscillations have a negligible effect on the flow evolution. This allows higher-order recovery of the accuracy on a Fourier basis by means of a post-processing procedure.
G. H. Keetels, H. J. H. Clercx, G. J. F. van Heijs