• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.894-899
• /
• 2000

An automatic mesh generation system with unstructured quadrilateral elements on trimmed NURBS surfaces has been developed.. In this paper, NURBS surface geometries in the IGES format have been used to represent model shape. NURBS surface is represented as parametric surface. So each surface could be mapped to a 2D parametric plane through the parametric domain. And then meshes with quadrilateral elements are constructed in this plane. Finally, the constructed meshes are mapped back to the original 3D surface through the parametric domain. In this paper, projection plane, quasi-expanded plane and parametric Plane are used as 2D mesh generation plane. For mapping 3D surface to parametric domain, Newton-Rhapson Method is employed. For unstructured mesh generation with quadrilateral elements on 2D plane, a domain decomposition algorithm using loop operators has been employed. Sample meshes are represented to demonstrate the effectiveness of the proposed algorithm.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.1057-1064
• /
• 2005

A conservative finite-volume numerical method using unstructured meshes, which is developed by the authors, is evaluated for its application to several 2-D benchmark problems using a variety of quadrilateral, triangular and hybrid meshes. The present pressure-based numerical method for unstructured mesh clearly demonstrates the same accuracy and robustness as that fur typical structured mesh.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.113-117
• /
• 1995

Viscous analysis on incompressible flows is performed using unstructured triangular meshes. A two-dimensional and axisymmetric incompressible Navier-Stokes equations are solved in time-marching form by artificial compressibility method. The governing equations are discretized by a cell-centered based finite-volume method. and a centered scheme is used for inviscid and viscous fluxes with fourth order artificial dissipation. An explicit multi-stage Runge-Kutta method is used for the time integration with local time stepping and implicit residual smoothing. Convergence properties are examined and solution accuracies are also validated with benchmark solution and experiment.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.57-60
• /
• 2006

This paper presents an assessment of vertex-centered and cell-centered finite-volume methods on unstructured meshes. The results indicate that the vertex-centered method is more reliable than the cell-centered method.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.8
• /
• pp.1-11
• /
• 2004

A comprehensive study has been made for the investigation of the convergence and stability characteristics of the LU scheme for solving the Navier-Stokes equations on unstructured meshes. For this purpose the characteristics of the LU scheme was initially studied for a scalar model equation. Then the analysis was extended to the Navier-Stokes equations. It was shown that the LU scheme has an inherent stiffness in the streamwise direction. This stiffness increases when the grid aspect ratio becomes high and the cell Reynolds number becomes small. It was also shown that the stiffness related to the grid aspect ratio can be effectively eliminated by performing proper subiteration. The results were validated for a flat-plate turbulent flow.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.37-45
• /
• 2002

A two-dimensional unsteady, inviscid flow solver has been developed for the simulation of airfoil-vortex interactions on unstructured dynamically adapted meshes. The Euler solver is based on a second-order accurate implicit time integration using a point Gauss-Seidel relaxation scheme and a dual time-step subiteration. A vertex-centered, finite-volume discretization is used in conjunction with the Roe's flux-difference splitting. An unsteady solution-adaptive dynamic mesh scheme is used by adding and deleting mesh points to take account of both spatial and temporal variations of the flow field. The effect of vortex interaction on the aeroelastic displacement of an airfoil attached to the idealized two degree-of-freedom spring system is investigated.

• Journal of computational fluids engineering
• /
• v.17 no.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.

• Journal of computational fluids engineering
• /
• v.19 no.1
• /
• pp.27-33
• /
• 2014

A new and efficient implicit scheme is proposed to obtain a steady-state solution in time integration and the comparison of characteristics with the approximation ways for the implicit method to solve the incompressible Navier-Stokes equations is provided. The conservative, finite-volume cell-vertex upwind scheme and artificial compressibility method using dual time stepping for time accuracy is applied in this paper. The numerical results obtained indicate that the direct application of Jacobian matrix to the Lower and upper sweeps of implicit LU-SGS leads to better performance as well as convergence regardless of CFL number and true time step than explicit scheme and approximation of Jacobian matrix. The flow simulation around box in uniform flow with unstructured meshes is demonstrated to check the validity of the current formulation.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.9
• /
• pp.1-11
• /
• 2004

A comprehensive study has been made for the investigation of the convergence and stability characteristics of the LU scheme for solving the Euler equations on unstructured meshes. The von Neumann stability analysis technique was initially applied to a scalar model equation, and then the analysis was extended to the Euler equations. The results indicated that the convergence rate is governed by a specific combination of flow parameters. Based on this insight, it was shown that the LU scheme does not suffer any convergence deterioration at all grid aspect ratios, as long as the local time step is defined using an appropriate parameter combination.

• The Bulletin of The Korean Astronomical Society
• /
• v.38 no.2
• /
• pp.101.2-101.2
• /
• 2013

We present a new computational fluid dynamic (CFD) simulation code. The code employs the moving and polyhedral unstructured mesh scheme, which is known as a superior approach to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for surfaces of every Voronoi cell to update the hydrodynamic states as well as to move former generated meshes. For the second-order accuracy, the MUSCL-Hancock scheme is implemented. To increase efficiency for generating Voronoi tessellation we also develop the incremental expanding method, by which the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We will discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.