CN113553716B - 一种多层固定的高温燃气通道设计方法 - Google Patents

一种多层固定的高温燃气通道设计方法 Download PDF

Info

Publication number
CN113553716B
CN113553716B CN202110849454.6A CN202110849454A CN113553716B CN 113553716 B CN113553716 B CN 113553716B CN 202110849454 A CN202110849454 A CN 202110849454A CN 113553716 B CN113553716 B CN 113553716B
Authority
CN
China
Prior art keywords
wall
temperature
layer
channel
ring cavity
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
CN202110849454.6A
Other languages
English (en)
Other versions
CN113553716A (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202110849454.6A priority Critical patent/CN113553716B/zh
Publication of CN113553716A publication Critical patent/CN113553716A/zh
Application granted granted Critical
Publication of CN113553716B publication Critical patent/CN113553716B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16LPIPES; JOINTS OR FITTINGS FOR PIPES; SUPPORTS FOR PIPES, CABLES OR PROTECTIVE TUBING; MEANS FOR THERMAL INSULATION IN GENERAL
    • F16L53/00Heating of pipes or pipe systems; Cooling of pipes or pipe systems
    • F16L53/70Cooling of pipes or pipe systems
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16LPIPES; JOINTS OR FITTINGS FOR PIPES; SUPPORTS FOR PIPES, CABLES OR PROTECTIVE TUBING; MEANS FOR THERMAL INSULATION IN GENERAL
    • F16L57/00Protection of pipes or objects of similar shape against external or internal damage or wear
    • F16L57/02Protection of pipes or objects of similar shape against external or internal damage or wear against cracking or buckling
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16LPIPES; JOINTS OR FITTINGS FOR PIPES; SUPPORTS FOR PIPES, CABLES OR PROTECTIVE TUBING; MEANS FOR THERMAL INSULATION IN GENERAL
    • F16L57/00Protection of pipes or objects of similar shape against external or internal damage or wear
    • F16L57/04Protection of pipes or objects of similar shape against external or internal damage or wear against fire or other external sources of extreme heat
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16LPIPES; JOINTS OR FITTINGS FOR PIPES; SUPPORTS FOR PIPES, CABLES OR PROTECTIVE TUBING; MEANS FOR THERMAL INSULATION IN GENERAL
    • F16L9/00Rigid pipes
    • F16L9/14Compound tubes, i.e. made of materials not wholly covered by any one of the preceding groups
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Fuel Cell (AREA)
  • Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)

Abstract

本发明公开了一种多层固定的高温燃气通道设计方法,具体包括:获取高温燃气通道沿程各位置所需传导的最小热流密度;基于最小热流密度确定高温燃气通道沿程各位置所需传导的实际热流密度的取值范围;基于冷却剂的流量对冷却槽的数目、槽宽与槽高进行匹配设计,并得到冷却结构的传热系数;基于高温燃气通道沿程各位置所需传导的实际热流密度、冷却结构的传热系数得到通道壁的导热热阻;选取设计参数取值以适应步骤4中通道壁的导热热阻,所述尺寸参数包括环腔的层数、各壁面的厚度、各环腔的高度、各环腔内相邻支撑肋之间的间距长度。通过对通道壁的具体尺寸进行优化设计,实现了高温流体通道在保证安全工作前提下耗费工质热沉最小。

Description

一种多层固定的高温燃气通道设计方法
技术领域
本发明涉及高温燃气通道设计技术领域,具体是一种多层固定的高温燃气通道设计方法。
背景技术
在高温燃气通道的应用过程中,考虑到通道壁材料的耐高温强度,高温燃气通道壁面一般有一个工作温度上限,当流体温度超过可用的材料温度上限时,必须采用合适的热疏导方案使得通道内壁温度工作在材料可用极限温度以下。
通常采用的可以长时间工作的热疏导方案一般为主动冷却方案,具体为采用一种冷却工质流过壁面内部合理设置的流动通道,将超出燃气通道内壁安全工作承受范围的热及时带走,保持通道内壁温度低于极限工作温度,在此过程中冷却工质温度增加。但是在实际的应用过程中,往往需要让高温燃气通道内壁尽可能工作在较为均匀且略低于其极限工作温度的条件下。这种情况下采用上述主动冷却方案,当冷却流量、允许温升、面积等不足时、冷却剂的冷却能力不足,会导致燃气通道内壁过热、冷却失败。但若是增加冷却工质质量、允许的温升、冷却面积等过大,则又会导致传导出来的热流过多、耗费了过多的工质热沉,造成能源浪费。可见在这种背景条件下,高温燃气通道的主动冷却方案需要经过优化设计,使高温燃气通道在冷却剂热沉消耗尽可能低的条件下即保证其可靠工作。
发明内容
针对上述现有技术中的不足,本发明提供一种多层固定的高温燃气通道设计方法,实现了高温燃气通道在保证安全工作前提下耗费工质热沉最小。
为实现上述目的,本发明提供一种多层固定的高温燃气通道设计方法,所述高温燃气通道包括第一连接座、第二连接座、冷却结构以及筒状结构的通道壁,所述通道壁包括内壁与套设在内壁上的外壁;所述内壁包括若干逐层套设的内壁面,所述外壁包括若干逐层套设的外壁面,所述内壁面与外壁面均为筒状结构;每相邻的两个内壁面之间、每相邻的两个外壁面之间分别围成筒状结构的内壁环腔与外壁环腔,所述内壁环腔与外壁环腔内设有支撑肋,各所述外壁面的一端与第一连接座相连,各所述外壁面的另一端与第二连接座相连;各所述内壁面的一端与第一连接座相连,各所述内壁面的另一端与第二连接座之间设有变形间隙,以用于为内壁面在受热应力变形延展时提供容纳空间;最外层内壁面与最内层的外壁面之间具有间隙环腔,该间隙环腔内设有限定结构,以使得内壁面在受热应力变形时偏向于轴向延展;所述冷却结构为沿周向间隔设在最外层的外壁面上的若干冷却槽,且最内层的内壁面上具有隔热涂层;
所述设计方法具体包括:
步骤1,获取高温燃气通道沿程各位置所需传导的最小热流密度;
步骤2,基于最小热流密度得到高温燃气通道沿程各位置所需传导的实际热流密度的取值范围;
步骤3,基于冷却剂的流量对冷却槽的数目、槽宽与槽高进行匹配设计,并得到冷却结构的传热系数;
步骤4,基于高温燃气通道沿程各位置所需传导的实际热流密度、冷却结构的传热系数得到通道壁的导热热阻;
步骤5,选取设计参数取值以适应步骤4中通道壁的导热热阻,所述设计参数包括内壁环腔的层数、各内壁面的厚度、各内壁环腔的高度、内壁环腔内相邻支撑肋之间的间距长度、外壁环腔的层数、各外壁面的厚度、各外壁环腔的高度、外壁环腔内相邻支撑肋之间的间距长度、间隙环腔的高度、间隙环腔的高度以及间隙环腔内相邻限定结构之间的间距长度;
步骤6,对步骤5中的设计结果进行有限元计算校核,如超出约束条件,则修改实际热流密度的取值范围,重复开展步骤3-5,直到满足约束条件。
在其中一个实施例中,步骤1中,
不考虑燃气辐射的情况下,所述最小热流密度的获取过程为:
qmin=hwg(Taw-Tmh)+qr (1)
式中,qmin为最小热流密度,Tmh为高温燃气通道内壁的可靠工作温度,Taw为高温燃气通道燃气表面的恢复温度,hwg为高温燃气通道内燃气与高温燃气通道内壁之间的对流传热系数,qr为燃气辐射热流;
Figure GDA0003546519490000021
式中,εw,ef为壁面有效黑度,σ0为斯蒂芬-玻尔兹曼常数,σ0=5.67×10-8W/(m2·K4),εg为高温燃气气体发射率,αw为壁面气体吸收比,Tg为高温燃气静温;
其中:
Figure GDA0003546519490000031
式中,Tc,ns为高温燃气总温,r为当地恢复系数,M为高温燃气通道沿程各位置的当地燃气流动马赫数,k为燃气混合物比热比;设某一截面A处高温燃气由液体火箭发动机燃烧产生,根据高温燃气流动马赫数M可以由迭代求解获得对应的假想发动机喉部面积At
Figure GDA0003546519490000032
进而高温燃气与壁面的对流换热系数hwg可以由Bartz公式获得:
Figure GDA0003546519490000033
式中,Dt当量喉部直径,
Figure GDA0003546519490000034
cpg为燃气定压比热,μg为燃气动力粘度,Prg为燃气的普朗特数,pc.ns为高温燃气总压,c*为特征速度,可由火箭发动机热力计算获得,R为喉部处喷管的曲率半径,该项影响很小,在此处可以取R=Dt,σ为修正系数;
若高温燃气通道内壁有涂层,则需要将上述公式中Tmh更换为涂层燃气侧温度Ttcg,进一步
Figure GDA0003546519490000035
式中,δtc与λtc分别为隔热涂层的厚度与热导率,此情况下,联立式(1)和(6),并采用迭代求解方法获得Ttcg,进而求得qmin
在其中一个实施例中,所述修正系数的计算过程为:
Figure GDA0003546519490000036
在其中一个实施例中,步骤3中,基于冷却剂的流量对冷却槽的数目、槽宽与槽高进行匹配设计,具体为:
获取冷却剂的流速VL,其中,冷却剂的流速VL的取值范围为1~10m/s,并基于冷却剂的流速VL得到各冷却槽的总流通面积AL
Figure GDA0003546519490000041
式中,ρL为冷却剂的密度,
Figure GDA0003546519490000042
为冷却剂的流量;
在此基础上:
AL=NLWLHL (9)
式中,NL、WL、HL为冷却槽的数目、槽宽、槽高,选取NL、WL、HL的具体取值使其满足总流通面积AL即可。
在其中一个实施例中,步骤3中,所述冷却结构的传热系数的计算过程为:
Figure GDA0003546519490000043
式中,hwL为冷却结构的传热系数,ρL为冷却剂密度,λL为冷却剂导热系数,PrL为冷却剂的普朗特数,VL为冷却剂的流速,μL为冷却剂的粘度,DL为冷却槽槽道的水力直径,对于矩形通道,则有:
Figure GDA0003546519490000044
在其中一个实施例中,步骤4,所述基于高温燃气通道沿程各位置所需传导的实际热流密度、冷却结构的传热系数得到通道壁的导热热阻,具体为:
总导热热阻、高温燃气通道内外壁温差与实际热流密度的关系为:
Figure GDA0003546519490000045
推导得到:
Figure GDA0003546519490000046
式中,q为高温燃气通道的实际热流密度,Rw表示通道壁的导热热阻,TL表示冷却剂温度,
Figure GDA0003546519490000047
表示冷却剂侧的对流换热热阻,Rg为燃气侧综合热阻:
Figure GDA0003546519490000048
其中,燃气侧综合热阻Rg包含了燃气侧的对流换热热阻、辐射热阻和涂层的热阻的综合作用。
在其中一个实施例中,步骤5中,选取设计参数取值以适应步骤4中通道壁的导热热阻,所述设计参数包括内壁环腔的层数、各内壁面的厚度、各内壁环腔的高度、内壁环腔内相邻支撑肋之间的间距长度、外壁环腔的层数、各外壁面的厚度、各外壁环腔的高度、外壁环腔内相邻支撑肋之间的间距长度,具体为:
在给定燃气侧参数、冷却剂侧参数、通道壁物性参数和壁面辐射参数等的情况下,可以通过有限元方法获得燃气通道壁温度分布和导热热阻Rw
在初步设计时,假设各内壁环腔、外壁环腔内外侧壁温一致,通过如下简化计算公式进行计算,待确定优化参数后,再细化校核,最终确定通道壁的各设计参数,简化计算方法具体为:
对于具有m个内壁环腔、n个外壁环腔的通道壁,其导热热阻主要包括m+1个内壁面的导热热阻,以及m个内壁环腔的辐射热阻、n+1个外壁面的导热热阻,以及n+1个内壁环腔的辐射热阻、最内层外壁面与最外层内壁面之间的辐射热阻,即间隙环腔的辐射热阻,因此,当m=1、n=1时,通道壁的导热热阻R1w需满足下式:
Figure GDA0003546519490000051
Figure GDA0003546519490000052
Figure GDA0003546519490000053
Figure GDA0003546519490000061
式中,λw表示壁面材料热导率;δa1表示第一层内壁面的壁厚,Ha1表示第一层内壁环腔的高度,La1表示第一层内壁环腔内相邻支撑肋之间的间距长度,δaf表示最后一层内壁面的壁厚,Ta1S表示第一层内壁环腔下表面的温度,Ta1X表示第一层内壁环腔上表面的温度;δb1表示第一层外壁面的壁厚,Hb1表示第一层外壁环腔的高度,Lb1表示第一层外壁环腔内相邻支撑肋之间的间距长度,δbf表示最后一层外壁面的壁厚,Tb1S表示第一层外壁环腔下表面的温度,Tb1X表示第一层内壁环腔上表面的温度;Hc表示间隙环腔的高度,Lc表示间隙环腔内相邻限定结构之间的间距长度,TcS表示间隙环腔下表面的温度,即最后一层内壁面的外壁温度;TcX表示间隙环腔上表面的温度,即第一层外壁面的内壁温度;ε为内壁面与外壁面材料的表面黑度;
此后,每当m的数量增加1时,在上式(14)的右边项增加如下参数:
Figure GDA0003546519490000062
Figure GDA0003546519490000063
式中,δai表示第i层内壁面的壁厚,Hai表示第i层内壁环腔的高度,Lai表示第i层内壁环腔内相邻支撑肋之间的间距长度,Tais表示第i层内壁环腔下表面的温度,Taix表示第i层内壁环腔上表面的温度;
每当n的数量增加1时,在上式(14)的右边项增加如下参数:
Figure GDA0003546519490000064
Figure GDA0003546519490000065
式中,δbj表示第j层外壁面的壁厚,Hbj表示第j层外壁环腔的高度,Lbj表示第j层外壁环腔内相邻支撑肋之间的间距长度,TbjS表示第j层外壁环腔下表面的温度,TbjX表示第j层外壁环腔上表面的温度;
最终:
Figure GDA0003546519490000071
并基于仿真或数值模拟,以结构质量最轻为优化目标、以热应力需求、壁面最大温度、冷却剂流量和温升为约束,选取m、δai、δaf、Hai、Lai、n、δbj、δbf、Hbj、Lbj、Hc、Lc的具体取值,其中,i=1,···,m、j=1,···,n。
本发明提供的一种多层固定的高温燃气通道设计方法,其通过将若干筒状结构的壁面逐层套设而形成高温燃气通道的通道壁,并在相邻的壁面之间形成环腔,使得由通道内壁传递至通道外壁的热流大幅度降低,并通过对通道壁的具体尺寸进行优化设计,实现了高温流体通道在保证安全工作前提下耗费工质热沉最小。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图示出的结构获得其他的附图。
图1为本发明实施例中高温燃气通道壁面的结构示意图;
图2为本发明实施例中支撑肋的第一种实施结构示意图;
图3为本发明实施例中支撑肋的第二种实施结构示意图;
图4为本发明实施例中支撑肋的第三种实施结构示意图;
图5为本发明实施例中支撑肋的第四种实施结构示意图;
图6为本发明实施例中支撑肋的第五种实施结构示意图;
图7为本发明实施例中支撑肋的第六种实施结构示意图;
图8为本发明实施例中多层固定的高温燃气通道设计方法的流程示意图。
附图标号:内壁面1、内壁环腔2、支撑肋3、冷却槽4、肋板5、外壁面6、外壁环腔7、间隙环腔8、限定结构9。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明,本发明实施例中所有方向性指示(诸如上、下、左、右、前、后……)仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等,如果该特定姿态发生改变时,则该方向性指示也相应地随之改变。
另外,在本发明中如涉及“第一”、“第二”等的描述仅用于描述目的,而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本发明中,除非另有明确的规定和限定,术语“连接”、“固定”等应做广义理解,例如,“固定”可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接,还可以是物理连接或无线通信连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
另外,本发明各个实施例之间的技术方案可以相互结合,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在,也不在本发明要求的保护范围之内。
如图1所示,本实施例公开了一种高温燃气通道,该高温燃气通道包括第一连接座、第二连接座、冷却结构以及筒状结构的通道壁,通道壁包括内壁与套设在内壁上的外壁;内壁包括若干逐层套设的内壁面1,外壁包括若干逐层套设的外壁面6,内壁面1与外壁面6均为筒状结构;每相邻的两个内壁面1之间、每相邻的两个外壁面6之间分别围成筒状结构的内壁环腔2与外壁环腔7,内壁环腔2与外壁环腔7内设有支撑肋3,各外壁面6的一端与第一连接座相连,各外壁面6的另一端与第二连接座相连;各内壁面1的一端与第一连接座相连,各内壁面1的另一端与第二连接座之间设有变形间隙,以用于为内壁面在受热应力变形延展时提供容纳空间;最外层内壁面与最内层的外壁面之间具有间隙环腔8,该间隙环腔内设有限定结构9,以使得内壁面在受热应力变形时偏向于轴向延展;冷却结构为沿周向间隔设在最外层的外壁面上的若干冷却槽4,且最内层的内壁面上具有隔热涂层。通过将通道壁设置为多层壁面1的结构,使得由通道壁内壁传递至通道壁外壁热流大幅度降低,实现了高温燃气通道在保证安全工作前提下耗费工质热沉最小。本实施例中,高温燃气通道内壁一般为可采用增材制造等先进制造技术制造的材料,其与燃气可以直接接触,也可以额外增加隔热涂层或者抗氧化涂层等。需要注意的是,上述的高温燃气通道中的支撑肋3既可以是图1所示的实心的肋条结构,也可以是图2-3所示的空心的圆管或方管结构,也可以是图4-7所示的实心的肋条与空心的圆管或方管的组合结构。
基于上述结构的高温燃气通道,参考图8,本实施例公开了一种多层固定的高温燃气通道设计方法,该设计方法具体包括:
步骤1,获取高温燃气通道沿程各位置所需传导的最小热流密度。
对于高温燃气通道沿程任意位置的所需传导的最小热流密度,在不考虑燃气辐射的情况下,其获取过程为:
qmin=hwg(Taw-Tmh)+qr (1)
式中,qmin为最小热流密度;Tmh为高温燃气通道内壁的可靠工作温度,由壁面的材料特性确定,单位为K;Taw为高温燃气通道燃气表面的恢复温度,单位为K;hwg为高温燃气通道内燃气与高温燃气通道内壁之间的对流传热系数,单位为W·m-2·K-1,qr为燃气辐射热流;
Figure GDA0003546519490000091
式中,εw,ef为壁面有效黑度,σ0为斯蒂芬-玻尔兹曼常数,σ0=5.67×10-8W/(m2·K4),εg为高温燃气气体发射率,αw为壁面气体吸收比,Tg为高温燃气静温;
其中:
Figure GDA0003546519490000092
式中,Tc,ns为高温燃气总温,M为高温燃气通道沿程各位置的当地燃气流动马赫数,k为燃气混合物比热比;r为当地恢复系数,表示为摩擦产生的温度增量与绝热压缩引起的温度增量之比,可由试验确定或者根据普朗特数的下列简化关系式估算,即分别为层流和湍流对应的关系式:
Figure GDA0003546519490000093
设某一截面A处高温燃气由液体火箭发动机燃烧产生,根据高温燃气流动马赫数M可以由迭代求解获得对应的假想发动机喉部面积At
Figure GDA0003546519490000101
进而高温燃气与壁面的对流换热系数hwg可以由Bartz公式获得:
Figure GDA0003546519490000102
式中,Dt当量喉部直径,
Figure GDA0003546519490000103
cpg为燃气定压比热,μg为燃气动力粘度,Prg为燃气的普朗特数,pc.ns为高温燃气总压,c*为特征速度,可由火箭发动机热力计算获得,R为喉部处喷管的曲率半径,该项影响很小,在此处可以取R=Dt,σ为修正系数;
若高温燃气通道内壁有涂层,则需要将上述公式中Tmh更换为涂层燃气侧温度Ttcg,进一步
Figure GDA0003546519490000104
式中,δtc与λtc分别为隔热涂层的厚度与热导率,此情况下,联立式(1)和(6),并采用迭代求解方法获得Ttcg,进而求得qmin
上述修正系数σ的计算过程为:
Figure GDA0003546519490000105
步骤2,基于最小热流密度确定高温燃气通道沿程各位置所需传导的实际热流密度的取值范围,具体为:
对于高温燃气通道沿程任一位置,其最小热流密度为qmin,其实际热流密度一般比qmin略大,例如,在得到最小热流密度为qmin,可以将实际热流密度q的初始取值范围定位:1.05qmin≤q≤1.5qmin,在后续步骤6的校核过程中,基于校核结构调整实际热流密度q的初始取值范围。
步骤3,基于冷却剂的流量对冷却槽的数目、槽宽与槽高进行匹配设计,并得到冷却结构的传热系数。
一般来说可选的流速范围较宽,流速初值可取1~10m/s。但从设计角度看,流速过高意味着冷却剂的流通面积较小、压力损失较大;而流速过低则意味着冷却剂的流通面积较大、增加了加工成本。因此出于优化设计的目的,通常建设选择中间流速VL=5m/s作为初值,并基于冷却剂的流速VL得到各冷却槽的总流通面积AL
Figure GDA0003546519490000111
式中,ρL为冷却剂的密度,
Figure GDA0003546519490000112
为冷却剂的流量;
在此基础上:
AL=NLWLHL(9)
式中,NL、WL、HL为冷却槽的数目、槽宽、槽高,选取NL、WL、HL的具体取值使其满足总流通面积AL即可。本实施例中,槽宽可选的范围较宽,但考虑到承压与加工,因此槽宽的取值为1~6mm。同时需要注意的是,槽与槽之间存在肋板,肋板宽度取值一般不超过槽宽、另外考虑到结构强度也不宜过窄,因此限定肋板宽度的取值为0.4~3mm,在确定槽宽与肋板宽度后,即能确定冷却槽的数目以及槽高。
而冷却结构的传热系数的计算过程为:
Figure GDA0003546519490000113
式中,hwL为冷却结构的传热系数,ρL为冷却剂密度,λL为冷却剂导热系数,PrL为冷却剂的普朗特数,VL为冷却剂的流速,μL为冷却剂的粘度,DL为冷却槽槽道的水力直径,对于矩形通道,则有:
Figure GDA0003546519490000114
步骤4,基于高温燃气通道沿程各位置所需传导的实际热流密度、冷却结构的传热系数得到通道壁的导热热阻,具体为:
总导热热阻、高温燃气通道内外壁温差与实际热流密度的关系为:
Figure GDA0003546519490000115
推导得到:
Figure GDA0003546519490000121
式中,q为高温燃气通道的实际热流密度,Rw表示通道壁的导热热阻,TL表示冷却剂温度,
Figure GDA0003546519490000122
表示冷却剂侧的对流换热热阻,Rg为燃气侧综合热阻:
Figure GDA0003546519490000123
其中,燃气侧综合热阻Rg包含了燃气侧的对流换热热阻、辐射热阻和涂层的热阻的综合作用。
步骤5,选取设计参数取值以适应步骤4中通道壁的导热热阻,设计参数包括环腔的层数、各壁面的厚度、各环腔的高度、各环腔内相邻支撑肋之间的间距长度,具体为:
在给定燃气侧参数、冷却剂侧参数、通道壁物性参数和壁面辐射参数等的情况下,可以通过有限元方法获得燃气通道壁温度分布和导热热阻Rw,其中,燃气侧参数包括对流换热系数与恢复温度,冷却剂侧参数包括对流换热系数与冷却剂温度;
在初步设计时,假设各环腔内外侧壁温一致,分别为Tns和TnX,通过如下简化计算公式进行计算,待确定优化参数后,再细化校核,最终确定通道壁的各设计参数,简化计算方法具体为:
对于具有m个内壁环腔、n个外壁环腔的通道壁,其导热热阻主要包括m+1个内壁面的导热热阻,以及m个内壁环腔的辐射热阻、n+1个外壁面的导热热阻,以及n+1个内壁环腔的辐射热阻、最内层外壁面与最外层内壁面之间的辐射热阻,即间隙环腔的辐射热阻,因此,当m=1、n=1时,通道壁的导热热阻R1w需满足下式:
Figure GDA0003546519490000131
Figure GDA0003546519490000132
Figure GDA0003546519490000133
Figure GDA0003546519490000134
式中,λaB、λbB表示表示壁面材料热导率;δa1表示第一层内壁面的壁厚,Ha1表示第一层内壁环腔的高度,La1表示第一层内壁环腔内相邻支撑肋之间的间距长度,δaf表示最后一层内壁面的壁厚,Ta1S表示第一层内壁环腔下表面的温度,Ta1X表示第一层内壁环腔上表面的温度;δb1表示第一层外壁面的壁厚,Hb1表示第一层外壁环腔的高度,Lb1表示第一层外壁环腔内相邻支撑肋之间的间距长度,δbf表示最后一层外壁面的壁厚,Tb1S表示第一层外壁环腔下表面的温度,Tb1X表示第一层内壁环腔上表面的温度;Hc表示间隙环腔的高度,Lc表示间隙环腔内相邻限定结构之间的间距长度,TcS表示间隙环腔下表面的温度,即最后一层内壁面的外壁温度;TcX表示间隙环腔上表面的温度,即第一层外壁面的内壁温度;ε为内壁面与外壁面材料的表面黑度;本实施例中,ε=0.8。
其中,对于第一层内壁环腔下表面的温度T1S,其计算过程为:
Figure GDA0003546519490000135
对于第一层内壁环腔上表面的温度Ta1X,则要通过迭代或者数值方法计算得到,其过程为常规技术手段,因此本实施例中不再赘述。对于,Tais与Taix的计算过程,采用与第一层内壁环腔中Ta1S与Ta1X通用的计算方式即可,只需将Tais的计算过程中,令Taw=T(i-1)x即可,同样,也可以采用同样的原理计算得到TcS、TcX、Tbjs与Tbjx
式(14)中的第一项为内壁面的导热热阻,第二项为外壁面的导热热阻,第三项为间隙环腔的辐射热阻;此后,每当m的数量增加1时,在上式(14)的右边项增加如下参数:
Figure GDA0003546519490000141
Figure GDA0003546519490000142
式中,δai表示第i层内壁面的壁厚,Hai表示第i层内壁环腔的高度,Lai表示第i层内壁环腔内相邻支撑肋之间的间距长度,Tais表示第i层内壁环腔下表面的温度,Taix表示第i层内壁环腔上表面的温度;
每当n的数量增加1时,在上式(14)的右边项增加如下参数:
Figure GDA0003546519490000143
Figure GDA0003546519490000144
式中,δbj表示第j层外壁面的壁厚,Hbj表示第j层外壁环腔的高度,Lbj表示第j层外壁环腔内相邻支撑肋之间的间距长度,TbjS表示第j层外壁环腔下表面的温度,TbjX表示第j层外壁环腔上表面的温度;
最终:
Figure GDA0003546519490000145
并基于仿真或数值模拟,以结构质量最轻为优化目标、以热应力需求、壁面最大温度、冷却剂流量和温升为约束,选取m、δai、δaf、Hai、Lai、n、δbj、δbf、Hbj、Lbj、Hc、Lc的具体取值,其中,i=1,···,m、j=1,···,n。至于变形间隙的长度,在m、δai、δaf、Hai、Lai、n、δbj、δbf、Hbj、Lbj设计完成后即能内壁面受到热应力变形的最大幅度,使得变形间隙的长度略大于该最大幅度值即可。
在仿真或数值模拟的过程中,最靠燃气侧这层的厚度可以适当取厚,即第一层内壁面的厚度以0.5~1.5mm为宜,而其它层内壁面与外壁面的厚度则可取0.3~1mm,对于各内壁环腔与外壁环腔的高度,其取值范围为1~4mm;各内壁环腔与外壁环腔内相邻支撑肋之间的间距长度,其取值范围为2~10mm。
步骤6,在完成上述各尺寸与数值的设计后,即完成高温燃气通道初步结构设计,在初步结构设计的基础上进行有限元计算校核,如超出约束条件,则修改实际热流密度的取值范围后重复步骤3-5,其中,约束条件为热应力需求、壁面最大温度、冷却剂流量和温升。若校核结构为超出约束条件,则扩大实际热流密度的取值范围,例如将实际热流密度的取值范围由1.05qmin~1.5qmin修改为1.5qmin~1.7qmin。若校核结果还具有余量,则缩小实际热流密度的取值范围,例如将实际热流密度的取值范围由1.05qmin~1.5qmin修改为1.05qmin~1.3qmin
以上所述仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是在本发明的发明构思下,利用本发明说明书及附图内容所作的等效结构变换,或直接/间接运用在其他相关的技术领域均包括在本发明的专利保护范围内。

Claims (7)

1.一种多层固定的高温燃气通道设计方法,其特征在于,所述高温燃气通道包括第一连接座、第二连接座、冷却结构以及筒状结构的通道壁,所述通道壁包括内壁与套设在内壁上的外壁;所述内壁包括若干逐层套设的内壁面,所述外壁包括若干逐层套设的外壁面,所述内壁面与外壁面均为筒状结构;每相邻的两个内壁面之间、每相邻的两个外壁面之间分别围成筒状结构的内壁环腔与外壁环腔,所述内壁环腔与外壁环腔内设有支撑肋,各所述外壁面的一端与第一连接座相连,各所述外壁面的另一端与第二连接座相连;各所述内壁面的一端与第一连接座相连,各所述内壁面的另一端与第二连接座之间设有变形间隙,以用于为内壁面在受热应力变形延展时提供容纳空间;最外层内壁面与最内层的外壁面之间具有间隙环腔,该间隙环腔内设有限定结构,以使得内壁面在受热应力变形时偏向于轴向延展;所述冷却结构为沿周向间隔设在最外层的外壁面上的若干冷却槽,且最内层的内壁面上具有隔热涂层;
所述设计方法具体包括:
步骤1,获取高温燃气通道沿程各位置所需传导的最小热流密度;
步骤2,基于最小热流密度得到高温燃气通道沿程各位置所需传导的实际热流密度的取值范围;
步骤3,基于冷却剂的流量对冷却槽的数目、槽宽与槽高进行匹配设计,并得到冷却结构的传热系数;
步骤4,基于高温燃气通道沿程各位置所需传导的实际热流密度、冷却结构的传热系数得到通道壁的导热热阻;
步骤5,选取设计参数取值以适应步骤4中通道壁的导热热阻,所述设计参数包括内壁环腔的层数、各内壁面的厚度、各内壁环腔的高度、内壁环腔内相邻支撑肋之间的间距长度、外壁环腔的层数、各外壁面的厚度、各外壁环腔的高度、外壁环腔内相邻支撑肋之间的间距长度、间隙环腔的高度、间隙环腔的高度以及间隙环腔内相邻限定结构之间的间距长度;
步骤6,对步骤5中的设计结果进行有限元计算校核,如超出约束条件,则修改实际热流密度的取值范围,重复开展步骤3-5,直到满足约束条件。
2.根据权利要求1所述多层固定的高温燃气通道设计方法,其特征在于,步骤1中,
不考虑燃气辐射的情况下,所述最小热流密度的获取过程为:
qmin=hwg(Taw-Tmh)+qr (1)
式中,qmin为最小热流密度,Tmh为高温燃气通道内壁的可靠工作温度,Taw为高温燃气通道燃气表面的恢复温度,hwg为高温燃气通道内燃气与高温燃气通道内壁之间的对流传热系数,qr为燃气辐射热流;
Figure FDA0003546519480000021
式中,εw,ef为壁面有效黑度,σ0为斯蒂芬-玻尔兹曼常数,σ0=5.67×10-8W/(m2·K4),εg为高温燃气气体发射率,αw为壁面气体吸收比,Tg为高温燃气静温;
其中:
Figure FDA0003546519480000022
式中,Tc,ns为高温燃气总温,r为当地恢复系数,M为高温燃气通道沿程各位置的当地燃气流动马赫数,k为燃气混合物比热比;设某一截面A处高温燃气由液体火箭发动机燃烧产生,根据高温燃气流动马赫数M可以由迭代求解获得对应的假想发动机喉部面积At
Figure FDA0003546519480000023
进而高温燃气与壁面的对流换热系数hwg可以由Bartz公式获得:
Figure FDA0003546519480000024
式中,Dt当量喉部直径,
Figure FDA0003546519480000025
cpg为燃气定压比热,μg为燃气动力粘度,Prg为燃气的普朗特数,pc.ns为高温燃气总压,c*为特征速度,可由火箭发动机热力计算获得,R为喉部处喷管的曲率半径,该项影响很小,在此处可以取R=Dt,σ为修正系数;
若高温燃气通道内壁有涂层,则需要将上述公式中Tmh更换为涂层燃气侧温度Ttcg,进一步
Figure FDA0003546519480000026
式中,δtc与λtc分别为隔热涂层的厚度与热导率,此情况下,联立式(1)和(6),并采用迭代求解方法获得Ttcg,进而求得qmin
3.根据权利要求2所述多层固定的高温燃气通道设计方法,其特征在于,所述修正系数的计算过程为:
Figure FDA0003546519480000031
4.根据权利要求3所述多层固定的高温燃气通道设计方法,其特征在于,步骤3中,基于冷却剂的流量对冷却槽的数目、槽宽与槽高进行匹配设计,具体为:
获取冷却剂的流速VL,其中,冷却剂的流速VL的取值范围为1~10m/s,并基于冷却剂的流速VL得到各冷却槽的总流通面积AL
Figure FDA0003546519480000032
式中,ρL为冷却剂的密度,
Figure FDA0003546519480000033
为冷却剂的流量;
在此基础上:
AL=NLWLHL (9)
式中,NL、WL、HL为冷却槽的数目、槽宽、槽高,选取NL、WL、HL的具体取值使其满足总流通面积AL即可。
5.根据权利要求3所述高温燃气通道设计方法,其特征在于,步骤3中,所述冷却结构的传热系数的计算过程为:
Figure FDA0003546519480000034
式中,hwL为冷却结构的传热系数,ρL为冷却剂密度,λL为冷却剂导热系数,PrL为冷却剂的普朗特数,VL为冷却剂的流速,μL为冷却剂的粘度,DL为冷却槽槽道的水力直径,对于矩形通道,则有:
Figure FDA0003546519480000035
6.根据权利要求5所述多层固定的高温燃气通道设计方法,其特征在于,步骤4,所述基于高温燃气通道沿程各位置所需传导的实际热流密度、冷却结构的传热系数得到通道壁的导热热阻,具体为:
总导热热阻、高温燃气通道内外壁温差与实际热流密度的关系为:
Figure FDA0003546519480000041
推导得到:
Figure FDA0003546519480000042
式中,q为高温燃气通道的实际热流密度,Rw表示通道壁的导热热阻,TL表示冷却剂温度,
Figure FDA0003546519480000043
表示冷却剂侧的对流换热热阻,Rg为燃气侧综合热阻:
Figure FDA0003546519480000044
其中,燃气侧综合热阻Rg包含了燃气侧的对流换热热阻、辐射热阻和涂层的热阻的综合作用。
7.根据权利要求1所述多层固定的高温燃气通道设计方法,其特征在于,步骤5中,选取设计参数取值以适应步骤4中通道壁的导热热阻,所述设计参数包括内壁环腔的层数、各内壁面的厚度、各内壁环腔的高度、内壁环腔内相邻支撑肋之间的间距长度、外壁环腔的层数、各外壁面的厚度、各外壁环腔的高度、外壁环腔内相邻支撑肋之间的间距长度,具体为:
在给定燃气侧参数、冷却剂侧参数、通道壁物性参数和壁面辐射参数的情况下,可以通过有限元方法获得燃气通道壁温度分布和导热热阻Rw
在初步设计时,假设各内壁环腔、外壁环腔内外侧壁温一致,通过如下简化计算公式进行计算,待确定优化参数后,再细化校核,最终确定通道壁的各设计参数,简化计算方法具体为:
对于具有m个内壁环腔、n个外壁环腔的通道壁,其导热热阻主要包括m+1个内壁面的导热热阻,以及m个内壁环腔的辐射热阻、n+1个外壁面的导热热阻,以及n+1个内壁环腔的辐射热阻、最内层外壁面与最外层内壁面之间的辐射热阻,即间隙环腔的辐射热阻,因此,当m=1、n=1时,通道壁的导热热阻R1w需满足下式:
Figure FDA0003546519480000051
Figure FDA0003546519480000052
Figure FDA0003546519480000053
Figure FDA0003546519480000054
式中,λw表示壁面材料热导率;δa1表示第一层内壁面的壁厚,Ha1表示第一层内壁环腔的高度,La1表示第一层内壁环腔内相邻支撑肋之间的间距长度,δaf表示最后一层内壁面的壁厚,Ta1S表示第一层内壁环腔下表面的温度,Ta1X表示第一层内壁环腔上表面的温度;δb1表示第一层外壁面的壁厚,Hb1表示第一层外壁环腔的高度,Lb1表示第一层外壁环腔内相邻支撑肋之间的间距长度,δbf表示最后一层外壁面的壁厚,Tb1S表示第一层外壁环腔下表面的温度,Tb1X表示第一层内壁环腔上表面的温度;Hc表示间隙环腔的高度,Lc表示间隙环腔内相邻限定结构之间的间距长度,TcS表示间隙环腔下表面的温度,即最后一层内壁面的外壁温度;TcX表示间隙环腔上表面的温度,即第一层外壁面的内壁温度;ε为内壁面与外壁面材料的表面黑度;
此后,每当m的数量增加1时,在上式(14)的右边项增加如下参数:
Figure FDA0003546519480000055
Figure FDA0003546519480000061
式中,δai表示第i层内壁面的壁厚,Hai表示第i层内壁环腔的高度,Lai表示第i层内壁环腔内相邻支撑肋之间的间距长度,Tais表示第i层内壁环腔下表面的温度,Taix表示第i层内壁环腔上表面的温度;
每当n的数量增加1时,在上式(14)的右边项增加如下参数:
Figure FDA0003546519480000062
Figure FDA0003546519480000063
式中,δbj表示第j层外壁面的壁厚,Hbj表示第j层外壁环腔的高度,Lbj表示第j层外壁环腔内相邻支撑肋之间的间距长度,TbjS表示第j层外壁环腔下表面的温度,TbjX表示第j层外壁环腔上表面的温度;
最终:
Figure FDA0003546519480000064
并基于仿真或数值模拟,以结构质量最轻为优化目标、以热应力需求、壁面最大温度、冷却剂流量和温升为约束,选取m、δai、δaf、Hai、Lai、n、δbj、δbf、Hbj、Lbj、Hc、Lc的具体取值,其中,i=1,…,m、j=1,…,n。
CN202110849454.6A 2021-07-27 2021-07-27 一种多层固定的高温燃气通道设计方法 Active CN113553716B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110849454.6A CN113553716B (zh) 2021-07-27 2021-07-27 一种多层固定的高温燃气通道设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110849454.6A CN113553716B (zh) 2021-07-27 2021-07-27 一种多层固定的高温燃气通道设计方法

Publications (2)

Publication Number Publication Date
CN113553716A CN113553716A (zh) 2021-10-26
CN113553716B true CN113553716B (zh) 2022-05-10

Family

ID=78132888

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110849454.6A Active CN113553716B (zh) 2021-07-27 2021-07-27 一种多层固定的高温燃气通道设计方法

Country Status (1)

Country Link
CN (1) CN113553716B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114580217B (zh) * 2022-05-07 2022-07-26 西南交通大学 一种用于准环对称仿星器的限制器及其设计方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010281525A (ja) * 2009-06-05 2010-12-16 Mitsubishi Heavy Industries Mechatronics Systems Ltd 排ガス熱回収システム
US10094285B2 (en) * 2011-12-08 2018-10-09 Siemens Aktiengesellschaft Gas turbine outer case active ambient cooling including air exhaust into sub-ambient cavity
WO2015006026A1 (en) * 2013-07-12 2015-01-15 United Technologies Corporation Gas turbine engine component cooling with resupply of cooling passage
ES2543366B1 (es) * 2014-02-15 2016-05-25 Technokontrol Global Ltd Procedimiento para transformar materiales inorgánicos, NFU, petróleo crudo en hidrocarburos y/o energías renovables limpias mediante un sistema de termólisis
KR101772837B1 (ko) * 2014-04-25 2017-08-29 미츠비시 히타치 파워 시스템즈 가부시키가이샤 가스터빈 연소기 및 해당 연소기를 구비한 가스터빈
KR101863779B1 (ko) * 2017-09-15 2018-06-01 두산중공업 주식회사 라이너 냉각을 촉진하는 나선형 구조 및 이를 포함하는 가스 터빈용 연소기
JP2019173600A (ja) * 2018-03-27 2019-10-10 株式会社Soken 内燃機関の排気浄化装置
CN108825404B (zh) * 2018-06-12 2019-06-25 中国人民解放军国防科技大学 一种组合发动机燃烧室内多股流动混合燃烧的计算方法
CN109752187B (zh) * 2018-12-12 2020-10-13 西安航天动力试验技术研究所 姿轨控发动机真空环境高速高温燃气快速升压降温系统
US11423198B2 (en) * 2020-01-02 2022-08-23 Viettel Group Aerodynamic designing method of the central section of small size gas turbine engines
CN213237519U (zh) * 2020-05-15 2021-05-18 苏州万硕联信息技术有限公司 一种用于高温低氧环境的燃烧模块以及燃烧器
CN112832928B (zh) * 2021-03-05 2022-04-22 中国科学院力学研究所 一种用于火箭发动机的等内壁强度的冷却结构设计方法

Also Published As

Publication number Publication date
CN113553716A (zh) 2021-10-26

Similar Documents

Publication Publication Date Title
CN113357054B (zh) 一种高温燃气通道设计方法
Bailey et al. Local heat transfer and flow distributions for impinging jet arrays of dense and sparse extent
Kumar et al. Numerical study of heat transfer in a finned double pipe heat exchanger
Şahin The effect of variable viscosity on the entropy generation and pumping power in a laminar fluid flow through a duct subjected to constant heat flux
EP0048021B1 (en) Heat transfer device having an augmented wall surface
Tang et al. Thermal performance enhancement of an ultra-thin flattened heat pipe with multiple wick structure
CN113553716B (zh) 一种多层固定的高温燃气通道设计方法
Afzal et al. Steady and transient state analyses on conjugate laminar forced convection heat transfer
CN112832929B (zh) 一种用于火箭发动机的等内壁面温度的冷却结构设计方法
CN112832928B (zh) 一种用于火箭发动机的等内壁强度的冷却结构设计方法
El-Jummah et al. Conjugate heat transfer computational fluid dynamic predictions of impingement heat transfer: the influence of hole pitch to diameter ratio X/D at constant impingement gap Z
WO2002055864A1 (en) Rocket engine member and a method for manufacturing a rocket engine member
US7148450B2 (en) Portable blackbody furnace
Sambamurthy et al. Array jet impingement on high porosity thin metal foams: effect of foam height, pore-density, and spent air crossflow scheme on flow distribution and heat transfer
CN113389659B (zh) 一种低热沉的高温燃气通道设计方法
CN111428184A (zh) 一种板翅式换热器芯体尺寸计算方法
Nayak et al. Natural convection heat transfer in heated verticaltubes with internal rings
Yao et al. Numerical and experimental investigations of the heat transfer and fatigue life of a new W-type mini microchannel heat exchanger
Rang et al. Effects of surface roughness on forced convection and friction in triangular ducts
Kilic A heat transfer analysis from a porous plate with transpiration cooling
WO2021242131A1 (en) A cooled prandtl probe assembly
Jassim et al. An experimental investigation in forced convective heat transfer and friction factor of air flow over aligned round and flattened tube banks
Kundu et al. Exact and approximate analytic methods to calculate maximum heat flow in annular fin arrays with a rectangular step profile
CN110030862A (zh) 热管理系统及方法
Liao et al. Augmentation of turbulent convective heat transfer in tubes with three-dimensional internal extended surfaces

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