• Particle and aerosol research
• /
• v.9 no.2
• /
• pp.103-110
• /
• 2013

The electrostatic precipitator (ESP) has been used for degrading atmospheric pollutants. These devices induce the electrical forces to facilitate the removal of particulate pollutants. The ions travel from the high voltage electrode to the grounded electrode by Coulomb force induced by the electric field when a high voltage is applied between two electrodes. The ions collide with gas molecules and exchange momentum with each other thus inducing fluid motion, electrohydrodynamic (EHD) flow. In this study, for the simulation of electric field and EHD flow in ESPs, an open source EHD solver, "espFoam", has been developed using open source CFD toolbox, OpenFOAM(R) (Open Field Operation and Manipulation). The electric potential distribution and ionic space charge density distribution were obtained with the developed solver, and validated with experimental results in the literature. The comparison results showed good agreement. Turbulence model is also incorporated to simulate turbulent flow; hence the developed solver can analyze laminar and turbulent flow. In distributions of electric potential and space charge, the distributions become distorted and asymmetric as the flow velocity increases. The effect of electrical drift flow was investigated for different flow velocities and the secondary flow in a flow of low velocity is successfully predicted.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.7
• /
• pp.32-40
• /
• 2003

An experimental research has been carried out for flow visualizations of the shock wave/turbulent boundary-layer interaction control utilizing the aeroelastic flaps, Spark shadowgraphs, kerosene-lampblack tracings for the surface streakline pattern, and interference fringe patterns over a thin oil-film applied at the downstream of the shock interactions have been obtained , Effects of variation of the shapes and thicknesses of the flaps are tested, and all the results are compared to the solid-wall reference case without flow-control mechanism , From the qualitative observation of the variation of skin friction utilizing the interference fringe patterns over the silicone oil-film, a strong spanwise variation of the skin friction with a narrow and long region of separation has been noticed near the centerline behind the shock structure, which phenomenon demonstrate a strong three-dimensionality of the shock interaction flows, Influence of the shape of the cavity under the flaps to the shock interaction is also tested, and it is observed that the shape of the cavity is not negligible.

• Journal of Korea Water Resources Association
• /
• v.52 no.2
• /
• pp.115-123
• /
• 2019

The effects of the Froude numbers on turbulent flow patterns downstream of a non-submerged spur dike were investigated in a laboratory flume. Three-dimensional velocities and water depths were measured using Acoustic Doppler Velocimetry and distance sensors under three Froude number conditions ($Fr_d=0.31$, 0.38, and 0.46). The results show that there are marginal differences in the velocity fields downstream of a spur dike due to the change of the Froude number. However, an increase of the Froude number was found to reduce cross-sectional area in the flow and to increase the strength of the jet-like flow. The jet-like flow was observed to displace the location of the maximum turbulence kinetic energy within a cross section toward the inner bank in the transverse direction.

• Advances in aircraft and spacecraft science
• /
• v.10 no.4
• /
• pp.303-315
• /
• 2023

The primary goal of this research is to investigate flow separation phenomena using various turbulence models. Also investigated are the effects of free-stream turbulence intensity on the flow over a NACA 0018 airfoil. The flow field around a NACA 0018 airfoil has been numerically simulated using RANS at Reynolds numbers ranging from 100,000 to 200,000 and angles of attack (AoA) ranging from 0° to 18° with various inflow conditions. A parametric study is conducted over a range of chord Reynolds numbers for free-stream turbulence intensities from 0.1 % to 0.5 % to understand the effects of each parameter on the suction side laminar separation bubble. The results showed that increasing the free-stream turbulence intensity reduces the length of the separation bubble formed over the suction side of the airfoil, as well as the flow prediction accuracy of each model. These models were used to compare the modeling accuracy and processing time improvements. The K- SST performs well in this simulation for estimating lift coefficients, with only small deviations at larger angles of attack. However, a stall was not predicted by the transition k-kl-omega. When predicting the location of flow reattachment over the airfoil, the transition k-kl-omega model also made some over-predictions. The Cp plots showed that the model generated results more in line with the experimental findings.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.6
• /
• pp.374-385
• /
• 2013

Sediment particle motions have been numerically simulated over a sinusoidal ripple. Turbulent boundary layer flows are generated by Large Eddy Simulation, and the sediment particle motions are simulated using Lagrangian particle tracking method. Two unsteady flow conditions are used in the experiment by employing two different wave amplitudes while keeping other conditions such as wave period same. As expected, the amount of suspended sediment particles is clearly dependent on the wave amplitude as it is increasing with increasing flow intensity. However, it is also observed that the pattern of suspension may be different as well due to the only different condition caused by wave amplitude. Specially, the time of maximum sediment suspension within the wave period is not coincident between the two cases because sediment suspension is strongly affected by the existence of turbulent eddies that are formed at different times over the ripple between the two cases as well. The role of these turbulent eddies on sediment suspension is important as it is also confirmed in previous researches. However, it is also found the time of these eddies' formation may also dependent on the wave amplitude over rippled beds. Therefore, it has been proved that various flow as well as geometric conditions under waves has to be considered in order to have better understanding on the sediment suspension process over ripples. In addition, it is found that high turbulent energy and strong upward flow velocities occur during the time of eddy formation, which also supports high suspension rate at these time steps. The results indicate that the relationship between the structure of flows and bedforms has to be carefully examined in studying sediment suspension at coastal regions.

• Journal of the Korean Society of Visualization
• /
• v.15 no.3
• /
• pp.20-26
• /
• 2017

The two-fluid equations are widely used to simulate two-phase flows in a nuclear reactor. For the two-fluid momentum equation, the wall and interfacial drag terms play an important role in predicting a two-phase flow behavior. Since the bubble density is much smaller than the water density, the bubble accelerates faster than the liquid in a nozzle. As a result, the bubble phase becomes faster than the liquid phase in the nozzle. In contrast, the opposite phenomena occur in the diffuser. The purpose of our study is to experimentally show these behaviors in an area-varying channel such as nozzle and diffuser. Experiments were made of turbulent bubbly flows in an area-varying horizontal channel. The velocities of the bubble and liquid phases were measured by the PIV technique. It was shown that the two-phase velocities were no longer close to each other in the area-varying regions. The bubble was faster than the liquid in the nozzle; in contrast, the bubble was slower than the liquid in the diffuser. Code simulations were also performed using the MARS code. By replacing the original wall drag model in the MARS code with Kim (1)'s wall drag partition model, we obtained the simulation results being consistent with experimental observations.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.4
• /
• pp.58-66
• /
• 1999

The flow characteristics of supersonic ejectors is often subject to compressibility, unsteadiness and shock wave systems. The numerical works carried out thus far have been of one-dimensional analyses or some Computational Fluid Dynamics(CFD) which has been applied to only a very simplified configuration. For the design of effective ejector-pump systems the effects of secondary mass flow on the supersonic ejector flow should be fully understood. In the present work the supersonic ejector-pump flows with a secondary mass flow were simulated using CFD. A fully implicit finite volume scheme was applied to axisymmetric compressible Navier-Stokes equations. The standard two-equation turbulence model was employed to predict turbulent stresses. The results obtained showed that the flow characteristics of constant area mixing tube types were nearly independent of the secondary flow rate, but the flow fields of ejector system with the second-throat were strongly dependent on the secondary flow rate due to the effect of the back pressure near the primary nozzle exit.

• Atmosphere
• /
• v.25 no.3
• /
• pp.501-510
• /
• 2015

In this study, the characteristics of flows around building groups are investigated using a computational fluid dynamics (CFD) model. For this, building groups with different volumetric ratios in a fixed area are considered. As the volumetric ratio of the building group increases, the region affected by the building group is widened. However, the wind-speed reduced area rather decreases with the volumetric ratio near the ground bottom (z ${\lesssim}$ 0.7H, here, H is the height of the building group) and, above 0.7H, it increases. As the volumetric ratio decreases (that is, space between buildings was widened), the size of recirculation region decreases but flow recovery is delayed, resulting in the wider wind-speed reduced area. The increase in the volumetric ratio results in larger drag force on the flow above the roof level, consequently reducing wind speed above the roof level. However, above z ${\gtrsim}$ 1.7H, wind speed increases with the volumetric ratio for satisfying mass conservation, resultantly increasing turbulent kinetic energy there. Inside the building groups, wind speed decreased with the volumetric ratio and averaged wind speed is parameterized in terms of the volumetric ratio and background flow speed. The parameterization method is applied to producing averaged wind speed for 80 urban areas in 7 cities in Korea, showing relatively good performance.

• Journal of computational fluids engineering
• /
• v.17 no.3
• /
• pp.68-74
• /
• 2012

Near-wall effect on wakes behind particles is one of the important factors in precise tracking of particles in turbulent flows. However, most aerodynamic force models for particles did not fully consider the wall effect. In the present study, we focused on changes of hydrodynamic forces acting on a particle depending on wall proximity. To this end, we developed an immersed boundary method with multi-direct forcing incorporated to a fully implicit decoupling procedure for incompressible flows. We validate the present immersed boundary method through two-dimensional simulations of flow over a circular cylinder. Comprehensive parametric studies on the effect of the wall proximity on the drag and lift forces acting on an immersed circular cylinder in a channel flow are performed in order to investigate general flow patterns behind the circular cylinder for a wide range of Reynolds number (0.01 ${\leq}$ Re ${\leq}$ 200). As the cylinder is closer to the wall, the drag coefficient decreases while the lift coefficient increases with a local maximum. Maximum drag and lift coefficients for different wall proximities decrease with increment of Reynolds number. Normalized drag and lift coefficients by their maximum values show universal correlations between the coefficients and wall proximity in a low Reynolds number regime (Re ${\leq}$ 1).

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.