Publication Cover
Numerical Heat Transfer, Part A: Applications
An International Journal of Computation and Methodology
Volume 69, 2016 - Issue 6
184
Views
23
CrossRef citations to date
0
Altmetric
Original Articles

Transient radiative transfer in two-dimensional graded index medium by Monte Carlo method combined with the time shift and superposition principle

, , &
Pages 574-588 | Received 07 May 2015, Accepted 20 Jul 2015, Published online: 04 Jan 2016
 

ABSTRACT

Transient radiative transfer (TRT) in a two-dimensional scattering medium with graded refractive index distribution subjected to a collimated short-pulse irradiation is solved by a modified Monte Carlo (MMC) method coupled with the time shift and superposition (TSS) principle. The boundaries are considered as Fresnel surfaces, the refractive index at the boundary mismatches with that of the surroundings, making the reflectivity at the boundary change with the incident directions. The incident pulse consists of two parts when it hits the boundary: bundles directly reflected by the outside boundary and bundles refracted into the medium. The accuracy of the present algorithm is confirmed first. Numerical results show that by using the TSS principle, the computational efficiency is greatly improved. Afterward, the TRT in the media with different graded refractive index distributions is investigated. The time-resolved reflectance and transmittance at different locations are given. Several trends on the time-resolved signals are observed and analyzed.

Nomenclature

a=

anisotropic factor

c0=

light speed in vacuum, ms−1

H(t)=

Heaviside step function

I=

radiation intensity, Wm−2 sr−1

L=

medium thickness in direction z, m

Nx=

number of sub-layers in the x direction

Ny=

number of sub-layers in the y direction

m0=

discrete time steps of the incident pulse

N=

total number of simulated bundles

n=

refractive index

r=

radial coordinate

r0=

radius of incident laser beam, mm

ri=

radius of the detector, mm

tp=

incident pulse width, s

=

dimensionless incident pulse width,

t*=

dimensionless time, t* = βc0t

Φ=

scattering phase function

μ=

direction cosine

β=

extinction coefficient, m−1

ω=

scattering albedo

Nomenclature

a=

anisotropic factor

c0=

light speed in vacuum, ms−1

H(t)=

Heaviside step function

I=

radiation intensity, Wm−2 sr−1

L=

medium thickness in direction z, m

Nx=

number of sub-layers in the x direction

Ny=

number of sub-layers in the y direction

m0=

discrete time steps of the incident pulse

N=

total number of simulated bundles

n=

refractive index

r=

radial coordinate

r0=

radius of incident laser beam, mm

ri=

radius of the detector, mm

tp=

incident pulse width, s

=

dimensionless incident pulse width,

t*=

dimensionless time, t* = βc0t

Φ=

scattering phase function

μ=

direction cosine

β=

extinction coefficient, m−1

ω=

scattering albedo

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.