忽略多次散射和发射的普遍不带坐标辐射传输方程求解

    科技2025-11-07  12

    d I λ k λ ρ d s = − I λ \frac{dI_{ \lambda }}{k_{\lambda}\rho ds}=-I_{\lambda} kλρdsdIλ=Iλ

    以上,就是在忽略多次散射和发射的增量贡献时的比尔-布格-朗伯定律, k λ k_{\lambda} kλ为质量消光截面,为每单位质量的消光面积,为质量吸收截面和质量散射截面之和,单位是 c m − 1 cm^{-1} cm1. ρ \rho ρ为物质密度,二者之积为消光系数。

    如果单个粒子的消光截面表示为 σ e \sigma_{e} σe,单位体积的粒子数密度为 N N N,则消光系数也可以表示为 σ e N \sigma_{e}N σeN,可替换上述公式中的 k λ ρ k_{\lambda}\rho kλρ.

    下面,求解上述微分方程:

    (1)上述公式变形为: d I λ d s = − k ρ I λ \frac{dI_{\lambda}}{ds}=-k\rho I_{\lambda} dsdIλ=kρIλ 形如 d y d x = f ( x ) ∗ g ( y ) \frac{dy}{dx}=f(x)*g(y) dxdy=f(x)g(y)的一阶微分方程,属于变量可分离的一阶微分方程。利用变量可分离的一阶微分方程解法进行求解。

    (2) d I λ I λ = − k ρ d s \frac{dI_{\lambda}}{I_{\lambda}}=-k\rho ds IλdIλ=kρds,对等式左右两边进行积分,积分下限和上限分别为0, s 1 s_{1} s1,则有: ∫ 0 s 1 1 I λ d I λ = ∫ 0 s 1 − k ρ d s \int_{0}^{s_{1}} \frac{1}{I_{\lambda}} dI_{\lambda}=\int_{0}^{s_{1}}-k\rho ds 0s1Iλ1dIλ=0s1kρds l n I λ ∣ 0 s 1 = ∫ 0 s 1 − k ρ d s lnI_{\lambda} \mid_{0}^{s_{1}}=\int_{0}^{s_{1}}-k\rho ds lnIλ0s1=0s1kρds l n I λ ( s 1 ) − l n I λ ( 0 ) = ∫ 0 s 1 − k ρ d s lnI_{\lambda}(s_{1})-lnI_{\lambda}(0)=\int_{0}^{s_{1}}-k\rho ds lnIλ(s1)lnIλ(0)=0s1kρds l n I λ ( s 1 ) I λ ( 0 ) = ∫ 0 s 1 − k ρ d s ln\frac{I_{\lambda}(s_{1})}{I_{\lambda}(0)}=\int_{0}^{s_{1}}-k\rho ds lnIλ(0)Iλ(s1)=0s1kρds

    (3)对上式左右两侧同取自然底数,化简,得: I λ ( s 1 ) I λ ( 0 ) = e ∫ 0 s 1 − k ρ d s \frac{I_{\lambda}(s_{1})}{I_{\lambda}(0)}=e^{\int_{0}^{s_{1}}-k\rho ds} Iλ(0)Iλ(s1)=e0s1kρds

    (4)上式化简,即得:

    I λ ( s 1 ) = I λ ( 0 ) e ∫ 0 s 1 − k ρ d s I_{\lambda}(s_{1})=I_{\lambda}(0) e^{\int_{0}^{s_{1}}-k\rho ds} Iλ(s1)=Iλ(0)e0s1kρds

    人 生 苦 短 , 得 用 L a T e X 人生苦短,得用LaTeX LaTeX.

    2020.10.0821 : 00 2020.10.08 21:00 2020.10.0821:00

    Processed: 0.013, SQL: 8