where q α determines the shape of the interpolation (e,g. 0 makes it linear). α min is zero in the current study, but the finite value of α max should be chosen carefully as it must be small enough for the stability of the employed numerical methods but also large enough to make the fluid flow negligible in the solid part. From a purely volumetric penalty approach, this value can be set as :