• Journal of the Korean Society of Combustion
• /
• v.6 no.2
• /
• pp.15-22
• /
• 2001

A pressure-based, unstructured finite volume method has been applied to couple the chemical kinetics and fluid dynamics and to capture effectively and accurately the steep gradient flame field. The pressure-velocity coupling is handled by two methodologies including the pressure-correction algorithm and the projection scheme. A stiff, operator-split projection scheme for the detailed nonequilibrium chemistry has been employed to treat the stiff reaction source terms. The conservative form of the governing equations are integrated over a cell-centered control volume with collocated storage for all transport variables. Computations using detailed chemistry and variable transport properties were performed for two laminar reacting flows: a counterflow hydrogen-air diffusion flame and a lifted methane-air triple flame. Numerical results favorably agree with measurements in terms of the detailed flame structure.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.5
• /
• pp.130-141
• /
• 1996

Numerical study of a confined, swirling, isothermal and spray-combusting flows has been presented. The pressure-velocity coupling in the Eulerian gas-phase equation is handled by the improved PISO algorithm. The droplet dispersion by turbulence is introduced by a Stochastic Separated Flow(SSF) model. The k-$\varepsilon$ turbulence model and the eddy dissipation model are employed to account for turbulence-combustion interaction. The detailed comparison with experimental data has been made for the isothermal jet swirling flows and the nearly monodisperse spray-combusting flow in the swirl combustor.

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

In this paper, some numerical methods recently developed for gas-liquid two-phase flows are reviewed. And then, a preconditioning method to solve cavitating flow by the author is introduced. This method employs a finite-difference Runge-Kutta method combined with MUSCL TVD scheme, and a homogeneous equilibrium cavitation model. So that it permits to treat simply the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Finally, numerical results such as detailed observations of the unsteady cavity flows, a sheet cavitation break-off phenomena and some data related to performance characteristics of hydrofoils are shown.

• Korean Journal of Remote Sensing
• /
• v.34 no.6_1
• /
• pp.879-891
• /
• 2018

In this study, the effects of urban redevelopment and building construction on the change of the detailed flows around the Pukyong National University (PKNU) campus located in the building-congested area was investigated using a CFD (computational fluid dynamics) model and GIS (geographic information system). For the analysis of the detailed flows before and after the constructions of the buildings around and within the campus, numerical simulations for the 16 inflow directions were performed before and after the construction. We used, as reference wind speeds at the inflow boundaries, the averaged wind speeds observed at the Gwangan light beacon (962) where there is no surrounding obstacle (i.e., building and terrain) acting as friction. We analyzed the area fractions in which wind speeds at z = 2.5 m changed after the construction for 16 inflow directions. The area fractions were relatively large in the east-south-easterly and southerly cases, because of the high-rise buildings constructed at the east and the apartment complex and the Engineering buildings constructed at the south of the PKNU campus. In the case of the easterly of which frequency is highest among the wind directions observed at the Daeyeon AWS (AWS 942) located inside the PKNU campus, the wind-speed change was not significant even after the constructions. It is shown that the building construction has affected the detailed flows around as well as even in the far downwind region of the constructed buildings. Also, it is shown that the GIS and CFD model are useful for analyzing the detailed flows in planning the urban redevelopment and/or building construction.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.81-85
• /
• 2005

Three-dimensional steady incompressible laminar entry flows in a square duct of $90^{\circ}$ bend are numerically simulated by a new solution code(PowerCFD) using unstructured cell-centered method. Solutions are obtained with three unstructured grid types of hexahedron, prism and hybrid at a Reynolds number, based on the hydraulic diameter and bulk velocity, of 790. Interesting features of the flow are presented in detail. Detailed comparisons between the computed solutions and the available experimental data are given mainly for the velocity distributions at cross-sections in a $90^{\circ}$ bend of a square duct with fully-developed entry flows. It is found that the code is capable of producing the nature of laminar flow in curved square duct with no grid type dependency.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.683-686
• /
• 2002

The hydrodynamic instability of the three-dimensional boundary-layer over a rotating disk has been numerically investigated for three cases flows using linear stability theory (i.e. Rossby number, Ro = -1, 0, and 1). Detailed numerical values of the disturbance wave number, wave frequency, azimuth angle, radius (Reynolds number, Re) and other characteristics have been calculated for $K{\acute{a}}rm{\acute{a}}n$, Ekman and $B{\"{o}}ewadt$ boundary-layer flows. Neutral curves for these flows are presented. Presented are the neutral stability results concerning the two instability modes (Type I and Type II) by using a two-point boundary value problem code COLUEW that was based upon the adaptive orthogonal collocation method using B-spline. The prediction from the present results on both instability modes among the three cases agrees with the previously known numerical and experimental data well.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.66-72
• /
• 1999

A numerical simulation method has been under development for solving turbulent flows around a ship model in maneuvering motion using the Reynolds Averaged Navier-Stokes equations. The method used second-order finite differences, collocated grids, pressure-Poisson equation and four-stage Runge-Kutta scheme as key components of the solution method. A modified Baldwin-Lomax model is used for the turbulence closure. This paper presents a preliminary result of the computational study on turbulent flows past a ship model in drift motion. Calculations are carried out for a Series 60 $C_B=0.6$ ship model, for which detailed experimental data are available. The results of the present calculations are compared with the experimental data for hydrodynamic forces acting on the model as well as velocity distributions at longitudinal sections. Only fair agreements has been achieved. The computational results show the complex asymmetrical shear flow patterns including three-dimensional separations followed by formation of bilge vortices both in bow and stern regions.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.281-283
• /
• 2001

As a companion to an adiabatic version developed by Ryu and his coworkers, we have built an isothermal magnetohydrodynamic code for astrophysical flows. It is suited for the dynamical simulations of flows where cooling timescale is much shorter than dynamical timescale, as well as for turbulence and dynamo simulations in which detailed energetics are unimportant. Since a simple isothermal equation of state substitutes the energy conservation equation, the numerical schemes for isothermal flows are simpler (no contact discontinuity) than those for adiabatic flows and the resulting code is faster. Tests for shock tubes and Alfven wave decay have shown that our isothermal code has not only a good shock capturing ability, but also numerical dissipation smaller than its adiabatic analogue. As a real astrophysical application of the code, we have simulated the nonlinear three-dimensional evolution of the Parker instability. A factor of two enhancement in vertical column density has been achieved at most, and the main structures formed are sheet-like and aligned with the mean field direction. We conclude that the Parker instability alone is not a viable formation mechanism of the giant molecular clouds.

• Journal of Power System Engineering
• /
• v.2 no.1
• /
• pp.31-38
• /
• 1998

A diffuser, an important equipment to change kinetic energy into pressure energy, has been studied for a long time. Though experimental and theoretical researches have been done, the understanding of energy transfer and detailed mechanism of energy dissipation is unclear. As far as numerical prediction of diffuser flows are concerned, various numerical studies have also been done. On the contrary, many turbulence models have constraint to the applicability of diffuser-like complex flows, because of anisotropy of turbulence near the wall and of local nonequilibrium induced by an adverse pressure gradient. The existing $k-{\epsilon}$ turbulence models have some problems in the case of being applied to complex turbulent flows. The purpose of this paper is to test the applicability of the nonlinear $k-{\epsilon}$ model concerning diffuser-like flows with expansion and streamline curvature. The results show that the nonlinear $k-{\epsilon}$ turbulence model predicted well the coefficient of pressure, velocity profiles and turbulent kinetic energy distributions, however the shear stress prediction was failed.

• Korean Journal of Remote Sensing
• /
• v.36 no.6_3
• /
• pp.1643-1652
• /
• 2020

To investigate the characteristics of detailed flows in a building-congested district, we coupled a computation fluid dynamics (CFD) model to the local data assimilation and prediction system (LDAPS), a current operational numerical weather prediction model of the Korea Meteorological Administration. For realistic numerical simulations, we used the meteorological variables such as wind speeds and directions and potential temperatures predicted by LDAPS as the initial and boundary conditions of the CFD model. We trilinearly interpolated the horizontal wind components of LDAPS to provide the initial and boudnary wind velocities to the CFD model. The trilinearly interpolated potential temperatures of LDAPS is converted to temperatures at each grid point of the CFD model. We linearly interpolated the horizontal wind components of LDAPS to provide the initial and boundary wind velocities to the CFD model. The linearly interpolated potential temperatures of LDAPS are converted to temperatures at each grid point of the CFD model. We validated the simulated wind speeds and directions against those measured at the PKNU-SONIC station. The LDAPS-CFD model reproduced similar wind directions and wind speeds measured at the PKNU-SONIC station. At 07 LST on 22 June 2020, the inflow was east-north-easterly. Flow distortion by buildings resulted in the east-south-easterly at the PKNU-SONIC station, which was the similar wind direction to the measured one. At 19 LST when the inflow was southeasterly, the LDAPS-CFD model simulated southeasterly (similar to the measured wind direction) at the PKNU-SONIC station.