天然气管道运行模拟及仿真技术研究.doc

天然气管道运行模拟及仿真技术研究 1011202045 蔡永军 科学计算选讲结课论文 为了预测天然气管道运行状态，制定合理的管输计划，更好的配置设备开机，天然气管道输送过程中需要进行工况模拟及仿真。实际工作中需要建立压缩机、阀门等设备的模型，确定管段的控制方程、气体的状态方程，针对给出的初始条件和边界条件，筛选确定天然气管网数学模型的离散方法与非线性方程组的求解算法寻找合理的非线性方程的求解算法，得到合理的数值解。 1天然气管道仿真数学模型 1.1管段的控制方程 对于管道中的任意管段，经过适当的简化可以用下列公式来描述： 连续性方程： （1） 运动方程： （2） 能量方程： （3） 式中：A——管道的横截面积，m2； ρ——流体密度，kg/m3； t——时间，s； x——坐标，m； u——速度，m/s； P——压力，Pa； θ——管道倾角，rad； λ——水力摩阻系数； D——管道内径，m； T——流体温度，k； k1——流体至管壁的换热系数； h——比焓； Tw——管壁的温度，k。 1.2 阀门控制方程 阀门控制方程如下： （4） 式中： Mup——阀门入口质量流量，kg/s； Mdw——阀门入口质量流量，kg/s Cg——阀门系数； Pup——阀的入口压力，Pa； Pdw——阀的出口压力，Pa。 1.3压缩机控制方程 简化后的压缩机控制方程如下 （5） 式中：——压缩机压比； m——多变压缩指数； n——压缩机的实际转速，rpm； n0——压缩机的额定转速，rpm； a1, b1, c1——系数； Q——给定状态下的体积流量，m3/s； 1.4 理想调节阀阀控制方程 理想调节阀控制方程如下： （6） 2气体的状态方程 采用BWRS气体状态方程，如下： （7） 式中：P——系统压力，KPa； T——系统温度，K； ρ——混合气体密度，Kmol/m3； R——气体常数，8.3143KJ/( Kmol.K)。 A0、 B0、C0、D0、E0、a、b、c、d、α、γ为方程的是一个参数，根据（8）确定。 （8） 式中：xi、xj——混合气体中i和j组分的摩尔分数； kij——为i、j组分间的交互作用系数。 3气体的焓方程 气体的焓方程如下： （9） 4 管道周边的热力模型 管道的有效土壤厚度采用等效圆筒法，传热半径由下式计算： （10） 式中：R2-R1——土壤厚度，m； R1——从管道中心至土壤层的半径，m； H——至管道中心的实际埋深，m； D——管道直径，m。 管道和周围环境的瞬态热力模型计算式如下： （11） 式中：k——周围环境导热系数； r——传热半径； Tr——r处的气体温度； Cp——气体定压比热； Tt——t时刻的气体温度。 单位管长热流量由下式表示。通过该公式计算管壁在任意节点的温度。 （12） 式中：k2——管壁至土壤换热系数； K1——流体至管壁换热系数； Tw——管壁温度； T0——R2处的温度； T——气体温度。 5 水力摩阻系数计算式 管段控制方程涉及的水力摩阻系数λ采用F.Colebrook－White公式计算，该公式表达如下： （13） 式中，e/D——管道粗糙度和内径的无因次比； Re——雷诺数。 6控制方程的离散化 由管道控制方程与气体状态方程组成的非线性偏微分方程组，一般不能得出管流气体基本变量的解析解，因此有必要应用计算数学的方法求解偏微分方程组的数值解。本专题中选用中心隐式差分法对控制方程进行离散化。确定采用的基本变量为气体的密度（ρ）、速度（u）和温度（T）。 6.1离散形式 引进变量φ，φ代表三个流动基本中的任意一个。在时间步长为Δx , 空间步长为Δt 的情况下，以空间i和时间网格点t采用中心隐式差分格式，则有以下离散形式： 对于基本流动变量： （14） 基本流动变量对时间的一阶偏导数： （15） 基本流动变量对空间的一阶偏导数： （16） 基本流动变量对时间的二阶偏导数： （17） 基本流动变量对空间的二阶偏导数： （18） 基本流动变量对空间及时间的二阶偏导数： （19） 6.2 离散后的控制方程 离散后的控制方程如下： 离散后的连续性方程： （20） 离散后的运动方程： （21） 离散后的能量方程： （22） 6.3 初始条件与边界条件 初始条件指