Computational Simulation of the Tethered Undersea Kites for Power Generation

Author(s):  
Amirmahdi Ghasemi ◽  
David J. Olinger ◽  
Gretar Tryggvason

The dynamic motion of tethered undersea kites (TUSK) is studied using numerical simulations. TUSK systems consist of a rigid-winged kite moving in an ocean current. The kite is connected by tethers to a platform on the ocean surface or anchored to the seabed. Hydrodynamic forces generated by the kite are transmitted through the tethers to a generator on the platform to produce electricity. TUSK systems are being considered as an alternative to marine turbines since the kite can move in high speed motions to increase power production compared to conventional marine turbines. The two-dimensional Navier-Stokes equations are solved on a regular structured grid that comprises the ocean current flow, and an immersed boundary method is used for the rigid kite. A two-step projection method along with Open Multi-Processing (OpenMP) is employed to solve the flow equations. The reel-out and reel-in velocities of the two tethers are adjusted to control the kite angle of attack and the resultant hydrodynamic forces. A baseline simulation was studied where a high net power output was achieved during successive kite power and retraction phases. System power output, vorticity flow fields, tether tensions, and hydrodynamic coefficients for the kite are determined. The power output results are in good agreement with established theoretical results for a kite moving in two dimensions.

2016 ◽  
Vol 138 (12) ◽  
Author(s):  
Amirmahdi Ghasemi ◽  
David J. Olinger ◽  
Gretar Tryggvason

The dynamic motion of tethered undersea kites (TUSK) is studied using numerical simulations. TUSK systems consist of a rigid winged-shaped kite moving in an ocean current. The kite is connected by tethers to a platform on the ocean surface or anchored to the seabed. Hydrodynamic forces generated by the kite are transmitted through the tethers to a generator on the platform to produce electricity. TUSK systems are being considered as an alternative to marine turbines since the kite can move at a high-speed, thereby increasing power production compared to conventional marine turbines. The two-dimensional Navier–Stokes equations are solved on a regular structured grid to resolve the ocean current flow, and a fictitious domain-immersed boundary method is used for the rigid kite. A projection method along with open multiprocessing (OpenMP) is employed to solve the flow equations. The reel-out and reel-in velocities of the two tethers are adjusted to control the kite angle of attack and the resultant hydrodynamic forces. A baseline simulation, where a high net power output was achieved during successive kite power and retraction phases, is examined in detail. The effects of different key design parameters in TUSK systems, such as the ratio of tether to current velocity, kite weight, current velocity, and the tether to kite chord length ratio, are then further studied. System power output, vorticity flow fields, tether tensions, and hydrodynamic coefficients for the kite are determined. The power output results are shown to be in good agreement with the established theoretical results for a kite moving in two dimensions.


Author(s):  
Amirmahdi Ghasemi ◽  
David J. Olinger ◽  
Gretar Tryggvason

In this paper, a numerical simulation of three dimensional motion of tether undersea kites (TUSK) for power generation is studied. TUSK systems consist of a rigid-winged kite, or glider, moving in an ocean current. One proposed concept uses a tethered kite which is connected by a flexible tether to a support structure with a generator on the ocean surface. The numerical simulation models the flow field in a three-dimensional domain near the rigid undersea kite wing by solving the full Navier-Stokes equations. A two-step projection method along with Open Multi-Processing (OpenMP) is employed to solve the flow equations. In order to track the rigid kite, an immersed boundary method is used. A NACA 0021 airfoil is used for the cross section shape of the kite, and the tension forces in the elastic tethers are modeled by a simple Hooke’s law. A grid refinement study has been carried out to ensure the independence of the numerical results on the grid mesh resolution. Also, the Reynolds number independency has been studied. PID control methods are used to adjust the kite pitch, roll and yaw angles during power (tether reel-out) and retraction (reel-in) phases to obtain desired kite trajectories. During the reel-out phase the kite moves in successive cross-current motions in a figure-8 pattern, the tether length increases and power is generated. During reel-in the kite motion is along the tether, and kite hydrodynamic forces are reduced so that net positive power is produced. Kite trajectories, hydrodynamic forces on the kite, kite tether tension and output power are determined and analyzed for a baseline TUSK simulation.


Author(s):  
Jiho You ◽  
Jinmo Lee ◽  
Donghyun You

A computational simulation methodology, which combines a computational fluid dynamics technique and a computational structural dynamics technique, is employed to design a deformable foil of which kinematics is inspired by the propulsive motion of a fin or a tail of fish and cetacean. The unsteady incompressible Navier-Stokes equations are solved using a second-order accurate finite-difference method and an immersed-boundary method to effectively impose boundary conditions on complex moving boundaries. A finite-element-based structural dynamics solver is employed to compute the deformation of the foil due to interaction with fluid. A phase angle between pitching and heaving motions as well as the flexibility of the foil, which is represented by the Youngs modulus are varied to find out how these factors affect the propulsion efficiency.


2017 ◽  
Vol 828 ◽  
pp. 5-32 ◽  
Author(s):  
T. Barker ◽  
J. M. N. T. Gray

In recent years considerable progress has been made in the continuum modelling of granular flows, in particular the $\unicode[STIX]{x1D707}(I)$-rheology, which links the local viscosity in a flow to the strain rate and pressure through the non-dimensional inertial number $I$. This formulation greatly benefits from its similarity to the incompressible Navier–Stokes equations as it allows many existing numerical methods to be used. Unfortunately, this system of equations is ill posed when the inertial number is too high or too low. The consequence of ill posedness is that the growth rate of small perturbations tends to infinity in the high wavenumber limit. Due to this, numerical solutions are grid dependent and cannot be taken as being physically realistic. In this paper changes to the functional form of the $\unicode[STIX]{x1D707}(I)$ curve are considered, in order to maximise the range of well-posed inertial numbers, while preserving the overall structure of the equations. It is found that when the inertial number is low there exist curves for which the equations are guaranteed to be well posed. However when the inertial number is very large the equations are found to be ill posed regardless of the functional dependence of $\unicode[STIX]{x1D707}$ on $I$. A new $\unicode[STIX]{x1D707}(I)$ curve, which is inspired by the analysis of the governing equations and by experimental data, is proposed here. In order to test this regularised rheology, transient granular flows on inclined planes are studied. It is found that simulations of flows, which show signs of ill posedness with unregularised models, are numerically stable and match key experimental observations when the regularised model is used. This paper details two-dimensional transient computations of decelerating flows where the inertial number tends to zero, high-speed flows that have large inertial numbers, and flows which develop into granular rollwaves. This is the first time that granular rollwaves have been simulated in two dimensions, which represents a major step towards the simulation of other complex granular flows.


Author(s):  
Amirmahdi Ghasemi ◽  
David J. Olinger ◽  
Gretar Tryggvason

In this paper, a numerical simulation of tether undersea kites (TUSK) used for power generation is undertaken. The effect of varying key design parameters in these systems is studied. TUSK systems consist of a rigid-winged kite, or glider, moving in an ocean current. One proposed TUSK concept uses a tethered kite which is connected by a flexible tether to a support structure with a generator on a surface buoy. The numerical simulation models the flow field in a three-dimensional domain near the rigid undersea kite wing by solving the full Navier-Stokes equations. A moving computational domain method is used to reduce the computational run times. A second-order corrector-predictor method, along with Open Multi-Processing (OpenMP), is employed to solve the flow equations. In order to track the rigid kite, which is a rectangular planform wing with a NACA 0021 airfoil, an immersed boundary method is used. The tension force in the elastic tether is modeled by a simple Hooke’s law, and the effect of tether damping is added. PID control methods are used to adjust the kite pitch, roll and yaw angles during power (tether reel-out) and retraction (reel-in) phases to obtain the desired kite trajectories. During the reel-out phase the kite moves in successive cross-current motions in a figure-8 pattern, the tether length increases and power is generated. During reel-in the kite motion is along the tether, and kite hydrodynamic forces are reduced so that net positive power is produced. The effects of different key design parameters in TUSK systems, such as the ratio of tether to current velocity, and tether retraction velocity, are then further studied. System power output, kite trajectories, and vorticity flow fields for the kite are also determined.


Fluids ◽  
2019 ◽  
Vol 4 (3) ◽  
pp. 164 ◽  
Author(s):  
Bourantas ◽  
Zwick ◽  
Joldes ◽  
Loukopoulos ◽  
Tavner ◽  
...  

We present a strong form, meshless point collocation explicit solver for the numerical solution of the transient, incompressible, viscous Navier-Stokes (N-S) equations in two dimensions. We numerically solve the governing flow equations in their stream function-vorticity formulation. We use a uniform Cartesian embedded grid to represent the flow domain. We discretize the governing equations using the Meshless Point Collocation (MPC) method. We compute the spatial derivatives that appear in the governing flow equations, using a novel interpolation meshless scheme, the Discretization Corrected Particle Strength Exchange (DC PSE). We verify the accuracy of the numerical scheme for commonly used benchmark problems including lid-driven cavity flow, flow over a backward-facing step and unbounded flow past a cylinder. We have examined the applicability of the proposed scheme by considering flow cases with complex geometries, such as flow in a duct with cylindrical obstacles, flow in a bifurcated geometry, and flow past complex-shaped obstacles. Our method offers high accuracy and excellent computational efficiency as demonstrated by the verification examples, while maintaining a stable time step comparable to that used in unconditionally stable implicit methods. We estimate the stable time step using the Gershgorin circle theorem. The stable time step can be increased through the increase of the support domain of the weight function used in the DC PSE method.


Water ◽  
2021 ◽  
Vol 13 (7) ◽  
pp. 902
Author(s):  
Zhi Li ◽  
Ben R. Hodges

A new high-performance numerical model (Frehg) is developed to simulate water flow in shallow coastal wetlands. Frehg solves the 2D depth-integrated, hydrostatic, Navier–Stokes equations (i.e., shallow-water equations) in the surface domain and the 3D variably-saturated Richards equation in the subsurface domain. The two domains are asynchronously coupled to model surface-subsurface exchange. The Frehg model is applied to evaluate model sensitivity to a variety of simplifications that are commonly adopted for shallow wetland models, especially the use of the diffusive wave approximation in place of the traditional Saint-Venant equations for surface flow. The results suggest that a dynamic model for momentum is preferred over diffusive wave model for shallow coastal wetlands and marshes because the latter fails to capture flow unsteadiness. Under the combined effects of evaporation and wetting/drying, using diffusive wave model leads to discrepancies in modeled surface-subsurface exchange flux in the intertidal zone where strong exchange processes occur. It indicates shallow wetland models should be built with (i) dynamic surface flow equations that capture the timing of inundation, (ii) complex topographic features that render accurate spatial extent of inundation, and (iii) variably-saturated subsurface flow solver that is capable of modeling moisture change in the subsurface due to evaporation and infiltration.


Author(s):  
Hassan F Ahmed ◽  
Hamayun Farooq ◽  
Imran Akhtar ◽  
Zafar Bangash

In this article, we introduce a machine learning–based reduced-order modeling (ML-ROM) framework through the integration of proper orthogonal decomposition (POD) and deep neural networks (DNNs), in addition to long short-term memory (LSTM) networks. The DNN is utilized to upscale POD temporal coefficients and their respective spatial modes to account for the dynamics represented by the truncated modes. In the second part of the algorithm, temporal evolution of the POD coefficients is obtained by recursively predicting their future states using an LSTM network. The proposed model (ML-ROM) is tested for flow past a circular cylinder characterized by the Navier–Stokes equations. We perform pressure mode decomposition analysis on the flow data using both POD and ML-ROM to predict hydrodynamic forces and demonstrate the accuracy of the proposed strategy for modeling lift and drag coefficients.


1991 ◽  
Vol 113 (4) ◽  
pp. 608-616 ◽  
Author(s):  
H. M. Jang ◽  
J. A. Ekaterinaris ◽  
M. F. Platzer ◽  
T. Cebeci

Two methods are described for calculating pressure distributions and boundary layers on blades subjected to low Reynolds numbers and ramp-type motion. The first is based on an interactive scheme in which the inviscid flow is computed by a panel method and the boundary layer flow by an inverse method that makes use of the Hilbert integral to couple the solutions of the inviscid and viscous flow equations. The second method is based on the solution of the compressible Navier–Stokes equations with an embedded grid technique that permits accurate calculation of boundary layer flows. Studies for the Eppler-387 and NACA-0012 airfoils indicate that both methods can be used to calculate the behavior of unsteady blade boundary layers at low Reynolds numbers provided that the location of transition is computed with the en method and the transitional region is modeled properly.


1994 ◽  
Vol 47 (6S) ◽  
pp. S3-S13 ◽  
Author(s):  
Parviz Moin ◽  
Thomas Bewley

A brief review of current approaches to active feedback control of the fluctuations arising in turbulent flows is presented, emphasizing the mathematical techniques involved. Active feedback control schemes are categorized and compared by examining the extent to which they are based on the governing flow equations. These schemes are broken down into the following categories: adaptive schemes, schemes based on heuristic physical arguments, schemes based on a dynamical systems approach, and schemes based on optimal control theory applied directly to the Navier-Stokes equations. Recent advances in methods of implementing small scale flow control ideas are also reviewed.


Sign in / Sign up

Export Citation Format

Share Document