Publication Cover
Numerical Heat Transfer, Part A: Applications
An International Journal of Computation and Methodology
Volume 71, 2017 - Issue 7
173
Views
6
CrossRef citations to date
0
Altmetric
Original Articles

Numerical analysis of combined-mode dual-phase-lag heat conduction and radiation in an absorbing, emitting, and scattering cylindrical medium

, &
Pages 769-788 | Received 04 Nov 2016, Accepted 16 Feb 2017, Published online: 27 Apr 2017
 

ABSTRACT

Combined-mode dual-phase-lag (DPL) heat conduction and radiation heat transfer is analyzed in a concentric cylindrical enclosure filled with a radiatively absorbing, emitting, and scattering medium. The governing energy equation is incorporated with volumetric radiation as a source term, essentially to take the effect of radiative heat flux into account. While the energy equation is solved using the lattice Boltzmann method (LBM), the finite volume method (FVM) is used to calculate the radiative information. To establish the accuracy of the proposed LBM formulation, the governing energy equation is also solved with the finite difference method (FDM). Thermal perturbation is caused by suddenly changing the temperature at the boundaries. Radial temperature distributions during transience as well as steady state (SS) are presented for a wide range of parameters such as lag ratio, extinction coefficient, scattering albedo, conduction–radiation (C-R) parameter, boundary emissivity, and radius ratio. Sample results are benchmarked with those available in the literature, and a good agreement between the present and reported results is found.

Nomenclature

A=

- area

B=

- lag ratio

b=

- number of propagation directions

C=

- speed of propagation

C1=

- coefficient

C2=

- coefficient

C3=

- coefficient 2ΓB

Cp=

- dimensionless specific heat at constant pressure

cp=

- specific heat at constant pressure

=

- dimensional propagation velocity in a lattice direction

fi=

- particle distribution function in the lattice direction

=

- equilibrium particle distribution function in particle direction

k=

- thermal conductivity

r=

- space variable

T=

- temperature

t=

- time variable

qc=

- conductive heat flux

qr=

- radiative heat flux

qt=

- total heat flux

Vs=

- speed of thermal wave

ΔV=

- volume of control volumes

TE=

- truncation error

N=

- conduction–radiation parameter

G=

- irradiation

I=

- intensity of radiation

D=

- directional weight

Nγ=

- number of angular divisions in polar space

Nδ=

- number of angular divisions in azimuthal space

Sp=

- source term at the node of the control volume

P=

- present control volume

=

- east, west, south, and north neighbor control volumes of P, respectively

=

- direction

=

- outward unit normal vector at face i

=

- unit normal vector at the wall toward the medium

=

- coefficient of discretized RTE in directions m and n at nodal point I

Greek=
α=

- thermal diffusivity

Γ=

- nondimensional thermal diffusivity

γ=

- parameter

λ=

- parameter

η=

- scaled nondimensional space variable

ζ=

- scaled nondimensional time

θ=

- scaled nondimensional temperature

ρ=

- density

τ=

- relaxation time in the BGK model

=

- dimensional phase lag for heat flow

=

- dimensional phase lag for temperature gradient

τq=

- dimensionless phase lag for heat flow

τT=

- dimensionless phase lag for temperature gradient

Δη=

- scaled nondimensional lattice size

Δζ=

- scaled nondimensional time step

σ=

- Stefan–Boltzmann constant

κa=

- absorption coefficient

β=

- extinction coefficient

ω=

- scattering albedo

ψ=

- source term

ε=

- emissivity

Ω=

- solid angle

γ=

- polar angle

δ=

- azimuthal angle

σs=

- scattering coefficient

Subscripts=
1=

- inner wall

2=

- outer wall

c=

- conductive

r=

- radiative

b=

- blackbody

i=

- direction

re=

- reference state of dimensionless parameter

Superscripts=
*=

- dimensional quantity

(1)=

- O(ξ) in the Chapman–Enskong expansion

(2)=

- O(ξ2) in the Chapman–Enskong expansion

mt=

- time level for FDM discretization

Nomenclature

A=

- area

B=

- lag ratio

b=

- number of propagation directions

C=

- speed of propagation

C1=

- coefficient

C2=

- coefficient

C3=

- coefficient 2ΓB

Cp=

- dimensionless specific heat at constant pressure

cp=

- specific heat at constant pressure

=

- dimensional propagation velocity in a lattice direction

fi=

- particle distribution function in the lattice direction

=

- equilibrium particle distribution function in particle direction

k=

- thermal conductivity

r=

- space variable

T=

- temperature

t=

- time variable

qc=

- conductive heat flux

qr=

- radiative heat flux

qt=

- total heat flux

Vs=

- speed of thermal wave

ΔV=

- volume of control volumes

TE=

- truncation error

N=

- conduction–radiation parameter

G=

- irradiation

I=

- intensity of radiation

D=

- directional weight

Nγ=

- number of angular divisions in polar space

Nδ=

- number of angular divisions in azimuthal space

Sp=

- source term at the node of the control volume

P=

- present control volume

=

- east, west, south, and north neighbor control volumes of P, respectively

=

- direction

=

- outward unit normal vector at face i

=

- unit normal vector at the wall toward the medium

=

- coefficient of discretized RTE in directions m and n at nodal point I

Greek=
α=

- thermal diffusivity

Γ=

- nondimensional thermal diffusivity

γ=

- parameter

λ=

- parameter

η=

- scaled nondimensional space variable

ζ=

- scaled nondimensional time

θ=

- scaled nondimensional temperature

ρ=

- density

τ=

- relaxation time in the BGK model

=

- dimensional phase lag for heat flow

=

- dimensional phase lag for temperature gradient

τq=

- dimensionless phase lag for heat flow

τT=

- dimensionless phase lag for temperature gradient

Δη=

- scaled nondimensional lattice size

Δζ=

- scaled nondimensional time step

σ=

- Stefan–Boltzmann constant

κa=

- absorption coefficient

β=

- extinction coefficient

ω=

- scattering albedo

ψ=

- source term

ε=

- emissivity

Ω=

- solid angle

γ=

- polar angle

δ=

- azimuthal angle

σs=

- scattering coefficient

Subscripts=
1=

- inner wall

2=

- outer wall

c=

- conductive

r=

- radiative

b=

- blackbody

i=

- direction

re=

- reference state of dimensionless parameter

Superscripts=
*=

- dimensional quantity

(1)=

- O(ξ) in the Chapman–Enskong expansion

(2)=

- O(ξ2) in the Chapman–Enskong expansion

mt=

- time level for FDM discretization

Log in via your institution

Log in to Taylor & Francis Online

PDF download + Online access

  • 48 hours access to article PDF & online version
  • Article PDF can be downloaded
  • Article PDF can be printed
USD 61.00 Add to cart

Issue Purchase

  • 30 days online access to complete issue
  • Article PDFs can be downloaded
  • Article PDFs can be printed
USD 716.00 Add to cart

* Local tax will be added as applicable

Related Research

People also read lists articles that other readers of this article have read.

Recommended articles lists articles that we recommend and is powered by our AI driven recommendation engine.

Cited by lists all citing articles based on Crossref citations.
Articles with the Crossref icon will open in a new tab.