• 대한전기학회논문지:전기물성ㆍ응용부문C
• /
• 제53권8호
• /
• pp.433-442
• /
• 2004

Various models and various boundary conditions have been suggested for fluid transport simulations of high density plasma discharges such as the inductively coupled plasma discharge. In this work, we compare the various models using one-dimensional simulations based on the FDM(finite difference method), the upwind scheme, the power-law scheme, and the dielectric relaxation scheme［l］ Comparing the exactness, the numerical stability and the efficiency of the various models. the most adoptable model is suggested.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 추계 학술대회논문집
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• Wind and Structures
• /
• 제23권6호
• /
• pp.505-521
• /
• 2016

The objective of this study is to investigate numerically the effect of building roof shaps on wind flow and pollutant dispersion in a street canyon with one row of trees of pore volume, $P_{vol}=96%$. A three-dimensional computational fluid dynamics (CFD) model is used to evaluate air flow and pollutant dispersion within an urban street canyon using Reynolds-averaged Navier-Stokes (RANS) equations and the Explicit Algebraic Reynolds Stress Models (EARSM) based on k-${\varepsilon}$ turbulence model to close the equation system. The numerical model is performed with ANSYS-CFX code. Vehicle emissions were simulated as double line sources along the street. The numerical model was validated by the wind tunnel experiment results. Having established this, the wind flow and pollutant dispersion in urban street canyons (with six roof shapes buildings) are simulated. The numerical simulation results agree reasonably with the wind tunnel data. The results obtained in this work, indicate that the flow in 3D domain is more complicated; this complexity is increased with the presence of trees and variability of the roof shapes. The results also indicated that the largest pollutant concentration level for two walls (leeward and windward wall) is observed with the upwind wedge-shaped roof. But the smallest pollutant concentration level is observed with the dome roof-shaped.

• 한국환경과학회지
• /
• 제14권1호
• /
• pp.41-51
• /
• 2005

Using three-dimensional non-hydrostatical numerical model with one way double nesting technique, atmo­spheric circulation in the mountainous coastal region in summer was investigated from August 13 through 15, 1995. During the day, synoptic westerly wind blows over Mt. Mishrung in the west of a coastal city, Sokcho toward the East Sea, while simultaneously, easterly upslope wind combined with both valley wind from plain (coast) toward mountain and sea-breeze from sea toward inland coast blows toward the top of the mountain. Two different directional wind systems confront each other in the mid of eastern slope of the mountain and the upslope wind goes up to the height over 2 km, becoming an easterly return flow in the upper level over the sea and making sea-breeze front with two kinds of sea-breeze circulations of a small one in the coast and a large one in the open sea. Convective boundary layer is developed with a thickness of about 1km over the ground in the upwind side of the mountain in the west and a thickness of thermal internal boundary layer from the coast along the eastern slope of the mountain is only confined to less than 200 m. On the other hand, after sunset, no prohibition of upslope wind generated during the day and downward wind combined with mountain wind from mountain towardplain and land-breeze from land toward under nocturnal radiative cooling of the ground surfaces should intensify westerly downslope wind, resulting in the formation of wind storm. As the wind storm moving down along the eastern slop causes the development of internal gravity waves with hydraulic jump motion in the coast, bounding up toward the upper level of the coastal sea, atmospheric circulation with both onshore and offshore winds like sea-breeze circulation forms in the coastal sea within 70 km until midnight and after that, westerly wind prevails in the coast and open seas.

• 한국전산유체공학회지
• /
• 제19권3호
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.