scholarly journals On the Calculation of Propulsive Characteristics of a Bulk-Carrier Moving in Head Seas

2020 ◽  
Vol 8 (10) ◽  
pp. 786
Author(s):  
S. Polyzos ◽  
G. Tzabiras

The present work describes a simplified Computational Fluid Dynamics (CFD) approach in order to calculate the propulsive performance of a ship moving at steady forward speed in head seas. The proposed method combines experimental data concerning the added resistance at model scale with full scale Reynolds Averages Navier–Stokes (RANS) computations, using an in-house solver. In order to simulate the propeller performance, the actuator disk concept is employed. The propeller thrust is calculated in the time domain, assuming that the total resistance of the ship is the sum of the still water resistance and the added component derived by the towing tank data. The unsteady RANS equations are solved until self-propulsion is achieved at a given time step. Then, the computed values of both the flow rate through the propeller and the thrust are stored and, after the end of the examined time period, they are processed for calculating the variation of Shaft Horsepower (SHP) and RPM of the ship’s engine. The method is applied for a bulk carrier which has been tested in model scale at the towing tank of the Laboratory for Ship and Marine Hydrodynamics (LSMH) of the National Technical University of Athens (NTUA).

2020 ◽  
Vol 8 (2) ◽  
pp. 141 ◽  
Author(s):  
Ville Viitanen ◽  
Timo Siikonen ◽  
Antonio Sánchez-Caja

In this paper, we conducted numerical simulations to investigate single and two-phase flows around marine propellers in open-water conditions at different Reynolds number regimes. The simulations were carried out using a homogeneous compressible two-phase flow model with RANS and hybrid RANS/LES turbulence modeling approaches. Transition was accounted for in the model-scale simulations by employing an LCTM transition model. In model scale, also an anisotropic RANS model was utilized. We investigated two types of marine propellers: a conventional and a tip-loaded one. We compared the results of the simulations to experimental results in terms of global propeller performance and cavitation observations. The propeller cavitation, near-blade flow phenomena, and propeller wake flow characteristics were investigated in model- and full-scale conditions. A grid and time step sensitivity studies were carried out with respect to the propeller performance and cavitation characteristics. The model-scale propeller performance and the cavitation patterns were captured well with the numerical simulations, with little difference between the utilized turbulence models. The global propeller performance and the cavitation patterns were similar between the model- and full-scale simulations. A tendency of increased cavitation extent was observed as the Reynolds number increases. At the same time, greater dissipation of the cavitating tip vortex was noted in the full-scale conditions.


Author(s):  
Karl W. Schulz ◽  
Trond S. Meling

A multi-strip numerical method, combining solution of the incompressible Reynolds Averaged Navier-Stokes (RANS) equations with a finite-element structural dynamics response, has been developed to analyze the flow-structure interaction of long, flexible risers. This solution methodology combines a number of individual hydrodynamic simulations corresponding to individual axial strips along the riser section with a full 3D structural analysis to predict overall VIV loads and displacements. The hydrodynamic loading for each riser strip is derived from a 2D finite-volume discretization of the governing RANS equations which is applicable to both single and multiple riser configurations. The entire flow-structure solution procedure is carried out in the time domain via a loose coupling strategy, such that the hydrodynamic loads from each riser strip are summed to obtain the overall loading along the span of each riser. This loading is then used to integrate forward a single time-step in the riser equations of motion to obtain an updated riser displacement profile. Closure of the coupled flow-structure method is achieved by updating the riser displacements for each of the corresponding hydrodynamic strips in the next time-step integration. The developed multi-strip method is applied to a single bare riser subjected to both uniform and shear current profiles. The flow conditions and riser configuration were chosen to match the Marintek rotating rig experiments, and comparisons between experimental and numerical results are presented for several flow configurations and axial tensions. In addition, a parametric study is presented using 16, 32, and 64 hydrodynamic strips for a given flow configuration to ascertain the sensitivity of the results to the number of strips chosen.


2010 ◽  
Vol 132 (4) ◽  
Author(s):  
Jixian Yao ◽  
Steven E. Gorrell ◽  
Aspi R. Wadia

Part I of this paper validated the ability of the unsteady Reynolds-Averaged Navier-Stokes (RANS) solver PTURBO to accurately simulate distortion transfer and generation through selected blade rows of two multistage fans. In this part, unsteady RANS calculations were successfully applied to predict the 1/rev inlet total pressure distortion transfer in the entirety of two differently designed multistage fans. This paper demonstrates that high-fidelity computational fluid dynamics (CFD) can be used early in the design process for verification purposes before hardware is built and can be used to reduce the number of distortion tests, hence reducing engine development cost. The unsteady RANS code PTURBO demonstrated remarkable agreement with the data, accurately capturing both the magnitude and the profile of total pressure and total temperature measurements. Detailed analysis of the flow physics identified from the CFD results has led to a thorough understanding of the total temperature distortion generation and transfer mechanism, especially for the spatial phase difference of total pressure and total temperature profiles. The analysis illustrates that the static parameters are more revealing than their stagnation counterpart and that pressure and temperature rise are more revealing while the pressure and temperature ratio could be misleading. The last stage is effectively throttled by the inlet distortion even though the overall engine throttle remains unchanged. The total temperature distortion generally grows as flow passes through the fan stages.


Author(s):  
Tiziano Ghisu ◽  
Francesco Cambuli ◽  
Pierpaolo Puddu ◽  
Irene Virdis ◽  
Mario Carta

The work by Hu and Li (2018, “Unsteady RANS Simulations of Wells Turbine Under Transient Flow Conditions,” ASME J. Offshore Mech. Arct. Eng., 140(1), p. 011901) presents the numerical simulation of a high-solidity Wells turbine by means of a computational fluid dynamics (CFD) (Reynolds-averaged Navier–Stokes (RANS)) approach. A key aspect highlighted by the authors is the presence of a hysteretic loop in the machine's performance curves, due (according to their explanation) to the interaction of vortices shed by the blade with the blade circulation, which is responsible for the aerodynamic forces. It is our opinion that this work contains some serious errors that invalidate the results. In this brief discussion, we aim to demonstrate how the hysteresis found and discussed by the authors should not be present in the turbine analyzed in Hu and Li (2018, “Unsteady RANS Simulations of Wells Turbine Under Transient Flow Conditions,” ASME J. Offshore Mech. Arct. Eng., 140(1), p. 011901), and it is unlikely to be present in any Wells turbine. The fact that Hu and Li find hysteresis in their simulations is most likely caused by numerical errors due to an insufficient temporal discretization. This and other inaccuracies could have been avoided with a more careful consideration of the available literature.


Author(s):  
Leo M. González ◽  
Alvaro Rodriguez ◽  
Carlos A. Garrido ◽  
Juan C. Suarez ◽  
Francisco Huera-Huarte

In this work, CFD computations showing the dynamic response of a long flexible cylinder subject to a stepped current immersed in the wake of another cylinder are presented. These two cylinders are placed upstream in tandem configuration, where the flexible cylinder is excited by vortex shedding mechanisms. This work completes from the computational point of view, the research started 2 years ago with experiments conducted at the E.T.S.I. Navales towing tank of the Technical University of Madrid. The flexible cylinder studied is 3 m long having an external diameter of 16 mm. A combination of two codes that simulate the fluid-structure interaction phenomenon was used to obtain the velocity and pressure fields and also to measure the deformation of the cylinder at the same points where the strain gauges where placed during the experiment. This code communicates a finite volume (FV) software that solves the Navier-Stokes equations and reports the shear and pressure fields on the flexible cylinder to a second finite element (FEM) code that is able to compute stresses and deformations. Deformations are reported back to the first fluid solver in order to compute the next time step. In the experiments, only the 65% length of the cylinders were under the water surface, consequently a VOF technique was used to simulate the free surface separation between air and water. The numerical stability of these two combined codes is one of the most delicate aspects of the simulation. Taking into account that the upstream cylinder was orders of magnitude more rigid than the downstream one, we considered the upstream cylinder as stationary and consequently having no role during the FEM calculation. Boundary conditions for the flexible cylinder where such that they should imitate the universal joints used in the experiments. The fundamental natural frequencies of oscillation were monitored and compared to the towing tank experiments.


2019 ◽  
Vol 63 (4) ◽  
pp. 219-234
Author(s):  
João Baltazar ◽  
José A. C. Falcão de Campos ◽  
Johan Bosschers ◽  
Douwe Rijpkema

This article presents an overview of the recent developments at Instituto Superior Técnico and Maritime Research Institute Netherlands in applying computational methods for the hydrodynamic analysis of ducted propellers. The developments focus on the propeller performance prediction in open water conditions using boundary element methods and Reynolds-averaged Navier-Stokes solvers. The article starts with an estimation of the numerical errors involved in both methods. Then, the different viscous mechanisms involved in the ducted propeller flow are discussed and numerical procedures for the potential flow solution proposed. Finally, the numerical predictions are compared with experimental measurements.


Algorithms ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 129
Author(s):  
Yuan Li ◽  
Ni Zhang ◽  
Yuejiao Gong ◽  
Wentao Mao ◽  
Shiguang Zhang

Compared with continuous elements, discontinuous elements advance in processing the discontinuity of physical variables at corner points and discretized models with complex boundaries. However, the computational accuracy of discontinuous elements is sensitive to the positions of element nodes. To reduce the side effect of the node position on the results, this paper proposes employing partially discontinuous elements to compute the time-domain boundary integral equation of 3D elastodynamics. Using the partially discontinuous element, the nodes located at the corner points will be shrunk into the element, whereas the nodes at the non-corner points remain unchanged. As such, a discrete model that is continuous on surfaces and discontinuous between adjacent surfaces can be generated. First, we present a numerical integration scheme of the partially discontinuous element. For the singular integral, an improved element subdivision method is proposed to reduce the side effect of the time step on the integral accuracy. Then, the effectiveness of the proposed method is verified by two numerical examples. Meanwhile, we study the influence of the positions of the nodes on the stability and accuracy of the computation results by cases. Finally, the recommended value range of the inward shrink ratio of the element nodes is provided.


Processes ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 1192
Author(s):  
Dong-Hyun Kim ◽  
Jong-Chun Park ◽  
Gyu-Mok Jeon ◽  
Myung-Soo Shin

In this paper, the efficiency of Propeller Boss Cap Fins (PBCF) installed at the bulk carrier was estimated under both Propeller Open Water (POW) and self-propulsion conditions. For this estimation, virtual model-basin tests (resistance, POW, and self-propulsion tests) were conducted through Computational Fluid Dynamics (CFDs) simulation. In the resistance test, the total resistance and the wake distribution according to ship speed were investigated. In the POW test, changes of thrust, torque coefficient, and open water efficiency on the propeller according to PBCF installation were investigated. Finally, the International Towing Tank Conference (ITTC) 1978 method was used to predict the effect of PBCF installation on self-propulsive coefficient and brake horsepower. For analyzing incompressible viscous flow field, the Reynolds-Averaged Navier–Stokes (RANS) equation with SST k-ω turbulence model was calculated using Star-CCM+ 11.06.010-R8. All simulation results were validated by comparing the results of model tests conducted at the Korea Research Institute of Ships and Ocean Engineering (KRISO). Consequently, for the self-propulsion test with the PBCF, a 1.5% reduction of brake horsepower was estimated in the simulation and a 0.5% reduction of the brake horsepower was estimated in the experiment.


2017 ◽  
Vol 2017 ◽  
pp. 1-12
Author(s):  
Jui-Hsiang Kao

This research develops an Advance-Tracing Boundary Element Method in the time domain to calculate the waves that radiate from an immersed obstacle moving with random acceleration. The moving velocity of the immersed obstacle is multifrequency and is projected along the normal direction of every element on the obstacle. The projected normal velocity of every element is presented by the Fourier series and includes the advance-tracing time, which is equal to a quarter period of the moving velocity. The moving velocity is treated as a known boundary condition. The computing scheme is based on the boundary integral equation in the time domain, and the approach process is carried forward in a loop from the first time step to the last. At each time step, the radiated pressure on each element is updated until obtaining a convergent result. The Advance-Tracing Boundary Element Method is suitable for calculating the radiating problem from an arbitrary obstacle moving with random acceleration in the time domain and can be widely applied to the shape design of an immersed obstacle in order to attain security and confidentiality.


2021 ◽  
Vol 11 (4) ◽  
pp. 1399
Author(s):  
Jure Oder ◽  
Cédric Flageul ◽  
Iztok Tiselj

In this paper, we present uncertainties of statistical quantities of direct numerical simulations (DNS) with small numerical errors. The uncertainties are analysed for channel flow and a flow separation case in a confined backward facing step (BFS) geometry. The infinite channel flow case has two homogeneous directions and this is usually exploited to speed-up the convergence of the results. As we show, such a procedure reduces statistical uncertainties of the results by up to an order of magnitude. This effect is strongest in the near wall regions. In the case of flow over a confined BFS, there are no such directions and thus very long integration times are required. The individual statistical quantities converge with the square root of time integration so, in order to improve the uncertainty by a factor of two, the simulation has to be prolonged by a factor of four. We provide an estimator that can be used to evaluate a priori the DNS relative statistical uncertainties from results obtained with a Reynolds Averaged Navier Stokes simulation. In the DNS, the estimator can be used to predict the averaging time and with it the simulation time required to achieve a certain relative statistical uncertainty of results. For accurate evaluation of averages and their uncertainties, it is not required to use every time step of the DNS. We observe that statistical uncertainty of the results is uninfluenced by reducing the number of samples to the point where the period between two consecutive samples measured in Courant–Friedrichss–Levy (CFL) condition units is below one. Nevertheless, crossing this limit, the estimates of uncertainties start to exhibit significant growth.


Sign in / Sign up

Export Citation Format

Share Document