Results are presented for the 3D numerical simulation of the water impact of a wave energy converter in free fall and subsequent heave motion. The solver, AMAZON-3D, employs a Riemann-based finite volume method on a Cartesian cut cell mesh. The computational domain includes both air and water regions with the air/water boundary captured automatically as a discontinuity in the density field thereby admitting break up and recombination of the free surface. Temporal discretisation uses the artificial compressibility method and a dual time stepping strategy. Cartesian cut cells are used to provide a boundary-fitted grid at all times. The code is validated by experimental data including the free fall of a cone and free decay of a single Manchester Bobber component.
A popular class of wave energy converter (WEC) consists of floating bodies which oscillate with one or more degrees of freedom and whose horizontal dimensions are small in comparison to the wave length. Such bodies are point absorbers and they essentially convert their heave motion into useful energy. Several such WECs composed of one or more point absorbers are currently under development. Examples include the Manchester Bobber [
In this paper, attention is focused on a particular WEC: the Manchester Bobber henceforth referred to as “the Bobber”. The Bobber was designed by researchers at the University of Manchester in the UK and consists of an array of novel heaving point absorbers which generate oscillatory shaft motion that is converted to unidirectional rotation through a freewheel/clutch system which in turn drives an electricity generator. In order to assess the survivability of the Bobber, a full set of flow variables is required e.g. pressure and velocity field, as well as integrated effects like forces and device response as outputs from laboratory experiments or CFD predictions. In this paper, the numerical study is focused on an isolated Bobber under water impacts but the same code could be used to simulate violent wave impacts.
Research work on the water impact problem has been mainly for 2D cases. Greenhow and Lin [
Our in-house AMAZON-3D finite volume code solves the incompressible Navier-Stokes equations in both air and water regions and treats the free surface as a contact surface in the density field that is captured automatically in a manner analogous to shock capturing in compressible flow. Meshing uses the Cartesian cut cell method which automatically produces a body fitted mesh for static or moving solid bodies. Solid objects are cut out of a background Cartesian mesh leaving a set of irregularly shaped cells whose surfaces coincide with the boundary of the solid. There is no requirement to re-mesh globally for the case of a moving body. All that is required is to update the cut cell data at the body surface for as long as the body motion continues. The background mesh needs not be uniform and in fact non-uniform Cartesisn meshes were used in our simulations with small cells around the body for improved resolution. We have used the code to study wave loading on a Bobber-type WEC device under extreme wave conditions in a numerical wave tank (Hu et al. [
This paper focuses on the numerical modelling of the vertical slamming on a floating Bobber and the associated pressures that might be expected when slamming occurs. The code has been validated by the free fall of a cone and free decay of a Bobber. Results including the vertical displacement will be presented for the water impact of bodies in free fall motion.
An essential component of a CFD simulation is the mesh generation and it is particularly important for the present case because of the movement of the solid body within computational domain.
It is well known that the finite volume method (FVM) involves discretization of the flow domain of interest and then integration of the flow equations over elemental cell volumes. The method enables correct flux balances across cell boundaries and conserves momentum throughout the grid. Therefore, advantage may be gained by using the FVM since the dependent variables remain at all times referenced to a Cartesian frame even when boundary-fitted using cut cells. The remainder (majority) of the cells are uncut flow cells that be treated in a straight forward manner. However, whilst the Cartesian cut cell algorithms can easily accommodate moving boundaries, there are pathological cases where the approach can sometimes provide relatively poor resolution of some particular geometric features. For example, a numerical instability may occur locally within the flow solver if cut cells become arbitrarily small. To overcome this problem, cell merging is implemented as in Clarke et al. [
The integral form of the Euler equations for 3D incompressible flow with variable density can be written as
where Q is the vector of conserved variables which encloses the time dependent domain of interest V, F is the flux vector function and n is the outward unit vector normal to the boundary S. B is a source term for body forces. Q, F and B are given by
where
and u, v and w are the flow velocity components and ub, vb and wb are the velocity components of the boundary S which are zero when the boundary is stationary. ρ is the density, p is the pressure, β is the coefficient of artificial compressibility and g is the gravitational acceleration.
We can then discretize Equation (1) over each cell within the flow domain using a finite volume formulation, this gives
where is the average value of Q in cell stored at the cell centre and denotes the volume of the cell. Fk is the numerical flux across the face k of the cell, ΔAk is the area of the face and m is the number of faces of the cell. The convective flux Fk is evaluated using Roe’s approximate Riemann solver, which assumes a 1D Riemann problem in the direction normal to the cell face and has the form
where and are the reconstructed values on the right and left at face k and A is the flux Jacobian evaluated by Roe’s average state. The quantities R, L and are the right and left eigenvectors of A and the eigenvalues of A respectively.
To achieve a time-accurate solution at each time step of the unsteady flow problem a first-order Euler implicit difference scheme is used to discretise Equation (5) as
Introducing a pseudo-time derivative into Equation (7), this gives
where τ is the pseudo-time and. The right-hand side of Equation (8) can be linearized using Newton’s method at the m+1 pseudo-time level and then can be written in the matrix form
where
and Im should be defined as
.
When is iterated to zero, the density and momentum equations are satisfied, and the divergence of the velocity at time level n + 1 is zero. The system of equations can be written in matrix form as
where D is a block diagonal matrix, L is a block lower triangular matrix and U is a block upper triangular matrix. Each of the elements in D, L and U is a matrix. An approximate LU factorization (ALU) scheme as proposed by Pan and Lomax [
Within each time step of the implicit integration the sub-iteration is terminated when the L2 norm of the iterates process
is less than a specified limit ε and in this paper.
For all test cases non-reflecting boundary conditions are applied at the top boundary allowing air to leave or enter the domain freely. The remaining boundaries are set as rigid walls.
At the interface between two immiscible fluids, the present method assumes that the system of equations for non-homogeneous, incompressible flow can treat the free surface numerically as a contact discontinuity in the density field. Special procedures to track the free surface are thus unnecessary since the free surface is captured automatically. It is asserted that the numerical solution of Equation (1) for a system containing one or more free surfaces will converge to the correct unique solution.
In the current approach, the pressure value p can be calculated from the term p/β in Equation (2) by multiplying by β. The total force is then obtained by integrating the pressure field along the body surface
, where Sb is the area of the body surface as defined by the cut cell surface fitting.
To simulate water impact on floating body, the body must be allowed to move freely according to the interaction between the fluid and the body. In this paper, our numerical simulation with body motion is used to heave only, which is main part of the body motion in this problem. In the following calculations, the density ratio of water to air is taken as 1000:1. The location of the free surface, i.e. the air/water interface, is defined as the density contour with the value 500 kg∙m−3. The value of the gravitational acceleration is taken as 9.8 ms−2. The typical physical time step used is and the pseudo time step is set as. The value of the artificial compressibility parameter is.
The first test involves a free falling of a cone with deadrise angle 20˚ into still water. The physical tank tests are described by Backer et al. [
In the physical experiments, the time history of pressure on the cone surface was recorded during the impact at two locations (S1 and S2) with a horizontal distance of 0.04 m and 0.09 m respectively from the symmetrical axis as shown in
ing part on the right side of the peak. However, there is a substantial difference in the peak pressure values between the experimental and numerical data. For example, the values of the peak pressure are 0.455 bar for the finer mesh and 0.404 bar for the coarser mesh but for the experimental this value is 0.719 bar. This indicates further mesh refinements especially for the region near the surface of the cone are still needed. To show this and to save the computation time, an axi-symmetric version of the free surface code has been applied to simulate the same flow problem. Three levels of grid (dx = dy = 0.005, 0.0025 and 0.00125 m respectively) have been used with a time step of . In
From these results, it can be concluded that for the test case of a free falling cone impacting on water surface the occurrence of peak surface pressures is an extremely localised phenomenon in both space and time, so ex-
tremely small mesh size and time step at least locally must be used in order to capture those moments accurately. This highlights the need for both spatially and temporally adaptive solutions for such flow problems. On the other hand, the integrated global values such as the impact force, body acceleration and moving speed are much less affected by the localised extreme pressure values. For the same interval time as above,
The small differences in velocity may be due to the laboratory model not being constrained to heave motion only. Other reasons which explain the differences may be due to the volume of cone (or associated buoyancy force) being slightly different from the laboratory model owing to small inaccuracies in the surface fitting of the geome-
try used in the cut cell model (see
The simulations were carried out using one processor of a 600 MHz NEC vector computer with the CPU times of 2 days for the coarse mesh and about 5 days for the finer mesh.
This test case simulates the free-decay of the Bobber as described in the physical tank tests by Thomas et al. [
where the total mass equals the Bobber mass of mf plus the counterweight mass mc, is the acceleration in the vertical direction positive upward and the integration on the right-hand side calculates the dynamic force from the pressure on the wetted area of the Bobber in the vertical direction only. According to the experimental tests by Thomas et al. [
In this test, the geometry of the Bobber is illustrated in
1.0 m with a water depth of 0.5 m. A non-uniform mesh is used with 88 × 88 × 48 = 371,712 cells with a mesh spacing of 0.016 m in the refined regions close to the body geometry.
In this test, the initial position of the Bobber in terms of its base is set at 0.16 m below the still water surface (see
In this paper, the motions of solid rigid bodies entering water through free fall moving in heave motion were simulated numerically. A number of tests were performed including the water entry of a rigid cone and a generic heave-type wave energy device called the Man-
chester Bobber. Additionally, free decay tests for the Bobber were carried out. Where available, the numerical results were compared to those from laboratory experiments, which show the potential of the present simulation method for dealing with 3D water/body interaction problems. In particular, the code is able to handle moving solid bodies of arbitrary complexity via the Cartesian cut cell approach and also to deal with complex free surfaces by capturing them automatically as contact discontinuities in the density field. The code has shown the capability for simulating heaving wave energy converters and is able to predict impact forces from slamming which will be useful for survivability assessment studies. Although not demonstrated here, the code is also able to deal with general wave/body interactions and so, given sufficient compute power, would be able to model the performance of generic wave energy converters.
This work was supported by the UK Engineering and Physical Sciences Research Council (EPSRC) under grant reference EP/D077621 for which the authors are most grateful. The authors would also like to thank T. Stallard, S. Weller and S. Thomas (Manchester University), for their experimental data for the free decay tests of the Manchester Bobber.