放射性气体浓度的预测模型.doc

2011年数学建模 承 诺 书 我们仔细阅读了数学建模的竞赛规则。 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与本队以外的任何人（包括指导教师）研究、讨论与赛题有关的问题。 我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们愿意承担由此引起的一切后果。 我们的参赛报名号为： 参赛组别（本科或专科）： 参赛队员 (签名) ： 队员1： 队员2： 队员3： 2011年建模 编 号 专 用 页 参赛队伍的参赛号码：（请各个参赛队提前填写好）： 竞赛统一编号（由竞赛组委会送至评委团前编号）： 竞赛评阅编号（由竞赛评委团评阅前进行编号）： 2011年数学建模 题 目 核电站时间内通过空间域的气体流量，再求出相同时间内气体增量，利用质量守恒定律和气体总量与浓度的关系式求出气体浓度表达式。 问题2：可以把气体的扩散想象成连续的，建立空间坐标系利用高架连续点源扩散的高斯烟羽模型，最后对于该方程进行分析，优化并求解。 问题3：要求求上风和下风L公里处的放射性物质浓度，主要考虑风速和放射性物质扩散速度在空间中的矢量运算，代入问题2求出的模型即可 问题4：需要搜集核辐射源距我国东南海以及美国西海岸的距离、风速、风向及相关数据，根据问题2的模型，预测核辐射是否会对我国及美国造成影响。 三、模型假设 污染物浓度在y、z风向上分布为正态分布； 全部高度风速均匀稳定； 扩散中污染物是守恒的（不考虑转化）； 模型只考虑平面风，且假设垂直地面的风的风速为零； 污染源的源强是连续且均匀的，初始时刻云团内部的浓度、温度呈均匀分布； 扩散过程中不考虑云团内部温度的变化，忽略热传递、热对流及热辐射； 泄漏气体是理想气体，遵守理想气体状态方程； 在水平方向，大气扩散系数呈各向同性； 把地面看作镜面 取x轴为平均风速方向，整个扩散过程中风速的大小、方向保持不变，不随地点、时间变化而变化； 地面对泄漏气体起全反射作用，不发生吸收或吸附作用； 整个过程中，泄漏气体不发生沉降、分解，不发生任何化学反应等。 四、符号说明 C：放射性物质浓度 m ：扩散系数 grad：梯度 ：空间域 ：流量（i=1、2） div ：散度记号 Q：泄漏源释放的气体总量 ：单位强度的电源函数 h ：核泄漏口德有效高度 ：核泄漏气团受到热力和自身动力抬升而产生的附加高度 H：核泄漏点源的有效高度 、：泄漏气体在y、z方向分布的标准差 ：平均风速 ：源强（即泄漏速度） 五、模型建立和求解 5．1 放射性物质浓度的预测模型 令气体泄漏时刻记作t=0，则t时刻时，空间任一点的放射性物质浓度为C（x，y，z）。有假设，单位时间通过单位法向面积的流量与浓度成正比，得到 ， 其中m是扩散系数，grad表示梯度，负号表示由浓度高相浓度低的地方扩散。 考察空间域，在无风情况下，气体匀速在大气层中向四周扩散，则空间域是一个规则球体，其体积为V，包围的曲面为S，S的外法线向量为，则在内通过的流量为： 内全体的增量为： 由质量守恒定律，得： 又根据曲面积分的奥氏公式 其中是散度记号，由以上公式，再利用积分中值定理不难得到： ，（） （1） 这是无界区域的抛物线型微分方程，由假设1，初始条件为作用在坐标原点的电源函数，可记作： （2） 其中，Q表示泄漏源释放的气体总量，是单位强度的电源函数，可解得方程（1）满足条件（2）的解为 （3） 只要知道x，y，z，t，就可以用来描述核电站周边不同距离地区、不同时段放射性物质浓度。 5．2 核电站周边放射性物质浓度的变化情况 针对在有风速时，选取进行分析x轴，y轴在水平面内垂直于x轴，正方向在x轴的左侧，z轴垂直于水平面，向上为正。此坐标系中，烟流中心与x轴重合或烟流在oxy平面的投影为x轴。 图 1 其中核泄漏口的有效高度，在实际核扩散中核泄漏气团从泄漏口排出时，由于受到热力抬升和本身动力抬升进而产生的一个附加高度。表示核泄漏点源的有效源高。 （4） 由概率统计理论可以写出方差的表达式为： （5） 由假设可以写出源强的积分公式： （6） 式中：、为泄漏气体在y、z方向分布的标准差，单位为 m；C（x,y,z