• Fire Science and Engineering
• /
• v.6 no.1
• /
• pp.17-22
• /
• 1992

We have derived the general transfer equation for governing the continuity, energy transfer, mass and momentum transfer, and turbulent energy dissipation rate within the fire compartment which has the 800t fire source at the center of the floor. The governing transfer equations have been descretized using the finite volume approach and numerically experimented under the SIMPLE algorithm. In order for the SIMPLE algorithm approach to be physically reliable, the test results are compared with those of Morita's SOR Method using Conjugate Residual Method and found to be close to physical values though the computational convergence time still remains to be upgraded. The treatment of source terms in the system of finite difference equations has been critical in order to converge the governing equations within the appropriate time steps. The criteria of convergence allowance for the whole domain have been checked and the sudden change of the non-linear effects from the source term have been avoided. The criteria has been allowed to be for 5$\times$10$^{-5}$ .

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

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.609-616
• /
• 2005

While the SIMPLE algorithm is most widely used for the simulations of flow phenomena that take place in the industrial equipments or the manufacturing processes, it is less adopted for the simulations of the free surface flow. Though the SIMPLE algorithm is free from the limitation of time step, the free surface behavior imposes the restriction on the time step. As a result, the explicit schemes are faster than the implicit scheme in terms of computation time when the same time step is applied to, since the implicit scheme includes the numerical method to solve the simultaneous equations in its procedure. If the computation time of SIMPLE algorithm can be reduced when it is applied to the unsteady free surface flow problems, the calculation can be carried out in the more stable way and, in the design process, the process variables can be controlled based on the more accurate data base. In this study, a modified SIMPLE algorithm is presented fur the free surface flow. The broken water column problem is adopted for the validation of the modified algorithm (MoSIMPLE) and for comparison to the conventional SIMPLE algorithm.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2000.11a
• /
• pp.665-672
• /
• 2000

The problems of discontinuous layer interface are very important in the algorithm and programming for the analysis of multi-layered consolidation using a numerical analysis, finite difference method(F.D,M.). Better results can be obtained by the process for discontinuous layer interface, since it can help consolidation analysis to model the actual ground Explicit method is simple for analysis algorithm and convenient for use except for applying the operator Crank-Nicolson method represents implicit method, which have different analysis method according to weighting factor. This method uses different algorithm according to dimension. And, this paper uses alternative direction implicit method. The purpose of this paper provides an efficient computer algorithm based on numerical analysis using finite difference method which account for multi-layered soils with confined aquifer to determine the degree of consolidation and excess pore pressures relative to time and positions more realistically.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.675-682
• /
• 2010

The flood impact pressure acting on a vertical wall resulting from a dam-breaking problem is simulated using a navier-Stokes(N-S) solver. The N-S solver uses Eulerian Finite Volume Method(FVM) along with Volume Of Fluid(VOF) method for 2-D incompressible free surface flows. A Split Lagrangian Advection(SLA) scheme for VOF method is implemented in this paper. The SLA scheme is developed based on an algorithm of Piecewise Linear Interface Calculation(PLIC). The coupling between the continuity and momentum equations is affected by using a well-known Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. Several two-dimensional numerical simulations of the dam-breaking problem are presented to validate the accuracy and demonstrate the capability of the present algorithm. The significance of the time step and grid resolution are also discussed. The computational results are compared with experimental data and with computations by other numerical methods. The results showed a favorable agreement of water impact pressure as well as the global fluid motion.

• Nuclear Engineering and Technology
• /
• v.32 no.2
• /
• pp.169-179
• /
• 2000

A detailed numerical analysis of the evolution of thermal stratification in a curved piping system in a nuclear power plant is performed. A finite volume based thermal-hydraulic computer code has been developed employing a body-fitted, non-orthogonal curvilinear coordinate for this purpose. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of momentum interpolation method. The SIMPLE algorithm is employed for the pressure and velocity coupling, and the convection terms are approximated by a higher-order bounded scheme. The thermal-hydraulic computer code developed in the present study has been applied to the analysis of thermal stratification in a curved duct and some of the predicted results are compared with the available experimental data. It is shown that the predicted results agree fairly well with the experimental measurements and the transient formation of thermal stratification in a curved duct is also well predicted.

• Fire Science and Engineering
• /
• v.19 no.4 s.60
• /
• pp.18-25
• /
• 2005

In this study, numerical analysis of two-dimensional unsteady natural convection of air in a square enclosure heated from below, was performed as a basic research of fire science. SIMPLE algorithm was used to the pressure term of momentum equations in the numerical analysis. The numerical analysis were studied for the two model cases and two heat conditions, respectively, which are different with insulation of enclosures and position of heat applied. Also, the ceiling temperatures of enclosure were measured to compare the accuracy of numerical analysis, and it is found that the temperature predicted by numerical analysis were agreed well with the measurements. Streamline and isotherm of the each model case were acquired for each time step.

• Structural Engineering and Mechanics
• /
• v.4 no.2
• /
• pp.177-190
• /
• 1996

A method to numerically evaluate the behaviour of single span beams, prestressed with external tendons and symmetrically loaded is presented. This algorithm, based on the Finite Difference Method, includes second order effects and large displacements in an attempt to more fully understand the behaviour of the beam up to collapse. The numerical technique discussed is particularly appropriate for the analysis of R.C. and P.C. beams rehabilitated or strengthened by means of external prestressing but it is reliable for the analysis of new beams as well.