Numerical simulation of viscoelastic fluid–structure interaction benchmarks and their application to the human eye

Abstract

We present a numerical solution method for time-dependent viscoelastic fluid–structure interaction employing the arbitrary Lagrangian Eulerian framework. The derived monolithic variational formulation is discretized in time using the shifted Crank–Nicolson scheme and in space using the finite element method. For the linearisation we employ Newton’s method with exact Jacobians. The viscoelastic fluid is modelled either using the Oldroyd-B or a Burgers-type model. The elastic structures are non-linear hyperelastic materials. We validate the implementation on benchmark problems and numerically analyse the convergence for global mesh refinement and adaptive mesh refinement using the dual-weighted residual method. Furthermore we numerically analyse the influence of the viscoelasticity of the fluid on typical goal functionals like the drag, the lift and the displacement. The derived numerical solution method is applied to ophthalmology where we analyse the interaction of the viscoelastic vitreous with its surrounding elastic structures.