Finite element simulation of fully non‐linear interaction between vertical cylinders and steep waves. Part 1: methodology and numerical procedure

A methodology for computing three‐dimensional interaction between waves and fixed bodies is developed based on a fully non‐linear potential flow theory. The associated boundary value problem is solved using a finite element method (FEM). A recovery technique has been implemented to improve the FEM solution. The velocity is calculated by a numerical differentiation technique. The corresponding algebraic equations are solved by the conjugate gradient method with a symmetric successive overrelaxation (SSOR) preconditioner. The radiation condition at a truncated boundary is imposed based on the combination of a damping zone and the Sommerfeld condition. This paper (Part 1) focuses on the technical procedure, while Part 2 [Finite element simulation of fully non‐linear interaction between vertical cylinders and steep waves. Part 2. Numerical results and validation. International Journal for Numerical Methods in Fluids 2001] gives detailed numerical results, including validation, for the cases of steep waves interacting with one or two vertical cylinders.