重力作用下密闭空间内非球形多分散气溶胶浓度计算方法技术

本发明专利技术公开了一种重力作用下密闭空间内非球形多分散气溶胶浓度计算方法，包括以下步骤：S1：对密闭空间内的多分散气溶胶的离散处理；S2：根据计算每个尺寸区间内的气溶胶粒子的体积形状因子；S3：基于气溶胶粒子的体积形状因子计算结果，以及气溶胶粒子在下降过程中受到的重力和阻力，计算气溶胶重力沉降速度；S4：对密闭空间进行空间离散处理；S5：计算空间内t时刻的气溶胶浓度。本发明专利技术能够提高重力作用下的密闭空间内非球形多分散气溶胶粒子浓度计算的准确性，提高核电厂事故情况下放射性源项释放评估的精度。事故情况下放射性源项释放评估的精度。事故情况下放射性源项释放评估的精度。

【技术实现步骤摘要】
重力作用下密闭空间内非球形多分散气溶胶浓度计算方法


[0001]本专利技术属于核电工程领域的模拟技术，尤其涉及一种重力作用下密闭空 间内非球形多分散气溶胶浓度计算方法。

技术介绍

[0002]核电厂严重事故后放射性物质裂变产物从堆芯释放，并以气态、蒸汽、 气溶胶等形式存在，其中气溶胶为分散并悬浮在气体介质中的固体或液体小 颗粒，是裂变产物的主要存在形式之一。事故情况下产生的气溶胶多为非球 形多分散粒子，其在安全壳内的行为直接影响最终释放到环境中的放射性源 项，重力沉降是一种有效的气溶胶自然去除机理，准确计算重力作用下密闭 空间内非球形多分散气溶胶浓度对澄清严重事故最终的释放源项具有重要意 义。
[0003]重力沉降模型考虑颗粒受到重力与阻力，当颗粒在气体中的运动达到受 力平衡时，将以一定的速度进行沉降，即重力沉降速度。模型假设气溶胶粒 子为标准球体进行重力和阻力的计算，但实际事故工况下，产生的气溶胶粒 子多为非标准球体，此种计算方法不适用。因此有必要提出一种密闭空间内 非标准球形多分散气溶胶浓度计算方法。

技术实现思路

[0004]本专利技术针对现有技术中忽略不规则形状的影响而导致的气溶胶重力沉降 计算不准确的问题，提供了一种重力作用下密闭空间内非球形多分散气溶胶 浓度计算方法。
[0005]为实现上述目的，本专利技术的技术方案为：
[0006]一种重力作用下的密闭空间内非球形多分散气溶胶浓度计算方法，包括 以下步骤：
[0007]S1：对密闭空间内的多分散气溶胶的离散处理；
[0008]S2：根据计算每个尺寸区间内的气溶胶粒子的体积形状因 子；其中L为气溶胶粒子的长轴长度(m)，D为颗粒的短轴长度(m)，h为颗 粒的厚度(m)，θ为气溶胶粒子的自由下落处于稳定位置时粒子长轴与水平面 的夹角，根据实验测量多处于55
°
～90
°
区间，R为气溶胶粒子的半径(m)， 下标i表示气溶胶尺寸区间的编号；所述的气溶胶粒子的体积形状因子是指： 粒子在重力作用下在空间中下降至稳定位置时规则粒子与不规则粒子所受重 力和阻力之比，与粒子的长轴、短轴、厚度，下降处于稳定位置时与水平面 的夹角和等效规则球体的半径相关；
[0009]S3：基于气溶胶粒子的体积形状因子计算结果，以及气溶胶粒子在下降 过程中受到的重力和阻力，计算气溶胶重力沉降速度；
[0010]S4：对密闭空间进行空间离散处理；
[0011]S5：计算空间内t时刻的气溶胶浓度。
[0012]所述的多分散气溶胶的离散处理是指：根据气溶胶粒子尺寸将多分散气 溶胶划
分为多个尺寸区间，确定不同尺寸区间的代表性尺寸和该区间内气溶 胶的量，针对每一尺寸区间内气溶胶进行重力沉降计算。
[0013]气溶胶粒子尺寸符合对数正态分布，概率分布函数 其中d
p
为粒子直径，d
g
为粒子的几 何平均直径，σ
g
为几何标准偏差，确定划分的粒子尺寸区间总数为N
tot
，尺寸 区间[d
qmin
,d
qmax
]内划分的尺寸区间数量为
[0014]所述的稳定位置是指：颗粒在气体中处于自由沉降状态时粒子受到的重 力和阻力相等时所处的位置。
[0015]所述的气溶胶粒子的重力是指：气溶胶粒子在气体中沉降时由于地球的 吸引而受到的力，与气溶胶粒子的体积和密度相关。
[0016]所述的气溶胶粒子的阻力是指：气溶胶粒子在气体中处于沉降状态时， 气体对运动气溶胶粒子的流体阻力，与气体粘度、气溶胶粒子尺寸和粒子运 动速度相关。
[0017]所述的重力沉降速度是指：当气溶胶粒子只受到重力沉降机制的作用， 在气体中垂直水平面方向运动的速度。
[0018]所述步骤S3中，气溶胶粒子达到受力平衡后，根据气溶胶粒子在重力 沉降过程中受到的重力等于气溶胶粒子所受阻力得 到气溶胶粒子沉降速度，同时采用体积形状因子修正后沉降速度得到，非球 形气溶胶粒子重力沉降速度
[0019]其中，μ
g
为气体的动力粘度(Pa*s)，R为气溶胶粒子的半径(m)，v
G
为气 溶胶粒子的重力沉降速度(m/s)，Cu(R
i
)为坎宁滑移系数，R为气溶胶粒子的半 径(m)，ρ
P
为气溶胶粒子的密度(kg/m3)，g为重力加速度(m/s2)，下标i表示气 溶胶尺寸区间的编号。
[0020]所述的重力沉降速度计算基于气溶胶粒子的体积形状因子计算结果是指： 对于不规则的球体，气溶胶粒子受到的重力并不等于根据测量直径计算得到 的体积与密度的乘积，阻力并不等于根据测量直径计算得到并与气体动力粘 度、粒子速度和坎宁滑移系数相关的值，因此用体积形状因子进行修正。
[0021]所述的坎宁滑移系数是用于表征当气溶胶颗粒尺寸与气体分子平均自由 程较接近时非连续效应的影响，其中λ为气体 分子平均自由程。
[0022]所述的空间离散是指：根据空间与沉积表面的距离，将大空间离散为多 个区间，并确定每个区间内不同粒径气溶胶粒子占总空间内气溶胶粒子的份 额，针对每个空间区间进行重力沉降计算。
[0023]所述的空间内气溶胶浓度是指：根据空间内初始气溶胶浓度、每个空间 区间内气
溶胶份额、每个空间区间距离沉积面的距离和重力沉降速度计算得 到的空间内气溶胶分布。
[0024]根据其中Co0为空间内初始气溶胶浓度 (kg/m3),Co(t)为t时刻空间内气溶胶浓度(kg/m3)，t为重力沉降持续的时 间(s)，S为气溶胶粒子与沉积面之间的距离，Frac为粒子占空间内总气溶胶 粒子的份额，下标i表示气溶胶尺寸区间的编号，下标j表示空间区间的编 号。
[0025]本专利技术由于采用以上技术方案，使其与现有技术相比具有以下的优点和 积极效果：
[0026]本专利技术对多分散气溶胶合理离散并考虑气溶胶粒子的体积形状因子时， 气溶胶粒子的重力沉降计算精度显著提高，即本专利技术解决了密闭空间内多分 散非球形气溶胶粒子的重力沉降计算不准确的问题，使用本专利技术所涉及的重 力作用下密闭空间内非球形多分散气溶胶浓度计算方法，能够提高气溶胶粒 子重力沉降计算的准确性，提高核电厂事故情况下放射性源项释放评估的精 度。
附图说明
[0027]图1为本专利技术计算密闭空间内多分散非球形气溶胶浓度计算流程图；
[0028]图2为非球形颗粒示意图；
[0029]图3为空间离散示意图；
[0030]图4为本专利技术实施例的实施效果图。
具体实施方式
[0031]以下结合附图和具体实施例对本专利技术提出的一种重力作用下密闭空间内 非球形多分散气溶胶浓度计算方法作进一步详细说明。根据下面说明，本发 明的优点和特征将更清楚。
[0032]参看图1，一种重力作用下的密闭空间内非球形多分散气溶胶粒子浓度 计算方法，包括以本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种重力作用下密闭空间内非球形多分散气溶胶浓度计算方法，其特征在于，包括以下步骤：S1：对密闭空间内的多分散气溶胶的离散处理；S2：根据计算每个尺寸区间内的气溶胶粒子的体积形状因子；其中L为气溶胶粒子的长轴长度，D为颗粒的短轴长度，h为颗粒的厚度，θ为气溶胶粒子的自由下落处于稳定位置时粒子长轴与水平面的夹角，R为气溶胶粒子的半径，下标i表示气溶胶尺寸区间的编号；S3：基于气溶胶粒子的体积形状因子计算结果，以及气溶胶粒子在下降过程中受到的重力和阻力，计算气溶胶重力沉降速度；S4：对密闭空间进行空间离散处理；S5：计算空间内t时刻的气溶胶浓度。2.根据权利要求1所述的重力作用下密闭空间内非球形多分散气溶胶浓度计算方法，其特征在于，所述多分散气溶胶的离散处理为根据气溶胶粒子尺寸将多分散气溶胶划分为多个尺寸区间，确定不同尺寸区间的代表性尺寸和该区间内气溶胶的量，针对每一尺寸区间内气溶胶进行重力沉降计算。3.根据权利要求1所述的重力作用下密闭空间内非球形多分散气溶胶浓度计算方法，其特征在于，气溶胶粒子尺寸符合对数正态分布，概率分布函数其中d
p
为粒子直径，d
g
为粒子的几何平均直径，σ
g
为几何标准偏差，确定划分的粒子尺寸区间总数为N
tot
，尺寸区间[d
qmin
,d
qmax
]内划分的尺寸区间数量为4.根据权利要求1所述的重力作用下密闭空间内非...

【专利技术属性】
技术研发人员：佟立丽，何丽雯，胡培政，王善普，曹学武，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：