• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.5
• /
• pp.580-586
• /
• 2004

The ray effects of finite volume method (FVM) or discrete ordinate method (DOM) are known to show a non-physical oscillation in solution of radiative heat transfer on a boundary. This wiggling behavior is caused by the finite discretization of the continuous control angle. This article proposes a combined procedure of the Monte-Carlo and finite-volume method (CMCFVM) for solving radiative heat transfer in absorbing, emitting, and an-isotropically scattering medium with an isolated boundary heat source. To tackle the problem, which is especially pronounced in a medium with an isolated heat source, the CMCFVM is suggested here and successfully applied to a two-dimensional circular geometry.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.347-352
• /
• 2000

In the present study, DOM and FVM have been used to analyze the radiative heat transfer in an axisymmetric cylindrical enclosure with obstacles. Heat flux distributions on the wall of enclosure form DOM and FVM are compared to those from simplified zone analysis for a nonparticipating medium. The comparison of DOM and FVM is also presented. Results show that there is a good agreement between FVM and simplified zone analysis. In addition, the effect of the thickness of the obstacle on the results is considered. Heat flux distribution on the surface of the obstacle is also presented.

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

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.10 no.1
• /
• pp.95-107
• /
• 1998

Natural convection and radiation heat transfer in a square enclosure containing absorbing, emitting, and isotopically scattering(participating) media is studied numerically using the finite volume method. Various numerical methods are employed to analyze the radiative heat transfer. However, it is very difficult to choose the proper method. In present study, a finite volume method(FVM) and a discrete ordinates method(DOM) are compared in rectangular enclosure. The SIMPLER algorithm is used to solve the momentum and energy equations. Thermal and flow characteristics are investigated according to the variation of radiation parameters such as optical thickness and scattering albedo. The result shows that the accuracy and the computing time of FVM are better than those of DOM in regular geometry.

• 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.35 no.6
• /
• pp.561-567
• /
• 2011

The inverse radiation analysis of a two-dimensional irregular configuration using unstructured triangular meshes is presented. In this study, an enclosure filled with an absorbing, emitting and scattering medium with diffusely emitting and reflecting opaque boundaries is considered. The finite volume method is applied to solve the radiative transfer equation in order to simulate the measured incident radiation values which are used as input data for the inverse analysis. The conjugate gradient method is adopted for the estimation of wall emissivities by minimizing the objective function at each iteration step. To verify the performance of the unstructured grid system, we compare the results with those using a structured grid system for the two-dimensional lopsided shape. The effect of measurement errors on the estimation accuracy is also investigated.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.3
• /
• pp.806-819
• /
• 1995

A study of radiative heat heat transfer has been done in the non-orthogonal coordinate system utilizing the finite volume method and the P.1 approximation. Radiation of absorbing, emitting and scattering media in a concentric annulus has been solved using the non-orthogonal coordinate and the calculations were compared with the existing results. The results obtained from the analysis using the finite volume method are in good agreement with the existing calculations for all optical thicknesses. It was also shown that for only optically thick cases, P-1 approximation can be used in a non-orthogonal coordinate. Convective heat transfer analysis has been carried out to obtain the temperature fields in a cross flow around a circular cylinder and the finite volume method was applied in the non-orthogonal coordinate system to analyze radiative heat transfer. Effects of the optical thickness, the ratio of the surface temperature of the cylinder tot he free stream temperature, and the scattering albedo on radiation have been presented.

• 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 the Korean Society of Propulsion Engineers
• /
• v.6 no.1
• /
• pp.21-29
• /
• 2002

Hybrid propulsion systems provide many advantages in terms of stable operation and safety. However, classical hybrid rocket motors have lower fuel regression rate and combustion efficiency compared to solid propellant rocket motor. The recent research efforts are focused on the improvement of volume limitation and regression rate in the hybrid rocket engine. The present study has numerically investigated the combustion processes in the hybrid rocket engine. The turbulent combustion is represented by the eddy breakup model and Hiroyasu and Nagle and Strickland-Constable model are used for soot formation and soot oxidation. Radiative heat transfer is modeled by finite volume method. To reduce the uncertainties for convective heat transfer near solid fuel surface having strong blowing effect, the Low Reynolds number $\kappa-\varepsilon$ turbulent model is employed. Based on numerical results, the detailed discussion has been made for the turbulent combustion processes in the vortex hybrid rocket engine.

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