基于GPU放射性核素海洋大气扩散研究.docVIP

基于GPU放射性核素海洋大气扩散研究 　　摘 要：目前对海洋环境中反应堆严重事故的研究较少，为给海上核应急提供参考，采用欧拉模型，针对海洋环境中放射性核素大气扩散问题，基于GPU编写了大气扩散模拟程序，选取133Xe和133I两种核素作为研究对象，分析核素在30公里范围内的扩散过程，并将计算数据实时显示。结果表明，133I的积分浓度在下风方向约为133Xe的一半，干沉降是133I总量减少的主导因素。基于GPU的程序较非GPU版本加速约15%。 　　关键词：反应堆事故；海洋大气扩散；欧拉模型；放射性核素；GPU 　　中图分类号：TL732 文献标志码：A 文章编号：2095-2945（2017）33-0001-04 　　引言 　　目前国内外对位于沿海或内陆核电站的放射性物质大气扩散研究已经比较成熟，建立了完善的应急处理机制。 　　近年，小型反应堆逐渐成为热门话题，我国研发出了多种应用于海洋核动力平台的堆型。但反应堆并不是绝对安全的，在严重事故的情况下，裂变产生的放射性核素有可能突破安全屏障大量进入到大气中。 　　本文假想位于海上的小型反应堆发生严重事故，对事故中放射性核素的大气扩散过程进行计算机模拟，期望为海上核应急提供参考。 　　为了能直观准确地反映大气扩散的详细过程，本文选用了污染物大气扩散模型中的欧拉模型。考虑到数值求解欧拉模型时计算量大、并行度高的特点，本文将GPU通用计算技术应用到大气扩散模拟研究中，一方面加速计算，降低计算成本；另一方面利于计算结果的实时显示，数据的直观性更强。 　　1 欧拉模型介绍 　　欧拉模型是一种网格模型，将需要研究的空间区域划分成小的网格，应用数值方法离散化求解关于浓度的对流-扩散方程。 　　相较于工业中较常使用的高斯模型，欧拉模型没有将烟羽视为一个整体，因此精度更高，能反映浓度随时间变化的整个过程，而不是直接通过公式得到时间积分浓度。 　　1.1 对流-扩散方程 　　欧拉模型中的对流-扩散方程是根据质量守恒定律和菲克定律导出的。这里定义在任意时刻、位置单位体积空气中某种核素的活度为该核素的浓度，记为c（x，y，z，t），单位Bq/m3。暂时不考虑核素活度的变化，并将空气视为不可压缩流体，则在直角坐标系中有以下方程描述核素的浓度变化[1]： 　　（1） 　　式中S为该点的源强，单位Bq/（m3?s）；u、v、w分别代表风速在x、y、z方向上的分量向量；Kx、Ky、Kz分别是核素在三个方向上的扩散参数，单位m2/s。 　　在网格节点上求解对流-扩散方程时，不直接将（1）式离散化，而是先忽略z方向上的物质交换，将Kx、Ky视为与横纵坐标无关，在高度相同的每一层内部做一次迭代： 　　（2） 　　之后计算第k层网格和上下两层网格在竖直方向的物质交换： 　　（3） 　　最底层和最上层网格的计算较为特殊，最底层网格需要考虑干沉降；最上层网格处于大气混合层顶部，通常认为污染物不会再向上传输，因此这两层网格节点的计算式为： 　　（2）～（5）式中下标i，j，k表示节点由东向西、由南向北、由低向高的序号，上标n表示经过n个时间步长Δt之后的值；上标n表示迭代过程的中间量；式（1）中的Kz由第k层与第k+1层之间的扩散参数Kz，k代替；vd为干沉降速度。 　　1.2 扩散参数的确定 　　x、y、z三个方向的扩散参数对于大气扩散模型是至关重要的，选取合适的计算公式能提高模拟的精度。本文采用Seinfeld和Pandis提出的公式计算x、y方向的扩散参数： 　　（6） 　　当大气稳定度为近中性时，Kz，k可由下式计算： 　　（7） 　　其中w？鄢为对流速度尺度，h为混合层高度，z为两层网格交界面的高度。 　　1.3 放射性核素的衰变 　　放射性核素在大气扩散过程中，每时每刻都在发生衰变，其浓度会进一步降低，因此需要在计算大气扩散的过程中加入对衰变量的计算。 　　对于衰变常数为λ的放射性核素，其在t时刻活度记为A（t），经过Δt时间的衰变后，其活度变为： 　　（8） 　　若该核素是另一核素衰变的子体，且母核素在严重事故中也被释放到大气中，情况就会变得复杂。设子核素的衰变常数为λ0，t时刻浓度为A0（t）；母核素衰变常数为λ1，t时刻浓度为A1（t）。若不考虑多级衰变，则母核素按照（8）式的规律衰变，而子核素经过Δt时间后浓度为： 　　（9） 　　可?子核素的浓度计算需要用到上一时刻母核素的浓度，因此在计算时要同时跟踪两种核素。另外数值计算中需要注意精度的问题，为了保证收敛，时间步长不能取得太长，而较短的时间步长内核素的衰变量及其微小，因此本文设定衰变的计算和扩散的计算不同步，在一种核素衰变了1%后计算一次衰变量，该时间间隔为： 　　（10）