Fourier spectral and wavelet solvers for the incompressible Navier-Stokes equations with volume-penalization: Convergence of a dipole-wall collision

  • Authors:

  • G. H. Keetels;U. D'Ortona;W. Kramer;H. J. H. Clercx;K. Schneider;G. J. F. van Heijst

  • Affiliations:

  • Department of Physics, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands;Laboratoire de Modèlisation et Simulation Numèrique en Mècanique et Gènie des Procèdès, CNRS and Universitès d'Aix-Marseille, 38 rue F. Joliot-Curie, 13451 Marse ...;Department of Physics, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands;Department of Physics, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands and Department of Applied Mathematics, University of Twente, P.O. Box 217, 7500 AE Ensch ...;Laboratoire de Modèlisation et Simulation Numèrique en Mècanique et Gènie des Procèdès, CNRS and Universitès d'Aix-Marseille, 38 rue F. Joliot-Curie, 13451 Marse ...;Department of Physics, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2007

Quantified Score

Hi-index 31.49

Visualization

Abstract

In this study, we use volume-penalization to mimic the presence of obstacles in a flow or a domain with no-slip boundaries. This allows in principle the use of fast Fourier spectral methods and coherent vortex simulation techniques (based on wavelet decomposition of the flow variables) to compute turbulent wall-bounded flow or flows around solid obstacles by simply adding one term in the equation. Convergence checks are reported using a recently revived, and unexpectedly difficult dipole-wall collision as a benchmark computation. Several quantities, like the vorticity isolines, truncation error, kinetic energy and enstrophy are inspected for a collision of a dipole with a no-slip wall and compared with available benchmark data obtained with a standard Chebyshev pseudospectral method. We quantify the possible deteriorating effects of the Gibbs phenomenon present in the Fourier based schemes due to continuity restrictions of the penalized Navier-Stokes equations on the wall. It is found that Gibbs oscillations have a negligible effect on the flow evolution allowing higher-order recovery of the accuracy on a Fourier basis by means of postprocessing. An advantage of coherent vortex simulations, on the other hand, is that the degrees of freedom of the flow computation can strongly be reduced. In this study, we quantify the possible reduction of degrees of freedom while keeping the accuracy. For an optimal convergence scenario the penalization parameter has to scale with the number of Fourier and wavelet modes. In addition, an implicit treatment of the Darcy drag term in the penalized Navier-Stokes equations is beneficial since this allows one to set the time step independent from the penalization parameter without additional computational or memory requirements.