• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.6
• /
• pp.373-381
• /
• 2016

In this study, we investigate the characteristics of the drag reduction of a circular cylinder having a detached splitter plate at the wake side. We measure the fluid force on a circular cylinder and visualize the field using particle image velocimetry (PIV) with a high Reynolds number, Re = 10,000. The experimental paraeters used were the width ratios (H/B = 0.5~1.5) of splitters to the prism width and the gap ratios (G/B = 0~2) between the prism and the splitter plate. The drag-reduction rate of the circular cylinder increased with H/B in the case of the same G/B, and it increased and then decreased with G/B in the case of the same H/B. The vortices of the opposite direction on the upper and lower sides of the detached splitter plate were generated by installing the plate. Reverse flow was caused by the vortices at the wake region of the circular cylinder, and the drag of the circular cylinder was decreased by the reverse flow.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

• Journal of computational fluids engineering
• /
• v.17 no.4
• /
• pp.81-86
• /
• 2012

The objective of the present study is to investigate the effects of various turbulence models on the aerodynamic noise of an air-conditioner (AC) indoor unit. The results from URANS (unsteady Reynolds-averaged Navier-Stokes) simulations with the standard k-$\varepsilon$, k-$\omega$ shear stress transport (SST) and Spalart-Allmaras (S-A) turbulence models were analyzed and compared with the noise data from the experiments. The frequency spectra of the far-field acoustic pressure were computed using the Farrasat equation derived from the Ffowcs Williams-Hawkings (FW-H) equation based on the acoustic analogy model. Two fixed fan casings and the rotating cross-flow fan were used as the source surfaces of the dipole noise in the Farrasat equation. The result with the standard k-$\epsilon$ model showed a much better agreement with the experimental data compared to the k-w SST and S-A models. The differences in the pressure spectra from the different turbulence models were discussed based on the instantaneous vorticity fields. It was found that the over-estimated power spectra with the k-w SST and S-A models are related to the emphasized small-scale vortices produced with these models.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.6
• /
• pp.16-22
• /
• 2004

Unsteady vortex shedding behind a rectangular cylinder near a wall influences both increasing of drag and dynamic stability of heavy vehicle, bridge or building. Incompressible Averaged Navier-Stokes equation with modified ${\varepsilon}-SST$ turbulence model is adapted for investigating the flow field between the rectangular cylinder and the wall. In case the vortex shedding happens, not only the averaged maximum velocity is higher than other cases, but the position of the maximum velocity is closer to the lower surface of rectangular cylinder. On this study, it is confirmed that the vortex shedding behind a rectangular cylinder can be suppressed by the passive control method added by horizontal and vertical fences to the lower surface of rectangular cylinder.

• Journal of Korean Society on Water Environment
• /
• v.25 no.1
• /
• pp.48-57
• /
• 2009

A two-dimensional (2D), laterally-averaged hydrodynamic and water quality model, CE-QUAL-W2 was applied to evaluate the performance on simulating the effect of flushing from Daecheong Reservoir on the downstream water quality variations during the flushing events held on November, 2003 and March, 2008. The hydraulic and water quality simulation results were compared with field measurement data, as well as a one-dimensional (1D), unsteady model (KORIV1) that revealed limited capability in the previous study due to missing the resuspension process of river bottom sediments. The results showed that although the 2D model made satisfactory performance in reproducing the temporal variations of dissolved matters including phosphate, ammonia and nitrate, it revealed poor performance in simulating the increase of biological oxygen demand and suspended sediment (SS) concentrations during the passage of the flushing flow. The reason of the error was that the resuspension process of the 2D model is only the function of shear stress induced by wind. In reality, however, as shown by significant correlation between bottom shear stress ($\tau$) and observed SS concentration, the resuspension process can be significantly influenced by current velocity in the riverine system, especially during flushing event. The results indicate that the resuspension of river bottom materials should be incorporated into the water quality modeling processes if $\tau$ is greater than a critical shear stress (${\tau}_c$) for better simulation of flushing effect.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.2
• /
• pp.147-157
• /
• 2011

A RANS(Reynolds averaged Navier-Stokes) based numerical method is developed for the evaluation of ship resistance and self-propulsion performances. In the usability aspect of CFD for the hull form design, the field grid around practical hull forms is generated by solving a grid Poisson equation based on the hull surface grid generated from station offsets and centerline profile. A body force technique is introduced to model the effects of the propeller in which the propeller loads are obtained from potential flow analysis using an unsteady lifting surface method. The free surface is captured by using a two-phase level-set method and the realizable $k-{\varepsilon}$ model is used for turbulence closure. The hull attitude in vertical plane, i.e., trim and sinkage, is calculated by using a quasi-steady method and then considered in the computation by translating and rotating the grid system according to the values. For the validation of the proposed method, the numerical results of resistance tests for KCS, KLNG, and KVLCC1 and of self-propulsion test for KCS are compared with experimental data.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.3
• /
• pp.245-259
• /
• 1994

This paper is concerned with the development of a three-dimensional numerical model for coastal circulation and heat transport with improved prediction ability. The model uses fully nonlinear, time-dependent three-dimensional, $\sigma$-transformed equations of motion and equation of heat transport The model was verified with experimental data for wind-driven current in a one-dimensional channel and thermal jets flowing into stagnant waters and applied for unsteady flow induced by tide and thermal jets in coastal waters around Kori nuclear power plant. The model results were in good agreements with experimental data sets for wind-driven current and thermal jet, and field observed data sets in coastal waters. This study has shown that the $\kappa$-$\varepsilon$ turbulence model is applicable to various coastal conditions without any modification of turbulence constants.

• Journal of Korea Water Resources Association
• /
• v.54 no.8
• /
• pp.629-641
• /
• 2021

In this study, a two-dimensional numerical model (Nays2DH) was applied to analyze the process of morphological changes in the river channel bed depending on the changes in the amount of flooding after fully opening the Sejong weir, which was constructed upstream of the Geum River. For this, numerical simulations were performed by assuming the flow conditions, such as a non-uniform flow (NF), unsteady flows (single flood event, SF), and a continuous flood event (CF). Here, in the cases of the SF and CF, the normalized hydrograph was calculated from real flood events, and then the hydrograph was reconfigured by the peak flow discharge according to the scenario, and then it was employed as the flow discharge at the upstream boundary condition. In this study, to quantitatively evaluate the morphological changes, we analyzed the time changes in the bed deformation the bed relief index (BRI), and we compared the aerial photographs of the study area and the numerical simulation results. As simulation results of the NF, when the steady flow discharge increases, the ratio of lower width to depth decreases and the speed of bar migration increases. The BRI initially increases, but the amount of change decreased with time. In addition, when the steady flow discharge increases, the BRI increased. In the case of SF, the speed of bar migration decreased with the change of the flow discharge. In terms of the morphological response to the peak flood discharge, the time lag also indicated. In other words, in the SF, the change of channel bed indicates a phase lag with respect to the hydraulic condition. In the result of numerical simulation of CF, the speed of bar migration depending on the peak flood discharges decreased exponentially despite the repeated flood occurrences. In addition, as in the result of SF, the phase lag indicated, and the speed of bar migration decreased exponentially. The BRI increased with time changes, but the rate of increase in the BRI was modest despite the continuous peak flooding. Through this study, the morphological changes based on the hydrological characteristics of the river were analyzed numerically, and the methodology suggested that a quantitative prediction for the river bed change according to the flow characteristic can be applied to the field.

• Journal of the Korean GEO-environmental Society
• /
• v.16 no.1
• /
• pp.45-53
• /
• 2015

This study has obtained Soil Water Retention Curve (SWRC) of the unsaturated soil from the volumetric pressure plate extractor test and the triaxial compression tests was also conducted. By using the rainfall data measured in the site the seepage analysis of unsteady flow was performed with the program of SEEP/W in Geostudio 2007 and stability of the slope was analyzed with SLOPE/W program. Results of analyses show that shear strength of the unsaturated soil increases with the increase of matric suction. And it was also found that the net volumetric stress and the apparent cohesion increased with the matric suction. The seepage analysis of rainfall represents that the increasing rate of negative pore pressure at the zone of large negative pore pressure is appeared to be high even though lower rainfall intensity, but this tendency declines with ground depth. The stability analysis of slope was carried out for the actual plane of failure with the data representing the field condition. The factor of safety thus calculated was about unity (1.0) or just below, which means that the adopted method of analysis is in good agreement with the field condition.

• Journal of Korea Water Resources Association
• /
• v.44 no.5
• /
• pp.363-376
• /
• 2011

1-D hydrodynamic numerical models have been most widely used in the field of flood analysis. The model's input data are upstream/downstream boundaries, roughness coefficients, cross-sections, and so on, and computational distance step and time step are the most important factors in order to guarantee the computational accuracy, stability, and efficiency. In this study, a theoretical explanation is presented for the basis of the previous empirical selection criteria of cross-section's location; also, the estimation technique of computationally variable distance step is proposed to reflect the properties of flow at every computational time step. Combining this technique with 1-D unsteady numerical model, it was applied to two events of Teton dam failure flood and the Han River flood. The numerical experimental results demonstrate that the accuracy and stability is increased when used more interpolated cross-sections and show that the proposed technique of computationally variable distance step has the same order of accuracy with smaller numbers of cross-section than previous empirical selection criteria. The practical use of this technique will be possible to analyze the river floods with high efficiency as well as accuracy and stability.