Handling geometric singularities efficiently with the immersed-boundary method: application to riblets

Abstract

We present an enhancement of classical immersed boundary methods for the direct numerical simulation of the incompressible Navier-Stokes equations, designed to efficiently cope with geometrical singularities. The method, while general, is described here in the context of the flow over a plane wall covered by small longitudinal riblets. The enhancement consists in prescribing that the solution near the boundary matches the steady Stokes solution. In this way, sharp corners, where velocity gradients are strongest and viscous effects dominate, become tractable with grids comparable to those necessary for a smooth wall. The method requires the availability of a local Stokes solution, that can be obtained either analytically or numerically. Two test cases are discussed: the laminar flow over one riblet, which allows to measure protrusion heights, and a fully turbulent channel flow with ribbed walls. In both cases, the method reproduces reference results at a substantially lower computational cost.