• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.69-78
• /
• 2008

The high-order numerical method based on the adaptive mesh refinement(AMR) on the quadrilateral unstructured grids has been developed in this paper. This adaptive-grid method, originally developed with MUSCL-TVD scheme, is now extended to the WENO (weighted essentially no-oscillatory) scheme with the Runge-Kutta time integration of fifth order in spatial and temporal accuracy. The multidimensional interpolation was studied in the preliminary research, which allows us to maintain the same order of accuracy for the computation of numerical flux between two adjacent cells of different levels. Some standard benchmark tests are done to validate this method for checking the overall capacity and efficiency of the present adaptive-grid technique.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.91-94
• /
• 2007

A GUI based pre/post processor program, which is based on the MFC and OpenGL library in the Windows O/S, hee been developed for NUFLEX Using this program, users are able to generate and modify structured or unstructured grid geometries, set all the parameters for the solver, and observe the results of the simulation in graphic view by vector or scalar plots. The mesh geometry data can be imported from or exported to other programs by supporting functions for reading from and writing to CGNS data format files.

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

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.98-105
• /
• 1999

The Pressure-based Unstructured-grid Navier-Stokes Solver PUNS-2/3D for incompressible steady and unsteady viscous flows has been developed. It is based on nonstaggered cell-centered finite volume method. Second-order upwind scheme with least-square reconstruction is used for convective fluxes. The SIMPLE method is implemented to couple the pressure and velocity fields. And the time derivatives in the momentum equations are discretised using a second-order Euler backward-differencing scheme. The discretised linear equations are solved by the preconditioned Biconjugate Gradient Stabilized method(Bi-CGSTAB). The developed solver is applied to validation problems using hybrid meshes.

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

• Transactions of the Korean Society of Automotive Engineers
• /
• v.6 no.2
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.333-334
• /
• 2006

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.86-95
• /
• 2008

For analyses of multi-phase flows in a water-cooled nuclear power plant, a three-dimensional SIMPLE-algorithm based hydrodynamic solver CUPID-S has been developed. As governing equations, it adopts a two-fluid three-field model for the two-phase flows. The three fields represent a continuous liquid, a dispersed droplets, and a vapour field. The governing equations are discretized by a finite volume method on an unstructured grid to handle the geometrical complexity of the nuclear reactors. The phasic momentum equations are coupled and solved with a sparse block Gauss-Seidel matrix solver to increase a numerical stability. The pressure correction equation derived by summing the phasic volume fraction equations is applied on the unstructured mesh in the context of a cell-centered co-located scheme. This paper presents the numerical method and the preliminary results of the calculations.