• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.

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

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved by successively in the moving least square sense. Some weighing functions were tested in order to investigate the up-winding effect for the convection term. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as FVM.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.41-50
• /
• 1994

An iterative time marching procedure for solving incompressible viscous flows has been applied to the flow around a propeller. This procedure solves three-dimensional Navier-Stokes equations on a moving, body-fitted, non-orthogonal grid using first-order accurate scheme for the time deivatives and second-and third-order accurate schemes for the spatial derivatives. To accelerate iterative process, a multigrid technique has been applied. This procedure is suitable for efficient execution on the current generation of vector or massively parallel computer architectures. Generally good agreement with published experimental and numerical data has been obtained. It was also found that the multigrid technique was efficient in reducing the CPU time needed for the simulation and improved the solution quality.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1998.05a
• /
• pp.469-474
• /
• 1998

이상 유동 지배 방정식에 두 상의 경계면 두께 개념을 도입하여 액체 내부에서의 기포 거동에 대한 동적 안정성을 이론적으로 해석하였다. 현재까지 기포의 안정성은 실험과 거시적인 물리적 이론에 의하여만 정적으로 연구되어 왔으나, 상간의 압력 불연속성을 표면장력 모델링을 이용하여 해석적으로 해결 함으로서 기포 안정성에 대한 정량적인 분석을 가능하게 하였다. 그 결과 기포의 안정성은 기존의 무차윈 Weber수 이외에 Reynolds수 그리고 유체의 점성력과 Capillary파와 관련된 무차원 수에 의해서도 영향을 받는 것으로 밝혀졌다.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.225-234
• /
• 1995

터빈익렬내부의 유동해석을 위해 비압축성 점성유동해석을 이용한 수치 해석 프로그램을 개발하였다. 지배방정식으로는 2차원의 비정상 비압축성 Navier-Stokes 방정식을 일반화된 곡선좌표계로 전환하여 암시적으로(implicitly) 반복적인 시간진행방법을 이용하여 유동해석을 하였다. 지배방정식의 각항들은 시간에 대해 1차의 정확도 그리고 영역에 대해서는 2차의 정확도, 대류항에 대해서는 3차의 정확도를 가지는 Upwind기법을 적용하였다. 특히, 실험적 접근이 매우 어려운 터빈의 정익과 회전하고 있는 동익과의 상호운동을 멀티블럭기법과 데이터 interface를 통해 보다 쉽게 해석할 수 있었다. 본 연구결과는 정익만을 계산한 타 연구자의 결과와의 비교시 매우 일치하였으며 물리적인 유동을 잘 파악할 수 있었다. 난류유동 해석을 위해서 Baldwin-Lomax 모델을 적용하였다.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.259-264
• /
• 2001

High-resolution simulations using vortex methods have been performed for simulating unsteady viscous flows around an impulsively started square cylinder. In order to investigate the phenomenon from laminar to transition flow, simulations are performed for Reynolds numbers 25, 50, 150 and 250. At extremely low Reynolds number, flow around a square cylinder is known to separate at the trailing edges rather than the leading edges. With an increase of Reynolds number, the flow separation at the leading edges will be developed. The main flow characteristics of developing recirculation region and separations from leading and trailing edges are studied with the unsteady behavior of the wake after the cylinder starts impulsively. A notable change in the flow evolution is found at Re=150, that is, it is shown that the flow separations begin at both leading and trailing edges of the square cylinder. On the other hand, when Re=250, the strong secondary vorticity from the rear surfaces of the square cylinder increases the drag coefficient as the primary vortex layer is pushed outwards. The comparisons between results of the present study and experimental data show a good consistency.

• Journal of Navigation and Port Research
• /
• v.27 no.5
• /
• pp.473-478
• /
• 2003

Numerical study was made on the flow characteristics around a circular pipeline between parallel walls. The incompressible Navier-Stokes equations were solved by using a third-order upwind differential scheme. When the distance near a wall is small enough, the vortex shedding is almost completely suppressed because of the interaction with the wall boundary layer separation. This study aims to clarify the characteristics of the vortex shedding regime as the body approaches a wall as Reynolds number varies. The feature of separated vorticity dynamics is analyzed at different conditions with particular attention to the interaction between the pipeline wake and the induced separation on the plane walls.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.35-43
• /
• 2008

A methodology for the simulation of compressible high Reynolds number flow over rigid and moving bodies on a structured Cartesian grid is described in this paper. The approach is based on a modified version of the Brinkman Penalization method. To avoid oscillations in the vicinity of the body and to simulate shcok-containing flows, a Weighted Essentially Non-Oscillatory scheme is used to discretize the spatial flux derivatives. For high Reynolds number viscous flow, two turbulence models of the two-equation Menter's SST URANS model and a two-equation Detached Eddy Simulation are implemented. Some simple flow examples are given to assess the accuracy of the technique. Finally, a moving grid capability is demonstrated.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.200-208
• /
• 2007

A methodology for the simulation of compressible high Reynolds number flow over rigid and moving bodies on a structured Cartesian grid is described in this paper. The approach is based on a modified version of the Brinkman Penalization method. To avoid oscillations in the vicinity of the body and to simulate shcok-containing flows, a Weighted Essentially Non-Oscillatory scheme is used to discretize the spatial flux derivatives. For high Reynolds number viscous flow, two turbulence models of the two-equation Menter's SST URANS model and a two-equation Detached Eddy Simulation are implemented. Some simple flow examples are given to assess the accuracy of the technique. Finally, a moving grid capability is demonstrated.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.1
• /
• pp.106-115
• /
• 1988

A numerical procedure for determining the airfoil shape that maximizes the lift is presented. The structure of the flow field is calculated by iteratively coupling potential flow and boundary analysis using the viscous-inviscid interaction method. The potential flow field is obtained by the vortex panel method and boundary layer flow is analyzed by means of integral approximation method which is capable of handling the laminar, transition and turbulent flow regimes. As the result of this study, it is found that the calculated flow regimes have good agreement with the existing experimented data. Davidon-Fletcher-Powell method and Augmented Lagrange Multiplier method are used for the optimal techniques. NACA 23012, NACA 65-3-21, NACA 64-2-415, NACA 64-2-A215 airfoils are used for determining the optimal airfoil shapes as a basic and compensate airfoils. Optimal design showed that the lift coefficients are increased by 17.4% at M$_{0}$=0.2 and 29% at M$_{0}$=0.3, compared with those of basic airfoil.oil.