摘要

Rhie-Chow interpolation is a commonly used method in CFD calculations on a co-located mesh in order to suppress non-physical pressure oscillations arising from chequerboard effects. A fully parallelized smoothed-interface immersed boundary method on a co-located grid is described in this paper. We discuss the necessity of modifications to the original Rhie-Chow interpolation in order to deal with a locally refined mesh. Numerical simulation with the modified scheme of Choi shows that numerical dissipation due to Rhie-Chow interpolation introduces significant errors at the immersed boundary. To address this issue, we develop an improved Rhie-Chow interpolation scheme that is shown to increase the accuracy in resolving the flow near the immersed boundary. We compare our improved scheme with the modified scheme of Choi by parallel simulations of benchmark flows: (i) flow past a stationary cylinder; (ii) flow past an oscillating cylinder; and (iii) flow past a stationary elliptical cylinder, where Reynolds numbers are tested in the range 10-200. Our improved scheme is significantly more accurate and compares favourably with a staggered grid algorithm. We also develop a scheme to compute the boundary force for the direct-forcing immersed boundary method efficiently.