• Journal of Korea Water Resources Association
• /
• v.38 no.11
• /
• pp.973-984
• /
• 2005

A three-dimensional thermal hydrodynamic model is developed for carrying out unsteady simulation of the selective withdrawal of the stably stratified flow in a geometrically complex, natural reservoir The governing equations are discretized on a non-staggered grid using a second-order accurate, finite-volume scheme. The numerical model is validated by applying it to simulate three-dimensional, turbulent, stratified, shear-layer flow case. The numerical predictions appear to capture reasonably well the general shape of velocity and temperature profiles observed in the laboratory experiments, while significant overestimation of the magnitude of velocity profiles is observed in the application to the flow in a natural reservoir. The physics of selective withdrawal as emerge from the numerical simulations are also discussed.

• Journal of Biomedical Engineering Research
• /
• v.13 no.4
• /
• pp.339-352
• /
• 1992

A numerical simulation of the steady and pulsatile flow across the end-to-side anastomosis was performed In order to understand the role of flow dynamics in the preferential bevel opment of distal anastomotic intimal hyperplasla. The finite element technique was employed to solve two-dimensional unsteady pulsatile flow in that region. The results of the steady flow revealed that low shear stresses occur at the proximally occluded host artery and at the recirculation region in the Inner wall just distal to the toe region of the anastomosis. The nor- mal;zed wall shear rate was increased, as was the recirculation zone size in the host artery of the by-pass graft anastomosis, with increased anastomotic junction angle. In order to min imize the size of the low wall shear region which might result in the intimal hyperplasia in the by-pass graft anastomosis, a smaller anastomotic junction angle is recommended. The pulsatile flow simulation revealed flow that regions of low and ascillating mali shear do exist near the anastomosis as In the steady simulation. The shift of stagnation point depends on the pulsation of the flow. As the flow was accelerated at systole, the stagnation point moved downstream, disappered at early diastole and reappeared during late diastole. Low shear stress was also found along both walls of the occluded proximal artery. However, the diastolic flow behavior is quite different from the steady results. The vortex near the occluded artery moved downstream and inwardly during late systole, and disappeared during diastole. Recirculations proximal to the toe and heel regions were significant during diastole. Shear stress oscillation was found along the opposite wall. The results of the present study revealed that tow shear occurs at the proximally occluded host artery aud the recirculation region in the inner wall Just dlstal to the toe region of the anastomosis. The present study suggested that the regions of fluctuated wall shear stress wit flow separation is correlated with the preferential developing regions of anastomosis neointial fibrous hyperplasia.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.5
• /
• pp.342-350
• /
• 2008

The effects of combustion instability on flow structure and flame dynamics with the inlet configurations in a model gas turbine combustor were investigated using large eddy simulation (LES). A G-equation flamelet model was employed to simulate the unsteady flame behaviors. As a result of mean flow field, the change of divergent half angle($\alpha$) at combustor inlet results in variations in the size and shape of the central toroidal recirculation (CTRZ) as well as the flame length by changing corner recirculation zone (CRZ). The case of ${\alpha}=45^{\circ}$ show smaller size and upstream location of CTRZ than those of $90^{\circ}$ and $30^{\circ}$ by the development of higher swirl velocity. The flame length in the case of ${\alpha}=45^{\circ}$ is shorter than other cases, while the case of ${\alpha}=30^{\circ}$ yields the longest flame length due to the decrease of effective reactive area with the absence of CRZ. Through the analysis of pressure fluctuation, it was identified that the case of ${\alpha}=45^{\circ}$ shows the largest damping effect of pressure oscillation in all configurations and brings in the noise reduction of 2.97dB, compared to that of ${\alpha}=30^{\circ}$ having the largest pressure oscillation. These reasons were discussed in detail through the analysis of unsteady phenomena related to recirculation zone and flame surface. Finally the effects of flame-acoustic interaction were evaluated using local Rayleigh parameter.

• Ecology and Resilient Infrastructure
• /
• v.10 no.4
• /
• pp.217-225
• /
• 2023

This study analyzed the effect of flood level reduction in the case of creating and operating offline storage for the Jangdong district, which can be used as a flood buffer space for the Geumgang River, through one-dimensional unsteady flow numerical simulation. In particular, the sensitivity analysis of changes in the height and width (length) of transverse weirs on flood level changes was performed to provide quantitative information necessary for flood control facility (embankment) design. As a result of analyzing the flood control effect of the offline storage based on the peak flood discharge and level, spatially, the flood control effect at the planned flood buffer space site and the downstream end was confirmed, and it was confirmed that the flood reduction effect at the downstream occurred the most. By design conditions of the transverse overflow weir, the greatest flood reduction effect was found under the condition that the overflow weir height based on the 50-year frequency flood level and the transverse overflow weir width (length) of 125 m were considered. The effect of delaying the time to reach the maximum flood due to the operation of the offline storage site was also presented based on unsteady flow modeling.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.6
• /
• pp.650-655
• /
• 2016

In this study, a computational simulation of the internal flow characteristics was carried out for a Urea-SCR Injector. A single hole swirl injector with a swirl disk and slanted nozzle was used in this simulation. The maximum needle lift and opening velocity were selected as the design parameters. To analyze the unsteady internal flow characteristics of the Urea-SCR injector, the moving grid technique was applied to simulate the delicate needle movement. According to the simulation results, the injected mass flow rate from the Urea-SCR injector decreased with increasing needle opening velocity and maximum needle lift. This is because the Urea-solution tends to fill the empty space that the needle previously occupied. The swirl flow is decreased as the flow goes through the injector nozzle, because of the friction with the nozzle wall. Also, during the maximum needle lift period, the swirl coefficient and mean swirl coefficient increase with increasing needle lift. The results of this study may be used as the basic design data of related injectors.

• Journal of Ship and Ocean Technology
• /
• v.5 no.2
• /
• pp.50-61
• /
• 2001

In this paper, nonlinear interactions between water waves and a horizontally submerged circular cylinder are numerically simulated. In this case, the nonlinear interactions between them generated a wave breaking phenomenon. The wave breaking phenomenon plays an important role in the wave farce. Negative drifting forces are raised at shallow submerged cylinders under waves because of the wave breaking phenomenon. For the numerical simulation, a finite difference method based on the unsteady incompressible Navier-Stokes equations and the continuity equation is adopted in the rectangular grid system. The free surface is simulated with a computational simulation method of two-layer flow by using marker density. The results are compared with some existing computational and experimental results.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.11a
• /
• pp.344-352
• /
• 2006

Evacuation performance, starting transient, and plume blowback at diffuser breakdown of a straight cylindrical supersonic exhaust diffuser with no externally supplied secondary flow are investigated. Pressure records in the transitional periods are measured by a small-scale cold-gas simulator. Flow-fields evolving in the diffuser-type ejector are solved by preconditioned Favre-averaged Navier-Stokes equations with a low-Reynolds number $k-{\varepsilon}$ turbulence model edited for turbulence compressibility effects. The present RANS method is properly validated with measured static wall pressure distributions and evacuation level at steady operation as well as the pressure records during the transition regime.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.557-560
• /
• 2002

In this study, simulation of weak shock waves are peformed by a two-dimensional thermal fluid or compressible fluid model of the lattice Boltzmann method. The shock wave represents an abrupt change in fluids properties, in which finite variations in pressure, internal energies, and density occur over the shock thickness. The characteristics of the proposed model with a simple distribution function is verified by calculation of the sound speeds, and the shock tube problem. The reflection of a weak shock wave by wedge propagating in a channel is performed. The results agree well with those by finite difference method or by experiment. In the simulation of unsteady shock wave diffraction around a sharp corner, we show a flow field of vortical structure near the comer.

• Journal of computational fluids engineering
• /
• v.17 no.4
• /
• pp.24-31
• /
• 2012

Effects of the subgrid length-scale in the Delayed-Detached Eddy Simulation(DDES) are investigated based on the Spalart-Allmaras(S-A) and the k-$\omega$ Shear Stress Transport(SST) turbulence models. Driver & Seegmiller's experimental results are used to validate numerical results. Grid convergence with grid resolution and subgrid length-scale is investigated. The simulation results show that the volume method for the subgrid length-scale is more resistant to unfavorable effects of the grid size in the periodic direction than the maximum method. Using a sufficient grid resolution and an appropriate subgrid length-scale, both S-A based DDES and SST based DDES methods can provide a good correlation with the experimental data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.1 s.256
• /
• pp.68-75
• /
• 2007

A three-dimensional computation is conducted to simulate a three-dimensional rotating stall in a low speed axial compressor. It is generally known that a tip leakage flow has an important role on a stall inception. However, almost of researchers have taken no interest in a role of the hub-comer-stall on the rotating stall even though it is a common feature of the flow in an axial compressor operating near stall and it has a large effect on the flows and loss characteristics. Using a time-accurate unsteady simulation, it is found that the hub-comer-stall may be a trigger to collapse the axisymmetric flows under high loads. An asymmetric disturbance is initially originated in the hub-comer-stall because separations are naturally unstable flow phenomena. Then this disturbance is transferred to the tip leakage flows from the hub-comer-stall and grows to be stationary stall cells, which adheres to blade passage and rotate at the same speed as the rotor. When stationary stall cells reach a critical size, these cells then move along the blade row and become a short-length-scale rotating stall. The rotational speed of stall cells quickly comes down to 79 percent of rotor so they rotate in the opposite direction to the rotor blades in the rotating frame.