Thermohydrodynamic characteristics of combined double-diffusive radiation convection heat transfer in a cavity

This work deals with the numerical analysis of a radiating gas flow caused by both temperature and buoyancy concentration gradients in a square cavity; in this regard, the set of governing equations, including conservation of mass, momentum, species, and energy are solved by a numerical technique. In terms of radiation, since the fluid is considered as a semitransparent medium, the radiative term in the energy equation appears and is calculated by numerical solving of the radiative transfer equation (RTE). Furthermore, all of the surrounding cavity walls are considered to be opaque, gray, and diffuse with constant emissivity. All of the flow equations are solved by the finite difference method (FDM) and the RTE by the discrete ordinate one (DOM). In the present study, an attempt is made to verify the optical thickness effects on flow, thermal behavior, and mass transform in a cavity flow, such that reciprocating trends were seen in this manner. Our numerical results show that the thermal field in double-diffusive convection flow reaches very fast its steady-state situation in comparison to the concentration distribution. Besides, it is found that the thermohydrodynamic characteristics of a double-diffusive convection flow of a radiating gas are much affected by optical thickness.     

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.     

Multiphysics Simulation Combining Large-Eddy Simulation,Wall Heat Conduction and Radiative Energy Transfer to Predict Wall Temperature Induced by a Confined Premixed Swirling Flame

A multi-physics simulation combining large-eddy simulation, conjugate heat transfer and radiative heat transfer is used to predict the wall temperature field of a confined premixed swirling flame operating under atmospheric pressure. The combustion model accounts for the effect of enthalpy defect on the flame structure whose stabilization is here sensitive to the wall heat losses. The conjugate heat transfer is accounted for by solving the heat conduction within the combustor walls and with the Hybrid-Cell Neumann-Dirichlet coupling method, enabling to dynamically adapt the coupling period. The latter coupling procedure is enhanced to determine statistics (mean, RMS, \(\ldots \)) in a permanent regime accurately and efficiently thanks to an acceleration technique which is derived and validated. The exact radiative heat transfer equation is solved with an advanced Monte Carlo method with a local control of the statistical error. The coupled simulation is carried out with or without accounting for radiation. Excellent results for the wall temperature are achieved by the fully coupled simulation which are then further analyzed in terms of radiative effects, global energy budget and fluctuations of wall heat flux and temperature.     

Group analysis for the effects of the temperature-dependent fluid viscosity and chemical reaction on free convective heat and mass transfer

A group analysis has been carried out to study heat and mass transfer characteristics of an incompressible Newtonian fluid having a temperature-dependent viscosity over a vertical stretching surface in the presence of thermal radiation and a chemical reaction. The Rosseland approximation is used to describe the radiative heat flux in the energy equation. The vertical surface is assumed to be permeable so as to allow for possible wall suction or injection. The governing differential equations are derived and transformed using the Lie group analysis. The transformed equations are solved numerically by applying the Runge—Kutta—Gill scheme with the shooting technique. Favorable comparisons with previously published works on various special cases of the problem are obtained     

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 transfer in absorbing, emitting, and anisotropically scattering boundary-layer flows with reflecting boundary

The heat transfer in absorbing, emitting, and anisotropically scattering boundary-layer flows with reflecting boundary over a flat plate, over a 90-deg wedge, and in stagnation flow is solved by application of the Galerkin method with the particular solution boundary condition I _p (τ₀,ξ,?μ) of the equation of radiative transfer for an inhomogeneous term and the Box method. The exact integral expressions for the radiation part of this problem are developed. The coupling between convective and radiative heat transfer in boundary-layer flows is described by a set of nonlinear simultaneous equations including differential equations and integrodifferential equations. The Galerkin method and the particular solution boundary condition I _p (τ₀,ξ,?μ) are used to analyze the radiation part of the problem. The nonsimilar boundary-layer equations are solved by the Box method. The present numerical procedure solutions are compared in tables with the other exact treating results, the P-3, and P-1 approximation methods for the case of isotropically scattering boundary-layer flows. The effects of linearly anistropically scattering and reflecting surface are taken into account. It is found that the present method is a reliable and efficient numerical procedure and scattering leads to a reduction in the total heat flux. The influence of the forward-backward scattering parameter on the total heat flux decreases with the increase of the surface reflectivity.     

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.     

Numerical heat transfer in a cavity with a solar control coating deposited to a vertical semitransparent wall

A transient two‐dimensional computational model of combined natural convection, conduction, and radiation in a cavity with an aspect ratio of one, containing air as a laminar and non‐participating fluid, is presented. The cavity has two opaque adiabatic horizontal walls, one opaque isothermal vertical wall, and an opposite semitransparent wall, which consists of a 6‐mm glass sheet with a solar control coating of SnS–Cu_xS facing the cavity. The semitransparent wall also exchanges heat by convection and radiation from its external surface to the surroundings and allows solar radiation pass through into the interior of the cavity. The momentum and energy equations in the transient state were solved by finite differences using the alternating direction implicit (ADI) technique. The transient conduction equation and the radiative energy flux boundary conditions are coupled to these equations. The results in this paper are limited to the following conditions: 10⁴≤Gr≤10⁶, an isothermal vertical cold wall of 21°C, outside air temperatures in the range 30°C≤T₀≤40°C and incident solar radiation of AM2 (750 W m⁻²) normal to the semitransparent wall. The model allows calculation of the redistribution of the absorbed component of solar radiation to the inside and outside of the cavity. The influences of the time step and mesh size were considered. Using arguments of energy balance in the cavity, it was found that the percentage difference was less than 4 per cent, showing a possible total numerical error less than this number. For Gr=10⁶ a wave appeared in the upper side of the cavity, suggesting the influence of the boundary walls over the air flow inside the cavity. A Nusselt number correlation as a function of the Rayleigh number is presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Enhancement of heat and mass transfer in metal hydride beds with the addition of Al plates

A numerical study is made of transient heat and mass transfer in metal hydride beds in the hydriding process with the addition of internal aluminum plates. The two-dimensional equations governing the hydration kinetics, hydrogen flow and heat transfer are solved by using the iterative method based on the finite-volume technique. It is found that the heat transfer is enhanced by the installation of aluminum plates, and the ratio of the gap distance between the aluminum plates (H) and the thickness of the bed (L) emerges to be an important parameter. The reaction process is also strongly influenced by H/L. An optimal value of H/L exists to yield the fastest reaction rate, which is also shown to depend on other relevant parameters, such as the thickness of hydride bed and the inlet pressure. Received on 17 July 1998     

Numerical simulation of steady-state and transient natural convection in an isothermal open cubic cavity

In this work the numeric results, of the steady-state and transient heat transfer by natural convection in a horizontal isothermal open cubic cavity are presented. The most important assumptions in the mathematical formulation are two, the flow is laminar and the Boussinesq approximation is valid. The conservation equations in primitive variables are solved using the finite volume method and the SIMPLEC algorithm. The advective terms are approximated by the SMART scheme and the diffusive terms are approximated using the central differencing scheme. The results are obtained for a Rayleigh number range from 10⁴ to 10⁷.The numerical model predicted flow instabilities and Nusselt number oscillations for high Rayleigh numbers.     

Oblique stagnation-point flow and heat transfer towards a shrinking sheet with thermal radiation

An analysis is made of steady two-dimensional oblique stagnation-point flow and radiative heat transfer of an incompressible viscous fluid towards a shrinking sheet which is shrunk in its own plane with a velocity proportional to the distance from a fixed point. Here the axis of the stagnation flow and that of the shrinking sheet are not aligned. A similarity transformation reduces the Navier-Stokes equations to a set of non-linear ordinary differential equations and are solved numerically using a shooting technique. The analysis of the results obtained shows that multiple solutions exist for a certain range of the ratio of the shrinking velocity to the free stream velocity. The effect of non-alignment for the wall shear stress and the horizontal velocity components are discussed. Streamline patterns are also shown for shrinking at the sheet with aligned and non-aligned cases. It is found that the temperature at a point in the fluid decreases with increase in effective Prandtl number (Pr _eff ). The results pertaining to the present study indicate that as Pr _eff increases, the rate of heat transfer also increases. The reported results are in good agreement with the available published work in the literature.     

Combined radiation and natural convection in a rectangular cavity with a transparent wall and containing a non-participating fluid

This paper describes a numerical method for the study of combined natural convection and radiation in a rectangular, two-dimensional cavity containing a non-participating (i.e. transparent) fluid. One wall of the cavity is isothermal, being heated either by solar radiation or independently. The opposite wall is partially transparent, permitting radiation exchanges between the cavity and its surroundings and/or the Sun; that wall also exchanges heat by convection from its external surface to the surroundings. The other two walls are adiabatic: convection and radiation there are balanced, so that there is no heat transfer through those walls. The equations of motion and energy are solved by finite difference methods. Coupled to these equations are the radiative flux boundary conditions which are used to determine the temperature distribution along the non-isothermal walls. A two-band radiation model has been employed. Results are presented for a square cavity with a vertical hot wall at 150 °C, the ambient at 20 °C and 10⁴ ? Ra ? 3 × 10⁵, in the absence of direct insolation. The effects on the flow and heat transfer in the cavity of radiation and external convection have been examined. More extensive results will be presented in subsequent papers.     

Numerical investigation of premixed combustion in a porous burner with integrated heat exchanger

In this paper, we perform a numerical analysis of a two-dimensional axisymmetric problem arising in premixed combustion in a porous burner with integrated heat exchanger. The physical domain consists of two zones, porous and heat exchanger zones. Two dimensional Navier–Stokes equations, gas and solid energy equations, and chemical species transport equations are solved and heat release is described by a multistep kinetics mechanism. The solid matrix is modeled as a gray medium, and the finite volume method is used to solve the radiative transfer equation to calculate the local radiation source/sink in the solid phase energy equation. Special attention is given to model heat transfer between the hot gas and the heat exchanger tube. Thus, the corresponding terms are added to the energy equations of the flow and the solid matrix. Gas and solid temperature profiles and species mole fractions on the burner centerline, predicted 2D temperature fields, species concentrations and streamlines are presented. Calculated results for temperature profiles are compared to experimental data. It is shown that there is good agreement between the numerical solutions and the experimental data and it is concluded that the developed numerical program is an excellent tool to investigate combustion in porous burner.     

Unsteady mixed convection on the stagnation-point flow adjacent to a vertical plate with a magnetic field

An analysis is performed to study the unsteady combined forced and free convection flow (mixed convection flow) of a viscous incompressible electrically conducting fluid in the vicinity of an axisymmetric stagnation point adjacent to a heated vertical surface. The unsteadiness in the flow and temperature fields is due to the free stream velocity, which varies arbitrarily with time. Both constant wall temperature and constant heat flux conditions are considered in this analysis. By using suitable transformations, the Navier–Stokes and energy equations with four independent variables (x, y, z, t) are reduced to a system of partial differential equations with two independent variables (, ). These transformations also uncouple the momentum and energy equations resulting in a primary axisymmetric flow, in an energy equation dependent on the primary flow and in a buoyancy-induced secondary flow dependent on both primary flow and energy. The resulting system of partial differential equations has been solved numerically by using both implicit finite-difference scheme and differential-difference method. An interesting result is that for a decelerating free stream velocity, flow reversal occurs in the primary flow after certain instant of time and the magnetic field delays or prevents the flow reversal. The surface heat transfer and the surface shear stress in the primary flow increase with the magnetic field, but the surface shear stress in the buoyancy-induced secondary flow decreases. Further the heat transfer increases with the Prandtl number, but the surface shear stress in the secondary flow decreases.     

Unsteady MHD free convective heat and mass transfer from a vertical porous plate with Hall current,thermal radiation and chemical reaction effects

An analysis is presented to investigate the effects of chemical reaction, thermal radiation and heat generation or absorption on unsteady free convective heat and mass transfer along an infinite vertical porous plate in the presence of a transverse magnetic field and Hall current. The governing partial differential equations are formulated and transformed by using a similarity transformation into a system of ordinary differential equations. The resulting equations are solved numerically using a fourth‐order Runge–Kutta scheme along with the shooting method. The Rosseland approximation is used to describe the radiative heat flux in the energy equation. Numerical results for the velocity, temperature and concentration distributions are shown graphically for different parametric values. The effects of parameters on the local friction coefficients, the Nusselt number and Sherwood numbers are depicted in tabulated form. Copyright © 2009 John Wiley & Sons, Ltd.     

MHD flow and heat transfer for the upper-convected Maxwell fluid over a stretching sheet

In the present work, the effect of MHD flow and heat transfer within a boundary layer flow on an upper-convected Maxwell (UCM) fluid over a stretching sheet is examined. The governing boundary layer equations of motion and heat transfer are non-dimensionalized using suitable similarity variables and the resulting transformed, ordinary differential equations are then solved numerically by shooting technique with fourth order Runge–Kutta method. For a UCM fluid, a thinning of the boundary layer and a drop in wall skin friction coefficient is predicted to occur for higher the elastic number. The objective of the present work is to investigate the effect of Maxwell parameter β, magnetic parameter Mn and Prandtl number Pr on the temperature field above the sheet.     

A least-squares finite element method for time-dependent incompressible flows with thermal convection

The time-dependent Navier–Stokes equations and the energy balance equation for an incompressible, constant property fluid in the Boussinesq approximation are solved by a least-squares finite element method based on a velocity–pressure–vorticity–temperature–heat-flux ( u –P–ω–T– q ) formulation discretized by backward finite differencing in time. The discretization scheme leads to the minimization of the residual in the l²-norm for each time step. Isoparametric bilinear quadrilateral elements and reduced integration are employed. Three examples, thermally driven cavity flow at Rayleigh numbers up to 10⁶, lid-driven cavity flow at Reynolds numbers up to 10⁴ and flow over a square obstacle at Reynolds number 200, are presented to validate the method.     

Convective drying analysis of three-dimensional porous solid by mass lumping finite element technique

A numerical analysis of convective drying of a 3D porous solid of brick material is carried out using the finite element method and mass lumping technique. The energy equation and moisture transport equations for the porous solid are derived based on continuum approach following Whitaker’s theory of drying. The governing equations are solved using the Galerkin’s weighted residual method, which convert the governing equations into discretized form of matrix equations. The resulting capacitance matrices are made diagonal matrices by following the classical row-sum mass lumping technique. Hence with the use of the Eulerian time marching scheme, the final equations are reduced to simple algebraic equations, which can be solved directly without using an equation solver. The proposed numerical scheme is initially validated with experimental results for 1D drying problem and then tested by application to convective drying of 3D porous solid of brick material for four different aspect ratios obtained by varying the cross section of the solid. The mass lumping technique could correctly predict the wet bulb temperature of the solid under evaporative drying conditions. A parametric study carried out for three different values of convective heat transfer coefficients, 15, 30 and 45 W/m² K shows an increased drying rate with increase in area of cross section and convective heat transfer coefficient. The proposed numerical scheme could correctly predict the drying behavior shown in the form of temperature and moisture evolutions.     

Thermocapillary convection in an absorbing-scattering medium subjected to irradiation

The criteria for the onset of thermocapillary convection in a horizontal radiating fluid layer heated by an incident thermal radiative energy source are determined. The fluid layer is an absorbing and isotropically scattering medium confined between a free upper surface and an insulated rigid lower surface. Linear analysis is performed on the continuity, momentum, energy, and approximate radiative equations. The resulting disturbance equations are solved using a numerical optimization technique to obtain the eigenvalues governing the onset of convective motion. The influence of thermal radiation on the critical Marangoni number is examined. Attention is drawn to the physical significance of the heat transfer mode, gravitational force, the scattering effect, and the surface radiative properties. The conditions leading to the onset of convection are presented as functions of the optical thickness, scattering albedo, Planck number, surface emissivities, and transmissivities.     

Viscous Dissipation and Thermal Radiation Effects on Mixed Convection from a Vertical Plate in a Non-Darcy Porous Medium

The paper presents an investigation of the influence of thermal radiation and viscous dissipation on the mixed convective flow due to a vertical plate immersed in a non-Darcy porous medium saturated with a Newtonian fluid. The physical properties of the fluid are assumed to be constant. The Rosseland approximation is used to describe the radiative heat flux in the energy equation. The governing partial differential equations are transformed into a system of ordinary differential equations and solved numerically using a shooting method. The results are analyzed for the effects of various physical parameters such as viscous dissipation, thermal radiation, mixed convection parameters, and the modified Reynolds number on dynamics. The heat transfer coefficient is also tabulated for different values of the physical parameters.