介质集有2个介质簇但只提供了1个_参与介质中的辐射传热和离散坐标法

93141a54d421b75c5066efc24b57b6ab.png

当热辐射穿过介质时,如果介质对辐射产生吸收、发射或散射等作用,则认为该介质参与了辐射过程。本文我们介绍了使用多个方向的离散集方法计算热辐射的散射模型,并深入讨论了离散化方法。

参与介质中的散射

当空气粒子与太阳辐照相互作用时,大气层就是参与介质。能量的吸收和释放是由于分子或原子尺度上的能量转移引起的,而散射则是因参与介质分子对入射能量的衍射、反射和折射(参考文献1)作用导致入射能量在空间中重新分布。下图示例显示了参与介质与入射辐射的相互作用。

0bba689fd5852150a80601d618f80ba2.png
参与介质与辐射相互作用的示例。

在 COMSOL Multiphysics® 软件的附加产品——传热模块中,我们可以通过离散坐标法 将辐射传递方程离散化。该方法基于辐射的离散方向数量以及用于数值积分的相关正交权重建立。一系列离散方向和关联的正交权重的组合称为正交集。下面,我们将解释和比较 COMSOL 软件中不同物理场接口的正交集,这些接口用于吸收介质中的辐射以及与参与介质中的辐射传热。

首先,以参与介质为例。辐射强度

equation?tex=I
equation?tex=%7B%5Cbf+s%7D 方向(这里也可称为单位矢量)传输的守恒方程,称为辐射传递方程(RTE),表示如下

(1)

equation?tex=%5Cbegin%7Bsplit%7D%7B%5Cbf+s%7D+%5Ccdot+%5Cnabla+I%28%7B%5Cbf+r%7D%2C%7B%5Cbf+s%7D%29+%26+%3D+%5Ckappa+I_%7B%5Crm+b%7D%28T%29+-+%28%5Ckappa+%2B+%5Csigma_%7B%5Crm+s%7D%29I%28%7B%5Cbf+r%7D%2C%7B%5Cbf+s%7D%29+%2B+%5Cmathcal%7BS%7D%28%7B%5Cbf+r%7D%2C%7B%5Cbf+s%7D%29%2C%5C%5C+%5Cmathcal%7BS%7D%28%7B%5Cbf+r%7D%2C%7B%5Cbf+s%7D%29+%26+%3D+%5Cfrac%7B%5Csigma_s%7D%7B4%5Cpi%7D+%5Cint+%5Climits_0%5E%7B4%5Cpi%7D%5Cphi%28%7B%5Cbf+s%7D%2C%7B%5Cbf+s%7D%27%29I%28%7B%5Cbf+r%7D%2C%7B%5Cbf+s%7D%27%29d%5COmega%27%2C+%5Cend%7Bsplit%7D%5C%5C 式中,
equation?tex=%5Ckappa 是吸收系数,
equation?tex=%5Csigma_%7B%5Crm+s%7D 是散射系数,
equation?tex=I_%7B%5Crm+b%7D 是黑体辐射强度。

equation?tex=d+%5COmega%3D%5Csin%5Cpsi+d+%5Cpsi+d+%5Cvarphi 为无穷小立体角,以球面度(sr)为单位,其中
equation?tex=%5Cpsi
equation?tex=%5Cvarphi 是球坐标系中
equation?tex=%7B%5Cbf+s%7D 方向(单位向量)的极角和方位角 (r,
equation?tex=%5Cpsi ,
equation?tex=%5Cvarphi ),如下图所示。在等式 1中,
equation?tex=%5Ckappa+I+%28%7B%5Cbf+r%7D%2C%7B%5Cbf+s%7D%29+ 数量的能量被吸收,
equation?tex=%5Ckappa+I_%7B%5Crm+b%7D+%28T%29+ 被发射;
equation?tex=%5Cmathcal%7BS%7D+%28%7B%5Cbf+r%7D%2C%7B%5Cbf+s%7D%29+ 是从任一方向散射到
equation?tex=+%7B%5Cbf+s%7D+ 方向的能量总量;
equation?tex=%5Csigma_%7B%5Crm+s%7D+I+%28%7B%5Cbf+r%7D%2C%7B%5Cbf+s%7D%29+ 是从
equation?tex=%7B%5Cbf+s%7D+ 方向散射至任一方向的能量总量。

对积分项

equation?tex=%5Cmathcal%7BS%7D+ 进行离散化是一项具有挑战性的任务。需要
equation?tex=n+ 对成组的离散方向
equation?tex=+%7B%5Cbf+s%7D_i 以及正交权重
equation?tex=%7B%5Cbf+%5Comega%7D_i ,即正交集
equation?tex=%5Cbig%28%28%5Comega_i%2C%7B%5Cbf+s%7D_i%29%5Cbig%29_%7B1+%5Cleq+i+%5Cleq+n%7D 。该正交集允许单位球面的数值离散化:

equation?tex=%5Cint+%5Climits_%7B0%7D%5E%7B%5Cpi%7D+%5Cint+%5Climits_%7B0%7D%5E%7B2%5Cpi%7D+%5Csin+%5Cpsi+d%5Cpsi+d%5Cvarphi+%3D+%5Cint+%5Climits_0%5E%7B4%5Cpi%7D+d+%5COmega+%5Capprox+%5Csum_i+%5CDelta+%5COmega_i+%5Capprox+%5Csum_i+%5Comega_i.%5C%5C


为了求解在参与介质中任一点沿离散方向

equation?tex=%28%7B%5Cbf+s%7D_i%29_%7B1%5Cleq+i+%5Cleq+n%7D+ 传输的入射能量或入射辐射
equation?tex=G ,求解等式(1)可以得到在离散方向
equation?tex=%7B%5Cbf+s%7D_i+ 传输的每个辐射强度
equation?tex=I_i%2C+1%5Cleq+i+%5Cleq+n 。入射辐射
equation?tex=G%28%7B%5Cbf+r%7D%29%3D%5Cint+I%28%7B%5Cbf+r%7D%2C%7B%5Cbf+s%7D%29+d%5COmega ,辐射热通量
equation?tex=%7B%5Cbf+q%7D_%7B%5Crm+r%7D%28%7B%5Cbf+r%7D%29%3D%5Cint+I%28%7B%5Cbf+r%7D%2C%7B%5Cbf+s%7D%29+%7B%5Cbf+s%7D+d%5COmega 以及散射项 mathcal{S} 的数值近似可由下式表示:

(2)

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值