Numerical analysis of interaction between non-gray radiation and forced convection flow over a recess using the full-spectrum k-distribution method

In the present work, the interaction between non-gray radiation and forced convection in a laminar radiating gas flow over a recess including two backward and forward facing steps in a duct is investigated numerically. Distributions of absorption coefficients across the spectrum (50 cm^?1 < η < 20,000 cm^?1) are obtained from the HITRAN2008 database. The full-spectrum k-distribution method is used to account for non-gray radiation properties, while the gray radiation calculations are carried out using the Planck mean absorption coefficient. To find the divergence of radiative heat flux distribution, the radiative transfer equation is solved by the discrete ordinates method. The effects of radiation–conduction parameter, wall emissivity, scattering coefficient and recess length on heat transfer behaviors of the convection–radiation system are investigated for both gray and non-gray mediums. In addition, the results of gray medium are compared with non-gray results in order to judge if the differences between these two approaches are significant enough to justify the usage of non-gray models. Results show that for air mixture with 10 % CO₂ and 20 % H₂O, use of gray model for the radiative properties may cause significant errors and should be avoided.     

Transient radiative and conductive heat transfer in non-gray semitransparent two-dimensional media with mixed boundary conditions

This paper is devoted to transient heat transfer involving radiation and conduction. Considering a non-gray purely absorbing media, the radiative heat transfer equation (RTE) is solved iteratively with the Discrete Ordinates Method (DOM) using an exponential differencing scheme. The energy balance equation is used to compute temperature at each time step with the Crank–Nicholson technique. Energy equation is coupled to the RTE through the radiative source term. Both equations are discretized with finite differencing schemes. The energy conservation leads to the sparse system of linear equations A× T=B which is solved with a bi-conjugate stabilized gradient technique (BCSG). Validation of the model with different test cases is achieved and application to transient heating of glass is also studied.     

Combined radiative and convective heat transfer in a three-dimensional rectangular channel at different wall temperatures

 This paper deals with a numerical study of combined convective and radiative heat transfer in a three-dimensional rectangular duct with hydrodynamically and thermally developing laminar flow. The gas is assumed to be an incompressible, absorbing, emitting, isotropically scattering, gray medium. Isothermal, gray, diffuse boundary walls at different temperatures are assumed. The finite-volume method (FVM) is adopted to describe both convective and radiative heat transfer. The coupled continuity and momentum equations are solved by means of SIMPLER algorithm. Numerical results for the radiative flux show very good agreement with the available data. The effects of aspect ratio, optical thickness, scattering albedo and wall emissivity on the mean bulk temperature are also investigated. By splitting the heat flux into convective and radiative contributions, the relative importance of these components is assessed for a typical range of values of the parameters. Received on 9 February 1999     

Heat transfer in open cell polyurethane foam insulation

This paper study systematic investigates the combined conductive and non-gray radiative heat transfer of open cell polyurethane (PU) foam in the pressure range between 760 and 0.02 Torr. Direct transmission measurements are also taken using Fourier transform infrared (FTIR) spectrometer. In doing so, experimental results are obtained for the spectral extinction coefficient from 2.5 to 25 μm. In addition, the P-3 approximation method along with the box model is employed to calculate the non-gray radiative heat flux. The diffusion approximation method is also applied to calculated the radiative conductivity. Also tested herein are three samples with different cell sizes ranging from 330 to 147 μm. According to those results, the spectral extinction coefficient increases with a decrease of cell size, leading to a decrease of thermal conductivity. Moreover, evacuating the gases in the foam cells can reduce the thermal conductivity of the PU foam by as much as 75%. Furthermore, radiative heat transfer accounts for about 4% of total heat transfer at 760 Torr and increases to 20% at 0.02 Torr. Received on 20 April 1998     

Heat transfer to laminar flow of nongray gases through a circular tube

Analyses are presented for infrared radiative energy transfer in gases when other modes of energy transfer simultaneously occur. Fully developed laminar flow of an absorbing emitting gas in a circular tube is considered under the conditions of uniform wall heat flux. Nongray as well as gray formulations are presented, and results are obtained for illustrative cases. Appropriate limiting solutions of the governing equations are obtained and conduction-radiation interaction parameters are evaluated. The influence of variable wall emittance (gray and nongray) upon radiative energy transfer in nongray gases is investigated. In particular, nongray results are obtained, in the large path length limit, for the flow of CO₂ through stainless steel tubes of various compositions. Finally, a correlation is presented which can be utilized to extend all nongray results for the parallel plate geometry, already available in literature, to yield results for the corresponding case of a circular tube. This work was supported by the National Science Foundation through Grant No. GK-16755.     

Transformational-zone method of calculation of complex heat exchange during flow of optically active medium inside tube of diffuse grey surface

The paper presents analytical and experimental investigations of influence of radiative heat transfer on complex heat exchange during flow of optically active gas inside a pipe of diffusegrey properties. It was assumed that the pipe is heated from the outside by a constant heat flux and gas flowing inside is both absorbing and emitting and of small optical density. The influence of length and radiative properties of the pipe surface and of the gas temperature distribution on the wall and in the gas were analysed. The influence of radiative energy transfer on overall heat transfer coefficient was estimated. Mathematical model of radiative convective heat exchange in a system of one-dimensional temperature field, based on zone division method of Hottel and surface transformation, was verified numerically and experimentally. The results of numerical calculations were compared with experimental results obtained during carbone dioxide (CO₂) flow inside electrically heated ceramic tube. The set of nonlinear differential equations was solved by Runge-Kutta method with Hamming modification and with the use of separable-kernel method.     

Composite heat transfer in case of steady laminar flow of a gray fluid with large optical density past a horizontal plate embedded in a saturated porous medium

The present paper discusses the problem of composite heat transfer and viscous friction of a moving gray medium with large optical density. Expressions for temperature and velocity distributions and the ratio of the radiative component to convective component of heat flux are obtained. It is observed that for a given value ofB the ratio of radiative heat flux to convective heat flux is maximum at the edge of the boundary layer and tends to an asymptotic value as the boundary is reached. However, for a given value ofK δ, the ratio of heat fluxes increases with increase inB (the porous parameter). The results also show that as the wall temperature approaches the value of free stream temperature, the ratio of heat fluxes decreases.     

Thermal radiation effects on MHD flow past a semi-infinite inclined plate in the presence of mass diffusion

MHD mixed free-forced heat and mass convective steady incompressible laminar boundary layer flow of a gray optically thick electrically conducting viscous fluid past a semi-infinite inclined plate for high temperature and concentration differences is studied. A uniform magnetic field is applied perpendicular to the plate. The density of the fluid is assumed to reduce exponentially with temperature and concentration. The usual Boussinesq approximation is neglected due to the high temperature and concentration differences between the plate and the ambient fluid. The Rosseland approximation is used to describe the radiative heat flux in the energy equation. The boundary layer equations governing the flow are reduced to ordinary differential equations, which are numerically solved by applying an efficient technique. The effects of the density/temperature parameter n, the density/concentration parameter m, the local magnetic parameter M_x and the radiation parameter R are examined on the velocity, temperature and concentration distributions as well as the coefficients of skin-friction, heat flux and mass flux.     

Solution to radiative transfer equation in a 3-D rectangular enclosure due to discontinuous heat flow divergence

In this paper a new model and computer code is presented by considering singular and discontinuous heat flow divergence. A hybrid model including Smith’s WSGG model and Coppale and Vervish’s model is used for calculating gas radiative properties. Energy equation is solved simultaneously to reach temperature field which specify gas radiative properties. S ₈ order of discrete ordinate method is used to solve RTE. It is assumed that walls of enclosure are gray, diffuse and opaque with specified temperature. Boundary conditions are corrected in each iteration that change temperature field.     

Radiative heat flux characteristics of methane flames in oxy-fuel atmospheres

Oxy-fuel combustion is a promising alternative for power generation with CO₂ capture, where the fuel is burned in an atmosphere enriched with oxygen and CO₂ is used as a diluent. This type of combustion is characterised by uncommon characteristics in terms of thermal heat transfer budget as compared to air supported systems. The study presents experimental results of radiative heat flux along the flame axis and radiant fractions of non-premixed jet methane flames developing in oxy-fuel environments with oxygen concentrations ranging from 35% to 70%, as well as in air. The flames investigated have inlet Reynolds numbers from 468 to 2340. The data collected have highlighted the effects of the flame structure and thermo-chemical properties of oxy-fuel combustion on the heat flux radiated by the flames. It was first observed that peak heat flux increases considerably with oxygen concentration. More generally the radiant fraction increases with both increasing Reynolds number in the laminar regime and oxygen concentration. It was found that despite a difference in flame temperature, the radiative characteristics of the flames (heat flux distributions and radiant fraction) in air were similar to those with 35% O₂ in CO₂. The radiative properties of flames in oxy-fuel atmosphere with CO₂ as diluents appear to be dominated by the flame temperature.     

Numerical analysis of conservative unstructured discretisations for low Mach flows

Unstructured meshes allow easily representing complex geometries and to refine in regions of interest without adding control volumes in unnecessary regions. However, numerical schemes used on unstructured grids have to be properly defined in order to minimise numerical errors. An assessment of a low Mach algorithm for laminar and turbulent flows on unstructured meshes using collocated and staggered formulations is presented. For staggered formulations using cell‐centred velocity reconstructions, the standard first‐order method is shown to be inaccurate in low Mach flows on unstructured grids. A recently proposed least squares procedure for incompressible flows is extended to the low Mach regime and shown to significantly improve the behaviour of the algorithm. Regarding collocated discretisations, the odd–even pressure decoupling is handled through a kinetic energy conserving flux interpolation scheme. This approach is shown to efficiently handle variable‐density flows. Besides, different face interpolations schemes for unstructured meshes are analysed. A kinetic energy‐preserving scheme is applied to the momentum equations, namely, the symmetry‐preserving scheme. Furthermore, a new approach to define the far‐neighbouring nodes of the quadratic upstream interpolation for convective kinematics scheme is presented and analysed. The method is suitable for both structured and unstructured grids, either uniform or not. The proposed algorithm and the spatial schemes are assessed against a function reconstruction, a differentially heated cavity and a turbulent self‐igniting diffusion flame. It is shown that the proposed algorithm accurately represents unsteady variable‐density flows. Furthermore, the quadratic upstream interpolation for convective kinematics scheme shows close to second‐order behaviour on unstructured meshes, and the symmetry‐preserving is reliably used in all computations. Copyright © 2016 John Wiley & Sons, Ltd.     

Three-dimensional fracture analysis using tetrahedral enriched elements and fully unstructured mesh

In the absence of automated and customized methods and tools, some of today’s existing methods for solving three-dimensional fracture problems require comprehensive finite element meshing, labor-intensive analysis and post-processing efforts. In this study, a tetrahedral enriched element method and related applications are presented that demonstrate employment of fully unstructured tetrahedral meshes for general mixed-mode three-dimensional fracture problems. As in the case of hexahedral enriched elements, the tetrahedral enriched elements also alleviate the needs of pre- and post-processing the finite element model, allowing direct computation of stress intensity factors in the solution phase. In addition, when tetrahedral enriched elements are used, the crack front region can also be meshed using unstructured elements allowing direct use of automatic free-meshing programs. The applications presented are plane-strain central crack problem, mode-I surface crack in a plate, inclined penny-shaped crack, edge-cracked bar under constant heat flux and lens-shaped crack embedded in a large elastic body. The results obtained are in good comparative agreement with those available in the literature. Thus, it is concluded that the enriched tetrahedral elements can be applied efficiently and accurately on a general three-dimensional fracture problem allowing usage of fully unstructured finite element meshes.     

Transient radiative heat transfer within a suspension of coal particles undergoing steam gasification

Transient radiative heat transfer in chemical reacting media is examined for a non-isothermal, non-gray, absorbing, emitting, and Mie-scattering suspension of coal particles, whose radiative properties vary with time as the particles undergo shrinking by endothermic gasification. A numerical model that incorporates parallel filtered collision-based Monte Carlo ray tracing, finite volume method, and explicit Euler time integration scheme is formulated for solving the unsteady energy equation that couples the radiative heat flux with the chemical kinetics. Variation of radiative properties, attenuation characteristics, temperature profiles, and extent of the chemical reaction are reported as a function of time. It is found that radiation in the visible and near IR spectrum incident on a cloud of coal particles greater than 2.5m is more likely to be forward scattered than absorbed, but the opposite is true as the particles shrink below 1.3m. The medium becomes optically thinner as the particles shrink and this effect is more pronounced for smaller initial coal particles because these offer higher volume fraction to particle diameter ratio and, consequently, attain higher temperatures, reaction rates, and shrinking rates.

---

Adaptive unstructured finite element method for two-dimensional detonation simulations

A non-equilibrium reacting flow methodology has been added to a conservative, monotonic, compressible flow solver to allow numerical simulations of gas detonations. This flow solver incorporates unstructured dynamically adaptive meshes with the Finite Element Method – Flux Corrected Transport (FEM-FCT) scheme, which has shown excellent predictive capability of various non-reacting compressible flows. A two-step induction parameter model was used to model the combustion of the gas phase coupled with an energy release equation which was simulated with a point implicit finite element scheme. This combustion model was then applied to a two-dimensional detonation test case of a hypothetical fuel:oxygen mixture. The detonation simulation yielded two transverse waves which continued to propagate. This feature and the detonation shock speed mean and fluctuations were found to be grid-independent based on a resolution of about twenty elements within the average induction length. The resolution was efficiently achieved with the unstructured dynamically adaptive finite elements, which were three orders of magnitude less in number then required for uniform discretization. Received 26 August 1996 / Accepted 31 March 1997     

基于非结构化同位网格的SIMPLE算法

通过基于非结构化网格的有限体积法对二维稳态Navier—Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散；扩散项的离散采用二阶中心差分格式；对于压力-速度耦合利用SIMPLE算法进行处理；计算节点的布置采用同位网格技术，界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算，讨论了非结构化同位网格有限体积法在实现SIMPLE算法时，迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知，本文的方法是准确和可信的。     

On radiative heat transfer in a plane layer of an absorbing medium

Recently there has arisen increased interest in the study of radiative heat transfer between geometrically simple systems, both as autonomous problems and as elements of more complex problems.Problems of this kind have been treated by many authors [1–111 who have considered gray, diffusely emitting and absorbing boundaries and gray nonscattering media. In most cases these investigations were restricted either to the derivation of approximate formulas for the net radiative flux, without an exact analysis of the temperature distribution in the layer [5–7], or to numerical computation [1–4], In the latter case, with the exception of [8], which contains a numerical analysis for the case of optical symmetry, no attempt was made to analyze the effect of the optical properties of the boundaries on the temperature field in the layer.These papers can be divided into two groups according to the method of analysis used. The first group includes papers based on the integral equations of radiative transfer, with the corresponding integral analytical methods [1, 2], Similar in nature are [3, 4] which use the slab method, applicable to electrical-analog computation, as well as a recent paper [8] based on probability methods.The second group of papers [5–7] is based on the so-called differential methods. Of particular interest is [7], which develops these methods to an advanced degree. In several papers the problem of radiative transfer is analyzed in conjunction with more complex problems (cf., e.g. [10, 11]).In the present work we shall attempt to carry out an approximate analytical study of problems connected with radiative heat transfer in a plane layer of an absorbing, emitting, nonscattering gray medium with temperature-independent optical properties. The layer is bounded by two parallel, diffusely emitting and diffusely reflecting, isothermal, gray planes.The paper presents the fundamental formulation of the problem, which consists in: (a) the determination of the net heat flux on the basis of given temperature distribution (direct formulation), and (b) the determination of the temperature distribution on the basis of given distribution of the net radiative heat source per unit volume and boundary temperatures (inverse formulation). The analysis is based on integral methods appropriate to the integral equations which represent the net total and hemispherical radiation flux densities [12].The author would like to thank S. S. Kutateladze for his interest in this work.     

A nongray analysis of radiating–convecting rectangular plate-fin arrays

A complete analysis of heat transfer in rectangular fin arrays has been conducted taking into consideration convection and radiation at all surfaces, as well as radiative exchange between the fins and all neighboring surfaces. The analysis assumes diffuse nongray surfaces and uses the spectral optical properties of stainless steel AISI 430. Three different gray models were proposed in order to assess the effect of the gray assumption on the heat-transfer characteristics. Results are presented for the temperature and radiative-flux distributions along the fin, the radiative flux along the base, the contribution of the radiative component to the overall heat transfer, and the effectiveness of the fin array. Significant deviations in some of these results were found between the nongray model and the gray models. In general, convection was found to be the more effective mode of heat transfer in fin arrays and the effectiveness of the array decreases as the contribution of the radiative component increases.     

A gas‐kinetic scheme for the simulation of turbulent flows on unstructured meshes

This paper presents a numerical method for simulating turbulent flows via coupling the Boltzmann BGK equation with Spalart–Allmaras one equation turbulence model. Both the Boltzmann BGK equation and the turbulence model equation are carried out using the finite volume method on unstructured meshes, which is different from previous works on structured grid. The application of the gas‐kinetic scheme is extended to the simulation of turbulent flows with arbitrary geometries. The adaptive mesh refinement technique is also adopted to reduce the computational cost and improve the efficiency of meshes. To organize the unstructured mesh data structure efficiently, a non‐manifold hybrid mesh data structure is extended for polygonal cells. Numerical experiments are performed on incompressible flow over a smooth flat plate and compressible turbulent flows around a NACA 0012 airfoil using unstructured hybrid meshes. These numerical results are found to be in good agreement with experimental data and/or other numerical solutions, demonstrating the applicability of the proposed method to simulate both subsonic and transonic turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.