• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.3
• /
• pp.575-583
• /
• 2017

Various types of flow patterns around the stepped weir and spillway, such as the skimming flow over such structures and the wave-type flow with a standing undular hydraulic jump and roller downstream of the structures, are developed in open channels. Unsteady three-dimensional numerical simulations are carried out using a hybrid RANS-LES turbulence modeling approach and the volume of fluid method for resolving free surface fluctuations to represent the turbulent flow including the skimming flow and wave-type flow over a stepped weir installed in a rectangular channel. The comparison of numerical results with an existing experimental measurement reveals that the present numerical simulations reasonably well reproduce the turbulent flow passing the stepped weir, in terms of time-averaged velocity profiles at selected locations downstream of the weir, flow topology characterized by the wave-type and skimming flows, the maximum height and length of the standing wave and the length of reattachment of recirculating zone. The numerical result further elucidates the distinct flow behaviors of the wave-type and skimming flow by presenting instantaneous intense variations of free surface and velocity vectors, the distributions of Reynolds shear stress and turbulent kinetic energy and three-dimensional complex features of coherent structures and total pressure distribution.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.12 no.2
• /
• pp.1005-1012
• /
• 2011

This research carried out a 3-dimensional simulation using computerized fluid dynamics (CFD) for the flow characteristics, temperature distribution, velocity distribution and residence time, etc. in a reactor in order to derive the optimal combustion conditions of an innovative combustion system. The area-weighted average temperature of the outlet of a furnace during combustion at a condition of fuel input rate 1.5 ton/hr, residence time 1.25 sec and air/fuel ratio 2.1 was $1,077^{\circ}C$, which is a suitable temperature for energy recovery and treatment of air pollutants. Exhaust gas is discharged through a duct at a 40~50 m/s maximum speed along strong vortexes at the center of a combustion chamber, so strong turbulence is created at the center of a combustion chamber to enhance the combustion speed and combustion efficiency. In this system, the optimum operation conditions to prevent incomplete combustion and suppress the formation of thermal NOx were air/fuel ratio 1.9~2.1 and fuel input rate 1.25~1.5 ton/hr.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.1
• /
• pp.92-100
• /
• 2018

Vortical systems are considered a main feature to sustain turbulence in a boundary layer through interaction. Such turbulent structures result in frictional drag and erosion or vibration in engineering applications. Research for controlling turbulent flow has been actively carried out, but in order to show the effect of vortices in a turbulent boundary layer, it is necessary to clarify the mechanism by which turbulent energy is transferred. For this purpose, it is convenient to demonstrate and capture phenomena in a laminar boundary layer. Therefore, in this study, the interactions of disturbed flow around a hemisphere on a flat plate in laminar flow were analyzed. In other words, a street of hairpin vortices was generated following a wake region formed after flow separation occurred over a hemisphere. Necklace vortices surrounding the hemisphere also appeared due to a strong adverse pressure gradient that brought high momentum fluid into the wake region thereby leading to an increase in the frequency of hairpin vortices. To mitigate the effect of these necklace vortices, local suction control was applied through a hole in front of the hemisphere. Flow visualization was recorded to qualitatively determine flow modifications, and hot-film measurements quantitatively supported conclusions on how much the power of the hairpin vortices was reduced by local wall suction.

• Wind and Structures
• /
• v.11 no.4
• /
• pp.275-289
• /
• 2008

Large eddy simulations have been performed within and over different types of urban building arrays. This paper adopted three dimensionless parameters, building frontal area density (${\lambda}_f$) the variation degree of building height (${\sigma}_h$), and the staggered degree of building range ($r_s$), to study the systematic influence of building spacing, height and layout on wind and turbulent characteristics. The following results have been achieved: (1) As ${\lambda}_f$ decrease from 0.25 to 0.18, the mean flow patterns transfer from "skimming" flow to "wake interference" flow, and as ${\lambda}_f$ decrease from 0.06 to 0.04, the mean flow patterns transfer from "wake interference" flow to "isolated roughness" flow. With increasing ${\lambda}_f$, wind velocity within arrays increases, and the vortexes in front of low buildings would break, even disappear, whereas the vortexes in front of tall buildings would strengthen and expand. Tall buildings have greater disturbance on wind than low buildings do. (2) All the wind velocity profiles and the upstream profile converge at the height of 2.5H approximately. The decay of wind velocity within the building canopy was in positive correlation with ${\lambda}_f$ and $r_s$. If the height of building arrays is variable, Macdonald's wind velocity model should be modified through introducing ${\sigma}_h$, because wind velocity decreases at the upper layers of the canopy and increases at the lower layers of the canopy. (3) The maximum of turbulence kinetic energy (TKE) always locates at 1.2 times as high as the buildings. TKE within the canopy decreases with increasing ${\lambda}_f$ and $r_s$ but the maximum of TKE are very close though ${\sigma}_h$ varies. (4) Wind velocity profile follows the logarithmic law approximately above the building canopy. The Zero-plane displacement $z_d$ heighten with increasing ${\lambda}_f$, whereas the maximum of and Roughness length $z_0$ occurs when ${\lambda}_f$ is about 0.14. $z_d$ and $z_0$ heighten linearly with ${\sigma}_h$ and $r_s$, If ${\sigma}_h$ is large enough, $z_d$ may become higher than the average height of buildings.

• Journal of Korea Water Resources Association
• /
• v.47 no.8
• /
• pp.703-715
• /
• 2014

This study aims to investigating the effect of Schmidt number (${\sigma}_c$) on sediment suspension and hydrodynamics calculation. The range of ${\sigma}_c$ is also studied based on the flux Richardson number ($Ri_f$) and gradient Richardson number ($Ri_g$). Numerical experiments are carried out by 1 dimensional vertical model. Both cohesive and non-cohesive sediments are tested under the conditions of pure current and oscillatory flow. The turbulence damping effect due to sediment suspension is examined considering ${\sigma}_c$ as a constant for the damping effect. The results of this study show the consistent effect of ${\sigma}_c$ on sediment suspension regardless of hydrodynamic condition. It is also found that the model overestimates the flow velocity and turbulent kinetic energy when the damping effect is not considered. Under the conditions of $Ri_f$ and $Ri_g$ causing density stratification, it is known that the vertical mixing of sediment is reasonably calculated in the range of ${\sigma}_c$ from 0.3 to 0.5.

• 한국신재생에너지학회:학술대회논문집
• /
• 2007.11a
• /
• pp.417-423
• /
• 2007

Oxy-gasification or oxygen-blown gasification, enables a clean and efficient use of coal and opens a promising way to CO2 capture. The coal gasification process of a slurry feed type, entrained-flow coal gasifier was numerically predicted in this paper. The purposes of this study are to develop an evaluation technique for design and performance optimization of coal gasifiers using a numerical simulation technique, and to confirm the validity of the model. By dividing the complicated coal gasification process into several simplified stages such as slurry evaporation, coal devolatilization, mixture fraction model and two-phase reactions coupled with turbulent flow and two-phase heat transfer, a comprehensive numerical model was constructed to simulate the coal gasification process. The influence of turbulence on the gas properties was taken into account by the PDF (Probability Density Function) model. A numerical simulation with the coal gasification model is performed on the Conoco-Philips type gasifier for IGCC plant. Gas temperature distribution and product gas composition are also presented. Numerical computations were performed to assess the effect of variation in oxygen to coal ratio and steam to coal ratio on reactive flow field. The concentration of major products, CO and H2 were calculated with varying oxygen to coal ratio (0.2-1.5) and steam to coal ratio(0.3-0.7). To verify the validity of predictions, predicted values of CO and H2 concentrations at the exit of the gasifier were compared with previous work of the same geometry and operating points. Predictions showed that the CO and H2 concentration increased gradually to its maximum value with increasing oxygen-coal and hydrogen-coal ratio and decreased. When the oxygen-coal ratio was between 0.8 and 1.2, and the steam-coal ratio was between 0.4 and 0.5, high values of CO and H2 were obtained. This study also deals with the comparison of CFD (Computational Flow Dynamics) and STATNJAN results which consider the objective gasifier as chemical equilibrium to know the effect of flow on objective gasifier compared to equilibrium. This study makes objective gasifier divided into a few ranges to study the evolution of the gasification locally. By this method, we can find that there are characteristics in the each scope divided.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.131-140
• /
• 2005

In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.

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

This paper proposes a novel airfoil named "KA2" for the blade of the wind turbine systems. Dynamic loads characteristics are analyzed and compared using aerodynamic data of ten airfoils including the proposed airfoil. The blade is divided into the sixteen elements in the longitudinal direction of the blade for applying the Blade Element Method Theory (BEMT) method, and in each element, torque, thrust, and pitching moment are calculated using turbulent time varying wind speed and aerodynamic data of each wing. Additionally, each force and torque is accumulated in the whole region of the blade for the estimation of representative values. The magnitude of such forces is comparatively analyzed for different airfoils. The angle of attack is constant below the rated wind speed due to the fact that the tip speed ratio is kept at the constant value, and it increases in the region of over rated wind speed as the tip speed ratio decreasing with constant rated rpm and increasing wind speed. Such increase in the angle of attack causes the changes of the force acting on the airfoil with different characteristics of lift and drag in the stall region of each different airfoil. Even though the mean wind speed is in the rated speed in a given time, because of the turbulence, it has either the over rated or under rated speed most of the time. Furthermore, the dynamic properties of each force are analyzed in this rated wind speed in order to objectively understand the dynamic properties of the blades which are designed based on the different airfoils. These dynamic properties are also compared by the standard deviation of time varying characteristics. Moreover, the output characteristics of the wind turbine are investigated with different airfoils and wind speeds. Based on these investigations, it was revealed that the proposed airfoil (KA2) is well applicable to the blade with passive pitch control system.

• Journal of Korean Society of Environmental Engineers
• /
• v.31 no.1
• /
• pp.29-34
• /
• 2009

This study was performed under four operational conditions for nitrogen removal in metal finishing wastewater. The conditions include electrode gap, reducing agent, the recycling of treated wastewater in 1st step and the simultaneous treatment of nitrate and other materials. Result showed that the removal efficiency of $NO_3{^-}-N$ was highest at the electrode gap of 10 mm. As the electrode gap was shorter than 10 mm, the removal efficiency of $NO_3{^-}-N$ decreased due to increasing in concentration polarization on electrode. And, in case that the electrode gap was longer than 10 mm, the removal efficiency of $NO_3{^-}-N$ increased with an increase in energy consumption. Because hydrogen ions are consumed when nitrate is reduced, reducing reaction of nitrate was effected more in acid solution. As 1.2 excess amount of zinc was injected, the removal efficiency of $NO_3{^-}-N$ increased due to increasing in amount of reaction with nitrate. As the effluent from 1st step in the reactor was recycled into the 1st step, the removal efficiency of $NO_3{^-}-N$ increased. Because the zinc were detached from the cathode and concentration-polarization was decreased due to formation of turbulence in the reactor. The presence of $NH_4{^+}-N$ did not affect the removal efficiency of $NO_3{^-}-N$ but the addition of heavy metal decreased the removal efficiency of $NO_3{^-}-N$. As chlorine is enough in wastewater, the simultaneous treatment of nitrate and ammonia nitrogen may be possible. The problem that heavy metal decrease the removal efficiency of $NO_3{^-}-N$ may be solved by increasing current density or using front step of electrochemical process for heavy metal removal.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.28 no.4
• /
• pp.648-656
• /
• 2022

As environmental regulations such as the International Maritime Organization (IMO)'s strategy to reduce greenhouse gases(GHG) are strengthened, technology development such as eco-friendly ships and alternative fuels is expanding. As part of this, ship propulsion technology using energy reduction and wind propulsion technology is emerging, especially in shipping companies and shipbuilders. By securing wind propulsion technology and introducing empirical research into shipbuilding and shipping, a high value-added market using eco-friendly technology can be created. Moreover, by reducing the fuel consumption rate of operating ships, GHG can be reduced by 6-8%. Rotor Sail (RS) technology is to generate a hydrodynamic lift in the vertical direction of the cylinder when the circular cylinder rotates at a constant speed and passes through the fluid. This is called the Magnus effect, and this study attempted to propose a plan to increase propulsion efficiency through a numerical analysis study on turbulence flow characteristics around RS, a wind power assistance propulsion system installed on a ship. Therefore, C_L and C_D values according to SR and AR changes were derived as parameters that affect the aerodynamic force of the RS, and the flow characteristics around the rotor sail were compared according to EP application.