漫射理论

漫射理论

  虽然解析理论从诸如Maxwell方程或波动方程出发，引入粒子的散射和吸收特性，并求出方差和相关函数这些统计量的适当的微分方程或积分方程．原则上说，这种理论考虑了多次散射、衍射和干涉效应，在数学上是严格的。但是，实际上这种理论不可能得到包含所有这些效应的解，真正能有用的解的各种理论也都是近似的，而且都有一定的参数使用范围．
  输运理论从能量的角度讨论波在包含散射粒子的介质中的传播问题。 输运理论的基本微分方程叫做输运方程，它等价于分子动力学理论以及中子输运理论中所用的Boltzmann（玻尔兹曼）方程.[1]漫射近似是输运理论一种较为成功的近似求懈方法．它是通过对介质内的漫射强度作Taylor展开并取其头两项的一种近似．[2] 在本章中，我们将从输运理论的基本概念出发，定义一些输运理论中重要的物理量，推导出输运理论基本的微分方程-输运方程，并推导了输运方程的漫射近似一漫射方程。

1. 输运理论中的基本物理量

1.1辐射强度

  输运理论中的基本物理量是辐射强度、平均辐射强度和辐射通量密度，其中最重要和最基本的量是辐射强度。 为了引进辐射强度，我们先引入一个新的物理量：辐射通量。辐射通量为单位时间内通过某一个曲面的辐射能。设辐射通量关于波长的分布函数为$P\lambda$，则波长属于区闻($\lambda$，$\lambda$+$d\lambda$)的辐射透量为$P\lambda$$d\lambda$。如果把整个频段上的辐射通量都积分起来，则总辐射通量为：辐射通量的度量单位为瓦．那么就有：
$$P={\int }_0^{\infty }P(\lambda )d\lambda \text{\hspace{1em}(1.1)}$$   辐射强度的度量单位是瓦$\cdot$米${}^{-}$${}^2$。球面度${}^{-}$${}^2$。它是波传播的输运理论中最基本的物理量之一。辐射强度的一种定义为单位立体角内，通过单位投影面积的辐射通量。 参考图(1.1)，辐射强度为：
$$I(\mathbf{r},\mathbf{s})=\frac{dP}{\cos \theta dad\omega d\mu }\text{\hspace{1em}(1.2)}$$  上式中，$dP$为在立体角$d$w内通过以单位矢量s${}_0$为法线的面元的辐射通量。辐射强度表示某个表面的辐射特征。 在空间r 点s 方向上的辐射强度可以用下述方式加以测量：取一个接收面$\Delta a$和接收立体角$\Delta w$都很小的接收器，把它对准方向 -s。如图(1.2)．于是，在$\Delta a$，$\Delta w$足够小的条件下。总的接收功率$P_r$除以$\Delta a$$\Delta w$则得到辐射强度。

图1.1：辐射强度I(r，辐射强度)及功率dP

图1.2：辐射强度的测量

   由于辐射强度是描述瞳面某一个方向上的辐射特征，而我们有时需要考察一下这个曲面综合的辐射特征，因此引入平均辐射强度 $U(r)$为：
$$U(\mathbf{r})=\frac{1}{4\pi }{\int }_{4\pi }I(\mathbf{r},\mathbf{s})d\omega \text{\hspace{1em}(1.3)}$$   另外，我们可以定义辐射通量密度 $F(r)$为单位面积上通过的所有方向上的辐射通量，根据辐射通量的定义，有数学表达式为：
$$F(r)={\int }_{4\pi }I(\mathbf{r},s)sd\omega \text{\hspace{1em}(1.4)}$$

1.2吸收系数、散射系数及相函数

   当波通过介质时，一部分入射能量被粒子吸收，另一部分入射能量被介质中的粒子所散射(一般来说散射并不是各向同性的)，余下来的部分按原来的传播方向继续前进。波的吸收、散射是波在介质中传播时所发生的普遍现象。在输运理论中用吸收系数、散射系数及相函数来描述这样的吸收和散射现象以及散射的各向异性。 下面我们给出一个简单的吸收系数、散射系数及相函数的定义。
   一、吸收系数${\mu }_a$
   除了真空，没有一种介质对任何波长的电磁波或其他类型的波是完全透明的。吸收是物质的普遍性质。下面我们来定义吸收系数：令平面波通过介质传播，由于介质对能量的吸收，经过薄层$dx$后，辐射强度从$I$减少到$I-d{I}_a$。Lambert[3]指出：$d{I}_a/I$应与吸收层的厚度$dx$成正比，即：
$$\frac{d{I}_a}{I}=-{\mu }_a\cdot dx\text{\hspace{1em}(1.5)}$$   其中比例系数便称为吸收系数。在这里，由于我们处理的是非均匀介质，因此吸收系数是空间坐标的函数。
   二、散射系数${\mu }_s$
   在混浊介质中，含有许多大的粒子，它们的大小数量级与入射波的波长可以比较。在这样的介质中，尽管粒子分布可能是均匀的，但当波通过介质时，一样会发生散射现象。 因此当波通过介质时不仅介质的吸收会使辐射强度光减弱，同时，介质的散射也是使辐射强度减弱的原因之一。设渡通过介质薄层$dx$后，由于散射而减弱的辐射强度为$-d{I}_s$，它与入射强度$I$，及薄层厚度$dx$成正比：$$\frac{d{I}_s}{I}=-{\mu }_s\cdot dx\text{\hspace{1em}(1.6)}$$   上式中的比倒系数${\mu }_s$称之为散射系数。同样，由于介质的非均匀性，散射系数也是空间坐标的函数。通常，我们把吸收系数和散射系数之和成为衰减系数或消光系数。$${\mu }_t={\mu }_a+{\mu }_s\text{\hspace{1em}(1.7)}$$  因此，吸收系数、散射系数和衰减系数是不互相独立的，通常我们只需要用其中的两个系数来描述介质的吸收和散射特性。
   三、相函数$p(\mathbf{o},\mathbf{i})$
  相函数来源于天文学的名称“相函数”(phase function)，那里指月相。在辐射传输理论中相函数$p(\mathbf{o},\mathbf{i})$描述的是介质中散射的非各向同性 ，通常是一个无量纲量。下面我们给出相函数的一个定义：假定有一单位振幅的线偏振平面电磁渡入射到介电常数为${\epsilon }_0$而磁导率为${\mu }_0$的介质中，它的电场强度为:
$$E_i(\mathbf{r})={\mathbf{e}}_i{e}^{ik\mathbf{i}\cdot \mathbf{r}}\text{\hspace{1em}(1.8)}$$   振幅$\left|E_i\right|=\left|{\mathbf{e}}_i\right|,k=\omega \sqrt{{\mu }_0{\epsilon }_0}=2\pi /\lambda$，$\lambda$是波在该介质中的波长，$i$是波传播方向上的单位矢量，$e_i$是波的偏振方向上的单位矢量。这个波入射到介质中的某一点上时，以该点为参考点，散射波以球面波的形式向四周传播散射，如图1.3所示，但在各个方向上的散射强度是不一样的。 在任一散射方向$\mathbf{o}$上离这个参考点距离为$R$的地方的电场强度$E_s$可以表示成下列形式：$$E_s(R)=\frac{f(\mathbf{o},\mathbf{i})}{R}{e}^{ik\mathbf{o}\cdot R}\text{\hspace{1em}(1.9)}$$  由于入射电场的振幅为单位振幅，因此，根据上式，我们可以把 ∣ f ( o , i ) ∣ 2 |f(\mathbf{o}, \mathbf{i})|^2