• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.1
• /
• pp.346-358
• /
• 1996

A finite volume radiation solution method was applied to a non- orthogonal coordinate system for the analysis of radiative-convective heat transfer about a circular cylinder in crossflow. The crossflow Reynolds number based on the cylinder radius was 20, and the fluid Prandtl number was 0.7. The radiative heat transfer coupled with convection was reasonably predicted by the finite volume radiation solution method. The investigation includes the effects of conduction- to-radiation parameter, optical thickness, scattering albedo and cylinder wall-emissivity on heat transfer about the cylinder. As the conduction- to-radiation parameter decreases, the radiative heat transfer rate increases and conduction rate as well due to the increase in temperature gradient on the cylinder wall which is caused by radiation enhancement. With an increase in the optical thickness, the Nusselt number increases significantly and the temperature gradient shows similar behavior. Though the radiative heat transfer increases with the scattering albedo, the total heat transfer decreases. This is because the decrease in the conduction heat transfer exceeds the increase in the radiation heat transfer. As the wall- emissivity increases, the radiation absorbed in the vicinity of the cylinder wall increases and thereby the total heat transfer increases, even though the conduction heat transfer decreases.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.1
• /
• pp.135-140
• /
• 2003

The convergence of finite volume method (FVM) or discrete ordinate method (DOM) is known to degrade for optical thickness greater than unity and large scattering albedo. The present article presents a convergence acceleration procedure for the FVM based on a full approximation storage (FAS) multigrid method. Among a variety of multigrid cycles, the V-cycle is used and the full multigrid algorithm (FMG) is applied to an analysis of radiation in irregular two-dimensional geometry. Solution convergence is discussed for the several cases of various optical thickness and scattering albedo. At small scattering albedo and optical thickness, there is no advantage to using the multigrid method for calculation CPU time. For large scattering albedo greater than 0.5 and optical thickness greater than unity, however, the multigrid method improves the convergence and the solution is rapidly obtained.

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

The finite volume method for radiation is applied to investigate a radiative heating of rocket base plane due to searchlight and plume emissions. Exhaust plume is assumed to absorb, emit and scatter the radiant energy isotropically as well as anisotropically, while the medium between plume boundary and base plane is cold and nonparticipating. Scattering phase function is modelled by a finite series of Legendre polynomials. After validating benchmark solution by comparison with that of previous works obtained by the Monte-Carlo method, further investigations have been done by changing such various parameters as plume cone angle, scattering albedo, scattering phase function, optical radius and nozzle exit temperature. The results show that the base plane is predominantly heated by the plume emission rather than the searchlight emission when the nozzle exit temperature is the same as that of plume.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2681-2692
• /
• 1996

The finite-volume method for radiation in a three-dimensional non-orthogonal gas turbine combustion chamber with absorbing, emitting and anisotropically scattering medium is presented. The governing radiative transfer equation and its discretization equation using the step scheme are examined, while geometric relations which transform the Cartesian coordinate to a general body-fitted coordinate are provided to close the finite-volume formulation. The scattering phase function is modeled by a Legendre polynomial series. After a benchmark solution for three-dimensional rectangular combustor is obtained to validate the present formulation, a problem in three-dimensional non-orthogonal gas turbine combustor is investigated by changing such parameters as scattering albedo, scattering phase function and optical thickness. Heat flux in case of isotropic scattering is the same as that of non-scattering with specified heat generation in the medium. Forward scattering is found to produce higher radiative heat flux at hot and cold wall than backward scattering and optical thickness is also shown to play an important role in the problem. Results show that finite-volume method for radiation works well in orthogonal and non-orthogonal systems.

• Journal of computational fluids engineering
• /
• v.9 no.2
• /
• pp.43-51
• /
• 2004

Interaction between fluid flow and thermal radiation has received considerable attention due to its numerous applications in engineering field. In this case the thermofluid properties of radiating fluid vary with the variation of temperature field caused by absorption and emission of radiant heat. To analyze the radiation heat transfer in radiating fluid, the simultaneous solution of the radiative transfer equation (RTE) and the fluid dynamics equations is required. This means that the numerical procedure used for the RTE must be computationally efficient to permit its inclusion in the other submodels, and must be compatible with the other transport equations. The finite volume method (FVM) and the discrete ordinates method (DOM) are usually employed to simulate radiation problems in generalized coordinates. These two representative methods are examined and compared, especially in view of the numerical integration of the radiation intensity over solid angle. The FVM shows better accuracy than the DOM owing to less constraints of the selection of control angle.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.62-70
• /
• 2004

Interaction between fluid flow and thermal radiation has received considerable attention due to its numerous applications in engineering field. In this case the thermofluid properties of radiating fluid vary with the variation of temperature field caused by absorption and emission of radiant heat. To analyze the radiation heat transfer in radiating fluid, the simultaneous solution of the radiative transfer equation (RTE) and the fluid dynamics equations is required. This means that the numerical procedure used for the RTE must be computationally efficient to permit its inclusion in the other submodels, and must be compatible with the other transport equations. The finite volume method (FVM) and the discrete ordinates method (DOM) are usually employed to simulate radiation problems in generalized coordinates. These two representative methods are examined and compared, especially in view of the numerical integration of the radiation intensity over solid angle. The FVM shows better accuracy than the DOM owing to less constraints of the selection of control angle.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.2
• /
• pp.213-220
• /
• 1998

The conventional discrete-ordinates method (DOM) is modified and developed for the analysis of two-dimensional axisymmetric cylindrical enclosure with curved wall. The objective of the present work is to extend the capability of the conventional DOM into a general axisymmetric geometry like nozzle-shaped enclosure, by adopting the arbitrary control angle as was done in the finite-volume method (FVM), while keeping the same two-dimensional solution procedure as in the conventional DOM. The present method is validated by applying it to three different benchmark problems of axisymmetric enclosure containing absorbing, emitting and scattering medium. Results presented in this work not only support the solution accuracy, but also moderate efficiency in the numerical calculation of axisymmetric radiation problem.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.146-153
• /
• 2004

The Numerical study has been carried out to investigate the effects of chemical reaction and thermal radiation on the rocket plume flow-field at various altitudes. The theoretical formulation is based on the Navier-Stokes equations for compressible flows along with the infinitely fast chemistry and thermal radiation. The governing equations were solved by a finite volume fully-implicit TVD(Total Variation Diminishing) code which uses Roe's approximate Riemann solver and MUSCL(Monotone Upstream-centered Schemes for Conservation Laws) scheme. LU-SGS (Lower Upper Symmetric Gauss Seidel) method is used for the implicit solution strategy. An equilibrium chemistry module for hydrocarbon mixture with detailed thermo-chemical properties and a thermal radiation module for optically thin media were incorporated with the fluid dynamics code. In this study, kerosene-fueled rocket was assumed operating at O/F ratio of 2.34 with a nozzle expansion ratio of 6.14. Flight conditions considered were Mach number zero at ground level, Mach number 1.16 at altitude 5.06km and Mach number 2.9 at altitude 17.34km. Numerical results gave the understandings on the detailed plume structures at different altitude conditions. The diffusive effect of the thermal radiation on temperature field and the effect of chemical recombination during the expansion process could be also understood. By comparing the results from frozen flow and infinitely fast chemistry assumptions, the excess temperature of the exhaust gas resulting from the chemical recombination seems to be significant and cannot be neglected in the view point of performance, thermal protection and flow physics.