scholarly journals A Nonlinear Computational Model for Floating Wind Turbines

Author(s):  
Ali Nematbakhsh ◽  
David J. Olinger ◽  
Gretar Tryggvason

The dynamic motion of floating wind turbines is studied using computational simulations. The full three-dimensional Navier-Stokes equations are solved on a regular structured grid, using a level set method for the free surface and modified immersed boundary method to model the turbine platform. The tethers, the tower, the nacelle and the rotor weight are include using reduced order dynamic models, resulting in an efficient numerical approach. Wind is modeled as a constant thrust force. Other aerodynamic loading, rotor gyroscopic effects, and aeroelastic effects are not considered in the current study. The response of a tension leg platform to moderate amplitude waves is examined. By using the current approach, nearly all the nonlinear and viscose effects can be considered while keeping the computational cost reasonable. The model is applied to a Tension Leg Platform (TLP) consisting of a ballasted cylindrical tank.

2013 ◽  
Vol 135 (12) ◽  
Author(s):  
Ali Nematbakhsh ◽  
David J. Olinger ◽  
Gretar Tryggvason

The dynamic motion of floating wind turbines is studied using numerical simulations. The full three-dimensional Navier–Stokes equations are solved on a regular structured grid using a level set method for the free surface and an immersed boundary method for the turbine platform. The tethers, the tower, the nacelle, and the rotor weight are included using reduced-order dynamic models, resulting in an efficient numerical approach that can handle nearly all the nonlinear hydrodynamic forces on the platform, while imposing no limitation on the platform motion. Wind speed is assumed constant, and rotor gyroscopic effects are accounted for. Other aerodynamic loadings and aeroelastic effects are not considered. Several tests, including comparison with other numerical, experimental, and grid study tests, have been done to validate and verify the numerical approach. The response of a tension leg platform (TLP) to different amplitude waves is examined, and for large waves, a nonlinear trend is seen. The nonlinearity limits the motion and shows that the linear assumption will lead to overprediction of the TLP response. Studying the flow field behind the TLP for moderate amplitude waves shows vortices during the transient response of the platform but not at the steady state, probably due to the small Keulegan–Carpenter number. The effects of changing the platform shape are considered, and finally, the nonlinear response of the platform to a large amplitude wave leading to slacking of the tethers is simulated.


Computation ◽  
2018 ◽  
Vol 6 (3) ◽  
pp. 50
Author(s):  
Jonatas Borges ◽  
Marcos Lourenço ◽  
Elie Padilla ◽  
Christopher Micallef

The immersed boundary method has attracted considerable interest in the last few years. The method is a computational cheap alternative to represent the boundaries of a geometrically complex body, while using a cartesian mesh, by adding a force term in the momentum equation. The advantage of this is that bodies of any arbitrary shape can be added without grid restructuring, a procedure which is often time-consuming. Furthermore, multiple bodies may be simulated, and relative motion of those bodies may be accomplished at reasonable computational cost. The numerical platform in development has a parallel distributed-memory implementation to solve the Navier-Stokes equations. The Finite Volume Method is used in the spatial discretization where the diffusive terms are approximated by the central difference method. The temporal discretization is accomplished using the Adams-Bashforth method. Both temporal and spatial discretizations are second-order accurate. The Velocity-pressure coupling is done using the fractional-step method of two steps. The present work applies the immersed boundary method to simulate a Newtonian laminar flow through a three-dimensional sudden contraction. Results are compared to published literature. Flow patterns upstream and downstream of the contraction region are analysed at various Reynolds number in the range 44 ≤ R e D ≤ 993 for the large tube and 87 ≤ R e D ≤ 1956 for the small tube, considerating a contraction ratio of β = 1 . 97 . Comparison between numerical and experimental velocity profiles has shown good agreement.


2021 ◽  
Author(s):  
Junshi Wang ◽  
Vadim Pavlov ◽  
Zhipeng Lou ◽  
Haibo Dong

Abstract Dolphins are known for their outstanding swimming performance. However, the difference in flow physics at different speeds remains elusive. In this work, the underlying mechanisms of dolphin swimming at three speeds, 2 m/s, 5 m/s, and 8 m/s, are explored using a combined experimental and numerical approach. Using the scanned CAD model of the Atlantic white-sided dolphin (Lagenorhynchus acutus) and virtual skeleton-based surface reconstruction method, a three-dimensional high-fidelity computational model is obtained with time-varying kinematics. A sharp-interface immersed-boundary-method (IBM) based direct numerical simulation (DNS) solver is employed to calculate the corresponding thrust production, wake structure, and surface pressure at different swimming speeds. It is found that the fluke keeps its effective angle of attack at high values for about 60% of each stroke. The total pressure force coefficient along the x-axis converges as the speed increase. The flow and surface pressure analysis both show considerable differences between lower (2 m/s) and higher (5 m/s and 8 m/s) speeds. The results from this work help to bring new insight into understanding the force generation mechanisms of the highly efficient dolphin swimming and offer potential suggestions to the future designs of unmanned underwater vehicles.


2006 ◽  
Vol 128 (6) ◽  
pp. 1394-1399 ◽  
Author(s):  
Donghyun You ◽  
Meng Wang ◽  
Rajat Mittal ◽  
Parviz Moin

A novel structured grid approach which provides an efficient way of treating a class of complex geometries is proposed. The incompressible Navier-Stokes equations are formulated in a two-dimensional, generalized curvilinear coordinate system complemented by a third quasi-curvilinear coordinate. By keeping all two-dimensional planes defined by constant third coordinate values parallel to one another, the proposed approach significantly reduces the memory requirement in fully three-dimensional geometries, and makes the computation more cost effective. The formulation can be easily adapted to an existing flow solver based on a two-dimensional generalized coordinate system coupled with a Cartesian third direction, with only a small increase in computational cost. The feasibility and efficiency of the present method have been assessed in a simulation of flow over a tapered cylinder.


1986 ◽  
Vol 163 ◽  
pp. 257-282 ◽  
Author(s):  
Philip Hall ◽  
Mujeeb R. Malik

The instability of a three-dimensional attachment-line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite-amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time-dependent Navier–Stokes equations for the attachment-line flow have been solved using a Fourier–Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite-amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment-line boundary layer is also investigated.


Author(s):  
Maurizio Collu ◽  
Michael Borg ◽  
Andrew Shires ◽  
Feargal P. Brennan

In the present article, progress on the development of an aero-hydro-servo-elastic coupled model of dynamics for floating Vertical Axis Wind Turbines (VAWTs) is presented, called FloVAWT (Floating Vertical Axis Wind Turbine). Aerodynamics is based on Paraschivoiu’s Double-Multiple Streamtube (DMST) model [1] [2], relying on blade element momentum (BEM) theory, but also taking into account three-dimensional effects, dynamic stall, and unsteady wind profiles and platform motions. Hydrodynamics is modelled with a time domain seakeeping model [3], based on hydrodynamic coefficients estimated with a frequency analysis potential method. In this first phase of the research program, the system is considered a rigid body. The mooring system is represented through a user defined force-displacement relationship. Due to the lack of experimental data on offshore floating VAWTs, the model has initially been validated by taking each module separately and comparing it against known experimental data, showing good agreement. The capabilities of the program are illustrated through a case study, giving an insight on the relative importance of aerodynamics loads and gyroscopic effects with respect to hydrodynamic load effects.


Author(s):  
Yannis Kallinderis ◽  
Hyung Taek Ahn

Numerical prediction of vortex-induced vibrations requires employment of the unsteady Navier-Stokes equations. Current Navier-Stokes solvers are quite expensive for three-dimensional flow-structure applications. Acceptance of Computational Fluid Dynamics as a design tool for the offshore industry requires improvements to current CFD methods in order to address the following important issues: (i) stability and computation cost of the numerical simulation process, (ii) restriction on the size of the allowable time-step due to the coupling of the flow and structure solution processes, (iii) excessive number of computational elements for 3-D applications, and (iv) accuracy and computational cost of turbulence models used for high Reynolds number flow. The above four problems are addressed via a new numerical method which employs strong coupling between the flow and the structure solutions. Special coupling is also employed between the Reynolds-averaged Navier-Stokes equations and the Spalart-Allmaras turbulence model. An element-type independent spatial discretization scheme is also presented which can handle general hybrid meshes consisting of hexahedra, prisms, pyramids, and tetrahedral.


Author(s):  
Junshi Wang ◽  
Huy Tran ◽  
Martha Christino ◽  
Carl White ◽  
Joseph Zhu ◽  
...  

Abstract A combined experimental and numerical approach is employed to study the hydrodynamic performance and characterize the flow features of thunniform swimming by using a tuna-inspired underwater vehicle in forward swimming. The three-dimensional, time-dependent kinematics of the body-fin system of the underwater vehicle is obtained via a stereo-videographic technique. A high-fidelity computational model is then directly reconstructed based on the experimental data. A sharp-interface immersed-boundary-method (IBM) based incompressible flow solver is employed to compute the flow. The primary objective of the computational effort is to quantify the thrust performance of the model. The body kinematics and hydrodynamic performances are quantified and the dynamics of the vortex wake are analyzed. Results have shown significant leading-edge vortex at the caudal fin and unique vortex ring structures in the wake. The results from this work help to bring insight into understanding the thrust producing mechanism of thunniform swimming and to provide potential suggestions in improving the hydrodynamic performance of swimming underwater vehicles.


1998 ◽  
Vol 185 ◽  
pp. 217-218
Author(s):  
Kwing L. Chan ◽  
Y.C. Kim

Currently, the most successful direct simulation of the solar granules (and the convection/radiation transition layer) is the three-dimensional (3D) model computed by Stein and Nordlund (1989). So far, there is no other similar 3D models available for comparison [however, see Ludwig et al. (1997) for a recent 2D calculation]. We are developing an alternative numerical approach to simulate the 3D radiation hydrodynamics of this layer. In this approach, the Eddington approximation is used to handle the radiation rather than solving the radiative transfer equations along rays, and the ADISM method (Chan and Wolff 1982) which solves the Navier Stokes equations in conservative forms is used to speed up the thermal relaxation of the fluid layer. We are in the process of testing the numerical accuracy of the codes. This paper summarizes the results of a test that illustrate the effects of vertical space resolution on the mean profiles of some important quantities.


2017 ◽  
Vol 21 (5) ◽  
pp. 1231-1257 ◽  
Author(s):  
Jian Cheng ◽  
Xiaodong Liu ◽  
Tiegang Liu ◽  
Hong Luo

AbstractA parallel, high-order direct Discontinuous Galerkin (DDG) method has been developed for solving the three dimensional compressible Navier-Stokes equations on 3D hybrid grids. The most distinguishing and attractive feature of DDG method lies in its simplicity in formulation and efficiency in computational cost. The formulation of the DDG discretization for 3D Navier-Stokes equations is detailed studied and the definition of characteristic length is also carefully examined and evaluated based on 3D hybrid grids. Accuracy studies are performed to numerically verify the order of accuracy using flow problems with analytical solutions. The capability in handling curved boundary geometry is also demonstrated. Furthermore, an SPMD (single program, multiple data) programming paradigm based on MPI is proposed to achieve parallelism. The numerical results obtained indicate that the DDG method can achieve the designed order of accuracy and is able to deliver comparable results as the widely used BR2 scheme, clearly demonstrating that the DDG method provides an attractive alternative for solving the 3D compressible Navier-Stokes equations.


Sign in / Sign up

Export Citation Format

Share Document