• International Journal of Air-Conditioning and Refrigeration
• /
• v.10 no.3
• /
• pp.162-175
• /
• 2002

To investigate smoke movement with radiation in a room fires, a numerical and experi-mental analysis were performed. In this paper, results from a field model based on a self-developed SMEP (Smoke Movement Estimating Program) were compared with Stockier's ex-periment and the experiments on various sized pool fires in a room with door The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k- $\varepsilon$ turbulence model with buoyancy term. Also it solves the radiation equation using the S-N discrete ordinates method (DOM). The result of the cal-culated smoke temperature considering radiation effect has shown good agreement compared with the experimental data, although there are large discrepancy in the hot smoke layer be-tween the temperature predicted by the SMEP with only convection effect and obtained by the experimental result. This large discrepancy is caused from the radiation effect of $H_2O$ and $CO_2$ gas under smoke productions. Hence the radiation effect under smoke in fire is the point to be specially considered in order to produce more realistic result.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.239-244
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved by successively in the moving least square sense. Some weighing functions were tested in order to investigate the up-winding effect for the convection term. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as FVM.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.84-90
• /
• 2001

Secondary instability of flow past a circular cylinder is examined using direct numerical simulation at Reynolds number 220 and 250. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme, and the pressure and viscosity terms are applied by the 4th order central difference. In spanwise, Navier-Stokes equation is distributed using Spectral Method. The critical Reynolds number for this instability is found to be about Re=190. The secondary instability leads re three-dimensionality with a spanwise wavelength about 4 cylinder diameters at onset (A-mode). Results of three-dimensional effect in wake of a circular cylinder are represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Bulletin of the Korean Chemical Society
• /
• v.31 no.9
• /
• pp.2613-2617
• /
• 2010

Chloroprene is a chlorine substituent of 1,3-butadiene. Butadiene rubber (BR) and chloroprene rubber (CR) composites were thermally aged at 60, 70, 80, and $90^{\circ}C$ for 2 - 185 days in a convection oven and changes of the crosslink densities by the accelerated thermal aging were investigated. The crosslink densities increased as the aging time elapsed and as the aging temperature became higher. Degrees of the crosslink density changes of the BR composite were on the whole larger than those of the CR one except the short-term thermal aging at 60 and $70^{\circ}C$. The crosslink densities abnormally increased after themal aging at high temperatures for a long time. Activation energies for the crosslink density changes of the rubber composites tended to increase with increase of the aging time and the variation showed a local minimum. The activation energies of the CR composite were lower than those of the BR one. The experimental results were explained with a role of ligand of chlorine atom of CR in a zinc complex, steric hindrance by chlorine atom of CR, and oxidation of rubber chain.

• Nuclear Engineering and Technology
• /
• v.32 no.4
• /
• pp.379-394
• /
• 2000

As an in-vessel retention (IVR) design concept in coping with a severe accident in the nuclear power plant during which time a considerable amount of core material may melt, external cooling of the reactor vessel has been suggested to protect the lower head from overheating due to relocated material from the core. The efficiency of the ex-vessel management may be estimated by the thermal margin defined as the ratio of the critical heat flux (CHF）to the actual heat flux from the reactor vessel. Principal factors affecting the thermal margin calculation are the amount of heat to be transferred downward from the molten pool, variation of heat flux with the angular position, and the amount of removable heat by external cooling In this paper a thorough literature survey is made and relevant models and correlations are critically reviewed and applied in terms of their capabilities and uncertainties in estimating the thermal margin to potential failure of the vessel on account of the CHF Results of the thermal margin calculation are statistically treated and the associated uncertainties are quantitatively evaluated to shed light on the issues requiring further attention and study in the near term. Our results indicated a higher thermal margin at the bottom than at the top of the vessel accounting for the natural convection within the hemispherical molten debris pool in the lower plenum. The information obtained from this study will serve as the backbone in identifying the maximum heat removal capability and limitations of the IVR technology called the Cerium Attack Syndrome Immunization Structures (COASISO) being developed for next generation reactors.

• Journal of computational fluids engineering
• /
• v.14 no.3
• /
• pp.86-94
• /
• 2009

The non-conservative form of momentum equations is often used for some two-phase flow codes instead of a conservative form because of numerical convenience. Another non-conservative form, so called, a semi-conservative form can improve the numerical solution of these codes maintaining the numerical convenience. It is close to the conservative form but still maintains the feature of the non-conservative form. A semi-conservative form of the momentum equations and a non-conservative form of the momentum equations are implemented in CUPID[1] code. The numerical results of the semi-conservative and the non-conservative forms are compared against analytical solutions and the solutions of the FLUENT code that uses the conservative form. The results clearly showed that the semi-conservative form of the momentum equations provides better solutions than the non-conservative form, especially for heterogeneous two-phase flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.11
• /
• pp.3589-3597
• /
• 1996

The Reynolds-averaged Navier-Stokes equations with the equations of the k-.epsilon. turbulence model are solved numerically in a general curvilinear system for a three-dimensional turbulent flow around an Ahmed body. The simulation is especially aimed at the evaluation of three finite differencing schemes for the convection term, which include the upwind differencing scheme(UDS), the second order upwind differencing scheme(SOU scheme) and the QUICK scheme. The drag coefficient, the velocity and pressure fields are found to be changed considerably with the adopted finite differencing schemes. It is clearly demonstrated that the large difference between computation and experiment in the drag coefficient is due to relatively high predicted values of pressure drag from both front part and vertical rear end base. The results also show that the simulation with the QUICK or SOU scheme predicts fairly well the flow field and gives more accurate drag coefficient than other finite differencing scheme.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.900-907
• /
• 2000

The three-dimensional turbulent flow in curved pipes susceptible to flow-accelerated corrosion has been analyzed numerically to predict the pressure and shear stress distributions on the inner surface of the pipes. The analysis employs the body-fitted non-orthogonal curvilinear coordinate system and a standard $ {\kappa}-{\varepsilon}$ turbulence model with wall function method. The finite volume method is used to discretize the governing equations. The convection term is approximated by a high-resolution and bounded discretization scheme. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of a modified version of momentum interpolation scheme. The SIMPLE algorithm is employed for the pressure and velocity coupling. The numerical calculations have been performed for two curved pipes with different bend angles and curvature radii, and discussions have been made on the distributions of the primary and secondary flow velocities, pressure and shear stress on the inner surface of the pipe to examine applicability of the present analysis method. As the result it is seen that the method is effective to predict the susceptible systems or their local areas where the fluid velocity or local turbulence is so high that the structural integrity can be threatened by wall thinning degradation due to flow-accelerated corrosion.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.12
• /
• pp.1657-1668
• /
• 1998

In the analysis of a mushy solidification system with natural convection using a fixed grid method, the enthalpy method has been used to account for the release of latent heat. The variable viscosity, Darcy source, and hybrid methods have been employed for the velocity suppression in a mushy region. The choice of the values of solid viscosity and permeability constant in conjunction with the Darcy source term plays an important role in forming the location and shape of the phase boundaries. In this work the effects of these major parameters related to steady-state behavior in the system of mushy solidification are investigated through a simple test problem. The effective specific heat based on the spatial gradients of the enthalpy and temperature is adopted for the treatment of the release of latent heat. The effects of the Prandtl and Rayleigh numbers on the shape of mushy region are examined using the hybrid method.