• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.474-481
• /
• 2001

This paper addresses a numerical method for predicting transient temperature distributions in the wall of a curved pipe subjected to internally thermal stratification flow. A simple and convenient numerical method of treating the unsteady conjugate heat transfer in the non-orthogonal coordinate systems is presented. The proposed method is implemented in a finite volume thermal-hydraulic computer code based on a cell-centered, non-staggered grid arrangement, the SIMPLEC algorithm, a higher-order bounded convection scheme, and the modified version of momentum interpolation method. Calculations are performed for the transient evolution of thermal stratification in two curved pipes, where the one has thick wall and the other has so thin wall that its presence can be negligible in the heat transfer analysis. The predicted results show that the thermally stratified flow and transient conjugate heat transfer in a curved pipe with a finite wall thickness can be satisfactorily analyzed by the present numerical method, and that the neglect of wall thickness in the prediction of pipe wall temperature distributions can provide unacceptably distorted results.

• 유체기계공업학회:학술대회논문집
• /
• 1998.12a
• /
• pp.42-48
• /
• 1998

Numerical study is presented for the analysis of three-dimensional incompressible turbulent flows in multiblade centrifugal fan. Reynolds-averaged Navier-Stokes equations with standard k - $\epsilon$ turbulence model are transformed to non-orthogonal curvilinear coordinates, and are discretized with finite volume approximations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. The computational area is divided into three blocks; core, impeller and scroll, which are linked by multi-block method. The flow inside of the fan is regarded as steady flow, and mathematical formula established from the cascade theory and empirical coefficient are employed to simulate tile flow through the impeller. From comparisons between the computational results and the experimental data, the validity of the mathematical formula for the blade forces was examined and good results were obtained qualitatively. Hence, we can get the flow characteristics of multi-blade centrifugal fan and it will be a corner stone of the development of the multiblade centrifugal fan.

• Journal of Korean Society for Atmospheric Environment
• /
• v.14 no.3
• /
• pp.219-228
• /
• 1998

In the present investigation, a numerical model developed for the prediction of the wind flow over complex terrain is validated by comparing with the field experiments. For the solution of the Reynolds - Averaged Clavier- stokes equations which are the governing equations of the microscale atmospheric flow, the model is constructed based on the finite-volume formulation and the SIMPLEC pressure-correction algorithm for the hydrodynamic computation. The boundary- fitted coordinate system is employed for the detailed depiction of topography. The boundary conditions and the modified turbulence constants suitable for an atmospheric boundary- layer are applied together with the k- s turbulence model. The full- scale experiments of Cooper's Ridge, Kettles Hill and Askervein Hill are chosen as the validation cases . Comparisons of the mean flow field between the field measurements and the predicted results show good agreement. In the simulation of the wind flow over Askervein Hill , the numerical model predicts the three dimensional flow separation in the downslope of the hill including the blockage effect due to neighboring hills . Such a flow behavior has not been simulated by the theoretical predictions. Therefore, the present model may offer the most accurate prediction of flow behavior in the leeside of the hill among the existing theoretical and numerical predictions.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.160-168
• /
• 1995

This paper describes computational efforts on the various energy and environmental problems using Patankar's SIMPLE method. The specific problems included in this study are : pollutant and flammable material dispersions in open and confined areas, aerator-induced flow in a lake for DO(dissolved oxygen) concentration, primary clarifier for water and waste water treatment, hood ventilation in workplace, cyclone and LNG combustors and Dow chlorination reactor. A control-volume based finite-difference method is employed together with the power-law scheme. The pressure-velocity coupling is resolved by the use of the revised version of SIMPLE, says SIMPLER and SIMPLEC. The Reynolds stresses are closed using the standard or the RNG $k-{\varepsilon}$ models. Turbulent reaction is modeled using two fast chemistry methods such as eddy breakup and conserved scalar models. Further, a nonequilibrium model is developed for the application of the chlorination process in the Dow reactor. Other important empirical models and physical insights appeared in this study are presented and discussed in a brief note. The computational method developed in this study is considered, in general, as a viable tool for the design and determination of the optimal condition of various engineering system of interest.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2157-2161
• /
• 2003

In this paper, the response surface method using three-dimensional Navier-Stokes analysis to optimize the shape of a forward-curved blades centrifugal fan, is described. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard k-e turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in forward-curved blades centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.93-108
• /
• 1997

This paper describes several computational results on the various energy and environmental problems using Patankar's SIMPLE method. The specific problems included in this study are : pollutant and flammable material dispersions in open and confined areas, buoyancy-driven flow in a lake, primary clarifier for water and waste water treatment, hood ventilation in workplace. cyclone combustor and Dow chlorination reactor. A control-volume based finite-difference method is employed together with the power-law scheme. The pressure-velocity coupling is resolved by the use of the revised version of SIMPLE, says SIMPLER and SIMPLEC. The Reynolds stresses are closed using the standard or RNG κ-ε models. A nonequilibrium turbulent reaction model is developed for the application of the chlorination process in the Dow thermal reactor. Other important empirical models and physical insights appeared in this study are presented and discussed in a brief note. The computational method developed in this study is considered, in general, as a viable tool for the design and determination of the optimal operating condition of various environmental engineering system of interest.

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

• The KSFM Journal of Fluid Machinery
• /
• v.2 no.1 s.2
• /
• pp.29-34
• /
• 1999

To develop the efficient numerical optimization method for the design of an airfoil, an evaluation of various methods coupled with two-dimensional Naviev-Stokes analysis is presented. Simplex method and Hook-Jeeves method we used as direct search methods, and steepest descent method, conjugate gradient method and DFP method are used as indirect search methods and are tested to determine the search direction. To determine the moving distance, the golden section method and cubic interpolation method are tested. The finite volume method is used to discretize two-dimensional Navier-Stokes equations, and SIMPLEC algorithm is used for a velocity-pressure correction method. For the optimal design of two-dimensional airfoil, maximum thickness, maximum ordinate of camber line and chordwise position of maximum ordinate are chosen as design variables, and the ratio of drag coefficient to lift coefficient is selected as an objective function. From the results, it is found that conjugate gradient method and cubic interpolation method are the most efficient for the determination of search direction and the moving distance, respectively.

• Journal of the Korean Society of Safety
• /
• v.10 no.3
• /
• pp.68-80
• /
• 1995

A series of parametric calculations have been performed in order to investigate the short-term and short-range plume and puff behavior of toxic gaseous and solid pollutant dispersion in an open atmosphere. The simulation is made by the use of the computer program developed by this laboratory, in which a control-volume based finite-difference method is used together with the SIMPLEC algorithm for the resolution of the pressure-velocity coupling appeared In Wavier-Stokes equation. The Reynolds stresses are solved by the standard two-equation k-$\varepsilon$ model modified for buoyancy together with the RNG(Renormalization Group) k-$\varepsilon$ model. The major parameters considered in this calculation are pollutant gas density and temperature, the relative velocity of pollutants to that of the surrounding atmospheric air, and particulate size and density together with the height released. The flow field is typically characterized by the formation of a strong recirculation region for the case of the low density gases such as $CH_4$ and air due to the strong buoyancy, while the flow is simply declining pattern toward the downstream ground for the case of heavy molecule like the $CH_2C1_2$and $CCl_4$, even for the high temperature, $200^{\circ}C$. The effect of gas temperature and velocity on the flow field together with the particle trajectory are presented and discussed in detail. In general, the results are physically acceptable and consistent.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.499-504
• /
• 2001

A validation of a 3D CFD model for predicting local subcooling of moderator in the vicinity of calandria tubes in a CANDU reactor is performed. The small scale moderator experiments performed at Sheridan Park Experimental Laboratory(SPEL) in Ontario, Canada[1] is used for the validation. Also a comparison is made between previous CFD analyses based on 2DMOTH and PHOENICS, and the current model analysis for the same SPEL experiment. For the current model, a set of grid structures for the same geometry as the experimental test section is generated and the momentum, heat and continuity equations are solved by CFX-4.3, a CFD code developed by AEA technology. The matrix of calandria tubes is simplified by the porous media approach. The standard $k-\varepsilon$ turbulence model associated with logarithmic wall treatment and SIMPLEC algorithm on the body fitted grid are used and buoyancy effects are accounted for by the Boussinesq approximation. For the test conditions simulated in this study, the flow pattern identified is a buoyancy-dominated flow, which is generated by the interaction between the dominant buoyancy force by heating and inertial momentum forces by the inlet jets. As a result, the current CFD moderator analysis model predicts the moderator temperature reasonably, and the maximum error against the experimental data is kept at less than $2.0^{\circ}C$ over the whole domain. The simulated velocity field matches with the visualization of SPEL experiments quite well.