An Upwinding Procedure for Numerical Radiation Heat Transfer

2000 ◽  
Author(s):  
Daniel R. Rousse ◽  
Guillaume Gautier ◽  
Jean-François Sacadura

Abstract This paper presents a skewed upwinding procedure for application to the Control Volume Finite Element Method (CVFEM) in the context of radiation heat transfer problems involving participating media. The proposed first order scheme is stable, economical, accurate and it inherently precludes the possibility of computing negative coefficients in the discretized algebraic equations while accounting for the direction of radiant propagation. The suggested first-order skew positive coefficients upwind scheme (SPCUS) is validated by application to several basic test problems, acknowledged by the radiative heat transfer community: its performance has proven to be excellent.

2011 ◽  
Vol 15 (3) ◽  
pp. 649-662 ◽  
Author(s):  
Raoudha Chaabane ◽  
Faouzi Askri ◽  
Ben Nasrallah

A new algorithm based on the lattice Boltzmann method (LBM) and the Control Volume Finite Element Method (CVFEM) is proposed as an hybrid solver for two dimensional transient conduction and radiation heat transfer problems in an optically emitting, absorbing and scattering medium. The LBM was used to solve the energy equation and the CVFEM was used to compute the radiative information. The advantages of the proposed methodology is to avoid problems that confronted when previous techniques are used to predict radiative heat transfer, essentially, in complex geometries and when there is scattering and/or non-black boundaries surfaces. This method combination, which is applied for the first time to solve this unsteady combined mode of heat transfer, has been found to accurately predict the effects of various thermo-physical parameters such as the scattering albedo, the conduction-radiation parameter and the extinction coefficient on temperature distribution. The results of the LBM-CVFEM combination were found to be in excellent agreement with the LBM-CDM (Collapsed Dimension Method)this proposed numerical approach include, among others, simple implementation on a computer, accurate CPU time, and capability of stable simulation.


Author(s):  
John Tencer ◽  
Kevin Carlberg ◽  
Roy Hogan ◽  
Marvin Larsen

Radiation heat transfer is an important phenomenon in many physical systems of practical interest. When participating media is important, the radiative transfer equation (RTE) must be solved for the radiative intensity as a function of location, time, direction, and wavelength. In many heat transfer applications, a quasi-steady assumption is valid. The dependence on wavelength is often treated through a weighted sum of gray gases type approach. The discrete ordinates method is the most common method for approximating the angular dependence. In the discrete ordinates method, the intensity is solved exactly for a finite number of discrete directions, and integrals over the angular space are accomplished through a quadrature rule. In this work, a projection-based model reduction approach is applied to the discrete ordinates method. A small number or ordinate directions are used to construct the reduced basis. The reduced model is then queried at the quadrature points for a high order quadrature in order to inexpensively approximate this highly accurate solution. This results in a much more accurate solution than can be achieved by the low-order quadrature alone. One-, two-, and three-dimensional test problems are presented.


Author(s):  
Singiresu S. Rao

A meshless local Petrov-Galerkin (MLPG) method is proposed to obtain the numerical solution of nonlinear heat transfer problems. The moving least squares scheme is generalized, to construct the field variable and its derivative continuously over the entire domain. The essential boundary conditions are enforced by the direct scheme. The radiation heat transfer coefficient is defined, and the nonlinear boundary value problem is solved as a sequence of linear problems each time updating the radiation heat transfer coefficient. The matrix formulation is used to drive the equations for a 3 dimensional nonlinear coupled radiation heat transfer problem. By using the MPLG method, along with the linearization of the nonlinear radiation problem, a new numerical approach is proposed to find the solution of the coupled heat transfer problem. A numerical study of the dimensionless size parameters for the quadrature and support domains is conducted to find the most appropriate values to ensure convergence of the nodal temperatures to the correct values quickly. Numerical examples are presented to illustrate the applicability and effectiveness of the proposed methodology for the solution of heat transfer problems involving radiation with different types of boundary conditions. In each case, the results obtained using the MLPG method are compared with those given by the FEM method for validation of the results.


2008 ◽  
Author(s):  
Mohammad Hadi Bordbar ◽  
Timo Hyppa¨nen

This paper describes the theoretical bases of the Radiative Exchange Method, a new numerical method for simulating radiation heat transfer. By considering radiative interaction between all points of the geometry and solving the radiation balance equation in a mesh structure coarser than the structure used in computational fluid flow calculation, this method is able to simulate radiative heat transfer in arbitrary 3D space with absorbing, emitting and scattering media surrounded by emitting, absorbing and reflecting surfaces. A new concept is introduced, that of the exchange factors between the different elements that are necessary for completing the radiative balance equation set. Using this method leads to a set of algebraic equations for the radiative outgoing power from each coarse cell being produced and the result of this set of equations was then used to calculate the volumetric radiative source term in the fine cell structure. The formulation of the exchange factor for a three-dimensional state and also a mesh size analysis that was conducted to optimize the accuracy and runtime are presented. The results of this model to simulate typical 3D furnace shape geometry, is verified by comparison with those of other numerical methods.


2014 ◽  
Vol 35 (2) ◽  
pp. 65-92 ◽  
Author(s):  
Paweł Kuczyński ◽  
Ryszard Białecki

Abstract The paper deals with a solution of radiation heat transfer problems in enclosures filled with nonparticipating medium using ray tracing on hierarchical ortho-Cartesian meshes. The idea behind the approach is that radiative heat transfer problems can be solved on much coarser grids than their counterparts from computational fluid dynamics (CFD). The resulting code is designed as an add-on to OpenFOAM, an open-source CFD program. Ortho-Cartesian mesh involving boundary elements is created based upon CFD mesh. Parametric non-uniform rational basis spline (NURBS) surfaces are used to define boundaries of the enclosure, allowing for dealing with domains of complex shapes. Algorithm for determining random, uniformly distributed locations of rays leaving NURBS surfaces is described. The paper presents results of test cases assuming gray diffusive walls. In the current version of the model the radiation is not absorbed within gases. However, the ultimate aim of the work is to upgrade the functionality of the model, to problems in absorbing, emitting and scattering medium projecting iteratively the results of radiative analysis on CFD mesh and CFD solution on radiative mesh.


Sign in / Sign up

Export Citation Format

Share Document