激活扩散模型

放射气体模型的预估模型

摘要

本文是以日本福岛核电站遭遇自然灾害发生核泄漏的背景而提出的。且结合了高斯烟羽模型、线性拟合，以及微分方程模型，运用MATLAB软件，分析泄漏源强度、风速、大气稳定度参数等因素对放射性气体扩散的影响，预测了放射性气体浓度在不同时间，不同地区的浓度变化，并且本文模型中数据可以根据不同的实际情况而加以改变，因而是本文的应用范围大大增加，可以适用于具有较强的应用型。

对于问题一，讨论在无风的情况下，放射性气体以s m/s的匀速在大气中向四周扩散。本问中由于不考虑风力的影响，且扩散出来的气体匀速向四周散开，这样经过任意时刻t，扩散的气体围成一个半径为st的球，且距球心位置不同的地方浓度值不同。采用列数列的表现方法，设定相同时间段t，把条件进行整理，并经过简单计算得出每段时间所预测得到的扩散距离r和浓度C。利用MATLAB软件对数据进行线性拟合，采用微分方程模型得到核电站周边放射性气体在不同地区，不同时间段的浓度变化，得出随着离泄漏源距离的延伸，最后放射性物质的浓度越来越小，趋近于零，即当x趋向无穷时，C(x,y,z,t)趋向于零；当时间趋向于无穷时，C(x,y,z,t)也趋于无穷。

对于问题二，要探究风速对放射性物质

放射气体模型的预估模型

摘要

本文是以日本福岛核电站遭遇自然灾害发生核泄漏的背景而提出的。且结合了高斯烟羽模型、线性拟合，以及微分方程模型，运用MATLAB软件，分析泄漏源强度、风速、大气稳定度参数等因素对放射性气体扩散的影响，预测了放射性气体浓度在不同时间，不同地区的浓度变化，并且本文模型中数据可以根据不同的实际情况而加以改变，因而是本文的应用范围大大增加，可以适用于具有较强的应用型。

对于问题一，讨论在无风的情况下，放射性气体以s m/s的匀速在大气中向四周扩散。本问中由于不考虑风力的影响，且扩散出来的气体匀速向四周散开，这样经过任意时刻t，扩散的气体围成一个半径为st的球，且距球心位置不同的地方浓度值不同。采用列数列的表现方法，设定相同时间段t，把条件进行整理，并经过简单计算得出每段时间所预测得到的扩散距离r和浓度C。利用MATLAB软件对数据进行线性拟合，采用微分方程模型得到核电站周边放射性气体在不同地区，不同时间段的浓度变化，得出随着离泄漏源距离的延伸，最后放射性物质的浓度越来越小，趋近于零，即当x趋向无穷时，C(x,y,z,t)趋向于零；当时间趋向于无穷时，C(x,y,z,t)也趋于无穷。

对于问题二，要探究风速对放射性物质

篇一：高斯烟羽模型

模型假设：

1、 坐标系

高斯模型的坐标系如图2.1所示，原点为排放点（若为高架源，原点为排放点在地面的投影），x轴正向为风速方向，y轴在水平面上垂直于x轴，正向在x轴的左侧，z轴垂直于水平面xoy，向上为正向。在此坐标系下烟流中心线或烟流中心线在xoy面的投影与x轴重合。

2、模型假设

(1)污染物的浓度在y、z轴上的分布是高斯分布（正态分布）的；

(2)污染源的源强是连续且均匀的，初始时刻云团内部的浓度、温度呈均匀分布；

(3)扩散过程中不考虑云团内部温度的变化，忽略热传递、热对流及热辐射；

(4)泄漏气体是理想气体，遵守理想气体状态方程；

(5)在水平方向，大气扩散系数呈各向同性；

(6)取x轴为平均风速方向，整个扩散过程中风速的大小、方向保持不变，不随地点、时间变化而变化；

(7)地面对泄漏气体起全反射作用，不发生吸收或吸附作用；

(8)整个过程中，泄漏气体不发生沉降、分解，不发生任何化学反应等。

3、模型公式推导

由正态分布假设可以导出下风向任意一点X（x,y,z）处泄漏气体浓度的函数为：

X(x,y,z)?A(x)e?ay2e?bz2

（1）

由概率统计理论可以写出方差的表达式为：

???2??y

???2??z???

Q??????0y2Xdy??

摘要

问题的重述

设有一座核电站遇自然灾害发生泄漏，浓度为p0的放射性气体以匀速排出，速度为m kg/s，在无风的情况下，匀速在大气中向四周扩散, 速度为s m/s. 1）请你建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。

2）当风速为k m/s时，给出核电站周边放射性物质浓度的变化情况。

3）当风速为k m/s时，分别给出上风和下风L公里处，放射性物质浓度的预测模型。 4）将你建立的模型应用于福岛核电站的泄漏，计算出福岛核电站的泄漏对我国东海岸，及美国西海岸的影响。

计算所用数据可以在网上搜索或根据具体情况自己模拟。

模型的假设

1. 放射性气体的传播服从菲克第一定律，即单位时间通过单位法向面积的流量与它的浓度成正比。

2. 假定无地形差异，不考虑逆温，对流对放射性气体浓度变化的影响。 3. 放射性气体在空气中扩散时没有损失，质量始终守恒，且在扩散过程中不发生化学反应。

4. 放射性气体由某一点源向各个方向等强度释放进行扩散，源强是连续并且均匀的。

主要基本符号说明

符号 说明 下风向距离 横风向距离 距地面高度 放射性气体在空间中的浓度 放射性污染物源强，即释放率 风速 分别为水平方向和垂直方向的扩

工业煤气基于MATLAB仿真的泄漏扩散影响研究

摘要：本文的研究目的是研究企业范围空间煤气泄漏的扩散规律和影响范围。采用matlab模拟煤气泄漏后CO 的浓度分布和扩散距离规律。通过建立煤气泄漏扩散数学模型, 对其影响煤气扩散的主要因素进行了分析、探讨了煤气毒性范围的划分, 然后在对煤气泄漏造成的危害和泄漏原因的基础上，运用扩散模型，计算煤气泄露扩散影响范围，然后用MATLAB对此进行模拟，得出不同的距离下煤气的浓度，并对其进行分析。因为大气稳定度、风速对煤气泄漏扩散的浓度影响起着非常重要的作用。大气稳定度和风速会显著改变有害气体的扩散状态。在风速和泄漏增大时, 煤气在开放空间扩散距离大, 影响范围广, 应合理布置煤气监控点, 预防煤气中毒。本文还鉴于煤气泄漏的危害之大，根据CO 的特性,对于煤气柜这种重大危险源的管理和控制可以得出一些经验，为采取措施预防其危害提供一定的依据。

关键词： 煤气泄漏； MATLAB；数值模拟； 扩散

一、前言

煤气泄漏的研究的背景及意义

我国当代工业以煤炭为主要能源的结构特点,决定了我国大多数工业企业的生产性气源以焦炉煤气和高炉煤气等为主,而煤气具有易燃易爆性!易散发性!剧毒性的特点,随着煤气在石油!化工!冶金

关于核电站泄漏放射性气体扩散的预估模型

摘要

由于核泄漏导致放射性气体扩散对经济和人身造成巨大损失的报道在国内外屡见不鲜，本文中日本福岛核泄漏事件更加使我们认识到对放射性气体扩散进行合理性的预估从而为以后类似于此的突发性事件作积极有效的补救措施的重要性。

对于问题一我们运用了点源烟羽扩散模型，用抛物型二阶偏微分方程解出理想状态下的不同时刻、不同地点的浓度表达式：C(x,y,z,t)?Qe(4?kt)32?x2?y2?z24kt。

此模型是建立在以泄漏点为圆心的一个无界球形区域内的。为了使模型更符合实际情况，能够被应用于现实生活中，我们在泄漏源有效高度的确定和考虑地面反射与吸收作用下对此模型进行了修正，最终得到问题一浓度的确定公式(14)C(x,y,z,t)的表达式。

对于问题二，我们采用高位连续点源烟羽扩散模式，其扩散服从正态分布，

并根据概率论的相关知识通过数学公式推导，得到理想状态下的高斯模型，由泄漏源有效高度，地面反射等因素的影响对其进行修正，又由于重力干沉积，雨洗湿沉积以及核衰变等因素对源强的影响，对高斯烟羽模型再次进行修正，最终得到泄漏源周边浓度变化情况即公式(32)，在风速为km/s的条件下浓度为

C(x,y,z,H)。

对于问

载流子的扩散运动

2010-03-09 09:19:17| 分类： 微电子物理 | 标签： |字号大中小 订阅

(什么是扩散系数？什么是扩散长度？) 作者：Xie M. X. (UESTC，成都市)

扩散是粒子在混乱热运动（布朗运动）基础之上的、在浓度梯度驱动之下的一种定向运动。半导体中载流子的扩散与原子的扩散不同；少数载流子与多数载流子的扩散也不相同。

（1）载流子扩散与原子扩散的区别： 半导体中载流子的扩散与原子的扩散，都是依靠浓度梯度所产生的一种定向运动；扩散流密度与浓度梯度成正比，其比例系数就是扩散系数（表征着扩散的快慢，单位是cm2/s）。但是载流子扩散与原子扩散的机理不同。

载流子扩散是在不断遭受散射的情况下所产生的定向运动；而原子扩散是在晶体热缺陷的帮助下所产生的定向运动。一般来说，温度越高，载流子遭受晶格振动的散射就越厉害，则扩散越慢；只有在低温下，晶格振动散射不大时，扩散才是随着温度的升高而加快（因为载流子的动能增加所致）。

而对于原子的扩散，温度越高，晶体热缺陷就越多（有指数关系），则扩散就越快；在低温下，几乎不产生热缺陷，则扩散也就慢得几乎无法进行。在热扩散技术中，就是通过高温来进行掺杂的（掺入施主或者受主杂质

7.2.2扩散系数

费克定律中的扩散系数Ｄ代表单位浓度梯度下的扩散通量，它表达某个组分在介质中扩散的快慢，是物质的一种传递性质。 一、气体中的扩散系数

气体中的扩散系数与系统、温度和压力有关，其量级为10m/s。通常对于二元气体Ａ、B的相互扩散，Ａ在B中的扩散系数和B在A中的扩散系数相等，因此可略去下标而用同一符号Ｄ表示，即DAB?DBA?D。

表７－１给出了某些二元气体在常压下（1.013?10Pa）的扩散系数。

对于二元气体扩散系数的估算，通常用较简单的由富勒（Fuller）等提出的公式：

5?520.0101T1.75D?11?MAMB （７－１９）

P[(?vA)1/3?(?vB)1/3]22式中，D－Ａ、B二元气体的扩散系数，m/s；

P－气体的总压，Pa； T－气体的温度，Ｋ；

MA、MB－组分Ａ、B的摩尔质量，kg/kmol；

、－组分Ａ、B分子扩散体积，cm/mol。

一般有机化合物可按分子式由表７－２查相应的原子扩散体积加和得到，某些简单物质则在表７－２种直接列出。 系统 H2－空气 He－空气 O2－空气 Cl2－空气 H2O－空气 NH3－空气 CO2－空气 SO2－空气 5表7－１ 某些二元气

扩散模拟软件 DICTRA

DICTRA 是一个专门用于模拟多元合金扩散型相变的通用软件包。 DICTRA 软件建立在对多组元扩散方程数值求解的基础上，其与 Thermo-Calc 软件配合使用， Thermo-Calc 用来处理必要的热力学计算。扩散模拟是基于储存于数据库内的已经评估的热 / 动力学数据数据。利用必要的热 / 动力学数据， DICTRA 可同时模拟包含 10 种元素以上合金的扩散过程。在开发 DICTRA 过程中，重点放在把基本算法与关键的热 / 动力学数据联系起来，从而使计算模拟更能接近于真实条件而发挥模拟的重要作用。 扩散模拟可应用于一维空间和三种不同的几何形状：圆柱，球形和平面， 利用这些条件可以模拟许多您感兴趣的扩散现象。例如，圆柱几何模型可以用于模拟管道中的扩散以及杆状析出物的溶解。

灵活方便的 、 带有强大专用模块的用户界面 DICTRA 开发了许多专业模块来完成 一些特殊的 功能，如数据读取、计算、后处理等。 用户界面包括 一个 可键入命 令的命令行界面 ，其中 所有命令都可采用缩写形式键入。使用宏文件可以保存复杂的 输入 命令序列， 该宏文件可以修 改、保存并 重复 使用 。 应用实例 ：

合金均匀化

Centrotherm扩散炉培训

培训内容

开机操作 生产操作 关机操作 维护保养 生产提示

开机操作(一)1、打开气柜压缩空 气的手动阀并调节 调压阀使其压力在 6Bar左右。 2、打开气柜氮气和 氧气手动阀，并调 节调压阀使其压力 在2Bar左右。

开机操作(二)3、打开电气柜上的 “Main Switch”在 “ON”位置，设备 上电；CMI电脑自 动启动。 4、打开设备冷却水， 先开出水阀，再开 进水阀。检查水压 在3-6kg /cm2 范围 内。

开机操作(三)5、输入XP操作系统用户“Operator”密码“6065”, 进入CMI计算机的XP系统，自动启动CMI软件。 6、检查并打开气柜上1-4管“CESAR”的开关，启 动“CESAR”系统。

开机操作(四)7、 设定各管“Paddle”温度为默认值，按下： “Manual/T_Paddle/All/Setpoint/Default” 。 8、设定升温速度为5℃/min,按下： “Manual/T_Paddle/All/RampPos”,输入“5”,并按确认键 “ENTER”。

开机操作(五)9、当温度到达500℃ 时，用手推动炉门 的气缸，使炉门关 闭。 10、顺时针打开TGA 上部调压