• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.201-210
• /
• 2009

A volume capturing method using unstructured grid system for numerical analysis of multiphase flows is introduced in the present paper. This method uses an interface capturing method (CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The novelty of CICSAM lies in the adaptive combination of high resolution discretization scheme which ensures the preservation of the sharpness and shape of the interface while retaining boundedness of the field, and no explicit interface reconstruction which is perceived to be difficult to implement on unstructured grid system. Several typical test cases for multiphase flows are presented, which are simulated by an in-house solution code(PowerCFD). This code employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with CICSAM. It is found that the present method simulates efficiently and accurately complex free surface flows such as multiphase flows.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.36-41
• /
• 1999

Shock wave propagating in the particle suspension has important applications. Examples are shock waves occurring in the solid rocket plume and detonation of dusty particles by shock waves. Experimental and numerical investigations on this subject have drawn much attention. More recently, Sivier et al. numerically simulated the experiment of Sommerfeld using the unstructured adaptive grid. They used the Eulerian-Eulerian approach based on the continuum assumption for both gas and particles. In the present paper, a new numerical method using the Lagrangian particle tracing technique and unstructured particle-adaptive grid for the polydisperse system is presented. It is explained why the existing numerical calculation has showed discrepancy with the experimental results by Sommerfeld.

• Environmental Sciences Bulletin of The Korean Environmental Sciences Society
• /
• v.10 no.S_1
• /
• pp.23-28
• /
• 2001

A new numerical weather prediction and dispersion model, the Operational Multi-scale Environment model with Grid Adaptivity(OMEGA) including an embedded Atmospheric Dispersion Model(ADM), is introduced as a next generation atmospheric simulation system for real-time hazard predictions, such as severe weather or the transport of hazardous release. OMEGA is based on an unstructured grid that can facilitate a continuously varying horizontal grid resolution ranging from 100 km down to 1 km and a vertical resolution from 20 -30 meters in the boundary layer to 1 km in the free atmosphere. OMEGA is also naturally scale spanning and time. In particular, the unstructured grid cells in the horizontal dimension can increase the local resolution to better capture the topography or important physical features of the atmospheric circulation and cloud dynamics. This means the OMEGA can readily adapt its grid to a stationary surface, terrain features, or dynamic features in an evolving weather pattern. While adaptive numerical techniques have yet to be extensively applied in atmospheric models, the OMEGA model is the first to exploit the adaptive nature of an unstructured gridding technique for atmospheric simulation and real-time hazard prediction. The purpose of this paper is to provide a detailed description of the OMEGA model, the OMEGA system, and a detailed comparison of OMEGA forecast results with observed data.

• Journal of computational fluids engineering
• /
• v.7 no.4
• /
• pp.35-41
• /
• 2002

The projectile afterbodies for zero-lift drag reduction has been analyzed using the Navier-Stokes equations with the κ-εturbidence model. The numerical method of a second order upwind scheme has been used on an unstructured adaptive grid system. Base drag reduction methods that have been found effective on axisymmetric bodies are boattailing, base bleed, base combustion, locked vortex afterbodies and multistep afterbodies. In this paper, turbulence flow and pressure charateristics have been studied for geometries of multistep afterbodies. The important geometrical and flow parameters relevant to the design of such afterbodies have been identified by step number, length and height. The flow over multistep aftoerbodies or base have many kinds of compressible flow characteristics including expansion waves at the trailing edge, recompression waves, separation and recirculating flow in the base region, shear flow and wake flow. The numerical results have been compared and analyzed with the experimental data. The flow characteristics have been clearly shown.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.8
• /
• pp.1-7
• /
• 2005

This paper proposes a local improvement method that combines extended topological clean up and optimization-based smoothing for homogenizing triangular grid system. First extended topological clean up procedures are applied to improve the connectivities of grid elements. Then, local optimization-based smoothing is performed for maximizing the distortion metric that measures grid quality. Using the local improvement strategy, we implement the grid homogenizations for two triangular grid examples. It is shown that the suggested algorithm improves the quality of the triangular grids to a great degree in an efficient manner and also can be easily applied to the remeshing algorithm in adaptive mesh refinement technique.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.39-47
• /
• 2006

In this paper, the authors developed a multidimensional interpolation method inside a finite volume cell in the computation of high-order accurate numerical flux such as the fifth order WEND (weighted essentially non-oscillatory) scheme. This numerical method starts from a simple Taylor series expansion in a proper spatial order of accuracy, and the WEND filter is used for the reconstruction of sharp nonlinear waves like shocks in the compressible flow. Two kinds of interpolations are developed: one is for the cell-averaged values of conservative variables divided in one mother cell (Type 1), and the other is for the vertex values in the individual cells (Type 2). The result of the present study can be directly used to the cell refinement as well as the convective flux between finer and coarser cells in the Cartesian adaptive grid system (Type 1) and to the post-processing as well as the viscous flux in the Navier-Stokes equations on any types of structured and unstructured grids (Type 2).

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.1
• /
• pp.182-189
• /
• 1993

To improve the resolution of complex flow field features, grid adaptation scheme of Anderson has been revised, which was based on the Poisson grid generator of Thompson. Anderson's original scheme adapts the grid to solution automatically, but if flow field is more or less complex, then the adaptivity is weak. So the technique of using threshold which is used in unstructured grid system is adopted. The regions of large variation in the solution are marked by marking function which has the property of total variation of the solution, and these regions have same values of weight but other regions are neglected. This updated method captures shocks clearly and sharpy. Four examples are demonstrated, (1) Hypersonic flow past a blunt body, (2) High speed inlet analysis, (3) Supersonic flow of M=1.4 over a 4% biconvex airfoil in a channel, (4) Hypersonic shock-on-shock interaction at M=8.03.