• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.48 no.3
• /
• pp.200-208
• /
• 1999

In this paper, the corona discharge is analyzed by Finite Element Method(FEM) combined with Flux-corrected Transport(FCT) algorithm. In the previous papers, Finite Difference Method(FDM) combined with FCT was used. Usually in the FDM, the regionof interest is discretized with structured grids. But to refine local regions with same resolution, much more grids are required for the structured grids than for unstructured grids than for unstructured grids. Therefore, we propose the FEM-FCT method to simulate the corona discharge. The proposed method has good flexibility in model shape and can reduce the computational cost by the local refinement where the physical quantities have steep gradients. Using the proposed method, we study the streamer growth of parallel plate electrodes which is initiated by the low and high perturbation density. We find that the varying the initial density of perturbation has very little effect on the streamer propagation. And the corona discharge of the rod-to-plane electrode is simulated. On the surface of the rod electrode, the high concentration of the electric field gives rise to many number of streamer seeds. The strong axial streamer propagate to the plane electrode. The weaker non-axial streamer repel each other and stop growing more. The results are very similar to those of the papers which used the FDM-FCT method on structured grids. Thus we can conclude that the proposed FEM-FCT method is more efficient than the conventional FDM-FCT method by virtue of the reduction in computational grids number.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.153-158
• /
• 2005

Numerical simulations of two-dimensional steady incompressible lid-driven flow in a square cavity are presented to verify the validity of a new solution code(PowerCFD) with unstructured grids. The code uses the non-staggered(collocated) grid approach which is very popular for incompressible flow analysis because of its numerical efficiency on the curvilinear or unstructured grids. Solutions are obtained for configurations with a Reynolds number as high as 10,000 with both rectangular and hybrid types of unstructured grid mesh. Interesting features of the flow are presented in detail and comparisons are made with benchmark solutions found in the literature. It is found that the code is capable of producing accurately the nature of the lid-driven cavity flow at high Reynolds numbers.

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

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

Three-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.910-915
• /
• 2003

A Navier-Stokes equation solver for incompressible viscous flows with free surface is developed and tested. This is based upon a fractional time step method and a non-staggered finite volume formulation for unstructured meshes. For time advancement scheme, Adams -Bashforth method for convective term and Crank-Nicolson method for diffusive term are applied. The interface between two fluids with different fluid properties is tracked with Piecewise Linear Interface Calculation(PLIC) Volume-of-Fluid(VOF) methods. Computational results are presented for some test problems: the broken dam, the sloshing in a rectangular tank, the filling of a cylindrical tank.

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

A field reconstruction scheme for a cell centered finite volume method on unstructured meshes is developed. Regardless of mesh quality, this method is exact within a machine accuracy if the solution is linear, which means it has full second order accuracy. It does not have any limitation on cell shape except convexity of the cells and recovers standard discretization stencils at structured orthogonal grids. Accuracy comparisons with other popular reconstruction schemes are performed on a simple example.

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

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.44-53
• /
• 2007

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed for a 3D component of a nuclear system code and a component-scale analysis tool. A two-fluid three-field model is used for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. An unstructured grid is adopted for realistic simulations of the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been applied to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the modified numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing very well.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09b
• /
• pp.1213-1214
• /
• 2005

• Journal of Korea Water Resources Association
• /
• v.31 no.3
• /
• pp.279-290
• /
• 1998

The height and speed of the shock wave are critical data in flood-control operations or in the design of channel walls and bridges along rivers with high flow velocities. Therefore, a numerical model is needed for simulating flow discontinuity over a wide range of conditions. In this study, a governing equation. As a Riemann solver Roe(1981)'s one is used. The model employs the modified MUSCL for handling the unstructured grids in this research. this model that adopts the explicit tradditional twl dimmensional dam break problems, two hydraulic dam break model is simulations, and a steady state simulation in a curved channel. Conclusions of this research are as follows : 1) the finite volume method can be combined with the Godonov-type method that is useful for modeling shocks. Hence, the finite volume method is suitable for modeling shocks. 2) The finite volume model combined with the modified MUSCL is successful in modeling shock. Therefore, modified MUSCL is proved to be valid.