• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.2
• /
• pp.122-132
• /
• 2013

In this research the development of unstructured grid discretization solution techniques is presented. The purpose is to describe such a conservative discretization scheme applied for experimental validation work. The objective of this paper is to better establish the effects of mesh generation techniques on velocity fields and particle deposition patterns to determine the optimal aerodynamic characteristics. In order to achieve the objective, the mesh surface discretization approaches used the VLA prototype manufacturing tolerance zone of the outer surface. There were 3 schemes for this discretization study implementation. They are solver validation, grid convergence study and surface tolerance study. A solver validation work was implemented for the simple 2D and 3D model to get the optimum solver for the VLA model. A grid convergence study was also conducted with a different growth factor and cell spacing, the amount of mesh can be controlled. With several amount of mesh we can get the converged amount of mesh compared to experimental data. The density around surface model can be calculated by controlling the number of element in every important and sensitive surface area of the model. The solver validation work result provided the optimum solver to employ in the VLA model analysis calculation. The convergence study approach result indicated that the aerodynamic trend characteristic was captured smooth enough compared with the experimental data. During the surface tolerance scheme, it could catch the aerodynamics data of the experiment data. The discretization studies made the validation work more efficient way to achieve the purpose of this paper.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.20 no.4
• /
• pp.419-425
• /
• 2014

The current work investigated the variation of numerical solutions according to the grid number, the distance of the first grid point off the ship surface, turbulence modeling and discretization. The subject vessel is KVLCC. A commercial code, Gridgen V15 and FLUENT were used the generation of the ship hull surface and spatial system and flow computation. The first part of examination, the effect of solutions were accessed depending on the grid number, turbulence modeling and discretization. The second part was focus on the suitable selection of the distance of the first grid point off the ship surface: $Y_P+$. When grid number and discretization were fixed the same value, the friction resistance showed differences within 1 % but the pressure resistance showed big differences 9 % depending on the turbulence modeling. When $Y_P+$ were set 30 and 50 for the same discretization, friction resistance showed almost same results within 1 % according to the turbulence modeling. However, when $Y_P+$ were fixed 100, friction resistance showed more differences of 3 % compared to $Y_P+$ of 30 and 50. Whereas pressure resistance showed big differences of 10 % regardless of turbulence modeling. When turbulence modeling and discretization were set the same value, friction, pressure and total resistance showed almost same result within 0.3 % depending on the grid number. Lastly, When turbulence modeling and discretization were fixed the same value, the friction resistance showed differences within 5~8 % but the pressure resistance showed small differences depending on the $Y_P+$.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.33-38
• /
• 2000

A Numerical study of laminar vortex-shedding past a circular cylinder has been performed widely by many researchers. Many factors, such as numerical technique and domain size, number and shape of grid, affected predicting vortex shedding and Strouhal number. In the present study, the effect of convection scheme, time discretization methods and grid dependence were investigated. The present paper presents the finite volume solution of unsteady flow past circular cylinder at Re=200, 400. The Strouhal number was predicted using UDS, CDS, Hybrid, Power-law, LUDS, QUICK scheme for convection term, implicit and crank-nicolson methods for time discretization. The grid dependence was investigated using H-type mesh and O-type mesh. It also studied that the effect of mesh size of the nearest adjacent grid of circular cylinder. The effect of convection scheme is greater than the effect of time discretization on predicting Strouhal. It has been found that the predicted Strouhal number changed with mesh size and shape.

• KIEAE Journal
• /
• v.16 no.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.

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

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

In order to reduce the excessive numerical dissipation, the new spatial discretization scheme is introduced. The present method in this paper has the formula that has an additional procedure of defining transferred properties at a cell-interface, based on AUSMPW+. The newly defined transferred property could eliminate numerical dissipation effectively in non-flow aligned grid system. In addition, the present method guarantees the monotonic characteristic in capturing a discontinuity. Through a stationary or moving contact discontinuity and a stationary or moving shock discontinuity, a vortex discontinuity and shock wave/ boundary layer interaction, it is verified that the accuracy of the present method is improved.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.96-101
• /
• 1998

Various discretiation methods of Laplacian operator in the Pressure-Poisson equation are investigated for the solution of incompressible Navier-Stokes equations using the non-staggered grid. Laplacian operators previously proposed by other researchers are applied to a Driven-Cavity problem. The computational results are compared with those of Ghia. The results show the characteristics of the discrete Laplacian operators.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.11
• /
• pp.1521-1530
• /
• 2002

This article proposes a pressure based method for predicting flows at all speeds. The compressible SIMPLE algorithm is extended to unstructured grid framework. Convection terms are discretized using second-order scheme with deferred correction approach. Diffusion term discretization is based on structured grid analogy that can be easily adopted to hybrid unstructured grid solver. This method also uses node centered scheme with edge based data structure for memory and computing time efficiency of arbitrary grid types. Both incompressible and compressible benchmark problems are solved using the above methodology. The demonstration of this method is extended to slip flow problem that has low Reynolds number but compressibility effect. It is shown that the proposed method can improve efficiency in memory usage and computing time without losing any accuracy.

• Journal of the Korean earth science society
• /
• v.39 no.4
• /
• pp.317-326
• /
• 2018

Effect of the nonuniform grid on the two-dimensional transport equation was investigated in terms of theoretical analysis and finite difference method (FDM). The nonuniform grid having a typical structure of the numerical weather forecast model was incorporated in the vertical direction, while the uniform grid was used in the zonal direction. The staggered and non-staggered grid were placed in the vertical and zonal direction, respectively. Time stepping was performed with the third-order Runge Kutta scheme. An error analysis of the spatial discretization on the nonuniform grid was carried out, which indicated that the combined effect of the nonuniform grid and advection velocity produced either numerical diffusion or numerical adverse-diffusion. An analytic function is used for the quantitative evaluation of the errors associated with the discretized transport equation. Numerical experiments with the non-uniformity of vertical grid were found to support the analysis.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.19-22
• /
• 1998

Automatic mesh generation procedures applied to industrial now problems lead to complex mesh topologies where usually no special considerations to mesh resolution are taken. In the present study a fast and flexible solution algorithm in combination with generalized higher order discretization schemes is presented and its application to intake port calculation is demonstrated.