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

• Nuclear Engineering and Technology
• /
• v.56 no.2
• /
• pp.520-525
• /
• 2024

The optimization design problem of nuclear reactor radiation shield is a typical multi-objective optimization problem with almost 10 sub-objectives and the sub-objectives are always demanded to be under tolerable limits. In this paper, a design method combining multi-objective optimization algorithms with paralleling discrete ordinate transportation code is developed and applied to shield design of the Savannah nuclear reactor. Three approaches are studied for light-weighted and compact design of radiation shield. Comparing with directly optimization with 10 objectives and the single-objective optimization, the approach by setting sub-objectives representing weight and volume as optimization objectives while treating other sub-objectives as constraints has the best performance, which is more suitable to reactor shield design.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.5
• /
• pp.1288-1300
• /
• 1994

Mathematical modeling and numerical calculation on the flow and thermal characteristics induced by fire in a partial enclosure are performed. The solution procedures include the Shvab-Zeldovich approximation for the physical transport equations, low Reynolds number k-.epsilon. model for the turbulent fluid flow and Discrete Ordinate method(DOM) to calculate the radiative heat transfer. PMMA(Polymethylmethacrylate) is adopted as a solid fuel. Two different cases are considered : combustions with and without gas radiation occuring in a open cavity for variable pyrolyzing location of PMMA. When the fire source is located at the left-wall, the flow region of flame gas is limited at the left-wall and ceiling and recirculation region of inlet gas is formulated at neat the floor. In case of neglecting the radiative heat transfer, more large flame size and higher temperature is predicted. It is essential to consider the radiative heat transfer for analysis of fire phenomenon.

• 한국연소학회:학술대회논문집
• /
• 1999.10a
• /
• pp.155-163
• /
• 1999

A numerical study for simulating a swirling pulverized coal combustion in axisymmetric geometry is done here by applying the weighted sum of gray gases model (WSGGM) approach with the discrete ordinate method (DOM) to model the radiative heat transfer equation. In the radiative transfer equation, the same polynomial equation and coefficients for weighting factors as those for gas are adopted for the coal/char particles as a function of partial pressure and particle temperature. The Eulerian balance equations for mass, momentum, energy, and species mass fractions are adopted with the standard ${\kappa}-{\varepsilon}$ turbulence model, whereas the Lagrangian approach is used for the particulate phase for soot. The eddydissipation model is employed for the reaction rate for gaseous mixture, and the single-step first-order reaction model for the devolatilization process for coal. By comparing the numerical results with experimental ones, the models used here are confirmed and found to be one of good alternatives for simulating the combustion as well as radiative characteristics.

• Journal of the Korean Society of Safety
• /
• v.15 no.1
• /
• pp.19-27
• /
• 2000

This paper numerically deals with combined heat transfer in a enclosure with a block. The block affected by hot wall is located centrally in the enclosure with a radiating gray gas. The discrete ordinate method(DOM) was used for solving the radiative transfer equation. Both laminar and turbulent cases were investigated for various Rayleigh number and standard k-$\varepsilon$ model was adopted to turbulent case. The effects of optical thickness, wall emissivity and fluid-solid thermal conductivity ratio are investigated on the flow and temperature fields. This study shows that as the wall emissivity decreases, the temperature distribution gradually becomes uniform and the heat transfer is reduced in enclosure. It is expected that this study can help to design the energy system related to the combined heat transfer and operate it safely.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.65-70
• /
• 1999

Radiative heating of a liquid rocket base plane due to plume emission is numerically investigated. Calculation of flow and temperature fields around rocket nozzle precedes and thereby realistic plume shape and temperature distribution inside the plume are obtained. Based on the calculated temperature field, radiative transfer equation is solved by discrete ordinate method. The averaged radiative heat flux reaching the base plane is about $5kW/m^2$ at the flight altitude of 10.9km. This value is small compared with radiative heat flux caused by constant-temperature (1500K) plume emission, but it is not negligibly small. At higher altitude (29.8km), view factor between the babe plane and the exhaust plume is increased due to the increased expansion angle of the plume. Nevertheless, the radiative heating disappears since the base plane is heated to high temperature (above 1000K) due to convective heat transfer.

• Laser Solutions
• /
• v.12 no.2
• /
• pp.17-25
• /
• 2009

Two methodologies for predicting thermal relaxation time of tissue subjected to pulsed laser irradiation is introduced by the calculation the optical penetration depth and by the investigation of the temperature diffusion behavior. First approach is that both x-axial and y-axial thermal relaxation times are predicted and they are superposed to achieve the thermal relaxation time (${\tau}_1$) for two-dimensional square tissue model. Another approach to achieve thermal relaxation time (${\tau}_2$) is measuring the time required for local temperature drop until $e^{-1}$ of the maximum laser induced heating.

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

• Environmental Sciences Bulletin of The Korean Environmental Sciences Society
• /
• v.3 no.1
• /
• pp.41-56
• /
• 1999

Since UV-B radiation measuring networks have not been established, numerical models which calculate the flux from other readily available meteorological measurements may play an important role. That is, such a problem can be solved by using parameterization models such as two stream approximation, the delta-Eddington method, doubling method, and discrete ordinate method. However, most UV-B radiative transfer models have not been validated with measurements, because such models are not intended as practical computational schemes for providing surface estimates of UV-B radiation. The main concern so far has been to demonstrate model sensitivity for cloudless skies. In particular, few have been concerned with real cloud information. Clouds and aerosols have generally been incorporated as constituents of particular atmospheric layers with specified optical depths and scattering properties. The parameterization model presented here is a combination of a detailed radiative transfer algorithm for a coludless sky radiative process and a more approximate scheme to handle cloud effects. The model input data requires a daily measurement of the total ozone amount plus a daily record of the amount and type of cloud in the atmosphere. Measurements for an examination of the models at the Department of Atmospheric Sciences, Pusan National University have been takenfrom February, 1995. These models can be used to calculate present and future fluxes where measurements have not been taken, and construct climatologies for the period before ozone depletion began.

• Atmosphere
• /
• v.23 no.1
• /
• pp.39-45
• /
• 2013

A parameterization for the scattering of longwave radiation by ice clouds has been developed based on spectral scattering property calculations with shapes and sizes of ice crystals. For this parameterization, the size distribution data by Fu (1996) and by Michell and Arnott (1994) are used. The shapes of ice crystal considered in this study are plate, solid column, hollow column, bullet-rosette, droxtal, aggregate, and spheroid. The properties of longwave scattering by ice crystals are presented as a function of the extinction coefficient, single-scattering albedo, and asymmetry factor. The heating rate and flux by the radiative parameterization model are calculated for wide range of ice crystal sizes, shapes, and optical thickness. The results are compared with the calculated results using a six-stream discrete ordinate scattering algorithm and Chou's method. The new method (with various habits and size distributions) provides a good simulation of the scattering properties and cooling rate in optically thin clouds (optical thickness < 5). Depending on the inclusion of scattering by ice clouds, the errors in the calculation of the cooling rates are significantly different.