高斯扩散模型公式

篇一：高斯烟羽模型

模型假设：

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??

工业煤气基于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)。

对于问

放射气体模型的预估模型

摘要

本文是以日本福岛核电站遭遇自然灾害发生核泄漏的背景而提出的。且结合了高斯烟羽模型、线性拟合，以及微分方程模型，运用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)也趋于无穷。

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

高斯投影坐标正反算

一、基本思想：

高斯投影正算公式就是由大地坐标（L，B）求解高斯平面坐标（x，y），而高斯投影反算公式则是由高斯平面坐标（x，y）求解大地坐标（L，B）。

二、计算模型：

基本椭球参数： 椭球长半轴a 椭球扁率f

椭球短半轴：b?a(1?f)

a2?b2椭球第一偏心率 ：e?

aa2?b2椭球第二偏心率 ：e?? b高斯投影正算公式：此公式换算的精度为0.001m

x?X??NN2??sinBcosB?l?simBcos3B(5?t2?9?2?4?4)l??4242???24???

N5246??sinBcosB(61?58t?t)l720???6y?NN3223??cosB?l???cosB(1?t??)l???6???3N5242225???cosB(5?18t?t?14??58?t)l120???5

其中：角度都为弧度

B为点的纬度，l???L?L0，L为点的经度，L0为中央子午线经度；

N为子午圈曲率半径，N?a(1?esinB)； t?tanB；

22?12?2?e?2cos2B

????180??3600

其中X为子午线弧长：

1616??X?a0B?sinBcosB?(a2?a4?a6)?(2a4?a6)sin

高斯投影坐标正反算

一、相关概念

大地坐标系由大地基准面和地图投影确定，由地图投影到特定椭圆柱面后在南北两极剪开展开而成，是对地球表面的逼近，各国或地区有各自的大地基准面，我国目前主要采用的基准面为：

1.WGS84基准面，为GPS基准面，17届国际大地测量协会上推荐，椭圆柱长半轴a=6378137m，短半轴b=6356752.3142451m；

2.西安80坐标系，1975年国际大地测量协会上推荐，椭圆柱长半轴a=6378140m，短半轴b=6356755.2881575m;

3.北京54坐标系，参照前苏联克拉索夫斯基椭球体建立,椭圆柱长半轴a=6378245m, 短半轴b=6356863.018773m;

通常所说的高斯投影有三种，即投影后：

a) 角度不变（正角投影），投影后经线和纬线仍然垂直； b) 长度不变； c) 面积不变；

大地坐标一般采用高斯正角投影，即在地球球心放一点光源，地图投影到过与中央经线相切的椭圆柱面上而成；可分带投影，按中央经线经度值分带，有每6度一带或每3度一带两种(起始带中央经线经度为均为3度，即：6度带1带位置0-6度，3度带1带位置1.5-4.5 度)，即所谓的高斯-克吕格投影。

图表 11高斯投影和分带

§8.3高斯投影坐标正反算公式

任何一种投影①坐标对应关系是最主要的；②如果是正形投影，除了满足正形投影的条件外（C-R偏微分方程），还有它本身的特殊条件。 8.3.1高斯投影坐标正算公式： B,l x,y

高斯投影必须满足以下三个条件：

①中央子午线投影后为直线；②中央子午线投影后长度不变；③投影具有正形性质，即正形投影条件。

由第一条件知中央子午线东西两侧的投影必然对称于中央子午线，即(8-10)式中，x为l的偶函数，y为l的奇函数；l 30

30 ，即l /

1/20，如展开为l的级数，收敛。

x m24

6

0 m2l m4l m6l y m3

5

（8-33）

1l m3l m5l

式中

m0,m1, 是待定系数，它们都是纬度B的函数。

由第三个条件知：

xx y

q

y l, l q

(8-33)式分别对l和q求偏导数并代入上式

m1 3m24

dm0dm22

dm44

3l 5m5l

dq

dql dql 2ml3

6m5

6l

dm1 (8-34)

dq

l

dm355

2l 4m4dq

l3

dmdq

l

上两式两边相等，其必要充分条件是同次幂l前的系数应相等，即

mdm01

dq

m 1dm1

22

dq

(8-35)

m1dm2

3 3

§8.3高斯投影坐标正反算公式

任何一种投影①坐标对应关系是最主要的；②如果是正形投影，除了满足正形投影的条件外（C-R偏微分方程），还有它本身的特殊条件。 8.3.1高斯投影坐标正算公式： B,l x,y

高斯投影必须满足以下三个条件：

①中央子午线投影后为直线；②中央子午线投影后长度不变；③投影具有正形性质，即正形投影条件。

由第一条件知中央子午线东西两侧的投影必然对称于中央子午线，即(8-10)式中，x为l的偶函数，y为l的奇函数；l 3030 ，即l / 1/20，如展开为l的级数，收敛。

x m0 m2l2 m4l4 m6l6 y m1l m3l3 m5l5

（8-33）

式中m0,m1, 是待定系数，它们都是纬度B的函数。 由第三个条件知：

x y x y

, q l l q

(8-33)式分别对l和q求偏导数并代入上式

dm0dm22dm44

m1 3m3l 5m5l l l

dqdqdq

(8-34) dm33dm55dm135

2m2l 4m4l 6m6l l l l

dqdqdq

2

4

上两式两边相等，其必要充分条件是同次幂l前的系数应相等，即

61

dm0

m1

dq

1dm1

m2

2dq

1dm2

Green公式、Stokes公式、Gauss公式在专业学科中

的应用

摘要

格林（Green）公式，斯托克斯（Stokes）公式和高斯（Gauss）公式是多元函数积分学中的三个基本公式，它们分别建立了曲线积分与二重积分、曲面积分与三重积分、曲线积分和曲面积分的联系。它们建立了向量的散度与通量、旋度与环量之间的关系，除了在数学上应用于计算多元函数积分，在其他领域也有很多重要的应用。本文将主要从这三个公式与物理学之间的联系展开介绍它们的其他应用，其中包括应用于GPS面积测量仪，确定外部扰动重力场，应用于保守场以及推证阿基米德定律和高斯定理等，帮助人们加深对格林公式、斯托克斯公式和高斯公式的理解，从而能够更准确地应用此三个公式。

关键词：格林公式斯托克斯公式高斯公式散度旋度应用

目录

一、引言 ......................................... 1 二、格林（Green）公式的应用 ...................... 1

（一）格林公式的定义 .............................. 1 1、单连通区域的概念 ..................