Radiative slab heating analysis for various fuel gas compositions in an axial-fired reheating furnace

A transient radiative slab heating analysis was performed to investigate the effect of various fuel mixtures on the performance of an axial-fired reheating furnace. The various fuel mixtures tested were assumed to be attained by mixing COG (coke oven gas) and BFG (blast furnace gas), which are the two main byproduct gases found in the integrated steel mill industry. The numerical prediction of radiative heat transfer was calculated using an FVM radiation solving method, which is a well-known and efficient method for curvilinear coordinates. The WSGGM (weighted sum of gray gas model) was also adopted to calculate the radiative heat transfer in composition dependent media. The entire furnace was divided into fourteen sub-zones to calculate the radiative thermal characteristics of the furnace without flow field calculations. Each sub-zone was assumed to have homogeneous media and wall temperatures. All of the medium and wall temperatures were computed by calculating the overall heat balance using some relevant assumptions. The overall heat balance was satisfied when the net heat input equaled the three sources of heat loss in each sub-zone, wall loss, skid loss, and slab heating loss.     

Transient radiative heating characteristics of slabs in a walking beam type reheating furnace

Transient radiative heating characteristics of slabs in a walking beam type reheating furnace is predicted by the finite-volume method (FVM) for radiation. The FVM can calculate the radiative intensity absorbed and emitted by hot gas as well as emitted by the wall with curvilinear geometry. The non-gray weighted sum of gray gas model (WSGGM) which is more realistic than the gray gas model is used for better accurate prediction of gas radiation. The block-off procedure is applied to the treatment of the slabs inside which intensity has no meaning. Entire domain is divided into eight sub-zones to specify temperature distribution, and each sub-zone has different temperatures and the same species composition. Temperature field of a slab is acquired by solving the transient 3D heat conduction equation. Incident radiation flux into a slab is used for the boundary condition of the heat conduction equation governing the slab temperature. The movement of the slabs is taken into account and calculation is performed during the residence time of a slab in the furnace. The slab heating characteristics is also investigated for the various slab residence times. Main interest of this study is the transient variation of the average temperature and temperature non-uniformity of the slabs.     

Numerical investigation of coupled natural convection and radiation in a differentially heated cubic cavity filled with humid air. Effects of the cavity size

A numerical study of natural convection with surface and air/H₂O mixture radiation in a differentially heated cubic square cavity is presented. The coupled flow and heat transfers in the cavity are predicted by coupling a finite volume method with a spectral line weighted sum of gray gase model to describe gas radiative properties. The radiative transfer equation is solved by means of the discrete ordinate method. Simulations are performed at Ra?=?10⁶, considering different combinations of passive wall and/or gas radiation properties and different cavity length. It was found that in presence of a participative medium representative of building, cavity length has a strong influence on temperature and velocity fields which affect the global circulation and heat transfers in the cavity. For each steady-state solution, the convective and radiative contributions to the global heat transfer are discussed. More specifically, boundary layer thickness, thermal stratification parameter, and three-dimensional effects are compared to pure convective case results. The results suggest that radiative effects, often considered as negligible in view of the relatively low optical thickness, may not be neglected when trying to predict regime transitions.     

Numerical investigation of combustion with non-gray thermal radiation and soot formation effect in a liquid rocket engine

A numerical analysis was carried out in order to investigate the combustion and heat transfer characteristics in a liquid rocket engine in terms of non-gray thermal radiation and soot formation. Governing gas and droplet phase equations with PSIC model, turbulent combustion model with liquid kerosene fuel, soot formation, and non-gray thermal radiative equations are introduced. A radiation model was implemented in a compressible flow solver in order to investigate the effects of thermal radiation. The finite-volume method (FVM) was employed to solve the radiative transfer equation, and the weighted-sum-of-gray-gases model (WSGGM) was applied to model the radiation effect by a mixture of non-gray gases and gray soot particulates. After confirming the two-phase combustion behavior with soot distribution, the effects of the O/F ratio, wall temperature, and wall emissivity on the wall heat flux were investigated. It was found that the effects of soot formation and radiation are significant; as the O/F ratio increases, the wall temperature decreases. In addition, as the wall emissivity increases, the radiative heat flux on the wall increases.     

Spectrally correlated radiation and laminar forced convection in the entrance region of a circular duct

The two-dimensional combined radiative and convective transfer in emitting and absorbing real gases in the entrance region of a duct with a jump of wall temperature is studied. The axial propagation of radiation is taken into account in the analysis. The flow field and the energy equations are solved simultaneously and the radiative properties of the flowing gases, CO₂ or H₂O, are modeled by using either the narrow-band correlated-k model or the global absorption distribution function (ADF) model. The results are presented in terms of temperature and radiative power fields, and of the evolution of bulk temperatures and of heat transfer coefficients. Due to the axial component of the radiative flux, the gas is preheated or precooled before the change in wall temperature and this induces a persistent difference between the results of 1-D and 2-D radiation analyses. Some differences between CO₂ and H₂O temperature and radiative power profiles, due to the different structures of their spectra, are put in evidence. The ADF model, only suitable for gray walls, is shown to be less accurate when the gas is heated than when it is cooled.     

Inverse Boundary Design Problems in Enclosures With Non-Gray Media

In this work, the inverse analysis is applied to radiative heat transfer boundary design problems with non-gray media. The objective of the inverse problem is to find the power of the heaters on the heater surface that produces the desired output, that is, temperature and heat flux distribution over the design surface. The inverse problem is formulated as an optimization problem for minimization of an objective function, which is defined by the sum of the squared difference between estimated and desired heat flux distributions over the design surface. The non-gray optimization problem is solved using the conjugate gradient method, which is a gradient-based optimization method. The spectral line weighted-sum-of-gray-gases model (SLW) is used to account for non-gray gas radiation properties. The radiative transfer equation is solved by the discrete ordinates method combined with two models for simulation of non-gray media. Enclosures with diffuse and gray walls are considered. Radiation is assumed the dominant mode of heat transfer. Example problems including homogeneous/nonhomogeneous, isothermal/nonisothermal media are considered. The results obtained using the SLW model and the gray model are compared.     

Simultaneous radiation and conduction heat transfer in a graded index semitransparent slab with gray boundaries

The simultaneous radiation and conduction heat transfer in a semitransparent slab of absorbing-emitting gray medium is solved in this paper. The refractive index of the medium spatially varies in a linear relationship, and the two boundary walls are diffuse and gray. A curved ray tracing technique in combination with a pseudo-source adding method is employed to deduce the radiative intensities on gray walls. Resorting to some of the results presented by Ben Abdallah and Le Dez, an exact expression of the radiative flux in medium is deduced. The influences on the temperature and radiative flux fields are examined, which are caused by the refractive index distribution, absorbing coefficient, thermal conductivity and the boundary wall emissivities. The results display the significant influences of the refractive index distribution and boundary wall emissivities on the radiative flux and temperature in medium.     

Numerical analysis of entropy generation through non-grey gas radiation in a cylindrical annulus

This study is devoted to analyze the radiative heat transfer of non-grey gas confined in a cylindrical annulus with isothermal walls. The radiative heat transfer equation is resolved through the Ray Tracing method, which is associated to the statistical narrow bands correlated–k (SNBcK) model to compute the medium radiative properties. Special focus is given on the components of radiative entropy generation and its dependency on geometrical and thermodynamic parameters. The results show that entropy generation is greatly affected by gas and wall temperatures. Moreover, the dominance between wall radiative entropy generation and the volumetric one depends mainly on differences between gas and wall temperatures.     

Assessment of gas radiative property models in the presence of nongray particles

In this study, a radiation code based on the method of lines solution of the discrete ordinates method for the prediction of radiative heat transfer in nongray gaseous media is developed by incorporation of two different spectral gas radiative property models, banded spectral line-based weighted sum of gray gases (banded SLW) and gray wide band (GWB) approximation in the presence of nongray absorbing–emitting–scattering particles. The aim is to introduce an accurate and CPU efficient spectral gas radiation model, which is compatible with spectral fuel/ash particle property models. Input data required for the radiation code and its validation are provided from two combustion tests previously performed in a 300 kWt atmospheric bubbling fluidized bed combustor test rig burning low calorific value Turkish lignite with high volatile matter/fixed carbon (VM/FC) ratio in its own ash. The agreement between wall heat fluxes and source term predictions obtained by global and banded SLW models reveal that global SLW model can be converted to an accurate wide band gas model (banded SLW) which can directly be coupled with spectral particle radiation. Furthermore, assessment of GWB approximation by benchmarking its predictions against banded SLW model shows that GWB gives reasonable agreement with a higher CPU efficiency when the particle absorption coefficient is at least in the same order of magnitude with the gas absorption coefficient.     

Nongray inverse boundary design problems in enclosures at radiative equilibrium

In this paper, an inverse analysis is used to find an appropriate heat flux distribution over the heater surface of radiant enclosures, filled with nongray media at radiative equilibrium from the knowledge of desired (prespecified) temperature and heat flux distributions over the given design surface. Regular and irregular 2D enclosures filled with nongray combustive gas products are considered. Radiation is considered the dominant mode of heat transfer and the medium temperature is obtained from the energy equation. To evaluate the nongray behavior of the participating gases properly, the spectral‐line weighted‐sum‐of‐gray‐gases (SLW) model with updated correlations is used. The dependence of absorption coefficients and the weights of the SLW model on the temperature of the medium makes the inverse problem nonlinear and difficult to handle. Here, the inverse problem is formulated as an optimization problem and the Levenberg‐Marquardt method has been used to solve it. The finite volume method is exploited for the discretization of the energy equation and the spatial discretization of the radiative transfer equation (RTE). The discrete ordinates method (TN quadrature) is used for the angular discretization of RTE. Five test cases, including homogeneous and inhomogeneous media, are investigated to prove the ability of the present methodology for achieving the desired conditions.     

Non-gray radiative convective conductive modeling of a double glass window with a cavity filled with a mixture of absorbing gases

Coupled radiation and natural convection heat transfer occurs in vertical enclosures with walls at different temperatures filled with gas media. In glass window thermal insulation applications in hot climates, infrared absorbing gases appear as an alternative to improve their thermal performance. The thermal modeling of glass windows filled with non-gray absorbing gases is somewhat difficult due to the spectral variation of the absorption coefficients of the gases and the phenomena of natural convection. In this work, the cumulative wavenumber (CW) model is used to treat the spectral properties of mixtures of absorbing gases and the radiative transport equation is solved using CW model and the discrete ordinates method. Due to the range of temperature variation, the mixture of gases is considered as homogeneous. The absorption coefficients were obtained from the database HITRAN. First, the natural convection in a cavity with high aspect ratio is modeled using a CFD code and the local and global Nusselt numbers are computed and compared with available empirical correlations. Also, the flow pattern for different Rayleigh numbers is analyzed. Then, the heat transfer in the gas domain is approximated by a radiative conductive model with specified heat flux at boundaries which is equivalent to convective transport at the walls surroundings. The energy equation in its two-dimensional form is solved by the finite volume technique. Three types of gas mixtures, highly absorbing, medium and transparent are investigated, to determinate their effectiveness in reducing heat gain by the gas ambient. Reflective glasses are also considered. The numerical method to solve radiative heat transport equation in gray and non-gray participant media was validated previously. The temperatures distributions in the gas and the glass domain are computed and the thermal performance of the gas mixtures is evaluated and discussed. Also, comparison with pure radiative conductive model is shown.     

Weighted sum of gray gases model optimization for numerical investigations of processes inside pulverized coal-fired furnaces

The effects of the number of significant figures (NSF) in the interpolation polynomial coefficients (IPCs) of the weighted sum of gray gases model (WSGM) on results of numerical investigations and WSGM optimization were investigated. The investigation was conducted using numerical simulations of the processes inside a pulverized coal-fired furnace. The radiative properties of the gas phase were determined using the simple gray gas model (SG), two-term WSGM (W2), and three-term WSGM (W3). Ten sets of the IPCs with the same NSF were formed for every weighting coefficient in both W2 and W3. The average and maximal relative difference values of the flame temperatures, wall temperatures, and wall heat fluxes were determined. The investigation showed that the results of numerical investigations were affected by the NSF unless it exceeded certain value. The increase in the NSF did not necessarily lead to WSGM optimization. The combination of the NSF (CNSF) was the necessary requirement for WSGM optimization.     

Efficiency analysis of radiative slab heating in a walking-beam-type reheating furnace

The thermal efficiency of a reheating furnace was predicted by considering radiative heat transfer to the slabs and the furnace wall. The entire furnace was divided into fourteen sub-zones, and each sub-zone was assumed to be homogeneous in temperature distribution with one medium temperature and wall temperature, which were computed on the basis of the overall heat balance for all of the sub-zones. The thermal energy inflow, thermal energy outflow, heat generation by fuel combustion, heat loss by the skid system, and heat loss by radiation through the boundary of each sub-zone were considered to give the two temperatures of each sub-zone. The radiative heat transfer was solved by the FVM radiation method, and a blocked-off procedure was applied to the treatment of the slabs. The temperature field of a slab was calculated by solving the transient heat conduction equation with the boundary condition of impinging radiation heat flux from the hot combustion gas and furnace wall. Additionally, the slab heating characteristics and thermal behavior of the furnace were analyzed for various fuel feed conditions.     

Reconstruction of Temperature Distribution in the Combustion Region of a Non-Gray Medium

Two-dimensional temperature distribution in the combustion region of a radiant furnace with non-gray gaseous products is reconstructed in an inverse manner. The discrete ordinate method is used to solve the radiative transfer equation, and the non-gray behavior of gaseous medium is modeled by the spectral line weighted sum of gray gas model. Conjugate gradient method is employed for the inverse analysis. The domain of interest is divided into two regions: the unknown- and known-temperature zones. The temperature distribution in the unknown-temperature zone, near the flame, is recovered in an inverse manner by the measurement of heat fluxes over the opposite walls in the direction of flame diffusion. Effects of the size of unknown-temperature zone, variable concentrations, and measurement errors in temperature reconstruction are discussed.     

Statistical modeling of thermal radiation transfer in buoyant turbulent diffusion flames

A statistical (Monte Carlo) method for radiative heat transfer has been incorporated in CFD modeling of buoyant turbulent diffusion flames in stagnant air and in a cross-wind. The model and the computational tool have been developed and applied to simulate both burner flames with controlled fuel supply rate and in self-sustained pool fires with burning rates coupled with flame radiation. The gas–soot mixture was treated either as gray (using the effective absorption coefficient derived from total emissivity data or the Planck mean absorption coefficient) or as non-gray (using the weighed sum of gray gases model). The comparison of predicted radiative heat fluxes indicates applicability of the gray media assumption in modeling of thermal radiation in case of high soot content. The effect of turbulence-radiation interaction is approximately taken into account in calculation of radiation emission, which is corrected to allow for temperature self-correlation and absorption-temperature correlation. In modeling buoyant propane flames in still air above 0.3 m diameter burner, extensive comparison is presented of the predictions with the measurements of gas species concentrations, temperature, velocity and their turbulent fluctuations, and radiative heat fluxes obtained in flames with different heat release rates. Similar to previously published experimental data, the predicted burning rate of flames above the acetone pools exposed to flame radiation increases with the pool diameter and approaches a constant level for large pool sizes. The magnitude of predicted burning rates is shown to be in agreement with the reported measurements. Augmentation of burning rate of the pool fire in a cross-wind because of increased net radiative heat flux received by the fuel surface and non-monotonic dependence of burning rate on cross-wind velocity, subject to the pool diameter, is predicted. The statistical treatment of thermal radiation transfer has been found to be robust and computationally efficient.     

Mixed Convection and Non-Gray Radiation in a Horizontal Rectangular Duct

Numerical studies for fluid flow and heat transfer in a horizontal rectangular duct are carried out. The flow is considered to be laminar, hydrodynamically and thermally developing. Heat transfer by both forced and natural convection is taken into account. The radiation from the gas is modeled with weighted sum of gray gases (WSGG) model. While considering non-gray radiation with WSGG, the fluid is considered to be a mixture of CO₂ and H₂O. Simulations are carried out with lower wall temperature than the inlet temperature of the gas. The effect of buoyancy and radiation on bulk mean temperature and Nusselt number are studied. The effects of temperature dependent properties are discussed. Comparative studies are carried out among forced convection, mixed convection, gray and non-gray gas radiation. It is found from the simulations that the assumption of gray gas can produce an error of ±10% over a non-gray model with WSGG for the cases studied.     

Scaling anisotropic scattering in radiative transfer in three-dimensional nonhomogeneous media

Radiative heat transfer in three-dimensional nonhomogeneous participating medium was investigated by using REM² method. The anisotropic scattering phase function was dealt with the scaling technique based on delta function approximation. The three-dimensional scaled isotropic results were compared with the published anisotropic scattering computations. A good agreement between the scaled isotropic approaches and the anisotropic solutions was found. The effects of scattering albedo, forward fraction of phase function, and wall emissivity were discussed. It was found that, with the increase of the scattering albedo, the radiative heat flux increases for forward scattering media, but decreases for backward scattering media. The radiative heat flux is increased with the increases of forward fraction of phase function and wall emissivity. The emissive power at the center of a cubical nonhomogeneous medium in radiative equilibrium with gray diffuse walls equals to the averaged blackbody emissive power of the six walls.     

Enhancement of Heat Transfer with Porous/Solid Insert for Laminar Flow of a Participating Gas in a 3-D Square Duct

In recent years, porous or solid insert has been used in a duct for enhancing heat transfer in high temperature thermal equipment, where both convective and radiative heat transfer play a major role. In the present work, the study of heat transfer enhancement is carried out for flow through a square duct with a porous or a solid insert. Most of the analyses are carried out for a porous insert. The hydrodynamically developing flow field is solved using the Navier–Stokes equation and the Darcy–Brinkman model is considered for solving the flow in the porous region. The radiative heat transfer is included in the analysis by coupling the radiative transfer equation to the energy equation. The fluid considered is CO₂ with temperature dependent thermophysical properties. Both the fluid and the porous medium are considered as gray participating medium. The increase in heat transfer is analyzed by comparing the bulk mean temperature, Nusselt number, and radiative heat flux for different porous size and orientation, Reyonlds number, and Darcy number.     

Evaluation of emissivity correlations for H2OCO2N2/air mixtures and coupling with solution methods of the radiative transfer equation

This study is directed towards the limitations of applying total emissivity correlations in computational fluid dynamics (CFD) computer codes for flame modeling. The predictions of nine widely applied total emissivity models for H₂O---CO₂ homogeneous mixtures are compared with the exponential wide band model (EWBM) calculations. The comparison covers a range of total pressures, temperatures and path lengths which are suitable for the use of fine numerical grids in CFD simulations of atmospheric and high pressure combustors.

Attention is paid to coupling of the property models with the radiative transfer equation (RTE) and their performance in non-homogeneous applications. In this respect both the total transmittance non-homogeneous (TTNH) model and the spectral group model (SGM) are used. The latter model is combined with five weighted sum of gray gases models (WSGGM), the single line based sum of gray gases model (SLW) and the k-distribution model. The non-homogeneous validation tests used in situ total radiance measurements in two non-luminous natural gas flames representing two industrial situations, a water cooled furnace and a refractory lined furnace. The main conclusions are as follows. The spectral group model provides an elegant and accurate method of coupling WSGGM, k-distribution and SLW property models to the equation of radiative transfer. Both homogeneous and non-homogeneous tests indicate the advantage of using the Smith, Shen and Friedman weighted sum of gray gases model over polynomial correlations and the SLW model. It has been shown that in the near burner region of a natural gas diffusion flame, the water vapor to carbon dioxide partial pressure ratio departs significantly from the value expected for the complete combustion of methane in air. This finding emphasizes the limitation of existing WSGGM to H₂O to CO₂ partial pressure ratios of one and two only.