CN114239429B - 船舶冷凝器壳侧热力性能数值模拟方法 - Google Patents
船舶冷凝器壳侧热力性能数值模拟方法 Download PDFInfo
- Publication number
- CN114239429B CN114239429B CN202111467409.0A CN202111467409A CN114239429B CN 114239429 B CN114239429 B CN 114239429B CN 202111467409 A CN202111467409 A CN 202111467409A CN 114239429 B CN114239429 B CN 114239429B
- Authority
- CN
- China
- Prior art keywords
- condenser
- unit
- outlet
- value
- flow
- 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/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
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明的目的在于提供船舶冷凝器壳侧热力性能数值模拟方法,包括如下步骤:将多结构冷凝器管侧分解为仅具有单进口、单出口特点的子单元,设定多结构冷凝器壳侧各工质在入口和出口处的流量焓值的初始值,求出第n个单元对应的管侧工质出口焓值以及管内平均速度和温度,离散冷凝器各单元的质量、动量、空气组分及温度传输方程,并生成系数矩阵,数值求解冷凝器各单元通道的系数矩阵,获得各单元us、Ts和Y值,求解Gsc,Gac和Hsc,判断冷凝器各单元出口处的Gsc,Gac和Hsc与上一步计算值或初始值的最大偏差是否小于收敛值,根据收敛情况输出结果。本发明可将船舶冷凝器复杂结构离散化为一系列仅具有单入口、出口结构的管壳式冷凝器子单元。
Description
技术领域
本发明涉及的是一种船舶模拟方法,具体地说是船舶二回路冷凝器模拟方法。
背景技术
目前已有的冷凝器壳侧热力性能数值模拟方法主要采用结构网格划分方式,数值计算方法采用基于有限差分法或有限体积法的二维计算流体力学模型,将复杂壳侧结构简化为多孔介质,采用质量源项描述冷凝过程,采用阻力源项修正复杂管结构对流动的影响,计算数据后处理采用文本导出方式。该方法对冷凝器边界条件往往采用定温、定热流或绝热边界条件,仅能计算简单结构的管壳、表面式冷凝器壳侧热力性能,网格前处理方法简单,后处理模块缺失,难以满足对船舶二回路系统多结构冷凝器的管束布置及壳侧流场进行精确、快速的热力性能模拟要求。
发明内容
本发明的目的在于提供一种基于流路质量、能量流分解和二维数值模拟方法相结合的船舶冷凝器壳侧热力性能数值模拟方法。
本发明的目的是这样实现的:
本发明船舶冷凝器壳侧热力性能数值模拟方法,其特征是:
(1)将多结构冷凝器管侧按照工质流动关系,分解为仅具有单进口、单出口特点的子单元,第n号子单元入口、出口变量包括循环水质量流量Gwr,n/Gwc,n和热流量Hwr,n/Hwc,n,其中下标r表示入口,下标c表示出口,n表示单元编号;
(2)设定多结构冷凝器壳侧各工质在入口和出口处的流量Gsr,n/Gsc,n和Hsn焓值的初始值;
(3)根据管、壳间换热能量守恒原则,第n个单元的管壳式换热器满足以下关系:
Gwr,nHwr,n–Gwc,nHwc,n=(Gsr-Gsc)nHsn
求出第n个单元对应的管侧工质出口焓值,以及管内平均速度uw和温度Tw,各子单元管侧之间按照工质流动流程具有串联、并联或独立的关系,当子单元1和子单元2具有串联关系时,两者间的变量满足:
Gwc,1=Gwr,2;Hwc,1=Hwr,2
当子单元1和子单元2具有并联关系时,各自入口、出口关系独立,但两者流量分配满足压降相等原则,即
△P1=△P2
且与各自上游、下游单元具有对应的流量关系,其Gw和Hw可依据质量守恒和能量守恒原理确定;
为冷凝器壳侧每个单元通道划分二维计算网格,根据壳侧蒸汽质量流量Gs、空气质量流量Ga、热流量Hs值、壁面边界及管侧传热特性设定各单元内部网格的us、Ts、空气质量分数Y变量初始值及边界条件;
(4)基于网格、边界条件和us/Ts/Y初始值信息,采用数值方法离散冷凝器各单元的质量、动量、空气组分及温度传输方程,并生成系数矩阵;其中,有限元离散方法采用大数法,对流项采用P1格式离散,扩散项采用中心节点法的P2格式离散,并采用SUPG方法考虑对流项占优数学问题;
(5)数值求解冷凝器各单元通道的系数矩阵,获得各单元us、Ts和Y值,并根据出口处的us、Ts、Y值求解Gsc,Gac和Hsc;
(6)判断冷凝器各单元出口处的Gsc,Gac和Hsc与上一步计算值或初始值的最大偏差是否小于收敛值;
(7)如果步骤(6)不收敛,则将第(6)步获得的冷凝器各单元的Gsc,Gac和Hsc作为已知量,重复步骤(3)-(6),直至收敛;
(8)如果步骤(6)收敛,则输出结果。
本发明的优势在于:本发明可将船舶冷凝器复杂结构离散化为一系列仅具有单入口、出口结构的管壳式冷凝器子单元,并根据工质流动的串、并联关系及压降(来自CFD计算结果)确定冷凝器各单元的CFD边界条件;采用迭代的方法取得最终收敛结果,通过各部件之间的关系确定各单元传热部件的边界条件,进而对各单元部件采用CFD方法求出数值解,再根据数值解更新各部件的质量流和能量流,进而更新各部件的CFD边界条件。重复上述步骤即可获得船舶多结构冷凝器参数设计及CFD计算方法。
附图说明
图1为本发明的流程图;
图2多结构冷凝器管侧子单元示意图。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
针对船舶二回路系统冷凝器结构紧凑、高度集成布置的设计需求,多介质套管式集成冷凝器被广泛发展应用,现有的管壳、表面式冷凝器壳侧热力性能数值模拟方法无法满足套管式集成冷凝器壳侧多种介质真实的流动和换热过程的模拟分析。
现有的冷凝器壳侧热力性能的数值模型方法多对壳体结构采用简单的网格划分方法,对于壁面传热特性采用绝热边界或定热流边界条件进行简化,边界条件取值精度差,导致数值计算方法精度有限。
结合图1-2,针对船舶二回路系统传统结构形式的冷凝器和多介质套管式结构的集成冷凝器,本发明提出一种基于流路质量、能量流分解和二维数值模拟方法相结合的数值仿真计算方法,程序流程如附图1所示,主要包括以下步骤:
(1)将多结构冷凝器管侧按照工质流动关系,分解为一系列仅具有单进口、单出口特点的子单元(附图2),第n号子单元入口、出口变量包括循环水质量流量Gwr,n/Gwc,n和热流量Hwr,n/Hwc,n,其中,下标r表示入口,下标c表示出口,n表示单元编号。
(2)根据换热器需求,设定多结构冷凝器壳侧各工质在入口和出口处的流量Gsr,n/Gsc,n和Hsn焓值的初始值。
(3)根据管、壳间换热能量守恒原则,第n个单元的管壳式换热器满足以下关系:
Gwr,nHwr,n–Gwc,nHwc,n=(Gsr-Gsc)nHsn
可求出第n个单元对应的管侧工质出口焓值,以及管内平均速度uw和温度Tw。各子单元管侧之间按照工质流动流程具有串联、并联或独立的关系,当子单元1和子单元2具有串联关系时,两者间的变量满足:
Gwc,1=Gwr,2;Hwc,1=Hwr,2
当子单元1和子单元2具有并联关系时,各自入口、出口关系独立。但两者流量分配满足压降相等原则,即
△P1=△P2
且与各自上游、下游单元具有对应的流量关系,其Gw和Hw可依据质量守恒和能量守恒原理确定。
为冷凝器壳侧每个单元通道划分二维计算网格,根据壳侧蒸汽质量流量Gs、空气质量流量Ga、热流量Hs值、壁面边界及管侧传热特性设定各单元内部网格的us、Ts、空气质量分数Y变量初始值及边界条件。
(4)基于网格、边界条件和us/Ts/Y初始值信息,采用数值方法离散冷凝器各单元的质量、动量、空气组分及温度传输方程,并生成系数矩阵;其中,有限元离散方法采用大数法,对流项采用P1格式离散,扩散项采用中心节点法的P2格式离散,并采用SUPG方法考虑对流项占优数学问题;
(5)数值求解冷凝器各单元通道的系数矩阵,获得各单元us、Ts和Y值,并根据出口处的us、Ts、Y值求解Gsc,Gac和Hsc;
(6)判断冷凝器各单元出口处的Gsc,Gac和Hsc与上一步计算值(或初始值)的最大偏差是否小于收敛值;
(7)如果步骤(6)不收敛,则将第(6)步获得的冷凝器各单元的Gsc,Gac和Hsc作为已知量,重复步骤(3)、(4)、(5)、(6),直至收敛;
(8)如果步骤(6)收敛,则打印计算结果。
Claims (1)
1.船舶冷凝器壳侧热力性能数值模拟方法,其特征是:
(1)将多结构冷凝器管侧按照工质流动关系,分解为仅具有单进口、单出口特点的子单元,第n号子单元入口、出口变量包括循环水质量流量Gwr,n/Gwc,n和热流量Hwr,n/Hwc,n,其中下标r表示入口,下标c表示出口,n表示单元编号;
(2)设定多结构冷凝器壳侧各工质在入口和出口处的流量Gsr,n/Gsc,n和Hsn焓值的初始值;
(3)根据管、壳间换热能量守恒原则,第n个单元的管壳式换热器满足以下关系:
Gwr,nHwr,n–Gwc,nHwc,n=(Gsr-Gsc)nHsn
求出第n个单元对应的管侧工质出口焓值,以及管内平均速度uw和温度Tw,各子单元管侧之间按照工质流动流程具有串联、并联或独立的关系,当子单元1和子单元2具有串联关系时,两者间的变量满足:
Gwc,1=Gwr,2;Hwc,1=Hwr,2
当子单元1和子单元2具有并联关系时,各自入口、出口关系独立,但两者流量分配满足压降相等原则,即
△P1=△P2
且与各自上游、下游单元具有对应的流量关系,其Gw和Hw可依据质量守恒和能量守恒原理确定;
为冷凝器壳侧每个单元通道划分二维计算网格,根据壳侧蒸汽质量流量Gs、空气质量流量Ga、热流量Hs值、壁面边界及管侧传热特性设定各单元内部网格的us、Ts、空气质量分数Y变量初始值及边界条件;
(4)基于网格、边界条件和us/Ts/Y初始值信息,采用数值方法离散冷凝器各单元的质量、动量、空气组分及温度传输方程,并生成系数矩阵;其中,有限元离散方法采用大数法,对流项采用P1格式离散,扩散项采用中心节点法的P2格式离散,并采用SUPG方法考虑对流项占优数学问题;
(5)数值求解冷凝器各单元通道的系数矩阵,获得各单元us、Ts和Y值,并根据出口处的us、Ts、Y值求解Gsc,Gac和Hsc;
(6)判断冷凝器各单元出口处的Gsc,Gac和Hsc与上一步计算值或初始值的最大偏差是否小于收敛值;
(7)如果步骤(6)不收敛,则将第(6)步获得的冷凝器各单元的Gsc,Gac和Hsc作为已知量,重复步骤(3)-(6),直至收敛;
(8)如果步骤(6)收敛,则输出结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111467409.0A CN114239429B (zh) | 2021-12-02 | 2021-12-02 | 船舶冷凝器壳侧热力性能数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111467409.0A CN114239429B (zh) | 2021-12-02 | 2021-12-02 | 船舶冷凝器壳侧热力性能数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114239429A CN114239429A (zh) | 2022-03-25 |
CN114239429B true CN114239429B (zh) | 2023-02-17 |
Family
ID=80753093
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111467409.0A Active CN114239429B (zh) | 2021-12-02 | 2021-12-02 | 船舶冷凝器壳侧热力性能数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114239429B (zh) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109657372B (zh) * | 2018-12-24 | 2021-07-13 | 成都安世亚太科技有限公司 | 一种新型的管壳式换热器多尺度耦合协同换热模拟方法 |
CN110610037B (zh) * | 2019-08-30 | 2023-04-18 | 中国科学院工程热物理研究所 | 逆流式回热器热力仿真方法 |
CN113095004B (zh) * | 2021-03-29 | 2022-08-05 | 西安交通大学 | 一种液态金属流动传热计算方法 |
-
2021
- 2021-12-02 CN CN202111467409.0A patent/CN114239429B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114239429A (zh) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Polley et al. | INTERFACING HEAT EXCHANGER NETWORK VNTHESIS AND DETAILED HEAT EXCHANGER DESIGN | |
CN111144054B (zh) | 一种氟盐冷却高温堆非能动余热排出系统自然循环特性模化方法 | |
CN114266171B (zh) | 一种u型管蒸汽发生器全耦合共轭传热计算方法 | |
Adelaja et al. | Computer aided analysis of thermal and mechanical design of shell and tube heat exchangers | |
CN105512402B (zh) | 空调热换器的仿真方法 | |
CN110610037B (zh) | 逆流式回热器热力仿真方法 | |
CN114491817B (zh) | 一种船舶核动力管网仿真中的母管模拟方法及系统 | |
CN115659908B (zh) | 一种印刷电路板换热器的三单元非平衡多孔介质方法 | |
CN114048572B (zh) | 一种大型变物性管壳式换热器的设计计算方法 | |
CN114239429B (zh) | 船舶冷凝器壳侧热力性能数值模拟方法 | |
CN115017674A (zh) | 一种基于Modelica的管壳式换热器建模仿真方法 | |
CN109959283B (zh) | 高温冷却器热力校核方法及系统 | |
Torresi et al. | CFD analysis of the flow through tube banks of HRSG | |
CN116150993A (zh) | 一种基于频域特性的换热器设计方法及换热器 | |
Teck et al. | Simulation of Shell and Tube Heat Exchanger: Influence of the Lower Flows and the Baffles on a Fluid Dynamics | |
Sodja et al. | Some aspects of the modeling of tube-and-shell heat-exchangers | |
Han et al. | Simplification method of thermal-fluid network with circulation reflux based on matrix operation | |
Ramakanth et al. | Prediction of channel diameter to reduce flow mal distribution in radiators using ANN | |
Sarmiento et al. | A numerical method for shell and thermosyphon heat exchanger analysis | |
Węglarz et al. | New calculation method for tube cross-flow heat exchangers | |
Cabello et al. | Validation of CFD models for double-pipe heat exchangers with empirical correlations | |
CN117350188A (zh) | 一种基于模型降阶的核反应堆热工回路系统快速计算方法 | |
Hur et al. | Numerical analysis of heat transfer in a plate heat exchanger | |
Libreros et al. | Critical Review of the Theoretical, Experimental and Computational Fluid Dynamics Methods for Designing Plate Fin Heat Exchangers. | |
CN113619810B (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 |