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

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.295-302
• /
• 2009

The flow in the draft tube cone of Francis turbines operated at partial discharge is a complex hydrodynamic phenomenon where an incoming steady axisymmetric swirling flow evolves into a three-dimensional unsteady flow field with precessing helical vortex (also called vortex rope) and associated pressure fluctuations. The paper addresses the following fundamental question: is it possible to compute the circumferentially averaged flow field induced by the precessing vortex rope by using an axisymmetric turbulent swirling flow model? In other words, instead of averaging the measured or computed 3D velocity and pressure fields we would like to solve directly the circumferentially averaged governing equations. As a result, one could use a 2D axi-symmetric model instead of the full 3D flow simulation, with huge savings in both computing time and resources. In order to answer this question we first compute the axisymmetric turbulent swirling flow using available solvers by introducing a stagnant region model (SRM), essentially enforcing a unidirectional circumferentially averaged meridian flow as suggested by the experimental data. Numerical results obtained with both models are compared against measured axial and circumferential velocity profiles, as well as for the vortex rope location. Although the circumferentially averaged flow field cannot capture the unsteadiness of the 3D flow, it can be reliably used for further stability analysis, as well as for assessing and optimizing various techniques to stabilize the swirling flow. In particular, the methodology presented and validated in this paper is particularly useful in optimizing the blade design in order to reduce the stagnant region extent, thus mitigating the vortex rope and expending the operating range for Francis turbines.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.75-80
• /
• 2009

The present paper focuses on the analysis of aero-acoustics characteristic by several boundary conditions. In this simulation, a high-order and high-resolution numerical schemes are used for the accurate computation of compressible flow with several boundary conditions including characteristic boundary conditions as well as extrapolation and zonal characteristic boundary condition. These boundary conditions are applied to the computation of two dimensional circular cylinder flows with Mach number of 0.3 and Reynolds number of 400. The computation results are validated with measurement datum and other computation results for the Strouhal frequency of vortex shedding, the mean drag coefficient and root-mean-square lift for the unsteady periodic flow regime. Secondary frequency is predicted by three kinds of boundary conditions characteristic.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.161-165
• /
• 2004

In the vortex particle method based on the vorticity-velocity formulation for solving the Wavier-Stokes equations, the unsteady, incompressible, viscous laminar flow over a NACA 0012 foil is simulated. By applying an operator-splitting method, the 'convection' and 'diffusion' equations are solved sequentially at each time step. The convection equation is solved using the vortex particle method, and the diffusion equation using the particle strength exchange(PSE) scheme which is modified to avoid a spurious vorticity flux. The scheme is improved for variety body shape using one image layer scheme. For a validation of the present method, we illustrate the early development of the viscous flow about an impulsively started NACA 0012 foil for Reynolds number 550.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1274-1284
• /
• 1999

The Regulatory Guide 1.82 recommends an analysis of hydraulic performance for sump of ECCS (Emergency Core Cooing System) when LOCA(Loss of Coolant Accident) occurs in a nuclear power plant. The present study deals with 3-dimensional, unsteady, turbulent and two-phase flow simulation to examine the behavior of mixture of reactor coolant and debris near the floor of containment building in conjunction with appropriate assumptions. The dispersed solid model has been adjusted to the interfacial momentum transfer between reactor coolant and debris. According to the results, the counterclockwiserecirculation zone had been formed in the region between sump and connection aisle about 376 second after LOCA occurs. The debris thickness accumulated on a sump screen periodically increases or decreases up to 2000 second, afterwards its peak decreases.

• Journal of computational fluids engineering
• /
• v.14 no.3
• /
• pp.49-55
• /
• 2009

This study focuses on the near-field vortical structure and dynamics of coaxial jets. The characteristics of laminar flow and mixing in coaxial jets are investigated using a unsteady flow simulation. In order to analyze the geometric effects on the vortical structure, several cases of different configurations are selected for various values of the velocity ratio of inner jet to outer jet. From the result, it is confirmed that the flow mixing is promoted by the development of vortical structure and the interaction between inner jet and outer jet. This feature is strongly related to the vortex frequency in the shear-layers. The vortex frequency depends on the velocity ratio and the lip thickness of inner nozzle, but the outer pipe length has no effect on the frequency variation.

• Journal of Advanced Marine Engineering and Technology
• /
• v.27 no.3
• /
• pp.421-428
• /
• 2003

The aerodynamic characteristics of the Wells turbine for wave energy conversion have been investigated by the numerical simulation to reproduce hysteretic behaviors. The pressure distributions on the suction surface of the blade were investigated to find out the cause of the hysteretic mechanism. The results have shown that the hysteretic behavior is associated with streamwise vertical flow appearing near the suction surface and become more obvious. as the tip clearance and solidity change. Also it has shown that such phenomena occur due to different behaviors of wakes in the accelerating and decelerating flow process. The CFD analysis shows a good agreement with experimental results.

• Journal of Navigation and Port Research
• /
• v.31 no.9
• /
• pp.793-799
• /
• 2007

The effect of seasonal wind on the tidal circulation in Jeju harbor was examined by using a numerical shallow water model. A finite element for analyzing shallow water flow is presented. The Galerkin method is employed for spatial discretization. Two step explicit finite element scheme is used to discretize the time function, which has advantage in problems treating large numbers of elements and unsteady state. The numerical simulation is compared with three cases; Case 1 does not consider the effect of wind, Case 2 and Case 3 consider the effect of summer and winter seasonal wind, respectively. According to result considering effect of seasonal wind, velocity of current vector shows slightly stronger than that of case 1 in the flow field. It can be concluded that the present method is a useful and effective tool in tidal current analysis.

• Journal of Power System Engineering
• /
• v.16 no.4
• /
• pp.51-59
• /
• 2012

The computational flow analysis using LES technique was carried out to investigate the flow characteristics of a sirocco fan under the maximum flowrate condition. The commercial SC/Tetra software was used for this unsteady and three-dimensional numerical analysis. In consequence, because a flow is unstable within the range of about 50% of a housing depth from a bellmouth around the cutoff region, the passing flow through the blade cascade occurred on the X-Y plane is a slow or a reverse with approaching to the housing inlet. Also, the secondary flow shows on the radial plane of a housing, and its vortex center exists within about 33% of a housing depth from a bellmouth except the cutoff region. Moreover, the flow occurring on the exit plane of a sirocco fan shows a complex secondary flow.

• International Journal of Fluid Machinery and Systems
• /
• v.9 no.1
• /
• pp.66-74
• /
• 2016

Concerning the numerical simulation of high-speed water jet with intensive cavitation this paper presents a practical compressible mixture flow method by coupling a simplified estimation of bubble cavitation and a compressible mixture flow computation. The mean flow of two-phase mixture is calculated by URANS for compressible fluid. The intensity of cavitation in a local field is evaluated by the volume fraction of gas phase varying with the mean flow, and the effect of cavitation on the flow turbulence is considered by applying a density correction to the evaluation of eddy viscosity. High-speed submerged water jets issuing from a sheathed sharp-edge orifice nozzle are treated when the cavitation number, ${\sigma}=0.1$, and the computation result is compared with experimental data The result reveals that cavitation occurs initially at the entrance of orifice and bubble cloud develops gradually while flowing downstream along the shear layer. Developed bubble cloud breaks up and then sheds downstream periodically near the sheath exit. The pattern of cavitation cloud shedding evaluated by simulation agrees experimental one, and the possibility to capture the unsteadily shedding of cavitation clouds is demonstrated. The decay of core velocity in cavitating jet is delayed greatly compared to that in no-activation jet, and the effect of the nozzle sheath is demonstrated.