A second-order cut-cell method for the numerical simulation of 2D flows past obstacles.

François Bouchon (1) , Thierry Dubois (1) , Nicolas James (1)

1 LMBP - Laboratoire de Mathématiques Blaise Pascal

We present a new second-order method, based on the MAC scheme on cartesian grids, for the numerical simulation of two-dimensional incompressible flows past obstacles. In this approach, the solid boundary is embedded in the cartesian computational mesh. Discretizations of the viscous and convective terms are formulated in the context of finite volume methods ensuring local conservation properties of the scheme. Classical second-order centered schemes are applied in mesh cells which are sufficiently far from the obstacle. In the mesh cells cut by the obstacle, fisrt-order approximations are proposed. The resulting linear system is nonsymetric but the stencil remains local as in the classical MAC scheme on cartesian grids. The linear systems are solved by a direct method based on the capacitance matrix method. The time integration is achieved with a second-order projection scheme. While in cut-cells the scheme is locally first-order, a global second-order accuracy is recovered. This property is assessed by computing analytical solutions for a Taylor-Couette problem. The efficiency and robustness of the method is supported by numerical simulations of 2D flows past a circular cylinder at Reynolds number up to 9 500. Good agreement with experimental and published numerical results are obtained.

Mots clés: Immersed boundary method; Cartesian grid method; Cut-cell method; Incompressible viscous flows; Complex geometry

Domaines: Analyse numérique [math.NA]; Equations aux dérivées partielles [math.AP]; Physique mathématique [math-ph]

François Bouchon, Thierry Dubois, Nicolas James. A second-order cut-cell method for the numerical simulation of 2D flows past obstacles.. 2011. ⟨hal-00570049⟩