CN112507631B - 一种窄通道流动失稳出口界限含气率限值测试方法及系统 - Google Patents

一种窄通道流动失稳出口界限含气率限值测试方法及系统 Download PDF

Info

Publication number
CN112507631B
CN112507631B CN202011301106.7A CN202011301106A CN112507631B CN 112507631 B CN112507631 B CN 112507631B CN 202011301106 A CN202011301106 A CN 202011301106A CN 112507631 B CN112507631 B CN 112507631B
Authority
CN
China
Prior art keywords
limit
dynamic self
condition
feedback system
outlet
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011301106.7A
Other languages
English (en)
Other versions
CN112507631A (zh
Inventor
黄彦平
王云
袁德文
刘文兴
岳倪娜
王艳林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nuclear Power Institute of China
Original Assignee
Nuclear Power Institute of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nuclear Power Institute of China filed Critical Nuclear Power Institute of China
Priority to CN202011301106.7A priority Critical patent/CN112507631B/zh
Publication of CN112507631A publication Critical patent/CN112507631A/zh
Application granted granted Critical
Publication of CN112507631B publication Critical patent/CN112507631B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/08—Probabilistic or stochastic CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Abstract

本发明公开了一种窄通道流动失稳出口界限含气率限值测试方法及系统，本发明首先分别获得竖直静止及在其他参数保持不变的情况下的典型运动条件下动态自反馈系统中并联窄通道流动失稳界限含气率；其次在此基础上在参数范围内获得竖直静止和典型运动条件下的并联窄通道流动失稳出口界限含气率限值的影响因子；最后利用条件概率分析方法获得耦合运动条件下动态自反馈条件下窄通道流动失稳出口界限含气率限值。本发明利用条件概率分析方法科学合理的确定耦合运动条件下动态自反馈条件下窄通道流动失稳出口界限含气率限值，解决了现行限值确定方法中的热工裕量保守性和安全评价的问题。

Description

一种窄通道流动失稳出口界限含气率限值测试方法及系统
技术领域
本发明属于热工力学技术领域，具体涉及一种耦合运动条件下动态自反馈系统中的并联窄通道流动失稳出口界限含气率限值测试方法及系统。
背景技术
在换热系统的热工水力设计过程中，不发生通道间的流道脉动(并联通道流动失稳)是重要的热工水力设计准则，其流动失稳的限值参数是相应换热系统热工水力设计的重要限值之一。这种限值的确定大多采用实验的方法获得的运行参数及附近工况参数范围内最极端条件下的限值，并在此基础上取一定的安全裕量作为最终设计准则限值，对于具有动态自反馈特性的换热系统而言，其设计准则限值的确定也沿用了上述方法。
然而，由于运动条件的复杂性使得上述方法在获取运动条件下的设计准则限值的时候难以适用。其主要原因在于在动态自反馈条件下，耦合的运动条件的作用难以实现与竖直静止时的点对点对应。在现行的运动条件下具有动态自反馈特性的窄通道流动失稳出口界限含气率限值往往在竖直静止的数据基础上采用过大的安全系数来保证设计值的安全，这种方式很大程度上掩藏了换热系统的热工裕量，不利于其热工性能的提升；同时，在运动条件下具有动态自反馈特性的系统中的并联通道流动失稳特性机制尚不清晰的情况下，以安全系数来保证设计的安全的科学性有待进一步论证，可靠性差，这也影响了相关换热系统的安全性评价。
发明内容
为了克服现有的窄通道流动失稳出口界限含气率限值测试方法存在掩藏了换热系统的热工裕量，不利于其热工性能的提升，同时并没有考虑耦合运动对其的影响，导致可靠性差的技术问题，本发明提供了解决上述问题的一种窄通道流动失稳出口界限含气率限值测试方法。
本发明通过下述技术方案实现：
一种窄通道流动失稳出口界限含气率限值测试方法，该方法包括以下步骤：
步骤S1，获取竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率；并获取典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率；且所述典型运动条件和所述竖直静止条件的其他参数条件相同；
步骤S2，基于获得的出口界限含气率得到单一典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；
步骤S3，基于单一典型运动条件下的影响因子得到耦合运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；
步骤S4，基于耦合运动条件下的影响因子优化得到耦合运动条件下动态自反馈系统中并联窄通道流道流动失稳出口界限含气率限值，用于优化换热系统的热工性能。
本发明是在获得竖直和典型运动条件下具有动态自反馈特性的系统中并联窄通道流动失稳界限含气率的基础上，分析计算得到综合作用下的耦合运动条件对该出口界限含气率的影响因子，从而优化得到运动条件下具有动态自反馈特性的系统中并联窄通道流动失稳出口界限含气率限值，优化了换热系统的热工性能，提升了相关换热系统的热工安全性。
优选的，本发明的步骤S2通过下式得到单一典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子：
Figure BDA0002786865390000021
式中，f(xi)表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，xc，i表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，xc，0表示竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，xi表示第i种典型运动条件。
优选的，本发明的步骤S3得到耦合运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子与单一典型运动条件下的影响因子的关系如下：
f(x1x2...xn)＝f(x1)[1-f(x2...xn)]+f(x2...xi)[1-f(x1)]
f(x2...xn)＝f(x1)f(x2)[1-f(x3...xn)]+f(x3...xi)[1-f(x1)f(x2)]
…
f(xn-1xn)＝f(x1)f(x2)...f(xn-1)[1-f(xn)]+f(xn)[1-f(x1)f(x2)...f(xn-1)]
式中，f(x1x2...xn)表示耦合n种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，f(xi)表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，i＝1,2,…,n。
优选的，本发明的步骤S4通过下式优化得到耦合运动条件下动态自反馈系统中并联窄通道流道流动失稳出口界限含气率限值为：
xc，mov＝xc，0f(x1x2...xn)
式中，xc，0表示竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，f(x1x2...xn)表示耦合n种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子。
优选的，本发明的步骤S4中，为保证换热系统的安全性，在优化得到耦合运动条件下的出口界限含气率限值xc，mov的过程中，竖直静止条件在参数范围内取最小出口界限含气率xc，min，各典型运动条件下的影响因子取最小值f(xi，min)；
其中，xc，min＝min xc，0，f(xi，min)＝min f(xi)。
另一方面，本发明还提出了一种窄通道流动失稳出口界限含气率限值测试系统，该系统包括数据获取模块、单一影响因子模块、耦合影响因子模块和优化模块；
所述数据获取模块用于获取竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率；并获取典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率；且所述典型运动条件和所述竖直静止条件的其他参数条件相同；
所述单一影响因子模块基于获得的出口界限含气率得到单一典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；
所述耦合影响因子模块基于单一典型运动条件下的影响因子得到耦合运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；
所述优化模块基于耦合运动条件下的影响因子优化得到耦合运动条件下动态自反馈系统中并联窄通道流道流动失稳出口界限含气率限值，用于优化换热系统的热工性能。
优选的，本发明的单一影响因子模块通过下式得到单一典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子：
Figure BDA0002786865390000041
式中，f(xi)表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，xc，i表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，xc，0表示竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，xi表示第i种典型运动条件。
优选的，本发明的耦合影响因子模块得到耦合运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子与单一典型运动条件下的影响因子的关系如下：
f(x1x2...xn)＝f(x1)[1-f(x2...xn)]+f(x2...xi)[1-f(x1)]
f(x2...xn)＝f(x1)f(x2)[1-f(x3...xn)]+f(x3...xi)[1-f(x1)f(x2)]
...
f(xn-1xn)＝f(x1)f(x2)...f(xn-1)[1-f(xn)]+f(xn)[1-f(x1)f(x2)...f(xn-1)]
式中，f(x1x2...xn)表示耦合n种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，f(xi)表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，i＝1,2,…,n。
优选的，本发明的优化模块通过下式优化得到耦合运动条件下动态自反馈系统中并联窄通道流道流动失稳出口界限含气率限值为：
xc，mov＝xc，0f(x1x2...xn)
式中，xc，0表示竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，f(x1x2...xn)表示耦合n种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子。
优选的，本发明为保证换热系统的安全性，所述优化模块在优化得到耦合运动条件下的出口界限含气率限值xc，mov的过程中，竖直静止条件在参数范围内取最小出口界限含气率xc，min，各典型运动条件下的影响因子取最小值f(xi，min)；
其中，xc，min＝minxc，0，f(xi，min)＝minf(xi)。
本发明具有如下的优点和有益效果：
1、本发明在分别获得竖直静止及典型运动条件在动态自反馈条件下对并联窄通道流动失稳出口界限含气率限值的基础上，利用条件概率分析方法分析确定耦合运动条件下动态自反馈条件下窄通道流动失稳出口界限含气率限值，解决了现行限值确定方法中的热工裕量保守性和安全评价的问题。
2、本发明可用为具有动态自反馈特性的换热系统，在运动条件下的热工水力设计与安全评价提供了重要的数据支撑，优化了换热系统的热工性能，提升了相关换热系统的热工安全性。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：
图1为本发明的方法流程示意图。
图2为本发明的系统原理框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
实施例1
相较于现有的测试方法，本实施例提出了一种窄通道流动失稳出口界限含气率限值测试方法，本实施例的方法在竖直静止及对应运动条件在动态自反馈条件下对并联窄通道流动失稳出口界限含气率限值的基础上，利用条件概率分析方法优化得到耦合运动条件下动态自反馈条件下窄通道流动失稳出口界限含气率限值。
具体如图1所示，本实施例的方法包括以下步骤：
1、获取竖直静止条件动态自反馈条件下并联窄通道流动失稳出口界限含气率
调节其他参数至既定值，台阶式提升并联通道的加热功率，构建冷热段的密度差驱动回路中的冷却工质运转，持续台阶式提升并联通道的加热功率直到流动失稳发生，并持续一段时间记录相关参数，获得竖直静止条件动态自反馈条件下并联窄通道流动失稳出口界限含气率。
2、获取典型运动条件动态自反馈条件下并联窄通道流动失稳出口界限含气率
在保持其他条件不变(即与1中相同)的情况下，启动运动平台，开展典型运动条件下的以并联通道为加热段的动态自反馈条件下的流动失稳特性研究。台阶式提升并联通道的加热功率，构建冷热段的密度差驱动回路中的冷却工质运转，持续台阶式提升并联通道的加热功率直到流动失稳发生，并持续一段时间记录相关参数，获得典型运动条件动态自反馈条件下并联窄通道流动失稳出口界限含气率。
3、获得典型运动条件对并联窄通道出口界限含气率的单因素影响因子。以xi表示第i种典型运动条件，则相应的影响因子f(xi)为第i种典型运动条件动态自反馈条件下窄通道流动失稳出口界限含气率xc，i和竖直静止条件下的动态自反馈条件下窄通道流动失稳出口界限含气率xc，0的比值。
Figure BDA0002786865390000071
其中，i＝1,2,…,n，其中n为表示典型运动条件数量。
4、采用条件概率分析方法得到耦合运动条件动态自反馈条件下窄通道流动失稳出口界限含气率限值影响因子，与单一典型运动条件的影响因子关系如下所示：
f(x1x2...xn)＝f(x1)[1-f(x2...xn)]+f(x2...xi)[1-f(x1)]
f(x2...xn)＝f(x1)f(x2)[1-f(x3...xn)]+f(x3...xi)[1-f(x1)f(x2)]
…
f(xn-1xn)＝
f(x1)f(x2)...f(xn-1)[1-f(xn)]+f(xn)[1-f(x1)f(x2)...f(xn-1)]
式中，f(x1x2...xn)表示耦合n种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子。
5、耦合运动条件动态自反馈条件下窄通道流动失稳出口界限含气率限值。通过下式优化获得的耦合运动条件对动态自反馈条件下窄通道流动失稳出口界限含气率限值为：
xc，mov＝xc，0f(x1x2...xn)
本实施例中，在优化得到耦合运动条件下的出口界限含气率限值xc，mov的过程中，竖直静止条件在参数范围内取最小出口界限含气率xc，min，各典型运动条件下的影响因子取最小值f(xi，min)；
其中，xc，min＝minxc，0，f(xi，mim)＝minf(xi)；
同时，为了兼顾性能优化，在确保安全性的前提下尽可能释放安全裕量。
实施例2
基于上述实施例1，本实施例还提出了一种窄通道流动失稳出口界限含气率限值测试系统，具体如图2所示，本实施例的系统包括数据获取模块、单一影响因子模块、耦合影响因子模块和优化模块。
本实施例的数据获取模块用于获取竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率；并获取典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率；且所述典型运动条件和所述竖直静止条件的其他参数条件相同；即本实施例的数据获取模块用于执行上述实施例1提出的方法步骤1和步骤2。
本实施例的单一影响因子模块基于获得的出口界限含气率得到单一典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；即本实施例的单一影响因子模块用于执行上述实施例1提出的方法步骤3。
本实施例的耦合影响因子模块基于单一典型运动条件下的影响因子得到耦合运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；即本实施例的耦合影响因子模块用于执行上述实施例1提出的方法步骤4。
本实施例的优化模块基于耦合运动条件下的影响因子优化得到耦合运动条件下动态自反馈系统中并联窄通道流道流动失稳出口界限含气率限值，用于优化换热系统的热工性能；即本实施例的优化模块用于执行上述实施例1提出的方法步骤5。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种窄通道流动失稳出口界限含气率限值测试方法，其特征在于，该方法包括以下步骤：
步骤S1，获取竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率；并获取典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率；且所述典型运动条件和所述竖直静止条件的其他参数条件相同；
步骤S2，基于获得的出口界限含气率得到单一典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；
步骤S3，基于单一典型运动条件下的影响因子得到耦合运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；
步骤S4，基于耦合运动条件下的影响因子优化得到耦合运动条件下动态自反馈系统中并联窄通道流道流动失稳出口界限含气率限值，用于优化换热系统的热工性能；
所述步骤S4通过下式优化得到耦合运动条件下动态自反馈系统中并联窄通道流道流动失稳出口界限含气率限值为：
xc，mov＝xc，0f(x1x2...xn)
式中，xc，0表示竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，f(x1x2...xn)表示耦合n种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；
所述步骤S4中，为保证换热系统的安全性，在优化得到耦合运动条件下的出口界限含气率限值xc，mov的过程中，竖直静止条件在参数范围内取最小出口界限含气率xc，min，各典型运动条件下的影响因子取最小值f(xi，min)；
其中，xc，min＝min xc，0，f(xi，min)＝min f(xi)；
f(xi)表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，i＝1，2，…，n。
2.根据权利要求1所述的一种窄通道流动失稳出口界限含气率限值测试方法，其特征在于，所述步骤S2通过下式得到单一典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子：
Figure FDA0003509985710000021
式中，f(xi)表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，xc，i表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，xc，0表示竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，xi表示第i种典型运动条件。
3.根据权利要求1所述的一种窄通道流动失稳出口界限含气率限值测试方法，其特征在于，所述步骤S3得到耦合运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子与单一典型运动条件下的影响因子的关系如下：
f(x1x2...xn)＝f(x1)[1-f(x2...xn)]+f(x2...xi)[1-f(x1)]
f(x2...xn)＝f(x1)f(x2)[1-f(x3…xn)]+f(x3…xi)[1-f(x1)f(x2)]
…
f(xn-1xn)＝f(x1)f(x2)...f(xn-1)[1-f(xn)]+f(xn)[1-f(x1)f(x2)...f(xn-1)]
式中，f(x1x2...xn)表示耦合n种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，f(xi)表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，i＝1，2，…，n。
4.一种窄通道流动失稳出口界限含气率限值测试系统，其特征在于，该系统包括数据获取模块、单一影响因子模块、耦合影响因子模块和优化模块；
所述数据获取模块用于获取竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率；并获取典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率；且所述典型运动条件和所述竖直静止条件的其他参数条件相同；
所述单一影响因子模块基于获得的出口界限含气率得到单一典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；
所述耦合影响因子模块基于单一典型运动条件下的影响因子得到耦合运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；
所述优化模块基于耦合运动条件下的影响因子优化得到耦合运动条件下动态自反馈系统中并联窄通道流道流动失稳出口界限含气率限值，用于优化换热系统的热工性能；
所述优化模块通过下式优化得到耦合运动条件下动态自反馈系统中并联窄通道流道流动失稳出口界限含气率限值为：
xc，mov＝xc，0f(x1x2...xn)
式中，xc，0表示竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，f(x1x2...xn)表示耦合n种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子；
为保证换热系统的安全性，所述优化模块在优化得到耦合运动条件下的出口界限含气率限值xc，mov的过程中，竖直静止条件在参数范围内取最小出口界限含气率xc，min，各典型运动条件下的影响因子取最小值f(xi，min)；
其中，xc，min＝min xc，0，f(xi，min)＝min f(xi)；
f(xi)表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，i＝1，2，…，n。
5.根据权利要求4所述的一种窄通道流动失稳出口界限含气率限值测试系统，其特征在于，所述单一影响因子模块通过下式得到单一典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子：
Figure FDA0003509985710000041
式中，f(xi)表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，xc，i表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，xc，0表示竖直静止条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率，xi表示第i种典型运动条件。
6.根据权利要求4所述的一种窄通道流动失稳出口界限含气率限值测试系统，其特征在于，所述耦合影响因子模块得到耦合运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子与单一典型运动条件下的影响因子的关系如下：
f(x1x2...xn)＝f(x1)[1-f(x2...xn)]+f(x2...xi)[1-f(x1)]
f(x2...xn)＝f(x1)f(x2)[1-f(x3…xn)]+f(x3…xi)[1-f(x1)f(x2)]
…
f(xn-1xn)＝f(x1)f(x2)...f(xn-1)[1-f(xn)]+f(xn)[1-f(x1)f(x2)...f(xn-1)]
式中，f(x1x2...xn)表示耦合n种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，f(xi)表示第i种典型运动条件下动态自反馈系统中并联窄通道流动失稳出口界限含气率限值影响因子，i＝1，2，…，n。
CN202011301106.7A 2020-11-19 2020-11-19 一种窄通道流动失稳出口界限含气率限值测试方法及系统 Active CN112507631B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011301106.7A CN112507631B (zh) 2020-11-19 2020-11-19 一种窄通道流动失稳出口界限含气率限值测试方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011301106.7A CN112507631B (zh) 2020-11-19 2020-11-19 一种窄通道流动失稳出口界限含气率限值测试方法及系统

Publications (2)

Publication Number Publication Date
CN112507631A CN112507631A (zh) 2021-03-16
CN112507631B true CN112507631B (zh) 2022-04-15

Family

ID=74959914

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011301106.7A Active CN112507631B (zh) 2020-11-19 2020-11-19 一种窄通道流动失稳出口界限含气率限值测试方法及系统

Country Status (1)

Country Link
CN (1) CN112507631B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113551870B (zh) * 2021-06-23 2022-09-20 中国核动力研究设计院 一种对并联通道流动失稳行为影响机制的表征方法及系统
CN116502470B (zh) * 2023-06-25 2023-09-12 中国核动力研究设计院 过冷度限值的确定方法、装置、设备、介质和程序产品

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110310710A (zh) * 2019-06-21 2019-10-08 西安交通大学 获得燃料元件矩形窄缝通道内干涸型临界热流密度的方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO330911B1 (no) * 2008-12-12 2011-08-15 Multi Phase Meters As Fremgangsmåte og apparat for måling av sammensetning og strømningsrater for en våtgass
CN102928026B (zh) * 2012-10-22 2015-04-08 哈尔滨工程大学 一种利用局部瞬时空泡份额获得整体瞬时空泡份额的方法
CN106291701B (zh) * 2016-10-13 2019-03-15 成都理工大学 储层检测方法及装置
CN111553022B (zh) * 2020-04-30 2020-12-08 中国核动力研究设计院 动态自反馈条件下堆芯dnbr限值修正系数的求解方法
CN111507016B (zh) * 2020-04-30 2020-12-15 中国核动力研究设计院 一种运动条件动态下并联窄通道流动失稳边界的确定方法
CN111859514B (zh) * 2020-07-30 2024-01-30 西安建筑科技大学 一种多工况运行下围护结构热工性能的优化方法及系统

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110310710A (zh) * 2019-06-21 2019-10-08 西安交通大学 获得燃料元件矩形窄缝通道内干涸型临界热流密度的方法

Also Published As

Publication number Publication date
CN112507631A (zh) 2021-03-16

Similar Documents

Publication Publication Date Title
CN112507631B (zh) 一种窄通道流动失稳出口界限含气率限值测试方法及系统
Feng et al. H/sub/spl infin//controller synthesis of fuzzy dynamic systems based on piecewise Lyapunov functions and bilinear matrix inequalities
Chen et al. Composite fast‐slow MPC design for nonlinear singularly perturbed systems
CN113294885B (zh) 空调器控制方法、装置、空调器及计算机可读存储介质
CN106899519B (zh) 支付渠道流量配置方法及装置
CN104792079A (zh) 一种蒸汽型溴化锂吸收式冷水机组的温度闭环控制方法
CN112149714A (zh) 基于数据挖掘确定煤电机组能效特征指标基准值的方法
CN109162856A (zh) 水轮发电机组的状态监控方法和系统
CN109114805A (zh) 一种设备能耗的确定方法及装置
CN113504765A (zh) 一种电厂循环流化床锅炉操作优化方法
CN112231909A (zh) 一种基于无机盐相图的分盐数据拟合方法
CN112000205A (zh) Bmc及其散热策略适配方法、装置、服务器及存储介质
CN111143967A (zh) 一种改进的电网限额决策方法及系统
CN115290218A (zh) 一种火力发电机组锅炉水冷壁壁温的软测量方法与系统
CN116028193A (zh) 一种混部集群的大数据任务动态高能效调度方法和系统
CN114970766A (zh) 一种基于线性拟合的电站机组运行参数基准值获取方法
CN112613158B (zh) 一种严重事故下安全壳内控制体热力学响应综合分析方法
CN115096033A (zh) 基于体温变化的制冷策略生成方法、系统、终端及介质
US10353987B2 (en) Identification of a breakpoint based on a correlation measurement
CN114066212A (zh) 一种基于历史工况的机组升降负荷工况寻优方法
CN114398235A (zh) 基于融合学习和假设检验的内存回收趋势预警装置及方法
CN114707435B (zh) 再生式换热器不同运行方式下流量转折点分析方法及系统
Shang et al. Covariance eigenpairs neighbour distance for fault detection in chemical processes
CN113568386B (zh) 一种基于区间估计的火电机组全工况大数据分析方法
CN117435887A (zh) 基于指令分析的新能源场站受限时段精准辨识系统及方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant