• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.1
• /
• pp.17-25
• /
• 2020

In this study, three-dimensional fluid flow analyses have been performed in order to investigate the performance characteristics of a horizontal axis tidal turbine (HATT) by solving three-dimensional Reynolds-averaged Navier-Stokes equations utilizing the shear-stress-transport turbulence model. The computational domain for the flow analysis has been composed of hexahedral grids, and the grid dependency test has been carried out so as to determine the optimum grid size. Performance characteristics of the HATT have been investigated in consideration of the effects of hub nose geometry, inflow angle, and the tower. It has been found that the power output can be enhanced along with an increase of the ratio of the length to the diameter of the turbine nose, and the power of HATT has been reduced by approximately 10% when the primary fluid flow had an inflow angle of 15°. The power output of downstream HATT is found to be lower than that of the upstream HATT by about 1%.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.3
• /
• pp.180-198
• /
• 2019

In recent years, countries like Europe and Japan have been involved in many researches on the Low-Crested Structure (LCS) which is the method to protect beach erosion and it is regarded as an alternative to the submerged breakwaters, and compiled its results and released the design manual. In the past, studies on LCS have focused on two-dimensional wave transmission and calculating required weight of armor units, and these were mainly examined and discussed based on experiments. In this study, three-dimensional numerical analysis is performed on permeable LCS. The open-source CFD code olaFlow based on the Navier-Stokes momentum equations is applied to the numerical analysis, which is a strongly nonlinear analysis method that enables breaking and turbulence analysis. As a result, the distribution characteristics of the LCS such as water level, water flow, and turbulent kinetic energy were examined and discussed, then they were carefully compared and examined in the case of submerged breakwaters. The study results indicate that there is a difference between the flow patterns of longshore current near the shoreline, the spatial distribution of longshore and on-offshore directions of mean turbulent kinetic energy in case of submerged breakwaters and LCS. It is predicted that the difference in these results leads to the difference in sand movement.

• Journal of Energy Engineering
• /
• v.8 no.2
• /
• pp.224-232
• /
• 1999

Latent heat storage system using micro-encapsuled phase change material is effective method for floor heating of house and building. The temperature profile in capsule block and flow rate of hot water are important parameters for the development of heat storage system. In the present study, a mathematical model based on 3-D, non-steady state, Navier-Stokes equations, scalar conservation equations and turbulence model ($\kappa$-$\varepsilon$), is used to predict the temperature profiles in capsule and the velocity vectors in hot water pipe. The multi-block grids and fine grids embedding are used to join the circle in hot water pipe and square in capsule block. The phase change process of the capsule is quite complex not only because the size of phase change material is very small, but also because phase change material is mixed with the cement to form thermal storage block. In calculation, it's assumed that the phenomena of phase change is limited only the thermal properties of phase change material and the change of boundary is not happened in capsule. The purpose of this study is to calculate the temperature profiles in capsule block and velocity vectors in hot water pipe using the numerical calculation. Two kinds of thermal boundary condition were considered, the first (case 1) is the adiabatic condition for the both outside surfaces of the wall, the second (case 2) is the case in which one surface is natural convection with atmosphere and another surface is adaibatic. Calculation results are shown that the temperature profile in capsule block for case 1 is higher than that for case 2 due to less heat loss in adaibatic surface. Specially, in the domain of near Y=0, the difference of temperature is greater in case 1 than in case 2. The detailed experimental data of capsule block on the temperature profile and the thermal properties such as specific heat and coefficient of heat transfer with the various temperature are required to predict more exact phenomena of heat transfer.

• Journal of the Korean earth science society
• /
• v.38 no.6
• /
• pp.395-404
• /
• 2017

In this study, we classified observed fog events in the Middle Nakdong River near Gumi and analyzed the meteorological characteristics before and after the fog formation. The observation was performed from 2013 to 2015 using visibility meter. A total of 74 fog events were observed and most of them were classified as steam fog. The duration of observed steam fogs was longer than that of typical inland fogs because the nocturnal evaporation from the water surface was enhanced by the topographical characteristics. In order to analyze the effect of evaporation from the stream on the fog duration, the evaporation was estimated using the Penman-Monteith and the Bulk aerodynamic methods. The estimated evaporation by the Bulk method was similar to the actual evaporation from the water surface. Therefore, the Bulk method is suitable for estimating the evaporation from water surface. The evaporation amount, estimated by using the Bulk method was higher on fog days than non-fog days at 06 LST and 07 LST. The added evaporation of fog days released latent heat to the atmosphere and provided energy to maintain the turbulence in the fog. This phenomenon was confirmed by the increase of wind speed, temperature and turbulent kinetic energy within the fog.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.14 no.2
• /
• pp.114-120
• /
• 2011

The characteristics of a helical turbine to be used for tidal stream energy conversion have been numerically studied with varying a few design parameters. The helical turbines were proposed aiming at mitgating the well known poor cut-in characteristics and the structural vibration caused by the fluctuating torque, and the basic concept is introducing some twisting angle of the vertical blade along the rotation axis of the turbine. Among many potential controling parameters, we focused, in this paper, on the twisting angle and the height to diameter ratio of the turbine, and, based on the numerical experiment, We tried to propose a configuration of such turbine for which better performance can be expected. The three-dimensional unsteady RANS equations were solved by using the commercial CFD software, FLUENT with k-${\omega}$ SST turbulence model, and the grid was generated by GAMBIT. It is shown that there are a range of the twisting angle producing better efficiency with less vibration and the minimum height to diameter ratio above which the efficiency does not improve considerably.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.30 no.6
• /
• pp.270-285
• /
• 2018

The behavior of wave-induced pore water pressure inside the rubble mound and seabed, and the resultant structure failure are investigated, which are used in design of the composite breakwater representing the coastal and harbor structures. Numerical simulation techniques have been widely used to assess these behaviors through linear and nonlinear methods in many researches. While the combination of strongly nonlinear analytical method and turbulence model have not been applied yet, which can simulate these characteristics more accurately. In this study, olaFlow model considering the wave-breaking and turbulent phenomena is applied through VOF and LES methods, which gives more exact solution by using the multiphase flow analytical method. The verification of olaFlow model is demonstrated by comparing the experimental and numerical results for the interactions of regular waves-seabed and regular waves-composite breakwater-seabed. The characteristics of the spatial distributions of horizontal wave pressure, excess-pore-water pressure, mean flow velocity and mean vorticity on the upright caisson, and inside the rubble mound and seabed are discussed, as well as the relation between the mean distribution of vorticity size and mean turbulent kinetic energy. And the stability of composite breakwater are also discussed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.635-642
• /
• 2008

The influence of unsubmerged resistance bodies in a channel turbulence flow on energy loss was investigated by hydraulic experiments. Square-shaped multi-piers were used for unsubmerged structure or rigid vegetation in an open channel. In experimental channel flows multi-piers were arranged in double or single row along the channel direction, and mean-concept uniform elevations were attained and measured with a set of discharges and channel slopes. Applying the experimental results to the Manning equation, the equivalent resistance coefficient n, which implicates flow resistance and energy loss due to bottom friction as well as drag, was evaluated with varying the interval of piers and the uniform water depth. And the experimentally evaluated n values were compared with the semi-theoretical formula of the equivalent resistance coefficient derived from momentum analysis including a drag interaction coefficient. From the comparisons it was found that the interaction effect of piers on flow resistance was significant for the overall energy losses in a channel flow. The n values decrease when the interval of piers in flow-direction is less than about 2.2 times of the pier width. And it was also found that the n values increase with the 2/3 power of water depth in the theoretical formula, since the drag interaction coefficient was found to be mostly dependent on the interval of piers.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.1
• /
• pp.30-45
• /
• 2013

The present work focuses on the unsteady aerodynamics and aeroelastic properties of a small-medium sized wind-turbine blade operating under ideal conditions. A tapered/twisted blade representative of commercial blades used in an experiment setup at the National Renewable Energy Laboratory is considered. The aerodynamic loads are computed using Computational Fluid Dynamics (CFD) techniques. For this purpose, FLUENT$^{(R)}$, a commercial finite-volume code that solves the Navier-Stokes and the Reynolds-Averaged Navier-Stokes (RANS) equations, is used. Turbulence effects in the 2D simulations are modeled using the Wilcox k-w model for validation of the CFD approach. For the 3D aerodynamic simulations, in a first approximation, and considering that the intent is to present a methodology and workflow philosophy more than highly accurate turbulent simulations, the unsteady laminar Navier-Stokes equations were used to determine the unsteady loads acting on the blades. Five different blade pitch angles were considered and their aerodynamic performance compared. The structural dynamics of the flexible wind-turbine blade undergoing significant elastic displacements has been described by a nonlinear flap-lag-torsion slender-beam differential model. The aerodynamic quasi-steady forcing terms needed for the aeroelastic governing equations have been predicted through a strip-theory based on a simple 2D model, and the pertinent aerodynamic coefficients and the distribution over the blade span of the induced velocity derived using CFD. The resulting unsteady hub loads are achieved by a first space integration of the aeroelastic equations by applying the Galerkin's approach and by a time integration using a harmonic balance scheme. Comparison among two- and three- dimensional computations for the unsteady aerodynamic load, the flap, lag and torsional deflections, forces and moments are presented in the paper. Results, discussions and pertinent conclusions are outlined.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.51 no.6
• /
• pp.273-284
• /
• 2002

In this paper, computer simulations of the physical Phenomena occurring in the arc region before and after current zero were carried out to evaluate the dielectric recovery characteristics of two types of double-flow nozzles. A commercial CFD Program "PHOENICS" is used for the simulation and the user-coded subroutines to consider the arcing phenomena were added to this program by the authors. The computed results were verified by the comparison with the test results presented by the research group of BBC. In order to investigate the state of the arc region after current zero, the simulation was carried out with four steps. They are cold gas flow analysis, steady state arc simulation, transient arc simulation before current zero, transient hot gas flow simulation after current zero. The semi-experimental arc radiation model is adapted to consider the radiation energy transport and Prandtl′s mixing length model is employed as the turbulence model. The electric field and the magnetic field were calculated with the same grid structure used for the simulation of the flow field. The streamer criterion was introduced to evaluate the dielectric recovery characteristics after current zero. Compared with the results obtained by assuming the current zero state in the former studies, it has been found that the results obtained by considering the state before current zerowere more accurate.

• Fire Science and Engineering
• /
• v.13 no.4
• /
• pp.20-29
• /
• 1999

This paper describes the smoke filling process of a fire field model based on a self-developed SMEP(Smoke Movement Estimating Program) code to the simulation of fire induced flows in the two types of atrium space containing a ceiling heat flux. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k- epsilon turbulence model with buoyancy term. Compressibility is assumed and the perfect gas law is used. The results of the calculated upper-layer average temperature and smoke layer interface height has shown reasonable agreement compared with the zone models. The zone models used are the CFAST developed at the Building and Fire Research Laboratory NIST U.S.A. and the NBTC one-room of FIRECALC developed at CSIRO, Australia. The smoke layer interface heights that are important in fire safety were not as sensitive as the smoke layer temperature to the nature of ceiling heat flux condition.