• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.483-487
• /
• 2011

In the present study, a CFD program based on a finite volume method was developed by using an unstructured polyhedral grid system for the accurate simulation with the complex geometry of computational domain. To simulate the transient flow induced by the moving solid object, the program used a fractional step method and a ALE (Algebric Lagrangian-Eulerian) method. The grid deformation for the moving of solid object were performed with a spring analogy based on the center coordinate of each computational grid. To verify the present program with these methodologies, the numerical results of the flow around the fixed and oscillating circular cylinder were compared with the previous numerical results.

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

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.380-388
• /
• 2010

A thermal-hydraulic code, named CUPID, has been developed for the analysis of transient two-phase flows in nuclear reactor components. A two-fluid three-field model was used for steam-water two-phase flows. To obtain numerical solutions, the finite volume method was applied over unstructured cell-centered meshes. In steam-water two-phase flows, a phase change, i.e., evaporation of condensation, results in a great change in the flow field because of substantial density difference between liquid and vapor phases. Thus, two-phase flows are very sensitive to the local pressure that determines the phase change. This in turn puts emphasis on the accurate evaluation of local pressure gradient. This paper presents a new numerical scheme to evaluate the pressure gradient at cell centers on unstructured meshes. The results of the new scheme for a simple test function a gravity-driven cavity, and a wall boiling two-phase flow are compared with those of the previous schemes in the cupid code.

• 한국전산유체공학회지
• /
• 제17권1호
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• 대한기계학회논문집B
• /
• 제26권5호
• /
• pp.648-657
• /
• 2002

This paper presents a numerical analysis on the solidification characteristics in twin-roll strip casting. Unstructured fixed-grid system was employed to deal with phase change. Melting of pure gallium was analyzed to confirm the validity of present program in both structured and unstructured grid systems. An algorithm for simultaneous calculation of the temperature in the roll and the molten metal pool was developed. The flow field in the pool and heat transfer features between pool and roll were shown. The effect of process parameters was also studied. Since the geometry of the molten metal Pool significantly deforms along the casting direction, unstructured grid system is more efficient. The unstructured grid system gives almost the same accuracy, even though the number of grids is only 60% of the structure done.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.311-317
• /
• 2011

The present paper deals with the continuous works of extending the multi-dimensional limiting process (MLP) for compressible flows, which has been quite successful in finite volume methods, into discontinuous Galerkin (DG) methods. From the series of the previous, it was observed that the MLP shows several superior characteristics, such as an efficient controlling of multi-dimensional oscillations and accurate capturing of both discontinuous and continuous flow features. Mathematically, fundamental mechanism of oscillation-control in multiple dimensions has been established by satisfaction of the maximum principle. The MLP limiting strategy is extended into DG framework, which takes advantage of higher-order reconstruction within compact stencil, to capture detailed flow structures very accurately. At the present, it is observed that the proposed approach yields outstanding performances in resolving non-compressive as well as compressive flaw features. In the presentation, further numerical analyses and results are going to be presented to validate that the newly developed DG-MLP methods provide quite desirable performances in controlling numerical oscillations as well as capturing key flow features.

• 한국전산유체공학회지
• /
• 제13권2호
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.

• 한국항공우주학회지
• /
• 제39권12호
• /
• pp.1097-1106
• /
• 2011

제자리 비행하는 헬리콥터 전기체 형상으로부터 발사되는 유도무기 운동을 모사하기 위해 비정렬 중첩격자기법을 6자유도 운동 기법과 연계하여 비정상 해석을 수행하였다. 주로 터, 꼬리로터, 그리고 유도무기를 포함하는 동체로 구성되는 헬리콥터 전기체 형상에 대한 비정상 해석을 수행하였으며, 6자유도 운동 기법을 이용하여 헬리콥터에서 발사되는 유도 무기의 초기 발사 거동을 예측하였다. 발사위치 및 유도무기의 추력크기에 따른 발사거동을 비교하였으며, 주로터에서 발생하는 비정상 내리흐름에 따른 유도무기의 비행 특성을 분석하였다.

• 한국항공우주학회지
• /
• 제31권2호
• /
• pp.1-9
• /
• 2003

본 연구에서는 비정렬 격자를 이용한 병렬화된 3차원 직접모사법을 개발하였다. 모사 입자의 추적에는 선형 형상 함수를 3차원으로 확장하여 적용하였다. 프로그램에서는 효율적인 병렬계산을 위하여 모사 입자의 수에 기초한 영역 재분할 방식을 이용하여 프로세서간의 작업 균형을 유지하였다. 밀도 차가 큰 두 개 이상의 유동이 동시에 존재할 경우를 효과적으로 모사하기 위해서 입자에 가중치를 적용하는 기법을 도입하였다. 프로그램의 검증을 위해서 3차원 삼각 날개 유동을 해석하여 실험치나 다른 계산 결과와 비교하였다. 또한 고도 100km의 희박한 환경헤서 우주발사체 탑재부 주의 유동을 해석하여 노즐 전방으로 고온의 풀룸이 역류하여 비행체 표면에 영향을 미칠 수 있는 가능성을 타진하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.524-529
• /
• 2010

In the present study, a CFD program is developed for the Fluid-Structure Interaction(FSI) analysis. The non-staggered, non-orthogonal, and unstructured grid system was also used to handle the complicated geometries in the program. In order to validate the capabilities of the developed CFD program, various models are investigated by using unstructured and nonorthogonal meshes. The predicted results are a good agreement with analytic solution, experimental data and commercial software. And also PISO algorithm is applied for transient flow analysis. The cyclic boundary condition and baffle cell are developed in order to improve the effectiveness of the calculation for complex geometry.