• Journal of Korean Society for Atmospheric Environment
• /
• v.21 no.5
• /
• pp.525-535
• /
• 2005

Using a three-dimensional computational fluid dynamics (CFD) model with the $k-{\varepsilon}$ turbulence closure scheme based on the renormalization group theory, flow regimes in urban street canyons are classified according to the building and street aspect ratios. The transition between skimming flow (SF) and wake interference flow (WIF) is determined with the size of double-eddy circulation generated behind the upwind building. The transition between WIF and isolated roughness flow (IRF) is determined with the flow reattachment distance from the upwind building. The critical aspect ratios at which the flow transition occurs are found and compared with those in previous studies. The results show that the flow-regime classification method used in this study is quite reasonable and that the values of the critical aspect ratios are generally consistent with those in fluid experiments or large-eddy simulation. The regression equation describing a relation between the building and street aspect ratios at the flow-regime transition is presented.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.2
• /
• pp.10-16
• /
• 2003

The unsteady supersonic flow over tandem cavities has been analyzed by the integration of Navier-Stokes equations with the k-$\varepsilon$ turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in cavities. The upwind TVD scheme based on the flux vector split with the van Leer limiters is used. The results show the principal frequency is very reasonable. The principal frequency of the rear cavity due to the front cavity has been analyzed by the combination of the several aspect ratios of cavities. In the case of the front cavity of low aspect ratio, the frequencies of tandem cavities are almost same, because two shear layers developed from each cavity are mixed and developed to one shear layer. However, in the case of the front cavity of high aspect ratio, the characteristis of frequency are very different, because the second shear layer is developed in the diffused first shear layer.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.3
• /
• pp.129-136
• /
• 2008

All measures to cope with flooding rely on flood predictions to some extent. To investigate these predictions such as maximum water level or inundation area, a numerical model has been developed. The governing equations of the model are the two-dimensional Saint-Venant equations. The governing equations are discretized explicitly by using the leap-frog scheme and upwind scheme based on quadtree grids. The predicted numerical results have been verified by comparing to those of a Thacker problem. As a result of verifications, the present model is not only nearly four times as efficient as uniform grids but also in close agreement with the previous models. Next, the developed model is applied to several flood events in the Uiryeong basin. A general tendency is found that as a frequency is increasing, overall water levels including peak water level are increasing. At only a 500 year frequency, maximum water level is higher than 18.5 m. Therefore, it can be predictable that inundation will be generated in a 500 year frequency.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.290-297
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second order scheme.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.311-316
• /
• 2001

The thin-layer Navier-Stokes equations are solved for the hypersonic flow over blunt cone configurations with applications to laminar as well as turbulent flows. The equations are expressed in the forms of flux-vector splitting and explicit algorithm. The upwind schemes of Steger-Warming and van Leer are investigated in their ability to accurately predict the heating loads along the surface of the body. A comparison with the second order extensions of these schemes is made and a hybrid scheme incorporating a combination of central differencing and flux-vector-splitting is presented. This scheme is also investigated in its ability to accurately predict heat transfer distributions.

• Journal of computational fluids engineering
• /
• v.4 no.3
• /
• pp.35-43
• /
• 1999

In this study, the shock wave focusing of an elliptic reflector is numerically simulated by solving the Euler equations. The numerical method is the second order upwind TVD scheme with a finite volume discretization. For the verification of the present method, we simulate the moving shock wave passing through a two-dimensional corner. The computed isopycnics are compared with the earlier experiment. Numerical results of the elliptic reflectors show that the density and pressure at the focusing point increase linearly as the aspect ratio of the reflector becomes deep. On the other hand, the gas dynamic focal length decreased with the increase of the reflector aspect ratio.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1995.10a
• /
• pp.480-486
• /
• 1995

Turbulent flow field in a subchannel of bare rod bundles has been numerically simulated using the control volume based finite element method. Launder & Ying model of Reynolds stress and Lam & Bremhorst low-Reynolds number model are implemented in k-$\varepsilon$ equations and momentum equations. Secondary flows are simulated using the stream function and vorticity approach. The control volume based finite element method enable to use the upwind scheme (donor cell scheme). Sensitivity of the constants in the models are studied, and proper values are found to get the close result to the measured flow distributions.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.91-97
• /
• 2001

The thin-layer Navier-Stokes equations are solved for the hypersonic flow over blunt cone configurations with applications to laminar as well as turbulent flows. The equations are expressed in the forms of flux-vector splitting and explicit algorithm. The upwind schemes of Steger-Warming and van Leer are investigated in their ability to accurately predict the heating loads along the surface of the body. A comparison with the second order extensions of these schemes is made and a hybrid scheme incorporating a combination of central differencing and flux-vector-splitting is presented. This scheme is also investigated in its ability to accurately predict heat transfer distributions.

• Proceedings of the KSR Conference
• /
• 2004.06a
• /
• pp.685-689
• /
• 2004

Numerical simulation is conducted to clarify the heat transfer and fluid flow characteristics of HVAC(Heating, Ventilating and Air-Conditioning) for double' deck train. The HVAC system is installed under the roof of carbody. In the lay-out of HVAC system, air duct must be installed to supply air to 1st and 2nd floor respectively. The standard k-$\epsilon$ and LES models for turbulence and SIMPLE algorithm for pressure equation hased on finite volume method are used to solve the physic a] HVAC model. To assure convergence, QUICK scheme for momentum equation and the 2nd order upwind scheme for turbulent equations arc used. From the results of simulation, the temperature and velocity magnitude are also distributed uniformly in the interior of double-deck passenger car.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.39-41
• /
• 2003

n this paper, direct numerical simulation of decaying compressible turbulence with passive scalar is performed by using 7th order upwind difference scheme or 8th order group velocity control scheme. The start Reynolds number (defined by Taylor scale) is 72 and turbulent Mach numbers are 0.2-0.9. The Schmidt numbers of passive scalar are 2-10. The Batchelor k-1 range are found in scalar spectra, and the high wavenumber spectra decays faster with increasing turbulent Mach number. The extend self-similarity (ESS) is found in the passive scalar in compressible turbulence.