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

• KIEAE Journal
• /
• 제16권1호
• /
• pp.81-87
• /
• 2016

Purpose: Heat transfers phenomena are described by the second order partial differential equation and its boundary conditions. In a three-dimensional structure of a building, the heat transfer phenomena generally include more than one material, and thus, become complicate. The analytic solutions are useful to understand heat transfer phenomena, but they can hardly be applied in engineering or design problems. Engineers and designers have generally been forced to use numerical methods providing reliable results. Finite volume methods with the unstructured grid system is only the suitable means of the analysis for the complex and arbitrary domains. Method: To obtain an numerical solution, a discretization method, which approximates the differential equations, and the interpolation methods for temperature and heat flux between two or more materials are required. The discretization methods are applied to small domains in space and time, and these numerical solutions form the descretized equations provide approximated solutions in both space and time. The accuracy of numerical solutions is dependent on the quality of discretizations and size of cells used. The higher accuracy, the higher numerical resources are required. The balance between the accuracy and difficulty of the numerical methods is critical for the success of the numerical analysis. A simple and easy interpolation methods among multiple materials are developed. The linear equations are solved with the BiCGSTAB being a effective matrix solver. Result: This study provides an overview of discretization methods, boundary interface, and matrix solver for the 3-dimensional numerical heat transfer including two materials.

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

본 연구에서는 비정렬격자계를 사용하여 덮개 꼬리 로터의 제자리 비행에 대한 압축성, 비정성 유동을 해석하였다. 유동계산을 위한 수치적 기법으로는 셀 중심에 기초한 유한체적법과 내재적인 시간적분법을 사용하였다. 계산은 로터의 한 블레이드에 대해 수행되었으며, 블레이드와 블레이드 사이에는 주기적 경계조건을 설정하였다. 덮개가 없는 로터 형상에 대한 성능은 실험 결과와 잘 일치함을 보였다. 덮개를 포함하는 로터 형상에 대한 계산은 비교된 실험 형상의 불확실성을 고려하여 추력이 일치하는 피치각을 가지는 경우에 대해 수행하였으며, 자세한 유동은 실험결과와 잘 일치함을 확인할 수 있었다. 그 결과로부터 본 방법이 블레이드 끝단간극을 포함하는 복잡한 3차원 덮개 꼬리 로터 형상 해석에 매우 유용하게 사용될 수 있음을 알 수 있었다.

• 대한기계학회논문집B
• /
• 제30권7호
• /
• pp.603-611
• /
• 2006

This paper presents the response surface optimization method using three-dimensional Wavier-Stokes analysis to optimize the blade shape of an axial flow fan. Reynolds-averaged Wavier-Stokes equations with $k-{\epsilon}$ turbulence model are discretized with finite volume approximations using the unstructured grid. Regression analysis is used for generating response surface, and it is validated by ANOVA and t-statistics. Four geometric variables, i.e., sweep and lean angles at mean and tip respectively were employed to improve the efficiency. The computational results are compared with experimental data and the comparisons show generally good agreements. As a main result of the optimization, the total efficiency was successfully improved. Also, detailed effects of sweep and lean on the axial flow fan are discussed.

• 한국전산유체공학회지
• /
• 제9권4호
• /
• pp.48-54
• /
• 2004

Vortex-shedding flows past a circular cylinder for 200≤ Re ≤ 5000 are numerically simulated with the PowerCFD code, using a finite volume method and an unstructured grid system, developed by the author. The simulation is peformed by solving the unsteady 2-D Wavier-Stokes equations with both no model and turbulence model. The resulting Reynolds number dependence of the Strouhal number and of the drag and lift coefficients is compared with both experiments and previous numerical results. It is found that, in the range of 200≤ Re ≤ 5000 the calculation method with a turbulence model is capable of producing reasonably more accurate results than that with no model for the main practically relevant parameters such as Strouhal number, drag and lift coefficients.

• 한국전산유체공학회지
• /
• 제9권4호
• /
• pp.41-47
• /
• 2004

The paper presents numerical simulations of laminar vortex-shedding flows past a circular cylinder for Re ≤ 500. The simulations are performed by solving the unsteady 2-D Navier-Stokes equations with a finite volume method using unstructured grid system. The resulting Reynolds number dependence of the Strouhal number and of the drag and lift coefficients is compared with experiments and with previous numerical results, showing good agreement. It is found that, for the truly laminar Reynolds number range the present calculation method described is capable of producing reasonably accurate results for the main practically relevant parameters such as Strouhal number, drag and lift coefficients.

• 한국전산유체공학회지
• /
• 제11권4호
• /
• pp.39-47
• /
• 2006

In this paper, the authors developed a multidimensional interpolation method inside a finite volume cell in the computation of high-order accurate numerical flux such as the fifth order WEND (weighted essentially non-oscillatory) scheme. This numerical method starts from a simple Taylor series expansion in a proper spatial order of accuracy, and the WEND filter is used for the reconstruction of sharp nonlinear waves like shocks in the compressible flow. Two kinds of interpolations are developed: one is for the cell-averaged values of conservative variables divided in one mother cell (Type 1), and the other is for the vertex values in the individual cells (Type 2). The result of the present study can be directly used to the cell refinement as well as the convective flux between finer and coarser cells in the Cartesian adaptive grid system (Type 1) and to the post-processing as well as the viscous flux in the Navier-Stokes equations on any types of structured and unstructured grids (Type 2).

• 한국방재학회 논문집
• /
• 제7권3호
• /
• pp.69-78
• /
• 2007

본 연구에서는 홍수가 범람하였을 때 제내지에서 발생하는 동역학적 거동을 정확히 모의하기 위해, 시간에 따른 제방붕괴 양상을 고려하여 제방붕괴시 제내지로 유입되는 월류량을 정확하게 산정하였다. 2차원 비선형 천수방정식을 지배방정식으로 사용하였으며, 비구조적 격자계가 적용된 유한체적법을 이용하여 제방붕괴를 모의하였다. 제방붕괴시 발생하는 충격파 흐름을 해석하기 위하여 HLLC approximate Rimann solver를 사용하였고, 수치진동을 제어하기 위해 TVD 제한자를 사용하는 WAF(Weighted Averaged Flux) 기법을 사용하였다. 또한 생성항은 연산자 분리기법을 이용하여 비 물리적인 결과가 나오지 않도록 하였다. 먼저 본 모형을 댐붕괴 문제에 적용하여 댐 붕괴시 발생하는 자유수면 변위를 계산하였으며, 경사식 방파제의 월류량을 산정하여 기존 실험결과와 비교 검증하였다. 그 결과 충격파를 잘 모의하고 있었으며, 월류량 또한 기존 실험결과와 일치하였다. 또한 제방 붕괴시 발생하는 흐름에 대해 높이와 폭을 각각 시간에 따른 함수로 가정하여 적용하였다. 제방붕괴 유형에 따른 월류량을 각각 비교한 결과, 제방이 갑작스럽게 붕괴된 경우에서의 월류량이 점진적으로 붕괴되는 조건에서의 월류량보다 크게 산정됨을 알 수 있었다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2007년도 춘계학술대회B
• /
• pp.2857-2861
• /
• 2007

Flow and aerodynamic characteristics were analyzed numerically for a commercial passenger airplane, Boeing 747-400, flying in the cruising condition. The model geometry with 100:1 in scale was obtained by the photo scanning measurement with the maximum error of 1.4% comparing with the real airplane dimension. The three-dimensional inviscid steady compressible governing equations were solved by the finite volume method in the unstructured grid system. The convective terms were treated by the Crank-Nicholson and first-order upwind schemes. In the computational results, the strong wing-tip vortices were clearly observed and the pressure contours on the airplane surface were suggested. The lift and drag forces in the wing with engines increase by 1.49% and 3.9%, respectively compared with the case without engines. The aerodynamic forces were estimated quantitatively for each element which consists of the airplane.

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

In order to realistically represent the complex turbulence-chemistry interaction at the partially premixed turbulent lifted flames encountered in the gas turbine combustors, the combined conserved-scalar/level-set flamelet approach has been adopted. The parallel unstructured-grid finite-volume method has been developed to maintain the geometric flexibility and computational efficiency for the solution of the physically and geometrically complex flows. Special emphasis is given to the swirl effects on the combustion characteristics of the lean-premixed gas turbine combustor. Numerical results suggest that the present approach is capable of realistically simulating the combustion characteristics for the lean-premixed gas turbine engines and the lifted turbulent jet flame with a vitiated coflow.