CN113268824B - 基于温控阀内流场的流体回路数值模拟方法、系统及介质 - Google Patents

基于温控阀内流场的流体回路数值模拟方法、系统及介质 Download PDF

Info

Publication number
CN113268824B
CN113268824B CN202110442932.1A CN202110442932A CN113268824B CN 113268824 B CN113268824 B CN 113268824B CN 202110442932 A CN202110442932 A CN 202110442932A CN 113268824 B CN113268824 B CN 113268824B
Authority
CN
China
Prior art keywords
temperature control
temperature
control valve
heat
flow field
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
CN202110442932.1A
Other languages
English (en)
Other versions
CN113268824A (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.)
Shanghai Institute of Satellite Engineering
Original Assignee
Shanghai Institute of Satellite Engineering
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 Shanghai Institute of Satellite Engineering filed Critical Shanghai Institute of Satellite Engineering
Priority to CN202110442932.1A priority Critical patent/CN113268824B/zh
Publication of CN113268824A publication Critical patent/CN113268824A/zh
Application granted granted Critical
Publication of CN113268824B publication Critical patent/CN113268824B/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/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/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft 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
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • 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

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)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Control Of Temperature (AREA)

Abstract

本发明提供了一种基于温控阀内流场的流体回路数值模拟方法、系统及介质，包括：对温控阀进行参数化建模；抽取温控阀腔内的流体域；对三维流场进行网格划分，并进行网格相关性验证；建立驱动泵、热沉、辐射器和管路模型，实时计算管路内各个节点的温度、压力和流速值；编译联动Fluent的UDF和Matlab/Simulinkd的M函数，实现不同温控阀开度下流量和压力的数据交换；迭代计算，调试PID参数直到控温点温度稳定在预设范围以内，根据控温效果选择合适的控温点；对仿真结果进行可视化分析。该方法能够精确模拟流体回路温控阀的实际工作流阻，得到更加合适的控温方案和PID控制参数，为温控阀的多参数设计提供平台指导。

Description

基于温控阀内流场的流体回路数值模拟方法、系统及介质
技术领域
本发明涉及航天器热控技术领域，具体地，涉及一种基于温控阀内流场的流体回路数值模拟方法、系统及介质。
背景技术
流体回路系统作为航天器主动热控技术中的一种重要方式，利用机械泵驱动工质强制对流循环来进行热量控制。由于其换热能力强、便于结构布局、易于组织内部换热、安全可靠等特点，在国内外的许多航天器上都得到了广泛的应用。流体回路一般由内外两个回路组成，温控阀通过调节两个回路的工质流量，实现航天器内部热量转移的控制，满足航天器热控目标多样化的需求。
国内外一般使用Flowmaster和Matlab等软件仿真模拟流体回路运行。地面试验发现，温控阀与流体回路的流阻系数存在不匹配的现象，温控阀的内流场结构与流体回路控制策略息息相关。原因在于流体回路驱动泵流量和扬程存在一定关系，温控阀开关运动时流体回路流阻发生变化，导致驱动泵输出流量发生变化。因此，目前的单软件仿真得到的流阻系数与实际运行差别较大，导致温控阀控制能力受限，甚至控温策略失效。
专利文献CN101052469A(申请号：CN200580037969.5)公开了一种生物化学领域微流体槽道盒定时协议的方法和结构创新，而本发明提供了航天器热控制领域流体回路仿真技术。
专利文献CN104375531A(申请号：CN201410602694.6)公开了一种适应单相流体回路的组合温度控制系统及方法。本发明与之相比的创新在于提出了新型交互式流体回路数值模拟方法，搭建了基于温控阀的流体回路多参数仿真平台。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种基于温控阀内流场的流体回路数值模拟方法、系统及介质。
根据本发明提供的基于温控阀内流场的流体回路数值模拟方法，包括：
步骤1：使用三维设计软件对温控阀进行参数化建模；
步骤2：抽取温控阀腔内的流体域；
步骤3：对三维流场进行网格划分，并进行网格相关性验证；
步骤4：建立驱动泵、热沉、辐射器和管路模型，实时计算管路内各个节点的温度、压力和流速值；
步骤5：编译联动仿真分析软件Fluent的UDF和Matlab/Simulinkd的M函数，实现不同温控阀开度下流量和压力的数据交换；
步骤6：迭代计算，调试PID参数直到控温点温度稳定在预设范围以内，根据控温效果选择合适的控温点；
步骤7：对仿真结果进行可视化分析。
优选的，所述步骤1包括：
步骤1.1：根据用户需求选择型谱化温控阀产品；
步骤1.2：对温控阀内流道结构进行参数化建模，选取出口阀腔宽度、进口阀腔宽度、进出口阀腔直径、V型槽槽深和V型槽开口角作为滑阀内流道建模的设计参数。
优选的，所述步骤2包括：
步骤2.1：将温控阀模型保存为.stp等文件格式，导入Spaceclaim软件中；
步骤2.2：选取温控阀的三个进出管边缘作为抽取的边界，使用体积抽取的命令抽取温控阀内流场；
步骤2.3：检查内流场模型是否存在缺失、畸形现象并修复，保留预设长度的温控阀进出管道，防止Fluent运算时发生回流。
优选的，所述步骤3包括：
步骤3.1：将温控阀内流场模型导入软件ICEM中，定义计算域的三个进出口表面和壁，采用四面体非结构化网格技术，对狭小射流区域进行局部加密；
步骤3.2：分别选取不同尺寸的网格进行仿真计算，经过对比验证网格无关性，确保计算结果的可靠性。
优选的，所述步骤4包括：
步骤4.1：工质首先通过热传导和对流换热将热量传输至蛇形管壁，管壁再通过热传导将热量传输至面板，最终由辐射面板通过辐射换热的方式将热量排散至空间环境；辐射器的动态特性为：
Figure BDA0003035654750000031
Figure BDA0003035654750000032
蛇形管路的动态特性为：
Figure BDA0003035654750000033
Figure BDA0003035654750000034
辐射板A、B进出口流体工质的温度关系为：
GAcfA(TA1-TA2)＝hA32AA32(TfA-θA2) (5)
GBcfB(TB1-TB2)＝hB32AB32(TfB-θB2) (6)
Figure BDA0003035654750000035
Figure BDA0003035654750000036
上式中mA1、mB1、cA1、cB1、θA1、θB1、AA1、AB1分别为辐射板A、B的质量、热容、温度和面积；αA、αB、εA、εB分别为辐射板A、B的吸收率和发射率；qA、qB分别为辐射板A、B的外热流密度；hA21、hB21、AA21、AB21分别为辐射板A、B和对应蛇形管路之间的等效换热系数和接触面积；θ0为空间环境温度，取4K；mA2、mB2、cA2、cB2、θA2、θB2、AA2、AB2分别为蛇形管路A、B的质量、热容、温度和面积；hA32、hB32、AA32、AB32为流体工质A、B与对应蛇形管路A、B的等效换热系数与换热面积；GA、GB、cfA、cfB分别为蛇形管路A、B内流体工质的质量流量和热容；TA1、TB1、TA2、TB2分别为辐射板A、B进出口处工质的温度；TfA、TfB为辐射板A、B内工质的平均温度；
步骤4.2：热沉是流体回路系统中负责收集载荷热耗的部件，热沉表面温度相同，被控对象的动态特性为：
Figure BDA0003035654750000037
热沉的动态特性方程为：
Figure BDA0003035654750000038
热沉进出口流体工质温度的关系如下：
h12A12(θ2-Tf1)＝Gfcf(T2-T1) (11)
上式中m1、c1、θ1分别为被控对象的质量、热容和温度；m2、c2、θ2为热沉的质量、热容和温度；Gf和cf为流体工质的质量流量和热容；T1为热沉进口工质的温度；T2为热沉出口工质的温度；Tf为热沉内工质的平均温度，且
Figure BDA0003035654750000041
Q1为被控对象的热耗；h11和A11为热沉与被控对象的等效换热系数与换热面积；h12和A12为热沉与流体工质的等效换热系数与换热面积；
步骤4.3：管路包覆多层隔热组件，忽略管路与周围环境的换热，此时管路的传递函数为一阶惯性加纯延时环节，惯性环节的时间常数Tpipe和延迟时间τ通过实验曲线整定，表达式为：
Figure BDA0003035654750000042
其中，s为复变量；
步骤4.4：驱动泵与流体回路的热量关系式满足：
Figure BDA0003035654750000043
Qpump＝(1-η)UI (14)
其中，Qpump为驱动泵产生的热量，η为驱动泵的工作效率，U为工作电压，I为工作电流。
优选的，所述步骤5包括：
步骤5.1：启动Simulink模型，计算得到温控阀进口处流量和表压，用M函数导入Fluent将作为内流场计算的初始边界条件；
步骤5.2：将网格化模型导入Fluent中，使用UDF将Simulink反馈的温控阀开度值传递到Fluent，计算后输出温控阀出口处流速和表压，该结果作为Simulink下一轮计算的输入。
优选的，所述步骤6包括：
步骤6.1：以温控阀进口处工质温度作为控温点，调试PID参数，输出控温点温度变化曲线；
步骤6.2：以目标单机表面作为控温点，调试PID参数，输出控温单机温度变化曲线；
步骤6.3：根据动态特性曲线比较控温效果，选择合适的控温点。
优选的，所述步骤7包括：
步骤7.1：Fluent输出温控阀内流场压力流线图，为结构优化提供理论依据；
步骤7.2：Simulink输出目标温度-温控阀开度曲线，为控温策略优化提供参考依据。
根据本发明提供的基于温控阀内流场的流体回路数值模拟系统，包括：
模块M1：使用三维设计软件对温控阀进行参数化建模；
模块M2：抽取温控阀腔内的流体域；
模块M3：对三维流场进行网格划分，并进行网格相关性验证；
模块M4：建立驱动泵、热沉、辐射器和管路模型，实时计算管路内各个节点的温度、压力和流速值；
模块M5：编译联动仿真分析软件Fluent的UDF和Matlab/Simulinkd的M函数，实现不同温控阀开度下流量和压力的数据交换；
模块M6：迭代计算，调试PID参数直到控温点温度稳定在预设范围以内，根据控温效果选择合适的控温点；
模块M7：对仿真结果进行可视化分析。
根据本发明提供的一种存储有计算机程序的计算机可读存储介质，所述计算机程序被处理器执行时实现上述的方法的步骤。
与现有技术相比，本发明具有如下的有益效果：
本发明方法能够精确模拟流体回路温控阀的实际工作流阻，得到更加合适的控温方案和PID控制参数，为温控阀的多参数设计提供平台指导。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为基于温控阀内流场的流体回路数值模拟方法流程图；
图2为Spaceclaim环境下温控阀腔内流场域结构图；
图3为流体回路辐射器Simulink模型图；
图4为流体回路热沉Simulink模型图；
图5为流体回路管路Simulink模型图；
图6为流体回路驱动泵Simulink模型图；
图7为流体回路控温算法Simulink模型图；
图8为Simulink环境下添加热扰动后温控阀开度和被控对象温度变化曲线图；
图9为驱动泵工作曲线图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
实施例1：
如图1，根据本发明提供的基于温控阀内流场的流体回路数值模拟方法，包括：步骤1：使用三维设计软件对温控阀进行参数化建模；步骤2：抽取温控阀腔内的流体域；步骤3：对三维流场进行网格划分，并进行网格相关性验证；步骤4：建立驱动泵、热沉、辐射器和管路模型，实时计算管路内各个节点的温度、压力和流速值；步骤5：编译联动仿真分析软件Fluent的UDF和Matlab/Simulinkd的M函数，实现不同温控阀开度下流量和压力的数据交换；步骤6：迭代计算，调试PID参数直到控温点温度稳定在预设范围以内，根据控温效果选择合适的控温点；步骤7：对仿真结果进行可视化分析。
所述步骤1包括：步骤1.1：根据用户需求选择型谱化温控阀产品；步骤1.2：对温控阀内流道结构进行参数化建模，选取出口阀腔宽度、进口阀腔宽度、进出口阀腔直径、V型槽槽深和V型槽开口角作为滑阀内流道建模的设计参数。
所述步骤2包括：步骤2.1：将温控阀模型保存为.stp等文件格式，导入Spaceclaim软件中；步骤2.2：选取温控阀的三个进出管边缘作为抽取的边界，使用体积抽取的命令抽取温控阀内流场；步骤2.3：检查内流场模型是否存在缺失、畸形现象并修复，保留预设长度的温控阀进出管道，防止Fluent运算时发生回流。
所述步骤3包括：步骤3.1：将温控阀内流场模型导入软件ICEM中，定义计算域的三个进出口表面和壁，采用四面体非结构化网格技术，对狭小射流区域进行局部加密；步骤3.2：分别选取不同尺寸的网格进行仿真计算，经过对比验证网格无关性，确保计算结果的可靠性。
所述步骤4包括：步骤4.1：工质首先通过热传导和对流换热将热量传输至蛇形管壁，管壁再通过热传导将热量传输至面板，最终由辐射面板通过辐射换热的方式将热量排散至空间环境；辐射器的动态特性为：
Figure BDA0003035654750000071
Figure BDA0003035654750000072
蛇形管路的动态特性为：
Figure BDA0003035654750000073
Figure BDA0003035654750000074
辐射板A、B进出口流体工质的温度关系为：GAcfA(TA1-TA2)＝hA32AA32(TfA-θA2) (5)
GBcfB(TB1-TB2)＝hB32AB32(TfB-θB2) (6)
Figure BDA0003035654750000075
Figure BDA0003035654750000076
上式中mA1、mB1、cA1、cB1、θA1、θB1、AA1、AB1分别为辐射板A、B的质量、热容、温度和面积；αA、αB、εA、εB分别为辐射板A、B的吸收率和发射率；qA、qB分别为辐射板A、B的外热流密度；hA21、hB21、AA21、AB21分别为辐射板A、B和对应蛇形管路之间的等效换热系数和接触面积；θ0为空间环境温度，取4K；mA2、mB2、cA2、cB2、θA2、θB2、AA2、AB2分别为蛇形管路A、B的质量、热容、温度和面积；hA32、hB32、AA32、AB32为流体工质A、B与对应蛇形管路A、B的等效换热系数与换热面积；GA、GB、cfA、cfB分别为蛇形管路A、B内流体工质的质量流量和热容；TA1、TB1、TA2、TB2分别为辐射板A、B进出口处工质的温度；TfA、TfB为辐射板A、B内工质的平均温度；
步骤4.2：热沉是流体回路系统中负责收集载荷热耗的部件，热沉表面温度相同，被控对象的动态特性为：
Figure BDA0003035654750000077
热沉的动态特性方程为：
Figure BDA0003035654750000078
热沉进出口流体工质温度的关系如下：h12A12(θ2-Tf1)＝Gfcf(T2-T1) (11)
上式中m1、c1、θ1分别为被控对象的质量、热容和温度；m2、c2、θ2为热沉的质量、热容和温度；Gf和cf为流体工质的质量流量和热容；T1为热沉进口工质的温度；T2为热沉出口工质的温度；Tf为热沉内工质的平均温度，且
Figure BDA0003035654750000079
Q1为被控对象的热耗；h11和A11为热沉与被控对象的等效换热系数与换热面积；h12和A12为热沉与流体工质的等效换热系数与换热面积；
步骤4.3：管路包覆多层隔热组件，忽略管路与周围环境的换热，此时管路的传递函数为一阶惯性加纯延时环节，惯性环节的时间常数Tpipe和延迟时间τ通过实验曲线整定，表达式为：
Figure BDA0003035654750000081
其中，s为复变量；
步骤4.4：驱动泵与流体回路的热量关系式满足：
Figure BDA0003035654750000082
Qpump＝(1-η)UI (14)
其中，Qpump为驱动泵产生的热量，η为驱动泵的工作效率，U为工作电压，I为工作电流。
所述步骤5包括：步骤5.1：启动Simulink模型，计算得到温控阀进口处流量和表压，用M函数导入Fluent将作为内流场计算的初始边界条件；步骤5.2：将网格化模型导入Fluent中，使用UDF将Simulink反馈的温控阀开度值传递到Fluent，计算后输出温控阀出口处流速和表压，该结果作为Simulink下一轮计算的输入。
所述步骤6包括：步骤6.1：以温控阀进口处工质温度作为控温点，调试PID参数，输出控温点温度变化曲线；步骤6.2：以目标单机表面作为控温点，调试PID参数，输出控温单机温度变化曲线；步骤6.3：根据动态特性曲线比较控温效果，选择合适的控温点。
所述步骤7包括：步骤7.1：Fluent输出温控阀内流场压力流线图，为结构优化提供理论依据；步骤7.2：Simulink输出目标温度-温控阀开度曲线，为控温策略优化提供参考依据。
根据本发明提供的基于温控阀内流场的流体回路数值模拟系统，包括：模块M1：使用三维设计软件对温控阀进行参数化建模；模块M2：抽取温控阀腔内的流体域；模块M3：对三维流场进行网格划分，并进行网格相关性验证；模块M4：建立驱动泵、热沉、辐射器和管路模型，实时计算管路内各个节点的温度、压力和流速值；模块M5：编译联动仿真分析软件Fluent的UDF和Matlab/Simulinkd的M函数，实现不同温控阀开度下流量和压力的数据交换；模块M6：迭代计算，调试PID参数直到控温点温度稳定在预设范围以内，根据控温效果选择合适的控温点；模块M7：对仿真结果进行可视化分析。
根据本发明提供的一种存储有计算机程序的计算机可读存储介质，所述计算机程序被处理器执行时实现上述的方法的步骤。
实施例2：
实施例2为实施例1的优选例。
本发明提供了基于参数化温控阀建模平台的流体回路数值模拟方法，具体实施步骤如下：
步骤1、使用Proe等三维设计软件对温控阀进行参数化建模。根据某温控阀结构形式，出口阀腔宽度a＝6mm、进口阀腔宽度b＝7mm、进出口阀腔直径c＝18mm、V型槽槽深d＝2mm和V型槽开口角e＝130°。
步骤2、使用Spaceclaim软件抽取温控阀腔内的流体域，见图2。将温控阀模型保存为.stp等文件格式，导入Spaceclaim软件中。选取温控阀的三个进出管边缘作为抽取的边界，使用“体积抽取”命令抽取温控阀内流场。检查内流场模型是否存在缺失、畸形现象并修复，保留一定的温控阀进出管道长度，防止Fluent运算时发生回流。
步骤3、使用ICEM等网格划分软件对三维流场进行网格划分，并进行网格相关性验证。网格细化后误差不超过2％即可认为网格无关。
步骤4、使用Matlab/Simulink建立驱动泵、热沉、辐射器和管路模型，实时计算管路内各个节点的温度、压力和流速值。图3为流体回路辐射器Simulink模型；图4为流体回路热沉Simulink模型；图5为流体回路管路Simulink模型；图6为流体回路驱动泵Simulink模型；图7为流体回路控温算法Simulink模型。
步骤5、编译联动Fluent的UDF和Matlab/Simulinkd的M函数，实现不同温控阀开度下流量和压力的数据交换。得到温控阀内流场压力流线图。添加外热流扰动后可得到图8所示的温控阀开度和被控对象温度变化曲线。基于对图8的分析，开展温控阀结构参数和流体回路控温策略的优化设计。
驱动泵输出参数扬程H和体积流量Q间满足一定的函数关系，其函数关系H＝F(Q)见图9。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以，本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件，而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构；也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。

Claims (9)

1.一种基于温控阀内流场的流体回路数值模拟方法，其特征在于，包括：
步骤1：使用三维设计软件对温控阀进行参数化建模；
步骤2：抽取温控阀腔内的流体域；
步骤3：对三维流场进行网格划分，并进行网格相关性验证；
步骤4：建立驱动泵、热沉、辐射器和管路模型，实时计算管路内各个节点的温度、压力和流速值；
步骤5：编译联动仿真分析软件Fluent的UDF和Matlab/Simulink的M函数，实现不同温控阀开度下流量和压力的数据交换；
步骤6：迭代计算，调试PID参数直到控温点温度稳定在预设范围以内，根据控温效果选择合适的控温点；
步骤7：对仿真结果进行可视化分析；
所述步骤4包括：
步骤4.1：工质首先通过热传导和对流换热将热量传输至蛇形管壁，管壁再通过热传导将热量传输至面板，最终由辐射面板通过辐射换热的方式将热量排散至空间环境；辐射器的动态特性为：
Figure FDA0004058404860000011
Figure FDA0004058404860000012
蛇形管路的动态特性为：
Figure FDA0004058404860000013
Figure FDA0004058404860000014
辐射板A、B进出口流体工质的温度关系为：
GAcfA(TA1-TA2)＝hA32AA32(TfA-θA2)(5)
GBcfB(TB1-TB2)＝hB32AB32(TfB-θB2)(6)
Figure FDA0004058404860000015
Figure FDA0004058404860000021
上式中mA1、mB1、cA1、cB1、θA1、θB1、AA1、AB1分别为辐射板A、B的质量、热容、温度和面积；αA、αB、εA、εB分别为辐射板A、B的吸收率和发射率；qA、qB分别为辐射板A、B的外热流密度；hA21、hB21、AA21、AB21分别为辐射板A、B和对应蛇形管路之间的等效换热系数和接触面积；θ0为空间环境温度，取4K；mA2、mB2、cA2、cB2、θA2、θB2、AA2、AB2分别为蛇形管路A、B的质量、热容、温度和面积；hA32、hB32、AA32、AB32为流体工质A、B与对应蛇形管路A、B的等效换热系数与换热面积；GA、GB、cfA、cfB分别为蛇形管路A、B内流体工质的质量流量和热容；TA1、TB1、TA2、TB2分别为辐射板A、B进出口处工质的温度；TfA、TfB为辐射板A、B内工质的平均温度；
步骤4.2：热沉是流体回路系统中负责收集载荷热耗的部件，热沉表面温度相同，被控对象的动态特性为：
Figure FDA0004058404860000022
热沉的动态特性方程为：
Figure FDA0004058404860000023
热沉进出口流体工质温度的关系如下：
h12A12(θ2-Tf1)＝Gfcf(T2-T1)(11)
上式中m1、c1、θ1分别为被控对象的质量、热容和温度；m2、c2、θ2为热沉的质量、热容和温度；Gf和cf为流体工质的质量流量和热容；T1为热沉进口工质的温度；T2为热沉出口工质的温度；Tf为热沉内工质的平均温度，且
Figure FDA0004058404860000024
Q1为被控对象的热耗；h11和A11为热沉与被控对象的等效换热系数与换热面积；h12和A12为热沉与流体工质的等效换热系数与换热面积；
步骤4.3：管路包覆多层隔热组件，忽略管路与周围环境的换热，此时管路的传递函数为一阶惯性加纯延时环节，惯性环节的时间常数Tpipe和延迟时间τ通过实验曲线整定，表达式为：
Figure FDA0004058404860000025
其中，s为复变量；
步骤4.4：驱动泵与流体回路的热量关系式满足：
Figure FDA0004058404860000031
Qpump＝(1-η)UI(14)
其中，Qpump为驱动泵产生的热量，η为驱动泵的工作效率，U为工作电压，I为工作电流。
2.根据权利要求1所述的基于温控阀内流场的流体回路数值模拟方法，其特征在于，所述步骤1包括：
步骤1.1：根据用户需求选择型谱化温控阀产品；
步骤1.2：对温控阀内流道结构进行参数化建模，选取出口阀腔宽度、进口阀腔宽度、进出口阀腔直径、V型槽槽深和V型槽开口角作为滑阀内流道建模的设计参数。
3.根据权利要求1所述的基于温控阀内流场的流体回路数值模拟方法，其特征在于，所述步骤2包括：
步骤2.1：将温控阀模型保存为.stp等文件格式，导入Spaceclaim软件中；
步骤2.2：选取温控阀的三个进出管边缘作为抽取的边界，使用体积抽取的命令抽取温控阀内流场；
步骤2.3：检查内流场模型是否存在缺失、畸形现象并修复，保留预设长度的温控阀进出管道，防止Fluent运算时发生回流。
4.根据权利要求1所述的基于温控阀内流场的流体回路数值模拟方法，其特征在于，所述步骤3包括：
步骤3.1：将温控阀内流场模型导入软件ICEM中，定义计算域的三个进出口表面和壁，采用四面体非结构化网格技术，对狭小射流区域进行局部加密；
步骤3.2：分别选取不同尺寸的网格进行仿真计算，经过对比验证网格无关性，确保计算结果的可靠性。
5.根据权利要求1所述的基于温控阀内流场的流体回路数值模拟方法，其特征在于，所述步骤5包括：
步骤5.1：启动Simulink模型，计算得到温控阀进口处流量和表压，用M函数导入Fluent将作为内流场计算的初始边界条件；
步骤5.2：将网格化模型导入Fluent中，使用UDF将Simulink反馈的温控阀开度值传递到Fluent，计算后输出温控阀出口处流速和表压，该结果作为Simulink下一轮计算的输入。
6.根据权利要求1所述的基于温控阀内流场的流体回路数值模拟方法，其特征在于，所述步骤6包括：
步骤6.1：以温控阀进口处工质温度作为控温点，调试PID参数，输出控温点温度变化曲线；
步骤6.2：以目标单机表面作为控温点，调试PID参数，输出控温单机温度变化曲线；
步骤6.3：根据动态特性曲线比较控温效果，选择合适的控温点。
7.根据权利要求1所述的基于温控阀内流场的流体回路数值模拟方法，其特征在于，所述步骤7包括：
步骤7.1：Fluent输出温控阀内流场压力流线图，为结构优化提供理论依据；
步骤7.2：Simulink输出目标温度-温控阀开度曲线，为控温策略优化提供参考依据。
8.一种基于温控阀内流场的流体回路数值模拟系统，其特征在于，采用权利要求1-7中任一项所述的基于温控阀内流场的流体回路数值模拟方法去执行，包括：
模块M1：使用三维设计软件对温控阀进行参数化建模；
模块M2：抽取温控阀腔内的流体域；
模块M3：对三维流场进行网格划分，并进行网格相关性验证；
模块M4：建立驱动泵、热沉、辐射器和管路模型，实时计算管路内各个节点的温度、压力和流速值；
模块M5：编译联动仿真分析软件Fluent的UDF和Matlab/Simulink的M函数，实现不同温控阀开度下流量和压力的数据交换；
模块M6：迭代计算，调试PID参数直到控温点温度稳定在预设范围以内，根据控温效果选择合适的控温点；
模块M7：对仿真结果进行可视化分析。
9.一种存储有计算机程序的计算机可读存储介质，其特征在于，所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。
CN202110442932.1A 2021-04-23 2021-04-23 基于温控阀内流场的流体回路数值模拟方法、系统及介质 Active CN113268824B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110442932.1A CN113268824B (zh) 2021-04-23 2021-04-23 基于温控阀内流场的流体回路数值模拟方法、系统及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110442932.1A CN113268824B (zh) 2021-04-23 2021-04-23 基于温控阀内流场的流体回路数值模拟方法、系统及介质

Publications (2)

Publication Number Publication Date
CN113268824A CN113268824A (zh) 2021-08-17
CN113268824B true CN113268824B (zh) 2023-07-11

Family

ID=77229259

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110442932.1A Active CN113268824B (zh) 2021-04-23 2021-04-23 基于温控阀内流场的流体回路数值模拟方法、系统及介质

Country Status (1)

Country Link
CN (1) CN113268824B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104932565A (zh) * 2015-06-04 2015-09-23 中国空气动力研究与发展中心高速空气动力研究所 一种高精度暂冲型引射式跨声速风洞流场控制结构
CN106527519A (zh) * 2016-12-07 2017-03-22 中国空气动力研究与发展中心高速空气动力研究所 一种大型直流暂冲式超声速风洞变速压流场控制方法
CN109002617A (zh) * 2018-07-23 2018-12-14 安徽江淮纳威司达柴油发动机有限公司 一种考虑热辐射的发动机排气歧管的温度场分析方法
WO2019152907A1 (en) * 2018-02-02 2019-08-08 Schlumberger Technology Corporation Flexible gradient-based reservoir simulation optimization
CN112528572A (zh) * 2020-11-16 2021-03-19 天华化工机械及自动化研究设计院有限公司 一种基于overset模型的低温碳化炉丝束加热过程三维模拟方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080275677A1 (en) * 2007-03-19 2008-11-06 Optimal Solutions Software, Llc System, methods, and computer readable media, for product design using coupled computer aided engineering models
CN203058203U (zh) * 2012-12-28 2013-07-17 北京天达京丰技术开发有限公司 一种蒸发汽回收循环利用的设备系统
US9400256B2 (en) * 2013-12-06 2016-07-26 Rolls-Royce Corporation Thermographic inspection techniques
CN106081126A (zh) * 2016-06-13 2016-11-09 王晨 仿生蜂窝状主动安全逃生舱嵌入航空飞行器的应用及设计
CN107016203A (zh) * 2017-04-14 2017-08-04 江苏大学 一种光伏水泵内部流动的数值模拟方法
CN110908414A (zh) * 2019-10-28 2020-03-24 沪东重机有限公司 一种管道汇流液体温度的控制系统及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104932565A (zh) * 2015-06-04 2015-09-23 中国空气动力研究与发展中心高速空气动力研究所 一种高精度暂冲型引射式跨声速风洞流场控制结构
CN106527519A (zh) * 2016-12-07 2017-03-22 中国空气动力研究与发展中心高速空气动力研究所 一种大型直流暂冲式超声速风洞变速压流场控制方法
WO2019152907A1 (en) * 2018-02-02 2019-08-08 Schlumberger Technology Corporation Flexible gradient-based reservoir simulation optimization
CN109002617A (zh) * 2018-07-23 2018-12-14 安徽江淮纳威司达柴油发动机有限公司 一种考虑热辐射的发动机排气歧管的温度场分析方法
CN112528572A (zh) * 2020-11-16 2021-03-19 天华化工机械及自动化研究设计院有限公司 一种基于overset模型的低温碳化炉丝束加热过程三维模拟方法

Also Published As

Publication number Publication date
CN113268824A (zh) 2021-08-17

Similar Documents

Publication Publication Date Title
Bonilla et al. Switching moving boundary models for two-phase flow evaporators and condensers
Tian et al. Coupling fast fluid dynamics and multizone airflow models in Modelica Buildings library to simulate the dynamics of HVAC systems
Li et al. Recent advances in dynamic modeling of HVAC equipment. Part 2: Modelica-based modeling
Matsushima et al. Dynamic simulation program with object-oriented formulation for absorption chillers (modelling, verification, and application to triple-effect absorption chiller)
Zupančič et al. Computer-aided physical multi-domain modelling: Some experiences from education and industrial applications
Hirsch et al. An efficient numerical solution method for detailed modelling of large 5th generation district heating and cooling networks
Picard Modeling, optimal control and HVAC design of large buildings using ground source heat pump systems
CN113268824B (zh) 基于温控阀内流场的流体回路数值模拟方法、系统及介质
Chandan Modeling and control of hydronic building HVAC systems
US11580281B2 (en) System and method for designing heating, ventilating, and air-conditioning (HVAC) systems
Jorissen et al. Integrated Modelica Model and Model Predictive Control of a Terraced House Using IDEAS.
Zhou et al. Dynamic modeling and thermodynamic analysis of lithium bromide absorption refrigeration system using Modelica
Morini et al. The simulation of transients in thermal plant. Part I: Mathematical model
Elarga et al. Numerical investigation of a CO2 cooling system connected to Spawn-of-energy-plus thermal zones
Filonenko et al. Modelica implementation of phase change material ventilation unit
Cui et al. Simulation modelling and design optimization of vertical ground heat exchanger-GEOSTAR program
Bonilla et al. Mean densities in dynamic mathematical two-phase flow models
CN112182905B (zh) 一种用于综合能源系统的供热管网仿真方法和装置
CN116227113A (zh) 热力系统仿真的求解方法、装置、电子设备及介质
Uren et al. State space model extraction of thermohydraulic systems–Part I: A linear graph approach
Sarmiento et al. A numerical method for shell and thermosyphon heat exchanger analysis
Chen et al. Smoothing techniques in dynamic building system simulation
Burhenne et al. Integration of HVAC models into a hygrothermal whole building simulation tool
CN111222230A (zh) 热泵类电采暖设备的建模处理方法及装置
Berrabah et al. Towards a quadrupole-based method for buildings simulation: Validation with ASHRAE 140 standard

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