CN111982213A - 一种核反应堆模拟燃料组件流量测量方法及系统 - Google Patents

一种核反应堆模拟燃料组件流量测量方法及系统 Download PDF

Info

Publication number
CN111982213A
CN111982213A CN202010843390.4A CN202010843390A CN111982213A CN 111982213 A CN111982213 A CN 111982213A CN 202010843390 A CN202010843390 A CN 202010843390A CN 111982213 A CN111982213 A CN 111982213A
Authority
CN
China
Prior art keywords
reynolds number
section
function
flow
venturi tube
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.)
Granted
Application number
CN202010843390.4A
Other languages
English (en)
Other versions
CN111982213B (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 CN202010843390.4A priority Critical patent/CN111982213B/zh
Publication of CN111982213A publication Critical patent/CN111982213A/zh
Application granted granted Critical
Publication of CN111982213B publication Critical patent/CN111982213B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/05—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
    • G01F1/34—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure
    • G01F1/36—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure the pressure or differential pressure being created by the use of flow constriction
    • G01F1/40—Details of construction of the flow constriction devices
    • G01F1/44—Venturi tubes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/05—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects
    • G01F1/34—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure
    • G01F1/36—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure the pressure or differential pressure being created by the use of flow constriction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F25/00—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume
    • G01F25/10—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume of flowmeters

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明公开了一种核反应堆模拟燃料组件流量测量方法。本发明还公开了一种核反应堆模拟燃料组件流量测量系统。本发明在不变更经典文丘里管几何结构的前提下，拓宽经典文丘里管的量程比，使经典文丘里管的可测最小流量进一步降低。基于拓宽后的Re范围内对应的流量系数α存在特殊的变化趋势，本发明拟提出一种宽量程经典文丘里管的流量计算方法，利用该计算方法无需迭代计算，通过直接求解方程的方法，实现高效高精度的宽量程文丘里管流量测量。

Description

一种核反应堆模拟燃料组件流量测量方法及系统
技术领域
本发明涉及核技术领域，具体涉及一种核反应堆模拟燃料组件流量测量方法及系统。
背景技术
《GB/T 2624-2006用安装在圆形截面管道中的差压装置测量满管流体流量》中的第四部分对经典文丘里管的使用介绍中，存在约束条件，要求雷诺数Re≥2×105，真实使用条件下，往往需要文丘里流量计有较高量程比，能够测量更小的流量，需要在不变更经典文丘里流量计几何结构前提下，拓宽流量计标定试验的流量范围，这种情况下雷诺数Re会小于2×105，在国标GB/T 2624-2006的附录B部分《超出GB/T 2624.4范围使用的经典文丘里管》有这样的描述，机械加工经典文丘里管当雷诺数Re减小到2×105以下时，往往会出现流量系数α值在随着Re的减小而稳定地减小之前略有增大。
而在核领域中，反应堆整体水力模拟试验对模拟燃料组件流量进行测量时，由于反应堆运行存在多种工况，不同工况下雷诺数变化幅度很大，这就使得通用技术中的流量测量方式不适用于核反应堆模拟燃料组件流量的测量，并且无相关文献可参考。
发明内容
本发明所要解决的技术问题是在核领域中，反应堆整体水力模拟试验对模拟燃料组件流量进行测量时，由于反应堆运行存在多种工况，不同工况下雷诺数变化幅度很大，这就使得通用技术中的流量测量方式不适用于核反应堆模拟燃料组件流量的测量，并且无相关文献可参考，目的在于提供一种核反应堆模拟燃料组件流量测量方法及系统，解决上述问题。
本发明通过下述技术方案实现：
一种核反应堆模拟燃料组件流量测量方法，包括以下步骤：
S1：使用文丘里管对模拟燃料组件中流体的流量系数α和雷诺数Re的关系进行标定生成标定关系函数；所述标定关系函数的自变量为雷诺数Re，应变量为流量系数α；
S2：对所述标定关系函数中上升段的进行拟合生成上升段拟合函数；对所述标定关系函数中下降段的进行拟合生成下降段拟合函数；
S3：将所述上升段拟合函数和文丘里管的理论关系函数联立生成上升段模型；将所述下降段拟合函数和文丘里管的理论关系函数联立生成下降段模型；所述文丘里管的理论关系函数根据文丘里管中雷诺数Re和流量系数α的理论关系获取；
S4：对模拟燃料组件中流体的流量进行测量时，将从文丘里管和测温仪表获取的数据输入所述下降段模型获取下降段输出雷诺数Re′；
当下降段输出雷诺数Re′位于所述标定关系函数的下降段内时，将所述下降段输出雷诺数Re′作为实际雷诺数输出；
当下降段输出雷诺数Re′未位于所述标定关系函数的下降段内时，将从文丘里管和测温仪表获取的数据输入所述上升段模型获取上升段输出雷诺数Re″；
如果所述上升段输出雷诺数Re″位于所述标定关系函数的上升段内时，将所述上升段输出雷诺数Re″作为实际雷诺数输出；
如果所述上升段输出雷诺数Re″未位于所述标定关系函数的上升段内时，扩大流量系数α和雷诺数Re的关系的标定范围并执行S1～S4；
S5：根据所述实际雷诺数获取模拟燃料组件中流体的流量。
现有技术中，核反应堆的运行包括大量工况，同时还包括有各种事故工况；在不同工况下，模拟燃料组件内的流量变化很大，尤其是在一些流量较小的工况下，单个模拟燃料组件内的流量会很低，普通测量流量的方式，如经典文丘里管测流量的方式会出现量程不足的问题。
本发明应用时，首先对经典文丘里管在一个较宽的流量范围开展实流标定试验，生成标定关系函数，而在进行标定过程中，由于在低雷诺数时，流量系数α和雷诺数之间关系会发生特殊变化，所以发明人对标定关系函数进行分段处理，发明人发现在一个宽量程的文丘里管中，标定关系函数会出现先上升，后下降的趋势；这里所谓的上升段是指斜率大于零的段，而下降段是指斜率小于零的段，为了便于后续判断，将上升段和下降段的交点并入下降段处理。
然后对上升段和下降段分别进行拟合，由于标定关系函数本身的特殊性，如果后续对标定关系函数进行直接求解，只能选择通过迭代手段求取数值解，运算速度慢并且求解精度低，而在本发明中采用分别拟合的方式，非常有利于后续求解，拟合后的多项式可以直接求出解析解，也便于后续进行结果判断。
然后生成上升段模型和下降段模型，这种模型是基于文丘里管和其他相关仪器采集的数据进行设计的，采用将拟合函数和规范中的文丘里管的理论关系函数进行联立获得，文丘里管的理论关系函数主要包括有文丘里管体积流量公式，这样就可以将文丘里管和其他相关仪器采集的数据直接带入模型获取对应的雷诺数，其他相关仪器主要为测温仪表。
在使用上升段模型和下降段模型进行流量测量时，需要进行分段判断，使用下降段模型进行流量判断时，需要最终得到的雷诺数落在下降段模型的雷诺数范围内，而使用上升段模型行流量判断时，需要最终得到的雷诺数落在上升段模型的雷诺数范围内，本身两个模型都是可以直接获取解析解，所以进行雷诺数范围的判断会非常便捷容易；如果无法满足上述条件，说明测量已经超量程了，此时需要重新进行标定并循环整个流程，这样本发明就具有非常强的适应性。本发明通过上述步骤，具备通用性，针对单相流体流量测量，只需在经典文丘里管实流标定阶段，适当拓宽流量标定范围，利用该流量计算方法，可高精度拟合雷诺数Re与流量系数α的函数关系；在后续实际流量测量阶段，可利用该雷诺数Re与流量系数α的函数关系式无需迭代计算，快速并高精度的求解得到实际流量，实现经典文丘里管的宽量程流量测量。
进一步的，步骤S2中：
所述上升段拟合函数为α＝DRe2+ERe+F，Re∈[Re1、Re2)；
所述下降段拟合函数为α＝ARe2+BRe+C，Re∈[Re2、Re3]；
式中α为流量系数，Re为雷诺数，A、B、C、D、E和F均为常数，Re1为上升段的雷诺数起点，Re2为上升段和下降段交点处的雷诺数，Re3为下降段的雷诺数终点。
进一步的，步骤S3中：
文丘里管的理论关系函数为A1*α-B1*Re＝0；其中
Figure BDA0002642230610000031
Figure BDA0002642230610000032
上式中，ρ为入口管取压点处流体密度，根据流体介质温度和压力确定；ΔP为文丘里管差压；d为文丘里管喉部直径；ε为膨胀系数，根据直径比、压力比和等熵指数确定；υ为入口管流体运动粘度，根据流体介质温度和压力确定；D为文丘里管入口管直径。
进一步的，步骤S3中：
下降段模型为A2Re2+B2Re+C2＝0；其中A2＝A1*A，B2＝A1*B-B1，C2＝A1*C；
上升段模型为A3Re2+B3Re+C3＝0；其中A3＝A1*D，B3＝A1*E-B1，C3＝A1*F；
进一步的，步骤S4包括以下子步骤：
下降段输出雷诺数Re′包括两个解：Re1′和Re2′；其中：
Figure BDA0002642230610000033
当Re1′和Re2′中任意一个位于[Re2、Re3]内时，将位于[Re2、Re3]内的下降段输出雷诺数Re′作为实际雷诺数输出。
进一步的，步骤S4包括以下子步骤：
上升段输出雷诺数Re″包括两个解：Re3′和Re4′；其中：
Figure BDA0002642230610000041
当Re3′和Re4′中任意一个位于[Re1、Re2)内时，将位于[Re1、Re2)内的上升段输出雷诺数Re″作为实际雷诺数输出。
进一步的，步骤S5包括以下子步骤：
根据下式获取模拟燃料组件中流体的流量：
Q＝B1*Re”'；
式中Re”'为实际雷诺数，Q为模拟燃料组件中流体的流量。
采用上述任意一项所述方法的一种核反应堆模拟燃料组件流量测量系统，包括：
标定单元：用于使用文丘里管对模拟燃料组件中流体的流量系数α和雷诺数Re的关系进行标定生成标定关系函数；所述标定关系函数的自变量为雷诺数Re，应变量为流量系数α；
拟合单元：用于对所述标定关系函数中上升段的进行拟合生成上升段拟合函数；对所述标定关系函数中下降段的进行拟合生成下降段拟合函数；
模型单元：用于将所述上升段拟合函数和文丘里管的理论关系函数联立生成上升段模型；将所述下降段拟合函数和文丘里管的理论关系函数联立生成下降段模型；所述文丘里管的理论关系函数根据文丘里管中雷诺数Re和流量系数α的理论关系获取；
测量单元：用于对模拟燃料组件中流体的流量进行测量时，将从文丘里管和测温仪表获取的数据输入所述下降段模型获取下降段输出雷诺数Re′；
当下降段输出雷诺数Re′位于所述标定关系函数的下降段内时，将所述下降段输出雷诺数Re′作为实际雷诺数输出；
当下降段输出雷诺数Re′未位于所述标定关系函数的下降段内时，将从文丘里管和测温仪表获取的数据输入所述上升段模型获取上升段输出雷诺数Re″；
如果所述上升段输出雷诺数Re″位于所述标定关系函数的上升段内时，将所述上升段输出雷诺数Re″作为实际雷诺数输出；
如果所述上升段输出雷诺数Re″未位于所述标定关系函数的上升段内时，扩大流量系数α和雷诺数Re的关系的标定范围，且标定单元、拟合单元、模型单元和测量单元循环工作；
所述测量单元还用于根据所述实际雷诺数获取模拟燃料组件中流体的流量。
进一步的，所述上升段拟合函数为α＝DRe2+ERe+F，Re∈[Re1、Re2)；
所述下降段拟合函数为α＝ARe2+BRe+C，Re∈[Re2、Re3]；
式中α为流量系数，Re为雷诺数，A、B、C、D、E和F均为常数，Re1为上升段的雷诺数起点，Re2为上升段和下降段交点处的雷诺数，Re3为下降段的雷诺数终点。
进一步的，文丘里管的理论关系函数为A1*α-B1*Re＝0；其中
Figure BDA0002642230610000051
上式中，ρ为入口管取压点处流体密度，根据流体介质温度和压力确定；ΔP为文丘里管差压；d为文丘里管喉部直径；ε为膨胀系数，根据直径比、压力比和等熵指数确定；υ为入口管流体运动粘度，根据流体介质温度和压力确定；D为文丘里管入口管直径。
本发明与现有技术相比，具有如下的优点和有益效果：
本发明一种核反应堆模拟燃料组件流量测量方法及系统，具备通用性，针对单相流体流量测量，只需在经典文丘里管实流标定阶段，适当拓宽流量标定范围，利用该流量计算方法，可高精度拟合雷诺数Re与流量系数α的函数关系；在后续实际流量测量阶段，可利用该雷诺数Re与流量系数α的函数关系式无需迭代计算，快速并高精度的求解得到实际流量，实现经典文丘里管的宽量程流量测量。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：
图1为本发明方法步骤示意图；
图2为宽量程经典文丘里管雷诺数Re与流量系数α关系图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
实施例
如图1所示，本发明一种核反应堆模拟燃料组件流量测量方法，包括以下步骤：
S1：使用文丘里管对模拟燃料组件中流体的流量系数α和雷诺数Re的关系进行标定生成标定关系函数；所述标定关系函数的自变量为雷诺数Re，应变量为流量系数α；
S2：对所述标定关系函数中上升段的进行拟合生成上升段拟合函数；对所述标定关系函数中下降段的进行拟合生成下降段拟合函数；
S3：将所述上升段拟合函数和文丘里管的理论关系函数联立生成上升段模型；将所述下降段拟合函数和文丘里管的理论关系函数联立生成下降段模型；所述文丘里管的理论关系函数根据文丘里管中雷诺数Re和流量系数α的理论关系获取；
S4：对模拟燃料组件中流体的流量进行测量时，将从文丘里管和测温仪表获取的数据输入所述下降段模型获取下降段输出雷诺数Re′；
当下降段输出雷诺数Re′位于所述标定关系函数的下降段内时，将所述下降段输出雷诺数Re′作为实际雷诺数输出；
当下降段输出雷诺数Re′未位于所述标定关系函数的下降段内时，将从文丘里管和测温仪表获取的数据输入所述上升段模型获取上升段输出雷诺数Re″；
如果所述上升段输出雷诺数Re″位于所述标定关系函数的上升段内时，将所述上升段输出雷诺数Re″作为实际雷诺数输出；
如果所述上升段输出雷诺数Re″未位于所述标定关系函数的上升段内时，扩大流量系数α和雷诺数Re的关系的标定范围并执行S1～S4；
S5：根据所述实际雷诺数获取模拟燃料组件中流体的流量。
现有技术中，核反应堆的运行包括大量工况，同时还包括有各种事故工况；在不同工况下，模拟燃料组件内的流量变化很大，尤其是在一些流量较小的工况下，单个模拟燃料组件内的流量会很低，普通测量流量的方式，如经典文丘里管测流量的方式会出现量程不足的问题。
本实施例实施时，首先对经典文丘里管在一个较宽的流量范围开展实流标定试验，生成标定关系函数，而在进行标定过程中，由于在低雷诺数时，流量系数α和雷诺数之间关系会发生特殊变化，所以发明人对标定关系函数进行分段处理，发明人发现在一个宽量程的文丘里管中，标定关系函数会出现先上升，后下降的趋势；这里所谓的上升段是指斜率大于零的段，而下降段是指斜率小于零的段，为了便于后续判断，将上升段和下降段的交点并入下降段处理。
然后对上升段和下降段分别进行拟合，由于标定关系函数本身的特殊性，如果后续对标定关系函数进行直接求解，只能选择通过迭代手段求取数值解，运算速度慢并且求解精度低，而在本发明中采用分别拟合的方式，非常有利于后续求解，拟合后的多项式可以直接求出解析解，也便于后续进行结果判断。
然后生成上升段模型和下降段模型，这种模型是基于文丘里管和其他相关仪器采集的数据进行设计的，采用将拟合函数和规范中的文丘里管的理论关系函数进行联立获得，文丘里管的理论关系函数主要包括有文丘里管体积流量公式，这样就可以将文丘里管和其他相关仪器采集的数据直接带入模型获取对应的雷诺数，其他相关仪器主要为测温仪表。
在使用上升段模型和下降段模型进行流量测量时，需要进行分段判断，使用下降段模型进行流量判断时，需要最终得到的雷诺数落在下降段模型的雷诺数范围内，而使用上升段模型行流量判断时，需要最终得到的雷诺数落在上升段模型的雷诺数范围内，本身两个模型都是可以直接获取解析解，所以进行雷诺数范围的判断会非常便捷容易；如果无法满足上述条件，说明测量已经超量程了，此时需要重新进行标定并循环整个流程，这样本发明就具有非常强的适应性。本发明通过上述步骤，具备通用性，针对单相流体流量测量，只需在经典文丘里管实流标定阶段，适当拓宽流量标定范围，利用该流量计算方法，可高精度拟合雷诺数Re与流量系数α的函数关系；在后续实际流量测量阶段，可利用该雷诺数Re与流量系数α的函数关系式无需迭代计算，快速并高精度的求解得到实际流量，实现经典文丘里管的宽量程流量测量。
为了进一步的说明本实施例的工作过程，步骤S2中：
所述上升段拟合函数为α＝DRe2+ERe+F，Re∈[Re1、Re2)；
所述下降段拟合函数为α＝ARe2+BRe+C，Re∈[Re2、Re3]；
式中α为流量系数，Re为雷诺数，A、B、C、D、E和F均为常数，Re1为上升段的雷诺数起点，Re2为上升段和下降段交点处的雷诺数，Re3为下降段的雷诺数终点。
为了进一步的说明本实施例的工作过程，步骤S3中：
文丘里管的理论关系函数为A1*α-B1*Re＝0；其中
Figure BDA0002642230610000071
Figure BDA0002642230610000072
上式中，ρ为入口管取压点处流体密度，根据流体介质温度和压力确定；ΔP为文丘里管差压；d为文丘里管喉部直径；ε为膨胀系数，根据直径比、压力比和等熵指数确定；υ为入口管流体运动粘度，根据流体介质温度和压力确定；D为文丘里管入口管直径。
为了进一步的说明本实施例的工作过程，步骤S3中：
下降段模型为A2Re2+B2Re+C2＝0；其中A2＝A1*A，B2＝A1*B-B1，C2＝A1*C；
上升段模型为A3Re2+B3Re+C3＝0；其中A3＝A1*D，B3＝A1*E-B1，C3＝A1*F；
为了进一步的说明本实施例的工作过程，步骤S4包括以下子步骤：
下降段输出雷诺数Re′包括两个解：Re1′和Re2′；其中：
Figure BDA0002642230610000081
当Re1′和Re2′中任意一个位于[Re2、Re3]内时，将位于[Re2、Re3]内的下降段输出雷诺数Re′作为实际雷诺数输出。
为了进一步的说明本实施例的工作过程，步骤S4包括以下子步骤：
上升段输出雷诺数Re″包括两个解：Re3′和Re4′；其中：
Figure BDA0002642230610000082
当Re3′和Re4′中任意一个位于[Re1、Re2)内时，将位于[Re1、Re2)内的上升段输出雷诺数Re″作为实际雷诺数输出。
为了进一步的说明本实施例的工作过程，步骤S5包括以下子步骤：
根据下式获取模拟燃料组件中流体的流量：
Q＝B1*Re”'；
式中Re”'为实际雷诺数，Q为模拟燃料组件中流体的流量。
本实施例采用上述任意一项所述方法的一种核反应堆模拟燃料组件流量测量系统，包括：
标定单元：用于使用文丘里管对模拟燃料组件中流体的流量系数α和雷诺数Re的关系进行标定生成标定关系函数；所述标定关系函数的自变量为雷诺数Re，应变量为流量系数α；
拟合单元：用于对所述标定关系函数中上升段的进行拟合生成上升段拟合函数；对所述标定关系函数中下降段的进行拟合生成下降段拟合函数；
模型单元：用于将所述上升段拟合函数和文丘里管的理论关系函数联立生成上升段模型；将所述下降段拟合函数和文丘里管的理论关系函数联立生成下降段模型；所述文丘里管的理论关系函数根据文丘里管中雷诺数Re和流量系数α的理论关系获取；
测量单元：用于对模拟燃料组件中流体的流量进行测量时，将从文丘里管和测温仪表获取的数据输入所述下降段模型获取下降段输出雷诺数Re′；
当下降段输出雷诺数Re′位于所述标定关系函数的下降段内时，将所述下降段输出雷诺数Re′作为实际雷诺数输出；
当下降段输出雷诺数Re′未位于所述标定关系函数的下降段内时，将从文丘里管和测温仪表获取的数据输入所述上升段模型获取上升段输出雷诺数Re″；
如果所述上升段输出雷诺数Re″位于所述标定关系函数的上升段内时，将所述上升段输出雷诺数Re″作为实际雷诺数输出；
如果所述上升段输出雷诺数Re″未位于所述标定关系函数的上升段内时，扩大流量系数α和雷诺数Re的关系的标定范围，且标定单元、拟合单元、模型单元和测量单元循环工作；
所述测量单元还用于根据所述实际雷诺数获取模拟燃料组件中流体的流量。
为了进一步的说明本实施例的工作过程，所述上升段拟合函数为α＝DRe2+ERe+F，Re∈[Re1、Re2)；
所述下降段拟合函数为α＝ARe2+BRe+C，Re∈[Re2、Re3]；
式中α为流量系数，Re为雷诺数，A、B、C、D、E和F均为常数，Re1为上升段的雷诺数起点，Re2为上升段和下降段交点处的雷诺数，Re3为下降段的雷诺数终点。
为了进一步的说明本实施例的工作过程，文丘里管的理论关系函数为A1*α-B1*Re＝0；其中
Figure BDA0002642230610000091
上式中，ρ为入口管取压点处流体密度，根据流体介质温度和压力确定；ΔP为文丘里管差压；d为文丘里管喉部直径；ε为膨胀系数，根据直径比、压力比和等熵指数确定；υ为入口管流体运动粘度，根据流体介质温度和压力确定；D为文丘里管入口管直径。
为了进一步的说明本实施例的工作过程，本实施例具体的实现过程如下：
如图2所示，第一步：对经典文丘里管在更宽流量范围开展实流标定试验，获得文丘里管入口段雷诺数Re与流量系数α的关系图；
第二步：将图2曲线分成两段，下降段和上升段，并分别利用二次项公式进行拟合，得到两段Re与α的二次项函数关系式；
第三步：结合两段Re与α的二次项函数关系式、流量系数α计算公式、雷诺数Re计算公式，通过求解方程的根及对方程根的判断，得到文丘里管Re，进而得到文丘里管实际流量。
为了进一步的说明本实施例的工作过程，本实施例每个步骤的实现过程如下：
第一步：对经典文丘里管在更宽流量范围开展实流标定试验，获得文丘里管入口段雷诺数Re与流量系数α的关系图，见图2，图2共有10个标定流量点；
第二步：将图2曲线分成两段，下降段和上升段，利用下降段的6个标定流量点进行二次项公式拟合，下降段Re∈[Re2、Re3]拟合公式为
α＝ARe2+BRe+C (1)
利用上升段的5个标定流量点进行二次项公式拟合，上升段Re∈[Re1、Re2)拟合公式为
α＝DRe2+ERe+F (2)
其中A、B、C、D、E、F均为常数；
第三步：在文丘里管实际使用中，结合两段Re与α的二次项函数关系式、流量系数α计算公式、雷诺数Re计算公式，即可求解文丘里管流量，流量具体求解过程如下：
《GB/T 2624-2006用安装在圆形截面管道中的差压装置测量满管流体流量》中规定：
文丘里管体积流量：
Figure BDA0002642230610000101
式中：ρ为入口管取压点处流体密度，根据流体介质温度、压力确定；ΔP为文丘里管差压；d为文丘里管喉部直径；ε为膨胀系数，根据直径比、压力比和等熵指数确定；α为文丘里管流量系数，是一个仅与雷诺数Re相关的函数。
为方便起见，令
Figure BDA0002642230610000102
则Q＝A1*α (3)
定义入口管雷诺数Re的计算公式为：
Figure BDA0002642230610000103
式中：V为入口管流体平均流速，De为入口管当量直径，D为文丘里管入口管直径；υ为入口管流体运动粘度，根据流体介质温度、压力确定；
将以上公式变形可得到
Figure BDA0002642230610000111
令
Figure BDA0002642230610000112
则Q＝B1*Re(4)
将公式(3)代入公式(4)可得到：
A1*α-B1*Re＝0 (5)
假设Re∈[Re2、Re3]
将公式(1)代入公式(5)得到：
A1*(ARe2+BRe+C)-B1*Re＝0
A1*ARe2+(A1*B-B1)Re+A1*C＝0
令A2＝A1*A
B2＝A1*B-B1
C2＝A1*C
则A2Re2+B2Re+C2＝0 (6)
如果方程(6)可解，其根为：
Figure BDA0002642230610000113
和
Figure BDA0002642230610000114
判断Re1、Re2是否在[Re2、Re3]内，如果有其中一个根位于[Re2、Re3]内，则求解得到真实的Re，之后利用公式(4)即可得到文丘里管流量Q＝B1*Re。
如果方程(6)不可解没有根或求解得到的根不在[Re2、Re3]内，则开展后续计算。
再假设Re∈[Re1、Re2)
将公式(2)入公式(5)得到：
A1*(DRe2+ERe+F)-B1*Re＝0
A1*DRe2+(A1*E-B1)Re+A1*F＝0
令A3＝A1*D
B3＝A1*E-B1
C3＝A1*F
则A3Re2+B3Re+C3＝0 (7)
如果方程(7)可解，其根为：
Figure BDA0002642230610000121
和
Figure BDA0002642230610000122
判断Re3、Re4是否在[Re1、Re2)内，如果有其中一个根位于[Re1、Re2)内，则求解得到真实的Re，之后利用公式(5)即可得到文丘里管流量Q＝B1*Re。
如果方程(7)不可解没有根或求解得到的根不在[Re1、Re2)内，则表明实际流体Re不在区间[Re1、Re3]内，超出文丘里管雷诺数标定范围，无法测得实际流量。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种核反应堆模拟燃料组件流量测量方法，其特征在于，包括以下步骤：
S1：使用文丘里管对模拟燃料组件中流体的流量系数α和雷诺数Re的关系进行标定生成标定关系函数；所述标定关系函数的自变量为雷诺数Re，应变量为流量系数α；
S2：对所述标定关系函数中上升段的进行拟合生成上升段拟合函数；对所述标定关系函数中下降段的进行拟合生成下降段拟合函数；
S3：将所述上升段拟合函数和文丘里管的理论关系函数联立生成上升段模型；将所述下降段拟合函数和文丘里管的理论关系函数联立生成下降段模型；所述文丘里管的理论关系函数根据文丘里管中雷诺数Re和流量系数α的理论关系获取；
S4：对模拟燃料组件中流体的流量进行测量时，将从文丘里管和测温仪表获取的数据输入所述下降段模型获取下降段输出雷诺数Re′；
当下降段输出雷诺数Re′位于所述标定关系函数的下降段内时，将所述下降段输出雷诺数Re′作为实际雷诺数输出；
当下降段输出雷诺数Re′未位于所述标定关系函数的下降段内时，将从文丘里管和测温仪表获取的数据输入所述上升段模型获取上升段输出雷诺数Re″；
如果所述上升段输出雷诺数Re″位于所述标定关系函数的上升段内时，将所述上升段输出雷诺数Re″作为实际雷诺数输出；
如果所述上升段输出雷诺数Re″未位于所述标定关系函数的上升段内时，扩大流量系数α和雷诺数Re的关系的标定范围并执行S1～S4；
S5：根据所述实际雷诺数获取模拟燃料组件中流体的流量。
2.根据权利要求1所述的一种核反应堆模拟燃料组件流量测量方法，其特征在于，步骤S2中：
所述上升段拟合函数为α＝DRe2+ERe+F，Re∈[Re1、Re2)；
所述下降段拟合函数为α＝ARe2+BRe+C，Re∈[Re2、Re3]；
式中α为流量系数，Re为雷诺数，A、B、C、D、E和F均为常数，Re1为上升段的雷诺数起点，Re2为上升段和下降段交点处的雷诺数，Re3为下降段的雷诺数终点。
3.根据权利要求2所述的一种核反应堆模拟燃料组件流量测量方法，其特征在于，步骤S3中：
文丘里管的理论关系函数为A1*α-B1*Re＝0；其中
Figure FDA0002642230600000011
Figure FDA0002642230600000021
上式中，ρ为入口管取压点处流体密度，根据流体介质温度和压力确定；ΔP为文丘里管差压；d为文丘里管喉部直径；ε为膨胀系数，根据直径比、压力比和等熵指数确定；υ为入口管流体运动粘度，根据流体介质温度和压力确定；D为文丘里管入口管直径。
4.根据权利要求3所述的一种核反应堆模拟燃料组件流量测量方法，其特征在于，步骤S3中：
下降段模型为A2Re2+B2Re+C2＝0；其中A2＝A1*A，B2＝A1*B-B1，C2＝A1*C；
上升段模型为A3Re2+B3Re+C3＝0；其中A3＝A1*D，B3＝A1*E-B1，C3＝A1*F。
5.根据权利要求4所述的一种核反应堆模拟燃料组件流量测量方法，其特征在于，步骤S4包括以下子步骤：
下降段输出雷诺数Re′包括两个解：Re1′和Re2′；其中：
Figure FDA0002642230600000022
当Re1′和Re2′中任意一个位于[Re2、Re3]内时，将位于[Re2、Re3]内的下降段输出雷诺数Re′作为实际雷诺数输出。
6.根据权利要求4所述的一种核反应堆模拟燃料组件流量测量方法，其特征在于，步骤S4包括以下子步骤：
上升段输出雷诺数Re″包括两个解：Re3′和Re4′；其中：
Figure FDA0002642230600000023
Figure FDA0002642230600000024
当Re3′和Re4′中任意一个位于[Re1、Re2)内时，将位于[Re1、Re2)内的上升段输出雷诺数Re″作为实际雷诺数输出。
7.根据权利要求4所述的一种核反应堆模拟燃料组件流量测量方法，其特征在于，步骤S5包括以下子步骤：
根据下式获取模拟燃料组件中流体的流量：
Q＝B1*Re”'；
式中Re”'为实际雷诺数，Q为模拟燃料组件中流体的流量。
8.采用权利要求1～7中任意一项所述方法的一种核反应堆模拟燃料组件流量测量系统，其特征在于，包括：
标定单元：用于使用文丘里管对模拟燃料组件中流体的流量系数α和雷诺数Re的关系进行标定生成标定关系函数；所述标定关系函数的自变量为雷诺数Re，应变量为流量系数α；
拟合单元：用于对所述标定关系函数中上升段的进行拟合生成上升段拟合函数；对所述标定关系函数中下降段的进行拟合生成下降段拟合函数；
模型单元：用于将所述上升段拟合函数和文丘里管的理论关系函数联立生成上升段模型；将所述下降段拟合函数和文丘里管的理论关系函数联立生成下降段模型；所述文丘里管的理论关系函数根据文丘里管中雷诺数Re和流量系数α的理论关系获取；
测量单元：用于对模拟燃料组件中流体的流量进行测量时，将从文丘里管和测温仪表获取的数据输入所述下降段模型获取下降段输出雷诺数Re′；
当下降段输出雷诺数Re′位于所述标定关系函数的下降段内时，将所述下降段输出雷诺数Re′作为实际雷诺数输出；
当下降段输出雷诺数Re′未位于所述标定关系函数的下降段内时，将从文丘里管和测温仪表获取的数据输入所述上升段模型获取上升段输出雷诺数Re″；
如果所述上升段输出雷诺数Re″位于所述标定关系函数的上升段内时，将所述上升段输出雷诺数Re″作为实际雷诺数输出；
如果所述上升段输出雷诺数Re″未位于所述标定关系函数的上升段内时，扩大流量系数α和雷诺数Re的关系的标定范围，且标定单元、拟合单元、模型单元和测量单元循环工作；
所述测量单元还用于根据所述实际雷诺数获取模拟燃料组件中流体的流量。
9.根据权利要求8所述的一种核反应堆模拟燃料组件流量测量系统，其特征在于，所述上升段拟合函数为α＝DRe2+ERe+F，Re∈[Re1、Re2)；
所述下降段拟合函数为α＝ARe2+BRe+C，Re∈[Re2、Re3]；
式中α为流量系数，Re为雷诺数，A、B、C、D、E和F均为常数，Re1为上升段的雷诺数起点，Re2为上升段和下降段交点处的雷诺数，Re3为下降段的雷诺数终点。
10.根据权利要求8所述的一种核反应堆模拟燃料组件流量测量系统，其特征在于，文丘里管的理论关系函数为A1*α-B1*Re＝0；其中
Figure FDA0002642230600000041
上式中，ρ为入口管取压点处流体密度，根据流体介质温度和压力确定；ΔP为文丘里管差压；d为文丘里管喉部直径；ε为膨胀系数，根据直径比、压力比和等熵指数确定；υ为入口管流体运动粘度，根据流体介质温度和压力确定；D为文丘里管入口管直径。
CN202010843390.4A 2020-08-20 2020-08-20 一种核反应堆模拟燃料组件流量测量方法及系统 Active CN111982213B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010843390.4A CN111982213B (zh) 2020-08-20 2020-08-20 一种核反应堆模拟燃料组件流量测量方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010843390.4A CN111982213B (zh) 2020-08-20 2020-08-20 一种核反应堆模拟燃料组件流量测量方法及系统

Publications (2)

Publication Number Publication Date
CN111982213A true CN111982213A (zh) 2020-11-24
CN111982213B CN111982213B (zh) 2022-02-11

Family

ID=73444179

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010843390.4A Active CN111982213B (zh) 2020-08-20 2020-08-20 一种核反应堆模拟燃料组件流量测量方法及系统

Country Status (1)

Country Link
CN (1) CN111982213B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113295395A (zh) * 2021-05-24 2021-08-24 无锡职业技术学院 一种电磁比例阀流量系数测试系统及方法

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5669515A (en) * 1979-11-10 1981-06-10 Toshiba Corp Measuring device of recycle flow rate in furnace core
US4425807A (en) * 1982-02-11 1984-01-17 Michael Victor Flow measuring device with constant flow coefficient
US4634559A (en) * 1984-02-29 1987-01-06 Aluminum Company Of America Fluid flow control process
CN1563940A (zh) * 2004-04-02 2005-01-12 浙江大学 基于希尔伯特－黄变换的水平管气液两相流流型识别方法
CN101149991A (zh) * 2007-10-19 2008-03-26 中国核动力研究设计院 用于燃料组件流量测量的模拟组件
EP2192391A1 (en) * 2008-12-01 2010-06-02 Services Pétroliers Schlumberger Apparatus and a method of measuring the flow of a fluid
RU2008150138A (ru) * 2008-12-19 2010-06-27 Иван Федорович Пивин (RU) Датчик
RU2010152725A (ru) * 2010-12-23 2012-06-27 Государственное образовательное учреждение высшего профессионального образования "Московский энергетический институт (технический у Способ работы тепловыделяющей сборки на входном участке и устройство для его осуществления
CN103353908A (zh) * 2013-06-20 2013-10-16 江苏大学 一种基于数值计算的管路阻力系数精确计算方法
US20140137663A1 (en) * 2012-11-20 2014-05-22 Fluid Handling Llc Valve ball for direct flow measurement
CN104897222A (zh) * 2014-03-07 2015-09-09 华东理工大学 密相气力输送过程的固相流量测量系统与测量方法
CN105973319A (zh) * 2016-07-15 2016-09-28 中国核动力研究设计院 一种控制棒驱动机构排污系统水力特性计算方法
EP3128212A1 (en) * 2015-08-04 2017-02-08 WATTS INDUSTRIES ITALIA S.r.l. Instrument for measuring the flowrate of a fluid
EP3322962A1 (en) * 2015-07-14 2018-05-23 Rabotti S.r.l. Method for the measurement of an unsteady liquid flow rate, in particular of a high pressure liquid flow
CN108983634A (zh) * 2017-06-02 2018-12-11 数学工程公司 用于重新缩放可执行仿真模型的系统和方法
CN110487338A (zh) * 2019-08-29 2019-11-22 东南大学 一种多孔平衡孔板设计方法及其评价方法

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5669515A (en) * 1979-11-10 1981-06-10 Toshiba Corp Measuring device of recycle flow rate in furnace core
US4425807A (en) * 1982-02-11 1984-01-17 Michael Victor Flow measuring device with constant flow coefficient
US4634559A (en) * 1984-02-29 1987-01-06 Aluminum Company Of America Fluid flow control process
CN1563940A (zh) * 2004-04-02 2005-01-12 浙江大学 基于希尔伯特－黄变换的水平管气液两相流流型识别方法
CN101149991A (zh) * 2007-10-19 2008-03-26 中国核动力研究设计院 用于燃料组件流量测量的模拟组件
EP2192391A1 (en) * 2008-12-01 2010-06-02 Services Pétroliers Schlumberger Apparatus and a method of measuring the flow of a fluid
RU2008150138A (ru) * 2008-12-19 2010-06-27 Иван Федорович Пивин (RU) Датчик
RU2010152725A (ru) * 2010-12-23 2012-06-27 Государственное образовательное учреждение высшего профессионального образования "Московский энергетический институт (технический у Способ работы тепловыделяющей сборки на входном участке и устройство для его осуществления
US20140137663A1 (en) * 2012-11-20 2014-05-22 Fluid Handling Llc Valve ball for direct flow measurement
CN103353908A (zh) * 2013-06-20 2013-10-16 江苏大学 一种基于数值计算的管路阻力系数精确计算方法
CN104897222A (zh) * 2014-03-07 2015-09-09 华东理工大学 密相气力输送过程的固相流量测量系统与测量方法
EP3322962A1 (en) * 2015-07-14 2018-05-23 Rabotti S.r.l. Method for the measurement of an unsteady liquid flow rate, in particular of a high pressure liquid flow
EP3128212A1 (en) * 2015-08-04 2017-02-08 WATTS INDUSTRIES ITALIA S.r.l. Instrument for measuring the flowrate of a fluid
CN105973319A (zh) * 2016-07-15 2016-09-28 中国核动力研究设计院 一种控制棒驱动机构排污系统水力特性计算方法
CN108983634A (zh) * 2017-06-02 2018-12-11 数学工程公司 用于重新缩放可执行仿真模型的系统和方法
CN110487338A (zh) * 2019-08-29 2019-11-22 东南大学 一种多孔平衡孔板设计方法及其评价方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
MIKUZ, BLAZ: "Numerical simulations of a turbulent flow in a fuel assembly", 《23RD INTERNATIONAL CONFERENCE NUCLEAR ENERGY FOR NEW EUROPE, (NENE 2014)》 *
杨红义: "堆芯燃料组件流阻特性模块式与整体式三维数值模拟方案比较", 《原子能科学技术》 *
王盛: "CNP1000反应堆旁漏流、流量分配及交混特性研究", 《中国核动力研究设计院科学技术年报（2009年）》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113295395A (zh) * 2021-05-24 2021-08-24 无锡职业技术学院 一种电磁比例阀流量系数测试系统及方法

Also Published As

Publication number Publication date
CN111982213B (zh) 2022-02-11

Similar Documents

Publication Publication Date Title
Stadtmüller et al. A test case for the numerical investigation of wake passing effects on a highly loaded LP turbine cascade blade
Aly et al. Experimental study of the pressure drop after fractal-shaped orifices in turbulent pipe flows
CN110411521B (zh) 一种基于双喷嘴的油井多相流分相含率在线计量方法
Grimshaw et al. Fast settling millimetre-scale five-hole probes
CN111982213B (zh) 一种核反应堆模拟燃料组件流量测量方法及系统
CN106096223A (zh) 一种五孔探针数据处理方法
CN111503025A (zh) 一种低压比轴流压缩机模型级级性能计算方法
Ismail et al. Experimental evaluation for the feasibility of test chamber in the open-loop wind tunnel
CN115824324B (zh) 一种基于旋进旋涡mems检测技术的湿气流量测量方法
CN115358026A (zh) 一种基于多元线性回归与曲面拟合的五孔探针数据处理方法
Crawford Design and Calibration of Seven Hole Probes for Flow Measurement
CN108593964B (zh) 一种七孔气流探针的标定方法
CN114993392B (zh) 基于加速度幅值建模的涡街湿气流量测量方法
Boufidi et al. A probabilistic uncertainty estimation method for turbulence parameters measured by hot-wire anemometry in short-duration wind tunnels
CN108280300B (zh) 基于计算流体力学的大数据流量计开发方法
CN106017582A (zh) 一种基于管径二分原理的皮托管流量测量方法
Toni et al. Unsteady flow field measurements in an industrial centrifugal compressor
Ku¨ nzelmann et al. Steady and unsteady flow field in a multistage low-speed axial compressor: a test case
Singh et al. Investigation the Effect on Coefficient of filled and unfilled port to port dimension with plaster of paris of 2-Holes offset probe
Chen et al. Performance studies on cavitation-resistance turbine flow sensor based on experiment and CFD simulation
Ajit Experimental work on the low Reynolds number behaviour of 2-hole offset probes
Burdett et al. An impulse response technique to improve the effective frequency response of pressure probes
Piradeepan An experimental and numerical investigation of a turbulent airfoil wake in a 90° curved duct
CN106909707B (zh) 利用分形维数测量风场地表粗糙度的方法
ÖZAHİ et al. DEVISED APPLICATION OF LABVIEW FOR AN AUTOMATIC TEST SYSTEM BASED ON GENERATION, CONTROL AND PROCESSING OF PULSATILE PIPE FLOWS.

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