• 한국항해항만학회:학술대회논문집
• /
• 한국항해항만학회 2011년도 춘계학술대회
• /
• pp.222-224
• /
• 2011

본 연구에서는 유한체적법과 함께 비정렬 격자(Unstructured Grid) 를 채택한 SELFE 모델을 기반으로 운영되는 유출유 확산모델을 개발하였다. 모델의 적용성을 검토하기 위해 미극 오레건주 콜롬비아 강의 유출유 모델링에 적용하였으며 양호한 결과를 얻을 수 있었다. 앞으로 이 모델은 하천, 호수 및 해양을 포함한 다양한 수계에 적용이 가능할 것으로 기대된다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 추계학술대회
• /
• pp.910-915
• /
• 2003

A Navier-Stokes equation solver for incompressible viscous flows with free surface is developed and tested. This is based upon a fractional time step method and a non-staggered finite volume formulation for unstructured meshes. For time advancement scheme, Adams -Bashforth method for convective term and Crank-Nicolson method for diffusive term are applied. The interface between two fluids with different fluid properties is tracked with Piecewise Linear Interface Calculation(PLIC) Volume-of-Fluid(VOF) methods. Computational results are presented for some test problems: the broken dam, the sloshing in a rectangular tank, the filling of a cylindrical tank.

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

A field reconstruction scheme for a cell centered finite volume method on unstructured meshes is developed. Regardless of mesh quality, this method is exact within a machine accuracy if the solution is linear, which means it has full second order accuracy. It does not have any limitation on cell shape except convexity of the cells and recovers standard discretization stencils at structured orthogonal grids. Accuracy comparisons with other popular reconstruction schemes are performed on a simple example.

• 한국항공우주학회지
• /
• 제35권8호
• /
• pp.677-684
• /
• 2007

본 연구에서는 삼차원 점성 유동을 효율적으로 해석하기 위해 사면체, 프리즘, 피라미드를 포함하는 비정렬 혼합격자계를 기반으로 하는 유동해석코드를 개발하였다. 유동의 지배방정식은 격자점 중심의 유한체적법을 사용하여 공간차분회었으며, 제어테적은 메디안 듀얼(median-dual)방법으로 구성하였다. 난류유동 해석은 Spalart-Allmaras 난류모형과 연계하여 계산되었다. 개발된 해석코드의 정상 유동 검증을 위해 삼차원 날개에 대한 층류, 난류유동을 해석하였으며, 비정상 유동 검증을 위해 조화운동에 의해 진동하는 삼차원 날개에 대한 유동해석을 수행하였다.

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

A new anisotropic mesh refinement method is proposed. The new method is based on a simple second order interpolation error indicator. Therefore, it is methodologically direct and intuitive as compared with traditional anisotropic refinement strategies. Moreover, it does not depend on the mesh type. The error indicator is face-wisely calculated for all faces in a mesh and the cell refinement type is determined by the configuration of face markings with a given threshold. For the sake of simplicity, an application for a poisson equation on a triangle mesh is considered. The error field and resultant mesh refinement pattern are compared and effects of the threshold selection are discussed. Applying anisotropic refinement with various thresholds, we observed higher convergence rates than those in the uniform refinement cases.

• 대한기계학회논문집B
• /
• 제35권6호
• /
• pp.561-567
• /
• 2011

2 차원 비직교 형상에 대해 비정렬 삼각 격자를 이용하여 복사 열전달의 역해석을 수행하였다. 본 논문에서는 확산적으로 방사 및 반사를 하는 벽면으로 이루어진 형상 내부에 흡수 및 방사, 산란하는 매질이 채워져 있는 문제를 고려하였다. 유한체적법을 사용하여 복사전달 방정식을 계산하였고 이 때 얻은 입사복사량을 역해석의 측정 데이터로 사용하였다. 벽면의 방사율을 추정하기 위해 켤레구배법을 적용하였으며, 목적 함수를 최소화하는 과정을 통해 해를 구하였다. 측정값의 측정 오차가 추정 정확도에 미치는 영향을 살펴보았고, 비정렬 격자계의 성능을 확인하기 위해 정렬 격자계를 이용하여 얻은 결과와 비교해 보았다.

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

An extension of our recently developed locally linear reconstruction scheme to 2 dimensional incompressible flow solver is presented. The solver is based on a semi-implicit fractional step method in which the convective term is discretized by Adams-Bashforth method and the diffusion term by Crank-Nicolson method. Several numerical examples are tested to demonstrate the mesh type independent accuracy of the solver, which include decaying vortex flow, square cavity flow, and flow around a circular cylinder. The above examples are solved on quadrilateral or hybrid meshes. For all numerical examples, we obtained reasonable results.

• 한국분무공학회지
• /
• 제5권1호
• /
• pp.13-22
• /
• 2000

The present study has numerically analyzed the vaporization characteristics of fuel droplets in the high temperature convective flow field. The axisymmetric governing equations for mass, momentum, energy, and species are solved by an iterative and implicite unstructured finite-volume method. The moving boundary due to vaporization is handled by the deformable unstructured grid technique. The pressure-velocity coupling in the density-variable flows is treated by the SIMPLEC algorithm. In terms of the matrix solver, Bi-CGSTAB is employed for the numerically efficient and stable convergence. The n-decane is used as a liquid fuel and the initial droplet temperature is 300K. Computations are performed for the nonevaporating and evaporating droplets with the relative interphase velocity(25m/s). The unsteady vaporization process has been simulated up to the nondimensional time, 25. Numerical results indicate that the mathematical model developed in this study succesfully simulates the main features of the droplet vaporization process in the convective environment.

• 한국전산유체공학회지
• /
• 제18권2호
• /
• pp.41-48
• /
• 2013

Ice accretion on the solid surface is an importance factor in assessing the performance of aircraft and wind turbine blade. Changes in the external shape due to ice accretion can greatly deteriorate the aerodynamic performance. In this study, a three-dimensional upwind-type second-order positivity-preserving finite volume CFD scheme based on the unstructured mesh topology is developed to simulate two-phase flow in atmospheric icing condition. The code is then validated by comparing with NASA IRT experimental data on the sphere. The present results of the collection efficiency are found to be in close agreement with experimental data and show improvement near the stagnation region.

• 대한기계학회논문집
• /
• 제19권11호
• /
• pp.3031-3040
• /
• 1995

Two-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with the unstructured triangular meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. The convective fluxes are obtained by the Roe's flux difference splitting scheme using edge-based connectivities and higher-order differences are achieved by a reconstruction procedure. The time integration is based on an explicit four-stage Runge-Kutta scheme. Numerical procedures with local time stepping and implicit residual smoothing have been implemented to accelerate the convergence for the steady-state solutions. Comparisons with experimental data and other numerical results have proven accuracy and efficiency of the present unstructured approach.