scholarly journals Numerical Method to Simulate Self-Propulsion of Aframax Tanker in Irregular Waves

2020 ◽  
Vol 2020 ◽  
pp. 1-10
Author(s):  
Gui-sheng Peng ◽  
Yang Gao ◽  
Wang Wen-hua ◽  
Lin Lin ◽  
Yi Huang

With the implement of ship energy efficiency design index (EEDI), computational fluid dynamics (CFD) technique has become an effective method to predict the ship performance and further guide the designers to optimize hull lines. However, due to the complexity of the propeller-hull interactions and the ship’s complex motions in waves, accurately predicting the speed-power performance of a self-propelled ship in actual seaway remains a challenge. In the present work, firstly, the resistance and self-propulsion experiments of Aframax model in waves are carried out at FORCE towing tank. Then, the CFD model and method are adopted to investigate the resistance and thrust under the conditions of regular and irregular waves in a three-dimensional numerical wave tank created by commercial software Star-CCM+. Therein, Reynolds-Averaged Navier–Stokes (RANS) equations and k-ε turbulent models were used for modeling the turbulent flow, and volume of fluid (VOF) method was applied to track the location and shape of transit-free surface. Based on the numerical method, the added resistance caused by regular waves was firstly investigated, and the self-propulsion of propeller in irregular waves was further performed. Furthermore, in order to simulate the rotation of the propeller, both the sliding mesh technique and overset mesh technique were discussed. Finally, compared with the experimental data, the numerical solutions have been validated, which shows potential to provide theoretical guidance and technical support for the self-propulsion performance of Aframax tanker in waves.

2007 ◽  
Author(s):  
Yutaka Masuyama ◽  
Yusuke Tahara ◽  
Toichi Fukasawa ◽  
Naotoshi Maeda

Database of full-scale three-dimensional sail shapes are presented with the aerodynamic coefficients for the upwind condition of IMS type sails. Three-dimensional shape data are used for the input of numerical calculations and the results are compared with the measured sail performance. The sail shapes and performance are measured using a sail dynamometer boat Fujin. The Fujin is a 34-foot LOA boat, in which load cells and charge coupled devices (CCD) cameras are installed to measure the sail forces and shapes simultaneously. The sailing conditions of the boat, such as boat speed, heel angle, wind speed, wind angle, and so on, are also measured. The tested sail configurations are as follows: mainsail with 130% jib, mainsail with 75% jib and mainsail alone. Sail shapes are measured at several height positions. The measured shape parameters are chord length, maximum draft, maximum draft position, entry angle at the luff and exit angle at the leech. From these parameters three-dimensional coordinates of the sails are calculated by interpolation. These three-dimensional coordinates are tabulated with the aerodynamic coefficients. Numerical calculations are performed using the measured sail shapes. The calculation methods are of two types; Reynolds-averaged Navier-Stokes (RANS)-based CFD and vortex lattice methods (VLM). A multi-block RANS-based CFD method was developed by one of the authors and is capable of predicting viscous flows and aerodynamic forces for complicated sail configuration for upwind as well as downwind conditions. Important features of the numerical method are summarized as follows: a Finite- Analytic scheme to discretize transport equations, a PISO type velocity-pressure coupling scheme, multi-block domain decomposition capability, and several choices of turbulence models depending on flows of interest. An automatic grid generation scheme is also included. Another calculation method, the vortex lattice method is also adopted. In this case, step-by-step calculations are conducted to attain the steady state of the sail in steady wind. Wake vortices are generated step-by-step, which flow in the direction of the local velocity vector. These calculated sail forces are compared with the measured one, and the validity of the numerical method is studied. The sail shape database and comparison with numerical calculations will provide a good benchmark for the sail performance analysis of the upwind condition of IMS type sails.


2021 ◽  
Author(s):  
Shan Wang ◽  
C. Guedes Soares

Abstract Three-dimensional effects on slamming loads predictions of a ship section are investigated numerically using the unsteady incompressible Reynolds-Average Navier-Stokes (RANS) equations and volume of fluid (VOF) method, which are implemented in interDyMFoam solver in open-source library OpenFoam. A convergence and uncertainty study is performed considering different resolutions and constant Courant number (CFL) following the ITTC guidelines. The numerical solutions are validated through comparisons of slamming loads and motions between the CFD simulations and the available experimental values. The total slamming force and slamming pressures on a 2D ship section and the 3D model are compared and discussed. Three-dimensional effects on the sectional force and the pressures are quantified both in transverse and longitudinal directions of the body considering various entry velocities. The non-dimensional pressure coefficient distribution on the 3D model is presented.


2006 ◽  
Vol 28 (3) ◽  
pp. 134-144
Author(s):  
Nguyen The Duc

The paper presents a numerical method to simulate two-phase turbulent cavitating flows in ducts of varying cross-section usually faced in engineering. The method is based on solution of two-phase Reynolds-averaged Navier-Stokes equations of two-phase mixture. The numerical method uses artificial compressibility algorithm extended to unsteady flows with dual-time technique. The discreted method employs an implicit, characteristic-based upwind differencing scheme in the curvilinear grid systems. Numerical simulation of an unsteady three-dimensional two-phase cavitating flow in a duct of varying cross-section with available experiment was performed. The unsteady important characteristics of the unsteady flow can be observed in results of numerical simulation. Comparison of predicted results with experimental data for time-averaged velocity and phase fraction are provided.


Author(s):  
Mustafa Koz ◽  
Serhat Yesilyurt

Microorganisms such as bacteria use their rotating helical flagella for propulsion speeds up to tens of tail lengths per second. The mechanism can be utilized for controlled pumping of liquids in microchannels. In this study, we aim to analyze the effects of control parameters such as axial span between helical rounds (wavelength), angular velocity of rotations (frequency), and the radius of the helix (amplitude) on the maximum time-averaged flow rate, maximum head, rate of energy transfer, and efficiency of the micropump. The analysis is based on simulations obtained from the three-dimensional time-dependent numerical model of the flow induced by the rotating spiral inside a rectangular-prism channel. The flow is governed by Navier-Stokes equations subject to continuity in time-varying domain due to moving boundaries of the spiral. Numerical solutions are obtained using a commercial finite-element package which uses arbitrary Lagrangian-Eulerian method for mesh deformations. Results are compared with asymptotic results obtained from the resistive-force-theory available in the literature.


2020 ◽  
Author(s):  
HaKun Jang ◽  
Celalettin Emre Ozdemir ◽  
Mayank Tyagi ◽  
Jun-Hong Liang

Abstract The purpose of this study is to numerically investigate the bed shear stress and near-bed mixing due to coherent vortex structures in the vicinity of a vertically wall-mounted circular cylinder subject to an imposed finite-depth oscillatory sinusoidal flow. Previous studies reveal that the Keulegan–Carpenter (KC) number influences the formation of lee-side wake vortex structures as well as the horseshoe vortex in front of a cylinder. Therefore, parametric studies in a moderately wide range of KC from 5 to 20 are numerically performed. In the present study, Direct Numerical Simulation (DNS) is conducted using the open-source software, OpenFOAM, that solves the three-dimensional unsteady incompressible Navier-Stokes equations using finite volume method. Nondimensional parameters used in the simulations are carefully chosen to represent the real physics. The numerical solutions are first validated using an analytical solution for the oscillating Stokes flow and the results are then systematically and quantitatively compared with the experimental measurements. The results show that the lee-side wake is significantly influenced by KC, and distinctive types of the lee-side wake are generated and classified based on KC. It is also found that both KC and the ratio of the thickness of the Stokes boundary layer to the water depth are heavily associated with the stability of the lee-side wake. In addition, the simulated size and lifespan of the horseshoe vortex agree well with the experimental data.


2011 ◽  
Vol 52-54 ◽  
pp. 267-272 ◽  
Author(s):  
Yong Hua Zhang ◽  
Jian Hui He ◽  
Guo Qing Zhang

This paper aims to understand influence of the obliquity of fin ray on its motion performance. An environment-friendly propulsion system mimicking undulating fins of stingray had been built. Investigations were presented by using three-dimensional unsteady Computational Fluid Dynamics (CFD) method. An unstructured, grid-based, unsteady Navier-Stokes solver with automatic adaptive remeshing was used to compute the unsteady flow around the fin through twenty complete cycles. The pressure distribution on fin surface was computed and integrated to provide fin forces which were decomposed into lift and thrust. Vortex contour maps of the fin with different obliquity of fin ray were displayed and compared. Finally, we draw a conclusion that the generated propulsion force of the biomimetic propulsor is gradually increase with the obliquity of the fin ray from 0 degree till a certain angle and then gradually decrease with the obliquity of the fin ray from the certain angle till 90 degree. The results provide valuable information for the optimization of robotic underwater undulating propulsor design.


2021 ◽  
Vol 143 (12) ◽  
Author(s):  
Sasuga Ito ◽  
Masato Furukawa ◽  
Kazutoyo Yamada ◽  
Kaito Manabe

Abstract Turbulence is one of the most important phenomena in fluid dynamics. Large eddy simulation (LES) generally allows us to analyze smaller eddies than when using simulations based on unsteady Reynolds-averaged Navier–Stokes equations (URANS). In addition, the numerical solutions of LES show good agreements with experiments and numerical solutions based on direct numerical simulation. URANS simulations are, however, frequently used in academia and industry because LES computations are much more expensive compared with URANS simulations. In this investigation, an optimization of unsolved coefficients of the k–ω two equations model is performed on the transonic flow around T106A low-pressure turbine cascade to improve the accuracy of turbulence prediction with URANS. For the optimization approach, two-dimensional URANS is combined with ensemble Kalman filter which is one of the data assimilation techniques. In the assimilation process, a time- and spanwise-averaged LES result is used as pseudo-experimental data. Three-dimensional URANS simulations are performed for the evaluation of the optimization effect. URANS simulations are also applied to a different turbine cascade flow for the evaluation of the robustness of the optimized coefficients. These URANS results confirmed that the optimized coefficients improve the accuracy of turbulence prediction.


1983 ◽  
Vol 105 (3) ◽  
pp. 639-645 ◽  
Author(s):  
Ching-Jen Chen ◽  
Young Hwan Yoon

Connective heat transfer for steady-state laminar flow in axisymmetric coordinates is considered. Numerical solutions for flow pattern and temperature distribution are obtained by the finite analytic numerical method applied to the Navier-Stokes equations expressed in terms of vorticity and stream function, and the energy equation. The finite analytic numerical method differs from other numerical methods in that it utilizes a local analytic solution in an element of the problem to construct the total numerical solution. Finite analytic solutions of vorticity, stream function, temperature, and heat transfer coefficients for flow with Reynolds numbers of 5, 100, 1000, and 2000, and Prandtl numbers of 0.1, 1.0, and 10.0 with uniform grid sizes, are reported for an axisymmetric pipe with a sudden expansion and contraction. The wall temperature is considered to be isothermal and differs from the inlet temperature. It is shown that the finite analytic is stable, converges rapidly, and simulates the convection of fluid flow accurately, since the local analytic solution is capable of simulating automatically the influence of skewed convection through the element boundary on the interior nodal values, thereby minimizing the false numerical diffusion.


Sign in / Sign up

Export Citation Format

Share Document