• Korean Journal of Computational Design and Engineering
• /
• v.1 no.1
• /
• pp.56-64
• /
• 1996

A quadrilateral-triangular mixed grid method for the solution of incompressible viscous flow is presented. The solution domain near the body surface is meshed using elliptic grid geneator to acculately simulate the viscous flow. On the other hand, we used unstructured triangular grid system generated by advancing front technique of a simple automatic grid generation algorithm in the rest of the computational domain. The present method thus is capable of not only handling complex geometries but providing accurate solutions near body surface. The numerical technique adopted here is PISO type finite element method which was developed by the present author. Investigations have been made of two-dimensional unsteady flow of Re=550 past a circular cylinder. In the case of use of the unstructured grid only, there exists a considerable amount of difference with the existing results in drag coefficient and vorticity at the cylinder surface; this may be because of the lack of the grid clustering to the surface that is a inevitable requirement to resolve the viscous flow. However, numerical results on the mixed grid show good agreements with the earlier computations and experimental data.

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

A conservative pressure-based finite-volume numerical method has been developed for computing flow and heat transfer by using an unstructured grid system. The method admits arbitrary convex polyhedra. Care is taken in the discretization and solution procedures to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are found by a novel second-order accurate spatial discretization. Momentum interpolation is used to prevent pressure checkerboarding and the SIMPLE algorithm is used for pressure-velocity coupling. The resulting set of coupled nonlinear algebraic equations is solved by employing a segregated approach, leading to a decoupled set of linear algebraic equations fer each dependent variable, with a sparse diagonally dominant coefficient matrix. These equations are solved by an iterative preconditioned conjugate gradient solver which retains the sparsity of the coefficient matrix, thus achieving a very efficient use of computer resources.

• ETRI Journal
• /
• v.42 no.4
• /
• pp.608-618
• /
• 2020

We present a novel GPU-based ray-casting algorithm for volume rendering of unstructured grid data. Our volume rendering system uses a ray-casting method that guarantees accurate rendering results. We also employ the per-pixel intersection list concept in the Bunyk algorithm to guarantee an accurate result for non-convex meshes. For efficient memory access for the lists on the GPU, we represent the intersection lists for all faces as an array with our novel construction algorithm. With the intersection lists, we perform ray-casting on a GPU, and a GPU thread handles each ray. To increase ray-coherency in a thread block and improve memory access efficiency, we extend a prior image-tile-based work distribution method to fit modern GPU architectures. We also show that a prior approach using a per-thread local buffer to reduce redundant computation is not appropriate for modern GPU architectures. Instead, we take an on-demand calculation strategy that achieves better performance even though it allows duplicate computations. We applied our method to three unstructured grid datasets with different characteristics. With a GPU, our method achieved up to 36.5 times higher performance for the ray-casting process and 19.7 times higher performance for the whole volume rendering process compared with the Bunyk algorithm using a CPU core. Also, our approach showed up to 8.2 times higher performance than a GPU-based cell projection method while generating more accurate rendering results. These results demonstrate the efficiency and accuracy of our method.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2218-2229
• /
• 2006

A conservative finite-volume numerical method for unstructured grids with the cell-centered method has been developed for computing flow and heat transfer by combining the attractive features of the existing pressure-based procedures with the advances made in unstructured grid techniques. This method uses an integral form of governing equations for arbitrary convex polyhedra. Care is taken in the discretization and solution procedure to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. For both convective and diffusive fluxes the forms superior to both accuracy and stability are particularly adopted and formulated through a systematic study on the existing approximation ones. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are computed by using a linear reconstruction based on the divergence theorem. Momentum interpolation is used to prevent the pressure checkerboarding and a segregated solution strategy is adopted to minimize the storage requirements with the pressure-velocity coupling by the SIMPLE algorithm. An algebraic solver using iterative preconditioned conjugate gradient method is used for the solution of linearized equations. The flow analysis code (PowerCFD) developed by the present method is evaluated for its application to several 2-D structured-mesh benchmark problems using a variety of unstructured quadrilateral and triangular meshes. The present flow analysis code by using unstructured grids with the cell-centered method clearly demonstrate the same accuracy and robustness as that for a typical structured mesh.

• International Journal of High-Rise Buildings
• /
• v.8 no.4
• /
• pp.255-264
• /
• 2019

This paper describes large eddy simulation of wind pressures on a square cylinder in a uniform flow and a high-rise building immersed in an atmospheric turbulent boundary layer. For the atmospheric boundary layer case, the inflow turbulence is generated by a numerical wind tunnel. In the numerical simulation, particular attention is devoted to the performance of an auto hexahedral non-structural mesh. Both simulations are performed for three grid systems: an auto hexahedral non-structured grid, a structured Cartesian grid and a non-structured triangular prism grid, and for three grid numbers. The present study shows that the auto hexahedral unstructured mesh achieves the best simulation results for wind pressures on the square cylinder and the high-rise building. When the grid number is sufficiently large, the differences among the results obtained from the three investigated grid systems are not significant. However, the advantage of the auto hexahedral unstructured mesh becomes clear when the grid number decreases, because it enables a balanced distribution of orthogonal grids. The results described in this paper demonstrate that the auto hexahedral non-structured mesh has good potential applicability to simulation of urban flows.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.2
• /
• pp.90-97
• /
• 2004

Flow mechanism of contractive and dilative motion is numerically investigated to obtain a propulsive force in highly viscous fluid. An axisymmetric code is developed with unstructured grid system based on cell-centered scheme. It is validated by comparing with the results of Stokes approximation for the problem of uniform flow past a sphere in low Reynolds number(R$_{n}$=1). The validated code is applied to the simulation of contractive and dilative periodic motion of body whose results are quantitatively compared with the two dimensional case. In order to investigate the grid dependency, two different grids are applied to the present computations. The present study provides key information for the development of an axisymmetric Micro-hydro-robot.t.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5B
• /
• pp.511-517
• /
• 2006

Flood damage is one of the most important and influential natural disaster which has an effect on human beings. Local concentrated heavy rainfall in urban area yields flood damage increase due to insufficient capacity of drainage system. When the excessive flood occurs in urban area, it yields huge property losses of public facilities involving roadway inundation to paralyze industrial and transportation system of the city. To prevent such flood damages in urban area, it is necessary to develop adequate inundation analysis model which can consider complicated geometry of urban area and artificial drainage system simultaneously. In this study, an urban flood analysis model adopting the unstructured computational grid was developed to simulate the urban flood characteristics such as inundation area, depth and integrated with subsurface drainage network systems. By the result, we can make use of these presented method to find a flood hazard area and to make a flodd evacuation map. The model can also establish flood-mitigation measures as a part of the decision support system for flood control authority.

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

The present paper deals with the accurate and robust limiting procedure for the multi-dimensional flow analysis on unstructured meshes. The multi-dimensional limiting process (MLP) which was successfully proposed on structured grid system is extended to unstructured meshes. Based on MUSCL-type framework on unstructured meshes, the new slope limiter is devised to satisfy the MLP condition, which is quite effective to regulate the unwanted oscillations, especially on multiple dimensions. Considering the neighborhood based on the vertex of the cell, as well as the edge, this limiting strategy captures the multi-dimensional flow features very accurately with the proper stencils. From the various numerical results, these desirable characteristics of the proposed limiting strategy are clearly shown.

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

The present paper deals with the accurate and robust limiting procedure for the multi-dimensional flow analysis on unstructured meshes. The multi-dimensional limiting process (MLP) which was successfully proposed on structured grid system is extended to unstructured meshes. Based on MUSCL-type framework on unstructured meshes, the new slope limiter is devised to satisfy the MLP condition, which is quite effective to regulate the unwanted oscillations, especially on multiple dimensions. Considering the neighborhood based on the vertex of the cell, as well as the edge, this limiting strategy captures the multi-dimensional flow features very accurately with the proper stencils. From the various numerical results, these desirable characteristics of the proposed limiting strategy are clearly shown.

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