CN107194076A - 一种管内超临界流体气化‑管外膜状冷凝传热计算方法 - Google Patents
一种管内超临界流体气化‑管外膜状冷凝传热计算方法 Download PDFInfo
- Publication number
- CN107194076A CN107194076A CN201710376075.3A CN201710376075A CN107194076A CN 107194076 A CN107194076 A CN 107194076A CN 201710376075 A CN201710376075 A CN 201710376075A CN 107194076 A CN107194076 A CN 107194076A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- heat transfer
- infinitesimal
- 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
Links
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
-
- 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)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
Abstract
本发明公开了一种管内超临界流体气化‑管外膜状冷凝传热计算方法,包括以下步骤:输入传热工艺与换热管参数;根据换热管长度划分换热微元,假设外壁过冷度初值;计算微元j外壁温度、管外冷凝传热系数及热负荷初值;计算微元j内壁温度、管内超临界传热系数;核算微元j热负荷;更新微元j外壁温度、管外传热系数;判断微元j传热计算是否收敛;计算微元j管内压降;计算微元j管内流体出口压力与温度;判断是否计算至最后一个微元;进入下一微元;完成计算。本发明可用于低温流体在超临界压力下气化时卧式换热器的传热校核计算,适用于流体进出口温差大、采用光滑换热管或不同类型强化传热管工况,计算精度高。
Description
技术领域
本发明属于低温流体超临界传热技术领域,具体涉及一种管内超临界流体气化-管外膜状冷凝传热计算方法。
背景技术
将工业气体液化以液态形式储运、并在用户终端使其再气化的技术被广泛应用于液化天然气(LNG)、液氧、液氮、液氢、液态二氧化碳等领域,可大大节省储运空间、提高储运效率。低温流体气化器是实现该过程的关键装备之一。
由于全球范围内的能源危机和环境恶化,LNG工业快速发展。为满足LNG接收站气化低温LNG的需要,专利EP0550845B1提出一种开架式气化器(ORV)结构,通过将海水自上而下喷淋在板型换热管束外表面上,把能量传递给LNG使其气化;专利US3818893提出一种浸没燃烧式气化器(SCV)型式,将天然气燃烧后产生的高温烟气排入水浴池中加热并剧烈搅拌海水,浸泡在海水中的换热盘管内LNG被池水加热气化;US6164247提出一种中间流体式气化器(IFV)结构,将蒸发器、凝结器与调温器三组管壳式换热器组合起来,在蒸发器中用海水加热壳程的低沸点中间介质使其气化,在凝结器中间介质管外冷凝提供热量气化管内的LNG,在调温器中LNG被进一步加热以达到要求的温度。目前,ORV、SCV与IFV已成为世界范围内广泛应用的三种基本负荷型LNG气化器。
IFV是一种具有广泛应用潜力的低温流体气化设备。IFV相比于ORV,它避免了用海水直接气化LNG时可能带来的结冰问题,且适用于海水含沙量较高的水域,因而适应性更强。IFV相比于SCV,它不需要额外燃烧天然气,因而能效更好,且没有CO2、NOx等环境污染问题。由于IFV的组合式气化结构设计,它可用于海水水质较差的陆上LNG气化终端、空间和重量受限的海上浮式储存与气化平台、以及冷能发电系统等不同场合。
为进一步提升IFV在工程应用中的技术优势,专利US6164247提出在蒸发器/冷凝器共用的壳体内设置分隔板和支撑板,以适应海上平台晃荡工况,保持稳定的气化效率;专利US6367429B2提出当采用工业回收热水作为热源时,可将蒸发器设计成两管程或多管程,以使结构更紧凑;专利US6367265B1提出在凝结器管箱中设置LNG喷雾冷却装置,以降低IFV启动或低负荷运行时的温差应力;专利US20160146403A1提出凝结器采用管外带翅片或沟槽的强化传热管,以强化中间介质的管外冷凝传热过程。
由于IFV采用组合式气化结构、设备紧凑、LNG温度及物性变化范围大,无法采用现有商业软件直接设计。文献中主要采用分布式参数模型(DPM)进行设计型和校核型计算。该方法是将IFV各传热段沿管长方向划分成多个微元,对每个微元进行能量衡算和总传热系数核算,其假设整个换热过程是稳态的、不存在向环境的漏热。文献中一般将各传热段均按单管程、单壳程处理,并忽略压降。例如,文献[大型LNG中间介质气化器换热面积计算方法,天然气与石油,31,31-35,2013]提出了超临界工况下IFV各传热段换热面积计算方法;文献[Thermal performance analysis of intermediate fluid vaporizer for liquefiednatural gas,Applied Thermal Engineering,65,564-574,2014]对超临界工况IFV进行了校核型计算,分析了海水进口温度与流量、LNG进口压力与流量等工艺参数对传热特性的影响。
但上述现有技术忽略了管内超临界流体流经换热管的压力降、不适用于采用强化传热管时IFV的传热计算、且采用的管内超临界流体传热系数关联式预测精度不高,因而传热计算精度低、适用范围受到限制。
发明内容
为了解决上述技术问题,本发明提供一种管内超临界流体气化-管外膜状冷凝传热计算方法。
为了实现本发明的目的,本发明采用了以下技术方案:
一种管内超临界流体气化-管外膜状冷凝传热计算方法,包括以下步骤:
S1、输入传热工艺参数与水平换热管参数,包括:
管内流体的组分及其含量、入口温度Tb (1)、入口压力Pb (1)、质量通量Gb、拟临界温度Tpc,其中管内流体的入口压力大于管内流体的临界压力;
管外流体的组分及其含量、饱和温度Tsat;
换热管的内直径di、外直径do、长度L、导热系数kw、单列换热管束中沿竖直方向的换热管排数Neff;
S2、根据换热管长度划分N个相等长度的传热微元,每个传热微元长度为dL=L/N;所述传热微元数目N满足:管内流体进出换热管温度差/N≤20K;
假设所述换热管外壁面过冷度ΔTsub的初值记为ΔTsub,0,所述ΔTsub,0取0.5~20K;
S3、计算传热微元j(j=1,2,…,N)外壁面温度、管外冷凝传热系数及热负荷初值
根据换热管外壁面过冷度初值ΔTsub,0,按公式(1)~(3)计算传热微元j外壁面温度Two、管外冷凝传热系数ho及热负荷dQ的初值,并分别记为Two,0 (j)、ho,0 (j)与dQ0 (j):
Two=Tsat-ΔTsub (1)
dQ=hoΔTsubπdodL (3)
式中:g为重力加速度;
ρL、ρV、kL、r、μL分别为管外流体在饱和温度Tsat下的液相密度、汽相密度、液相导热系数、蒸发焓、液相粘度;
m为管束效应修正系数,按如下方式进行取值:
在换热器样机中选取具有X排换热管的单列换热管束,按照GB/T27698.2-2011汽-液冷凝器用热交换器的测试程序,测试每根换热管内流体的进出口温度、压力和流量,管外流体的进出口温度、压力和流量,并计算沿竖直方向自上而下布置的第1排换热管外的冷凝传热系数ho (1)与第M排换热管外的冷凝传热系数ho (M),其中X>M,按公式(4)回归得到m值:
ho (M)=ho (1)·[M1-m-(M-1)1-m] (4)
S4、计算传热微元j内壁面温度、管内超临界流体传热系数
根据传热微元j外壁面温度初值Two,0 (j)和热负荷初值dQ0 (j),按公式(5)计算内壁面温度Twi (j);根据传热微元j的入口温度Tb (j)、入口压力Pb (j),按公式(6)计算管内超临界流体传热系数hi (j):
式中:
Tb以管内流体在传热微元j入口处的温度Tb (j)代入;
ρb、Cpb、kb、μb、Hb分别以管内流体在传热微元j入口温度Tb (j)、入口压力Pb (j)下的密度、定压比热、导热系数、粘度和焓值代入;
ρwi、Hwi分别以管内流体在传热微元j内壁面温度Twi (j)、入口压力Pb (j)下的密度和焓值代入;
S5、核算传热微元j热负荷
根据传热微元j内壁面温度Twi (j)、管内流体传热系数hi (j),按公式(7)核算传热微元j的热负荷dQ(j):
dQ=hi(Twi-Tb)πdidL (7)
式中,Tb以管内流体在传热微元j入口处的温度Tb (j)代入;
S6、更新传热微元j外壁面温度、管外冷凝传热系数
以核算后的传热微元j热负荷dQ(j)代入公式(1)~(3),求解更新后的传热微元j外壁面温度Two (j)、管外冷凝传热系数ho (j);
S7、判断传热微元j传热计算是否收敛
根据传热微元j热负荷的初值dQ0 (j)与核算值dQ(j),按公式(8)计算相对误差ε:
若ε≤0.001,则传热微元j传热计算收敛,进入S8;
否则,返回S4,依据更新后的外壁面温度和热负荷,重新计算内壁面温度、管内超临界流体传热系数;
S8、计算传热微元j管内流体压力降
根据传热微元j的入口温度Tb (j)、入口压力Pb (j),按公式(9)~(10)计算传热微元j管内流体压力降ΔPb,并记为ΔPb (j):
式中:f为直管摩擦阻力系数;
R为U形换热管的弯曲半径;当采用直换热管时,取R=0;
S9、计算传热微元j管内流体出口压力与温度
根据传热微元j入口处管内流体的压力Pb (j)和流经传热微元j的压力降ΔPb (j),按公式(11)计算传热微元j出口处管内流体的压力Pb (j+1);根据传热微元j入口处管内流体的焓值Hb (j),按公式(12)计算传热微元j出口处管内流体的焓值Hb (j+1):
根据传热微元j出口处管内流体的焓值Hb (j+1)和压力Pb (j+1),通过查询物质理化性能手册或采用热力学方程计算方法,得到传热微元j出口处管内流体的温度Tb (j+1);
S10、判断是否计算至最后一个传热微元,即当前j值是否等于N
若是,则进入步骤S12;若否,则进入步骤S11;
S11、进行下一传热微元的传热计算
进入下一传热微元j+1,并执行步骤S3;
S12、完成传热计算。
进一步,所述步骤S3中m值根据换热管的类型按照如下方式进行取值:
对于普通光滑管,m=0.1667~0.25;
对于二维外翅片管,m=0.01~0.1;
对于三维非连续翅片管,m=0.23~0.36。
进一步,所述换热器由多列、多排换热管组成时,所述步骤S3中的管外冷凝传热系数ho按公式(13)计算平均值:
式中,y表示第y列换热管束;Y表示换热器共有Y列换热管束;h0,y表示第y列换热管束的管外冷凝传热系数;Neff,y表示第y列换热管束中沿竖直方向的换热管排数。
本发明的有益效果在于:
本发明提出了一种管内超临界流体气化-管外膜状冷凝传热计算方法,尤其是一种卧式换热器中管内超临界流体气化-管外膜状冷凝传热过程的校核型计算方法,考虑了管内超临界流体的压力降、管外膜状冷凝传热的管束效应影响,采用了高精度的传热系数关联式,适用于流体进出口温差大、采用普通光滑换热管或不同类型强化传热管工况,计算精度高。
附图说明
图1为本发明中管内超临界流体气化-管外膜状冷凝传热计算流程示意图。
图2为本发明实施例1计算的某LNG气化过程中沿管长方向的温度分布与管内超临界流体传热系数分布结果。
图3为本发明实施例2计算的某LNG气化过程中沿管长方向的温度分布与管内超临界流体传热系数分布结果。
图4为本发明实施例3计算的某LNG气化过程中沿管长方向的温度分布与管内超临界流体传热系数分布结果。
具体实施方式
下面结合结合附图和实施例对本发明技术方案做出更为具体的说明:
如附图1所示,一种管内超临界流体气化-管外膜状冷凝传热计算方法,适用于低温流体在超临界压力下气化时卧式换热器的传热校核计算,包括以下步骤:
S1、输入传热工艺参数与水平换热管参数,包括:
管内流体的组分及其含量、入口温度Tb (1)、入口压力Pb (1)、质量通量Gb、拟临界温度Tpc,其中管内流体的入口压力大于管内流体的临界压力;
管外流体的组分及其含量、饱和温度Tsat;
换热管的内直径di、外直径do、长度L、导热系数kw、单列换热管束中沿竖直方向的换热管排数Neff;
所述管内流体的组分及其含量、管外流体的组分及其含量主要用于计算相应流体的物理性质参数,比如密度、粘度、导热系数、焓等。
所述含量可以是质量含量或体积含量或摩尔含量。
S2、根据换热管长度划分N个相等长度的传热微元,每个传热微元长度为dL=L/N;所述传热微元数目N满足:管内流体进出换热管温度差/N≤20K;
假设所述换热管外壁面过冷度ΔTsub的初值记为ΔTsub,0,所述ΔTsub,0取0.5~20K;
S3、计算传热微元j(j=1,2,…,N)外壁面温度、管外冷凝传热系数及热负荷初值
根据换热管外壁面过冷度初值ΔTsub,0,按公式(1)~(3)计算传热微元j外壁面温度Two、管外冷凝传热系数ho及热负荷dQ的初值,并分别记为Two,0 (j)、ho,0 (j)与dQ0 (j):
Two=Tsat-ΔTsub (1)
dQ=hoΔTsubπdodL (3)
式中:g为重力加速度;
ρL、ρV、kL、r、μL分别为管外流体在饱和温度Tsat下的液相密度、汽相密度、液相导热系数、蒸发焓、液相粘度;
m为管束效应修正系数,按如下方式进行取值:
在换热器样机中选取具有X排换热管的单列换热管束,按照GB/T27698.2-2011汽-液冷凝器用热交换器的测试程序,测试每根换热管内流体的进出口温度、压力和流量,管外流体的进出口温度、压力和流量,并计算沿竖直方向自上而下布置的第1排换热管外的冷凝传热系数ho (1)与第M排换热管外的冷凝传热系数ho (M),按公式(4)回归得到m值:
ho (M)=ho (1)·[M1-m-(M-1)1-m] (4)
换热器管束外冷凝传热过程存在边界效应,即末排换热管外冷凝传热系数不降反升。为避免边界效应对冷凝传热试验过程的影响,所述换热器样机中单列换热管束的换热管总排数X需大于待测换热管所在排数M。
S4、计算传热微元j内壁面温度、管内超临界流体传热系数
根据传热微元j外壁面温度初值Two,0 (j)和热负荷初值dQ0 (j),按公式(5)计算内壁面温度Twi (j);根据传热微元j的入口温度Tb (j)、入口压力Pb (j),按公式(6)计算管内超临界流体传热系数hi (j):
式中:
Tb以管内流体在传热微元j入口处的温度Tb (j)代入;
ρb、Cpb、kb、μb、Hb分别以管内流体在传热微元j入口温度Tb (j)、入口压力Pb (j)下的密度、定压比热、导热系数、粘度和焓值代入;
ρwi、Hwi分别以管内流体在传热微元j内壁面温度Twi (j)、入口压力Pb (j)下的密度和焓值代入;
所述传热微元j的入口温度Tb (j)、入口压力Pb (j),是由步骤S9计算得到的传热微元(j-1)的出口温度和出口压力;对于传热微元1即是指整根换热管的入口温度Tb (1)、入口压力Pb (1);
S5、核算传热微元j热负荷
根据传热微元j内壁面温度Twi (j)、管内流体传热系数hi (j),按公式(7)核算传热微元j的热负荷dQ(j):
dQ=hi(Twi-Tb)πdidL (7)
式中,Tb以管内流体在传热微元j入口处的温度Tb (j)代入;
S6、更新传热微元j外壁面温度、管外冷凝传热系数
以核算后的传热微元j热负荷dQ(j)代入公式(1)~(3),求解更新后的传热微元j外壁面温度Two (j)、管外冷凝传热系数ho (j);
S7、判断传热微元j传热计算是否收敛
根据传热微元j热负荷的初值dQ0 (j)与核算值dQ(j),按公式(8)计算相对误差ε:
若ε≤0.001,则传热微元j传热计算收敛,进入S8;
否则,返回S4,依据更新后的外壁面温度和热负荷,重新计算内壁面温度、管内超临界流体传热系数;
S8、计算传热微元j管内流体压力降
根据传热微元j的入口温度Tb (j)、入口压力Pb (j),按公式(9)~(10)计算传热微元j管内流体压力降ΔPb,并记为ΔPb (j):
式中:f为直管摩擦阻力系数;
R为U形换热管的弯曲半径;当采用直换热管时,取R=0;
S9、计算传热微元j管内流体出口压力与温度
根据传热微元j入口处管内流体的压力Pb (j)和流经传热微元j的压力降ΔPb (j),按公式(11)计算传热微元j出口处管内流体的压力Pb (j+1);根据传热微元j入口处管内流体的焓值Hb (j),按公式(12)计算传热微元j出口处管内流体的焓值Hb (j+1):
根据传热微元j出口处管内流体的焓值Hb (j+1)和压力Pb (j+1),通过查询物质理化性能手册或采用热力学方程计算方法,得到传热微元j出口处管内流体的温度Tb (j+1);
S10、判断是否计算至最后一个传热微元,即当前j值是否等于N
若是,则进入步骤S12;若否,则进入步骤S11;
S11、进行下一传热微元的传热计算
进入下一传热微元j+1,并执行步骤S3;
S12、完成传热计算。
所述步骤S3中m值根据换热管的类型按照如下方式进行取值:
对于普通光滑管,m=0.1667~0.25;
对于二维外翅片管,m=0.01~0.1;
对于三维非连续翅片管,m=0.23~0.36。
换热管通常分为普通光滑管和强化传热管两种类型,所述强化传热管按表面结构的不同,又分为二维外翅片管、异形强化管、三维非连续翅片管。其中所述普通光滑管为不带翅片且壁面光滑的换热管;所述二维外翅片管是将铝带、铜带或钢带构成的翅片缠绕在换热管表面,所述二维外翅片管的翅片表面规则;所述三维非连续翅片管的翅片比二维外翅片管的翅片具有更不规则的扩展表面。
所述换热器由多列、多排换热管组成时,所述步骤S3中的管外冷凝传热系数ho按公式(13)计算平均值:
式中,y表示第y列换热管束;Y表示换热器共有Y列换热管束;h0,y表示第y列换热管束的管外冷凝传热系数;Neff,y表示第y列换热管束中沿竖直方向的换热管排数。
实施例1
以某LNG气化过程为例,依据本发明计算方法,完成该气化过程中换热管的校核型传热计算,具体步骤如下:
S1、输入传热工艺参数与换热管参数
LNG气化工况1中,LNG组成见表1,传热工艺参数见表2,换热管参数见表3。
表1 LNG气化工况1的LNG组分及其含量
组分 | 甲烷 | 乙烷 | 丙烷 | 正丁烷 | 氮 |
摩尔含量/% | 84.5 | 12.9 | 1.5 | 0.5 | 0.6 |
表2 LNG气化工况1的传热工艺参数
表3 LNG气化工况1的换热管参数
类型 | 普通光滑管、直管 |
材质 | 不锈钢 |
内直径di/m | 0.015 |
外直径do/m | 0.019 |
长度L/m | 9 |
导热系数kw/W·m-1·K-1 | 16.2 |
管排深度Neff | 1 |
U形管弯曲半径R/m | 0 |
S2、根据换热管长度划分N个相等长度的传热微元
传热微元数目N=200;每个传热微元长度为dL=9/200=0.045m;
假设换热管外壁面过冷度ΔTsub的初值,记为ΔTsub,0,取10K;
S3、计算传热微元j(j=1,2,…,N)外壁面温度、管外冷凝传热系数及热负荷初值
以j=1为例:
按公式(1)计算,Two,0 (1)=260.15K;
对于光滑换热管,取m=0.1667,按公式(2)计算,ho,0 (1)=2126.95W·m-2·K-1;
按公式(3)计算,dQ0 (1)=57.13W;
S4、计算传热微元j内壁面温度、管内超临界流体传热系数
以j=1为例:
按公式(5)计算,Twi (1)=257.20K;
传热微元1入口温度为Tb (1)=111.15,入口压力Pb (1)=12000kPa,按公式(6)计算,此时Tb (1)<Tpc<Twi (1),因此n=0.4+0.2(Twi/Tpc-1)=0.4093,hi (1)=957.14W·m-2·K-1;
S5、核算传热微元j热负荷
以j=1为例:
按公式(7)核算,dQ(1)=296.44W;
S6、更新传热微元j外壁面温度、管外冷凝传热系数
以j=1为例:
将dQ(1)=296.44W带入公式(1)~(3),求解得更新后的Two (1)=180.32K,ho (1)=1228.58W·m-2·K-1;
S7、判断传热微元j传热计算是否收敛
以j=1为例:
按公式(8)计算,相对误差ε=0.807>0.001,传热计算不收敛,返回S4,按Two,0 (1)=180.32K、dQ0 (1)=296.44W,重新计算内壁面温度、管内超临界流体传热系数;
重复上述迭代过程,最终传热计算收敛时结果为:
Two (1)=216.14K,ho (1)=1477.59W·m-2·K-1,dQ(1)=214.34W,Twi (1)=205.04K,hi (1)=1076.60W·m-2·K-1;
S8、计算传热微元j管内流体压力降
以j=1为例:
按公式(9)~(10)计算,ΔPb (1)=0.0034kPa;
S9、计算传热微元j管内流体出口压力与温度
以j=1为例:
按公式(11)计算,Pb (2)=11999.997kPa;
按公式(12)计算,Hb (2)=-18310.97J/kg;
根据传热微元1出口处管内流体的焓值Hb (2)和压力Pb (2),通过调用美国国家标准与技术研究院开发的工质物性计算软件(NIST Refprop),得到传热微元1出口处管内流体的温度Tb (2)=113.15K;
S10、判断是否计算至最后一个传热微元
当前j值为1,不等于N,即未计算至最后一个传热微元,进入S11;
S11、进入下一传热微元
令j=2,进入下一传热微元,并执行S3;
如此依序计算每个传热微元,直到完成所有传热微元的传热计算。
S12、完成传热计算。
本实施例中,利用本发明方法计算的管内LNG沿管长方向的温度分布与管内超临界流体传热系数分布结果如图2所示。LNG最终的出口温度为255.12K,LNG进、出换热管总的压力降为0.79kPa,管内超临界流体传热系数为1076.60~1558.19W·m-2·K-1。
实施例2
LNG气化工况2,LNG组成与传热工艺参数同表1与表2,换热管参数见表4。
表4 LNG气化工况2的换热管参数
类型 | 普通光滑管、U形管 |
材质 | 不锈钢 |
内直径di/m | 0.015 |
外直径do/m | 0.019 |
长度L/m | 9 |
导热系数kw/W·m-1·K-1 | 16.2 |
管排深度Neff | 10 |
U形管弯曲半径R/m | 0.04 |
按本发明的传热计算方法,划分传热微元数目N=200,取外壁面过冷度初值ΔTsub,0=0.5K,取管束效应修正系数m=0.1667。
本实施例中,利用本发明方法计算的管内LNG沿管长方向的温度分布与管内超临界流体传热系数分布结果如图3所示。LNG最终的出口温度为248.84K,LNG进、出换热管总的压力降为7.08kPa,管内超临界流体传热系数为1094.65~1577.07W·m-2·K-1。
实施例3
LNG气化工况3,LNG组成如表5所示,传热工艺参数见表6,换热管参数见表7。
表5 LNG气化工况3的LNG组分及其含量
组分 | 甲烷 | 乙烷 | 丙烷 | 正丁烷 | 氮 |
摩尔含量/% | 95 | 4.6 | 0.38 | 0 | 0.02 |
表6 LNG气化工况3的传热工艺参数
表7 LNG气化工况3的换热管参数
类型 | 二维外翅片管、直管 |
材质 | 不锈钢 |
内直径di/m | 0.015 |
外直径do/m | 0.019 |
长度L/m | 9 |
导热系数kw/W·m-1·K-1 | 16.2 |
管排深度Neff | 1 |
U形管弯曲半径R/m | 0 |
按本发明的传热计算方法,划分传热微元数目N=200,取外壁面过冷度初值ΔTsub,0=1K,取管束效应修正系数m=0.04。
本实施例中,利用本发明方法计算的管内LNG沿管长方向的温度分布与管内超临界流体传热系数分布结果如图4所示。LNG最终的出口温度为252.66K,LNG进、出换热管总的压力降为0.89kPa,管内超临界流体传热系数为1176.25~1567.55W·m-2·K-1。
Claims (3)
1.一种管内超临界流体气化-管外膜状冷凝传热计算方法,其特征在于包括以下步骤:
S1、输入传热工艺参数与水平换热管参数,包括:
管内流体的组分及其含量、入口温度Tb (1)、入口压力Pb (1)、质量通量Gb、拟临界温度Tpc,其中管内流体的入口压力大于管内流体的临界压力;
管外流体的组分及其含量、饱和温度Tsat;
换热管的内直径di、外直径do、长度L、导热系数kw、单列换热管束中沿竖直方向的换热管排数Neff;
S2、根据换热管长度划分N个相等长度的传热微元,每个传热微元长度为dL=L/N;所述传热微元数目N满足:管内流体进出换热管温度差/N≤20K;
假设所述换热管外壁面过冷度ΔTsub的初值记为ΔTsub,0,所述ΔTsub,0取0.5~20K;
S3、计算传热微元j(j=1,2,…,N)外壁面温度、管外冷凝传热系数及热负荷初值
根据换热管外壁面过冷度初值ΔTsub,0,按公式(1)~(3)计算传热微元j外壁面温度Two、管外冷凝传热系数ho及热负荷dQ的初值,并分别记为Two,0 (j)、ho,0 (j)与dQ0 (j):
Two=Tsat-ΔTsub (1)
<mrow>
<msub>
<mi>h</mi>
<mi>o</mi>
</msub>
<mo>=</mo>
<mn>0.79</mn>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<msub>
<mi>g&rho;</mi>
<mi>L</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&rho;</mi>
<mi>V</mi>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>k</mi>
<mi>L</mi>
<mn>3</mn>
</msubsup>
<mi>r</mi>
</mrow>
<mrow>
<msub>
<mi>&mu;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>d</mi>
<mi>o</mi>
</msub>
<msub>
<mi>&Delta;T</mi>
<mrow>
<mi>s</mi>
<mi>u</mi>
<mi>b</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mn>0.25</mn>
</msup>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>e</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mi>m</mi>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
dQ=hoΔTsubπdodL (3)
式中:g为重力加速度;
ρL、ρV、kL、r、μL分别为管外流体在饱和温度Tsat下的液相密度、汽相密度、液相导热系数、蒸发焓、液相粘度;
m为管束效应修正系数,按如下方式进行取值:
在换热器样机中选取具有X排换热管的单列换热管束,按照GB/T27698.2-2011汽-液冷凝器用热交换器的测试程序,测试每根换热管内流体的进出口温度、压力和流量,管外流体的进出口温度、压力和流量,并计算沿竖直方向自上而下布置的第1排换热管外的冷凝传热系数ho (1)与第M排换热管外的冷凝传热系数ho (M),其中X>M,按公式(4)回归得到m值:
ho (M)=ho (1)·[M1-m-(M-1)1-m] (4)
S4、计算传热微元j内壁面温度、管内超临界流体传热系数
根据传热微元j外壁面温度初值Two,0 (j)和热负荷初值dQ0 (j),按公式(5)计算内壁面温度Twi (j);根据传热微元j的入口温度Tb (j)、入口压力Pb (j),按公式(6)计算管内超临界流体传热系数hi (j):
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>Q</mi>
</mrow>
<msub>
<mi>k</mi>
<mi>w</mi>
</msub>
</mfrac>
<mfrac>
<mrow>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>o</mi>
</msub>
<mo>/</mo>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>d</mi>
<mi>L</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mn>0.0183</mn>
<msubsup>
<mi>Re</mi>
<mi>b</mi>
<mn>0.82</mn>
</msubsup>
<msubsup>
<mi>Pr</mi>
<mi>b</mi>
<mn>0.5</mn>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>w</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mi>b</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>0.3</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<mi>p</mi>
<mi>b</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<mfrac>
<msub>
<mi>k</mi>
<mi>b</mi>
</msub>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
1
式中:
Tb以管内流体在传热微元j入口处的温度Tb (j)代入;
ρb、Cpb、kb、μb、Hb分别以管内流体在传热微元j入口温度Tb (j)、入口压力Pb (j)下的密度、定压比热、导热系数、粘度和焓值代入;
ρwi、Hwi分别以管内流体在传热微元j内壁面温度Twi (j)、入口压力Pb (j)下的密度和焓值代入;
S5、核算传热微元j热负荷
根据传热微元j内壁面温度Twi (j)、管内流体传热系数hi (j),按公式(7)核算传热微元j的热负荷dQ(j):
dQ=hi(Twi-Tb)πdidL (7)
式中,Tb以管内流体在传热微元j入口处的温度Tb (j)代入;
S6、更新传热微元j外壁面温度、管外冷凝传热系数
以核算后的传热微元j热负荷dQ(j)代入公式(1)~(3),求解更新后的传热微元j外壁面温度Two (j)、管外冷凝传热系数ho (j);
S7、判断传热微元j传热计算是否收敛
根据传热微元j热负荷的初值dQ0 (j)与核算值dQ(j),按公式(8)计算相对误差ε:
<mrow>
<mi>&epsiv;</mi>
<mo>=</mo>
<mo>|</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>Q</mi>
<mo>-</mo>
<msub>
<mi>dQ</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>Q</mi>
</mrow>
</mfrac>
<mo>|</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
若ε≤0.001,则传热微元j传热计算收敛,进入S8;
否则,返回S4,依据更新后的外壁面温度和热负荷,重新计算内壁面温度、管内超临界流体传热系数;
S8、计算传热微元j管内流体压力降
根据传热微元j的入口温度Tb (j)、入口压力Pb (j),按公式(9)~(10)计算传热微元j管内流体压力降ΔPb,并记为ΔPb (j):
<mrow>
<mi>f</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>64</mn>
<mo>/</mo>
<msub>
<mi>Re</mi>
<mi>b</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Re</mi>
<mi>b</mi>
</msub>
<mo><</mo>
<mn>2300</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.316</mn>
<mo>/</mo>
<msubsup>
<mi>Re</mi>
<mi>b</mi>
<mn>0.25</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>2300</mn>
<mo>&le;</mo>
<msub>
<mi>Re</mi>
<mi>b</mi>
</msub>
<mo><</mo>
<mn>8000</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.0054</mn>
<mo>+</mo>
<mn>0.3964</mn>
<mo>/</mo>
<msubsup>
<mi>Re</mi>
<mi>b</mi>
<mn>0.3</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Re</mi>
<mi>b</mi>
</msub>
<mo>&GreaterEqual;</mo>
<mn>8000</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mi>b</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mi>f</mi>
<mfrac>
<mi>L</mi>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>+</mo>
<mn>0.294</mn>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>R</mi>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&CenterDot;</mo>
<mfrac>
<msubsup>
<mi>G</mi>
<mi>b</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mn>2</mn>
<msub>
<mi>&rho;</mi>
<mi>b</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:f为直管摩擦阻力系数;
R为U形换热管的弯曲半径;当采用直换热管时,取R=0;
S9、计算传热微元j管内流体出口压力与温度
根据传热微元j入口处管内流体的压力Pb (j)和流经传热微元j的压力降ΔPb (j),按公式(11)计算传热微元j出口处管内流体的压力Pb (j+1);根据传热微元j入口处管内流体的焓值Hb (j),按公式(12)计算传热微元j出口处管内流体的焓值Hb (j+1):
<mrow>
<msubsup>
<mi>P</mi>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&Delta;P</mi>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>H</mi>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mn>4</mn>
<mo>&CenterDot;</mo>
<mi>d</mi>
<mi>Q</mi>
</mrow>
<mrow>
<msubsup>
<mi>&pi;d</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<msub>
<mi>G</mi>
<mi>b</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msubsup>
<mi>H</mi>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
根据传热微元j出口处管内流体的焓值Hb (j+1)和压力Pb (j+1),通过查询物质理化性能手册或采用热力学方程计算方法,得到传热微元j出口处管内流体的温度Tb (j+1);
S10、判断是否计算至最后一个传热微元,即当前j值是否等于N
若是,则进入步骤S12;若否,则进入步骤S11;
S11、进行下一传热微元的传热计算
进入下一传热微元j+1,并执行步骤S3;
S12、完成传热计算。
2.如权利要求1所述的管内超临界流体气化-管外膜状冷凝传热计算方法,其特征在于:所述步骤S3中m值根据换热管的类型按照如下方式进行取值:
对于普通光滑管,m=0.1667~0.25;
对于二维外翅片管,m=0.01~0.1;
对于三维非连续翅片管,m=0.23~0.36。
3.如权利要求1所述的管内超临界流体气化-管外膜状冷凝传热计算方法,其特征在于:所述换热器由多列、多排换热管组成时,所述步骤S3中的管外冷凝传热系数ho按公式(13)计算平均值:
<mrow>
<msub>
<mi>h</mi>
<mi>o</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Y</mi>
</munderover>
<msub>
<mi>N</mi>
<mrow>
<mi>e</mi>
<mi>f</mi>
<mi>f</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
</msub>
<msub>
<mi>h</mi>
<mrow>
<mi>o</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Y</mi>
</munderover>
<msub>
<mi>N</mi>
<mrow>
<mi>e</mi>
<mi>f</mi>
<mi>f</mi>
<mo>,</mo>
<mi>y</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,y表示第y列换热管束;Y表示换热器共有Y列换热管束;h0,y表示第y列换热管束的管外冷凝传热系数;Neff,y表示第y列换热管束中沿竖直方向的换热管排数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710376075.3A CN107194076B (zh) | 2017-05-25 | 2017-05-25 | 一种管内超临界流体气化-管外膜状冷凝传热计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710376075.3A CN107194076B (zh) | 2017-05-25 | 2017-05-25 | 一种管内超临界流体气化-管外膜状冷凝传热计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107194076A true CN107194076A (zh) | 2017-09-22 |
CN107194076B CN107194076B (zh) | 2019-12-10 |
Family
ID=59875740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710376075.3A Active CN107194076B (zh) | 2017-05-25 | 2017-05-25 | 一种管内超临界流体气化-管外膜状冷凝传热计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107194076B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109298018A (zh) * | 2018-11-16 | 2019-02-01 | 南京工业大学 | 一种可应用多种冷却介质和模拟不同重力加速度环境的喷雾冷却试验台 |
CN111581584A (zh) * | 2020-05-21 | 2020-08-25 | 西安石油大学 | 一种地热开发过程中的压降换热定量计算方法 |
CN111767636A (zh) * | 2020-05-21 | 2020-10-13 | 浙江工业大学 | 单通错流换热器的显式传热计算方法 |
CN111814351A (zh) * | 2020-07-23 | 2020-10-23 | 清华大学 | 一种超临界压力流体换热器制备方法及系统 |
CN113361158A (zh) * | 2021-05-31 | 2021-09-07 | 江西理工大学 | 一种连续铸造结晶器内腔金属温度的测量方法及系统 |
CN113378105A (zh) * | 2021-06-24 | 2021-09-10 | 东北电力大学 | 一种超临界逆流式省煤器过程传热系数的计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0550845B1 (en) * | 1991-12-12 | 1996-09-18 | KABUSHIKI KAISHA KOBE SEIKO SHO also known as Kobe Steel Ltd. | Vaporizer for liquefied natural gas |
CN104406426A (zh) * | 2014-10-14 | 2015-03-11 | 西安倍更能源技术有限公司 | 一种带有中间传热介质的气化器设计方法 |
-
2017
- 2017-05-25 CN CN201710376075.3A patent/CN107194076B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0550845B1 (en) * | 1991-12-12 | 1996-09-18 | KABUSHIKI KAISHA KOBE SEIKO SHO also known as Kobe Steel Ltd. | Vaporizer for liquefied natural gas |
CN104406426A (zh) * | 2014-10-14 | 2015-03-11 | 西安倍更能源技术有限公司 | 一种带有中间传热介质的气化器设计方法 |
Non-Patent Citations (2)
Title |
---|
白宇恒 等: "大型LNG中间介质气化器换热面积计算方法", 《天然气与石油》 * |
陈永东 等: "LNG成套装置换热器关键技术分析", 《天然气工业》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109298018A (zh) * | 2018-11-16 | 2019-02-01 | 南京工业大学 | 一种可应用多种冷却介质和模拟不同重力加速度环境的喷雾冷却试验台 |
CN109298018B (zh) * | 2018-11-16 | 2024-04-02 | 南京工业大学 | 一种可应用多种冷却介质和模拟不同重力加速度环境的喷雾冷却试验台 |
CN111581584A (zh) * | 2020-05-21 | 2020-08-25 | 西安石油大学 | 一种地热开发过程中的压降换热定量计算方法 |
CN111767636A (zh) * | 2020-05-21 | 2020-10-13 | 浙江工业大学 | 单通错流换热器的显式传热计算方法 |
CN111581584B (zh) * | 2020-05-21 | 2023-03-24 | 西安石油大学 | 一种地热开发过程中的压降换热定量计算方法 |
CN111767636B (zh) * | 2020-05-21 | 2024-04-12 | 浙江工业大学 | 单通错流换热器的显式传热计算方法 |
CN111814351A (zh) * | 2020-07-23 | 2020-10-23 | 清华大学 | 一种超临界压力流体换热器制备方法及系统 |
CN111814351B (zh) * | 2020-07-23 | 2023-03-21 | 清华大学 | 一种超临界压力流体换热器制备方法及系统 |
CN113361158A (zh) * | 2021-05-31 | 2021-09-07 | 江西理工大学 | 一种连续铸造结晶器内腔金属温度的测量方法及系统 |
CN113378105A (zh) * | 2021-06-24 | 2021-09-10 | 东北电力大学 | 一种超临界逆流式省煤器过程传热系数的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107194076B (zh) | 2019-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107194076A (zh) | 一种管内超临界流体气化‑管外膜状冷凝传热计算方法 | |
Pu et al. | Thermal performance analysis of intermediate fluid vaporizer for liquefied natural gas | |
Liu et al. | Thermal-hydraulic performance analysis of printed circuit heat exchanger precooler in the Brayton cycle for supercritical CO2 waste heat recovery | |
Zhao et al. | Experimental and numerical investigation of thermal-hydraulic performance of supercritical nitrogen in airfoil fin printed circuit heat exchanger | |
Kwon et al. | Experimental investigation of heat transfer coefficient of mini-channel PCHE (printed circuit heat exchanger) | |
CN103258097B (zh) | 考虑流型的用于优化设计有机朗肯循环换热器管道长度的方法 | |
Qi et al. | Performance analysis of submerged combustion vaporizer | |
Dong et al. | Experimental investigation on the heat transfer performance of molten salt flowing in an annular tube | |
Yoo et al. | Experimental study of propane condensation heat transfer and pressure drop in semicircular channel printed circuit heat exchanger | |
Wu et al. | Numerical simulation and experimental research on the comprehensive performance of the shell side of the spiral wound heat exchanger | |
Wang et al. | Multi-objective optimization design and performance evaluation of a novel multi-stream intermediate fluid vaporizer with cold energy recovery | |
Cai et al. | Numerical investigation on thermal–hydraulic performance of supercritical LNG in a Zigzag mini-channel of printed circuit heat exchanger | |
Wang et al. | Turbulent convective heat transfer of methane at supercritical pressure in a helical coiled tube | |
Jian et al. | Experimental study on the tube-side thermal-hydraulic performance of spiral wound heat exchangers | |
Zhao et al. | Experimental and numerical study on thermal‐hydraulic performance of printed circuit heat exchanger for liquefied gas vaporization | |
Moghaddam et al. | Flow pattern maps, pressure drop and performance assessment of horizontal tubes with coiled wire inserts during condensation of R-600a | |
Wang et al. | Simulation of performance of intermediate fluid vaporizer under wide operation conditions | |
Li et al. | Review of the LNG intermediate fluid vaporizer and its heat transfer characteristics | |
Wang et al. | Analysis on the heat transfer performance of supercritical liquified natural gas in horizontal tubes during regasification process | |
Jiang et al. | Adaptive design methodology of segmented non-uniform fin arrangements for trans-critical natural gas in the printed circuit heat exchanger | |
Ruan et al. | Experimental study on heat transfer in a model of submerged combustion vaporizer | |
Kubo et al. | Correlation between pressure loss and heat transfer coefficient in boiling flows in printed circuit heat exchangers with semicircular and circular mini-channels | |
Sun et al. | Research on the falling film flow and heat transfer characteristics of FLNG spiral wound heat exchanger under sea conditions | |
Wang et al. | Investigation on flow distribution of supercritical natural gas in a printed circuit heat exchanger | |
Jian et al. | Experimental investigation of the condensation mechanisms in the shell side of spiral wound heat exchangers |
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 | ||
CB02 | Change of applicant information |
Address after: 230031 Hefei Changjiang Road, Shushan District, Anhui, No. 888 Applicant after: Hefei General Machinery Research Institute Co., Ltd. Applicant after: Hefei General Machinery Research Institute special equipment inspection station Co., Ltd. Address before: 230031 Hefei Changjiang Road, Shushan District, Anhui, No. 888 Applicant before: Hefei Universal Machinery Inst. Applicant before: Special Equipment Inspection Station of Hefei General Machinery Research Institute |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |