• 한국전산유체공학회지
• /
• 제4권1호
• /
• pp.8-18
• /
• 1999

In devising a numerical approximation for the convective spatial transport of a fluid mechanical quantity, it is noted that the convective motion of a scalar quantity occurs in one-way, or from upstream to downstream. This consideration leads to a new scheme termed a pure upwind difference scheme (PUDS) in which an estimated value for a fluid mechanical quantity at a control surface is not influenced from downstream values. The formal accuracy of the proposed scheme is third order accurate. Two typical benchmark problems of a wall-driven fluid flow in a square cavity and a buoyancy-driven natural convection in a tall cavity are computed to evaluate performance of the proposed method. for comparison, the widely used simple upwind scheme, power-law scheme, and QUICK methods are also considered. Computation results are encouraging: the proposed PUDS sensitized to the convection direction produces the least numerical diffusion among tested convection schemes, and, notable improvements in representing recirculation of fluid stream and spatial change of a scalar. Although the formal accuracy of PUDS and QUICK are the same, the accuracy difference of approximately a single order is observed from the revealed results.

• 대한기계학회논문집B
• /
• 제24권9호
• /
• pp.1227-1238
• /
• 2000

This study examines the effects of the discretization schemes on numerical solutions of viscoelastic fluid flows. For this purpose, a temporally evolving mixing layer, a two-dimensional vortex pair interacting with a wall, and a turbulent channel flow are selected as the test cases. We adopt a fourth-order compact scheme (COM4) for polymeric stress derivatives in the momentum equations. For convective derivatives in the constitutive equations, the first-order upwind difference scheme (UD) and artificial diffusion scheme (AD), which are commonly used in the literature, show most stable and smooth solutions even for highly extensional flows. However, the stress fields are smeared too much and the flow fields are quite different from those obtained by higher-order upwind difference schemes for the same flow parameters. Among higher-order upwind difference schemes, a third-order compact upwind difference scheme (CUD3) shows most stable and accurate solutions. Therefore, a combination of CUD3 for the convective derivatives in the constitutive equations and COM4 for the polymeric stress derivatives in the momentum equations is recommended to be used for numerical simulation of highly extensional flows.

• Journal of Hydrospace Technology
• /
• 제2권2호
• /
• pp.57-67
• /
• 1996

The formulation for hybrid-QUICK scheme of convective transport terms in finite-volume calculation procedure is presented. Source terms are modified to apply the hybrid-QUICK scheme. Test calculations are performed for wall-driven cavity flow at Re=$10_2$, $10_3$, and $10_4$. These include the evaluation of boundary conditions approximated by third-order finite difference scheme. The stable and converged solutions are obtained without unsteady terms in the momentum equations. The results using hybrid-QUICK scheme show no difference with those using hybrid scheme at low Re ($=10_2$) and are better at higher Re ($10_3$, and $10_4$).

• 대한토목학회논문집
• /
• 제7권4호
• /
• pp.147-157
• /
• 1987

대류 분산 모형의 유한 차분 방법으로 양해법, Bresler방법, 음해법, upstream차분법과 Chaudhari방법등을 선택하여 각 차분법들의 특성을 규명하고 수치실험을 통하여 이들의 효율적인 사용 방안을 제시하였다. 비교 분석 결과 Chaudhari방법은 수치 분산 현상에 가장 둔감한 반면 조건부 안정이고, Bresler방법은 overshooting에 민감한 반면 무조건 안정이라는 특성이 있다. 분산이 지배적인 흐름에서는 양해법이 가장 정확하고, 대류가 지배적인 흐름에서는 Chaudhari방법이 가장 정확하다. 계산 시간(CPU)은 양해법 또는 Chaudhari방법이 비슷하게 가장 작고 Bresler방법이 항상 가장 크다.

• Journal of Advanced Marine Engineering and Technology
• /
• 제20권3호
• /
• pp.101-109
• /
• 1996

Two-dimensional lid-driven closed flows within square cavity were studied numerically for four Reynolds numbers : $10^4$, 3$\times10^4$, 5$\times10^4$ and 7.5$\times10^4$. A convective difference scheme to maintain the same spatial accurary by irregular grid correction is adopted by applying the interior division principle. Grid number is $80\times80$and its minimum size is about 1/400 of the cavity height. At Re=$10^4$, periodic migration of small eddies appearing in corner separation region and its temporal sinusoidal fluctuation are represented. At three higher Reynolds numbers(3$\times10^4$, 5$\times10^4$ and 7.5$\times10^4$), an organizing structure of four consecutive vorticles at two lower corners is revealed from time-mean flow patterns. But, instantaneous flow characteristics show very random unsteady fluctuation mainly due to the interaction between rotating shed vortices and stationary eddies within the corners.

• 대한기계학회논문집
• /
• 제17권2호
• /
• pp.467-474
• /
• 1993

Viscous flow and heat transfer phenomena in an axisymmetric cylinder which models a diesel engine have been numerically studied. In order to search for a way to minimize numerical diffusion, the effectiveness and the appropriateness of two selected numerical schemes for convective terms in the governing equations have been tested. They are Linear Upwind Difference Scheme and Hybrid Scheme. Using a standard k-.epsilon. turbulence model, the calculation has been carried out basically up to 180.deg. of crank angle. As a result, it was shown from comparison with previous experimental data that Linear Upwind Difference Scheme is less influenced than Hybrid Scheme by the numerical diffusion and it was suggested that these effects of numerical diffusion can be more significant than those due to turbulence modeling.

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

• 한국자동차공학회논문집
• /
• 제6권5호
• /
• pp.86-96
• /
• 1998

A numerical study has been carried out of three-dimensional turbulent flows around a MIRA reference vehicle model both with and without wheels in computation. Two convective difference schemes with two k-$\varepsilon$ turbulence models are evaluated for the performance such as drag coefficient, velocity and pressure fields. Pressure coefficients along the surfaces of the model are compared with experimental data. The drag coefficient, the velocity and pressure fields are found to change considerably with the adopted finite difference schemes. Drag forces computed in the various regions of the model indicate that design change decisions should not rely just on the total drag and that local flow structures are important. The results also indicate that the RNG model with the QUICK scheme predicts fairly well the tendency of velocity and pressure fields and gives more reliable drag coefficient rather than the other cases.

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

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

• 한국전산유체공학회지
• /
• 제3권2호
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.