Improving convergence of volume penalised fluid-solid interactions

From MaRDI portal
Publication:6316331

DOI10.1016/J.JCP.2020.110043arXiv1903.11914MaRDI QIDQ6316331FDOQ6316331

Eric W. Hester, Keaton J. Burns, G. M. Vasil

Publication date: 28 March 2019

Abstract: We analyse and improve the volume-penalty method, a simple and versatile way to model objects in fluid flows. The volume-penalty method is a kind of fictitious-domain method that approximates no-slip boundary conditions with rapid linear damping inside the object. The method can then simulate complex, moving objects in general numerical solvers without specialised algorithms or boundary-conforming grids. Volume penalisation pays for this simplicity by introducing an equation-level error, the extitmodelerror, that is related to the damping time etall1. While the model error has been proven to vanish as the damping time tends to zero, previous work suggests convergence at a slow rate of mathcalO(eta1/2). The stiffness of the damping implies conventional volume penalisation only achieves first order numerical accuracy. We analyse the volume-penalty method using multiple-scales matched-asymptotics with a signed-distance coordinate system valid for arbitrary smooth geometries. We show the dominant model error stems from a displacement length that is proportional to a Reynolds number extRe dependent boundary layer of size mathcalO(eta1/2extRe1/2). The relative size of the displacement length and damping time leads to multiple error regimes. Our key finding derives a simple smoothing prescription for the damping that eliminates the displacement length and reduces the model error to mathcalO(eta) in all regimes. This translates to second order numerical accuracy. We validate our findings in several comprehensive benchmark problems and finally combine Richardson extrapolation of the model error with our correction to further improve convergence to mathcalO(eta2).












This page was built for publication: Improving convergence of volume penalised fluid-solid interactions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6316331)