买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

龙图腾网恭喜哈尔滨工程大学申请的专利一种管道瞬变流的计算方法获国家发明授权专利权，本发明授权专利权由国家知识产权局授予，授权公告号为：CN115758926B 。

龙图腾网通过国家知识产权局官网在2025-05-06发布的发明授权授权公告中获悉：该发明授权的专利申请号/专利号为：202211385031.4，技术领域涉及：G06F30/28；该发明授权一种管道瞬变流的计算方法是由田兆斐;张乾平设计研发完成，并于2022-11-07向国家知识产权局提交的专利申请。

本一种管道瞬变流的计算方法在说明书摘要公布了：本发明公开一种管道瞬变流的计算方法，包括：步骤S1：设计基于非恒定流阻力模型的高效波追踪方程组；步骤S2、确定短管划分标准，进行管道节点划分以及时间步长的确定；步骤S3：给定初始条件和边界条件，进行稳态计算；步骤S4：具体进行高效波追踪法的计算。本发明技术方案在进行管道节点划分中，短管划分比例系数和波速调整系数会影响计算精度和计算效率，两个系数越小，计算精度越高但计算效率将会降低。本发明中采用高效波追踪法能够解决该问题，因为高效波追踪法的计算时间对管道分段数不敏感，所以可以尽可能的满足精度需求而选取小的短管划分比例系数和波速调整系数。

本发明授权一种管道瞬变流的计算方法在权利要求书中公布了：1.一种管道瞬变流的计算方法，其特征在于，所述方法包括：步骤S1、设计基于非恒定流阻力模型的高效波追踪方程组，所述高效波追踪方程组的表达式具体包括： 其中，符号fA表示所述管道的端点A处的下游波，fB表示所述管道的端点B处的下游波，符号LAB表示管道的端点A和所述管道的端点B间的距离；符号g表示重力加速度；符号J表示所述管道的当前摩擦损失，符号J0表示所述管道的初始摩擦损失；下角标A、B表示单根管道的上下游端点，J-J0A表示所述端点A到所述端点B的压力波损失，J-J0B表示端点B到端点A的压力波损失，式1和式2中的各项可由下式计算得到： 其中，符号f表示初始流速正方向的下游波；符号F表示初始流速负方向的上游波；符号H表示静压头，符号Q表示体积流量；符号a表示波速；符号D表示所述管道的直径；符号Ac表示所述管道的横截面积；符号λ表示所述管道的沿程阻力系数，符号C*表示取值与流态有关的Brunone阻力系数；符号Q0表示管道稳态时的流速；符号QA表示所述端点A处的流速，符号QB所述端点B处的流速，符号λA表示所述端点A处的沿程阻力系数，符号λB表示所述端点B处的沿程阻力系数，QA,0所述端点A处的初始稳态流速，QB,0所述端点B处的初始稳态流速，符号λA,0表示所述端点A处的初始稳态沿程阻力系数，符号λB,0表示所述端点B处的初始稳态沿程阻力系数；步骤S2、确定短管划分标准，进行管道节点划分以及时间步长的确定；步骤S3：给定初始条件和边界条件，进行稳态计算；步骤S4：具体进行高效波追踪法的计算，具体包括：步骤S41、每一根管道都建立一个基于空间域的压力波库，压力波库由两个分别存储管道上游波F和下游波f的数组组成，数组大小为N+1；下游波数组中的元素0,1…N与管道空间节点位置一一对应，上游波数组元素与管道节点位置的对应关系相反，数组元素存储管道相应位置处的压力波幅值；压力波N表示传入管道端点的压力波，压力波0表示管道端点传出的压力波；步骤S42、在每一个时间步长下，调用传播到边界节点处的压力波N，进行边界节点方程的求解，由边界处产生的压力波与入射波的方向是相反的，上游边界节点A处进入的压力波为上游压力波库F中的压力波N，产生的压力波为下游压力波库中的压力波0，新产生的压力波0在插入到数组前先进行从节点A到节点B的压力波衰减计算，新产生的压力波记为F0，f0，由式6和式7来进行衰减计算： 对于新产生的压力波0直接从数组元素0处插入数组中，使数组中的元素朝着压力波传播的方向前进一个位置；数组元素N被调用后直接删除；步骤S43、在计算过程中追踪压力波f和F，通过下式获取边界节点的压头和体积流量，更新管道端点的瞬时水力状态，在下一个时间步长下返回步骤S2，直至完成瞬态计算；

如需购买、转让、实施、许可或投资类似专利技术，可联系本专利的申请人或专利权人哈尔滨工程大学，其通讯地址为：150000 黑龙江省哈尔滨市南岗区南通大街145号；或者联系龙图腾网官方客服，联系龙图腾网可拨打电话0551-65771310或微信搜索“龙图腾网”。