使用VB.net实现的天然气物理性质计算BWRS方程

Public Class Gas_Hysical_Property
    Private Pressure As Single '定义压力,MPa
    Private Temperature As Single '定义温度,℃

    Private Relative_Molecular_Mass(17) As Single '纯组分相对分子质量
    Private Percentage_Gas_Temperament(17) As Single '输入的气体组分
    Private Total_Relative_Molecular_Mass As Single '计算得到的气体分子量

    Private Temperature_Critical(17) As Single '临界温度,K
    Private Density_Critical(17) As Single '临界密度,kmol/m3
    Private Acentric_Factor(17) As Single '偏心因子
    Private Interaction_Coefficient(17, 17) As Single 'BWRS交互作用系数

    Private Const_A(10) As Single '通用常数Ai
    Private Const_B(10) As Single '通用常数Bi

    Private Const R = 8.3143 '气体常数，kJ/（kmol*K）
    ''''纯组分状态方程的11个参数,以及混合气体的11个参数''''
    Private B0_Pure(17), B0_Mix As Single
    Private A0_Pure(17), A0_Mix As Single
    Private C0_Pure(17), C0_Mix As Single
    Private Gama_Pure(17), Gama_Mix As Single
    Private b_Pure(17), b_Mix As Single
    Private a_Pure(17), a_Mix As Single
    Private Alpha_Pure(17), Alpha_Mix As Single
    Private c_Pure(17), c_Mix As Single
    Private D0_Pure(17), D0_Mix As Single
    Private d_Pure(17), d_Mix As Single
    Private E0_Pure(17), E0_Mix As Single
    ''''纯组分状态方程的11个参数,以及混合气体的11个参数''''

    Private Density_mol As Single '定义摩尔密度,kmol/m3
    Private Density_mol_Temp(200) As Single '求密度时的的临时变量
    Private Density As Single '定义密度，kg/m3
    Private Z As Single '定义压缩因子
    Private Partial_Differential_Coefficient_Temperture As Single '偏微分系数，偏p比偏T
    Private Partial_Differential_Coefficient_Density As Single '偏微分系数，偏p比偏T
    Private SHCPC_B(17) As Single, SHCPC_C(17) As Single, SHCPC_D(17) As Single, SHCPC_E(17) As Single, SHCPC_F(17) As Single
    Private Cp_Ideal(17) As Single
    Private Cp0_Mix As Single, Cv0_Mix As Single '理想气体纯组分的比热
    Private Cp_Mix As Single, Cv_Mix As Single '理想气体混合物的比热
    Private Di As Single 'J-T系数（Di），℃/MPa
    Private Adiabatic_Exponent_Temperature As Single '定温绝热指数
    Private Adiabatic_Exponent_Volume As Single '定容绝热指数
    Private LHV_Pure(17) As Single '纯组分的低热值
    Private LHV_Mix As Single '天然气混合物的低热值
    Private Relative_Density As Single '相对密度
    Private X_Temp As Single : Private Y_Temp As Single : Private C_Temp As Single
    Private Dynamic_Viscosity As Single '天然气的动力粘度(mPa*s)


    Private Sub Gas_Hysical_Property_Load(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles MyBase.Load

    End Sub

    Private Sub Button1_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Button1.Click
        Call ClearData()
        Pressure = Val(TextBox1.Text) * 1000 '将压力单位转化为kPa
        Temperature = Val(TextBox2.Text) + 273.1 '将温度转化为K
        LHV_Pure = {35.9, 59.477, 64.4, 87.667, 93.24, 122.85, 123.65, 155.72, 156.73, 156.73, 156.73, 156.73, 156.73, 156.73, 156.73, 0, 0, 23.36} '各组分的低热值