Finite Element Simulation of Fully Non-Linear Interaction Between Vertical Cylinders and Steep Waves

The document summarizes a finite element method for simulating the interaction between steep ocean waves and fixed vertical cylinders. An iterative method with preconditioning is used to solve the algebraic equations arising from the discretization. Linear shape functions are used in the finite element formulation. The fluid domain is divided into hexahedral elements to discretize the governing Laplace equation and boundary conditions at the free surface, rigid surfaces, and far-field boundary.

Uploaded by

yaseromidi67

Available Formats

Download as PDF, TXT or read online on Scribd

0% found this document useful (0 votes)

28 views

Finite Element Simulation of Fully Non-Linear Interaction Between Vertical Cylinders and Steep Waves

Uploaded by

yaseromidi67

Available Formats

Download as PDF, TXT or read online on Scribd

You are on page 1/ 3

82011209 Subject : Finite element simulation of fully non-linear interaction between vertical cylinders and steep waves 1.

1. Introduction : In the present paper we apply the FEM to analyse the fully non-linear three dimensional interaction between waves and fixed bodies. An iterative method with a preconditioner is used to solve the algebraic equations. A postprocessor based on the recovery technique is adopted in order to improve the accuracy of the finite element solution without increasing the number of elements. A radiation condition is imposed, which is to minimize the effect of the reflection of waves from the truncated boundary of the fluid domain.

Figure 1. Sketch of the problem

2. Governing equation and boundary conditions of fluid motion Assumptions: The viscosity, surface tension and compressibility are neglected.( ) Governing equation: The fluid flow is then irrotational so that a velocity potential exists. Boundary conditions:

where n is the normal vector of the surface pointing out of the fluid domain. 3. Finite Element Formulation At each time step, a mixed boundary value problem has to be solved in which the unknown functions satisfy the Laplace equation in the fluid domain, the Dirichlet condition on the free surface and on the far end, and the Neumann condition on the rigid surface. The FEM is used in which the fluid domain is discretized into a set of small elements and the velocity potential is expressed in terms of a shape function.

Where

is the velocity potential at the nodes.

In this work, linear shape functions are used because this choice makes the calculation of influence coefficients relatively simple.Using the Galerkin method, the Laplace equation and the boundary conditions satisfied by the velocity potential can be discretized in the following form:

This Equation can further be written in the matrix form:

Where

Discretization of the fluid domain: The discretization of the fluid domain into a mesh of small elements is a fundamental part of the numerical procedure. The work presented here is primarily directed towards wave interactions with vertical cylinders. This allows some relatively simple techniques to be used. The fluid domain is first divided into a number of small hexahedra.