没有合适的资源?快使用搜索试试~ 我知道了~
在本文中,我们引入了一种有效的数值方法来表征由任意入射贝塞尔光束照射的随机离散粒子的多重散射。 具体地,完美地满足麦克斯韦方程组的贝塞尔光束的矢量表达式结合旋转欧拉角被用来表示任意入射的贝塞尔光束。 利用混合矢量有限元边界积分特征基函数方法来求解涉及多个具有随机分布的离散粒子的散射问题。 由于有限元方法的灵活性,所采用的方法可以方便地解决由均匀分布的均质颗粒,非均质颗粒和各向异性颗粒引起的多重散射问题。 包括一些数值结果,以说明所提方法的有效性和功能,并显示当随机离散粒子被贝塞尔光束照射时的散射行为。
资源推荐
资源详情
资源评论
Multiple scattering of arbitrarily incident Bessel beams
by random discrete particles
Zhiwei Cui,* Yiping Han, and Xia Ai
School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, China
*Corresponding author: zwcui@mail.xidian.edu.cn
Received August 29, 2013; revised September 27, 2013; accepted September 27, 2013;
posted October 1, 2013 (Doc. ID 196605); published October 24, 2013
In this paper, we introduce an efficient numerical method to characterize the multiple scattering by random
discrete particles illuminated by Bessel beams with arbitrary incidence. Specifically, the vector expressions of
Bessel beams that perfectly satisfy Maxwell’s equations in combination with rotation Euler angles are used to
represent the arbitrarily incident Bessel beams. A hybrid vector finite element–boundary integral–characteristic-
basis function method is utilized to formulate the scattering problems involving multiple discrete particles with a
random distribution. Due to the flexibility of the finite element method, the adopted method can conveniently deal
with the problems of multiple scattering by randomly distributed homogeneous particles, inhomogeneous
particles, and anisotropic particles. Some numerical results are included to illustrate the validity and capability
of the proposed method and to show the scattering behaviors of random discrete particles when they are
illuminated by Bessel beams. © 2013 Optical Society of America
OCIS codes: (290.5850) Scattering, particles; (290.4210) Multiple scattering; (140.0140) Lasers and laser
optics; (050.1755) Computational electromagnetic methods.
http://dx.doi.org/10.1364/JOSAA.30.002320
1. INTRODUCTION
The problem of multiple scattering by random media com-
posed of many discrete particles is an important subject of
research owing to the wide range of possible applications in
academic research and industry. Over the past few decades,
the multiple scattering of plane waves by random discrete
particles has been extensively investigated. Also, many meth-
odologies have been developed to analyze this problem, for
instance, multiple scattering theory [
1–4], radiative transfer
theory [
5,6], the T-matrix method [7–10], the sparse-matrix
canonical-grid method [
11,12], the characteristic basis function
method (CBFM) [
13,14], and the hybrid finite element–
boundary integral–characteristic-basis function method (FE-
BI-CBFM) [
15].
In recent years, with the development of laser sources and
the expansion of their applications, there has been a growing
interest in the study of multiple scattering by random discrete
particles illuminated by laser beams. For the case of an inci-
dent focused Gaussian beam, an early study was carried out
by Mackowski and Mishchenko [
16]. In that paper, they ap-
plied the T-matrix method to simulate the multiple scattering
of a Gaussian beam by multiple discrete particles with a ran-
dom distribution. However, they focused their attention on
the restricted case of an on-axis normally incident Gaussian
beam. Later, we generalized the CBFM based on surface
integral equations to include the case of illumination by an
arbitrarily incident focused Gaussian beam and applied it to
examine the scattering behavior of random discrete particles
[
14]. In our work, the incident Gaussian beam was described
by utilizing the Davis beam model in combination with rota-
tion Euler angles. Nevertheless, the scattering of arbitrarily
incident Bessel beams by randomly distributed particles
has not been reported. In fact, Bessel beams, as another kind
of laser beams, have attracted widespread attention in various
fields ever since their first introduction by Durnin [
17], be-
cause of their special characteristics of nondiffraction and
self-reconstruction. That is, these beams can maintain the
same intensity profile and the intensively localized intensity
distribution. Therefore, they will not suffer diffraction during
the wave propagation process within the Rayleigh distance.
Moreover, they are able to wholly reform at some distance
beyond the obstruction as long as the whole beam is not
blocked if such beams encounter an obstruction. It is thus im-
portant to understand their behavior upon interaction with
random media composed of many discrete particles.
In this work, we introduce an efficient numerical method to
characterize the multiple scattering by random discrete par-
ticles illuminated by Bessel beams. Because the zeroth-order
Bessel beam has the typical characteristics of nondiffraction
and self-reconstruction, and can be easily obtained in the lab-
oratory, here we only consider the zeroth-order Bessel beam.
Specifically, the arbitrarily incident beams are described by
the vector expressions of the zeroth-order Bessel beam
[
18,19] in combination with the rotation Euler angles [20].
The scattering problems involving multiple particles with
a random distribution are analyzed by adopting the hybrid
FE-BI-CBFM presented in the authors’ previous paper [
15].
This method exploits the unique features of the hybrid
finite element–boundary integral (FE-BI) method [
21–24] and,
more importantly, the unique features of random discrete par-
ticles. It is designed in such a manner that it first decomposes
the original problem into many subregions, where each sep-
arate particle is regarded as a subregion, and then it employs
the finite element method (FEM) to deal with each subregion.
The subregions are coupled to each other through the boun-
dary integral equations based on Green’s function. To reduce
2320 J. Opt. Soc. Am. A / Vol. 30, No. 11 / November 2013 Cui et al.
1084-7529/13/112320-08$15.00/0 © 2013 Optical Society of America
computational burdens, the resultant FE-BI matrix equation is
solved by utilizing the CBFM described in [
13]. Compared with
other methods, the adopted hybrid FE-BI-CBFM technique
can conveniently deal with the problems of multiple scatter-
ing by randomly distributed homogeneous particles, inhomo-
geneous particles, and anisotropic particles owing to the
flexibility of the FEM.
This paper is structured as follows. In Section
2, we first
give a brief description of the Monte Carlo method that is used
to generate random discrete particles. The mathematical ex-
pressions for describing the arbitrarily incident Bessel beams
are then given. Finally, we briefly recall the fundamental
aspects of the FE-BI-CBFM that is used to characterize the
multiple scattering by random discrete particles. Section
3
shows the numerical results of this work. Finally, Section
4
concludes the paper.
2. FORMULATION
A. Generation of Random Discrete Particles
For the purpose of the following numerical study, the posi-
tions of multiple discrete particles with a random distribution
need to be generated. In this work, we adopt the Monte Carlo
method [
25] to generate the random discrete particles and ob-
tain the positions of particles that are needed for computation
of the unknown fields on the surfaces of the particles. Without
any loss of generality, particle positions are generated
randomly in a cubic box and the particles are assumed to
be uniform spheres that may contain inhomogeneities. The
procedure of the Monte Carlo method is described as follows.
Consider a cubic box of side length l containing M uniform
spherical particles. Let r be the radius of the primary particle.
Then the fractional volume occupied by the particles is
given by
f
4πM
3l
3
r
3
: (1)
Initially, these M particles are placed randomly inside the
primary cell with no overlap. In each cycle every particle is
subject to be randomly displaced once. The acceptance of
its new position is according to whether it overlaps another
particle or not. The displacement is random and is not gov-
erned by any interparticle force except that it cannot pen-
etrate other particles. Detailed steps can be found in [
25]. As
an illustration, we use the Monte Carlo method to generate
125 randomly distributed particles, as shown in Fig.
1. These
particles are uniform spheres, and the fractional volume
is f 10%.
B. Description of the Incident Bessel Beams
To launch the arbitrarily incident Bessel beams in the adopted
FE-BI-CBFM, the key point is to use appropriate mathematical
expressions. For this work, we utilize the vector expressions
of the zeroth-order Bessel beam presented by Mishra [
18]in
combination with the rotation Euler angles [
20] to describe
the arbitrarily incident beams. For the sake of convenient
description, we define two Cartesian coordinate systems, as
illustrated in Fig.
2. The particles are attached to a global
coordinate system Oxyz, and the incident Bessel beam
propagates from the negative w to the positive w in the beam
coordinate system O
B
uvw, with the leading electric field
polarized in the positive u direction. The beam center is
located at the point O
B
, and the coordinates of O
B
in the global
coordinate system are x
0
;y
0
;z
0
. The frame system Oxyz can
be obtained from the beam system O
B
uvw by rotations
through Euler angles α; β; γ followed by a translation of
x
0
;y
0
;z
0
. Specifically, we first translate the beam center
x
0
;y
0
;z
0
to the origin of particle system Oxyz, then rotate
around the z axis an angle α, and then around the line of nodes
(temporary u axis) an angle β such that the temporary w axis
is in the direction of the beam axis, and finally around the w
axis an angle γ. To be specific, the relationship between these
two coordinate systems can be written as
"
x − x
0
y − y
0
z − z
0
#
A
"
u
v
w
#
; (2)
where the elements of the transformation matrix A,
A
"
a
11
a
12
a
13
a
21
a
22
a
23
a
31
a
32
a
33
#
; (3)
are defined by
Fig. 1. Illustration of 125 randomly distributed uniform spherical
particles generated with the Monte Carlo method.
Fig. 2. Geometry of Cartesian coordinates of the beam and particle.
Cui et al. Vol. 30, No. 11 / November 2013 / J. Opt. Soc. Am. A 2321
剩余7页未读,继续阅读
资源评论
weixin_38633897
- 粉丝: 11
- 资源: 972
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功