• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.6
• /
• pp.1179-1194
• /
• 1992

QUICKER scheme has several attractive properties. However, under highly convective conditions, it produces overshoots and possibly some oscillations on each side of steps in the dependent variable when the flow is convected at an angle oblique to the grid line. Fortunately, it is possible to modify the QUICKER scheme using non-linear and linear functional relationship. Details of the development of polynomial upwinding scheme are given in this paper, where it is seen that this non-linear scheme has also third order accuracy. This polynomial upwinding scheme is used as the basis for the SHARPER and SMARTER schemes. Another revised scheme was developed by partial modification of QUICKER scheme using CDS and UPWIND schemes(QUICKUP). These revised schemes are tested at the well known bench mark flows, Two-Dimensional Pure Convection Flows in Oblique-Step, Lid Driven Cavity Flows and Buoyancy Driven Cavity Flows. For pure convection oblique step flow test problem, QUICKUP, SMARTER and SHARPER schemes remain absolutely monotonic without overshoot and oscillation. QUICKUP scheme is more accurate than any other scheme in their relative accuracy. In high Reynolds number Lid Driven Cavity Flow, SMARTER and SHARPER schemes retain lower computational cost than QUICKER and QUICKUP schemes, but computed velocity values in the revised schemes produced less predicted values than QUICKER scheme which is strongly effected by overshoot and undershoot values. Also, in Buoyancy Driven Cavity Flow, SMARTER, SHARPER and QUICKUP schemes give acceptable results.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.80-86
• /
• 1994

A numerical method for simulations of inviscid incompressible flow fields around a ship advancing on the free surface is developed. A body fitted coordinate system, generated by numerically solving elliptic type partial differential equations is used to conform the ship and free surface configurations. Three dimensional Euler equations transformed to the non-staggered body fitted coordinate system are discretised by finite difference method. Time and spatial derivatives are discretised by forward and centered differencings, respectively, and artificial dissipations are added to discretised convection terms for improvements of numerical stability. At each time steps, free surface elevations are recomputed to satisfy nonlinear free surface conditions. Poisson equations for pressure field are solved iteratively and the velocity field for next time step is extrapolated. To verify the developed numerical method, flow fields around a Wigley model are simulated(Fn=0.250-0.408) and compared with experimental data to show good agreements.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.19 no.4
• /
• pp.108-116
• /
• 2005

This paper presents finite element analyses solution in the travelling magnetic field problem. The travelling magnetic field problem is subject to convective-diffusion equation. Therefore, the solution derived from Galerkin-FEM with linear interpolation function may oscillate between the adjacent nodes. A simple model with Dirichlet, Neumann and Periodic boundary condition respectively, have been analyzed to investigate stabilities of solutions. It is concluded that the solution of Galerkin-FEM may oscillate according to boundary condition and element type, but that of Upwind-FFM is stable regardless boundary condition.

• Journal of the Korea Computer Graphics Society
• /
• v.15 no.1
• /
• pp.1-6
• /
• 2009

This paper presents a new CIP method for unstructured mesh to reduce the numerical dissipation. To reflect precise physical characteristics, CIP method updates both the physical quantity and the derivative information. The proposed method uses the Finite Volume Method(FVM) to solve the non-advection term of CIP equation. And we performed several experiments to improve the accuracy of third-order interpolation. Our result shows that our algorithm has less numerical dissipation than that of linear advection solver.

• Journal of Biosystems Engineering
• /
• v.10 no.1
• /
• pp.24-38
• /
• 1985

경사진 밀폐 공간에서 마주 보는 두 벽면의 온도 차로 인하여 발생되는 자연 대류 현상은 여러 공학 분야에서 볼 수 있는 중요한 열전달 현상으로서, 최근 들어 평판형 태양열 집열기를 설계하려는 사람들에게 많은 관심의 대상이 되고 있다. 평판형 태양열 집열기의 경우 덮개판으로 부터의 대류 열손실을 감소시킴으로서 집열 효율을 높일 수 있을 뿐만 아니라 사용목적에 따라 소형 집열기를 제작할 수 있어 경제적으로 유리하게 될 것이다. 밀폐된 공간에서 최초에 정지 상태에 있는 얇은 유체층을 하부에서 가열시켜 주면 열팽창 현상이 일어나고, 이것에 의한 부력이 점도나 열전도도 등의 안정화 요인을 극복할 수 있을 정도로 커지면 System이 불안정하게 되어 자연 대류 현상이 수반되며 이 때문에 열전달율이 급격히 증가하게 된다. 이러한 현상의 지배 방정식은 연립 비선형 편미분 방정식으로 특수한 경계 조건외에는 일반적으로 해석적 해를 구하기가 어렵기 때문에 실험적 연구가 많이 실시되어 왔고 이들 결과의 대부분은 전반적인 열전달 특성치만을 구하는데 집중되어 왔다. 본 연구에서는 수치 해석법의 하나인 유한 차분법을 도입하여 이차원으로 가정한 경사진 평판형 밀폐 공간에서의 자연 대류 현상의 지배 방정식을 유한 차분화시켜, $$2.74{\times}10^3\leq_-Gr\leq_-2.0{\times}10^6$$, Pr=0.73, $$15^{\circ}\leq_-a\leq_-150^{\circ}$$, 종횡비는 1, 2, 3, 5, 9에 대하여 정상 상태에서의 해를 구하면서 수치적으로 실험하였다. 본 연구에서 얻어진 결론을 요약하면 다음과 같다. (1) 해석적으로 구하기 어려운 경사진 밀폐 공간에서 자연대류현상의 지배 방정식을 유한 차분법으로 해결할 수 있으며, 대류항과 확산항을 따로 고려한 유한차분법이 효과적임을 확인하였다. (2) 저온과 고온 벽면에서의 온도를 각각 균일하게 놓고 단변으로 이루어진 벽면은 완전히 절연되어 있는 경우에 대하여 수치해를 구한결과, 이전의 해석적 및 실험적 결과와 일치하였으며, 시간의 경과에 따른 온도 및 유선의 변화를 현상학적으로 관찰할 수 있었다. (3) 평균 열전달 계수에 미치는 경사각의 효과를 살펴본 결과 종횡비가 1 인 경우 경사각이 $45^{\circ}$에서, 종횡비가 2, 3, 5, 9인 경우 경사각이 $60^{\circ}$에서 각각 평균 열전달 계수 최대치가 나타났다. (4) Ra수(Rayleigh number) 가 증가될수록, 경사각에 상관없이 평균 열전달 계수도 증가되었다. 한편 Ra수 및 경사각의 변화에 따라 종횡비가 증가될수록 평균 열전달 계수는 경사각이 $90^{\circ}$인 경우를 제외하고는 감소됨을 볼 수 있었다. 경사각이 $90^{\circ}$인 경우, 평균 열전달 계수는 종횡비가 2 인 곳에서 최대치를 얻을 수 있었으며, 종횡비가 계속 증가될수록 평균 열전달 계수는 점차 감소되어짐을 볼 수 있었다.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.4
• /
• pp.19-30
• /
• 2000

A finite-volume method is developed to solve turbulent flows around modern commercial hull forms with bow and stern bulbs. The RANS equations are solved. The cell-centered finite-volume method employs QUICK and central difference scheme for convective and diffusive flux discretization, respectively. The SIMPLEC method is adopted for the velocity-pressure coupling. The developed numerical methods are applied to calculate turbulent flow around KRISO 3600TEU container ship. Surface meshes are generated into five blocks: bow and stern bulbs, overhang, fore and afterbody. 3-D field grid system with O-H topology is generated using elliptic grid generation method. Surface friction lines and wake distribution at propeller plane is compared with experiment. The calculated results show that the present method can be used to predict flow around a modern commercial hull forms with bulbs.