A Phase Field Formulation of Frictional Sliding Contact for 3D Fully Eulerian Fluid Structure Interactions

Abstract: Frictional sliding contact in hydrodynamic environments can be found in a range of engineering applications. Accurate modeling requires an integrated numerical framework capable of resolving large relative motions, multiphase interactions, and nonlinear contact responses. Building on our previously developed fully Eulerian fluid structure formulation, we introduce a phase field based formulation for dynamic frictional contact in 3D. Contact detection is achieved via the overlap of diffuse interfaces of colliding solids. The normal contact response is defined as a volumetric body force proportional to the overlap parameter, while the tangential response is computed using the Coulomb friction model. The direction of the friction forces are derived by projecting phase-averaged relative velocities onto the local tangent plane of colliding bodies. This proposed unified treatment enables the computation of both normal and frictional forces within a single momentum balance equation, avoiding separate velocity fields for individual solids. We present several test cases with increasing complexity to verify and demonstrate our proposed frictional contact model. Verification against the Hertzian contact problem shows excellent agreement with the analytical solution, with errors below $3\%$ in the traction profile. In the sliding block benchmark, the computed displacement profiles closely follow the analytical solution for point-mass systems across multiple friction coefficients. The ironing problem demonstrates stable force predictions under finite deformation, with normal and tangential forces matching kinetic friction laws. The robustness and scalability of the proposed formulation are further demonstrated through a representative ship ice interaction scenario with free surface and frictional sliding effects.