CN109977598B - 针对阀下游排放管的载荷分析方法 - Google Patents
针对阀下游排放管的载荷分析方法 Download PDFInfo
- Publication number
- CN109977598B CN109977598B CN201910280875.4A CN201910280875A CN109977598B CN 109977598 B CN109977598 B CN 109977598B CN 201910280875 A CN201910280875 A CN 201910280875A CN 109977598 B CN109977598 B CN 109977598B
- Authority
- CN
- China
- Prior art keywords
- valve
- discharge pipe
- model
- downstream
- thermal hydraulic
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- 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/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了针对阀下游排放管的载荷分析方法,为了计算获得排放载荷,需要首先进行初始和边界条件设定、然后进行热工水力分析、最后将热工水力分析结果通过高效便捷的方法转换为载荷。其中,比较关键的技术点在于:热工水力分析过程中关键模型的选取以及怎样将热工水力分析结果通过一定的手段转换为载荷结果。
Description
技术领域
本发明涉及核电站排放管线载荷分析领域,具体涉及针对阀下游排放管的载荷分析方法。
背景技术
如图1所示,水封体排放管线是由依次连接的稳压器2、U型管状的水封体管5、安全阀3、阀下游排放管1、卸压箱4组成;
在安全阀(或卸压阀)的上游,由于考虑到稳压器安全阀的不凝性气体泄漏等因素,安全阀(或卸压阀)与稳压器汽空间连接管线上有一段底部充水的U型管状的水封体管5,起到水封的作用以密封不凝性气体、保持稳压器安全阀体的温度等。
当任一安全阀(或卸压阀)开启时,水封体管5内的水封体会在压力差的作用下、在很短的时间内(不到1秒)通过阀下游排放管,从而对阀下游排放管造成一定的冲击力。为了确保该瞬态对阀下游排放管造成的冲击力在阀下游排放管的承受范围之内,需要通过热工水力计算及力学处理过程来研究瞬态过程中产生的排放载荷,从而为阀下游排放管设计(包括材料、布置、尺寸等)、支撑设计提供技术支持。
目前,国内没有能够计算水封体排放管线中阀下游排放管载荷的分析方法。
发明内容
本发明目的在于提供针对阀下游排放管的载荷分析方法。
本发明通过下述技术方案实现:
针对水封体排放管线的载荷分析模型构建方法,包括以下步骤:
在能够模拟两相传热传质过程、临界流、汽液界面追踪的平台中,构建阀下游排放管模型,同时构建稳压器模型、水封体管模型、安全阀模型、卸压箱模型,将稳压器模型、水封体管模型、安全阀模型、阀下游排放管模型、卸压箱模型通过边界连接和数据传递连接的方式连通构成水封体排放管线模型;其中,阀下游排放管模型是由n个节块模型通过边界连接和数据传递连接的方式连通构成,n个节块模型一一对应于沿阀下游排放管轴向、对阀下游排放管分段划分为形成的n个控制体节块,每段控制体节块的长度控制在10cm-20cm的范围内;
所述稳压器模型包括与稳压器对应的热工水力分析模型B,安全阀模型包括与安全阀对应的热工水力分析模型,卸压箱模型包括与卸压箱对应的热工水力分析模型D,水封体管模型包括与水封体管对应的热工水力分析模型E,节块模型包括与控制体节块对应的热工水力分析模型a,阀下游排放管模型包括与阀下游排放管对应的热工水力分析模型A,热工水力分析模型A由n个热工水力分析模型a通过边界连接和数据传递连接的方式连通构成,水封体排放管线模型包括由热工水力分析模型A、热工水力分析模型B、热工水力分析模型C、热工水力分析模型D、热工水力分析模型E通过边界连接和数据传递连接的方式连通构成与水封体排放管线对应的热工水力分析模型Z;
构建阀下游排放管载荷分析模型:以热工水力分析模型a输出的热工水力分析结果信息集a作为输入变量和热工水力分析模型A输出的热工水力分析结果信息集A作为输入变量,以阀下游排放管载荷为解释变量,从而构建起阀下游排放管载荷分析模型。
本发明的设计原理为:
在现有技术中,对于阀下游排放管载荷的分析一般采用的方法是:直接采用软件建立力学模型,其做法是将稳压器排放***分为稳压器和管道***两个部分,分别采用不同的简化方法进行模拟,最后将2个有限元模型耦合在一起,从而对整个稳压器排放***进行力学分析。
但是,本发明的针对对象是具有水封体的管道***,当安全阀(或卸压阀)开启时,由于阀门(阀门上游高达17.0MPa左右,阀门下游0.1MPa左右)两侧压力差巨大,因此水封会以临界流速通过阀门,并以高达几十甚至几百米每秒的流速通过下游各排放管道。因此水封体在运动的过程中会汽化,因此阀下游排放管内部是复杂的两相流,包括:泡状流、弹状流、环状流、弥散流,并且由于压力传递、汽液相间的阻力和换热,管道内部的压力、汽液相密度、汽液相温度、汽液相流速、空泡份额等都随时在变化,每个管道所受的动态冲击载荷也是实时变化的。由于管道内部流体的压力、流速、加速度、空泡份额、热物性等随时间一直在变化,因此管道会受到随时间变化的冲击载荷,冲击载荷的大小与内部流体的运动状态有关。因此,传统直接采用力学法分析的方法将会导致分析结果与实际情况差别较大,且由于两相态的阀下游排放管所受动态冲击载荷,传统整体管道一体化的力学分析,并未考虑这种影响,也是影响分析结果与实际结果差异较大的重要因素。本发明采用将复杂两相流动热工水力参量和管道所受载荷联系起来的方法,同时在对阀下游排放管进行热工水力分析时采用多控制体节块的方法进行分析,这样使得最终的热工水力参量、以及每个控制体节块热工水力参量都能与管道所受载荷关联,这样就使得整个分析结果更加考虑到实际的内部载荷情况,也使得分析结果更为精准。
在本发明中,由于考虑管道不同位置的情况处于两相变化态的,为了分析更加可靠,因此,本发明建立分析模型时,与常规的以大部件为分析节点的方式不同,且并不是直接利用力学角度出发进行的载荷分析,对于上述内容,本发明的主要技术贡献为:1、提出了先建立两相下热工水力分析模型,然后对热工水力分析模型输出的参量用于阀下游排放管载荷分析的模型建立,即本发明的阀下游排放管载荷分析模型时在热工水力分析模型输出的参量的基础上建立的,这样在整个模型的过程中,并不是一个静态数据的简单分析,而是在一个动态的热工水力过程中,不断获得各种动态数据,然后依赖这些动态数据输出构建具有可动态分析的载荷分析模型。2、由于本发明是针对的是具有水封体的管线,因此阀下游排放管是两相环境,且这种两相状态是变化的,因此为了更加准确的分析结果,提出了针对阀下游排放管多分控制体节块的方法进行动态热工水力分析,即对阀下游排放管分为多个段,每段形成1个控制体节块,每段都独立、但又互相关联影响,因此,每个控制体节块都能输出对应的热工水力分析参数,那么在此基础上,本发明在建立阀下游排放管载荷分析模型时,就可以依靠每段控制体节块的热工水力分析参数和整体阀下游排放管的总体参数进行对应的载荷力学分析。
优选的,上述阀下游排放管载荷分析模型优选通过对管道以及流体进行动量守恒分析以及受力分析,从而获得管道所受载荷的计算。因此所述阀下游排放管载荷分析模型是由加速项模型、动能变化模型、压降模型、流体体积力模型求和处理组成;其中加速项模型是指对各个控制体节块内的加速项进行积分处理;动能变化模型是指对整个阀下游排放管的入口和出口的动能变化进行处理;压降模型是指对整个阀下游排放管的入口和出口的压降变化进行处理;流体体积力模型是指对各个控制体节块内的体积力进行积分处理。
优选的,所述热工水力分析结果信息集a包括:控制体节块内的液相平均流速、控制体节块内的汽相平均流速、控制体节块内的液相密度、控制体节块内汽相密度、制体节块内的空泡份额;所述热工水力分析结果信息集A包括:阀下游排放管的入口截面处压力、阀下游排放管的出口截面处压力;阀下游排放管的入口液相流速、阀下游排放管的出口液相流速、阀下游排放管的入口汽相流速、阀下游排放管的出口汽相流速、阀下游排放管的入口空泡份额、阀下游排放管的出口空泡份额;阀下游排放管的入口液相密度、阀下游排放管的出口液相密度;阀下游排放管的入口汽相密度、阀下游排放管的入口汽相密度。
优选的,所述阀下游排放管载荷分析模型是由加速项模型、动能变化模型、压降模型、流体体积力模型求和处理组成;其中加速项模型是指对各个控制体节块内的加速项进行积分处理;动能变化模型是指对整个阀下游排放管的入口和出口的动能变化进行处理;压降模型是指对整个阀下游排放管的入口和出口的压降变化进行处理;流体体积力模型是指对各个控制体节块内的体积力进行积分处理;阀下游排放管载荷分析模型具体为:
其中,
表示阀下游排放管载荷;t为时间;S入表示阀下游排放管的入口截面积、S出表示阀下游排放管的出口截面积,S入和S出为基础信息、由配置获得;P入表示阀下游排放管的入口截面处压力、P出表示阀下游排放管的出口截面处压力;ul入表示阀下游排放管的入口液相流速、ul出表示阀下游排放管的出口液相流速、ug入表示阀下游排放管的入口汽相流速、ug出表示阀下游排放管的出口汽相流速、a入表示阀下游排放管的入口空泡份额、a出表示阀下游排放管的出口空泡份额;ρl入表示阀下游排放管的入口液相密度、ρl出表示阀下游排放管的出口液相密度;ρg入表示阀下游排放管的入口汽相密度、ρg出表示阀下游排放管的出口汽相密度;
g为重力加速度,dVi表示第i个控制体节块的体积;
优选的,所述平台为RELAP5。
优选的,所述热工水力分析模型包括临界流模型,所述临界流模型包括Henry-Fauske模型、非均匀流模型、突扩突缩模型。
针对阀下游排放管的载荷分析方法:包括以下步骤:
采集目标水封体排放管线的基础信息,将基础信息分配给热工水力分析模型Z,基础信息包括水封体的信息、安全阀的信息、阀下游排放管的几何信息、卸压箱的几何信息、水封体排放管线的内部流体的初始条件和边界条件;
设置热工水力分析模型Z的分析时间步长;
通过n个热工水力分析模型a获得n个随时间实时变化的热工水力分析结果信息集a,通过热工水力分析模型A获得随时间实时变化的热工水力分析结果信息集A;热工水力分析结果信息集a包括:控制体节块内的液相平均流速、控制体节块内的汽相平均流速、控制体节块内的液相密度、控制体节块内汽相密度、制体节块内的空泡份额;所述热工水力分析结果信息集A包括:阀下游排放管的入口截面处压力、阀下游排放管的出口截面处压力;阀下游排放管的入口液相流速、阀下游排放管的出口液相流速、阀下游排放管的入口汽相流速、阀下游排放管的出口汽相流速、阀下游排放管的入口空泡份额、阀下游排放管的出口空泡份额;阀下游排放管的入口液相密度、阀下游排放管的出口液相密度;阀下游排放管的入口汽相密度、阀下游排放管的入口汽相密度;
通过阀下游排放管载荷分析模型对n个随时间实时变化的热工水力分析结果信息集a和随时间实时变化的热工水力分析结果信息集A进行分析,得到随时间实时变化的阀下游排放管载荷值。
优选的,得到随时间实时变化的阀下游排放管载荷值后,将离散的阀下游排放管载荷值拟合成随时间实时变化的载荷波。
优选的,所述分析时间步长控制在10-4秒的量级。
优选的,水封体的信息包括水封体体积、水封温度分布,安全阀的信息包括阀门开启压力、阀门排量、阀门完全开启时间,阀下游排放管的几何信息包括管道尺寸、管道布置、弯头角度及弯曲半径,卸压箱的几何信息包括卸压箱尺寸,水封体排放管线的内部流体的初始条件和边界条件包括气体种类、初始压力、初始温度。
本发明与现有技术相比,具有如下的优点和有益效果:
本发明的总体设想是:
为了计算获得排放载荷,需要首先进行初始和边界条件设定、然后进行热工水力分析、最后将热工水力分析结果通过高效便捷的方法转换为载荷。其中,比较关键的技术点在于:热工水力分析过程中关键模型的选取以及怎样将热工水力分析结果通过一定的手段转换为载荷结果。
上述方法准备用于中国三代压水堆核电厂(站)的设计。该研究方法的分析成果将直接应用于管道应力分析,从而为管道材料设计、管道尺寸及布置、管道支撑设计提供技术支持。此外,该分析方法还可以用于论证稳压器安全阀(或卸压阀)定值设计的合理性问题。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1为水封体排放管线示意图;
图2为本发明流程示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
针对水封体排放管线的载荷分析模型构建方法,包括以下步骤:
在能够模拟两相传热传质过程、临界流、汽液界面追踪的平台中,构建阀下游排放管模型,同时构建稳压器模型、水封体管模型、安全阀模型、卸压箱模型,将稳压器模型、水封体管模型、安全阀模型、阀下游排放管模型、卸压箱模型通过边界连接和数据传递连接的方式连通构成水封体排放管线模型;其中,阀下游排放管模型是由n个节块模型通过边界连接和数据传递连接的方式连通构成,n个节块模型一一对应于沿阀下游排放管轴向、对阀下游排放管分段划分为形成的n个控制体节块,每段控制体节块的长度控制在10cm-20cm的范围内;
所述稳压器模型包括与稳压器2对应的热工水力分析模型B,安全阀模型包括与安全阀3对应的热工水力分析模型C,卸压箱模型包括与卸压箱4对应的热工水力分析模型D,水封体管模型5包括与水封体管对应的热工水力分析模型E,节块模型包括与控制体节块对应的热工水力分析模型a,阀下游排放管模型包括与阀下游排放管1对应的热工水力分析模型A,热工水力分析模型A由n个热工水力分析模型a通过边界连接和数据传递连接的方式连通构成,水封体排放管线模型包括由热工水力分析模型A、热工水力分析模型B、热工水力分析模型C、热工水力分析模型D、热工水力分析模型E通过边界连接和数据传递连接的方式连通构成与水封体排放管线对应的热工水力分析模型Z;
构建阀下游排放管载荷分析模型:以热工水力分析模型a输出的热工水力分析结果信息集a作为输入变量和热工水力分析模型A输出的热工水力分析结果信息集A作为输入变量,以阀下游排放管载荷为解释变量,从而构建起阀下游排放管载荷分析模型。
在本发明中,例如如图1中,对于阀下游排放管,我们对其在平台中构建起对于的模型,对于一段,阀下游排放管假设其为1m长,对于在建立模型时,则将其划分为10cm一段,则有10个控制体节块,因此***输出10组热工水力分析结果信息集a,热工水力分析结果信息集A则是以整个阀下游排放管的入口和出口作为考虑对象得出的。
本发明的设计原理为:
在现有技术中,对于阀下游排放管载荷的分析一般采用的方法是:直接采用软件建立力学模型,其做法是将稳压器排放***分为稳压器和管道***两个部分,分别采用不同的简化方法进行模拟,最后将2个有限元模型耦合在一起,从而对整个稳压器排放***进行力学分析。
但是,本发明的针对对象是具有水封体的管道***,当安全阀(或卸压阀)开启时,由于阀门(阀门上游高达17.0MPa左右,阀门下游0.1MPa左右)两侧压力差巨大,因此水封会以临界流速通过阀门,并以高达几十甚至几百米每秒的流速通过下游各排放管道。因此水封体在运动的过程中会汽化,因此阀下游排放管内部是复杂的两相流,包括:泡状流、弹状流、环状流、弥散流,并且由于压力传递、汽液相间的阻力和换热,管道内部的压力、汽液相密度、汽液相温度、汽液相流速、空泡份额等都随时在变化,每个管道所受的动态冲击载荷也是实时变化的。由于管道内部流体的压力、流速、加速度、空泡份额、热物性等随时间一直在变化,因此管道会受到随时间变化的冲击载荷,冲击载荷的大小与内部流体的运动状态有关。因此,传统直接采用力学法分析的方法将会导致分析结果与实际情况差别较大,且由于两相态的阀下游排放管所受动态冲击载荷,传统整体管道一体化的力学分析,并未考虑这种影响,也是影响分析结果与实际结果差异较大的重要因素。本发明采用将复杂两相流动热工水力参量和管道所受载荷联系起来的方法,同时在对阀下游排放管进行热工水力分析时采用多控制体节块的方法进行分析,这样使得最终的热工水力参量、以及每个控制体节块热工水力参量都能与管道所受载荷关联,这样就使得整个分析结果更加考虑到实际的内部载荷情况,也使得分析结果更为精准。
在本发明中,由于考虑管道不同位置的情况处于两相变化态的,为了分析更加可靠,因此,本发明建立分析模型时,与常规的以大部件为分析节点的方式不同,且并不是直接利用力学角度出发进行的载荷分析,对于上述内容,本发明的主要技术贡献为:1、提出了先建立两相下热工水力分析模型,然后对热工水力分析模型输出的参量用于阀下游排放管载荷分析的模型建立,即本发明的阀下游排放管载荷分析模型时在热工水力分析模型输出的参量的基础上建立的,这样在整个模型的过程中,并不是一个静态数据的简单分析,而是在一个动态的热工水力过程中,不断获得各种动态数据,然后依赖这些动态数据输出构建具有可动态分析的载荷分析模型。2、由于本发明是针对的是具有水封体的管线,因此阀下游排放管是两相环境,且这种两相状态是变化的,因此为了更加准确的分析结果,提出了针对阀下游排放管多分控制体节块的方法进行动态热工水力分析,即对阀下游排放管分为多个段,每段形成1个控制体节块,每段都独立、但又互相关联影响,因此,每个控制体节块都能输出对应的热工水力分析参数,那么在此基础上,本发明在建立阀下游排放管载荷分析模型时,就可以依靠每段控制体节块的热工水力分析参数和整体阀下游排放管的总体参数进行对应的载荷力学分析。
优选的,上述阀下游排放管载荷分析模型优选通过对管道以及流体进行动量守恒分析以及受力分析,从而获得管道所受载荷的计算。因此所述阀下游排放管载荷分析模型是由加速项模型、动能变化模型、压降模型、流体体积力模型求和处理组成;其中加速项模型是指对各个控制体节块内的加速项进行积分处理;动能变化模型是指对整个阀下游排放管的入口和出口的动能变化进行处理;压降模型是指对整个阀下游排放管的入口和出口的压降变化进行处理;流体体积力模型是指对各个控制体节块内的体积力进行积分处理。
优选的,所述热工水力分析结果信息集a包括:控制体节块内的液相平均流速、控制体节块内的汽相平均流速、控制体节块内的液相密度、控制体节块内汽相密度、制体节块内的空泡份额;所述热工水力分析结果信息集A包括:阀下游排放管的入口截面处压力、阀下游排放管的出口截面处压力;阀下游排放管的入口液相流速、阀下游排放管的出口液相流速、阀下游排放管的入口汽相流速、阀下游排放管的出口汽相流速、阀下游排放管的入口空泡份额、阀下游排放管的出口空泡份额;阀下游排放管的入口液相密度、阀下游排放管的出口液相密度;阀下游排放管的入口汽相密度、阀下游排放管的入口汽相密度。
优选的,所述阀下游排放管载荷分析模型是由加速项模型、动能变化模型、压降模型、流体体积力模型求和处理组成;其中加速项模型是指对各个控制体节块内的加速项进行积分处理;动能变化模型是指对整个阀下游排放管的入口和出口的动能变化进行处理;压降模型是指对整个阀下游排放管的入口和出口的压降变化进行处理;流体体积力模型是指对各个控制体节块内的体积力进行积分处理;阀下游排放管载荷分析模型具体为:
其中,
表示阀下游排放管载荷;t为时间;S入表示阀下游排放管的入口截面积、S出表示阀下游排放管的出口截面积,S入和S出为基础信息、由配置获得;P入表示阀下游排放管的入口截面处压力、P出表示阀下游排放管的出口截面处压力;ul入表示阀下游排放管的入口液相流速、ul出表示阀下游排放管的出口液相流速、ug入表示阀下游排放管的入口汽相流速、ug出表示阀下游排放管的出口汽相流速、a入表示阀下游排放管的入口空泡份额、a出表示阀下游排放管的出口空泡份额;ρl入表示阀下游排放管的入口液相密度、ρl出表示阀下游排放管的出口液相密度;ρg入表示阀下游排放管的入口汽相密度、ρg出表示阀下游排放管的出口汽相密度;
g为重力加速度,dVi表示第i个控制体节块的体积。
优选的,所述平台为RELAP5。
优选的,所述热工水力分析模型包括临界流模型,所述临界流模型包括Henry-Fauske模型、非均匀流模型、突扩突缩模型。
实施例2
如图2所示
针对阀下游排放管的载荷分析方法:包括以下步骤:
对具体的目标管线机械能分析,确定其需要分析的工况,从而采集目标水封体排放管线的基础信息,将基础信息分配给热工水力分析模型Z(瞬态分析),基础信息包括水封体的信息、安全阀的信息、阀下游排放管的几何信息、卸压箱的几何信息、水封体排放管线的内部流体的初始条件和边界条件;
设置热工水力分析模型Z的分析时间步长;
通过n个热工水力分析模型a获得n个随时间实时变化的热工水力分析结果信息集a,通过热工水力分析模型A获得随时间实时变化的热工水力分析结果信息集A;热工水力分析结果信息集a包括:控制体节块内的液相平均流速、控制体节块内的汽相平均流速、控制体节块内的液相密度、控制体节块内汽相密度、制体节块内的空泡份额;所述热工水力分析结果信息集A包括:阀下游排放管的入口截面处压力、阀下游排放管的出口截面处压力;阀下游排放管的入口液相流速、阀下游排放管的出口液相流速、阀下游排放管的入口汽相流速、阀下游排放管的出口汽相流速、阀下游排放管的入口空泡份额、阀下游排放管的出口空泡份额;阀下游排放管的入口液相密度、阀下游排放管的出口液相密度;阀下游排放管的入口汽相密度、阀下游排放管的入口汽相密度;
通过阀下游排放管载荷分析模型对n个随时间实时变化的热工水力分析结果信息集a和随时间实时变化的热工水力分析结果信息集A进行分析,得到随时间实时变化的阀下游排放管载荷值。阀下游排放管载荷分析模型采用实施例1方案中的公式。
优选的,得到随时间实时变化的阀下游排放管载荷值后,将离散的阀下游排放管载荷值拟合成随时间实时变化的载荷波。
优选的,所述分析时间步长控制在10-4秒的量级。对于时间步长的设置,由于本发明是一个瞬态过程,若时间步长设置的过长(例如10-3)会影响计算精度,经过我们测试,发现控制在10-4秒的量级,对于整个分析能起到非常好的精度控制。对于10-4秒的量级,则会形成1万个样点,而若步长设置的过短(例如10-5)则会导致失真,且对于计算效率也是非常差的。
优选的,水封体的信息包括水封体体积、水封温度分布,安全阀的信息包括阀门开启压力、阀门排量、阀门完全开启时间,阀下游排放管的几何信息包括管道尺寸、管道布置、弯头角度及弯曲半径,卸压箱的几何信息包括卸压箱尺寸,水封体排放管线的内部流体的初始条件和边界条件包括气体种类、初始压力、初始温度。
其中若有三组安全阀,且阀门定值不同,可以这样选定分析工况:
①若通过超压分析得出不同定值的阀门开启时间间隔较大(大于0.5s),则分析最大定值安全阀单独开启的工况;
②若通过超压分析得出不同定值的阀门开启时间间隔较小(小于0.5s),则分析不同定值安全阀相继开启的工况。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.针对阀下游排放管的载荷分析方法,其特征在于,包括以下步骤:
采集目标水封体排放管线的基础信息,将基础信息分配给热工水力分析模型Z,基础信息包括水封体的信息、安全阀的信息、阀下游排放管的几何信息、卸压箱的几何信息、水封体排放管线的内部流体的初始条件和边界条件;
设置热工水力分析模型Z的分析时间步长;
通过n个热工水力分析模型a获得n个随时间实时变化的热工水力分析结果信息集a,通过热工水力分析模型A获得随时间实时变化的热工水力分析结果信息集A;热工水力分析结果信息集a包括:控制体节块内的液相平均流速、控制体节块内的汽相平均流速、控制体节块内的液相密度、控制体节块内汽相密度、制体节块内的空泡份额;所述热工水力分析结果信息集A包括:阀下游排放管的入口截面处压力、阀下游排放管的出口截面处压力;阀下游排放管的入口液相流速、阀下游排放管的出口液相流速、阀下游排放管的入口汽相流速、阀下游排放管的出口汽相流速、阀下游排放管的入口空泡份额、阀下游排放管的出口空泡份额;阀下游排放管的入口液相密度、阀下游排放管的出口液相密度;阀下游排放管的入口汽相密度、阀下游排放管的出口汽相密度;
通过阀下游排放管载荷分析模型对n个随时间实时变化的热工水力分析结果信息集a和随时间实时变化的热工水力分析结果信息集A进行分析,得到随时间实时变化的阀下游排放管载荷值;
其中,
在能够模拟两相传热传质过程、临界流、汽液界面追踪的平台中,构建阀下游排放管模型,同时构建稳压器模型、水封体管模型、安全阀模型、卸压箱模型,将稳压器模型、水封体管模型、安全阀模型、阀下游排放管模型、卸压箱模型通过边界连接和数据传递连接的方式连通构成水封体排放管线模型;其中,阀下游排放管模型是由n个节块模型通过边界连接和数据传递连接的方式连通构成,n个节块模型一一对应于沿阀下游排放管轴向、对阀下游排放管分段划分为形成的n个控制体节块,每段控制体节块的长度控制在10cm-20cm的范围内;
所述稳压器模型包括与稳压器(2)对应的热工水力分析模型B,安全阀模型包括与安全阀(3)对应的热工水力分析模型C,卸压箱模型包括与卸压箱(4)对应的热工水力分析模型D,水封体管模型包括与水封体管(5)对应的热工水力分析模型E,节块模型包括与控制体节块对应的热工水力分析模型a,阀下游排放管模型包括与阀下游排放管(1)对应的热工水力分析模型A,热工水力分析模型A由n个热工水力分析模型a通过边界连接和数据传递连接的方式连通构成,水封体排放管线模型包括由热工水力分析模型A、热工水力分析模型B、热工水力分析模型C、热工水力分析模型D、热工水力分析模型E通过边界连接和数据传递连接的方式连通构成与水封体排放管线对应的热工水力分析模型Z;
构建阀下游排放管载荷分析模型:以热工水力分析模型a输出的热工水力分析结果信息集a作为输入变量和热工水力分析模型A输出的热工水力分析结果信息集A作为输入变量,以阀下游排放管载荷为解释变量,从而构建起阀下游排放管载荷分析模型;
所述阀下游排放管载荷分析模型是由加速项模型、动能变化模型、压降模型、流体体积力模型求和处理组成;其中加速项模型是指对各个控制体节块内的加速项进行积分处理;动能变化模型是指对整个阀下游排放管的入口和出口的动能变化进行处理;压降模型是指对整个阀下游排放管的入口和出口的压降变化进行处理;流体体积力模型是指对各个控制体节块内的体积力进行积分处理;
水封体的信息包括水封体体积、水封温度分布,安全阀的信息包括阀门开启压力、阀门排量、阀门完全开启时间,阀下游排放管的几何信息包括管道尺寸、管道布置、弯头角度及弯曲半径,卸压箱的几何信息包括卸压箱尺寸,水封体排放管线的内部流体的初始条件和边界条件包括气体种类、初始压力、初始温度。
2.根据权利要求1所述的针对具有阀下游排放管的载荷分析方法,其特征在于,得到随时间实时变化的阀下游排放管载荷值后,将离散的阀下游排放管载荷值拟合成随时间实时变化的载荷波。
3.根据权利要求1所述的针对具有阀下游排放管的载荷分析方法,其特征在于,所述分析时间步长控制在10-4秒的量级。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910280875.4A CN109977598B (zh) | 2019-04-09 | 2019-04-09 | 针对阀下游排放管的载荷分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910280875.4A CN109977598B (zh) | 2019-04-09 | 2019-04-09 | 针对阀下游排放管的载荷分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109977598A CN109977598A (zh) | 2019-07-05 |
CN109977598B true CN109977598B (zh) | 2022-02-18 |
Family
ID=67083654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910280875.4A Active CN109977598B (zh) | 2019-04-09 | 2019-04-09 | 针对阀下游排放管的载荷分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109977598B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111125972B (zh) * | 2019-12-26 | 2021-10-19 | 西安交通大学 | 核电厂破口失水事故水力载荷分析方法 |
CN111524623B (zh) * | 2020-04-30 | 2022-02-22 | 中国核动力研究设计院 | 一种稳压器安全阀定值和排布方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2339984Y (zh) * | 1998-03-10 | 1999-09-22 | 申景泉 | 全量型安全阀在线定压仪 |
CN107895095A (zh) * | 2017-12-15 | 2018-04-10 | 西安交通大学 | 棒束子通道与临界热流密度机理模型的耦合分析方法 |
CN109018287A (zh) * | 2018-08-01 | 2018-12-18 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 用于深海核动力水下平台的稳压器卸压*** |
CN109299536A (zh) * | 2018-09-20 | 2019-02-01 | 西安交通大学 | 一种大型压水堆核电厂稳压器水封形成过程计算方法 |
CN109453559A (zh) * | 2018-11-29 | 2019-03-12 | 合肥科锐特环保工程有限公司 | 一种乙醇洗脱液回收装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7139219B2 (en) * | 2004-02-12 | 2006-11-21 | Tempress Technologies, Inc. | Hydraulic impulse generator and frequency sweep mechanism for borehole applications |
US8196395B2 (en) * | 2009-06-29 | 2012-06-12 | Lightsail Energy, Inc. | Compressed air energy storage system utilizing two-phase flow to facilitate heat exchange |
CN105443857B (zh) * | 2015-12-13 | 2018-06-29 | 北京化工大学 | 一种用于阀门的变阻尼减振降噪装置 |
US10807681B2 (en) * | 2017-05-25 | 2020-10-20 | Sumathi Paturu | Fire-escape devices of the off-shore rigs with emphasis on a detachable island rig |
-
2019
- 2019-04-09 CN CN201910280875.4A patent/CN109977598B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2339984Y (zh) * | 1998-03-10 | 1999-09-22 | 申景泉 | 全量型安全阀在线定压仪 |
CN107895095A (zh) * | 2017-12-15 | 2018-04-10 | 西安交通大学 | 棒束子通道与临界热流密度机理模型的耦合分析方法 |
CN109018287A (zh) * | 2018-08-01 | 2018-12-18 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 用于深海核动力水下平台的稳压器卸压*** |
CN109299536A (zh) * | 2018-09-20 | 2019-02-01 | 西安交通大学 | 一种大型压水堆核电厂稳压器水封形成过程计算方法 |
CN109453559A (zh) * | 2018-11-29 | 2019-03-12 | 合肥科锐特环保工程有限公司 | 一种乙醇洗脱液回收装置 |
Non-Patent Citations (2)
Title |
---|
两相排放载荷分析方法研究;吴丹等;《核动力工程》;20150415;第36卷(第2期);160-164 * |
金飞.AP1000自动卸压***的特点和设计优化.《 中国核科学技术进展报告(第二卷)——中国核学会2011年学术年会论文集第2册(铀矿冶分卷、核能动力分卷(上))》.2011, * |
Also Published As
Publication number | Publication date |
---|---|
CN109977598A (zh) | 2019-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Simonetti et al. | Evaluation of air compressibility effects on the performance of fixed OWC wave energy converters using CFD modelling | |
Zhou et al. | Investigation of hydraulic transients of two entrapped air pockets in a water pipeline | |
Greyvenstein | An implicit method for the analysis of transient flows in pipe networks | |
Zhou et al. | Influence of entrapped air pockets on hydraulic transients in water pipelines | |
Wang et al. | Water hammer simulation using explicit–implicit coupling methods | |
CN110197049B (zh) | 一种基于瞬变反问题的非金属管道泄漏定位方法 | |
CN106844913B (zh) | 一种基于三维cfd的滞留气团热力学特性模拟方法 | |
CN108763683B (zh) | 一种三角函数框架下新weno格式构造方法 | |
CN109918787A (zh) | 基于有限体积法的输水管道内水气两相均质流的模拟方法 | |
CN109977598B (zh) | 针对阀下游排放管的载荷分析方法 | |
Torres et al. | Study of the turbine power output of an oscillating water column device by using a hydrodynamic–Aerodynamic coupled model | |
Torres et al. | A methodology to determine the optimal size of a wells turbine in an oscillating water column device by using coupled hydro-aerodynamic models | |
CN105184395A (zh) | 含余热利用***的火电机组的初参数确定方法 | |
Tomor et al. | Validation of a discrete model for flow distribution in dividing-flow manifolds: Numerical and experimental studies | |
CN111125972B (zh) | 核电厂破口失水事故水力载荷分析方法 | |
Li et al. | 1D-3D coupling investigation of hydraulic transient for power-supply failure of centrifugal pump-pipe system | |
Zolfaghary Azizi et al. | Comparison of numerical methods for two-fluid model for gas–liquid transient flow regime and its application in slug modeling initiation | |
Adamkowski et al. | Improved discharge measurement using the pressure-time method in a hydropower plant curved penstock | |
Simonetti | Optimization of oscillating water column wave energy converters-a numerical study | |
CN113657049A (zh) | 一种池式钠冷快堆主冷却剂***传热与流动快速仿真方法 | |
Abouri et al. | A stable Fluid-Structure-Interaction algorithm: Application to industrial problems | |
Mahood et al. | Analytical and numerical investigation of transient gas blow down | |
Adamkowski et al. | A new approach to using the classic Gibson method to measure discharge | |
Xu et al. | Theoretical investigation on the throttle pressure reducing valve through CFD simulation and validating experiments | |
Greyvenstein et al. | Dynamic modelling of heat, mass and momentum transfer in the Pebble Bed Modular Reactor |
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 |