基于Fluent石蜡相变材料模拟
- 格式:docx
- 大小:334.30 KB
- 文档页数:19
石蜡/蛭石复合相变储能材料的强化传热研究边鹏旭打宋春燕",崔肖打李金洪2,刘利洁2,黄凯越2(1.石河子大学理学院低维材料物性及器件物理实验室,新疆石河子832003;2.中国地质大学(北京)材料科学与工程学院,北京100083)摘要:石蜡/蛭石复合相变储能材料(PEVC)具有潜在储热优势和应用前景。
文章采用数值模拟对PEVC进行强化传热特性分析,获得相变传热过程。
理论分析发现,PEVC整体处于吸热过程,当未达到石蜡熔点,PEVC未相变时,温度、热流变化较快;达到石蜡熔点时,材料呈糊状,温度、热流变化慢;当继续升高温度,石蜡发生完全相变,温度波动变小,热流恒定。
符合相变材料潜热高、导热系数大、密度大等特点。
总体上数值模拟结果与实验结果相一致,本研究有望为实现石蜡类复合相变储能材料在建筑中的应用提供物理依据。
关键词:石蜡/蛭石复合相变材料;数值模拟;传热中图分类号:TB34文献标志码:A文章编号:2095-2945(2021)16-0028-05Abstract:The paraffin/vermiculite composite phase change energy storage material(PEVC)has potential heat storage advantages and application prospects,and has attracted increasing attention.In this paper,numerical simulation software is used to analyze the enhanced heat transfer characteristics of PEVC to obtain the phase change heat transfer process.Theoretical analysis found that PEVC is in an endothermic process.When the melting point of paraffin is not reached and the phase transition of PEVC does not change,the temperature and heat flow change rapidly;when the melting point of paraffin is reached,the material is pasty,and the temperature and heat flow change slowly;when it continues to rise at high temperatures,the paraffin phase changes completely,the temperature fluctuation becomes smaller,and the heat flow is constant.It is consistent with the characteristics of high latent heat,large thermal conductivity and high density of phase change materials.The simulation results of numerical simulation software are consistent with the experimental results.This study is expected to provide a physical basis for the application of paraffin-based composite phase change energy storage materials in construction.Keywords:paraffin/vermiculite composite phase change material;numerical simulation;heat transfer随着社会的发展,化石资源储量的减少和生态环境的破坏等能源问题变得日益显著。
计算流体力学作业FLUENT 模拟燃烧问题描述:长为2m、直径为的圆筒形燃烧器结构如图1所示,燃烧筒壁上嵌有三块厚为 m,高 m的薄板,以利于甲烷与空气的混合。
燃烧火焰为湍流扩散火焰。
在燃烧器中心有一个直径为 m、长为 m、壁厚为 m的小喷嘴,甲烷以60 m/s的速度从小喷嘴注入燃烧器。
空气从喷嘴周围以 m/s的速度进入燃烧器。
总当量比大约是(甲烷含量超过空气约28%),甲烷气体在燃烧器中高速流动,并与低速流动的空气混合,基于甲烷喷嘴直径的雷诺数约为×103。
假定燃料完全燃烧并转换为:CH4+2O2→CO2+2H2O反应过程是通过化学计量系数、形成焓和控制化学反应率的相应参数来定义的。
利用FLUENT的finite-rate化学反应模型对一个圆筒形燃烧器内的甲烷和空气的混合物的流动和燃烧过程进行研究。
1、建立物理模型,选择材料属性,定义带化学组分混合与反应的湍流流动边界条件2、使用非耦合求解器求解燃烧问题3、对燃烧组分的比热分别为常量和变量的情况进行计算,并比较其结果4、利用分布云图检查反应流的计算结果5、预测热力型和快速型的NO X含量6、使用场函数计算器进行NO含量计算一、利用GAMBIT建立计算模型第1步启动GAMBIT,建立基本结构分析:圆筒燃烧器是一个轴对称的结构,可简化为二维流动,故只要建立轴对称面上的二维结构就可以了,几何结构如图2所示。
(1)建立新文件夹在F盘根目录下建立一个名为combustion的文件夹。
(2)启动GAMBIT(3)创建对称轴①创建两端点。
A(0,0,0),B(2,0,0)②将两端点连成线(4)创建小喷嘴及空气进口边界①创建C、D、E、F、G点②连接AC、CD、DE、DF、FG。
(5)创建燃烧筒壁面、隔板和出口①创建H、I、J、K、L、M、N点(y轴为,z轴为0)。
②将H、I、J、K、L、M、N向Y轴负方向复制,距离为板高度。
③连接GH、HO、OP、PI、IJ、JQ、QR、RK、KL、LS、ST、TM、MN、NB。
fluent凝固过程模拟
Fluent是一种流体动力学(CFD)软件,可以用于模拟各种凝固过程。
以下是一些常见的凝固过程模拟:
1.熔体凝固。
熔体凝固是固体制备和合金制备中的重要过程。
在Fluent中,可以
使用相变模型或传热模型来模拟熔体凝固。
相变模型可以描述相变过程中
的温度和浓度变化,而传热模型可以描述传热过程中的温度和相变。
2.熔铸凝固。
熔铸凝固是一种制备单晶或多晶的方法。
在Fluent中,可以使用多
相流模型来模拟熔铸凝固。
多相流模型可以描述不同相之间的相互作用,
包括气-液、气-固、液-固之间的传热和质量传递。
3.熔凝法制备陶瓷。
熔凝法是一种制备陶瓷和玻璃的方法。
在Fluent中,可以使用多相
流模型来模拟熔凝法。
多相流模型可以描述熔体和固体之间的相互作用,
以及气体的传输和反应。
4.铸造过程。
铸造是一种制备金属件的方法。
在Fluent中,可以使用多相流模型
来模拟铸造过程。
多相流模型可以描述金属液体和气体之间的相互作用,
包括浸润、气泡形成和气体排放等。
5.涂层制备。
涂层制备是一种将功能性材料涂覆在基材上的方法。
在Fluent中,可以使用多相流模型来模拟涂层制备过程。
多相流模型可以描述涂料和基材之间的相互作用,以及气体的传输和反应。
总之,Fluent可以用于模拟各种凝固过程,可以根据不同的物理和化学过程选择不同的模型进行模拟,并可以对流场、温度场、浓度场等进行分析和优化。
- 20 -高 新 技 术0 引言相变材料(Phase Change Material,PCM)通过吸收/释放潜热以储存能量,可实现能量在时空上的转换,缓解能源危机。
目前,将胶囊封装技术应用于相变材料领域是国内外的研究热点。
与普通相变材料相比,相变材料微胶囊化可以增加传热面积、减少与外界环境的反应。
相变微胶囊具有良好的储热调温特性,可应用于建筑节能、控温包装等领域。
相变过程包含多种复杂的热物理现象,相变微胶囊尺寸很小,而传统实验基本难以观测到微胶囊内部的相变,采用数值方法可以有效解决这一难题。
热焓法是被广泛应用的一种数值方法,它是以温度与焓值作为变量,在整个区域内建立统一的能量方程求解。
叶会文等运用热焓法研究了相变微胶囊的蓄热特性,发现相变微胶囊的换热效率高,适用于小温差的储热系统[1]。
戴晓丽建立了单个相变微胶囊凝固过程的模型,得到胶囊粒径、ste 数等对胶囊相变过程的影响规律[2]。
郝睿研究了相变微胶囊的壁材对相变过程的影响,发现微胶囊壁材的厚度和导热系数对相变过程有着很大影响[3]。
在相变微胶囊的传热过程研究中,目前较少有在微米尺度上预测相变微胶囊内的温度分布、固/液界面等瞬态现象,系统地分析了胶囊粒径、壁材的厚度、壁材的导热系数对相变过程的影响。
该文以芯材是石蜡、壁材是密胺树脂的相变微胶囊单体为研究对象,运用热焓法对相变微胶囊的熔化过程进行了数值模拟,预测了微胶囊内的物相分布情况,研究了胶囊粒径、囊壁厚度、壁材的导热系数对相变传热过程的影响,为相变微胶囊在控温包装领域的设计和应用提供了理论指导。
1 模型的建立采用Fluent 软件进行模拟仿真。
下面介绍其物理模型和计算过程。
1.1 物理模型图1为石蜡相变微胶囊的电镜扫描微观结构图,选取其中的一个单体进行研究,其示意图如图2所示,微胶囊外壁面温度恒定。
为了简化计算,对相变微胶囊模型进行5条假设。
1)相变微胶囊的囊壁、相变芯材具有均质和各向同性。
fluent熔化凝固模型液相分数1. 引言熔化凝固模型是一种在工程和科学领域中常用的模型,用于研究物质的相变过程。
在熔化凝固模型中,液相分数是一个重要的参数,用于描述物质在熔化和凝固过程中液态物质的比例。
在本文中,我们将深入探讨使用Fluent软件来模拟和分析熔化凝固过程,并重点研究液相分数对物质性质的影响。
2.熔化凝固模型熔化凝固模型是通过建立数学方程来描述物质在熔化和凝固过程中的温度变化和相变过程。
其中液相分数是一个重要的参数,可以通过质量分数或体积分数来表示。
在熔化过程中,液相分数会逐渐增加,而在凝固过程中,液相分数会逐渐减少。
3. Fluent软件介绍Fluent是一种流体动力学模拟软件,常用于模拟和分析流体和气体的行为。
它可以模拟多种物理现象,包括热传导、流动、传热和相变过程。
在熔化凝固模型中,Fluent可以用来模拟和分析物质的温度变化、相变过程以及液相分数的变化。
4.使用Fluent模拟熔化凝固过程4.1 设置模拟参数在Fluent中,首先需要设置模拟的参数。
这包括物质的热物性、初始温度、凝固温度等。
还需要定义模拟域的边界条件,如热通量、边界温度等。
4.2 定义模型在Fluent中,可以选择不同的模型来模拟熔化凝固过程。
常用的模型包括体积平衡法、相场法、多相流模型等。
这些模型可以根据物质的性质和模拟需求进行选择。
4.3 求解方程在Fluent中,需要定义一组方程来描述熔化凝固过程。
这包括质量守恒方程、能量守恒方程、相变方程等。
这些方程会随着时间的推移进行求解,得到物质的温度分布、液相分数等相关数据。
5.分析液相分数对物质性质的影响使用Fluent模拟熔化凝固过程后,我们可以通过分析液相分数对物质性质的影响来深入了解这一过程。
具体的分析方法包括:5.1 温度变化分析通过对模拟结果中的温度数据进行分析,可以了解液相分数对物质的温度变化有何影响。
例如,在熔化过程中,液相分数的增加是否会导致温度的降低?5.2 相变速率分析通过对模拟结果中的相变速率数据进行分析,可以了解液相分数对物质的相变速率有何影响。
基于相变流体的船用柴油机缸套冷却的数值模拟
颉利东;邹得球;王军建
【期刊名称】《中国水运(下半月)》
【年(卷),期】2014(014)008
【摘要】利用CFD软件对6,300船用柴油机缸套冷却进行数值模拟,建立了柴油机缸套冷却腔的三维模型,分别以水和相变石蜡乳状液作为冷却介质,获得柴油机冷却腔温度场数据.对以上两种冷却介质的数值计算比较分析,发现相变石蜡乳状液作为冷却介质时,其冷却腔温度分布较均匀,且整体温度较低,缸套冷却效果更好.文中可为船用大功率柴油机相变冷却系统的优化设计提供技术指导.
【总页数】4页(P95-97,112)
【作者】颉利东;邹得球;王军建
【作者单位】宁波大学海运学院,浙江宁波315211;宁波大学海运学院,浙江宁波315211;宁波大学海运学院,浙江宁波315211
【正文语种】中文
【中图分类】TK424
【相关文献】
1.船用低速柴油机缸套冷却水及处理 [J], 刘飞
2.船用柴油机排气冷却消声器流体与消声特性 [J], 陈涛;胡霖;黎南
3.高强化柴油机气缸套周围冷却水流动的数值模拟和试验研究 [J], 屈盛官;夏伟;王颖;黄荣华;成晓北
4.船用GY8300型柴油机缸套冷却水腔腐蚀 [J], 陈世兴
5.柴油机气缸套冷却水空化流的三维数值模拟 [J], 夏冬生;张会臣;于彦
因版权原因,仅展示原文概要,查看原文内容请购买。
石油工程应用实例7.1 热油管路温降分析7.1.1 案例简介与思路分析我国所产原油80%以上为凝点较高的含蜡原油。
目前含蜡原油多采用加热管输的方式。
在热油管道运行过程中,不可避免地会因计划性或事故性的原因而发生停输。
原油管道的停输再启动,一直是管道安全运行的核心技术问题。
由于热油管道大部分都采用埋地敷设,因此在实际生产、运行中对于停输后热油管道以及土壤温度场分布的研究,具有重要的意义和作用。
【案例简介】本例的输油管为Φ500×6mm钢管,管中心埋深为1.5m,管外保温层厚度10mm,沥青防腐层厚度6mm。
由于涉及到与土壤环境的传热问题,需要考虑大地恒温层及管道热力影响区的水平边界来划定土壤区域范围,本例中选择土壤垂深10m,水平方向的热力半径为10m,如图7-126所示。
图7-1 管路温降几何示意图【思路分析】在模拟的过程中分别划分出流体区域与固体区域,然后在选择求解器时采用非耦合式隐式算法,打开能量方程,由于存在温度梯度,因此需要编写UDF来定义边界温度。
整个模拟计算的流程如图7-127所示。
图7-2 井下节流阀流动模拟流程图【光盘文件】——参见附带光盘中的“START\Ch7\7-6.msh”文件。
——参见附带光盘中的“END\Ch7\7-6.cas,6-6.dat”文件。
——参见附带光盘中的“AVI\Ch7\7-6.avi”文件。
[1]启动GAMBIT,单击→→,依次创建坐标(0,0,0)、(0.25,0,0)、(0.256,0,0)、(0.266,0,0)、(0.272,0,0)的5个点,如图7-128所示。
图7-3 创建的点[2]单击→→,打开Move/Copy Vertices对话框,选中之前创建的5个点,选择Copy选项,在Operation选项卡中选择Rotate,在旋转角度Angle中输入90,旋转轴Axis 选项保持默认,即围绕z轴旋转,单击按钮即可。
在旋转角度Angle中输入-90,其他选项保持不变,单击按钮,得到的点如图7-129所示。
计算流体力学作业FLUENT 模拟燃烧问题描述:长为2m、直径为0.45m的圆筒形燃烧器结构如图1所示,燃烧筒壁上嵌有三块厚为0.0005 m,高0.05 m的薄板,以利于甲烷与空气的混合。
燃烧火焰为湍流扩散火焰。
在燃烧器中心有一个直径为0.01 m、长为0.01 m、壁厚为0.002 m的小喷嘴,甲烷以60 m/s的速度从小喷嘴注入燃烧器。
空气从喷嘴周围以0.5 m/s的速度进入燃烧器。
总当量比大约是0.76(甲烷含量超过空气约28%),甲烷气体在燃烧器中高速流动,并与低速流动的空气混合,基于甲烷喷嘴直径的雷诺数约为5.7×103。
假定燃料完全燃烧并转换为:CH4+2O2→CO2+2H2O反应过程是通过化学计量系数、形成焓和控制化学反应率的相应参数来定义的。
利用FLUENT的finite-rate化学反应模型对一个圆筒形燃烧器内的甲烷和空气的混合物的流动和燃烧过程进行研究。
1、建立物理模型,选择材料属性,定义带化学组分混合与反应的湍流流动边界条件2、使用非耦合求解器求解燃烧问题3、对燃烧组分的比热分别为常量和变量的情况进行计算,并比较其结果4、利用分布云图检查反应流的计算结果5、预测热力型和快速型的NO X含量6、使用场函数计算器进行NO含量计算一、利用GAMBIT建立计算模型第1步启动GAMBIT,建立基本结构分析:圆筒燃烧器是一个轴对称的结构,可简化为二维流动,故只要建立轴对称面上的二维结构就可以了,几何结构如图2所示。
(1)建立新文件夹在F盘根目录下建立一个名为combustion的文件夹。
(2)启动GAMBIT(3)创建对称轴①创建两端点。
A(0,0,0),B(2,0,0)②将两端点连成线(4)创建小喷嘴及空气进口边界①创建C、D、E、F、G点C D E F Gx 0 0.01 0.01 0 0y 0.005 0.005 0.007 0.007 0.225②连接AC、CD、DE、DF、FG。
相变材料系统流程模拟与优化摘要:本文采用PROII模拟软件,以某炼厂相变材料生产装置减压蒸馏系统为基础,分别研究了回流比、不同塔压、进料位置、进料组成变化等参数对该系统的影响,对减压蒸馏系统在进料组成变化的情况进行了适应性评价。
通过合理优化操作条件,在保证相变材料产品质量控制不变的条件下,探索了分离相变材料的合理化工艺路线。
关键词:减压蒸馏;相变材料;操作参数; PROII储热技术的开发和利用能够有效提高能源综合利用水平,主要应用于热电联产、集中供热、风电消纳以及移动供热等。
储热技术的核心和关键是相变材料。
近十年来,储能相变材料一直备受关注,人们对相变材料作为热能储存的兴趣几乎呈指数增长,是国际研究的热点。
相变材料是指在保持温度不变的情况下发生相变并提供潜热的物质,可用于热量贮存和温度控制领域,对于提高能源利用效率、改善能源结构具有重要意义。
烃类物质用作相变材料时,相变潜热与烃类物质的化学组成直接相关,其中正构烷烃含量越高、碳数分布越窄,相变潜热大,优先选用纯正构烷烃作为相变材料几乎是本领域公知的技术[1]。
快速开发出n-C12~n-C22等正构烷烃的生产技术以替代进口产品是国内企业不断开发的目标。
1 相变材料系统分离工艺石油基相变材料的分离可采用减压精馏方式,利用精确馏程切割正构烷烃产品。
减压蒸馏在炼化工业主要应用于原油减压装置深拔处理过程,石油基相变材料减压蒸馏可以参考深拔处理流程精确切割产品。
相变材料减压精馏系统主要由减压塔构成,并附加换热器、分离罐、泵等其他设备。
蜡液自液蜡原料罐经进料泵转送至减压塔T-01,通过T-01脱除轻组份。
T-01塔底物料由泵直接打入减压塔T-02,塔顶产品为NC18。
T-02塔底物料由泵打入T-03,T-03分离塔顶产品为NC19。
T-03塔底物料由泵打入T-04,T-04分离塔顶产品为NC20。
T-04塔底物料送至罐区储存完成切割任务。
2 减压精馏系统流程模拟与分析2.1原料性质本文减压精馏系统的原料来自含蜡原油经常减压装置粗分后再通过糠醛精制和酮苯脱油脱蜡得到的蜡油,原料流量为4 t/h,塔顶压力为2~4kPa。
金属翅片、骨架强化石蜡相变传热的数值模拟
常钊;陈宝明;王惠临;罗丹
【期刊名称】《煤气与热力》
【年(卷),期】2022(42)3
【摘要】在充满相变材料(石蜡)的矩形腔体内分别添加金属翅片、金属翅片和金属骨架,采用模拟方法研究纯石蜡方腔、含翅片方腔、含翅片(即金属翅片)-骨架(即金属骨架)方腔对腔体内石蜡相变传热的影响。
3种方腔内的相变传热过程均为导热与自然对流传热共同作用。
含翅片-骨架方腔中的石蜡熔化速率最快,且完全熔化时间最短。
含翅片方腔中的石蜡熔化速率最慢,完全熔化时间最长。
纯石蜡方腔中的石蜡熔化速率与完全熔化时间居中。
相同加热时间,含翅片-骨架方腔的中心线温度分布最均匀,且高于其他两种方腔。
含翅片-骨架方腔的中心线平均温度最先达到稳定,然后是纯石蜡方腔,最后是含翅片方腔。
在3种方腔中,含翅片-骨架方腔有利于加速石蜡熔化速率,缩短熔化时间。
含翅片方腔阻碍自然对流传热,不利于石蜡的相变传热。
【总页数】6页(P22-27)
【作者】常钊;陈宝明;王惠临;罗丹
【作者单位】山东建筑大学热能工程学院
【正文语种】中文
【中图分类】TU832.1
【相关文献】
1.铝泡沫和石墨泡沫强化石蜡相变传热的数值模拟
2.带翅片的小型平板CPL蒸发器相变传热的数值模拟
3.金属骨架断面对石蜡相变传热影响的模拟分析
4.不同结构金属翅片对石蜡相变传热的影响
5.含不同结构金属骨架石蜡相变传热数值模拟
因版权原因,仅展示原文概要,查看原文内容请购买。
含蜡原油在模型环道中温度场的数值模拟
李男;陈艺彤;李莹
【期刊名称】《北京石油化工学院学报》
【年(卷),期】2018(026)003
【摘要】结合室内模型环道装置进行蜡沉积实验,应用Fluent软件对环道测试段温度场进行数值模拟计算.模拟了不同壁温、油温条件下管道轴向不同位置内管壁的径向温度场,利用等效导热系数的方法对不同蜡沉积层厚度下管道内管壁的温度分布进行模拟.模拟结果表明,油温与壁温间温度差、蜡沉积层厚度均会对内管壁处温度梯度产生影响;油温与壁温间的温度差越大,内管壁处的温度梯度越大;随着沉积层厚度增加,保温效果变得明显,内管壁处的温度梯度减小.随着轴向位置距离的增加,内管壁的温度梯度值呈现减小的趋势.
【总页数】6页(P46-50,64)
【作者】李男;陈艺彤;李莹
【作者单位】中国石油大学(北京),北京 102249;中国石油大学(北京),北京 102249;中国石油大学(北京),北京 102249
【正文语种】中文
【中图分类】TE832
【相关文献】
1.已加工表面热源模型研究及磨削温度场数值模拟 [J], 王德祥;孙树峰;颜丙亮;刘新福;江京亮
2.单晶连铸凝固过程温度场数值模拟(Ⅰ)——温度场的物理模型及数学模型 [J], 范新会
3.寒区引水隧洞低温相变温度场-渗流场耦合模型及其数值模拟研究 [J], 姜海波;貊祖国;后雄斌;金瑾
4.模型环道蜡沉积装置测试段流场温度场数值模拟 [J], 孙涛;蔡均猛
5.基于CEL模型的搅拌摩擦焊温度场及材料流动数值模拟 [J], 武晓燕;罗巍;曹志明;王怡嵩;江海涛
因版权原因,仅展示原文概要,查看原文内容请购买。
Liquid-Solid Phase Change Modeling Using FluentAnirudh and Joseph Lam 2002 Fluent Users’ Group MeetingSolidification Modelu uFLUENT can be used to solve fluid flow problems involving both solidification and melting for pure materials and alloys Instead of tracking the liquid-solid front explicitly, an enthalpyporosity formulation is used For pure materials, (Tsolidus and Tliquidus are equal), a method based on specific heat is used The energy equation is written in terms of enthalpy:∂ r (ρH )+ ∇.(ρ uH ) = ∇.(k∇T ) + S ∂tuuuThe liquid fraction, can be defined as:u uβ = 0 if T ≤ Tsolidus β = 1 if T ≥ Tliquidusβ= T − Tsolidus Tliquidus − Tsolidusuif Tsolidus< T <Tliquidus ,2UGM 2002ConfidentialSolidification ModeluuuuOther relationships between the liquid fraction and temperature (and species concentrations) are possible, but are not implemented yet The latent heat content/release can be written in terms of the latent heat and state of melting of the material, L: ∆ H = β L Instantaneous latent heat content/release can vary between zero (for a solid) and L (for a liquid) Solution for temperature is essentially an iteration between the energy and the liquid fraction equationUGM 2002Confidential3Treatment of Mushy Zone in Alloysu uuIn Mushy zone, liquid fraction lies between 0 and 1 Mushy zone is modeled as a ''pseudo'' porous medium: porosity decreases from 1 to 0 as the material solidifies Appropriate momentum sink terms are added to account for the pressure drop caused by the solid material(1 − β ) 2 r r S= 3 A mush ( v − v p ) (β + ε )uβ : ε : A r mush : vP :liquid volume fraction small number to prevent division by 0 is the mushy zone constant and the solid pull velocityuPorosity in a solidified cell is zero: velocities takes on the solid (pull) velocity Similar sinks are also added to the turbulence equations (1 − β ) 2 S= 3 A φ (β + ε ) mush4UGM 2002ConfidentialContact Resistance Due to ShrinkageuWall T Tw Near wall cellulTw Rc l/k TPresence of an air gap between the walls and the solidified material causes an additional heat transfer resistance between the cooling walls and adjacent fluid with liquid fraction less than 1 This contact resistance is accounted for by adding appropriate resistance in the wall heat flux calculation:Contact Thermal Resistance Near the Wall( T − Tw ) q= ( l / k + R c ( 1 − β ))k : thermal conductivity β : liquid volume fraction, and Rc : contact resistance[m 2K / W]UGM 2002 Confidential5Pull Velocity for Continuous CastinguIn continuous casting processes, the solidified material is “continuously” pulled out from the computational domainu uuSolid material has a finite velocity that influences solidification front The exact computation of the pull velocity for the solid material is dependent on the E and σ of the solid and the forces acting on it In FLUENT, an approximated Laplacian equation is used:r ∇ vp = 02uAppropriate Dirichlet boundary conditions are used on walls and outlets and the fluid velocity is picked up at the interface mold Solidified Shell PullHorizontal Continuous Casting of IronUGM 2002 Confidential6Solidification/Melting Model InputsuDefine/models/ solidification and melting ...uDefine/materials...UGM 2002Confidential7Continuous Casting of Round BilletSlag (oil) InletInput Parameters : Parameters:u uMoldLiquid Symmetry axis Moving wallu u u u u u u uDiameter of the mold : 0.115 m mold: Diameter of the nozzle : 0.03 m nozzle: Mold length : 0.50 m length: Casting speed: 0.03 m/s speed: Solidus temperature : 1490 °C temperature: Liquidus temperature : 1530 °C temperature: Melt superheat: 20 °C superheat: Mold heat transfer coefficient : 1270W/m2/°C coefficient: Spray heat transfer coefficient: 1080W/m2/°C coefficient: Simulated length : 3.0 m length:Casting speedSchematic configurationUGM 2002Confidential8Continuous Casting of Round BilletInput Parameters (continued): (continued):Density: Density: 7020 kg/m3 u Viscosity: Viscosity: 6.2x10 -3 J/kg u Latent heat of fusion: 270x103 J/kg fusion: u Specific heat : 680 J/kg/ °C heat: u Thermal conductivity : 34 W/m/ °C conductivity: u Thermal expansion coefficient: 1x10-4/K coefficient: u Permeability coefficient: 5.0x10 -11 m 2 coefficient: u Mesh size : 227x48 cells (10,896 total) size: u Near wall treatment: enhanced wall treatment treatment:Solidified shell thickness (mm)0 0.0 Expt. of Ushijima (1962) Computed profile 0 20 30 40 50uPressure-outletuComputational mesh in the upper regionSolidification profile is obtained based on 30% liquid fraction (after Aboutalebi et al, 1995) Comparison of predicted solidified shell thickness with the experiments of Ushijima (1962) shows good agreement9UGM 2002ConfidentialDistance from the meniscus (m)0.51.01.52.02.53.03.5Melting of Pure Material: GalliumInsulated wall gCold wallt = 5 min (current) t = 8 min (current) t = 12 min (current) t = 5 min (Exp.*) t = 8 min (Exp.*) t = 12 min (Exp.*) t = 5 min (Rep. Comp.) t = 8 min (Rep. Comp.) t = 12 min (Rep. Comp.) *(Webb and Viskanta, NHT Vol. 9, pp. 539-558, 1986)Hot wall4.5 cmInsulated wall 9.0 cmMaterial: Gallium(Pure) Property Units Density kg/m3 Specific Heat j/kg-k Therm. Cond. w/m-k Viscosity kg/m-s Therm. Exp. Coeff. 1/k Melting Heat j/kg Solidus Temperature Liquidus Temperature1.2Validation of Pure Metal Melting (Gallium)Value 5904 371.5 29 0.001639 9.9999997e-05 80000 k 302 k 3021.0 0.8 0.6 0.4 0.2 0.0 0Y-CoordinateSte = 0.042, Pr = 0.021, Ra = 2.2 x 105UGM 2002 Confidential0.20.40.6 0.8 1 X-Coordinate1.21.410Solidification with VOFu uuuA liquid metal droplet spreads and solidifies on a cooled substrate The VOF (free surface) and phase change models in FLUENT are used together to simulate the process Results are useful for applications of liquid metal jetting, used in electronics manufacturing Also useful for printing, painting, spray gluing applicationsConfidentialUGM 200211Liquid Metal Solidificationt = 50 µs solidification (black)uVOF model uses single set of fluid equationsutracks liquid/gas volume fraction and location and shape of gasliquid interface tracks liquid-solid interface within the liquid phase source terms in energy equation depend on liquid fraction and available latent heatuPhase change model uses enthalpy-porosity methodu uUGM 2002Confidential12Liquid Metal Solidification ...t = 250 µsu usolidification (black)uSolid layer is colored black Spread factor, dfinal/dinitial = 3.2 as predicted by FLUENT is in good agreement with published numerical finding, 2.971 Simulation shows that FLUENT can capture combined phenomena of free surface flows and phase changet = 750 µssolidification (black)t = 1.9 mssolidification (complete)R. Holt and Albert Y. Tong, “The normal incidence impact and solidification phenomena of a liquid metal droplet onto a rigid substrate”, ASME FED-Vol. 234, 1995 IMECE, 215-224.UGM 2002 Confidential1Brendon13Vibrating Mold with SolidificationuA 3D round billet case is studied with mold oscillationsu uThe mold stroke: 0.025m and frequency: 1Hz The average free-board is 0.15mu uSolidification in the mold is analyzed Free-surface of the steel is assumed shearfree and not dynamically recalculated as with a VOF analysislThis is *not* a limitation but a simplification adopted hereuAll other process parameters are same as the 2D studyConfidentialUGM 200214Vibrating Mold with SolidificationSome instantaneous solution fields Velocity vectors PathlinesuTemperatureLiquid FractionUGM 2002Confidential15Vibrating Mold with Solidificationu uLiquid FractionTemperatureuThese animations on a diameter-plane The mold water region contours are showing velocity magnitude in both cases The dynamics of the flow-fields are slower than the mold-oscillation rateuuIn a couple of mold-cycles, no significant change in the flow-field should be expected Oscillation frequency in the actual mold can be slowerUGM 2002Confidential16Mold Flow with Meniscus Atmospheric Pressure and SolidificationuSEN (V-in)A slab caster moldu u uuuFlow and heat transfer Solidification Mold top is open to atmosphere Free surface at the meniscus No chemistry is consideredSEN-jet (internal) Mold wall GeometryOpen slab Wall Pressure OutletUGM 2002Confidential17Mold Flow with Meniscus and SolidificationuMold top is open to atmosphereu uuNo mold powder is considered Air on top of the mold maintains a mild recirculating flow Outlet pressure reflects the force of the support-rolls which would determine the free-board in the molduSolidified shell on mold wall develops depending on the process parametersu u u uCasting speed Incoming superheat of liquid steel The jet penetration and SEN design Mold water flow rateConfidentialVelocity Vectors Liquid FractionUGM 200218Mold Flow with Meniscus and SolidificationuuuContact resistance is automatically introduced based on the local state of solidification Similarly, pull velocity is assigned on the solidified material only Figure shows shell thickness, free surface (meniscus) and the velocity field in a slab caster mold on the wide, symmetry faceUGM 2002Confidential19Summary of Melting/Solidification ModeluFeaturesu u u u u u uEnthalpy-Porosity formulation Accessible through GUI interface Contact resistance at walls (customizable) Marangoni convection at "walls" Solid pull velocities (both user-specified and computed) VOF-compatible Sink terms for momentum & turbulent quantities in mush & solid zones Lack of species transports model User-specified H-T curveuLimitationu uuBoth species model and custom-H-T curve should be available in the next release - macro-segregation can be modeledUGM 2002Confidential20。
利用FLUENT软件模拟计算含蜡原油管道的停输降温过程张园园
【期刊名称】《管道技术与设备》
【年(卷),期】2009(000)006
【摘要】研究含蜡原油管道的停输降温过程,对确定安全停输时间、提出再启动方案以及制订停输检修计划具有重要的指导作用.介绍了利用FLUENT软件模拟求解含蜡原油管道停输降温问题的方法,将双边壁面设为耦合壁面,不需要再单独设定边界条件;将析蜡潜热转化为附加原油比热容,不需要跟踪固液相界面的移动.通过与实验结果进行对比,表明该方法能够准确地计算出任一时刻管道横截面上各点的温度,而且能表现出降温过程中自然对流的变化和固液相界面的移动,使得此类问题的求解更加简单.
【总页数】4页(P1-4)
【作者】张园园
【作者单位】中国石油化工股份有限公司石油勘探开发研究院,北京,100083【正文语种】中文
【中图分类】TE832
【相关文献】
1.热油管道停输降温过程研究进展 [J], 刘刚;张园园
2.海底含蜡原油管道停输降温过程的数值模拟 [J], 高艳波;马贵阳;曹先慧;齐浩;杜义朋
3.含蜡原油管道停输降温的数值计算 [J], 李传宪;纪冰;魏国庆
4.基于有限元法的含蜡原油管道停输过程温降数值模拟 [J], 李玉春;李中义;张弼;王姝萌
5.热油管道停输与再启动过程模拟计算软件 [J], 安家荣;史秀敏;张国忠
因版权原因,仅展示原文概要,查看原文内容请购买。
fluent相变沸腾内温度过高
Fluent是一款流体动力学模拟软件,常用于模拟流体流动、传热、传质等过程。
相变沸腾是指物质从液态变为气态的过程,通常伴随着热量的传递和相变潜热的释放。
在Fluent中模拟相变沸腾时,如果内温度过高,可能有以下原因:
1.热通量过大:这可能是由于流体进入管道时流速过高,或者管道内部存在太多的热源,导致壁面温度过高。
2.热边界条件设置不当:在模拟过程中,如果给定的热边界条件不合适,也可能导致温度过高。
3.物理模型选择不当:例如,如果选择了不合适的相变模型或传热模型,也可能影响模拟结果。
为了解决这个问题,可以尝试以下方法:
1.检查并调整热通量:根据实际情况调整流体的流速或减少管道内部的热源,以降低壁面温度。
2.重新设置热边界条件:根据实际情况重新设置合适的热边界条件,以更准确地模拟实际情况。
3.尝试不同的物理模型:根据实际情况选择更合适的相变模型或传热模型,以获得更准确的模拟结果。
总之,对于Fluent中相变沸腾内温度过高的问题,需要仔细检查模拟设置和物理模型,并根据实际情况进行调整。
fluent多相流模拟温度变化使用Fluent多相流模拟温度变化引言:在工程和科学领域中,了解物质在不同温度下的行为是非常重要的。
为了研究和预测温度变化对不同物体的影响,科学家和工程师使用了各种模拟方法。
其中一种常用的方法是使用Fluent多相流模拟软件。
本文将探讨使用Fluent多相流模拟温度变化的原理和应用。
一、Fluent多相流模拟的基本原理Fluent是一种基于计算流体力学(CFD)的软件,它使用数值方法来解决流体流动和传热问题。
多相流模拟是Fluent的一个重要功能,它可以模拟多种物质在不同温度下的相互作用和传热过程。
在Fluent中,多相流模拟是通过将模拟区域分为离散的网格单元来实现的。
每个网格单元内的物质被假设为均匀的,并且可以具有不同的热传导系数和热容量。
通过在不同网格单元之间建立质量、能量和动量的平衡方程,可以模拟物质在不同温度下的传热和流动行为。
二、Fluent多相流模拟温度变化的应用1. 工业领域Fluent多相流模拟在工业领域中具有广泛的应用。
例如,在石油和化工行业,科学家和工程师可以使用Fluent来模拟化工过程中的温度变化。
通过分析和优化温度分布,可以提高生产效率和产品质量。
2. 环境工程在环境工程中,Fluent多相流模拟也发挥着重要作用。
例如,在城市规划中,可以使用Fluent来模拟建筑物和街道上的温度分布。
通过分析不同材料的热传导性能和建筑物的排热能力,可以减少城市热岛效应,改善城市的舒适度。
3. 生物医学在生物医学研究中,Fluent多相流模拟可以用于模拟人体内部的温度变化。
例如,在热疗领域,科学家可以使用Fluent来模拟热疗过程中的温度分布,以优化治疗效果并减少副作用。
此外,Fluent还可以用于模拟血液流动和热传导,帮助医生更好地理解和治疗心脑血管疾病。
三、Fluent多相流模拟温度变化的优势1. 准确性Fluent多相流模拟使用基于数值方法的计算模型,可以准确地预测物质在不同温度下的行为。
一、 问题背景: 为了解决日益严重的能源短缺问题,如何更加充分地利用现有的 化石能源,开发利用绿色能源成为世界各国关注的重要议题。随着现 有化石能源的逐步开采利用,世界各国已经普遍认识到降低对传统能 源(如煤炭、石油、天然气能源等)的依赖性,以及对绿色能源(如 太阳能、风能、地热能等)实现充分开发利用的重要性,使用再生类 能源并通过提高能源利用效率的方式成为应对能源枯竭现状的重要 手段。蓄热技术就是这类能够提高能源利用效率的典型技术手段, 蓄 热技术通过将间歇性或者不稳定的热量通过蓄能介质暂时储存, 在有 使用要求时释放能量,解决能源利用高峰阶段造成的能源匹配不足的 问题。经过多年应用发展,蓄热技术已经在太阳能、地热能、风能、 工业废热、电网系统的“移峰填谷”等领域有了一定程度的应用,并 表现出强劲的发展势头。
二、 蓄热技术 蓄热技术一般通过利用蓄热介质的比热容、潜热等物理特性实现 对采集能源多余热量的暂时储存,主要分为显热蓄热、潜热蓄热和化 学热反应蓄热,其中前两种技术属于物理蓄热范畴。 显热蓄热 显热蓄热通过提升蓄热介质材料的温度进而提高物质内能的方 式实现储热,储热能力取决于材料的比热容(提升物质单位温升所需 要的外部能量)等物理参数。显热蓄热的突出弱点在于由于蓄能介质
(例如水)在多余能量的储存过程中伴随着物质温度的上升, 不能满 足部分 设备对于恒温放热的要求。同时,显热蓄热材料蓄热能力有限, 储能密度较低,往往需要较大容积的容器提供储能保证, 限制了显热 蓄热技术的大规模推广应用。 潜热蓄热 潜热蓄热技术充分利用了相变介质在相态改变时会吸收或释放 巨大的能量,并藉此实现对多余能量的储存和释放。 这类利用相变过 程实现能量吸收释放的材料被称为相变材料。 相比显热蓄热材料,相 变材料的相变潜热与蓄热材料的比热特性相比在储热能力方面有了 极大的提升,同等质量的储热介质能够实现对更多能量的储存, 降低 了对储热设备容积、质量等方面的要求,降低了整体设备成本。同时, 潜热蓄热材料在相变吸热放热过程中近似等温过程, 方便了实际工程 控制。尽管如此,由于相变蓄热材质在导热、传热特性方面的不足, 以及某些无机盐相变材料存在的相分离和过冷现象以及某些固有的 化学反应,会严重影响储热设备的使用寿命。 化学热反应蓄热 化学热反应在反应过程中会释放大量热量,与外部环境进行热交 换,化学热反应蓄热即是利用这一热交互过程实现能量和吸收和释 放。基本原理如下:利用正化学反应吸热,将能量转化为化学能储存, 之后利用负反应放热,将通过正反应储存的能量通过热量的方式释放 出来。化学热反应蓄热与相变潜热蓄热技术有相似的优势特点, 但是 由于化学热反应存在反应过程复杂、反应速度过快、反应过程复杂难 于控制等问题,对于热反应发生装置要求较高,目前仅仅在一些特殊 要求应用场合得到了一定程度的应用。 综上所述,在热能综合利用和蓄能技术对比中, 采用相变材质的 潜热蓄热技术具有更为突出的应用前景,如图1所示为一种典型的通 过相变材料实现蓄热的蓄热系统构成(相变蓄能热水器系统) 。该系 统首先通过吸收外部太阳辐射热量对水进行加热,并导入热水箱结构 进行储存,对于多余的能量,系统通过热水将热能再次转化为相变材 料的相变能量,实现对多余能量的储存。在夜间缺少太阳能的时候, 可以通过相变蓄热模块对冷水进行初步加热,并将初步加热的温水导 入热水箱通过传统的辅助加热设备加热到固定的使用温度。 上述系统 能够在连续多天没有阳光或者日照不足的情况下, 实现对所储存能量 的更高效的利用,降低对传统能源的依赖程度。
图1相变材料实现蓄热的蓄热系统构成 三、相变材料介绍
根据蓄热材料的组成分为无机类、有机类和复合相变材料三个大 的类别。通过相变材料的相变过程温度变化范围分为高温相变材料 (相变温度》250C),中低温相变材料(相变温度:100C〜250C), 低温相变材料(相变温度W 100C)。如图2所示,对相关相变材料分 类及类别内的物质种类进行了简要划分。
无机相变蓄热材料:熔融盐类、结晶水合盐类、金属合金类无机物
混合相变蓄热材料:通过将有机相变材料和无机相变材 料进行一定比例的混合。
潜热蓄热 材料
低温蓄热材料:石蜡、脂肪酸、醇类、盐类水合物;
高低温蓄热材料:单纯盐类、金属与合金、金属氧化物、混合盐类; 图2相变材料划分及典型代表
本次毕业设计采用ANS YS Flue nt计算流体力学仿真软件对考虑 对流状态下的石蜡相变材料融化和凝固过程进行了数值模拟, 研究石 蜡材料在热交互过程中的储热机理和传热效果, 为基于潜热蓄热激励 的储能设备的设计提供参考数据。
4.1固-液相变传热相关机理 相变材料的凝固和融化过程实际上完全是一个物理变化过程, 在 热力学中往往利用潜热的概念(KJ/kg )对单位质量的相变物质相态 转化所需要的能量进行定义,并称相变材料在相态转化过程中的潜热 为溶解热(或者凝固热)。相变材料相态转变过程主要分为三个阶段: 1. 相变材料的显热存储阶段:相变材料的温度低于外部热源温 度,材料与热源发生热传导实现材料之间的能量交互, 相变材料温度 逐步上升并达到物质熔点温度。
2. 相变材料的潜热存储阶段: 相变材料达到材料熔点后, 依旧持 续从外部热源吸收热量并转化为自身能量, 但在存储阶段相变材料的 温度不变, 吸
依据材料 类别 有机相变蓄热材料:石蜡、脂肪酸类、多元醇等有机物
依据材料 相变温度 中低温蓄热材料;
四、基于 Flue nt的石蜡相变材料蓄热模拟 收的热量转化为融化热能量, 并促使相变材料由固态转 变为液态。 3. 相变材料的第二显热存储阶段: 当全部相变材料转换为液态形 式,持续从外部热源吸收的热量将促使液态相变材料温度进一步升 高。 相变传热过程是一个复杂的固 -液两相转化过程,相变过程中固 - 液交接面的移动依赖于外部边界条件, 界面移动的速度取决于相变材 料结构、属性以及外部热力学边界条件,涉及了固 - 液两相界面的移 动,在数学上属于高度非线性的复杂问题, 无法单独对系统方程进行 求解,需要采用单独的求解器对系统进行分别求解。 4.2 相变传热的数学模型和求解机理 由于相变材料的相变过程中, 对于单一的相变材料, 由于有确定 的相变温度, 存在不同状态条件下的相变材料共存的现象。 不同状态 的相变材料由于具有不同的物理属性(导热性质和比热容) ,并且在 固态相变材料中仅仅具有热传导方式的热量转化而在流体状态的相 变材料中还包括了热对流的传热方式。目前,对于固 - 液相变材料传 热模拟的数值求解思路主要包括以下两种: 1.温度法模型 温度法模型适合于单一成分有固定相变温度的相变材料传热过 程模拟,以温度为系统因变量,分别对不同状态下的相变材料进行控 制能量方程的建立,并进一步进行温度求解。 2.焓法模型 相比前述温度法模型,焓法模型不仅适用于单一组分的相变传热 模拟过程,对于模糊的边界移动等传热问题同样适用,该方法将焓作 为系统因变量,在固、液以及混合模糊区域建立统一的能量守恒控制 万程。
图3固-液两相系统示意图 如图3所示,为典型的相变材料固-液两相系统示意图,选择对 应控制体
V作为研究对象,经过相界面的移动,在时刻 t,系统分为 固体相"s和液体相人,将相间由于温度变化产生的循环流动速度定义 为v,控制体内材料的变化以及相界面移动取决于交界面 S处的传热边 界条件,初始材料温度以及材料的热力学物理属性。 对于,焓法模型, 应用积分形式(如公式1):
式中,以温度和焓参数作为系统因变量,相互关系可以通过公式phvdA =f kS7TdA + jqdV
(1)
相界 2实现转换: fh~\
Cl
公式2中,参数入和勺分别表示相变材料固相比热和液相比热, 参数仪和也分别表示相变材料的液体相与固体相焓值,单位为 J/kg。 将上述方程带入分别表示固体相和液体相的能量控制方程(公式3和 公式4): 叭
dTl 刊味丽+ -也)刊•(勺叫)+切⑷
式中,’「和"分别表示固体相和液体相的密度参数,[和"表示 对应相态的导热系数,°和八分别表示对应相态的温度,灯表示梯度 算子。想要求解温度控制方程,需要首先对流体场的速度方程及相关 系数进行求解。将公式1,2带入公式3,4后,可以得到如下方程:
p- = fcv'r(5)
\p. {kh"二匚:⑹ P _
5 通过上述公式可以将活动区域以及相态界面区域的能量控制方 程统一为一个共同的控制方程组,实现了相变传热模拟的数值计算。
对于上述问题的求解方法,对于一维问题求解主要包括:Neuma n法, Lightfoot 积分法, Paterson 法,对于多维问题求解主要包括有限差 分法,有限元法和有限体积法。 4.3 基于 Fluent 的石蜡材料相变传热模拟 4.3.1 Fluent 计算流体力学数值模拟软件介绍 Fluent 是一款较为常用的集成流体、热力学领域数值模拟技术 的数值模拟软件,采用了多重网格技术,具有求解收敛稳定、收敛速 度快的特点,软件提供了丰富的物理模型。依托ANSYS有限元仿真软 件多场耦合平台,ANSYS Mes模块为Flue nt提供了较好的前处理功 能,可以高效集成多种CAD辅助建模软件,能够生成二维数值模拟所 需要的三角形、四边形网格,生成三维数值模拟所需要的四面体、六 面体网格,并能够方便实现对模型细节网格的处理,包括网格细化、 节点耦合等处理方式。同时,ANS YS平台为Flue nt软件提供了更人 性化的结果后处理窗口 CFD-Post。ANS YSF台下基于Flue nt的流体 数值模拟一般流程如图 4 所示。