• Journal of computational fluids engineering
• /
• v.14 no.3
• /
• pp.33-48
• /
• 2009

In the present study, a fast grid deformation technique has been incorporated into the unsteady compressible and incompressible viscous flow solvers on unstructured hybrid meshes. An algebraic method based on the basis decomposition of normal edge vector was used for the deformation of viscous elements, and a ball-vertex spring analogy was adopted for inviscid elements among several spring analogy methods due to its robustness. The present method was validated by comparing the results obtained from the grid deformation and the rigid motion of entire grids. Fish swimming motion of an NACA0012 airfoil and flapping wing motion of a generic fighter were also simulated to demonstrate the robustness of the present grid deformation technique.

• Journal of computational fluids engineering
• /
• v.18 no.3
• /
• pp.79-86
• /
• 2013

Preconditioned compressible Navier-Stokes equations are solved for almost incompressible flows. Unstructured meshes are utilized for spatial discretization of complex flow domain. Effectiveness of the current preconditioning algorithm, with respect to various Reynolds numbers and Mach numbers, is demonstrated by the solution of canonical problems for incompressible flows, e.g. driven cavity flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.311-317
• /
• 2011

The present paper deals with the continuous works of extending the multi-dimensional limiting process (MLP) for compressible flows, which has been quite successful in finite volume methods, into discontinuous Galerkin (DG) methods. From the series of the previous, it was observed that the MLP shows several superior characteristics, such as an efficient controlling of multi-dimensional oscillations and accurate capturing of both discontinuous and continuous flow features. Mathematically, fundamental mechanism of oscillation-control in multiple dimensions has been established by satisfaction of the maximum principle. The MLP limiting strategy is extended into DG framework, which takes advantage of higher-order reconstruction within compact stencil, to capture detailed flow structures very accurately. At the present, it is observed that the proposed approach yields outstanding performances in resolving non-compressive as well as compressive flaw features. In the presentation, further numerical analyses and results are going to be presented to validate that the newly developed DG-MLP methods provide quite desirable performances in controlling numerical oscillations as well as capturing key flow features.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.13-22
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second-order scheme.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1752-1756
• /
• 2008

A simple, accurate and efficient mesh generation technique, the cut-cell method, is able to represent an arbitrarily complex geometry. Both structured and unstructured grid meshes are used in this method. First, the numerical domain is constructed with regular Cartesian grids as a background grid and then the solid boundaries or bodies are cut out of the background Cartesian grids. As a result, some boundary cells can be contained two numerical conditions such as the flow and solid conditions, where the special treatment is needed to simulate such physical characteristics. The HLLC approximate Riemann solver, a Godunov-type finite volume method, is employed to discretize the advection terms in the governing equations. Also, the TVD-WAF method is applied on the Cartesian cut-cell grids to stabilize numerical results. Present method is validated for the rectangular dam break problems. Initially, a conventional grid is constructed with the Cartesian regular mesh only and then applied to the dam-break flow simulation. As a comparative simulation, a cut-cell grids are applied to represent the flow domain rotated with arbitrary angles. Numerical results from this study are compared with the results from the case of the Cartesian regular mesh only. A good agreement is achieved with other numerical results presented in the literature.

• Journal of Energy Engineering
• /
• v.17 no.4
• /
• pp.218-226
• /
• 2008

A component-scale two-phase analysis code has been developed for a realistic simulation of two-phase flow transients in a light water nuclear reactor component. In the code, a two-fluid three-field model is adopted and the governing equations are solved on an unstructured mesh. For the numerical solution scheme, the semi-implicit method used in the RELAP5 code was selected, which has been proved to be very stable and accurate for most of practical applications. However, some modifications were needed for its application to an unstructured non-staggered grid. This paper presents the modified semi-implicit numerical method for unstructured grid and the preliminary results of the calculations.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.93-98
• /
• 2001

An adaptive Cartesian grid method having the best elements of structured, unstructured, and Cartesian grids is developed to solve the steady two-dimensional Euler equations. The solver is based on a cell-centered finite-volume method with Roe's flux-difference splitting and implicit point Gauss-seidel time integration method. Calculations of several compressible flows are carried out to show the efficiency of the developed computer code. The results were generally in good agreements with existing data in the literature and the developed code has the good ability to capture important feature of the flows.

• Journal of computational fluids engineering
• /
• v.15 no.1
• /
• pp.81-87
• /
• 2010

In this study, we conducted resistance and propulsion performance test of ship composed of the Resistance Test, Propeller Open Water Test and Self Propulsion Test using the CFD(Computational Fluid Dynamics). We used commercial RANS(Reynolds Averaged Navier Stokes equation) solver, as a calculating tool. The unstructured grids were used in a bow and stern of ship, having complex shape, for a convenience of generating grids, and the structured grids were adopted in a central hull and rest of hull having a relatively simple shape which is called hybrid grid method. In addition, The sliding mesh method was adopted to rotate a propeller directly in the Propeller Open Water and Self Propulsion Test. The Resistance Test and Self Propulsion Test were calculated using Volume of Fluid (VOF) model and considering a free surface. And all The three cases were applied realizable k-epsilon model as the turbulence model. The results of calculations were verified for the suitability of calculations by comparing MOERI's EFD results.

• Journal of computational fluids engineering
• /
• v.10 no.3 s.30
• /
• pp.70-76
• /
• 2005

A two-layer incompressible time-accurate Euler solver is applied to analyze flow fields around a submerged body moving at a critical speed near a pycnocline. Discontinuities in the dependent variables across the material interface are captured without any dissipation or oscillation using the ghost fluid method on an unstructured grid. It is shown that the material interlace has significant effects on forces acting on a submerged body moving near a pycnocline regardless of the small difference in densities of two layers. Contrary to the shallow water waves, a submerged body can reach a critical speed at very low Froude number due to the small difference in the densities of the two layers.

• Journal of Mechanical Science and Technology
• /
• v.15 no.5
• /
• pp.665-670
• /
• 2001

The interaction of a travelling shock with the shear layer of a flat plate is studied computationally. The Euler and Navier-Stokes equations are solved numerically on quadrilateral unstructured adaptive grids. The flat plate is installed horizontally on the central axis of a shock tube. The shear layer is first created by two shock waves at different speeds splitted by a flat plate. A series of small vortices is developed as a consequence in the shear layer. The shock wave reflected at the end wall impinges the shear layer. The complicated shock dynamics in the evolution to the pseudo-steady state is represented with the morphological transformation of a planar shock into an oblique shock.