Computational framework for monolithic coupling for thin fluid flow in contact interfaces

TL;DR

This paper introduces a comprehensive computational framework for simulating thin fluid flow in contact interfaces, effectively coupling fluid and solid mechanics, handling complex geometries, contact constraints, and trapped fluid zones, with demonstrated robustness and applicability.

Contribution

The paper presents a novel monolithic finite-element framework for coupled fluid-solid contact problems, including fluid entrapment and non-linear pressurisation, suitable for both one- and two-way coupling approaches.

Findings

Robust convergence demonstrated on model geometries.

Effective simulation of fluid flow and contact pressure evolution.

Quantification of trapped fluid pools impact on interface properties.

Abstract

We developed a computational framework for simulating thin fluid flow in narrow interfaces between contacting solids, which is relevant for a range of engineering, biological and geophysical applications. The treatment of this problem requires coupling between fluid and solid mechanics equations, further complicated by contact constraints and potentially complex geometrical features of contacting surfaces. We developed a monolithic finite-element framework for handling mechanical contact, thin incompressible viscous flow and fluid-induced tractions on the surface of the solid, suitable for both one- and two-way coupling approaches. Additionally, we consider the possibility of fluid entrapment in "pools" delimited by contact patches and its pressurisation following a non-linear compressibility constitutive law. Furthermore, image analysis algorithms were adapted to identify the local status of each interface element within the Newton-Raphson loop. First, an application of the proposed framework for a problem with a model geometry is given, and the robustness is demonstrated by the residual-wise and status-wise convergence. The full capability of the developed two-way coupling framework is demonstrated on a problem of a fluid flow in contact interface between a solid with representative rough surface and a rigid flat. The evolution of the contact pressure, fluid flow pattern and the morphology of trapped fluid zones until the complete sealing of the interface is displayed. Additionally, we demonstrated an almost mesh-independent result of a refined post-processing approach to the real contact-area computation. The developed framework permits not only to study the evolution of effective properties of contact interfaces, but also to highlight the difference between one- and two-way coupling approaches and to quantify the effect of multiple trapped fluid "pools" on the coupled problem.