一种制冷剂气液两相可压缩流体充注质量流量解算方法技术

本发明专利技术公开了一种制冷剂气液两相可压缩流体充注质量流量解算方法。本发明专利技术通过制冷剂在变截面管道中入口段、收缩段及喉道处的压强值，根据连续性方程和等熵流方程，建立实际质量流量方程，根据伯努利方程和连续性方程得到入口段制冷剂的密度和入口段制冷剂流速之间的关系，利用训练好的PINN模型，得到入口段制冷剂密度，从而解算出制冷剂的实际质量流量；本发明专利技术仅使用压力传感器并结合物理神经网络即可解算制冷剂充注过程中的质量流量，能够非接触、低成本、快速准确地控制制冷剂充注的质量，能够实现制冷剂的高效利用和绿色节能减排。排。排。

【技术实现步骤摘要】
一种制冷剂气液两相可压缩流体充注质量流量解算方法


[0001]本专利技术涉及冷剂充注流量检测技术，具体涉及一种制冷剂气液两相可压缩流体充注质量流量解算方法。

技术介绍

[0002]气液两相流普遍存在于石油、化工、电力、冶金等领域，与其相关的各种参数的测量也是不可或缺的。比如，在锅炉、热交换器、冷凝器、油气运输管路等设备的状态监控、安装管理等的过程中，都需要对其中涉及到的气液两相流参数(如流量、流型、相含率等)进行可靠的测量。然而，由于气液两相流固有的复杂性，两相流参数的检测难度非常大，尤其是流量的在线不分离的测量。
[0003]其中制冷剂是一种在制冷系统中不断循环的工作物质，通过各种热机改变自身的状态(气液态转换)来实现能量转换，因此将制冷剂看作一种可压缩气液两相流体。制冷设备都有其制冷剂充注量的标准，过多或过少的充注都会影响设备的制冷效果，因此实现制冷剂气液混合物的精确充注以及制冷过程中临界流量参数的准确测量对提高制冷效率、降低排放和成本起着至关重要的作用。尤其是医药领域中的贮藏及运输中转仓储对超低温制冷系统的苛刻要求，因此急需一种实现制冷剂的准确充注的方法。由于制冷剂在填充和流动过程中的挥发性，其充注过程流量检测非常困难。

技术实现思路

[0004]针对以上现有技术存在的问题，本专利技术提出了一种制冷剂气液两相可压缩流体充注质量流量解算方法。
[0005]本专利技术的制冷剂气液两相可压缩流体充注质量流量解算方法，包括以下步骤：
[0006]1)装置连接：
[0007]a)变截面管道包括入口段、收缩段、喉道、扩散段和出口段，其中，入口段的截面积为A1，喉道的截面积为A2，收缩段的锥角为冷媒罐通过导管连接至变截面管道的入口段；变截面管道的出口段通过导管连接至待充注制冷外机；
[0008]b)在入口段、收缩段和喉道的内壁分别安装第一、第二和第三压力传感器，第一、第二和第三压力传感器分别经过集成处理电路连接至上位机；
[0009]c)打开冷媒罐，使制冷剂通过变截面管道，通过第一、第二和第三压力传感器分别测得入口段、收缩段和喉道处的压强P1、P2和P3；
[0010]2)根据连续性方程和等熵流方程，建立实际质量流量方程：
[0011]a)在变截面管道中，制冷剂充注过程宏观上看作无粘等熵可压缩流体流动过程，制冷剂充注过程中无摩擦效应，满足以下质量连续性方程：
[0012][0013]其中，g为重力加速度，为质量流量，G0为理想情况下质量流量，γ1＝ρ1g为入口段制冷剂比重，ρ1为入口段制冷剂密度，γ2为喉道处制冷剂比重，V1为入口段制冷剂流速，
V2为喉道制冷剂流速；
[0014]b)制冷剂充注过程中，通过变截面管道的制冷剂满足以下等熵流方程为：
[0015][0016]其中，k为制冷剂的比热比；
[0017]c)将(1)式和(2)式结合，得到理想情况下质量流量G0为：
[0018][0019]d)引入变截面管道的流量系数C和制冷剂的膨胀系数Y，得到实际质量流量G的方程为：
[0020][0021]其中，G为制冷剂的实际质量流量，D1和D2分别为入口段和喉道的截面圆的直径，膨胀系数Y表达为：
[0022][0023]其中，k为制冷剂的比热比，
[0024]3)根据伯努利方程及连续性方程，计算得到入口段制冷剂的流速和密度之间的关系：
[0025]入口段和收缩段的压力采集处满足以下伯努利方程：
[0026][0027]解得：
[0028][0029]从而通过入口段的压强P1和收缩段的压强P3推导出入口段制冷剂流速V1与入口段制冷剂密度ρ1之间的关系；
[0030]4)利用物理神经网络PINN模型得到入口段制冷剂的密度ρ1：
[0031]a)设定物理神经网络PINN模型的物理约束：
[0032]对于制冷剂气液混合流的瞬态一维流动，通过将质量守恒定律和动量定律应用于变截面管道获得以下质量守恒方程和动量守恒方程：
[0033][0034][0035]其中，ρ为制冷剂气液混合流密度，A为变截面管道横截面积，V为速度，P为压强，λ为摩擦系数，t为时间，x为沿变截面管道流动的距离；
[0036]制冷剂气液混合流微观上看作由含有气泡的不可压缩液体组成，气泡均匀分布并
根据多变过程可得到：
[0037][0038]其中，ρ
g
为气相密度，ρ
g0
为初始热力学条件下的气相密度，P0为初始热力学条件下的压强，n为多变指数；
[0039]根据气体质量比定义制冷剂气液混合流密度ρ：
[0040][0041]其中，气体质量比M
g
和M
l
分别为气相和液相流体的分子量，ρ
l
为液相密度；
[0042]声速c表示为：
[0043][0044]质量守恒方程(8)引入声速c，得到引入声速的质量守恒方程为：
[0045][0046]引入声速的质量守恒方程(13)和动量守恒方程(9)统称为偏微分方程；
[0047]b)构建物理神经网络PINN模型：
[0048]构建多层前馈神经网络f(x,t；θ)作为物理神经网络PINN模型，其中θ为构建多层前馈神经网络过程中可调参数权重和偏差的集合，物理神经网络输入x为沿变截面管道流动的距离，t为时间，V(x,t；θ)、P(x,t；θ)和ρ(x,t；θ)分别为时空坐标(x,t)下所对应的制冷剂速度V、压强P和密度ρ的物理神经网络的输出的估计值，其中
[0049]V、P和ρ都是关于(x,t)的函数，都能够用时空坐标(x,t)以及中间的各项参数解析地表达出来；
[0050]多层前馈神经网络f(x,t；θ)用以逼近两个偏微分方程映射：
[0051][0052]c)构建总损失函数：
[0053]定义物理约束残差项：
[0054][0055][0056][0057]其中，r1(x,t；θ)、r2(x,t；θ)和r3(x,t；θ)分别表示引入声速的质量守恒方程(13)的残差项、动量守恒方程(9)的残差项和方程(7)的残差项，统称为物理约束残差项，作为构造物理损失函数中的算子；
[0058]利用物理约束残差项构造物理损失函数
[0059][0060]其中，为时空坐标下选取残差训练点的集合，x为选取残差点坐标(x,t)；
[0061]从而通过偏微分方程(13)和(9)及方程(7)的残差项得到的物理损失函数作为物理约束；
[0062]将入口段的压强P1作为先验物理数据，并作为数据约束放到损失函数中去作为函数拟合中的约束，构造数据损失函数
[0063][0064]其中，为时空坐标下先验物理数据点的集合，P(x,t；θ)为时空坐标下物理神经网络的输出的估计值，P1(x,t)为时空坐标对应的入口段的压强P1的真实值即先验物理数据，从而以通过先验物理本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种制冷剂气液两相可压缩流体充注质量流量解算方法，其特征在于，所述方法包括以下步骤：1)装置连接：a)变截面管道包括入口段、收缩段、喉道、扩散段和出口段，其中，入口段的截面积为A1，喉道的截面积为A2，收缩段的锥角为冷媒罐通过导管连接至变截面管道的入口段；变截面管道的出口段通过导管连接至待充注制冷外机；b)在入口段、收缩段和喉道的内壁分别安装第一、第二和第三压力传感器，第一、第二和第三压力传感器分别经过集成处理电路连接至上位机；c)打开冷媒罐，使制冷剂通过变截面管道，通过第一、第二和第三压力传感器分别测得入口段、收缩段和喉道处的压强P1、P2和P3；2)根据连续性方程和等熵流方程，建立实际质量流量方程：a)在变截面管道中，制冷剂充注过程宏观上看作无粘等熵可压缩流体流动过程，制冷剂充注过程中无摩擦效应，满足以下质量连续性方程：其中，g为重力加速度，为质量流量，G0为理想情况下质量流量，γ1＝ρ1g为入口段制冷剂比重，ρ1为入口段制冷剂密度，γ2为喉道处制冷剂比重，V1为入口段制冷剂流速，V2为喉道制冷剂流速；b)制冷剂充注过程中，通过变截面管道的制冷剂满足以下等熵流方程为：其中，k为制冷剂的比热比；c)将(1)式和(2)式结合，得到理想情况下质量流量G0为：d)引入变截面管道的流量系数C和制冷剂的膨胀系数Y，得到实际质量流量G的方程为：其中，G为制冷剂的实际质量流量，D1和D2分别为入口段和喉道的截面圆的直径，膨胀系数Y表达为：其中，k为制冷剂的比热比，3)根据伯努利方程及连续性方程，计算得到入口段制冷剂的流速和密度之间的关系：入口段和收缩段的压力采集处满足以下伯努利方程：解得：
从而通过入口段的压强P1和收缩段的压强P3推导出入口段制冷剂流速V1与入口段制冷剂密度ρ1之间的关系；4)利用物理神经网络PINN模型得到入口段制冷剂的密度ρ1：a)设定物理神经网络PINN模型的物理约束：对于制冷剂气液混合流的瞬态一维流动，通过将质量守恒定律和动量定律应用于变截面管道获得以下质量守恒方程和动量守恒方程：面管道获得以下质量守恒方程和动量守恒方程：其中，ρ为制冷剂气液混合流密度，A为变截面管道横截面积，V为速度，P为压强，λ为摩擦系数，t为时间，x为沿变截面管道流动的距离；制冷剂气液混合流微观上看作由含有气泡的不可压缩液体组成，气泡均匀分布并根据多变过程可得到：其中，ρ
g
为气相密度，ρ
g0
为初始热力学条件下的气相密度，P0为初始热力学条件下的压强，n为多变指数；根据气体质量比定义制冷剂气液混合流密度ρ：其中，气体质量比M
g
和M
l
分别为气相和液相流体的分子量，ρ
l
为液相密度；声速c表示为：质量守恒方程(8)引入声速c，得到引入声速的质量守恒方程为：引入声速的质量守恒方程(13)和动量守恒方程(9)统称为偏微分方程；b)构建物理神经网络PINN模型：构建多层前馈神经网络f(x,t；θ)作为物理神经网络PINN模型，其中θ为构建多层前馈神经网络过程中可调参数权重和偏差的集合，物理神经网络输入x为沿变截面管道流动的距离，t为时间，V(x,t；θ)、P(x,t；θ)和ρ(x,t；θ)分别为时空坐标(x,t)下...

【专利技术属性】
技术研发人员：娄文忠，宣炜琨，张正谦，丁男希，付胜华，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：