• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.28 no.6
• /
• pp.332-349
• /
• 2016

OLAFOAM is the powerful CFD code and is an expanded version of $OpenFOAM^{(R)}$, for wave mechanics simulation. The $OpenFOAM^{(R)}$ does provide many solvers to correspond to each object of the numerical calculation in a variety of fields. OLAFOAM's governing equation bases on VARANS (Volume-Averaged Reynolds-Averaged Navier-Stokes) equation, and the finite volume method is applied to numerical techniques. The program is coded in C++ and run on the Linux operating system. First of all, in this study, OLAFOAM was validated for 1) wave transformation inside porous structure under bore and regular wave conditions, 2) wave transformation by submerged breakwater under regular wave condition, and 3) regular wave transformation and resultant vertical velocity distribution under current by comparison with existing laboratory measurements. Hereafter, this study, which is almost no examination carried out until now, analyzed closely variation characteristics of water surface level, wave height, frequency spectrum, breaking waves, averaged velocity and turbulent kinetic energy around porous submerged breakwater in the wave and current coexisting field for the case of permeable or impermeable rear beach. It was revealed that the wave height fluctuation according to current direction(following or opposing) was closely related to the turbulent kinetic energy, and others.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.6
• /
• pp.369-381
• /
• 2017

It has been presumed that highly nonlinear skewed waves frequently observed in a surf zone could significantly influence the transmission behaviour via a porous wave breaker due to its larger inertia force than its nonlinear counterparts of zero skewness [Cnoidal waves]. In this study, in order to confirm this perception, a numerical simulation has been implemented for 6 waves the skewness of that range from 1.02 to 1.032. A numerical simulation are based on the Tool Box called as the ihFoam that has its roots on the OpenFoam. Skewed waves are guided by the shoal of 1:30 slope, and the flow in the porous media are analyzed by adding the additional damping term into the RANS (Reynolds Averaged Navier-Stokes equation). Numerical results show that the highly nonlinear skewed waves are of higher transmitted ratio than its counterparts due to its stronger inertia force. In this study, in order to see whether or not the damping at the porous structure has an effect on the wave celerity, we also derived the dispersive relationships of Nonlinear Shallow Water Eq. [NSW] with damping at the porous structure being accounted. The newly derived dispersive relationships shows that the phase lag between the damping friction and the free surface elevation due to waves significantly influence the wave celerity.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.11-23
• /
• 1997

This paper contains a new approach to blade section design method for marine propellers. The hydrodynamic characteristics of 2-D section are highly influenced by its geometrical parameters i.e., thickness and camber distributions and leading edge radius etc. To consider fully turbulent flow field near 2-D section. the finite volume method with k-${\varepsilon}$ turbulent model which solve Reynolds time averaged Navier-Stokes(RANS) equation is applied. In this study, O-type grid system that can provide many calculation points on blade surface is used. The results were compared with those of the experiment of NACA0012 to confirm the accuracy of the developed codes. The goal of this study is the development of a blade section with high efficiency and low drag. To achieve this, we carried out the tests of lift, drag and cavitation characteristics in cavitation tunnel. The results of experiment were compared with numerical results in order to validate the proposed blades design method. By comparing the numerical results with the experiments, we found that the new blade section, KH28 allows superior performance in efficiency and cavitation avoidance characteristics. We further investigated the blade section design method and an application study of this section, KH28 to apply to the marine propeller. In order to improve the accuracy of numerical results on prediction of lift and drag, we conclude here that the 2-layer boundary model must be used.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.6
• /
• pp.470-475
• /
• 2017

Since the intake air of gas turbine engine of marine purpose contains water particles, inertial separator for separating the air and water particles are provided. Saw type and wave type separator are now used to separate inflow water particle from the gas. In this paper, the design parameters of saw type separator are studied by numerical simulations. Using the commercial CFD program, Star-CCM+, Lagrangian-Eulerian method was used to perform the analysis of two phase flow of the mist in the air. This method solves Reynolds-Averaged Navier-Stokes equations in Eulerian framework for the continuous phase, while solves equation of motion for individual particles in Lagrangian framework. Lagrangian multiphase method was applied to monitor the particles of different sizes and shapes and to verify collision between particles by chasing particles. Water particles were injected through injectors located at the inlet of the separator and escape mode was used which assumes that the particles attached on the surface of inertial separator were removed from the simulation, effectively escaping the solution domain. Through the numerical computations with the inlet condition of constant water particle size in the wetness fraction of 85%, efficiency of eliminating the water particle and the pressure drop between the inlet and outlet were examined.

• 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.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 Korean Solar Energy Society
• /
• v.35 no.2
• /
• pp.1-10
• /
• 2015

Promising onshore wind farm sites in Gangwon province of Korea were investigated in this study. Gangwon province was divided into twenty five simulation regions and a commercial program based on Reynolds averaged Navier-Stokes equation was used to find out wind resource maps of the regions. The national wind atlas with a period 2007-2009 developed by Korea institute of energy research was used as climatologies. The wind resource maps were combined to construct a wind resource map of Gangwon province with a horizontal spatial resolution of 100m. In addition to the wind resource, national environmental zoning map, distance from substation, residence and automobile road, Beakdudaegan mountain range, terrain slope, airport and military reservation district were considered to find out promising wind farm sites. A commercial wind farm design program was used to find out developable wind farm capacities in promising wind farm site with and without excluding environmental protection regions. The total wind farm capacities with and without excluding the protection regions were estimated to be 46MW and 598MW, respectively, when a 2MW commercial wind turbine was employed.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.65-65
• /
• 2016

Various types of flow conditions are developed in the region just downstream of hydraulic structures such as weir and drop structures. One of distinct flow conditions occurred downstream of drop structures is the wave type flow with undular hydraulic jump formation. We present three-dimensional numerical simulations of a wave type flow formed downstream of a stepped weir which were experimentally investigated by Kang et al. (2010). The turbulent flow over the weir structure is modeling using the unsteady Reynolds-averaged Navier-Stokes (URANS) simulation employing the Spalart-Allmaras one equation model and the detached eddy simulation. Numerical modeling and the performance of turbulence modeling approaches are evaluated by comparing with the experimental measurements in terms of the free surface variation, the shapes and sizes of undular wave, roller near at free surface, recirculation zone near the channel bottom downstream of the structures, and streamwise velocity profiles at selected longitudinal locations.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.1
• /
• pp.17-27
• /
• 2019

We carry out the numerical simulation to test a hypothesis that freak waves can be triggered by the infragravity waves of bound mode underlying the ever-present swells and its constructive interaction with swells using the Tool Box called the ihFoam that has its roots on the OpenFoam, and Bi-spectrum. Numerical simulation is implemented for the SamChcuk LNG Plant where freak waves have been reported in front of the private wharf during its construction phase due to the uncompleted northern breakwater. Infra-gravity waves of bound mode is generated using the difference wave-wave interaction between the local wind waves of 7 s and a swell of 11.4 s based on the Bi-spectrum. For the sake of comparison, numerical simulation for infra-gravity waves of free mode is also carried out. Numerical results show that stem waves along the private wharf for SamChcuk LNG Plant can be triggered by the infra-gravity waves of bound mode coming from the north, which eventually leads to freak waves when encounters the reflected waves from the south jetty.