• 한국전산유체공학회지
• /
• 제3권1호
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

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

The flow field around a high-speed train including cross-wind effects has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations in the inertial frame using the iterative time marching scheme. The governing equations are differenced with 1st-order accurate backward difference scheme for the time derivatives, 3th-order accurate QUICK scheme for the convective terms and 2nd-order accurate central difference scheme for the viscous terms. The Marker-and-Cell concept was applied to efficiently solve continuity equation, which is differenced with 2nd-order accurate central difference scheme. The 4th-order artificial damping is added to the continuity equation for numerical stability. A C-H type of elliptic grid system is generated around a high-speed train including ground. The Baldwin-Lomax turbulent model was implemented to simulate the turbulent flows. To validate the present procedure, the flow around a high speed train at constant yaw angle of $45^{\circ}\;and\;90^{\circ}$ has been simulated. The simulation shows 3-D vortex generation in the lee corner. The flow separation is also observed around the rear of the train. It has concluded that the results of present study properly agree with physical flow phenomena.

• 대한기계학회논문집B
• /
• 제20권5호
• /
• pp.1613-1623
• /
• 1996

In this paper, two discretization methods, hybrid and QUICK, are tested for the Navier-Stokes equations written in general nonorthogonal body fitted coordinates. Comparison is made by calculating two laminar flows at low Reynolds numbers of 10 - 100. One is a two-dimensional channel of gradually expanding cross section and the other is an axisymmetric flow through a circular tube having a circular constriction. Results show that the QUICK scheme results in a numerical solution more accurate than that of hybrid. The QUICK scheme also shows faster convergence for both test cases. As the number of grid points increases, all numerical solutions converge with more oscillation. The number of grid points in the y-direction(cross stream direction) is also shown to play a significant role in the approximation of convection term within separated flow zone.

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.46-53
• /
• 1997

The flow field around a three dimensional minivan-like body has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using second-order accurate schemes for the time derivatives, and third/second-order scheme for the spatial derivatives. The Marker-and-Cell concept is applied to efficiently solve continuity equation. A H-H type of multi-block grid system is generated around a three dimensional minivan-like body. Turbulent flows have been modeled by the Baldwin-Lomax turbulent model. To validate present procedure, the flows around the Ahmed body with 12.5° of slant angle are simulated. A good agreement with other numerical results is achived. After code validation, the flows around a mimivan-like body are simulated. The simulation shows three dimensional vortex-pair just behind body. The flow separation is also observed on the rear of the body. It has concluded that the results of present study properly agreed with physical flow phenomena.

• 대한기계학회논문집B
• /
• 제29권5호
• /
• pp.568-576
• /
• 2005

A 3-D finite element model combined with a volume tracking method is presented in this work to simulate the mold filling for casting processes. Especially, the analysis involves an adaptive grid method that is created under a criterion of element categorization of filling states and locations in the total region at each time step. By using an adaptive grid wherein the elements, finer than those in internal and external regions, are distributed at the surface region through refinement and coarsening procedures, a more efficient analysis of transient fluid flow with free surface is achieved. Adaptive grid based on VOF method is developed in tetrahedral element system. Through a 3-D analysis of the benchmark test of the casting process, the efficiency of the proposed adaptive grid method is verified. Developed FE code is applied to a typical industrial part of the casting process such as aluminum road wheel.

• 동력기계공학회지
• /
• 제19권4호
• /
• pp.36-42
• /
• 2015

TMixed convective flow in a bottom heated and top cooled rectangular channel can be significantly affected by the channel aspect ratio, Prandtl number, Reynolds number, Rayleigh number and angle of inclination. In such a mixed convection, the flow pattern plays an important role in various technological processes. In this study, a numerical investigation is carried out to explore mixed convection in a three-dimensional rectangular channel with bottom heated and top cooled uniformly. The three-dimensional governing equations are discretized using the finite volume method. In the range of low Reynolds number($0{\leq}Re{\leq}9.6{\times}10^{-2}$), the effects of the aspect ratio($2{\leq}AR{\leq}12$) and Gr/Re are presented and discussed. The longitudinal roll number in the channel is increased with increasing aspect ratio, and the roll number induced, regardless of the aspect ratio number, is even in the range of aspect ratios between 2 and 12, New vortex flow structure containing inclined longitudinal rolls is found, which is affected by aspect ratio and Reynolds number. The ratio Gr/Re is used to check the relative magnitudes of forced and natural convection in the mixed convective flow of high viscous fluid.

• 동력기계공학회지
• /
• 제18권3호
• /
• pp.79-86
• /
• 2014

Mixed convective vortex flow in the three-dimensional rectangular channel filled with high viscous fluid(Pr=909) is investigated computationally under various operating conditions. The Reynolds number is varied from 0 to $5{\times}10^{-1}$, the Rayleigh number from $10^3$ to $5{\times}10^4$. The three-dimensional governing equations are discretized using the finite volume method. The effects of Reynolds number and Rayleigh number are presented and discussed. From a parametric study, it is found that vortex flow pattern of mixed convection in rectangular channels can be classified into three flow patterns basically, but the new vortex flow structures containing wave rolls are found, which are affected by Rayleigh number and Reynolds number. From this results, we can draw a flow regime map to delineate various vortex flow patterns in the high viscosity fluid mixed convective flow.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2007년도 학술발표회 논문집
• /
• pp.438-443
• /
• 2007

개수로 흐름에 대한 연구는 1700년대 중엽 Chezy에 의해 이론적인 기초를 다졌으며, 광범위하고 조직적인 관측연구는 Darcy(1803-1858)에 의해 약 150년 전에 폭 2m, 길이 600m에 이르는 수로에서 실험관측을 수행하고자 시도하였다. Darcy의 동료이자 후계자인 Bazin은 제한된 조건의 현장관측뿐 아니라 다양한 조건의 수로를 제작하여 실험관측을 수행하였으며 그의 실험자료는 Bazin 자신 뿐만 아니라 Manning이나 Ganguillet와 Kutter 등 여러 연구자들의 경험식 개발에도 이용되었다. Nikuradse(1933)의 균일조도 원형관수로 실험결과로부터 관로 흐름은 층류, 천이층류, 완난류, 천이난류, 전난류 등 다섯 종류의 흐름특성을 가지고 있음을 확인하였고 Bazin(1829-1917)과 varwick(1945)의 실험결과로부터 개수로 흐름에서도 관로 흐름과 마찬가지로 층류, 천이층류, 완난류, 천이난류, 전난류 등 다섯가지의 흐름특성이 존재함을 알 수 있다. 본 연구에서는 기존의 지수형 마찰계수 산정식에서 단순히 조고만의 함수였던 ${\alpha}$에 물의 기본적인 성질인 표면장력, 점성력, 밀도와 자연 하천의 경사, 수심, 수면, 폭, 조고의 영향을 고려한 수심에 관한 무차원수 $Y_h$, 수면 폭에 관한 무차원수 $Y_b$, 조고에 관한 무차원수 $Y_k$를 도입하였다. 따라서 개수로 흐름 해석에 있어서 기존의 마찰계수 산정치보다 여러 영향을 반영하여 정확한 값을 산정할 수 있을 것으로 판단된다.

• 대한조선학회논문집
• /
• 제44권3호
• /
• pp.248-257
• /
• 2007

Turbulent flows around a sailing yacht hull were calculated using CFD technique. Grid system was generated by using Gridgen package and Fluent package was used to calculate flows around a 30ft-class yacht hull designed by MOERI. The drafts at starboard and port sides of a yacht. when the hull was heeled, were determined by adjusting the same displacement in the generated grid system. Pressure distribution on the sailing yacht hull was obtained and the changes of drag and side force by heel and leeway variation were shown. The flow calculations have been carried out both with and without free surface, and keel and rudder were included for both cases. It was found that the calculated results with free surface gives better agreement with experiments than those without free surface, although there are still a room for the improvement in correct prediction of forces.

• 한국해양공학회지
• /
• 제20권4호
• /
• pp.76-82
• /
• 2006

The large-eddy simulation(LES) technique, based an a message passing interface method(MPI), was applied to investigate the turbulent flaw phenomena around a ship. The Smagorinski model was used in the present LES simulation to simulate the turbulent flaw around a ship. The SPMD(sidsngle program multiple data) technique was used for parallelization of the program using MPI. All computations were performed an a 24-node PC cluster parallel machine, composed of 2.6 GHz CPU, which had been installed in the Advanced Ship Engineering Research Center(ASERC). Numerical simulations were performed for the Wigley hull, and the Series 60 hull(CB=0.6) using 1/4-, 1/2-, 1- and 2-million grid systems and the computational results had been compared to the experimental ones.