Revisiting the Slip Boundary Condition: Surface Roughness as a Hidden Tuning Parameter

TL;DR

This study explores how surface roughness influences slip boundary conditions in fluid flow simulations, revealing that numerical discretization and mesh quality significantly affect flow stability and force predictions, challenging existing wall modeling assumptions.

Contribution

It demonstrates that numerical surface roughness and mesh distortion critically impact flow stability and force calculations, highlighting the importance of high-order geometry mappings in simulations.

Findings

Numerical surface roughness can destabilize potential flow solutions.

High-order geometry mappings stabilize flow simulations.

Manipulating roughness and mesh distortion controls drag and lift forces.

Abstract

In this paper, we investigate the effect of boundary surface roughness on numerical simulations of incompressible fluid flow past a cylinder in two and three spatial dimensions furnished with slip boundary conditions. The governing equations are approximated using a continuous finite element method, stabilized with a Galerkin least-squares approach. Through a series of numerical experiments, we demonstrate that: $(i)$ the introduction of surface roughness through numerical discretization error, or mesh distortion, makes the potential flow solution unstable; $(ii)$ when numerical surface roughness and mesh distortion are minimized by using high-order isoparametric geometry mappings, a stable potential flow is obtained in both two and three dimensions; $(iii)$ numerical surface roughness, mesh distortion and refinement level can be used as control parameters to manipulate drag and lift forces resulting in numerical values spanning more than an order of magnitude. Our results cast some doubt on the predictive capability of the slip boundary condition for wall modeling in turbulent simulations of incompressible flow.