• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.49-54
• /
• 1997

Rice's monotone streamline upwind finite element method, which was proposed to treat convection-dominated flows, is applied to the linear triangular element. An alignment technique of unstructured grids with given velocity fields is used to prevent the interpolation error produced in evaluating the convection term in the upwind method. The alignment of grids is accomplished by optimizing a target function defined with the inner-product of a properly chosen side vector in the element with the velocity field. Two pure advection problems are considered to demonstrate the superiorities of the present approach in solving the convection-dominated flow on the unstructured grid. Solutions obtained with aligned grids are much closer to the exact solutions than those with initial regular grids. The capability of the present approach in predicting the appearance of the secondary vortex in the laminar confined jet impingement is shown by comparing streamlines to those produced by SIMPLE on a highly stretched grid toward the impingement plate.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.29-37
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geometries on unstructured triangular meshes. The flux terms are discretized based on a cell-centered finite-volume formulation with the Roe's flux-difference splitting. The solution is advanced in time using an implicit backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with the Gauss-Seidel relaxation scheme. The effect of turbulence effects is approximated with a standard $k-{\varepsilon}$ two-equation model which is solved separately from the mean flow equations using the same backward-Euler time integration scheme. The triangular meshes are generated using an advancing-front/layer technique. Validations are made for flows over the NACA0012 airfoil and the Douglas 3-element airfoil. Good agreements are obtained between the numerical results and the experiment.

• Proceeding of EDISON Challenge
• /
• 2012.04a
• /
• pp.57-60
• /
• 2012

본 연구에서는 현재 국내외에서 연구 중인 Blended Wing Body(BWB) 항공기의 엔진흡기 유동을 해석하기 위해 익형과 비행조건을 변화시켜 가며 전산유체해석을 수행하였다. 엔진의 위치에 따라 엔진이 효율적으로 동작하기 위한 조건인 흡기에서의 유동 속도와 그 분산을 중심으로 해석한 결과 익형 표면에서는 경계층의 영향으로 엔진흡기에서 유동속도가 낮고, 속도분산이 높음을 확인할 수 있었다. 한편, 익형 아랫면에서는 높은 비행속도에서 속도분산이 급격히 증가하였다. 이를 통해, 해석에 사용한 익형이 BWB의 동체로 활용하기에 적합한 엔진흡기조건을 갖는지 판별하였다.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.4
• /
• pp.179-189
• /
• 1996

An unstructured grid finite-volume method has been applied to predict the linear and nonlinear attenuation characteristics of the expansion chamber silencer system. In order to achieve a grid flexibility and a solution adaptation for geometrically silencer system. In order to achieve a grid flexibility and a solution adaptation for geometrically complex flow regions associated with the actual silencers, the unstructured mesh algorithm in context with the node-centered finite volume method has been employed. The present numerical model has been validated by comparison with the analytical solutions and the experimental data for the acoustic field of the concentric expansion chamber with and without pulsating flows, as well as the axisymmetric blast flowfield with open end. Effects of the chamber geometry on the nonlinear wave attenuation characteristics is discussed in detail.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.17-29
• /
• 2004

Internal waves are computed using a ghost fluid method on an unstructured grid. Discontinuities in density and dynamic pressure are captured in one cell without smearing or oscillations along a multimaterial interface. A time-accurate incompressible Navier-Stokes/Euler solver is developed based on a three-point backward difference formula for the physical time marching. Artificial compressibility is introduced with respect to pseudotime and an implicit method is used for the pseudotime iteration. To track evolution of an interface, a level set function is coupled with the governing equations. Roe's flux difference splitting method is used to calculate numerical fluxes of the coupled equations. To get higher order accuracy, dependent variables are reconstructed based on gradients which are calculated using Gauss theorem. For each edge crossing an interface, dynamic pressure is assigned for a ghost node to enforce the continuity of total pressure along the interface. Solitary internal waves are computed and the results are compared with other computational and experimental results.

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

A three dimensional inviscid parallel flow solver has been developed for the simulation of rotor blades in forward flight. The computational domain is divided into stationary and rotating zones for the more efficient mesh adaptation. The conservative mesh treatment algorithm is used for the convection of flow variables and fluxes across the sliding boundary. A deforming mesh algorithm using modified spring analogy is used for the blade motion. In the present paper, detail descriptions of numerical analysis for forward flight are introduced. Some results are presented for a two bladed AH-1G rotor and compared with experimental data.

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

A comprehensive study has been made for the investigation of the convergence characteristics of the LU scheme for the Euler equations using von Neumann stability analysis. The stability results indicate that the convergence rate is governed by a specific parameter combination. Based on this insight it is shown that the LU scheme will not suffer convergence deterioration at any grid aspect ration if the local time step is defined using appropriate parameter combination. The numerical results demonstrate that this time step definition gives uniform convergence for grid aspect ratios from one to $1\times10^4$.

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

In this study, a multi-block structured grid deformation code based on a hybrid of a transfinite interpolation algorithm and spring analogy was developed. The configuration was modeled by a Bezier surface. A combination of the spring analogy for block vertices and the transfinite interpolation for interior grid points helps to increase the robustness and makes it suitable for distributed computing. An elliptic smoothing operator was applied to the block faces with sub-faces in order to maintain the grid smoothness and skewness. The capability of this code was demonstrated on a range of simple and complex configurations including an airfoil and a wing-body configuration.

• 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.