• Journal of Mechanical Science and Technology
• /
• v.19 no.4
• /
• pp.1006-1017
• /
• 2005

In this work, we developed an accurate and efficient radiative finite volume method applicable for the complex 2D planar and 3D geometries using an unstructured-grid finite volume method. The present numerical model has fully been validated by several benchmark cases including the radiative heat transfer in quadrilateral enclosure with isothermal medium, tetrahedral enclosure, a three-dimensional idealized furnace, as well as convection-coupled radiative heat transfer in a square enclosure. The numerical results for all cases are well agreed with the previous results. Special emphasis is given to the parallelization of the unstructured-grid radiative FVM using the domain decomposition approach. Numerical results indicate that the present parallel unstruc­tured-grid FVM has the good performance in terms of accuracy, geometric flexibility, and computational efficiency.

• The KSFM Journal of Fluid Machinery
• /
• v.8 no.3 s.30
• /
• pp.7-15
• /
• 2005

A Navier-Stokes code has been developed to simulate the flow through a cross flow fan. It is based on an unstructured finite volume method and uses moving grid technique to model the rotation of the fan. A low Reynolds number turbulence model is used to calculate eddy viscosity. The basic algorithm is SIMPLE. Numerical simulations over a wide range of flow rate aye carried out to validate the code. Comparison of all numerical solutions with experimental data confirms the validity of the present code. Present numerical solutions show a noticeable improvement over a previous numerical method which is based on a model of body force to simulate the rotation of the impeller.

• The KSFM Journal of Fluid Machinery
• /
• v.9 no.4 s.37
• /
• pp.27-35
• /
• 2006

A Navier-Stokes code has been developed to simulate the flow through a cross flow fan. It is based on an unstructured finite volume method and uses moving grid technique to model the rotation of the fan. A low Reynolds number turbulence model is used to calculate eddy viscosity. The basic algorithm is SIMPLE. Numerical simulations over a wide range of flow rate are carried out to validate the code. Comparison of all numerical solutions with experimental data confirms the validity of the present code. Present numerical solutions show a noticeable improvement over a previous numerical method which is based on a model of body force to simulate the rotation of the impeller.

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

• Journal of applied mathematics & informatics
• /
• v.9 no.2
• /
• pp.561-577
• /
• 2002

In this paper, a two-dimensional finite volume unstructured mesh method (FVUM) based on a triangular background interpolation mesh is developed for analysing the evolution of the saltwater intrusion into single and multiple coastal aquifer systems. The model formulation consists of a ground-water flow equation and a salt transport equation. These coupled and non-linear partial differential equations are transformed by FVUM into a system of differential/algebraic equations, which is solved using backward differentiation formulas of order one through five. Simulation results are compared with previously published solutions where good agreement is observed.

• Journal of ILASS-Korea
• /
• v.5 no.1
• /
• pp.13-22
• /
• 2000

The present study has numerically analyzed the vaporization characteristics of fuel droplets in the high temperature convective flow field. The axisymmetric governing equations for mass, momentum, energy, and species are solved by an iterative and implicite unstructured finite-volume method. The moving boundary due to vaporization is handled by the deformable unstructured grid technique. The pressure-velocity coupling in the density-variable flows is treated by the SIMPLEC algorithm. In terms of the matrix solver, Bi-CGSTAB is employed for the numerically efficient and stable convergence. The n-decane is used as a liquid fuel and the initial droplet temperature is 300K. Computations are performed for the nonevaporating and evaporating droplets with the relative interphase velocity(25m/s). The unsteady vaporization process has been simulated up to the nondimensional time, 25. Numerical results indicate that the mathematical model developed in this study succesfully simulates the main features of the droplet vaporization process in the convective environment.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.4
• /
• pp.179-189
• /
• 1996

An unstructured grid finite-volume method has been applied to predict the linear and nonlinear attenuation characteristics of the expansion chamber silencer system. In order to achieve a grid flexibility and a solution adaptation for geometrically silencer system. In order to achieve a grid flexibility and a solution adaptation for geometrically complex flow regions associated with the actual silencers, the unstructured mesh algorithm in context with the node-centered finite volume method has been employed. The present numerical model has been validated by comparison with the analytical solutions and the experimental data for the acoustic field of the concentric expansion chamber with and without pulsating flows, as well as the axisymmetric blast flowfield with open end. Effects of the chamber geometry on the nonlinear wave attenuation characteristics is discussed in detail.

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

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

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