• 대한기계학회논문집B
• /
• 제20권11호
• /
• pp.3589-3597
• /
• 1996

The Reynolds-averaged Navier-Stokes equations with the equations of the k-.epsilon. turbulence model are solved numerically in a general curvilinear system for a three-dimensional turbulent flow around an Ahmed body. The simulation is especially aimed at the evaluation of three finite differencing schemes for the convection term, which include the upwind differencing scheme(UDS), the second order upwind differencing scheme(SOU scheme) and the QUICK scheme. The drag coefficient, the velocity and pressure fields are found to be changed considerably with the adopted finite differencing schemes. It is clearly demonstrated that the large difference between computation and experiment in the drag coefficient is due to relatively high predicted values of pressure drag from both front part and vertical rear end base. The results also show that the simulation with the QUICK or SOU scheme predicts fairly well the flow field and gives more accurate drag coefficient than other finite differencing scheme.

• 설비공학논문집
• /
• 제11권1호
• /
• pp.73-80
• /
• 1999

This paper describes a numerical study of three-dimensional buoyant turbulent flow in a stairwell model with three convective differencing schemes, which include the upwind differencing scheme, the hybrid scheme and QUICK scheme. The Reynolds-averaged Navier-Stokes and energy equations are solved with a two-equation turbulence model. The Boussinesq approximation is used to model buoyancy terms in the governing equations. Three-dimensional predictions of the velocity and temperature fields are presented and are compared with experimental data. Three-dimensional simulations with each scheme have predicted the overall features of the flow fairly satisfactorily. A better agreement with experimental is achieved with QUICK scheme.

• 대한기계학회논문집B
• /
• 제27권7호
• /
• pp.995-1006
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes for large eddy simulation is evaluated by a dynamic analysis. Large eddy simulation of isotropic turbulence is performed with various dissipative and non-dissipative schemes to investigate the effect of numerical dissipation on the resolved solutions. It is shown by the present dynamic analysis that upwind schemes reduce the aliasing error and increase the finite differencing error. The existence of optimal upwind scheme that minimizes total numerical error is verified. It is also shown that the finite differencing error from numerical dissipation is the leading source of numerical errors by upwind schemes. Simulations of a turbulent channel flow are conducted to show the existence of the optimal upwind scheme.

• 대한기계학회논문집B
• /
• 제27권7호
• /
• pp.973-983
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes is evaluated for large eddy simulation of turbulent flow. Two turbulent flows are considered: turbulent channel flow at Re = 23000 and flow over a circular cylinder at Re = 3900. The effects of numerical dissipation on the finite differencing and aliasing errors and the subgrid-scale stress are investigated. It is shown through the simulations that compact upwind schemes are not suitable for LES, whereas the fourth order-compact centered scheme is a good candidate for LES provided that proper dealiasing of nonlinear terms is performed. The classical issue on the aliasing error and the treatment of nonlinear terms is revisited with compact difference schemes.

• 한국유체기계학회 논문집
• /
• 제1권1호
• /
• pp.72-80
• /
• 1998

Computational analysis of incompressible flows by numerically solving Navier-Stokes equations using multi-block finite volume method is conducted on a parallel computing system. Numerical algorithms adopted in this study $include^{(1)}$ QUICK upwinding scheme for convective $terms,^{(2)}$ central differencing for other terms $and^{(3)}$ the second-order Euler differencing for time-marching procedure. Structured grids are used on the body-fitted coordinate with multi-block concept which uses overlaid grids on the block-interfacing boundaries. Computational code is parallelized on the MPI environment. Numerical accuracy of the computational method is verified by solving a benchmark test case of the flow inside two-dimensional rectangular cavity. Computation in the axial compressor cascade is conducted by using 4 PE's md, as results, no numerical instabilities are observed and it is expected that the present computational method can be applied to the turbomachinery flow problems without major difficulties.

• 한국전산유체공학회지
• /
• 제1권1호
• /
• pp.64-70
• /
• 1996

Three dimensional turbulent flow fields around ships are simulated by a numerical method. Reynolds Averaged Navier-Stokes equations are used where Reynolds stresses are approximated by Baldwin-Lomax and Sub-Grid Scale(SGS) turbulence models. Body-fitted coordinate system is introduced to conform three dimensional ship geometries. The governing equations are discretized by a finite volume method. Temporal derivatives are approximated by the forward differencing and the convection terms are approximated by the QUICK or Kawamura scheme. The 2nd-order centered differencing is used for other spatial derivatives. Pressure and velocity fields are simultaneously iterated by the Highly Simplified Marker-And-Cell method. To verify the numerical method and turbulence models, flow fields around ships are simulated and compared to the experiments.

• 한국해양공학회지
• /
• 제11권4호
• /
• pp.141-151
• /
• 1997

This paper deals with a numerical study of flow-guider applied to controlling current in a bay. Two dimensional numerical model for tidal currents based on the depth averaged equation is developed and standard k-.epsilon. model is adopted to determine the turbulence diffusion. Equations are described in a generalized coordinate system to be implemented by non-staggered grid system and discretized by using finite volume method. Unsteady flow is simulated by fully implicit scheme. Hybrid scheme and central differencing are used to compute the convective terms and source terms, respectively. The tidal current in a rectangular bay is simulated and it gives satisfactory results. The realistic and distinct models of a large structure placed in bay are also exemplified with or without flow-guiders. The simulation results show that the flow-guider gives the residual tidal current in the bay by the different flux with respect to the direction of tidal current.

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

Three dimensional turbulent flow fields around ships are simulated by a numerical method. Reynolds Averaged Navier-Stokes equations are used where Reynolds stresses are approximated by Baldwin-Lomax and Sub-Grid Scale(SGS) turbulence models. Body-fitted coordinate system is introduced to conform three dimensional ship geometries. The governing equations are discretized by a finite volume method. Temporal derivatives are approximated by the forward differencing and the convection terms are approximated by the QUICK or Kawamura scheme. The 2nd-order centered differencing is used for other spatial derivatives. Pressure and velocity fields are simultaneously iterated by the Highly Simplified Marker-And-Cell method. To verity the numerical method and turbulence models, flow fields around ships are simulated and compared to the experiments.

• 대한기계학회논문집
• /
• 제12권5호
• /
• pp.1158-1174
• /
• 1988

본 연구에서는 계단형 벽면조건을 없게 하기 위해서 비직교 좌표계(non-orth- ogonal coordinate system)를 사용하여 수치해석하였다. 비직교 좌표계를 이용한 수 치해석의 예는 Thompson등이 Laplace방정식 혹은 Poisson방정식을 해석함으로써 비직 교 격자망을 구성한 바 있고, Fahgri와 Asako는 대수적 비직교 좌표변환으로 유한차분 방정식을 유도하여 비정규경계면을 갖는 관로에서의 유동특성을 해석하였으며 이재헌 과 이상렬은 Fahgri와 Asako의 방법을 비정규경계면을 갖는 밀폐공간내에서의 자연대 류의 수치해석에 적용한 바 있다. 본 해석에서도 Fahgri와 Asako의 변환법으로 유한 차분방정식을 유도하였는데, 이 방법을 사용할 경우 확대관의 경사벽면을 계단형으로 만들지 않고 유한차분방정식을 유도할 수 있어서 계단형 벽면으로 인한 해의 오차를 제거할 수 있다. Fig.2는 본 해석에서 사용한 비직교 격자망을 나타낸다.

• 한국전산유체공학회지
• /
• 제8권2호
• /
• pp.49-56
• /
• 2003

Three-dimensional endwall flow within a linear cascade passage of high performance turbine blade is simulated with a 3-D Navier-Stokes CFD code (MOSA3D), which is based on body-fitted coordinate system, pressure-correction and finite volume method. The endwall flow characteristics, including the development and generation of horseshoe vortex, passage vortex, etc. are clearly simulated, consistent with the generally known tendency. The effects of both turbulence model and convective differencing scheme on the prediction performance of endwall flow are systematically analyzed in the present paper. The convective scheme is found to have stronger effect than the turbulence model on the prediction performance of endwall flow. The present simulation result also indicates that the suction leg of the horseshoe vortex continues on the suction side until it reaches the trailing edge.