• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.57-70
• /
• 2008

A high-order accurate Euler flow solver based on a discontinuous Galerkin finite-element method has been developed for the numerical simulations of blade-vortex interaction phenomena on unstructured meshes. A free vortex in freestream was investigated to assess the vortex-preserving property and the accuracy of the present flow solver. Blade-vortex interaction problems in subsonic and transonic freestreams were simulated by adopting a multi-level solution-adaptive dynamic mesh refinement/coarsening technique. The results were compared with those of other numerical and experimental methods. It was shown that the present discontinuous Galerkin flow solver can preserve the vortex structure for significantly longer vortex convection time and can accurately capture the complex unsteady blade-vortex interaction flows, including generation and propagation of acoustic waves.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.57-70
• /
• 2008

A high-order accurate Euler flow solver based on a discontinuous Galerkin finite-element method has been developed for the numerical simulations of blade-vortex interaction phenomena on unstructured meshes. A free vortex in freestream was investigated to assess the vortex-preserving property and the accuracy of the present flow solver. Blade-vortex interaction problems in subsonic and transonic freestreams were simulated by adopting a multi-level solution-adaptive dynamic mesh refinement/coarsening technique. The results were compared with those of other numerical and experimental methods. It was shown that the present discontinuous Galerkin flow solver can preserve the vortex structure for significantly longer vortex convection time and can accurately capture the complex unsteady blade-vortex interaction flows, including generation and propagation of acoustic waves.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.25-30
• /
• 2003

The numerical simulations with unstructured mesh by cell-centered and vertex-centered approaches were performed for the quadrilateral and triangular meshes. For the 2-D incompressible supersonic vortex flow, the simulation results and the analytic solution were compared and the accuracy was assessed. The calculation efficiency was measured by the parameter defined by the consumed CPU time multiplied by absolute error, As a results, equilateral triangular mesh yielded the best accuracy and efficiency among the tested meshes.

• Proceeding of EDISON Challenge
• /
• 2013.04a
• /
• pp.354-359
• /
• 2013

일반적으로 건축물의 설계시 풍동 실험을 통한 풍환경의 평가를 수행하고 있으며, 이는 환경 영향 평가법에서 정한 건축 사업 시행 시 수반되어야 할 자연환경, 생활환경 그리고 사회경제환경의 영향 평가의 일환으로 실시되고 있다. 그러나, 풍동 실험의 경우 여러 가지 현실적 제약조건으로 설계와 실험의 피드백 (Feedback)이 원활하지 못하며, 특히 대상 건축물이 공장과 같이 대기 오염원이 되는 경우 실험은 더욱 어려운 형편이다. 이에 대한 보완책으로 전산 유체 역학을 이용한 건축물의 풍압 해석에 의한 풍하중 추정이나 인접 지형-지물의 영향을 고려한 건축물 주위의 풍환경 평가가 있다. 전산 모사에 의해 풍동 실험의 미비점을 보완하고, 보다 상세한 정보를 확보함으로써 건축물의 구조적 안전성의 증대와 환경 피해 감소를 기할 수 있다. 그러나 복잡한 지형-지물이나 건축물 주위의 풍환경에 대한 전산 모사는 주로 두 가지의 기술적 어려움을 수반하게 되다. 그 중 하나는 고정 경계면을 이루는 형상의 복잡성으로 인해 기존에 많이 이용하고 있는 Body-fitted 격자계를 이용하는 경우, 격자 생성 과정이 매우 복잡하고 어려울 뿐 만 아니라 생성된 격자가 주로 비정렬 (unstructured) 특성을 갖게 되어 수치해석 과정의 효율을 저하시키는 요인이 되며, 격자의 형상도 수치해석의 수렴성을 저하시키는 예가 많다. 다른 어려움으로 풍환경은 전형적인 난류 유동장으로서 난류의 전산 해석은 아직도 해결하지 못한 부분이 많다는 점이다. 이에 본 논문에서는 복잡한 지형-지물이나 건축물의 풍하중과 풍환경의 전산 모사 기술 확보를 위하여 수행중인 연구의 일환으로 물체 형상의 기하학적 복잡성의 극복을 위한 가상경계법 (Immersed Boundary Method)과 난류 유동장의 물리적 엄밀성을 높이기 위한 다와동 모사 (Large Eddy Simulation)을 이용한 물체 형상과 무관한 유동장 해석 기술 개발에 대하여 다루고자 한다. 먼저 최근에 유동 해석에 이용되는 방법인 가상경계법(IBM)은 물체를 포함한 전체 전산 영역을 직교 좌표계에 의해 이산화하고, 유동장내 존재하는 물체의 표면에서의 점착 조건을 만족시키기 위하여 지배 방정식에 적절한 외력을 추가로 고려하는 방법이다. 본 연구에서는 가상경계법을 이용하여 경계층에 위치한 건물 형상의 각진 물체 주위 사이에 형성되는 공동 내부의 비정상 유속 및 압력에 대한 전산 해석을 수행하고, 풍상측 전면에 형성되는 경계층에 의한 영향을 분석하였다.

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.68-73
• /
• 2007

Computational Fluid Dynamics (CFD) approach is now playing an important role in the engineering process in these days. Generating proper grid system in time for the region of interest is prerequisite for the efficient numerical calculation of flow physics using CFD approach. Grid generation is, however, usually considered as a major obstacle for a routine and successful application of numerical approaches in the engineering process. CFD approach based on the unstructured grid system is gaining popularity due to its simplicity and efficiency for generating grid system compared to the structured grid approaches, especially for complex geometries. In this paper an automated triangular surface grid generation using CAD(Computer Aided Design) surface data is proposed. According to the present method, the CAD surface data imported in the STL(Stereo-lithography) format is processed to identify feature edges defining the topology and geometry of the surface shape first. When the feature edges are identified, node points along the edges are distributed. The initial fronts which connect those feature edge nodes are constructed and then they are advanced along the CAD surface data inward until the surface is fully covered by triangular surface grid cells using Advancing Front Method. It is found that this approach can be implemented in an automated way successfully saving man-hours and reducing human-errors in generating triangular surface grid system.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.5
• /
• pp.1-11
• /
• 2006

The effects of continuous blowing on flow control and stall suppression for flows over a NACA 0015 airfoil at low Reynolds numbers were numerically investigated through its parameter variation on unstructured meshes. The aerodynamic force and moment variations due to flow control were examined, along with the stall angle-of-attack change for stall suppression. The results showed that blowing with relatively strong jet increases lift at the cost of drag increment below stall angle. Continuous blowing delays flow stall when it is implemented near the leading edge. When the blowing jet was aligned along the flow direction on the airfoil, the favorable flow control effect was most significant below the stall angle of attack.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.211-217
• /
• 2009

The numerical program has been developed for the purpose of the complicate geometries application using CIP method. The non-staggered, non-orthogonal, and unstructured grid system can be also used for the various geometries in the program. For validating CIP solver, the lid-driven cavity flow and solitary wave propagation flow are carried out. Test results show a good agreement with the verified results. The dynamic solver was used for the behavior of moving body. Interface process between the two solvers is introduced. The research was performed on the flow problem around torpedo and log and the flow problem in a tank in order to analyze the three phase flow problem Although the comparison to the verified results was not quantitatively performed, the trend of the results was reasonable.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1581-1593
• /
• 2001

An efficient variable-reordering method for finite element meshes is used and the effect of variable-reordering is investigated. For the element renumbering of unstructured meshes, Cuthill-McKee ordering is adopted. The newsy reordered global matrix has a much narrower bandwidth than the original one, making the ILU preconditioner perform bolter. The effect of variable reordering on the convergence behaviour of saddle point type matrix it studied, which results from P2/P1 element discretization of the Navier-Stokes equations. We also propose and test 'level(0) preconditioner'and 'level(2) ILU preconditioner', which are another versions of the existing 'level(1) ILU preconditioner', for the global matrix generated by P2/P1 finite element method of incompressible Navier-Stokes equations. We show that 'level(2) ILU preconditioner'performs much better than the others only with a little extra computations.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.18-24
• /
• 2003

In the present study, solution algorithms for the connotation of unsteady flows on an unstructured mesh me presented Dual time stepping is incorporated to achieve the 2-nd order temporal accuracy while reducing the linearization and the factorization errors associated with a linear solver. Hence, any time step can be used by only considering physical phenomena. Gauss-Seidel scheme is used to solve linear system of equations. Rigid motion and suing analogy method for moving mesh are all considered and compared. Special treatments of suing analogy for high aspect ratio cells are presented. Finally, numerical results for oscillating ing are compared with experimental data.

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.44-53
• /
• 2007

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed for a 3D component of a nuclear system code and a component-scale analysis tool. A two-fluid three-field model is used for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. An unstructured grid is adopted for realistic simulations of the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been applied to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the modified numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing very well.