• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.883-898
• /
• 2019

This paper employs computational tools to predict power increase (or speed loss) and propulsion performances in waves of KVLCC2. Two-phase unsteady Reynolds averaged Navier-Stokes equations have been solved using finite volume method; and a realizable k-ε model has been applied for the turbulent closure. The free-surface is obtained by solving a VOF equation. Sliding mesh method is applied to simulate the flow around an operating propeller. Towing and self-propulsion computations in calm water are carried out to obtain the towing force, propeller rotating speed, thrust and torque at the self-propulsion point. Towing computations in waves are performed to obtain the added resistance. The regular short head waves of λ/LPP = 0.6 with 4 wave steepness of H/λ = 0.007, 0.017, 0.023 and 0.033 are taken into account. Four methods to predict speed-power relationship in waves are discussed; Taylor expansion, direct powering, load variation, resistance and thrust identity methods. In the load variation method, the revised ITTC-78 method based on the 'thrust identity' is utilized to predict propulsive performances in full scale. The propulsion performances in waves including propeller rotating speed, thrust, torque, thrust deduction and wake fraction, propeller advance coefficient, hull, propeller open water, relative rotative and propulsive efficiencies, and delivered power are investigated.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.736-749
• /
• 2019

This paper employs computational tools to investigate the diffraction effects in regular head waves on the added resistance and wake on the propeller plane. The objective ships are a 66,000 DWT bulk carrier and a 3,600 TEU container ship. Fixed and free to heave and pitch conditions at design speed have been taken into account. Two-phase unsteady Reynolds averaged Navier-Stokes equations have been solved using the finite volume method; and a realizable k-ε model has been applied for the turbulent closure. The free surface is obtained by solving a VOF equation. The computations are carried out at the same scale of the model tests. Grid and numerical wave damping zones are applied to remove unwanted wave reflection at the boundaries. The computational results are analyzed using the Fourier series. The added resistances in waves at the free condition are higher than those at the fixed condition, which are nearly constant for all wavelengths. The wake velocity in waves is higher than that in calm water, and is accelerated where the wave crest locates on the propeller plane. When the vertical motion at the stern goes upward, the wake velocity also accelerated.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.2
• /
• pp.147-157
• /
• 2011

A RANS(Reynolds averaged Navier-Stokes) based numerical method is developed for the evaluation of ship resistance and self-propulsion performances. In the usability aspect of CFD for the hull form design, the field grid around practical hull forms is generated by solving a grid Poisson equation based on the hull surface grid generated from station offsets and centerline profile. A body force technique is introduced to model the effects of the propeller in which the propeller loads are obtained from potential flow analysis using an unsteady lifting surface method. The free surface is captured by using a two-phase level-set method and the realizable $k-{\varepsilon}$ model is used for turbulence closure. The hull attitude in vertical plane, i.e., trim and sinkage, is calculated by using a quasi-steady method and then considered in the computation by translating and rotating the grid system according to the values. For the validation of the proposed method, the numerical results of resistance tests for KCS, KLNG, and KVLCC1 and of self-propulsion test for KCS are compared with experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.2
• /
• pp.163-177
• /
• 2010

This article examines hydrodynamic characteristics and speed performances of a ship attached with a full spade and a twisted rudder based on a computational method. For this study, a 13,100 TEU container carrier is selected. The turbulent flows around a ship are analyzed by solving the Reynolds-averaged Navier-Stokes equation together with the application of Reynolds stress turbulence model. The computations are carried out at the conditions of rudder, bare hull, hull-rudder and hull-propeller-rudder. An asymmetric body-force propeller is applied. The speed performance is predicted by the model-ship performance analysis method of the revised ITTC'78 method. The hydrodynamic forces are compared in both rudder-open-water and self-propulsion conditions. The flow characteristics, the speed performance including propulsion factors and the rudder-cavitation performance are also compared. The model tests are conducted at a deep-water towing tank to validate the computational predictions. The computational predictions show that the twisted rudder is superior to the full spade rudder in the respect of the speed and the cavitation performances.

• Journal of the Korea Society for Simulation
• /
• v.20 no.1
• /
• pp.97-105
• /
• 2011

In rivers and streams, biofilms are thin layers of greenish-brown slime attached to rocks, plants, and other surfaces. Biofilms play key roles in primary production and cycling of nutrients, water quality remediation, suspended sediment removal, and energy flow to higher trophic levels. In the present study, we developed a two-dimensional cellular automata model to simulate mixed biofilms of toxin-sensitive and toxin-producing species in hydrodynamic flow. The flow was generated by a stochastic process for uniform flow and by using the Navier-Stokes equation for non-uniform flow. Minimized local rules governing reproduction and mortality of the species were executed in the self-organizing processes to elucidate interactions between toxin-producing and toxin-sensitive species in competition over nutrients. We briefly discuss the morphology of the simulated biofilm under different flow conditions.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.28-40
• /
• 2002

Turbulent flows around two-dimensional wing sections in ground effect are analysed by incompressible RANS equations and a finite difference method. The Baldwin-Lomax algebraic turbulence model is used to simulate high Reynolds number flows. The main purpose of this study is to clarify the two-dimensional ground effect and its flow characteristics due to different ground boundary conditions, i.e., moving and fixed bottom boundary. As a first step, to validate the present numerical code, the computational result of Clark-Y(t/C 11.7%) is compared with published numerical results and experimental data. Then, NACA4412 section in ground effect is calculated for various ground clearances with two bottom boundary conditions. According to the computational results, the difference in the lift and moment simulated with the two bottom boundary conditions is negligible, but the drag force simulated by the fixed bottom is to some extent smaller than that by the moving bottom. Therefore, it can be concluded that the drag force measured in a wind tunnel with the fixed bottom could be smaller than that with the moving bottom.

• The KSFM Journal of Fluid Machinery
• /
• v.18 no.6
• /
• pp.45-56
• /
• 2015

This paper aims to develop a submerged propeller turbine for micro hydropower plant which allows to sustain high values of efficiency in a broad range of hydrological conditions (H=2~6 m, $Q=0.15{\sim}0.39m^3/s$). The two aspects to be considered in this development are mechanical simplicity and high-efficiency operation. Unlike conventional turbines that have spiral casing and gear box, this is directing driving and no spiral casing. A 10 kW class turbine which has the most high potential of the power generation has been developed. The most important element in the design of turbine is the runner blade. The initial blade is designed using inverse design method and then the runner geometry is modified by classical hydraulic method. The design process is carried out in two steps. First, the blade shape is fix and then other components of submerged propeller turbine are designed. Computational fluid dynamics analyses based on the Navier-Stokes equations have been used to obtain overall performance data for the blade and the full turbine, respectively. The results generated by performance parameters(head, guide vane opening angle and rotational speed) variations are theoretically analysed. The evaluation criteria for the blade and the turbine performances are the pressure distribution and flow's behavior on the runner blades and turbine. The results of simulation reveals an efficiency of 91.5% and power generation of 10.5kW at the best efficiency point at the head of 4m and a discharge of $0.3m^3/s$.

• Structural Engineering and Mechanics
• /
• v.66 no.3
• /
• pp.411-421
• /
• 2018

Although a rubber isolation cushion can reduce the dynamic response of a structure itself, it has little influence on the height of a sloshing wave and even may induce magnification action. Vertical baffles are set into a base-isolated Concrete Rectangular Liquid Storage Structure (CRLSS), and baffles are opened as holes to increase the energy dissipation of the damping. Problems of liquid nonlinear motion caused by baffles are described using the Navier-Stokes equation, and the space model of CRLSS is established considering the Fluid-Solid Interaction (FSI) based on the Finite Element Method (FEM). The dynamic response of an isolated CRLSS with various baffles under an earthquake is analyzed, and the results are compared. The results show that when the baffle number is certain, the greater the number of holes in baffles, the worse the damping effects; when a single baffle with holes is set in juxtaposition and double baffles with holes are formed, although some of the dynamic response will slightly increase, the wallboard strain and the height of the sloshing wave evidently decrease. A configuration with fewer holes in the baffles and a greater number of baffles is more helpful to prevent the occurrence of two failure modes: wallboard leakage and excessive sloshing height.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.10
• /
• pp.1053-1060
• /
• 2011

This study is performed to determine the drug concentration profiles of drug-eluting stents (DES) for an ideal circular ring stent and intertwined stent models for various Reynolds numbers (Re = 200, 400, and 800). The Navier.Stokes equations coupled with the advection-diffusion equation are solved numerically in order to determine how the flow patterns and drug deposition are affected in the in-stent and post-stent regions where flow separation and recirculation occur. The presence of DES within the arterial segment affects the local drug distribution in the flow field. As a result, the drug concentration for the intertwined stent is higher over the in-stent region in comparison with the ideal stents. For a given stent geometry, the local drug concentration in the in-stent region decreases with Reynolds number, while for a given Reynolds number, the local drug concentration is relatively insensitive to the stent geometry. The results show that drug concentration along the arterial wall is significantly higher within the in-stent and post-stent regions for the intertwined stent geometry than for the ideal stent geometries.