A numerical scheme to simulate time-domain motion responses of multibody floating systems has been successfully proposed. This scheme is integrated into a time-domain simulation tool, with fully coupled hydrodynamic coefficients obtained from the hydrodynamic software - WAMIT which solves the Boundary Value Problem (BVP). The equations of motion are transformed into standard state-space format, using the constant coefficient approximation and the impulse response function method. Thus the Ordinary Differential Equation (ODE) solvers in MATLAB can be directly employed. The time-domain responses of a single spar at sea are initially obtained. The optimal Linear Quadratic Regulator (LQR) controller is further applied to this single spar, by assuming that the Dynamic Positioning (DP) system can provide the optimized thruster forces. Various factors that affect the controlling efficiency, e.g., the time steps ∆τ and ∆t, the weighting factors(Q,R), are further investigated in detail. Next, a two-body floating system is studied. The response amplitude operators (RAOs) of each body are calculated and compared with the single body case. Then the effects of the body-to-body interaction coefficients on the time-domain responses are further investigated. Moreover, the mean drift force is incorporated in the DP system to further mitigate the motion responses of each body. Finally, this tool is extended to a three-body floating system, with the relative motions between them derived.

The hydrodynamic properties of a horizontal cylinder which is free to pitch about an off-centred axis are studied and used to derive the equations of motion of a wave energy converter which extracts energy from incoming sea waves with a linear power-take-off mechanism. The present work follows from a recent study which compared the performance of an off-centred cylinder with those of the Edinburgh Duck wave energy converter. The small decrease in performance found is offset by a reduction in the likely costs associated with the manufacturing of the cylindrical cam compared with those of the asymmetric profile. As part of the survivability strategy in very energetic seas-states it had been planned to completely submerge the device so as to reduce the mooring forces. However, experiments with scale models show that a good absorption capacity is retained even when fully-submerged. The hydrodynamic properties of a horizontal cylinder that pierces the free-surface and of one that is fully submerged are therefore of central concern in this study. These properties are well known for the case of very long cylinders but they are now found for cylinders with different widths, drafts, submergence levels and water-depths. The hydrodynamic forces and moments at the off-centred axis are, furthermore, derived through the application of transformation formulae. The equation of motion of the off-centred cylinder is derived for one degree of freedom and its performance as a wave energy converter is analysed. A relationship which relates the resonance of the device with the location of the off-centred axis and its mass distribution is derived and used to optimize the design for average sea conditions attained at a real location. Design cases associated with three diameters of the cylinder are looked into detail for both a fully-submerged and free-surface piercing cylinder. The one degree of freedom model is extended to include a multi-body which has three degrees of freedom in order to describe the dynamics of a proposed wave powered desalination system based on a cylindrical Duck device. This mathematical model is derived through linearised Lagrangian equations of motion in which the hydrodynamic forces are included as generalised external forces. The advantage of such approach is to reduce the number of equations associated with multi-body systems by removing the reaction forces of holomonic constraints from the system of equations to solve. This model is validated through experiments with a scale model performed in the curved tank of the University of Edinburgh with both regular waves and mixed seas.

Wave scattering by arrays of cylinders has received special attention by many authors and analytical solutions have been derived. The investigation of optical and acoustical analogies to the problem of interaction of water waves with rigid and flexible cylinder arrays is the main focus of this thesis. In acoustics, a sound may be attenuated while it propagates through a layer of bubbly liquid. In fact, if the natural frequency of the bubbles is in the range of the wave periods, the attenuation becomes more evident. The ultimate objective of the research described herein is to determine if this phenomenon may also be found in the interaction between water waves and arrays of flexible cylinders. In a first approach, arrays of rigid cylinders are studied in shallow water. The array is treated as an effective medium, which allows for the definition of reflection and transmission coefficients for the array, and theories from Hu and Chan (2005) associated with the Fabry-Perot interferometer are compared against direct computations of wave scattering using the commercial code WAMIT. Reflection and transmission coefficients from WAMIT are evaluated by applying a Maximum Likelihood Method. The results from WAMIT were found to be in good agreement with those obtained from the effective medium theory. Due to observed inconsistencies for short wave periods and small incident angles, the effective width of the medium is defined and corrected. For the case of a flexible cylinder, generalized modes corresponding to deformations of the cylinder's surface are formulated and added to WAMIT's subroutine. Equations of motion are derived from the theory of vibration for thin shells and mass and stiffness matrices are defined. The objective is to maximize wave attenuation from the array of flexible cylinders. Therefore, the natural periods of the "breathing" mode for these cylinders is set in the range of the studied wave periods. Then, material properties, as well as mass and stiffness matrices, are chosen to achieve this effect.

The vast and rich body of literature covering the numerical modelling of hydrodynamic floating body systems has demonstrated their great power and versatility when applied to offshore marine energy systems. It is possible to model almost any type of physical phenomenon which could be expected within such a system, however, limitations of computing power continue to restrict the usage of the most comprehensive models to very narrow and focused design applications. Despite the continued evolution of parallel computing, one major issue that users of computational tools invariably face is how to simplify their modelled systems in order to achieve practically the necessary computations, whilst capturing enough of the pertinent physics, with great enough ‘resolution’, to give robust results. The challenge is, in particular, to accurately deliver a complete spectrum of results, that account for all of the anticipated sea conditions and allow for the optimisation of different control scenarios. This thesis examines the uncertainty associated with the effects of viscosity and nonlinear behaviour on a small scale model of an oscillating system. There are a wide range of Computational Fluid Dynamics (CFD) methods which capture viscous effects. In general however, the oscillating, six degree-of-freedom floating body problem is best approached using a linear potential flow based Boundary Element Method (BEM), as the time taken to process an equivalent model will differ by several orders of magnitude. For modelling control scenarios and investigating the effects of different sea states, CFD is highly impractical. As potential flows are inviscid by definition, it is therefore important to know how much of an impact viscosity has on the solution, particularly when different scales are of interest during device development. The first aim was to develop verified and validated solutions for a generic type decaying system. The arrangement studied was adapted from an array tank test experiment which was undertaken in 2013 by an external consortium (Stratigaki et al., 2014). Solutions were delivered for various configurations and gave relatively close approximations of the experimental measurements, with the modelling uncertainties attributed to transient nonlinear effects and to dissipative effects. It was not possible however to discern the independent damping processes. A set of CFD models was then developed in order to investigate the above discrepancies, by numerically capturing the nonlinear effects, and the effects of viscosity. The uncontrolled mechanical effects of the experiment could then be deduced by elimination, using known response patterns from the measurements and derived results from the CFD simulations. The numerical uncertainty however posed a significant challenge, with the outcomes supported by verification evidence, and detailed discussions relating to the model configuration. Finally, the impact of viscous and nonlinear effects were examined for two different interacting systems – for two neighbouring devices, and an in-line array of five devices. The importance of interaction behaviour was tested by considering the transfer of radiation forces between the model wave energy converters, due to the widely accepted notion that array effects can impact on energy production yields. As there are only very limited examples of multi-body interaction analysis of wave energy devices using CFD, the results with this work provide important evidence to substantiate the use of CFD for power production evaluations of wave energy arrays. An effective methodology has been outlined in this thesis for delivering specific tests to examine the effects of viscosity and nonlinear processes on a particular shape of floating device. By evaluating both the inviscid and viscous solutions using a nonlinear model, the extraction of systematic mechanical effects from experimental measurements can be achieved. As these uncontrolled frictional effects can be related to the device motion in a relatively straightforward manner, they can be accommodated within efficient potential flow model, even if it transpires that they are nonlinear. The viscous effects are more complex; however, by decomposing into shear and pressure components, it may in some situations be possible to capture partially the dynamics as a further damping term in the efficient time-domain type solver. This is an area of further work.