• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• Transactions of the Korean hydrogen and new energy society
• /
• v.31 no.4
• /
• pp.387-394
• /
• 2020

A parametric study of a 2.5 kW class propeller type micro hydraulic turbine was performed. In order to analyze the internal flow characteristics in the hydraulic turbine, three dimensional Reynolds-averaged Navier-Stokes equations with shear stress transport turbulence model were used and the hexahedral grid system was used to construct computational domain. To secure the reliability of the numerical analysis, the grid dependency test was performed using the grid convergence index method based on the Richardson extrapolation, and the grid dependency was removed when about 1.7 million nodes were used. For the parametric study, the axial distance at shroud span (L) between the inlet guide vane and the runner, and the inlet and outlet blade angles (β₁, β₂) of the runner were selected as the geometric parameters. The inlet and outlet angles of the runner were defined in the 3 spans from the hub to tip, and a total of 7 geometric parameters were investigated. It was confirmed that the outlet angles of the runner had the most sensitive effect on the power and efficiency of the micro hydraulic turbine.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.1
• /
• pp.72-78
• /
• 2008

In this study, advanced (fluid-structure interaction (FSI)) analysis system has been developed in order to predict turbine cascade performance with blade deformation effect due to aerodynamic loads. Intereference effects due to the relative movement of the rotor cascade with respect to the stator cascade are also considered. Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation k-ω SST turbulence models are solved to accurately predict fluid dynamic loads considering flow separation effects. A fully implicit time marching scheme based on the (coupled Newmark time-integration method) with high artificial damping is efficiently used to compute the complex fluid-structure interaction problem. Predicted aerodynamic performance considering structural deformation effect of the blade shows somewhat different results compared to the case of rigid blade model. Cascade performance evaluations for different elastic axis positions are importantly presented and its aeroelastic effects are investigated.

• Asian Journal of Atmospheric Environment
• /
• v.10 no.3
• /
• pp.125-136
• /
• 2016

In this study, the effects of the mesh barrier on the free dispersion of ammonia were numerically investigated under different atmospheric conditions. This study presents the detail and flow feature of the dispersion of ammonia through the mesh barrier on various free stream conditions to decline and limit the toxic danger of the ammonia. It is assumed that the dispersion of the ammonia occurred through the leakage in the pipeline. Parametric studies were conducted on the performance of the mesh barrier by using the Reynolds-averaged Navier-Stokes equations with realizable k-${\varepsilon}$ turbulence model. Numerical simulations of ammonia dispersion in the presence of mesh barrier revealed significant results in a fully turbulent free stream condition. The results clearly show that the flow behavior was found to be a direct result of mesh size and ammonia dispersion is highly influenced by these changes in flow patterns in downstream. In fact, the flow regime becomes laminar as flow passes through mesh barrier. According to the results, the mesh barrier decreased the maximum concentration of the ammonia gas and limited the risk zone (more than 500 ppm) lower than 2 m height. Furthermore, a significant reduction occurs in the slope of the upper boundary of $NH_3$ risk zone distribution at downstream when a mesh barrier is presented. Thus, this device highly restricts the leak distribution of ammonia in the industrial plan.

• 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.

• Wind and Structures
• /
• v.17 no.6
• /
• pp.647-670
• /
• 2013

The paper presents a numerical approach to study of fluid flow characteristics and to predict performance of wind turbines. The numerical model is based on Finite-volume method (FVM) discretization of unsteady Reynolds-averaged Navier-Stokes (URANS) equations. The movement of turbine blades is modeled using moving mesh technique. The turbulence is modeled using commonly used turbulence models: Renormalization Group (RNG) k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ and k-${\omega}$ turbulence models. The model is validated with the experimental data over a large range of tip-speed to wind ratio (TSR) and blade pitch angles. In order to demonstrate the use of numerical method as a tool for designing wind turbines, two dimensional (2-D) and three-dimensional (3-D) simulations are carried out to study the flow through a small scale Darrieus type H-rotor Vertical Axis Wind Turbine (VAWT). The flows predictions are used to determine the performance of the turbine. The turbine consists of 3-symmetrical NACA0022 blades. A number of simulations are performed for a range of approaching angles and wind speeds. This numerical study highlights the concerns with the self-starting capabilities of the present VAWT turbine. However results also indicate that self-starting capabilities of the turbine can be increased when the mounted angle of attack of the blades is increased. The 2-D simulations using the presented model can successfully be used at preliminary stage of turbine design to compare performance of the turbine for different design and operating parameters, whereas 3-D studies are preferred for the final design.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.12 no.1
• /
• pp.9-14
• /
• 2009

A numerical investigation was performed based on the Reynolds-Averaged Navier-Stokes(RANS) equations for the two-dimensional unsteady flow around a vertical axis turbine(VAT) with three or four blades. VAT is one of the promising devices for tidal current energy conversion. The geometry of the turbine blade was $NACA65_3$-018 airfoil, for which CFD analysis using Fluent was carried out at several angles of attack and the results were compared with the corresponding experimental data for validation and calibration. Then CFD simulations were carried out for the whole vertical axis turbine with a two-dimensional setup. The CFD simulation demonstrated the usefulness of the method to study the typical unsteady flows around VATs and the results showed that the optimum turbine efficiency could be achieved for carefully selected combinations of the number of blade and Tip-Speed Ratio(TSR).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.11
• /
• pp.999-1007
• /
• 2013

The present work reports numerical results of the pressure drop and heat transfer characteristics of pipes with various shapes such as circular, elliptical, circumferential wavy and twisted using a three-dimensional simulation. Numerical simulations are calculated for laminar to turbulent flows. The fully developed flow in pipes was modeled using steady incompressible Reynolds-averaged Navier-Stokes (RANS) equations. The friction and Colburn factor of each pipe are compared with those of a circular tube. The overall flow and heat transfer calculations are evaluated by the volume and area goodness factor. Finally, the objective of the investigation is to find a pipe shape that decreases the pressure loss and increases the heat transfer coefficient.

• 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.