• The KSFM Journal of Fluid Machinery
• /
• v.13 no.6
• /
• pp.42-50
• /
• 2010

The present analysis deals with a numerical procedure for optimizing the shape of elliptical dimples in a cooling channel. The three-dimensional Reynolds-averaged Navier-Stokes (RANS) analysis is employed in conjunction with the SST model for predictions of the turbulent flow and the heat transfer. Three non-dimensional geometric design variables, such as the ellipse dimple diameter ratio, ratio of the dimple depth to the average diameter, and ratio of the distance between dimples to the pitch are considered in the optimization. Twenty-one experimental points within design space are selected by Latin Hypercube Sampling. Each objective function values at these points are evaluated by RANS analysis and producing optimal point using surrogate model. The linear combination of heat transfer coefficient and friction loss related terms with a weighting factor is defined as the objective function. The results show that the optimized elliptical dimple shape improves considerably the heat transfer performance than the circular dimple shape.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.11 no.2
• /
• pp.95-106
• /
• 1999

A numerical model is represented to calculate the wave fields such as the reflected waves, the transmitted waves and the depth-averaged velocities over submerged breakwaters for the normally incident wave trains of nonlinear mono-chromatic wave and solitary wave. The finite amplitude shallow water equations with the effects of bottom friction are solved numerically in time domain using an explicit dissipative Lax-Wendroff finite difference method. The numerical model is verified by comparisons with the other numerical results and the measured data. It is found that the submerged breakwater may be more useful for protecting the energies of monochromatic waves rather than solitary waves. Finally, the armor stability on submerged breakwater is indirectly analyzed using the hydrodynamic characteristics of flow fields.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.317-324
• /
• 1994

The circulation of Lake Paldang is analysed numerically as an ultimate goal to develop a vehicle predicting the dispersion and concentration of pollutants and sediment flowed into the lake. In finite difference formulation of 2-D depth averaged governing equations. Abbott's 3-time level scheme is employed and for nonlinear terms time centering iteration technique in time and space is used. Model parameters for shear stresses and eddy diffusivities are determined through measured data in and near the lake. Predicted velocities for steady flow are shown to be close to the measured velocities and further improved by taking into account of wind effect. This indicates that the wind effect is needed for proper circulation analysis and it calls for the inclusion of the wind effect. Simulated results of unsteady flow caused by flood inflows and release through Paldang dam are found to characterize the flow features quite well as expected. This implies that the developed model can be used as a tool to analysing the circulation in the lake.

• Journal of Korea Water Resources Association
• /
• v.37 no.12
• /
• pp.979-991
• /
• 2004

In this study, RMA2 and RMA4, the 2-D depth-averaged models, were employed to simulate the two-dimensional mixing characteristics of the pollutants in the natural streams. The velocity and depth were first calculated using RMA2, 2-D hydrodynamic model, and then the resulting flow field was inputted to RMA4, 2-D water quality model, to compute the concentration field. RMA models were verified using the velocity and concentration data measured in S-curved meandering channel. The results showed that the RMA2 model simulated well the phenomenon that the maximum velocity line is located at the Inner bank of meandering channel, and the RMA4 model was well adapted to reproduce the general mixing behavior and the separation of tracer clouds. Comparing model simulations with measured data in the field experiments, RMA2 model simulated well general flow field and tendency that the maximum velocity line skewed toward the outer bank which were found in field experiments. The simulations of RMA4 model showed that the center of the tracer cloud tends to follow the path in which the maximum velocity occurs. In this study, the dispersion coefficients are fine-tuned based on the measured coefficients calculated using field concentration data, and the results show reasonable agreement with predictive equations.

• Water for future
• /
• v.21 no.1
• /
• pp.67-76
• /
• 1988

Analyzed was the circulation phenomena in the open channel with sudden expansion, by applying the Galerkin's finite element method to the depth-averaged 2-dimensional continuity and momentum equations. Wave tests were done in the simplified channel in order to review the validity of this newly developed model and the computed results were within 0.5% of $L_2$-norm error, and application of this model to the simulation of simplified dam-break gave very close results compared with the analytical solution, thus, it can be concluded that this model is valid and efficient. The main flow in the expanded channel was defined as a new initial condition with given velocity and the flow in the expanded portion was at rest in simulating the circulation, and besides the Neumann's condition the slip boundary condition for lateral wall was found to be proper condition than the no-slip condition. It can be concluded, from the numerical tests in the sudden expension, that the circulating phenomena depend mainly on the convective inertia and the effect of turbulence due to bottom shear and lateral shear is insignificant.

• Journal of computational fluids engineering
• /
• v.18 no.3
• /
• pp.72-78
• /
• 2013

The turbulent flow behavior of air supply and exhaustion in the Shin-gum-ho subway station is analyzed for ordinary and emergency state. The depth of Shin-gum-ho station is 43.6m which consists of the island-type platform(8th floor in underground) and a two-story lobby (first & second floor in underground). An emergency stairway connects between the platform and the lobby. Ventilation operation mode for ordinary state is set up as a combination of air supply and exhaustion in the lobby and platform, while for emergency state it is set up as a full air supply in the lobby and a full exhaustion in the platform. The entire station is covered for simulation. The ventilation diffusers are modeled as 95 square shapes of $0.6m{\times}0.6m$ in the lobby and as 222 square shapes of $0.6m{\times}0.6m$ and 4 rectangular shapes of $1.2m{\times}0.8m$ in the platform. The total of 7.5million grids are generated and whole domain is divided to 22 blocks for MPI efficiency of calculation. Large eddy simulation(LES) is applied to solve the momentum equation and Smagorinsky model($C_s$=0.2) is used as SGS(subgrid scale) model. The time-averaged velocity fields are compared to experimental data and show a good agreement with it.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1177-1189
• /
• 1994

This paper aims to develop the numerical model for prediction of the channel migration by analyzing of sediment and flow characteristics with patterns of channel in alluvial rivers. Flow in rivers constitutes to be the meandering or the braided form and rarely straight channel through morphologically stable patterns with mutual actions between the flowing water and bed materials. In order to develop the model for simulation of the channel migration, the channels are divided into two types with positive or negative sign by the direction of curvature radius of the centerline channel ($r_c$). That is, the single bend-channel consists of only one curvature of positive or negative sign and the multi-bend channel consists of two more curvatures of positive or negative sign, respectively. The model analyzes the sediment and flow characteristics under the influence of superelevation, spiral motion, irregularity in bed topography and depth-averaged velocity of channels. For reliability of this model, the single bend-channel and the multi bend channel are compared with experiment data in other models and the measured field data in the Keum-River, respectively. As a result, the both com parisians turn out to be excellent.

• Journal of Wetlands Research
• /
• v.13 no.1
• /
• pp.79-94
• /
• 2011

As demonstrated in study for non-submerged groynes, the flow field is predominantly two-dimensional, with mainly horizontal eddies. The eddies shed form the tips of the groynes and migrate in the flow direction. These eddies have horizontal dimensions in the order of tens of meters and time-scales in the order of minutes. In the standard flow simulations, these motions are usually not resolved, due to a too coarse grid, too large time steps and, more importantly, the use of inadequate turbulence modelling. using for example a k-${\varepsilon}$ model, it is necessary to introduce substantial modifications. Therefore simulation resolved in this study, were carried out using the DELFT-3D-MOR programme, which is part of the DELFT3D software package of WL/Delft Hydraulics and In this study, apply a two-dimensional depth-averaged model, taking an horizontal large eddy simulation(HLES). The bed morphology computed when using HLES, as well as the associated time-scale, is similar to what has been obseved in a field case. When using a mean-flow model with-out HELS, the bed morphology is less realistic and the morphological time-scale is much larger. This slow development is the result of neglecting(or averaging). the strong velocity fluctuations associated with the time-varying eddy formation.

• Journal of Korea Water Resources Association
• /
• v.48 no.11
• /
• pp.891-903
• /
• 2015

This paper presents numerical simulations of flow around a circular patch of vegetation using a depth-averaged two-dimensional numerical model which is capable of simulating flow structure in vegetated open channel. In order to account for vegetation effect, drag force terms are included in governing equations. Numerical simulations are conducted with various solid volume fractions (SVF). Flow passes through a circular patch and low velocity region, which is called wake region, is formed downstream of the patch. When SVF is larger than 0.08, a recirculation is observed. The location of the recirculation is moved further downstream as SVF decreases. Von-$K{\acute{a}}rm{\acute{a}}n$ vortex street is developed beyond the wake region due to interaction between two shear layers induced by a circular patch of vegetation. The vortex is developed as SVF is larger than 0.08, and the location of the vortex is consistent with the maximum of turbulence kinetic energy. The location of the peak of turbulence kinetic energy is moved further downstream as SVF decreases.

• Journal of Wetlands Research
• /
• v.14 no.4
• /
• pp.439-455
• /
• 2012

Coupled wave induced current modeling system(WIC) was developed from combining with the nearshore spectral wave model, SWAN, the wave induced force model, WIF, and the flow model, EFDC. The reasonable results were obtained from WIC modeling system. The ratio of the wave height calculated with respect to refraction and diffraction effects over submerged spherical shoal was occurred approximately 1~5 % errors compared to Goda(2000)'s result. The radiation stress suggested by Longuet-Higgins and Stewart(1960), the stresses due to rollers in breaking waves proposed by Dally and Osiecki(1994), and Kim(2004)'s new spreading approach instead of the previous lateral mixing approach were added to calculate wave induced force. The results of the WIC modeling system show good agreement with Nishimura et al.(1985)'s laboratory measurements and better than Kim(2004)'s 2 dimensional depth averaged numerical computations for a plane beach with detached breakwater. The present flow field computed agrees reasonably well with the measured flow field. The relative merit of WIF model in WIC modeling system is unconditional stable for time increment.