• Journal of Energy Engineering
• /
• v.24 no.2
• /
• pp.110-114
• /
• 2015

A numerical investigation has been carried out on recirculation zone formed downstream of a compressible flow over a backward facing step. The study has been performed by solving Two-Dimensional Navier-Stokes equations. The system of governing equations has been solved, using an explicit Harten-Yee Non- MUSCL Modified flux type TVD scheme and a zero-equation algebraic turbulence model to calculate the eddy viscosity coefficient. The recirculation region dimensions are characterized over a range of Mach numbers of fixed step height 5mmThe detail of recirculation zone such as pressure, temperature, recirculation length, strength etc are reported. The variations of these characteristics due to change of Mach number are also presented.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.33 no.4
• /
• pp.360-369
• /
• 1997

We calculated the residual current forced by buoyancy, wind stress, and tidal stress in the Chinhae Bay using a three-dimensional diagnostic model. The calculated current was also compared with the observation. The flow directs outward from the central area of the Bay in the upper layer, and also forms eddy-shape stucture in the upper and middle layers. The flow of bottom layer shows an opposite pattern compared to those of top and middle layers. The maximum speed was 6.05 em/see (September) and 3.49 cm/sec (November) in the upper layer, and 4.39 cm/sec on both month in the middle layer. The Kinetic energy of November (8.39xlO' W) was larger than that of September (1.24xlO 'W), mainly resulting from larger buoyancy effect in September. In spite of the roughness of the grid size(1 km) and wind date, the calculated flow shows eorrelation(r=0.71) with the observation. We expect that the correlation be increased by increased by adopting the fine grid and the variable coefficients of diffusion and viscosity.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.30 no.2
• /
• pp.109-115
• /
• 1994

This paper presents a new concept to reduce turbulent frictional drag by injecting micro-bubble into near the buffer layer of turbulent boundary layer on flat plate. The concentrations of micro bubble distribution in the boundary was calculater by eddy viscosity equations in the governing equations. When near region of the buffer layer of turbulent boundary layer is filled with micro-bulle of air and viscous of the region is kept low, the velocity profile in the near region should be changed substantially. Then the Reynolds stress in the region becomes less, which guide to lower velocity gradient there. It results in reduction of velocity gradient at the viscous sublayer, which gives the reduction of shear stress at the wall.

• 한국방재학회:학술대회논문집
• /
• 2010.02a
• /
• pp.97.1-97.1
• /
• 2010

본 연구에서는 천수방정식의 매개변수인 조도계수와 난류점성계수에 따른 사행수로에서의 유속분포를 수치모의하였다. 벽면조도계수와 텐서형 난류점성계수를 변화시켜 8가지 수치모의 조건을 구성하였으며, 만곡부 정점에서의 횡방향 유속분포를 수리실험에 의한 측정값과 비교하였다. 벽면조도계수가 커질수록 벽면에서의 유속이 감소하여 실측값과 가까워지지만 벽면 이외의 영역에서는 오히려 유속이 증가하여 실측값보다 커졌다. 난류점성계수가 커질수록 만곡부 내외측의 유속차는 작게 나타났으며 xx와 yy방향으로 가중된 난류점성계수를 가지는 경우가 xy방향으로 가중된 난류점성계수를 가지는 경우에 비해 내외측 유속차가 약간 작게 나타났다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.5 no.3
• /
• pp.240-249
• /
• 1993

The suspended sediment mechanism was systematically investigated by analyzing existing data together with new data obtained under several wave folds. A simple empirical formula has been developed to predict the time-averaged suspended sediment concentration above movable bed by the effect of wave and current interaction. The present study shows linear dependency of the diffusion coefficient of suspended sediment on the eddy viscosity coefficient obtained by considering the ripple geometry. Furthermore, it is proved that reference sediment concentration is reduced by the effect of ripple asymmetry.

• The KSFM Journal of Fluid Machinery
• /
• v.9 no.4 s.37
• /
• pp.27-35
• /
• 2006

A Navier-Stokes code has been developed to simulate the flow through a cross flow fan. It is based on an unstructured finite volume method and uses moving grid technique to model the rotation of the fan. A low Reynolds number turbulence model is used to calculate eddy viscosity. The basic algorithm is SIMPLE. Numerical simulations over a wide range of flow rate are carried out to validate the code. Comparison of all numerical solutions with experimental data confirms the validity of the present code. Present numerical solutions show a noticeable improvement over a previous numerical method which is based on a model of body force to simulate the rotation of the impeller.

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.20-27
• /
• 2007

Two-dimensional stagnation flow toward a plane wall coated with magnetic fluid of uniform thickness is investigated. The flow field is represented as a similarity solution of the Navier-Stokes equation for this incompressible laminar flow. The resulting third order ordinary differential equation is solved numerically by using the shooting method and by determining two shooting parameters so as to satisfy the boundary and interface conditions. Features of the flow including streamline patterns are investigated for the varying values of density ratio, viscosity ratio, and Reynolds number. An adverse flow with double eddy pair in magnetic fluid region is found to emerge as the Reynolds number becomes higher than a threshold value. The results for the interface velocity, interface and wall shear stress, and boundary layer and displacement thickness are also presented.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.2
• /
• pp.157-166
• /
• 1999

Numerical calculations of turbulent natural convection in an enclosure of the 20 kYA oil-immersed transformer model are presented. The transformer is modelled as two concentric cylinders with different heights and diameters. The thermal boundary layers are well represented in the temperature distributions along the wall of the transformer model. The flow stratification between the hot and cold walls can not be seen in the transformer model. The turbulence eddy viscosity has its maximum at the center of the core and its maximum values at the top of the core are larger than those at the bottom of the core.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.18 no.2
• /
• pp.7-20
• /
• 1981

A method to predict the three-dimensional turbulent boundary layer on ship forms is introduced. The present differential method is in the scope of thin boundary layer theory and adopting the eddy-viscosity turbulence model. Two different numerical schemes are taken in this paper to handle the sign-changing cross-flows. The method is applied to predict the boundary layer development on real ship forms; SSPA Model 720($C_B$=0.675) and HSVA Tanker Model($C_B$=0.85). The results are qualitatively in good agreements with measurements except at the very stern. Therefore the method seems to be very promising if further developments are accomplished to handle the thick stern boundary layer effectively.

• 한국항공운항학회:학술대회논문집
• /
• 2004.11a
• /
• pp.240-246
• /
• 2004

For the prediction on the onset of oblique shock wave-induced vortex breakdown, computational studies on the Oblique Shock wave/Vortex Interaction (OSVI) are conducted and compared with both experimental results and analytic model. A Shock-stable numerical scheme, the Roe scheme with Mach number-based function (RoeM), and a two-equation eddy viscosity-transport approach are used for three-dimensional turbulent flow computations. The computational configuration is identical to available experiment, and we attempt to ascertain the effect of parameters such as a vertex strength, streamwise velocity deficit, and shock strength at a freestream Mach number of 2.49. Numerical simulations using the ${\kappa}-{\omega}SST$ turbulence model and suitably modeled vortex profiles are able to accurately reproduce many fine features through a direct comparison with experimental observations. The present computational approach to determine the criterion on the onset of oblique shock wave-induced vortex breakdown is found to be in good agreement with both the experimental result and the analytic prediction.