• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.397-403
• /
• 2010

The Immersed boundary method(IBM) is one of CFD techniques which can simulate flow field around complex objectives using simple Cartesian grid system. In the previous studies the IBM has mostly been implemented to fractional step method based Navier-Stokes solvers. In these cases, pressure buildup near IB was found to occur when linear interpolation and stadard mass conservation is used and the interpolation scheme became complicated when higher order of interpolation is adopted. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. Bi-linear and quadratic interpolation equations were formulated by using only geometric information of boundary to reconstruct velocities near IB. Flow around 2D circular cylinder at Re=40 and 100 was solved by using these formulations. It was found that the pressure buildup was not observed even when the bi-linear interpolation was adopted. The use of quadratic interpolation made the predicted aerodynamic forces in good agreement with those of previous studies.

• Proceeding of EDISON Challenge
• /
• 2016.11a
• /
• pp.1-3
• /
• 2016

대형 화물 트럭(heavy-duty truck)은 화물 적재에 용이하지만, 공기역학적으로 불리한 형상을 가진다. 이러한 단점을 극복하고자 대형 화물 트럭에는 공기저항력(aerodynamic drag)을 줄일 수 있는 여러 가지 장치가 달려있다. 본 논문에서는 디플렉터(deflector) 형상이 항력 감소에 어떠한 영향을 주며, 평판 형태와 굴곡진 형태의 디플렉터 형상에 대한 항력 계수 비교를 EDISON-CFD를 활용하여 비교하였다. 해석 결과, 측풍(side-wind)의 영향을 무시하며 차량 속도 95 km/h로 등속을 유지하는 조건에서 평판 형태의 Model 1과 바깥으로 굴곡진 Model 2에서 전체 항력 계수가 낮게 나타났다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.2
• /
• pp.89-107
• /
• 2017

In this paper, we have proposed a modified Marker-And-Cell (MAC) method to investigate the problem of an unsteady 2-D incompressible flow with heat and mass transfer at low, moderate, and high Reynolds numbers with no-slip and slip boundary conditions. We have used this method to solve the governing equations along with the boundary conditions and thereby to compute the flow variables, viz. u-velocity, v-velocity, P, T, and C. We have used the staggered grid approach of this method to discretize the governing equations of the problem. A modified MAC algorithm was proposed and used to compute the numerical solutions of the flow variables for Reynolds numbers Re = 10, 500, and 50000 in consonance with low, moderate, and high Reynolds numbers. We have also used appropriate Prandtl (Pr) and Schmidt (Sc) numbers in consistence with relevancy of the physical problem considered. We have executed this modified MAC algorithm with the aid of a computer program developed and run in C compiler. We have also computed numerical solutions of local Nusselt (Nu) and Sherwood (Sh) numbers along the horizontal line through the geometric center at low, moderate, and high Reynolds numbers for fixed Pr = 6.62 and Sc = 340 for two grid systems at time t = 0.0001s. Our numerical solutions for u and v velocities along the vertical and horizontal line through the geometric center of the square cavity for Re = 100 has been compared with benchmark solutions available in the literature and it has been found that they are in good agreement. The present numerical results indicate that, as we move along the horizontal line through the geometric center of the domain, we observed that, the heat and mass transfer decreases up to the geometric center. It, then, increases symmetrically.

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

Numerical analysis of two dimensional incompressible laminar free jet flow was carried out by using finite difference SMAC scheme. Flow characteristics of free jet flow such as jet width, similarity of jet velocity and hypothetical origin were investigated for different Reynolds numbers of Re=30 and 100. The reliability of predictions were confirmed by comparison with exact solution. Non-dimensional velocity distribution showed similarity of jet flow velocity after initial region. In the region of laminar flow, the location of hypothetical origin becomes more distant with Reynolds number.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.8
• /
• pp.653-662
• /
• 2007

Due to the difficulty raised from the coupling of cavitation modeling with turbulent flow, numerical simulation for two phase flow remains one of the challenging issues in the society. This research focuses on the development of numerical code to deal with incompressible two phase flow around 2D hydrofoil by combing the cavitation model suggested by Kunz et al. with $k-{\varepsilon}$ turbulent model. The simulation results are compared to experimental data to verify the validity of the developed code. Also, the comparison of the calculation results is made with LES results to evaluate the capability of $k-{\varepsilon}$ turbulence model. The calculation results show very good agreement with experimental observations even though this code can not grasp the small scaled bubbles in the calculation wheres LES can hold the real physics. This code will be extended to 3D compressible two phase flow for the study on the fluid dynamics in the inducers and impellers.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.44-53
• /
• 2012

Immersed boundary method(IBM) is a numerical scheme proposed to simulate flow field around complex objectives using simple Cartesian grid system. In the previous studies, the IBM has mostly been implemented to fractional step method based Navier-Stokes solvers. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. The weight coefficients of the bi-linear and quadratic interpolation equations were formulated by using only geometric information of boundary to reconstruct velocities near IB. Flow around 2D circular cylinder at Re=40 and 100 was solved by using these formulations. It was found that the pressure buildup was not observed even when the bi-linear interpolation was adopted. The use of quadratic interpolation made the predicted aerodynamic forces in good agreement with those of previous studies. For an analysis of moving boundary, we smulated an oscillating circular cylinder with Re=100 and KC(Keulegan-Carpenter) number of 5. The predicted flow fields were compared with experimental data and they also showed good agreements.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.10
• /
• pp.3272-3281
• /
• 1996

The flow analysis method which had been developed for the numerical calculation of 3-dimensional, incompressible and turbulent flow within an axial compressor was extended to the flow field within centrifugal impeller. In this method based on the SIMPLE(Semi Implicit Method Pressure Linked Equations) algorithm, the coordinate transformation was adopted and the standard k-.epsilon. model using wall function was used for turbulent flow analysis. The calculated flow fields have agreed very well with measurement results. Especially, 3-dimensional and viscous flow characteristics including secondary flows, jet-wake flow and decreased pressure rise along impeller passage, which can't be predicted by inviscid Q3D calculation were predicted very reasonably.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.55-61
• /
• 2010

A dependence of weighting parameter in preconditioning method for solving low Mach number flow with incompressible flow nature is investigated. The present preconditioning method employs a finite-difference method applied Roe‘s flux difference splitting approximation with the MUSCL-TVD scheme and 4th-order Runge-Kutta method in curvilinear coordinates. From the computational results of benchmark flows through a 2-D backward-facing step duct it is confirmed that there exists a suitable value of the weighting parameter for accurate and stable computation. A useful method to determine the weighting parameter is introduced. With this method, high accuracy and stable computational results were obtained for the flow with low Mach number in the range of Mach number less than 0.3.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.13 no.1
• /
• pp.164-169
• /
• 2010

Proper Orthogonal Decomposition (POD) is applied to the analysis of 2-dimensional cylinder wake flow. Time histories of flow variables were obtained by the incompressible CFD analysis. By using the method of snapshots the correlation matrix was constructed, and then eigenvalues, POD modes and time coefficients were calculated. Finally the flow field was reconstructed by using a few of the lower POD modes, and compared to the original ones.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.825-830
• /
• 2004

A numerical analysis based on two-dimensional and three-dimensional incompressible Navier-Stokes equations has been carried out for double-circular-arc (DCA) compressor cascades. Two types of double-circular-arc cascades were used in this analysis. The appropriate turbulence model for compressor analysis was selected among the conventional turbulence models such as Baldwin-Lomax, k-$\varepsilon$ and k-$\varepsilon$ models. The results of current study were compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE/PWIM algorithm for collocated grid and hybrid scheme for the convective terms were the main features of numerical tools. As commonly known, turbulence modeling is very important for the prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. For selection of turbulence model, 2-D analysis was performed. And then, k-$\varepsilon$ turbulence model with wall function chosen as the reasonable turbulence model for 3-D calculation was used to increase the efficiency of computation times. A reasonable result of 3-D flow pattern passing through the double-circular-arc cascade was obtained.