• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 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.

• 한국건축시공학회지
• /
• 제19권1호
• /
• pp.39-46
• /
• 2019

본 연구에서는 비압축성 Navier-Stokes 방정식을 적용한 ISPH 기법을 이용하여 3차원 유동 수치해석 모델을 개발하였다. 수치해석을 위해 MATLAB을 사용하여 ISPH 프로그램을 구현하였다. ISPH의 커널 함수로 piecewise cubic spline 함수를 사용하였다. 벽 경계조건으로 고정 가상 입자를 사용하였으며, 가상 밀도를 적용하여 자유 표면 경계 부근의 입자들을 결정하였다. 수치해석 모델과 코드의 정도를 확인하기 위해 $T_{500}$ 시험, 슬럼프 플로우 시험, L-box 시험의 수치해석 결과와 실험 결과를 비교하였다. 수치해석 결과 고유동 콘크리트의 점성계수 및 항복응력 변화에 따른 유동 현상의 특성을 잘 묘사하였으며, 기존의 실험값과 비교적 잘 일치함을 확인할 수 있었다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.299-300
• /
• 2003

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제16권1호
• /
• pp.65-84
• /
• 2012

We present an application of the variational multiscale methodology to the computation of concentric annular pipe flow. Isogeometric analysis is utilized for higher order approximation of the solution using Non-Uniform Rational B-Splines (NURBS) functions. The ability of NURBS to exactly represent curved geometries makes NURBS-based isogeometric analysis attractive for the application to the flow through the curved channel.

• 대한기계학회논문집B
• /
• 제21권3호
• /
• pp.341-349
• /
• 1997

The characteristics of defect correction method are discussed in a sample heat conduction problem showing the numerical solution of the error correction equation can predict the error of the numerical solution of the original governing equation. A way of using defect correction method combined with the existing algorithm for the incompressible fluid flow, is proposed and subsequently tested for the driven square cavity problem. The error correction equations for the continuity equation and the momentum equations are considered to estimate the errors of the numerical solutions of the original governing equations. With this new approach, better velocity and pressure fields can be obtained by correcting the original numerical solutions using the estimated errors. These calculated errors also can be used to estimate the orders of magnitude of the errors of the original numerical solutions.

• Journal of applied mathematics & informatics
• /
• 제22권1_2호
• /
• pp.67-82
• /
• 2006

In this paper, we present a finite difference method for the implementation of the rotation of a circular cylinder in the incompressible flow field by solving the two-dimensional unsteady Navier-Stokes equations. The approach is to use staggered grid method so that the accuracy and order of convergence of the associated algorithms can be maintained. The proposed method is easy to be implemented and is effective. A set of simulations for the flow dynamics is provided to show the computational results.

• 한국CDE학회논문집
• /
• 제1권1호
• /
• pp.56-64
• /
• 1996

A quadrilateral-triangular mixed grid method for the solution of incompressible viscous flow is presented. The solution domain near the body surface is meshed using elliptic grid geneator to acculately simulate the viscous flow. On the other hand, we used unstructured triangular grid system generated by advancing front technique of a simple automatic grid generation algorithm in the rest of the computational domain. The present method thus is capable of not only handling complex geometries but providing accurate solutions near body surface. The numerical technique adopted here is PISO type finite element method which was developed by the present author. Investigations have been made of two-dimensional unsteady flow of Re=550 past a circular cylinder. In the case of use of the unstructured grid only, there exists a considerable amount of difference with the existing results in drag coefficient and vorticity at the cylinder surface; this may be because of the lack of the grid clustering to the surface that is a inevitable requirement to resolve the viscous flow. However, numerical results on the mixed grid show good agreements with the earlier computations and experimental data.

• 한국추진공학회지
• /
• 제21권5호
• /
• pp.71-79
• /
• 2017

액체-기체의 2상 유동에 대한 수치해석 기법을 연구하였다. 비압축성 방정식에는 가상 압축성 기법을 적용하였으며 LS와 VOF를 합친 CLSVOF 기법을 적용하여 액체-기체 경계면을 추적하였다. CLSVOF의 격자 의존도를 파악하기 위해 h=1/64, 1/128, 1/128, 그리고 1/160의 격자로 Zalesak's disk 문제와 액체의 3차원 변형 문제의 수치해석을 실시했으며 격자가 최대 보존 오차에 미치는 영향을 조사하였다. 비압축성 2상 유동 방정식을 적용하여 Rayleigh-Taylor 불안정성에 대한 수치해석을 실시하였고 밀도 차에 의한 액체 표면 불안정성이 나타난다는 것을 알 수 있었다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.389-396
• /
• 2010

In this paper, a new computational code was developed using Chorin's artificial compressibility method to solve the two-dimensional incompressible Navier-Stokes equations. In spatial derivatives, Roe's flux difference splitting was used for the inviscid flux, while central differencing was used for the viscous flux. Furthermore, AF-ADI with dual time stepping method was implemented for accurate unsteady computations. Two-equation turbulence models, Menter's $k-{\omega}$ SST model and Coakley's $q-{\omega}$ model, hae been adopted to solve high-Reynolds number flows. A number of numerical simulations were carried out for steady laminar and turbulent flow problems as well as unsteady flow problem. The code was verified and validated by comparing the results with other computational results and experimental results. The results of numerical simulations showed that the present developed code with the artificial compressibility method can be applied to slve steady and unsteady incompressible flows.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 춘계 학술대회논문집
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.