Search

Computational fluid dynamics (CFD) is used in this thesis to model wind and tidal stream turbines and to investigate tidal turbine fence performance. There are two primary objectives of this work. The first is to develop and validate an actuator line method for the simulation of wind and tidal turbines which applies the blade forces to the flow field without the need for a regularisation kernel. The second is to examine tidal fences using, in part, the newly developed actuator line method. A potential flow equivalence method for determining the relative velocity to the blade chord and flow angle at the rotor blades in the actuator line method is proposed and validated. Results for simulations using this method compare favourably with those from both experiments and alternative computational methods, although the present model’s results deviate from experimental results in the vicinity of the blade tips. A CFD-embedded blade element-momentum tool is used to design rotors for operation in infinitely wide tidal fences spanning a tidal channel. Rotors are designed for fences with several different blockage ratios, with those designed for high blockage flows having greater solidity than those designed for operation in fences with lower blockage. It is found that designing rotors for operational blockage conditions can significantly improve the power output achieved by a tidal fence. Improved power output for higher blockage conditions is achieved by the application of greater thrust to the flow. Actuator line simulations of short (up to 8 turbines) fences with varying intra-rotor spacing and number of rotors confirm that hydrodynamic performance of the rotors improves as the spacing is reduced and as rotors are added to a fence. The position of a rotor within the fence impacts its performance; rotors at the ends of a fence extract reduced power compared to those at the centre of the fence, particularly for tip speed ratios greater than the design tip speed ratio.

A three dimensional computational model of a full scale axial flow tidal turbine has been used to investigate the effects of a range of realistic environmental conditions on turbine performance. The model, which is based on the Reynolds averaged Navier-Stokes equations, has been developed using the commercial flow solver ANSYS Fluent. A 1:30 scale tidal turbine is simulated in an open channel for comparison to existing experimental data. The rotor blades are directly resolved using a body-fitted, unstructured computational grid. Rotor motion is enabled through a sliding mesh interface between the rotor and the channel boundaries. Reasonably good agreement in thrust and power is observed. The computed performance curves are offset from the measured performance curves by a small increment in rotor speed. Subsequently, a full scale axial flow turbine is modelled in a variety of conditions representative of tidal channel flows. A parametric study is carried out to investigate the effects of flow shear, confinement and alignment on turbine performance, structural loading, and wake recovery. Mean power and thrust are found to be higher in sheared flow, relative to uniform flow of equivalent volumetric flow rate. Large fluctuations in blade thrust and torque occur in sheared flow as the blade passes through the high velocity freestream flow in the upper portion of the profile and the lower velocity flow near the channel bed. A stronger shear layer is formed around the upper portion of the wake in sheared flow, leading to enhanced wake mixing. Mean power and thrust are reduced when the turbine is simulated at a lower position in a sheared velocity profile. However, fluctuations in blade loading are increased due to the higher velocity gradient. The opposite effects are observed when the turbine operates at greater heights in sheared flow. Flow misalignment has a negative impact on mean rotor thrust and power, as well as on unsteady blade loading. Although the range of unsteady loading is not increased significantly, additional perturbations are introduced due to interactions between the blade and the nacelle. A deforming surface is introduced using the volume-of-fluid method. Linear wave theory is combined with the existing free surface model to develop an unsteady inflow boundary condition prescribing combined sheared flow and free surface waves. The relative effects of the sheared profile and wave-induced velocities on turbine loading are identified through frequency analysis. Rotor and blade load fluctuations are found to increase with wave height and wave length. In a separate study, the performance of bi-directional ducted tidal turbines is investigated through a parametric study of a range of duct profiles. A two dimensional axi-symmetric computational model is developed to compare the ducted geometries with an unducted device under consistent blockage conditions. The best-performing ducted device achieves a peak power coefficient of approximately 45% of that of the unducted device. Comparisons of streamtube area, velocity and pressure for the flow through the ducted device shows that the duct limits the pressure drop across the rotor and the mass flow through the rotor, resulting in lower device power.

Offshore structures are required to survive in extreme wave environments. Historically, the design of these offshore structures and vessels has relied on wave-tank experiments and linear theory. Today, with advances in computing power, it is becoming feasible to supplement these methods of analysis with fully nonlinear numerical simulation. This thesis is concerned with the development of an efficient method to perform this numerical modelling, in the context of potential flow theory. The interaction of a steep ocean wave with a floating body involves a moving free surface and a wide range of length scales. Attempts to reduce the size of the simulation domain cause problems with wave reflection from the domain edge and with the accurate creation of incident waves. A method of controlling the wave field around a structure is presented. The ability to effectively damp an outgoing wave in a very short distance is demonstrated. Steep incident waves are generated without the requirement for the wave to evolve over a large time or distance before interaction with the body. This enables a general wave-structure interaction problem to be modelled in a small tank, and behave as if it were surrounded by a large expanse of water. The suitability of the boundary element method for performing this modelling is analysed. Potential improvements are presented with respect to accuracy, robustness, and computational complexity. Evidence of third order diffraction is found for an FPSO model.

This thesis details the development, verification and validation of an unsteady unstructured high order (≥ 3) h/p Discontinuous Galerkin - Fourier solver for the incompressible Navier-Stokes equations on static and rotating meshes in two and three dimensions. This general purpose solver is used to provide insight into cross-flow (wind or tidal) turbine physical phenomena.Simulation of this type of turbine for renewable energy generation needs to account for the rotational motion of the blades with respect to the fixed environment. This rotational motion implies azimuthal changes in blade aero/hydro-dynamics that result in complex flow phenomena such as stalled flows, vortex shedding and blade-vortex interactions. Simulation of these flow features necessitates the use of a high order code exhibiting low numerical errors. This thesis presents the development of such a high order solver, which has been conceived and implemented from scratch by the author during his doctoral work.To account for the relative mesh motion, the incompressible Navier-Stokes equations are written in arbitrary Lagrangian-Eulerian form and a non-conformal Discontinuous Galerkin (DG) formulation (i.e. Symmetric Interior Penalty Galerkin) is used for spatial discretisation. The DG method, together with a novel sliding mesh technique, allows direct linking of rotating and static meshes through the numerical fluxes. This technique shows spectral accuracy and no degradation of temporal convergence rates if rotational motion is applied to a region of the mesh. In addition, analytical mappings are introduced to account for curved external boundaries representing circular shapes and NACA foils.To simulate 3D flows, the 2D DG solver is parallelised and extended using Fourier series. This extension allows for laminar and turbulent regimes to be simulated through Direct Numerical Simulation and Large Eddy Simulation (LES) type approaches. Two LES methodologies are proposed.Various 2D and 3D cases are presented for laminar and turbulent regimes. Among others, solutions for: Stokes flows, the Taylor vortex problem, flows around square and circular cylinders, flows around static and rotating NACA foils and flows through rotating cross-flow turbines, are presented.

When traditional linearised theory is used to study free-surface flows past a surface-piercing object or over an obstruction in a stream, the geometry of the object is usually lost, having been assumed small in one or several of its dimensions. In order to preserve the nonlinear nature of the geometry, asymptotic expansions in the low-Froude or low-Bond limits can be derived, but here, the solution invariably predicts a waveless free-surface at every order. This is because the waves are in fact, exponentially small, and thus beyond-all-orders of regular asymptotics; their formation is a consequence of the divergence of the asymptotic series and the associated Stokes Phenomenon.In this thesis, we will apply exponential asymptotics to the study of two new problems involving nonlinear geometries. In the first, we examine the case of free-surface flow over a step including the effects of both gravity and surface tension. Here, we shall see that the availability of multiple singularities in the geometry, coupled with the interplay of gravitational and cohesive effects, leads to the discovery of a remarkable new set of solutions.In the second problem, we study the waves produced by bluff-bodied ships in low-Froude flows. We will derive the analytical form of the exponentially small waves for a wide range of hull geometries, including single-cornered and multi-cornered ships, and then provide comparisons with numerical computations. A particularly significant result is our confirmation of the thirty-year old conjecture by Vanden-Broeck & Tuck (1977) regarding the impossibility of waveless single-cornered ships.

Many candidate sites for tidal stream energy extraction can be classified as headland sites. Understanding the tidal resource of such sites is of fundamental importance to the industry. This paper examines an upper bound for the power might can be generated from an idealised headland. The dependence of this on length of turbine rows, number of turbine rows and turbine blockage ratio are examined. Conclusions are drawn from this which are applicable to real headland sites.

If some solid bulk cargoes such as concentrates, unprocessed nickel ores and iron ore fines contain sufficient moisture, then cargo compaction during a voyage may cause liquefaction. This has been known to result in major cargo displacement ('cargo shift'), causing the vessel to capsize. This has resulted in a number of fatalities. To provide context to this interdisciplinary problem, an overview of the key considerations is presented, which includes: (i) a summary of the wide range of factors that influence a cargo's likelihood to shift; (ii) the soil mechanics principles that can be used to explain cargo shift; and (iii) the regulatory controls that are used to prevent its occurrence during shipment. Under the requirements of international maritime regulation, primarily the International Maritime Bulk Solid Cargoes (IMSBC) Code, a cargo that may liquefy must be shipped below a threshold moisture content, the 'Transportable Moisture Limit' (TML)—as determined from one of three laboratory test methods specified within the IMSBC Code. However, there are recognised problems with these tests, namely that for any given sample tested: (i) the tests give variable TML results, particularly when conducted in different laboratories; and (ii) each of the three test methods give a different TML result. As the TML determines a cargo's transportability, these issues have commercial and safety implications. To address these issues for the two most popular test methods, the Flow Table Test (FTT) and Proctor/Fagerberg Test (PFT), tests were performed on sand and iron ore fines. For each test, major sources of variability were identified, along with issues inherent in the test methodology. From these findings, a number of recommendations and several modifications have been proposed that, if implemented within the IMSBC Code, will improve the tests' reliability. The results were then compared to establish whether the different TML results obtained from the two test methods could be aligned—this showed that calibration between the methods was not possible on the basis of their TML. The test method assessment also incorporated observations at TML test facilities, interviews with leading researchers and an assessment of archive research material. Whilst the findings related to variability will help improve the reliability of the TML test methods, the accuracy of these methods is still not understood (i.e. how well do they indicate the likelihood for an actual cargo to liquefy?). This thesis therefore gives consideration to the likely accuracy of the TML test methods. It presents an improved representation of the cargo shift problem, with a view of providing future researchers with a basis to further assess this accuracy and thus the suitability of the TML as a criteria to detect cargoes which may liquefy.

The problem of flow separation around islands is investigated using a dynamically adaptive finite element model to allow for resolution of the shear layers that form in the advent of separation. The changes in secondary circulation and vertical motion that occur in both attached and separated flows are documented, as is the degree of closure of the wake eddies. In the numerical experiments presented, the strongest motion always takes place at the sides of the idealised island, where flow curvature and shear act together to induce ascent. In contrast, it is the slower motion within the wake eddies that allow streamlines to extend from the bottom to the surface. We find no evidence for closure of the wake eddies. Rather, all of our separated experiments show that streamlines that pass through the eddies originate outside of the shear layers and frictional boundary layers on the upstream side of the idealised island. The numerical experiments demonstrate the potential for dynamically adaptive, unstructured meshes to resolve the separated shear layers that occur downstream of the idealised island, as well as the narrow boundary layers that form on the island itself. © 2010 Springer-Verlag.

Shallow offshore foundations, which achieve their stability through the foundation bearing on the seabed, can in most applications be idealised as large rigid circular footings subjected to vertical, horizontal and moment loading. A small strain linear-elastic perfectly-plastic three- dimensional finite element program was developed to analyse this combined loading problem. The selection of a suitable three-dimensional finite element for accurate and computationally efficient analysis was based on the element's ability to model incompressible soil conditions, using exact numerical integration. A new approach for evaluating element suitability is developed and is quantified in terms of the parameter free degrees-of-freedom (equal to the degrees-of-freedom minus the incompressibility constraints). The tetrahedron family of three- dimensional elements is found to be in general more suitable and computationally efficient than the Serendipity and Lagrangian cube families. The 20-node quadratic strain tetrahedron is adopted for all analyses presented in this thesis. For combined loading, most of the available elastic analytical solutions are for rigid circular footings placed at the surface. Finite element analyses reveal certain inadequacies in some of these solutions. The elastic numerical analyses also examine the effect of footing embedment for three cases of embedment geometry. This demonstrates that the increase in footing stiffness (reduction in displacement) due to embedment for horizontal and moment loading is developed at shallower depths, and has a greater magnitude, than for vertical loading. The stability of a rough rigid circular footing placed on the surface of an undrained clay is examined. Zero thickness interface elements, despite experiencing some numerical instability problems, were found to model the footing-soil interaction better than conventional continuum elements. However, application of interface elements to footings which can lose contact with the soil, under moment loading conditions, resulted in numerical instability of the solution. A simpler model was therefore used to define the shape of the three-dimensional failure envelope at high vertical load. Comparisons with the semi-empirical inclination factors of bearing capacity solutions are included.

This study presents a numerical investigation of ducted tidal turbines, employing three-dimensional Reynolds-averaged Navier-Stokes simulations. Bidirectional ducted turbines are modelled with and without aperture, referred to as ducted and open- centre turbines respectively. The work consists of two investigations. In the first, the turbine rotors are represented by actuator discs, a simplification which captures changes in linear momentum and thus the primary interaction of the turbine with the flow through and around the duct, while greatly reducing computational complexity. In the second investigation, the turbine rotors are represented through a CFD-integrated blade element momentum model, employing realistic rotor data, capturing swirl and blade drag in addition to the extraction of linear momentum. Both modelling techniques were employed to investigate the performances of bare, ducted, and open-centre turbines, relating these to the flow fields exhibited. For axial flow, substantial decreases in power generated by the ducted and open-centre turbines were found, relative to a bare turbine of equal total device diameter. For open-centre turbines, an increase in aperture size leads to a further reduction in power generated. Increased blockage was shown to positively affect the performance of all devices. Two further measures of performance were employed: power density, normalising the power by the rotor area, and basin efficiency, relating the power generated to the overall power removed from the flow. Moderate increases in power density can be achieved for the ducted and open-centre devices, while their basin efficiencies are of similar value to that of the bare turbine. For yawed inflow, the performance of the bare turbine decreases, whilst that of the ducted and open-centre turbines increases. This is due to an increased flow velocity following flow acceleration around the inlet lip of the duct and also an increase in effective blockage as ducts present greater projected frontal area when approached non-axially.

This study presents a numerical investigation of a generic horizontal axis cross-flow marine turbine. The numerical tool used is the commercial Computational Fluid Dynamics package ANSYS FLUENT 12.0. The numerical model, using the SST k-w turbulence model, is validated against static, dynamic pitching blade and rotating turbine data. The work embodies two main investigations. The first is concerned with the influence of turbine solidity (ratio of net blade chord to circumference) on turbine performance, and the second with the influence of blockage (ratio of device frontal area to channel crosssection area) and free surface deformation on the hydrodynamics of energy extraction in a constrained channel.Turbine solidity was investigated by simulating flows through two-, three- and four-bladed turbines, resulting in solidities of 0.019, 0.029 and 0.038, respectively. The investigation was conducted for two Reynolds numbers, Re = O(10^5) & O(10^6), to reflect laboratory and field scales. Increasing the number of blades from two to four led to an increase in the maximum power coefficient from 0.43 to 0.53 for the lower Re and from 0.49 to 0.56 for the higher Re computations. Furthermore, the power curve was found to shift to a lower range of tip speed ratios when increasing solidity.The effects of flow confinement and free surface deformation were investigated by simulating flows through a three-bladed turbine with solidity 0.125 at Re = O(10^6) for channels that resulted in cross-stream blockages of 12.5% to 50%. Increasing the blockage led to a substantial increase in the power and basin efficiency; when approximating the free surface as a rigid lid, the highest power coefficient and basin efficiency computed were 1.18 and 0.54, respectively. Comparisons between the corresponding rigid lid and free surface simulations, where Froude number, Fr = 0.082, rendered similar results at the lower blockages, but at the highest blockage an increase in power and basin efficiency of up to 7% for the free surface simulations over that achieved with a rigid lid boundary condition. For the free surface simulations with Fr = 0.082, the energy extraction resulted in a drop in water depth of up to 0.7%. An increase in Fr from 0.082 to 0.131 resulted in an increase maximum power of 3%, but a drop in basin efficiency of 21%.

This thesis describes experiments that were carried out using focused wave groups in the UK Coastal Research Facility (UKCRF). Considerable effort was put into calibrating the UKCRF to determine the relationship between the input signals sent to the paddles and the waves generated in the facility. Focused wave groups of various sizes and phases, based on NewWave theory were generated, and measurements were made of the resulting surface elevation data, water particle kinematics, wave runup and overtopping volumes. NewWave theory models the profile of extreme waves in a Gaussian (random) sea. The thesis describes the first time this model has been applied in the context of coastal wave transformation. A method for the separation of the underlying harmonic structure of a focused wave group is described and results presented. This technique has been used in relatively deep water but is shown to work successfully in the coastal zone until wave overturning. A method has been devised to provide a theoretical Stokes-like expansion of the free and bound waves to model the surface elevation and water particle kinematics of the focused wave groups. Satisfactory agreement is achieved between the theoretical predictions of UKCRF measurements. Suggestions are made for an improved model. The underlying harmonic structure of the focused wave groups is presented as stacked time histories that give insight into the wave transformation process from deep to shallow water. Particular attention is paid to the low frequency wave generated as the wave group interacts with the beach. This is compared to the low frequency wave that is generated by a solitary wave in the UKCRF. Runup and overtopping measurements are in reasonable agreement with predictions based on certain empirical formulae, but not others. These comparisons are useful in identifying those formulae able to predict runup and overtopping of extreme waves in the coastal zone.

An understanding of extreme waves is important in the design and analysis of offshore structures, such as oil and gas platforms. With the increase of interest in the shipping of LNG, the design of import and export terminals in coastal water of slowly varying intermediate depth requires accurate analysis of steep wave shoaling. In this thesis, data from laboratory experiments involving random wave simulations on very gentle slopes have been analysed in terms of a model of large wave events, and the results interpreted by observation of the shape and magnitude of the large wave events. The auto-correlation function of the free surface elevation time histories, called NewWave, has been calculated from the wave spectrum and shown to fit very well up to the point where waves start to break (when compared to the ‘linear’ surface elevation time history). It has been shown that NewWave is an appropriate model for the shape of the ‘linear’ part of large waves provided kd > 0.5. A Stokes-like expansion for NewWave analysis has been demonstrated to match the average shape of the largest waves, accounting for the dominant vertical asymmetry. Furthermore, an appropriate local wave period derived from NewWave has been inserted into a Miche-based limiting criterion, using the linear dispersion equation, to obtain estimates for the limiting wave height. Overall, the analysis confirms the Miche-type criterion applies to limiting wave height for waves passing over very mild bed slopes. A derivation of general Green-Naghdi (GN) theory, which incorporates non-linear terms in its formulation, is also presented. This approach satisfies the boundary conditions exactly and approximates the field equations. The derived 2-dimensional vertical GN Level 1 model, capable of simulating steep waves on varying water depth, is validated against solitary waves and their interactions, and solitary waves on varying water depth and gives good qualitative agreement against the KdV equation. The developed and validated numerical model is used to simulate focussed wave groups on both constant depth and gentle slope. In general, the behaviour of waves simulated by the numerical model is very similar to that observed in the experimental data. There is evidence of vertical asymmetry as the water depth is reduced, owing to the non-linearity. Although the main physics is still controlled by linear dispersion, the higher order harmonics become increasingly important for shoaling waves. The numerical results also show a slope-induced wave set-up that keeps on increasing in amplitude as the wave group travels on the gentle slope.

Tidal stream energy has gained attention as a source of predictable and renewable energy. Devices resembling underwater wind turbines, placed in fast tidal streams, have been proposed to extract this energy. Arrays of many such devices will need to be deployed to deliver a significant amount of energy to the electricity grid. One consequence of energy extraction is that the array provides a resistance to the tidal stream, which may change the local and far field hydrodynamics, which in turn affects the power available to the array. Array-scale hydrodynamic changes affect the flow presented to the devices, which in turn affects the total resistance the array provides to the flow. This thesis is concerned with the interactions between device, array, and the tidal stream resource, to better understand the power potential of turbine arrays. Linear momentum actuator disc theory is employed to describe the operation of an idealised turbine array partially spanning a wide channel. The model is comprised of two quasi-independent sub-models, an array-scale model, describing flow phenomena around the array, which provides the upstream boundary condition to the device-scale model, describing the flow around a device. The thrust applied by the array is the sum of the thrust applied by the devices, completing the sub-model coupling. The numerical simulation of arrays in depth-averaged simulations is then investigated using the two-scale concept developed in the analytic partial-array model. It is shown that the device-scale flow must be modelled with a sub-grid scale model in order to correctly describe the unresolved device-scale flow and hence estimate the power available to an idealised array. Turbulence modelling in depth-averaged simulations of turbine arrays is also discussed. Temporal variations in tidal amplitude and strength mean that generator capacity will need to be economically matched to the available resource. As device performance may consequently depart from the relationship derived in idealised models when power capping is employed, blade element momentum theory is modified to parameterise tidal turbine performance during power capping. The array-scale effect of power capping is studied in depth-averaged simulations, in which it is shown that a significant reduction in device thrust may occur during power capping, reducing the impact of energy extraction from the tidal stream.

The increased use of sonar surveys for both industrial and leisure activities has motivated the research for cost effective, automated processed for seabed classification. Seabed classification is essential for many fields including dredging, environmental studies, fisheries research, pipeline and cable route surveys, marine archaeology and automated underwater vehicles. The advancement in both sonar technology and sonar data storage has led to large quantities of sonar data being collected per survey. The challenge, however, is to derive relevant features that can summarise these large amounts of data and provide discrimination between several seabed types present in each survey.The main aim of this work is to classify sidescan bathymetric datasets. However, in most sidescan bathymetric surveys, only a few ground-truthed areas (if any) are available. Since sidescan ‘ground-truthed’ areas were also provided for this work, they were used to test feature extraction, selection and classification algorithms. Backscattering amplitude, after using bathymetric data to correct for variations, did not provide enough discrimination between sediment classes in this work which lead to the investigation of other features. The variation of backscattering amplitude at different scales corresponds to variations in both micro bathymetry and large scale bathymetry. A method that can derive multiscale features from signals was needed, and the wavelet method proved to be an efficient method of doing so. Wavelets are used for feature extraction in 1D sidescan bathymetry survey data and both the feature selection and classification stages are automated. The method is tested on areas of known types and in general, the features show good correlation with sediment types in both types of survey.The main disadvantage of this method, however, is that signal futures are calculated per swathe (or received signal). Thus, sediment boundaries within the same swathe are not detected. To solve this problem, information present in consecutive pings of data can be used, leading to 2-D feature extraction.Several textural classification methods are investigated for the segmentation of sidescan sonar images. The method includes 2D wavelets and Gabor filters. Effects of filter orientation filter scale and window size are observed in both cases, and validated on given sonar images.For sidescan bathymetric datasets, a novel method of classification using both sidescan images and depth maps is investigated. Backscattering amplitude and bathymetry images are both used for feature extraction. Features include amplitude-dependent features, textural features and bathymetric variation features. The method makes use of grab samples available in given areas of the survey for training the classifiers. Alternatively, clustering techniques are used to group the data. The results of applying the method on sidescan bathymetric surveys correlate with the grab samples available as well as the user-classified areas.An automatic method for sidescan bathymetric classification offers a cost effective approach to classify large areas of seabed with a fewer number of grab samples. This work sheds light on areas of feature extraction, selection and classification of sonar data.

The energy resource in the Bristol Channel is of national strategic significance to meet the future demand for low carbon energy, and is the single largest resource area for tidal energy in the UK. However, the complex tidal dynamics of the Bristol Channel are not yet fully understood, in particular the sensitivity of the tidal dynamics to changes such as those caused by energy extraction. This is the motivation for the development of a simplified two-dimensional model to simulate the tidal flows in Bristol Channel. The tidal resonance in Bristol Channel has been investigated by exciting the model with a single tidal component and a series of sensitivity tests have been carried out on tuneable model parameters.

This abstract uses numerical models of several candidate sites around Great Britain to consider the predictability and temporal variations in tidal stream power.

The Transverse Horizontal Axis Water Turbine (THAWT) has been proposed as a tidal device which can be easily scaled and requires fewer foundations, bearings seals and generators than a more conventional axial-flow device. The THAWT device is a horizontally deployed variant of the Darrieus cross-flow turbine, in which the blades can be oriented into a truss configuration to produce long, stiff multi-bay rotors.This thesis establishes and combines a set of numerical models, which predict the hydrodynamic and structural performance of the THAWT device, with sufficient confidence to assess the feasibility of such a device at a full scale installation and to optimise its performance.Tests of 1/20th scale experimental models of the THAWT device have demonstrated that the truss configured device is capable of producing power with an efficiency close to that of the parallel configured turbine. In addition, variations in the configuration of the scale models have indicated how several design parameters affect the hydrodynamic performance of the device.A two-dimensional Navier-Stokes blade element model has been developed, in which the THAWT device is represented using an actuator cylinder. The addition of a hydrostatic free surface deformation correction has resulted in a model which is capable of matching experimental results with sufficient fidelity and accuracy.Blade loads from the numerical hydrodynamic model have been applied to a beam finite element analysis, to predict the stresses induced in the hydrofoils of the THAWT device.The numerical and hydrodynamic models are combined with a Linear Channel Momentum model to predict the performance of the THAWT device at a full scale tidal location. The effect that the device has on the channel flow indicates how much energy is available for extraction and how this energy might be most efficiently obtained. When considering material fatigue the analysis suggests that the structural design considerations dominate over the hydrodynamic considerations.

A significant proportion of the world's population and physical assets are located in low lying coastal zones. Accurate prediction of wave induced run-up and overtopping of sea defences are important in defining the extent and severity of wave action, and in assessing risk to people and property from severe storms and tsunamis.This thesis describes a one-dimensional numerical model based on the Boussinesq equations of Madsen and Sorensen (1992) and the non-linear shallow water equations. The model is suitable for simulating propagation of weakly non-linear and weakly dispersive waves from intermediate to zero depth, such that any inundation and/or overtopping caused by the incoming waves is also calculated as part of the simulation. Wave breaking is approximated by locally switching to the non-linear shallow water equations, which can model broken waves as bores. A piston paddle wavemaker is incorporated into the model for complete reproduction of laboratory experiments. A domain mapping technique is used in the vicinity of the paddle to transform a time-varying domain into a fixed domain, so that the governing equations can be more readily solved.First, various aspects of the numerical model are verified against known analytical and newly derived semi-analytical solutions. The complete model is then validated with laboratory measurements of run-up and overtopping involving solitary waves. NewWave focused wave groups, which give the expected shape of extreme wave events in a linear random sea, are used for further validation. Simulations of experiments of wave group run-up on a plane beach yield very good agreement with the measured run-up distances and free surface time series. Wave-by-wave overtopping induced by focused wave groups is also successfully simulated with the model, with satisfactory agreement between the experimental and the predicted overtopping volumes. Repeated simulations, now driven by second order paddle displacement signals, give insight into second order error waves spuriously generated by using paddle signals derived from linear theory. Separation of harmonics reveals that the long error wave is significantly affecting the wave group shape and leading to enhanced runu-up distances and overtopping volumes. An extensive parameter study is carried out using the numerical model investigating the influence on wave group run-up of linear wave amplitude at focus, linear focus location, and wave group phase at focus. For a given amplitude, both the phase and the focus location significantly affect the wave group run-up. It is also found that the peak optimised run-up increases with the wave amplitude, but wave breaking becomes an inhibiting factor for larger waves. This methodology is proposed for extreme storm wave induced run-up analysis.

Wave drift responses at low frequencies and springing at high frequencies can both be modelled as second order processes involving terms proportional to the square of the wave elevation in irregular seas. Although the models of both these phenomena share many common features, springing introduces several additional complexities. These concern both the second order hydrodynamic analysis, and the probabilistic modelling of the resulting first plus second order response. This article highlights these in the context of springing, which is particularly relevant to the design of Tension Leg Platforms and offshore structures susceptible to lightly damped resonant responses, at periods in the range from around 2 to 5 seconds. The three closely linked aspects of hydrodynamics, spectral analysis and probabilistic modelling are considered, and sample numerical results discussed.

The problem of second-order water wave diffraction of an incident monochromatic wave field by an array of bottom-mounted circular cylinders is solved by a semi-analytical approach. The solution for the second-order potential is obtained by combining eigenfunction expansions with an integral representation. Unlike the indirect approach for second-order forces (Lighthill 1979; Molin 1979), this approach gives complete information about local flow characteristics (pressure, velocities, wave elevation, etc.) thus providing a basis for solving the third-order problem. The results obtained are compared with other published data, and new detailed results, useful for benchmarking purposes, are given. Finally the influences of wave incidence, cylinder radius and cylinder configuration are considered. This leads to the suggestion that there exists a near-trapping phenomenon for the second-order wave in an array of cylinders, at half the wave frequency at which the corresponding linear near-trapped mode occurs.

The interaction between fully nonlinear water waves in a wave tank and fixed or floating structures with vertical and flared side surfaces is investigated. The higher-order boundary element method is used to solve the mixed boundary value problem in an Eulerian formulation at each time step. The time stepping scheme updates the boundary conditions on the free water surface and body surface, based on a Lagrangian description. In order to increase the efficiency of the calculation, the domain decomposition technique is implemented, with continuity conditions enforced on the interface between adjacent subdomains by an iterative procedure. For the calculation of wave forces acting on structures, some auxiliary functions are used, instead of predicting the time derivative of the potential directly. By means of these auxiliary functions, the coupled fluid–structure interaction can be decoupled easily, so that the calculation of body motions becomes independent of the hydrodynamic force. In addition, mesh regridding using the Laplace smoothing technique and interpolation are applied on the free surface to avoid possible numerical instabilities. Numerical results are obtained for the wave interactions with vertical cylinders and flared structures at different wave numbers, in which the bodies are fixed, freely floating in one direction and totally freely floating, respectively. The effect of different flared shapes on the forces, wave run-up and body motions is also studied.

Tidal stream energy has gained attention as a source of predictable and renewable energy. Devices resembling underwater wind turbines, placed in fast tidal streams, have been proposed to extract this energy. Arrays of many such devices will need to be deployed to deliver a significant amount of energy to the electricity grid. One consequence of energy extraction is that the array provides a resistance to the tidal stream, which may change the local and far field hydrodynamics, which in turn affects the power available to the array. Array-scale hydrodynamic changes affect the flow presented to the devices, which in turn affects the total resistance the array provides to the flow. This thesis is concerned with the interactions between device, array, and the tidal stream resource, to better understand the power potential of turbine arrays. Linear momentum actuator disc theory is employed to describe the operation of an idealised turbine array partially spanning a wide channel. The model is comprised of two quasi-independent sub-models, an array-scale model, describing flow phenomena around the array, which provides the upstream boundary condition to the device-scale model, describing the flow around a device. The thrust applied by the array is the sum of the thrust applied by the devices, completing the sub-model coupling. The numerical simulation of arrays in depth-averaged simulations is then investigated using the two-scale concept developed in the analytic partial-array model. It is shown that the device-scale flow must be modelled with a sub-grid scale model in order to correctly describe the unresolved device-scale flow and hence estimate the power available to an idealised array. Turbulence modelling in depth-averaged simulations of turbine arrays is also discussed. Temporal variations in tidal amplitude and strength mean that generator capacity will need to be economically matched to the available resource. As device performance may consequently depart from the relationship derived in idealised models when power capping is employed, blade element momentum theory is modified to parameterise tidal turbine performance during power capping. The array-scale effect of power capping is studied in depth-averaged simulations, in which it is shown that a significant reduction in device thrust may occur during power capping, reducing the impact of energy extraction from the tidal stream.

A three dimensional computational model of a full scale axial flow tidal turbine has been used to investigate the effects of a range of realistic environmental conditions on turbine performance. The model, which is based on the Reynolds averaged Navier-Stokes equations, has been developed using the commercial flow solver ANSYS Fluent. A 1:30 scale tidal turbine is simulated in an open channel for comparison to existing experimental data. The rotor blades are directly resolved using a body-fitted, unstructured computational grid. Rotor motion is enabled through a sliding mesh interface between the rotor and the channel boundaries. Reasonably good agreement in thrust and power is observed. The computed performance curves are offset from the measured performance curves by a small increment in rotor speed. Subsequently, a full scale axial flow turbine is modelled in a variety of conditions representative of tidal channel flows. A parametric study is carried out to investigate the effects of flow shear, confinement and alignment on turbine performance, structural loading, and wake recovery. Mean power and thrust are found to be higher in sheared flow, relative to uniform flow of equivalent volumetric flow rate. Large fluctuations in blade thrust and torque occur in sheared flow as the blade passes through the high velocity freestream flow in the upper portion of the profile and the lower velocity flow near the channel bed. A stronger shear layer is formed around the upper portion of the wake in sheared flow, leading to enhanced wake mixing. Mean power and thrust are reduced when the turbine is simulated at a lower position in a sheared velocity profile. However, fluctuations in blade loading are increased due to the higher velocity gradient. The opposite effects are observed when the turbine operates at greater heights in sheared flow. Flow misalignment has a negative impact on mean rotor thrust and power, as well as on unsteady blade loading. Although the range of unsteady loading is not increased significantly, additional perturbations are introduced due to interactions between the blade and the nacelle. A deforming surface is introduced using the volume-of-fluid method. Linear wave theory is combined with the existing free surface model to develop an unsteady inflow boundary condition prescribing combined sheared flow and free surface waves. The relative effects of the sheared profile and wave-induced velocities on turbine loading are identified through frequency analysis. Rotor and blade load fluctuations are found to increase with wave height and wave length. In a separate study, the performance of bi-directional ducted tidal turbines is investigated through a parametric study of a range of duct profiles. A two dimensional axi-symmetric computational model is developed to compare the ducted geometries with an unducted device under consistent blockage conditions. The best-performing ducted device achieves a peak power coefficient of approximately 45% of that of the unducted device. Comparisons of streamtube area, velocity and pressure for the flow through the ducted device shows that the duct limits the pressure drop across the rotor and the mass flow through the rotor, resulting in lower device power.

Computational fluid dynamics (CFD) is used in this thesis to model wind and tidal stream turbines and to investigate tidal turbine fence performance. There are two primary objectives of this work. The first is to develop and validate an actuator line method for the simulation of wind and tidal turbines which applies the blade forces to the flow field without the need for a regularisation kernel. The second is to examine tidal fences using, in part, the newly developed actuator line method. A potential flow equivalence method for determining the relative velocity to the blade chord and flow angle at the rotor blades in the actuator line method is proposed and validated. Results for simulations using this method compare favourably with those from both experiments and alternative computational methods, although the present model’s results deviate from experimental results in the vicinity of the blade tips. A CFD-embedded blade element-momentum tool is used to design rotors for operation in infinitely wide tidal fences spanning a tidal channel. Rotors are designed for fences with several different blockage ratios, with those designed for high blockage flows having greater solidity than those designed for operation in fences with lower blockage. It is found that designing rotors for operational blockage conditions can significantly improve the power output achieved by a tidal fence. Improved power output for higher blockage conditions is achieved by the application of greater thrust to the flow. Actuator line simulations of short (up to 8 turbines) fences with varying intra-rotor spacing and number of rotors confirm that hydrodynamic performance of the rotors improves as the spacing is reduced and as rotors are added to a fence. The position of a rotor within the fence impacts its performance; rotors at the ends of a fence extract reduced power compared to those at the centre of the fence, particularly for tip speed ratios greater than the design tip speed ratio.

If some solid bulk cargoes such as concentrates, unprocessed nickel ores and iron ore fines contain sufficient moisture, then cargo compaction during a voyage may cause liquefaction. This has been known to result in major cargo displacement ('cargo shift'), causing the vessel to capsize. This has resulted in a number of fatalities. To provide context to this interdisciplinary problem, an overview of the key considerations is presented, which includes: (i) a summary of the wide range of factors that influence a cargo's likelihood to shift; (ii) the soil mechanics principles that can be used to explain cargo shift; and (iii) the regulatory controls that are used to prevent its occurrence during shipment. Under the requirements of international maritime regulation, primarily the International Maritime Bulk Solid Cargoes (IMSBC) Code, a cargo that may liquefy must be shipped below a threshold moisture content, the 'Transportable Moisture Limit' (TML)—as determined from one of three laboratory test methods specified within the IMSBC Code. However, there are recognised problems with these tests, namely that for any given sample tested: (i) the tests give variable TML results, particularly when conducted in different laboratories; and (ii) each of the three test methods give a different TML result. As the TML determines a cargo's transportability, these issues have commercial and safety implications. To address these issues for the two most popular test methods, the Flow Table Test (FTT) and Proctor/Fagerberg Test (PFT), tests were performed on sand and iron ore fines. For each test, major sources of variability were identified, along with issues inherent in the test methodology. From these findings, a number of recommendations and several modifications have been proposed that, if implemented within the IMSBC Code, will improve the tests' reliability. The results were then compared to establish whether the different TML results obtained from the two test methods could be aligned—this showed that calibration between the methods was not possible on the basis of their TML. The test method assessment also incorporated observations at TML test facilities, interviews with leading researchers and an assessment of archive research material. Whilst the findings related to variability will help improve the reliability of the TML test methods, the accuracy of these methods is still not understood (i.e. how well do they indicate the likelihood for an actual cargo to liquefy?). This thesis therefore gives consideration to the likely accuracy of the TML test methods. It presents an improved representation of the cargo shift problem, with a view of providing future researchers with a basis to further assess this accuracy and thus the suitability of the TML as a criteria to detect cargoes which may liquefy.

Trapped modes around a row of bottom-mounted vertical circular cylinders in a channel are examined. The cylinders are identical, and their axes equally spaced in a plane perpendicular to the channel walls. The analysis has been made by employing the multipole expansion method under the assumption of linear water wave theory. At least the same number of trapped modes is shown to exist as the number of cylinders for both Neumann and Dirichlet trapped modes, with the exception that for cylinders having large radius the mode corresponding to the Dirichlet trapped mode for one cylinder will disappear. Close similarities between the Dirichlet trapped modes around a row of cylinders in a channel and the near-resonant phenomenon in the wave diffraction around a long array of cylinders in the open sea are discussed. An analogy with a mass–spring oscillating system is also presented.

An understanding of extreme waves is important in the design and analysis of offshore structures, such as oil and gas platforms. With the increase of interest in the shipping of LNG, the design of import and export terminals in coastal water of slowly varying intermediate depth requires accurate analysis of steep wave shoaling. In this thesis, data from laboratory experiments involving random wave simulations on very gentle slopes have been analysed in terms of a model of large wave events, and the results interpreted by observation of the shape and magnitude of the large wave events. The auto-correlation function of the free surface elevation time histories, called NewWave, has been calculated from the wave spectrum and shown to fit very well up to the point where waves start to break (when compared to the ‘linear’ surface elevation time history). It has been shown that NewWave is an appropriate model for the shape of the ‘linear’ part of large waves provided kd > 0.5. A Stokes-like expansion for NewWave analysis has been demonstrated to match the average shape of the largest waves, accounting for the dominant vertical asymmetry. Furthermore, an appropriate local wave period derived from NewWave has been inserted into a Miche-based limiting criterion, using the linear dispersion equation, to obtain estimates for the limiting wave height. Overall, the analysis confirms the Miche-type criterion applies to limiting wave height for waves passing over very mild bed slopes. A derivation of general Green-Naghdi (GN) theory, which incorporates non-linear terms in its formulation, is also presented. This approach satisfies the boundary conditions exactly and approximates the field equations. The derived 2-dimensional vertical GN Level 1 model, capable of simulating steep waves on varying water depth, is validated against solitary waves and their interactions, and solitary waves on varying water depth and gives good qualitative agreement against the KdV equation. The developed and validated numerical model is used to simulate focussed wave groups on both constant depth and gentle slope. In general, the behaviour of waves simulated by the numerical model is very similar to that observed in the experimental data. There is evidence of vertical asymmetry as the water depth is reduced, owing to the non-linearity. Although the main physics is still controlled by linear dispersion, the higher order harmonics become increasingly important for shoaling waves. The numerical results also show a slope-induced wave set-up that keeps on increasing in amplitude as the wave group travels on the gentle slope.

This study presents a numerical investigation of ducted tidal turbines, employing three-dimensional Reynolds-averaged Navier-Stokes simulations. Bidirectional ducted turbines are modelled with and without aperture, referred to as ducted and open- centre turbines respectively. The work consists of two investigations. In the first, the turbine rotors are represented by actuator discs, a simplification which captures changes in linear momentum and thus the primary interaction of the turbine with the flow through and around the duct, while greatly reducing computational complexity. In the second investigation, the turbine rotors are represented through a CFD-integrated blade element momentum model, employing realistic rotor data, capturing swirl and blade drag in addition to the extraction of linear momentum. Both modelling techniques were employed to investigate the performances of bare, ducted, and open-centre turbines, relating these to the flow fields exhibited. For axial flow, substantial decreases in power generated by the ducted and open-centre turbines were found, relative to a bare turbine of equal total device diameter. For open-centre turbines, an increase in aperture size leads to a further reduction in power generated. Increased blockage was shown to positively affect the performance of all devices. Two further measures of performance were employed: power density, normalising the power by the rotor area, and basin efficiency, relating the power generated to the overall power removed from the flow. Moderate increases in power density can be achieved for the ducted and open-centre devices, while their basin efficiencies are of similar value to that of the bare turbine. For yawed inflow, the performance of the bare turbine decreases, whilst that of the ducted and open-centre turbines increases. This is due to an increased flow velocity following flow acceleration around the inlet lip of the duct and also an increase in effective blockage as ducts present greater projected frontal area when approached non-axially.

This thesis details the development, verification and validation of an unsteady unstructured high order (≥ 3) h/p Discontinuous Galerkin - Fourier solver for the incompressible Navier-Stokes equations on static and rotating meshes in two and three dimensions. This general purpose solver is used to provide insight into cross-flow (wind or tidal) turbine physical phenomena.Simulation of this type of turbine for renewable energy generation needs to account for the rotational motion of the blades with respect to the fixed environment. This rotational motion implies azimuthal changes in blade aero/hydro-dynamics that result in complex flow phenomena such as stalled flows, vortex shedding and blade-vortex interactions. Simulation of these flow features necessitates the use of a high order code exhibiting low numerical errors. This thesis presents the development of such a high order solver, which has been conceived and implemented from scratch by the author during his doctoral work.To account for the relative mesh motion, the incompressible Navier-Stokes equations are written in arbitrary Lagrangian-Eulerian form and a non-conformal Discontinuous Galerkin (DG) formulation (i.e. Symmetric Interior Penalty Galerkin) is used for spatial discretisation. The DG method, together with a novel sliding mesh technique, allows direct linking of rotating and static meshes through the numerical fluxes. This technique shows spectral accuracy and no degradation of temporal convergence rates if rotational motion is applied to a region of the mesh. In addition, analytical mappings are introduced to account for curved external boundaries representing circular shapes and NACA foils.To simulate 3D flows, the 2D DG solver is parallelised and extended using Fourier series. This extension allows for laminar and turbulent regimes to be simulated through Direct Numerical Simulation and Large Eddy Simulation (LES) type approaches. Two LES methodologies are proposed.Various 2D and 3D cases are presented for laminar and turbulent regimes. Among others, solutions for: Stokes flows, the Taylor vortex problem, flows around square and circular cylinders, flows around static and rotating NACA foils and flows through rotating cross-flow turbines, are presented.

Offshore structures are required to survive in extreme wave environments. Historically, the design of these offshore structures and vessels has relied on wave-tank experiments and linear theory. Today, with advances in computing power, it is becoming feasible to supplement these methods of analysis with fully nonlinear numerical simulation. This thesis is concerned with the development of an efficient method to perform this numerical modelling, in the context of potential flow theory. The interaction of a steep ocean wave with a floating body involves a moving free surface and a wide range of length scales. Attempts to reduce the size of the simulation domain cause problems with wave reflection from the domain edge and with the accurate creation of incident waves. A method of controlling the wave field around a structure is presented. The ability to effectively damp an outgoing wave in a very short distance is demonstrated. Steep incident waves are generated without the requirement for the wave to evolve over a large time or distance before interaction with the body. This enables a general wave-structure interaction problem to be modelled in a small tank, and behave as if it were surrounded by a large expanse of water. The suitability of the boundary element method for performing this modelling is analysed. Potential improvements are presented with respect to accuracy, robustness, and computational complexity. Evidence of third order diffraction is found for an FPSO model.

The Transverse Horizontal Axis Water Turbine (THAWT) has been proposed as a tidal device which can be easily scaled and requires fewer foundations, bearings seals and generators than a more conventional axial-flow device. The THAWT device is a horizontally deployed variant of the Darrieus cross-flow turbine, in which the blades can be oriented into a truss configuration to produce long, stiff multi-bay rotors.This thesis establishes and combines a set of numerical models, which predict the hydrodynamic and structural performance of the THAWT device, with sufficient confidence to assess the feasibility of such a device at a full scale installation and to optimise its performance.Tests of 1/20th scale experimental models of the THAWT device have demonstrated that the truss configured device is capable of producing power with an efficiency close to that of the parallel configured turbine. In addition, variations in the configuration of the scale models have indicated how several design parameters affect the hydrodynamic performance of the device.A two-dimensional Navier-Stokes blade element model has been developed, in which the THAWT device is represented using an actuator cylinder. The addition of a hydrostatic free surface deformation correction has resulted in a model which is capable of matching experimental results with sufficient fidelity and accuracy.Blade loads from the numerical hydrodynamic model have been applied to a beam finite element analysis, to predict the stresses induced in the hydrofoils of the THAWT device.The numerical and hydrodynamic models are combined with a Linear Channel Momentum model to predict the performance of the THAWT device at a full scale tidal location. The effect that the device has on the channel flow indicates how much energy is available for extraction and how this energy might be most efficiently obtained. When considering material fatigue the analysis suggests that the structural design considerations dominate over the hydrodynamic considerations.

A significant proportion of the world's population and physical assets are located in low lying coastal zones. Accurate prediction of wave induced run-up and overtopping of sea defences are important in defining the extent and severity of wave action, and in assessing risk to people and property from severe storms and tsunamis.This thesis describes a one-dimensional numerical model based on the Boussinesq equations of Madsen and Sorensen (1992) and the non-linear shallow water equations. The model is suitable for simulating propagation of weakly non-linear and weakly dispersive waves from intermediate to zero depth, such that any inundation and/or overtopping caused by the incoming waves is also calculated as part of the simulation. Wave breaking is approximated by locally switching to the non-linear shallow water equations, which can model broken waves as bores. A piston paddle wavemaker is incorporated into the model for complete reproduction of laboratory experiments. A domain mapping technique is used in the vicinity of the paddle to transform a time-varying domain into a fixed domain, so that the governing equations can be more readily solved.First, various aspects of the numerical model are verified against known analytical and newly derived semi-analytical solutions. The complete model is then validated with laboratory measurements of run-up and overtopping involving solitary waves. NewWave focused wave groups, which give the expected shape of extreme wave events in a linear random sea, are used for further validation. Simulations of experiments of wave group run-up on a plane beach yield very good agreement with the measured run-up distances and free surface time series. Wave-by-wave overtopping induced by focused wave groups is also successfully simulated with the model, with satisfactory agreement between the experimental and the predicted overtopping volumes. Repeated simulations, now driven by second order paddle displacement signals, give insight into second order error waves spuriously generated by using paddle signals derived from linear theory. Separation of harmonics reveals that the long error wave is significantly affecting the wave group shape and leading to enhanced runu-up distances and overtopping volumes. An extensive parameter study is carried out using the numerical model investigating the influence on wave group run-up of linear wave amplitude at focus, linear focus location, and wave group phase at focus. For a given amplitude, both the phase and the focus location significantly affect the wave group run-up. It is also found that the peak optimised run-up increases with the wave amplitude, but wave breaking becomes an inhibiting factor for larger waves. This methodology is proposed for extreme storm wave induced run-up analysis.

When traditional linearised theory is used to study free-surface flows past a surface-piercing object or over an obstruction in a stream, the geometry of the object is usually lost, having been assumed small in one or several of its dimensions. In order to preserve the nonlinear nature of the geometry, asymptotic expansions in the low-Froude or low-Bond limits can be derived, but here, the solution invariably predicts a waveless free-surface at every order. This is because the waves are in fact, exponentially small, and thus beyond-all-orders of regular asymptotics; their formation is a consequence of the divergence of the asymptotic series and the associated Stokes Phenomenon.In this thesis, we will apply exponential asymptotics to the study of two new problems involving nonlinear geometries. In the first, we examine the case of free-surface flow over a step including the effects of both gravity and surface tension. Here, we shall see that the availability of multiple singularities in the geometry, coupled with the interplay of gravitational and cohesive effects, leads to the discovery of a remarkable new set of solutions.In the second problem, we study the waves produced by bluff-bodied ships in low-Froude flows. We will derive the analytical form of the exponentially small waves for a wide range of hull geometries, including single-cornered and multi-cornered ships, and then provide comparisons with numerical computations. A particularly significant result is our confirmation of the thirty-year old conjecture by Vanden-Broeck & Tuck (1977) regarding the impossibility of waveless single-cornered ships.

Applying dimples on spherical bluff bodies have been found to be a major technique to control the surface flow and reduce the stream-wise force. However, not much work has been done on dimple applications on cylindrical Objects . This research project aims to investigate, firstly, whether the stream-wise force on a circular cylinder could be reduced by using dimples on its outer surface. Secondly, if it could, whether there is an optimum design that gives the biggest force reduction. Thirdly, what the roles of different approaches in this project might be? In order to find answers to these questions, experimental and computational methods have been used. One smooth and three dimpled cylinders have been tested. It was found that applying dimples on the outer surface of a cylinder is very effective in reducing stream-wise force. Drag coefficient significantly reduced by 17% to 48% depending on different Reynolds numbers. However, the dimpled design did not seem to be the key factor of the reduction, and the level of drag reduction was found to be insensitive to different dimpled patterns. Furthermore, it has also been found that CFD method would not be suitable as the research tool for this project because of its demanding requirement in computer power. Therefore, rig tests remained the primary method for this project. In addition, flow visualization technique was also used to show the wake region behind both the smooth and dimpled cylinders. It was clear that a smaller wake was associated with the dimpled case, which resulted in the pressure drag on the dimpled cylinder being smaller than for the smooth cylinder.