1流动、传热及传质的控制方程
第五章 多组分反应流体守恒方程
2. 混合物分数
假设流量为 1 公斤/秒的混合物由两种成分混合而成,燃料的流量为 f 公斤/秒,空气的
流量为(1-f)公斤/秒。定义混合物分数 f 为燃料中所含元素的质量除以混合物的质量。于
是混合物分数 f 可以用流动中任一点的燃料、氧化剂和燃烧产物的质量分数来表示:
f
= YF
+
ν
1 +
1 YPr
下面将证明 f 是守恒量。
=
m&
''
⎡ ⎢(h
+
⎣
u2 2
) x+Δx
− (h +
u2 2
⎤ )x ⎥
⎦
上式除以 Δx ,并取极限 Δx → 0 ,得
− dQ& x'' = m& '' ( dh + u du )
dx
dx dx
(5-22) (5-23)
热流通量包括热传导产生的热通量和由于组分扩散引起的附加焓通量,如果不考虑热 辐射,热流通量的一般矢量表达形式如下:
t?xuxxpxuxxuuu???????????????520四能量守恒方程四能量守恒方程根据热力学第一定律控制体内能量变化率等于获得的外热的总和与对外做功的总和
第五章 多组分反应流体守恒方程
燃烧现象包含流体流动、传热、传质和化学反应以及它们之间的相互作用。燃烧过程是 一种综合的物理化学过程。本章我们将介绍控制燃烧过程的基本方程组:混合物质量守恒方 程、组分质量守恒方程、动量守恒方程以及能量守恒方程。着重介绍以后各章需要用到的多 组分反应流体一维和二维流动的守恒方程组,以便为分析各类火焰现象奠定基础。在附录 G 比较详细地推导了多组分反应流体多维流动的守恒方程和 Shvab-Zeldovich 变换公式,供读 者参考。
1 流体流动与传热
流体静压力的计算
如图所示,采用一复式U形压差计测量容器中O点处的压力,两段U形 管A和B中水银柱读数分别为60、70cm,中间一段充满水,求O点处 的压力
pO p2 0.8 w g
200
O
pa
RB=700
' p2 p2 p1' 0.6 Hg g
1 1’ 1’’
RA=600
P+dp
dp
dz
gdz 0
p
z
对于不可压缩流体,密度为常数
p
gz 常数
0
1 2
0
p1
gz1
p2
gz2
h z1 z2
p2 p1 g z1 z2
或
p2 p1 gh
注:上式只适用于重力场中静止的不可压缩的单一连续流体; 静压力只与各点的垂直位置有关,而与水平位置无关; 只有在压力变化不大时,气体才可适用上式。
聚合反应器与聚合反应操作
课程的教学内容
化工原理(流体流动与传热) 化学反应工程基础
聚合反应工程分析
搅拌聚合釜内的传递过程 搅拌聚合釜的放大 聚合过程与聚合反应器
第1章
流体流动与传热
1.1 流体流动
泵
水
水 池
水 封
连续介质假定:将流体视为由无数质点组成的、彼 此间没有空隙的连续介质。
R
d
2 1 u1
D
2
u2
1
管道两测点间连接压差计读数代表什么意义?
B
p1 pA gh1 p2
n
lg m xi lg i
i 1
基于Fick_扩散模型的陶瓷坯体干燥数值模拟
第42卷第8期2023年8月硅㊀酸㊀盐㊀通㊀报BULLETIN OF THE CHINESE CERAMIC SOCIETY Vol.42㊀No.8August,2023基于Fick 扩散模型的陶瓷坯体干燥数值模拟霍㊀平1,2,刘俊帅1,2,王虹凯1,2(1.华北理工大学机械工程学院,唐山㊀063210;2.河北省工业机器人产业技术研究院,唐山㊀063210)摘要:本文使用Fick 扩散模型对陶瓷坯体的干燥过程进行了模拟,在实际干燥工艺参数数值范围内,分析了热风速度㊁温度和相对湿度三个因素对陶瓷坯体内部温度㊁含水率和干燥速率的影响㊂结果表明,增加热风速度能够加大坯体表面区域的水分散失速率,但对坯体内部含水率变化影响较小;提高温度能够显著增加坯体内部的干燥速率,当温度从35ħ增加至75ħ时,最大干燥速率的变化幅度为46.34%;相对湿度对坯体平衡含水率影响较大,当相对湿度从5%增大至85%时,平衡含水率从0.8%增大至5.1%(均为质量分数),提高相对湿度能够改善坯体干燥均匀性,保证坯体干燥质量㊂模拟和试验数据基本吻合,计算结果将为进一步深入研究陶瓷坯体干燥的传热传质过程,以及后续干燥曲线的优化提供理论依据㊂关键词:陶瓷坯体;传热传质;干燥参数;温度分布;含水率分布;干燥速率中图分类号:TQ174㊀㊀文献标志码:A ㊀㊀文章编号:1001-1625(2023)08-2904-11Numerical Simulation of Ceramic Green Body Drying Process Based on Fick Diffusion ModelHUO Ping 1,2,LIU Junshuai 1,2,WANG Hongkai 1,2(1.College of Mechanical Engineering,North China University of Science and Technology,Tangshan 063210,China;2.Hebei Industrial Robot Industrial Technology Research Institute,Tangshan 063210,China)Abstract :The numerical simulation of drying process in ceramic green body was carried out with Fick diffusion model.The influences of hot wind speed,temperature and relative humidity on internal temperature,water content and drying rate of ceramic green body were analyzed within the numerical range of actual drying process parameters.The results show that water loss rate of surface area of green body increases with the increase of hot wind speed,and the wind speed has little effect on water content inside green body.The drying rate inside green body is significantly increases by increasing the temperature,and the maximum drying rate changes by 46.34%when temperature increases from 35ħto 75ħ.When relative humidity increases from 5%to 85%,the equilibrium water content increases from 0.8%to 5.1%(all are mass fraction),and the increase of relative humidity can improve drying uniformity of green body and ensure the drying quality of green body.The numerical simulation results are consistent with the experimental results.The numerical simulation data provide a theoretical foundation for further research on the heat and mass transfer process in the drying process of ceramic green body and the optimization of the drying curve of ceramic green body.Key words :ceramic green body;heat and mass transfer;drying parameter;temperature distribution;moisture content distribution;drying rate㊀收稿日期:2023-04-11;修订日期:2023-06-13作者简介:霍㊀平(1968 ),女,教授㊂主要从事机器人技术方面的研究㊂E-mail:huop@0㊀引㊀言陶瓷坯体干燥是陶瓷生产工艺中的重要环节,现有的干燥曲线是经过生产实践拟合的曲线,具有普遍性,但陶瓷新产品开发过程中具有个性化差异,用实验的方法优化干燥曲线时其实验周期长,可重复性较差㊂因此,采用数值模拟方法对陶瓷坯体传热传质的干燥过程进行研究,对于精准预测干燥参数㊁缩短产品生产周期具有非常重要的意义[1-5]㊂第8期霍㊀平等:基于Fick 扩散模型的陶瓷坯体干燥数值模拟2905㊀国内外学者主要从陶瓷坯体内部参数对坯体热风干燥的影响进行研究㊂Nait-Ali 等[6]研究分析了干燥过程中高岭土坯体的导热系数与含水量的关系,结果表明当含水率较高时,导热系数随含水率降低而不断增大,当含水率下降到一定值后,导热系数随着含水率下降而减小㊂李杰等[7]面临移动坯体气固耦合传热的问题时,利用FLUENT 软件模拟了辊道窑排烟段含有坯体在内的流场和温度场,模拟结果得出:可通过适当减小排烟口与坯体间的距离,以及合理调节烟道和阀门开度等来实现节能㊂Lauro 等[8]采用向后欧拉法对氧化铝浆料的干燥过程进行有限元模拟,采用与含水量有关的水分活度㊁有效水分扩散率㊁导热率和比热容等干燥参数,研究干燥初期收缩率对干燥速率的影响,模拟得到陶瓷坯体内部的平均含水率㊁干燥速率和表面平均温度随时间的变化规律㊂赖日东[9]采用欧拉薄壁模型对进风口干空气和湿空气两种工况下陶瓷坯体的干燥过程进行了模拟,得出干空气与湿空气干燥效果不同的结论㊂朱庆霞等[10]以Whitaker 体积平均方程为基础,推导出多孔介质内部热质传递的等效耦合扩散模型,探讨了干燥介质参数对干燥过程的影响,结果表明,风速对干燥速率的影响很小,介质的温度和湿度对干燥速率影响较大㊂综上所述,陶瓷坯体干燥的传热传质过程本质是坯体内部水分逐渐向表面扩散的过程,所以分析陶瓷坯体内部含水率㊁干燥速率的影响因素对陶瓷产品干燥过程具有重要意义㊂本文基于Fick 扩散模型,使用COMSOL Multiphysics 6.0对陶瓷坯体干燥过程进行数值模拟,分析干燥参数(速度㊁温度㊁相对湿度)对陶瓷坯体内部温度㊁含水率和干燥速率的影响,并探究其规律,为精准预测干燥介质参数㊁降低试验成本,以及优化后续陶瓷坯体干燥过程提供理论依据㊂1㊀物理模型及网格划分1.1㊀物理模型本文以卫生陶瓷坯体方板为例,参考实际干燥箱尺寸,建立流体域和坯体几何模型,其尺寸为400mm ˑ400mm ˑ500mm,陶瓷坯体方板尺寸为150mm ˑ150mm ˑ20mm,陶瓷坯体方板位于干燥箱的中心位置,干燥箱流体域及坯体几何模型如图1(a)所示,坯体中心截面如图1(b)所示(v 为风速)㊂图1㊀干燥箱流体域及坯体几何模型和坯体中心截面Fig.1㊀Fluid domain and green body geometry model of drying oven and central section of thebody 图2㊀网格划分情况Fig.2㊀Grid division conditions2906㊀陶㊀瓷硅酸盐通报㊀㊀㊀㊀㊀㊀第42卷1.2㊀网格划分为节约计算资源,提升网格质量,本文固体区域使用规则六面体网格,流体域为自由四面体网格㊂为提高计算精度,对固体域与流体域的分界处进行局部加密处理,网格划分情况如图2所示㊂2㊀数学模型2.1㊀模型假设考虑到模型的复杂性,对模型作如下假设:1)在干燥初始阶段,坯体内部的温度及含水率分布均匀;2)坯体内部各向同性,均匀且连续;3)坯体内部热量只通过热传导的方式传递,其内部水分只通过水分扩散的方式传递;4)壁面是绝热壁面㊂2.2㊀流体动力学控制方程2.2.1㊀控制方程干燥箱内的流体流动采用流体控制方程表示㊂1)连续性方程根据质量守恒定律,可以将连续性方程表达为[11-12](ρa)τ+Δ㊃(ρa u)=0(1)式中:ρa为湿空气的密度,kg/m3;τ为时间,s;u为空气的表观速度,m/s;Δ为向量微分算符㊂2)动量守恒方程动量守恒方程又称为Navier-Stokes方程[13],如式(2)所示㊂(ρa u)∂τ+Δ(ρa uu)+Δp=Δ㊃μeff(Δu+(Δu)T)(2)式中:Δp为压力梯度,Pa;μeff为有效流体黏度,Pa㊃s;T为向量矩阵的转置㊂3)能量守恒方程能量守恒定律即热力学第一定律[14],如式(3)所示㊂Tτ(ρa C a)+Δ㊃(ρa uC a T)=Δ(λaΔT)(3)式中:C a为气体比热容,kJ/(kg㊃K);λa为热空气导热系数,W/(m㊃K)㊂2.2.2㊀传热传质控制方程1)传热控制方程本文假设坯体的水分只在表面蒸发,因此其内部导热方程[15]如式(4)所示㊂ρC p T τ=λΔ2T(4)式中:ρ为坯体的密度,kg/m3;C p为坯体的比热容,kJ/(kg㊃K);λ为坯体导热系数,W/(m㊃K);T为坯体内部的温度,ħ㊂2)传质控制方程坯体内的水分传递遵循Fick第二定律[16-17],如式(5)所示㊂Mτ=DΔ2M(5)式中:M为物料内部的含水率,kg/kg;D为固体中有效水分扩散系数,m2/s㊂2.3㊀初始条件及边界条件试验测得坯体初始含水率为21.7%(含水率均为质量分数),初始温度为25ħ,进风口设置为速度入口㊂COMSOL Multiphysics6.0计算环境下边界条件的设置如表1所示㊂第8期霍㊀平等:基于Fick 扩散模型的陶瓷坯体干燥数值模拟2907㊀表1㊀边界条件的设置Table 1㊀Setting of boundary conditionBoundary condition Values Inlet speed /(m㊃s -1)0.2,0.4,0.6,0.8,1.0Inlet temperature /ħ35,45,55,65,75Inlet relative humidity /%5,25,45,65,85Export borders Pressure outlet Equivalent diameter /m 0.02Wall boundaries Thermal insulation and no slippage Fluid domain materials Wet air Thermal conductivity of green body /(W㊃m -1㊃K -1)1.675Specific heat capacity /(J㊃kg -1㊃K -1)(4183M t +1200)/(1+M t )Gas thermal conductivity /(W㊃m -1㊃K -1)0.0261Green body density /(kg㊃m -3)1985.82.4㊀主要评价指标2.4.1㊀干基含水率干基含水率是坯体内部的水分质量与绝干坯体质量的比值,其计算公式如式(6)所示[18-19]㊂M t =m t -m d m d (6)式中:M t 为t 时刻坯体的干基含水率;m t 为t 时刻坯体的质量,kg;m d 为绝干坯体的质量,kg㊂2.4.2㊀干燥速率干燥速率为单位时间坯体内所散失的水分质量,其计算公式如式(7)所示[20-23]㊂D r =M i +1-M i t i +1-t i (7)式中:M i 为t i 时刻物料的干基含水率,%;M i +1为t i +1时刻物料的干基含水率,%㊂2.5㊀模拟验证2.5.1㊀网格无关性验证图3㊀干燥1h 时平均含水率和平均温度随不同网格数量的变化情况Fig.3㊀Change conditions of average water content and average temperature with number of different grids drying for 1h 选用了八种不同数量的网格在相同条件下进行网格无关性验证㊂参考唐山某陶瓷厂的干燥介质曲线,选取热风速度为0.4m /s,热风温度为55ħ,热风相对湿度为45%,干燥1h 时平均含水率和平均温度随不同网格数量的变化情况如图3所示㊂由图3可知,当网格数目超过17万时,模拟结果开始保持稳定,在保证计算精度的前提下,考虑到计算资源问题,选择网格数量为21万㊂2.5.2㊀模型可靠性验证采用唐山某陶瓷厂的坯料制备试验用干燥陶瓷坯体样品,由于坯体结构简单,呈薄板状,可通过挤压成型的方式进行制备㊂将制备好的坯体样品放入干燥箱中,在风速为0.4m /s㊁温度为55ħ㊁相对湿度为45%的条件下,对坯体进行干燥,在坯体迎风端中心截面处的红线上等分取17个点㊂当干燥1h 时,用温湿度检测仪(三量WM810,温度量程为-10~60ħ,水分量程为2%~70%,测量误差范围为ʃ1%)测量该17个点处的温度和含水率数据,坯体中心截面上17个点的位置如图4(a)所示㊂在基本工况下,对坯体干燥进行模拟与实验,得到坯体内部1h 时各点的温度和含水率,模拟与实验结果进行了对比验证,温度和含水率的最大误差分别在ʃ8%㊁ʃ10%㊂同一干燥条件下,干燥1h 时的温度云2908㊀陶㊀瓷硅酸盐通报㊀㊀㊀㊀㊀㊀第42卷图和水分云图如图4(b)和4(c)所示,模拟值与实验值对比如图5所示㊂图4㊀干燥1h 时坯体中心截面的实物图与其温度分布和水分分布Fig.4㊀Physical diagrams of center section of green body and its temperature distribution and water distribution drying for 1h 图5㊀干燥1h 时坯体各点温度和含水率的模拟值与实验值对比Fig.5㊀Comparison between simulated and experimental values of temperature and water content at various points of green body drying for 1h 3㊀结果与讨论3.1㊀风速对陶瓷坯体干燥过程的影响不同风速对坯体内温度分布和含水率分布的影响有一定的规律性,为简化图形,本文中只分析了风速为0.2㊁0.6和1.0m /s 下的温度分布和含水率分布情况㊂3.1.1㊀风速对坯体内部温度和含水率影响图6为干燥1h 时不同风速下坯体表面和内部温度的分布㊂从图6中的坯体表面㊁x-y 截面和x-z 截面的温度云图分布可以看出,坯体迎风面温度较高,坯体背风面温度较低,且坯体表面温度高于内部温度,形成温度梯度,特别是风速在0.2m /s 时坯体内温度梯度更大;随着风速增加,坯体内温度逐渐升高,风速由0.2m /s 提高到0.6m /s 的坯体升温速率大于0.6m /s 提高到1m /s 时的坯体升温速率㊂垂直于流速方向的y-z 截面上的温度对称分布,且坯体表面温度高于其内部温度㊂坯体迎风面热空气湍动程度较高,对流传热速率较大,且风速能带走坯体迎风面的一部分水分,表面水分蒸发所吸收的热量相对较少,致使坯体迎风面温度较高㊂图7为干燥1h 时不同风速下坯体表面和内部含水率分布㊂从图7中的表面云图可以看出,不同风速下坯体表面含水率较低,且表面含水率无变化,其含水率均低于4%;从坯体各个中心截面可以看出,坯体外第8期霍㊀平等:基于Fick 扩散模型的陶瓷坯体干燥数值模拟2909㊀表面含水率较低,其内部含水率较高,形成明显的湿度梯度,这是因为坯体外表面直接与热空气接触,坯体表面水分受热蒸发,在传热㊁传质的推动力下,坯体内部水分向其表面迁移,并在表面蒸发㊂当风速为0.2m /s 时,坯体中心截面内部的干基含水率明显高于风速为0.6m /s 时的含水率,这是因为热风速度较低时,坯体表面气流湍动能程度较弱,坯体表面水蒸气浓度较大,不利于表面水分蒸发,坯体内部与表面间的含水率梯度小,传质驱动力小,干燥速率较低,导致坯体内部含水率较高㊂图6㊀干燥1h 时不同风速下坯体表面和内部温度的分布Fig.6㊀Temperature distribution of surface and interior of green body at different wind speeds drying for 1h 图7㊀干燥1h 时不同风速下坯体表面和内部含水率分布Fig.7㊀Water content distribution of surface and interior of green body at different speeds drying for 1h 图8为热风温度㊁相对湿度恒定时风速对含水率的影响㊂由图8可知,当坯体干燥1~4h 时,随着风速增加,含水率变化较小,当风速从0.2m /s 增加到0.4m /s 时,在2h 时干基含水率从10.8%减小至10.1%,仅减小了6.48%,这是因为风速的增加增强了坯体表面气体流动的湍动程度,削弱了边界层效应,坯体与空2910㊀陶㊀瓷硅酸盐通报㊀㊀㊀㊀㊀㊀第42卷气域界面处的对流传质速率增大,使坯体表面的区域水分浓度差变大,临近坯体表面的内部区域水分扩散速率增大㊂但是坯体表面对流传质速率远大于坯体内部扩散速率,因此热风速度对坯体内部含水率的影响较小㊂3.1.2㊀风速对坯体内部干燥速率影响图9为热风温度㊁相对湿度恒定时风速对干燥速率的影响㊂由图9可知,当热风温度㊁相对湿度相同且风速从0.2m/s增加到1.0m/s时,最大干燥速率提高了6.72%㊂当风速较低时,风速的增加对干燥速率的影响更显著,当风速较低时,风速的增加对干燥速率的影响更显著㊂当风速为0.2m/s时,坯体表面与热空气间的对流传质速率较低,坯体内部与表面的含水率梯度较小,随风速升高,坯体内部温度梯度增大,增强了传质驱动力㊂当风速从0.2m/s增加到1.0m/s时,干燥速率变化趋势相似,坯体表面水分已经蒸发,内部水分扩散速率依然较慢,在截面的同一含水率下干燥速率不会发生显著变化㊂图8㊀当热风温度㊁相对湿度恒定时风速对含水率的影响Fig.8㊀Influence of wind speed on water content with constant hot wind temperature and relativehumidity 图9㊀当热风温度㊁相对湿度恒定时风速对干燥速率的影响Fig.9㊀Influence of wind speed on drying rate with constant hot wind temperature and relative humidity3.2㊀热风温度对陶瓷坯体干燥过程的影响3.2.1㊀热风温度对坯体内部温度和含水率影响图10为干燥1h时不同热风温度下坯体表面和内部温度分布㊂从图10中坯体表面㊁x-y截面和x-z截面的温度云图分布可以看出:坯体迎风面温度较高,坯体背风面温度较低,且坯体表面温度高于内部温度,形成温度梯度,特别是在75ħ时坯体内温度梯度更大,35ħ时坯体内温度分布较均匀;随着热风温度升高,坯体内温度逐渐升高㊂垂直于流速方向的y-z截面上的温度对称分布,且坯体表面温度较高㊂坯体迎风面直接与热空气接触,在热风的风力下,对流换热率较大,坯体迎风面温度升高较快;热风温度越高,对流传热系数越大,坯体表面与热空气间的对流传热率越高,坯体表面温度越高,内部温度较低,形成较明显的温度梯度㊂图11为干燥2h时不同风温下坯体表面和内部含水率分布㊂由图11可知,在同一热风温度下坯体表面含水率与其内部含水率形成含水率梯度,且坯体迎风面含水率相对于背风面稍微较高㊂不同热风温度下,热风温度越高,坯体内部干基含水率越低,且当热风温度为75ħ时,坯体内部含水率较低,与热风温度从35ħ增加至55ħ相比,热风温度从55ħ增加至75ħ时相坯体内含水率降低相对较快㊂这是因为当热风温度较低时,热风温度的提高使得水分扩散速率得到较大的提高,水分活度提高较大,坯体内部水分迁移速率较大,干燥速率较快,坯体内含水率下降较快㊂图12为当热风速度㊁相对湿度恒定时风温对含水率影响㊂由图12可知,在干燥初期,随着温度升高,含水率下降越快,当温度为75ħ时,含水率下降最快,达到平衡所需时间最短,当温度从35ħ增加到75ħ时,平衡含水率减小了72.22%,这是因为温度增加使坯体内水分扩散系数增大,干基含水率下降速度增大㊂随着干燥时间延长,坯体内干基含水率缓慢下降,且干燥后期热风温度越高,干基含水率下降越缓慢,这是因为干燥后期坯体内水分较少,水分迁移量少,干燥速率较低㊂第8期霍㊀平等:基于Fick 扩散模型的陶瓷坯体干燥数值模拟2911㊀图10㊀干燥1h 时不同风温下坯体表面和内部温度分布Fig.10㊀Temperature distribution of surface and interior of green body at different temperatures drying for 1h 图11㊀干燥2h 时不同风温下坯体表面和内部含水率分布Fig.11㊀Water content distribution of surface and interior of green body at different temperatures drying for 2h 3.2.2㊀热风温度对干燥速率影响图13为当风速㊁相对湿度恒定时热风温度对干燥速率的影响㊂由图13可知,当热风速度㊁相对湿度相同,以及热风温度从35ħ增加到75ħ时,最大干燥速率变化幅度为46.34%㊂在相同干基含水率下,温度越高,干燥速率越大,随着干基含水率减小,温度越高,干燥速率下降速度越快,这是由于温度越高,水分扩散率越大,干燥速率随含水率的下降幅度较大㊂2912㊀陶㊀瓷硅酸盐通报㊀㊀㊀㊀㊀㊀第42卷图12㊀当风速㊁相对湿度恒定时热风温度对含水率的影响Fig.12㊀Influence of hot wind temperature on water content with constant wind speed and relativehumidity 图13㊀当风速㊁相对湿度恒定时热风温度对干燥速率的影响Fig.13㊀Influence of hot wind temperature on drying rate with constant wind speed and relative humidity3.3㊀热风相对湿度对陶瓷坯体干燥过程的影响3.3.1㊀相对湿度对坯体内部温度和含水率影响通过模拟可知,热风相对湿度对温度的影响较小,可以忽略不计㊂图14为不同相对湿度下5h 时坯体表面和内部含水率分布㊂由图14可知:在同一相对湿度下,坯体内部含水率较高,其表面含水率较低,热风相对湿度较大,坯体表面及内部含水率较高;随着热风相对湿度增加,坯体表面及内部含水率较大,这是因为坯体表面的饱和水蒸气分压增大,其表面水分蒸发率小,使得坯体表面含水率较高,坯体内部水分与表面的水分浓度梯度较小,其内部水分迁移较慢,则坯体内含水率也较高㊂图14㊀不同相对湿度下5h 时坯体表面和内部含水率分布Fig.14㊀Water content of distribution of surface and interior of green body at different relative humidity for 5h 当热风速度㊁温度恒定时相对湿度对含水率的影响如图15所示㊂由图15可知,当热风速度㊁热风温度相同时,热风相对湿度对陶瓷坯体内部干基含水率的影响较大;当热风速度㊁温度一定时,随着热空气相对湿度增大,含水率下降变慢,达到的平衡含水率最大,当相对湿度从5%增大到85%时,平衡含水率从0.8%增大至5.1%,这是因为坯体周围空气与陶瓷坯体表面的饱和水蒸气分压的压差减小,当空气中水蒸气浓度与坯体表面的水分浓度相等时,坯体表面水分不再蒸发,坯体内部水分不再向表面迁移,坯体内水分达到平衡㊂第8期霍㊀平等:基于Fick 扩散模型的陶瓷坯体干燥数值模拟2913㊀3.3.2㊀热风相对湿度对干燥速率影响图16为当风速㊁温度恒定时相对湿度对干燥速率的影响㊂由图16可知,当热风速度㊁温度相同,以及热风相对湿度从5%升高到85%时,最大干燥速率从11.1%/h 减小至8.8%/h,变化幅度为20.72%,干燥速率达到最大时的干基含水率最小,且干燥速率随干基含水率减小短暂升高后缓慢减小㊂在相同干基含水率下,相对湿度越高,湿度饱和差越小,干燥速率越小㊂在相同干燥速率下,相对湿度越高,空气中的绝对湿度越大,坯体干基含水率越高㊂图15㊀当风速㊁温度恒定时相对湿度对含水率的影响Fig.15㊀Influence of relative humidity on water content with constant wind speed andtemperature 图16㊀当风速㊁温度恒定时相对湿度对干燥速率的影响Fig.16㊀Influence of relative humidity on drying rate with constant wind speed and temperature4㊀结㊀论为揭示热风速度㊁温度㊁相对湿度对干燥过程中陶瓷坯体内部含水率和干燥速率的影响规律,为后续坯体干燥过程优化和质量提高提供理论参考,本研究基于Fick 扩散第二定律,对陶瓷坯体干燥过程进行了数值模拟,分析研究了热风速度㊁温度㊁相对湿度对陶瓷坯体热风干燥过程的影响,主要结论如下:1)在热风温度为55ħ㊁相对湿度为45%的条件下,随着风速增加,坯体内温度梯度增大,温度升高㊂提高风速能够加快坯体表面区域水分散失,当风速大于等于0.4m /s 时,风速对同一时间点下含水率影响较小㊂在生产中,将热风速度控制在合理的范围内能够防止坯体表面开裂㊂2)在热风速度为0.4m /s㊁相对湿度为45%的条件下,当温度从35ħ增加至75ħ时,温度梯度逐渐增大,最大干燥速率提高了46.34%,平衡含水率减小了72.22%㊂温度升高能够加快坯体内部干基含水率的变化,提高干燥速率,同时降低平衡含水率㊂3)在热风速度0.4m /s㊁温度为55ħ条件下,相对湿度对坯体内温度影响较小,相对湿度增大能够提高平衡含水率,当相对湿度从5%增大至85%时,平衡含水率从0.8%增大至5.1%㊂合理控制相对湿度可缩小坯体内部水分梯度,使坯体内水分分布更加均匀,有利于提高产品的良品率㊂参考文献[1]㊀左㊀元.论陶瓷坯体开裂的原因[J].蚌埠学院学报,2015,4(3):60-62.ZUO Y.On the reasons for cracking of the ceramic body[J].Journal of Bengbu University,2015,4(3):60-62(in Chinese).[2]㊀曾令可,税安泽.陶瓷工业实用干燥技术与实例[M].北京:化学工业出版社,2008:57-58.ZENG L K,SUI A Z.Practical drying techniques and examples in the ceramic industry[M].Beijing:Chemical Industry Press,2008:57-58(inChinese).[3]㊀吴海虹,朱道正,卞㊀欢,等.农产品干燥技术发展现状[J].现代农业科技,2016(14):279-281.WU H H,ZHU D Z,BIAN H,et al.Development status of agricultural products drying technology [J].Modern Agricultural Science andTechnology,2016(14):279-281(in Chinese).[4]㊀申常胜,李㊀镔,韦中华,等.陶瓷注塑成型技术的研究与进展[J].中国陶瓷工业,2022,29(5):24-29.SHEN C S,LI B,WEI Z H,et al.Research and development of ceramic injection molding[J].China Ceramic Industry,2022,29(5):24-29(in Chinese).2914㊀陶㊀瓷硅酸盐通报㊀㊀㊀㊀㊀㊀第42卷[5]㊀李㊀杰,朱永红,罗㊀添,等.辊道窑烧成带含坯体流场和温度场的数值模拟[J].中国陶瓷,2022,58(3):45-50.LI J,ZHU Y H,LUO T,et al.Numerical simulation of flow field and temperature field of gas and porcelain blank in the firing zone of a roller kiln[J].China Ceramics,2022,58(3):45-50(in Chinese).[6]㊀NAIT-ALI B,OUMMADI S,PORTUGUEZ E,et al.Thermal conductivity of ceramic green bodies during drying[J].Journal of the EuropeanCeramic Society,2017,37(4):1839-1846.[7]㊀李㊀杰,赖日东,吴志豪,等.辊道窑排烟段含坯体流场和温度场特性的数值研究[J].陶瓷学报,2022,43(2):317-322.LI J,LAI R D,WU Z H,et al.Numerical study on characteristics of flow field and temperature field of roller kiln containing porcelain blank in the smoke exhaust section[J].Journal of Ceramics,2022,43(2):317-322(in Chinese).[8]㊀LAURO N,OUMMADI S,ALZINA A,et puter model of drying behaviour of ceramic green bodies with particular reference to moisturecontent dependent properties[J].Journal of the European Ceramic Society,2021,41(14):7321-7329.[9]㊀赖日东.陶瓷坯体干燥过程温湿度变化数值模拟与实验研究[D].景德镇:景德镇陶瓷学院,2015:27-37.LAI R D.Numerical simulation and experimental investigation of ceramic body temperature and humidity variation in drying process[D].Jingdezhen:Jingdezhen Institute of Ceramics,2015:27-37(in Chinese).[10]㊀朱庆霞,罗民华,梁华银.大规格建陶坯体对流干燥热质耦合传递的数值模拟[J].陶瓷学报,2009,30(4):503-507.ZHU Q X,LUO M H,LIANG H Y.Heat and mass transfer inside porcelain tiles during convective drying process[J].Journal of Ceramics, 2009,30(4):503-507(in Chinese).[11]㊀DEFRAEYE T,BLOCKEN B,CARMELIET J.Analysis of convective heat and mass transfer coefficients for convective drying of a porous flatplate by conjugate modelling[J].International Journal of Heat and Mass Transfer,2012,55(1/2/3):112-124.[12]㊀王振华.仓储粮堆湿热传递过程的数值模拟与试验研究[D].北京:中国农业大学,2014:26-39.WANG Z H.Numerical simulation and experimental study on heat and mass transfer of stored grain pile[D].Beijing:China Agricultural University,2014:26-39.(in Chinese)[13]㊀RANJBARAN M,EMADI B,ZARE D.CFD simulation of deep-bed paddy drying process and performance[J].Drying Technology,2014,32(8):919-934.[14]㊀DA SILVA W P,DE OLIVEIRA F V S,DE ARAÚJO N G,et al.Modeling of water transport in roof tiles by removal of moisture at isothermalconditions[J].Heat and Mass Transfer,2012,48(5):809-821.[15]㊀毕清跃.多场耦合的多层玉米籽粒热风干燥特性研究[D].吉林:东北电力大学,2021:20-28.BI Q Y.Multi-field coupled hot air drying characteristics of multilayer corn seeds[D].Jilin:Northeast Dianli University,2021:20-28(in Chinese).[16]㊀ARAÚJO M V,SANTOS R S,DA SILVA R M,et al.Drying of industrial hollow ceramic brick:analysis of the moisture content and temperatureparameters[J].Defect and Diffusion Forum,2017,380:72-78.[17]㊀TEIXEIRA D B M K,TEIXEIRA D A D B,BARBOSA D L A G,et al.Heat and mass transfer during drying of clay ceramic materials:a three-dimensional analytical study[J].Diffusion Foundations,2017,10:93-106.[18]㊀陈衍男,王㊀晓,穆㊀岩,等.天麻蒸制后红外干燥特性及失水动力学研究[J].食品工业科技,2018,39(22):30-34+40.CHEN Y N,WANG X,MU Y,et al.Drying characteristics and kinetics research of gastrodia elata blume under infrared blast drying after steaming[J].Science and Technology of Food Industry,2018,39(22):30-34+40(in Chinese).[19]㊀欧阳梦云,王㊀燕,罗凤莲,等.预处理及干燥方法对果蔬干制品VC含量影响的研究进展[J].食品与机械,2018,34(7):179-182.OUYANG M Y,WANG Y,LUO F L,et al.Advances in vitamin C content in dried fruits and vegetables by using pretreatment and drying methods[J].Food&Machinery,2018,34(7):179-182(in Chinese).[20]㊀ZHU A S,SHEN X Q.The model and mass transfer characteristics of convection drying of peach slices[J].International Journal of Heat andMass Transfer,2014,72:345-351.[21]㊀ELMIZADEH A,SHAHEDI M,HAMDAMI parison of electrohydrodynamic and hot-air drying of the quince slices[J].Innovative FoodScience&Emerging Technologies,2017,43:130-135.[22]㊀OJEDIRAN J O,OKONKWO C E,ADEYI A J,et al.Drying characteristics of yam slices(Dioscorea rotundata)in a convective hot air dryer:application of ANFIS in the prediction of drying kinetics[J].Heliyon,2020,6(3):e03555.[23]㊀ZHAO J H,FENG S,GRUNEWALD J,et al.Drying characteristics of two capillary porous building materials:calcium silicate and ceramic brick[J].Building and Environment,2022,216:109006.Copyright©博看网. All Rights Reserved.。
液态金属凝固中的传热、传质及液体流动
t R2 K2
K为凝固系数。
在实际的生产中,通常不需计算出铸件的凝固时间, 只需通过比较它们的相对厚度或模数就可制定生产工艺。
铸件温度场及凝固时间的精确计算——计算机数值模拟
4、焊接温度场
准稳定温度场的概念
图4-4 “厚板”表面运动点热源的温度场
图4-5 薄板焊接时的温度场分布
(图b是否有误?)
3.数值计算法 数值计算法是把所研究的物体从时间和位置上分割成许多小
单元,对于这些小单元用差分方程式近似地代替微分方程式, 给出初始条件和边界条件,逐个计算各单元温度的一种方法。 即使铸件形状很复杂,也只是计算式和程序烦杂而已,在原则 上都是可以计算的。
数值计算法比其它近似计算法准确性高,当单元选得足够小
无限长圆棒试样 测温及结果处理
2.铸件的两种凝固方式
图4-3 合金成分和温度梯度对凝固方式的影响 a)、b)为层状凝固,c)、d)为体积凝固 影响因素:(1)化学成分(液-固相线距离)
(2)温度梯度
层状凝固过程 层状凝固缩孔特点
体积凝固过程 体积凝固方式的缩松
3、铸件凝固时间计算
——与铸件厚度及温度场(凝固速度)相关
1)铸型的蓄热系数 铸型的蓄热系数越大,对铸 件的冷却能力就越大,铸件是的温度梯度就越大。 铸型的导热系数越大,能把铸型内表面吸收的热迅 速传至外表面,使铸型内表面保持强的吸热能力, 铸件内的温度梯度也就大。如金属型、涂料等的影 响。
2)铸型的预热温度的影响 铸型预热温度越高,对铸件的冷却作用就越小,铸件断面上的温度梯度也 就越小。 3.浇注条件的影响 过热热量加热了铸型,所以过热度越大,相当于铸型预热温度越高。铸件 内的温度场越平坦。 4.铸件结构的影响 1)铸件的壁厚 厚壁铸件比薄壁铸件含有更多的热量,当凝固层向中心推 进时,把铸型加热到更高温度,所以铸件内温度场较平坦。 2)铸件的形状 铸件的棱角和弯曲表面,与平面的散热条件不同,向外凸 出的部分,散出的热量为较大何种的铸型所吸收,铸件的冷速较大,如果铸 件内凹的表面,则相反。
传热和传质基本原理 第四章 三传类比
4.2.2 柯尔本类似律
雷诺类似律或忽略了层流底层的存在,普朗特正 对此进行改进,推导出普朗特类似律:
冯卡门认为紊流核心与层流底层之间还存在一个 过渡层,于是又推导出了卡门类似律:
契尔顿和柯尔本根据许多层流和紊流传质的实验结果, 在1933年和1934年,得出:
简明适用,引入了流体的 重要物性Sc数。
24
根据薄膜理论,通过静止气层扩散过程的传质系数可定义为:
25
在紧贴壁面处,湍动渐渐消失,分子扩散起主导 作用,在湍流核心区,湍流扩散起主导,传质系 数与扩散系数成下列关系
另外,δ的数值决定于流体的流动状态,即雷诺 数。
26
4.4.2
同一表面上传质过程对传热过程的影响
设有一股温度为t2 的流体流经温度为t1的壁面。传递过程 中,组分A、B从壁面向流体主流方向进行传递,传递速 率分别为NA、NB。可以认为在靠近壁面处有一层滞留薄 层,假设其厚度为y0 ,求壁面与流体之间的热交换量。
边界层厚度
1904年普朗特首先提出
39
4.5.1 边界层理论的基本概念
边界层的定义
流体在绕过固体壁面流动时,紧 靠固体壁面形成速度梯度较大的 流体薄层称为流动边界层 流速相当于主流区速度的0.99处到固 体壁面间的距离定义为边界层的厚度
边界层的形成与特点
Re x 2 10 5
vl Re
以此两式计算管内流动质交换系数结果很接近。
17
18
紊流
19
例题: 试计算空气沿水面流动时的对流质交换系数hm和每 小时从水面上蒸发的水量。已知空气的流速 u=3m/s,沿气流方向的水面长度l=0.3m,水面的温 度为15 ℃ ,空气的温度为20 ℃ , 空气的总压力 1.013*105Pa,其中水蒸汽分压力p2=701Pa,相当 于空气的相对湿度为30%。
金属凝固理论 第4章 液态金属凝固过程中的传热、传质及液体流动
金属的凝固温度越高,在凝固过程中铸件表面和铸型内表 面的温度越高,铸型内外表面的温差就越大,致使铸件断 面温度场出现较大的梯度。如有色金属与钢铁相比,其温 度场较平坦。
2020/10/15
11/56
(2)铸型性质的影响 1)铸型的蓄热系数
铸型的蓄热系数越大,对铸件的冷却能力就越大, 铸件内的温度梯度就越大。铸型的导热系数越大, 能把铸型内表面吸收的热迅速传至外表面,使铸 型内表面保持强的吸热能力,铸件内的温度梯度 也就大。
向中心推进时,把铸型加热到更高温度,所以铸 件内温度场较平坦。
2)铸件的形状 铸件的棱角和弯曲表面,与平面的散热条件不同。
向外凸出的部分,散出的热量被较大体积的铸型 所吸收,铸件的冷速较大,如果铸件内凹的表面, 则相反。
2020/10/15
13/56
三、铸件的凝固方式及影响因素
TL
TS
L
S+L S
第四章 液态金属凝固过程 中的传热、传质及液体流动
2020/10/15
1/56
第一节 凝固过程中的传热
在材料成形过程中,液态金属的过热热量和 凝固潜热主要是通过传导而释放的。
2020/10/15
2/56
一、凝固过程中的热传导及傅里叶方程
温度场基本概念: 稳定温度场: 不随时间而变的温度场(即温度只是坐标的函 数),其表达式为:
2020/10/15
T n
Tw Tf
4/56
凝固过程中,热量传递有三种形式:传导、辐射、 对流。
以热传导为主。 热传导过程取决于温度的分布——温度场:温度
随空间和时间的变化。 T = f(x,y,z,t) Fourier热传导方程:
CFD技术在“数值传热学”教学中的应用
CFD技术在“数值传热学”教学中的应用作者:杜敏来源:《中国电力教育》2012年第33期摘要:提出将先进的CFD技术引入“数值传热学”的教学中,生动、形象地展示各种数值方法的计算结果,将抽象的概念、理论变成形象的画面,并结合基础理论进行讲解,便于学生对所学内容的深入理解,从而达到激发学生的学习兴趣,改善教学效果的目的。
关键词:数值传热学;CFD技术;Fluent软件;教学改革作者简介:杜敏(1981-),男,山东泰安人,江苏大学能源与动力工程学院,讲师;王助良(1964-),男,湖北鄂州人,江苏大学能源与动力工程学院,教授。
(江苏镇江212013)基金项目:本文系江苏大学高级专业人才科研启动基金(项目编号:11JDG152)的研究成果。
中图分类号:G642.0 文献标识码:A 文章编号:1007-0079(2012)33-0058-02数值传热学,[1]又称计算传热学,是指对描写流动与传热问题的控制方程采用数值方法,通过计算机求解的一门传热学与数值方法相结合的交叉学科。
本课程的目的是使学生掌握一种能够预测传热与传质、流体流动过程的数值方法,用以解决工程实际中大量存在的且用解析方法难以解决的传热与流体流动问题。
但由于数值传热学比较抽象、枯燥、难懂,应用的数学知识较多,大部分学生对该课程都有某种程度的畏惧感。
计算流体动力学(Computational Fluid Dynamics,CFD)技术[2]是利用数值方法通过计算机求解描述流体流动的数学方程,获得空间和时间离散位置处的数值解,揭示流动的物理规律和研究流动的物理特性的学科,目前已成为国际上一个强有力的研究领域,[3]是进行传热、传质、动量传递及燃烧、多相流和化学反应研究的核心和重要技术,广泛应用于航天设计、汽车设计、生物医学工业、化工处理工业、涡轮机设计、半导体设计等诸多工程领域。
随着计算机技术不断发展与进步,一些教学工作者已尝试将CFD技术应用到“工程流体力学”、“传热学”、“空气调节”等课程的教学中,[4-6]并取得了良好的教学效果。
传热和传质基本原理-----第四章-三传类比
相当于空气的相对湿度为30%。
38
4.5 边界层类比
流体流动的控制方程是非线性的偏微分方程组,处理 非线性偏微分方程依然是当今科学界的一大难题
实际工程问题:靠近固体 壁面的一薄层流体速度变 化较大,而其余部分速度 梯度很小
➢ 远离固体壁面,视为理想流 体--欧拉方程、伯努利方程
➢ 靠近固体壁面的一薄层流体, 进行控制方程的简化--流动 边界层
27
❖ 在薄层内取一微元体,那么进入微元体的热流为 由温度梯度引起的导热热流、由进入微元体的传 递组分本身具有的焓。
稳定状态时,微元体处于热平衡,满足下列关系式:
令
无因次数为传质阿克曼修正
系数,表示传质速率的大小、
方向对传热的影响。
28
得 边界条件为
令
得方程的解为:
代入边界条件,最后得到流体在薄层内的温度分别为:
水蒸 汽的汽化潜热r=2463.1kJ/kg,Sc=0.6.,Pr=0.7。 试计算干空气的温度。
2.试计算空气沿水面流动时的对流质交换系数hm和每小时从 水面上蒸发的水量。已知空气的流速u=3m/s,沿气流方向
的
水面长度l=0.3m,水面的温度为15 ℃,空气的温度20℃,
空气的总压力1.013*105Pa,其中水蒸汽分压力p2=701Pa,
➢边界层厚度
1904年普朗特首先提出
39
4.5.1 边界层理论的基本概念
边界层的定义
流体在绕过固体壁面流动时,紧 靠固体壁面形成速度梯度较大的 流体薄层称为流动边界层
流速相当于主流区速度的0.99处到固 体壁面间的距离定义为边界层的厚度
边界层的形成与特点 Re vl
平板绕流
Re x
v0 x
流动、传热及传质的控制方程
对流传质方程
01
对流传质方程是描述流体流动过程中物质传递的方程,它基于Fick第二定律和 Darcy定律。
02
对流传质方程的一般形式为:ρSc▽·vc = -▽P/ρ + ν▽²vc + (1/ρ)▽·(ρD▽c),其中 Sc是斯密特数,v是速度矢量,P是压力,ν是动力粘度,D是扩散系数,c是浓度。
有限元法在结构分析、固体力学、流体力学等领域有广泛应用。
有限体积法
01
02
03
有限体积法是一种求解偏微分方 程的数值方法,它将连续的求解 区域离散化为有限个小的体积单 元,并对每个体积单元构造近似 函数。
有限体积法特别适合处理流体动 力学问题,因为它能够很好地捕 捉到流体运动的特性,如速度和 压力的连续性。
熵守恒方程
总结词
描述流体熵在空间中的变化。
详细描述
熵守恒方程是热力学的基本方程之一,它表明在封闭系统中,流体的熵不会凭空产生或 消失。该方程基于熵增原理,表示流体在流场中单位时间内熵的增加等于流入该控制体
的净熵流量。
02 传热的控制方程
热传导方程
总结词
描述了物体内部热量的传递过程。
详细描述
热传导方程,也称为傅里叶定律,表 示在物体内部,温度梯度导致热量从 高温区域流向低温区域。该方程基于 能量守恒原理,并考虑了导热系数的 影响。
03
对流传质方程适用于描述流体流动过程中的物质传递过程,如化工、环境、食 品等领域中的流动和传递过程。
化学反应动力学方程
1
化学反应动力学方程是描述化学反应速率和反应 机理的数学模型。
2 3
1-1 流体静力学基本方程式
第一章流体流动本章学习指导1. 本章学习的目的通过本章学习,掌握流体流动过程的基本原理、管内流动的规律,并运用这些原理和规律去分析和计算流体流动过程的有关问题,诸如:(1)流体输送:流速的选择,管径的计算,输送机械选型。
(2)流动参数的测量:压强(压力)、流速(流量)等。
(3)不互溶液体(非均相物系)的分离和分散(混合)。
(4)选择适宜的流体流动参数,以适应传热、传质和化学反应的最佳条件。
2. 本章重点掌握的内容(1)静力学基本方程的应用(2)连续性方程、柏努力方程的物理意义、适用条件、应用柏努力方程解题的要点和注意事项。
(3)管路系统总能量损失方程(包括数据的获得)本章应掌握的内容(1)两种流型(层流和湍流)的本质区别,处理两种流型的工程方法(解析法和实验研究方法)(2)流量测量(3)管路计算本章一般了解的内容(1)边界层的基本概念(边界层的形成和发展,边界层分离)(2)牛顿型流体和非牛顿型流体3. 本章学习应注意的问题(1)流体力学是传热和传质的基础,它们之间又存在着密切的联系和相似性,从开始学习流体流动就要学扎实,打好基础。
(2)应用柏努力方程、静力学方程解题要绘图,正确选取衡算范围。
解题步骤要规范。
4. 本章教学时数分配知识点1-1 授课学时数2 自学学时数4知识点1-2 授课学时数3 自学学时数6知识点1-3 授课学时数1 自学学时数2知识点1-4 授课学时数3 自学学时数6知识点1-5 授课学时数1 自学学时数2知识点1-6 授课学时数2 自学学时数45. 本章学习资料必读书籍姚玉英主编. 化工原理(上册) (第一章"流体流动")·天津:天津大学出版社.1999参考书籍1.陈敏恒等.化工原理,上册.北京:化学工业出版社.19992.谭天恩等.化工原理,上册.北京:化学工业出版社.19903.蒋维钧.化工原理,上册.北京:清华大学出版社.19924.姚玉英.化工原理例题与习题,第三版.北京:化学工业出版社.19985.柴诚敬等.化工原理学习指导.天津:天津科技出版社.19926.柴诚敬,张国亮.化工流体流动和传热.北京:化学工业出版社.20007.张言文.化工原理60讲,上册.北京:轻工业出版社.19978.J.M.Coulson and J.F.Richrdson.Chemical Engineering Vol2.3rd ed.-oxford:Pergamon,19949.C.J Geankoplis. Transport Processes and Unit Operations, 2rd ed.Boston: Allyn and Baccon, Inc. 199310.W. L. McCabe and J. C. Smith.Unit Operations of Chemical Engineering, 5th ed. New York: McGraw. Hill Inc., 1993.通过本章的学习,掌握气体吸收的基本概念和气体吸收过程的基本计算方法。
第一章 1[1].1流体流动静力学基本方程
第一章 1[1].1流体流动静力学基本方程第一章流体流淌1-0 概述一学习本章的意义:1.流体存在的广泛性。
在化工厂中,管道和设备中绝大多数物质都是流体(包括气体、液体或气液混合物)。
只是到最后,有些产品才是固体。
2 .通过讨论流体流淌逻辑,可以正确设计管路和合理挑选泵、压缩机、风机等流体输送设备,并且计算其所需的功率。
3 .流体流淌是化工原理各种单元操作的基础,对强化传热、传质具有重要的实践意义。
由于热量传递,质量传递,以及化学反应都在流淌状态下举行,与流体流淌密切相关。
所以大家要仔细学习这一章,充分打好基础。
二流体流淌的讨论范畴1 流体定义:具有流淌性的液体和蔼体统称为流体。
2 延续性介质假定:流体是由大量的单个分子组成,而每个分子之间彼此有一定的间隙,它们将随时都在作无规章随机的运动。
所以,若把流体分子作为讨论对象,则流体将是一种不延续介质,这将使讨论十分困难。
好在在化工生产过程中,我们对流体流淌逻辑的讨论感爱好的并非是单个分子的微观运动,而是流体宏观的机械运动。
所以我们不取单个分子作为考察对象,而取比分子平均自由程大得多,比设备尺寸小得多的这样一个流体质点作为最小考察对象,质点是由大量分子组成的微团,它可以代表流体的性质。
流体可以看成是由大量微团组成的,质点间无空隙,而是弥漫所占空间的延续介质,从而可以使用延续函数的数学工具对流体的性质加以描述。
提高:延续性介质假定如图1所示,考虑一个微元体积内流体平均密度的变化状况:取包含P(x,y,z)点在内的微元体积⊿V,其中包含流体的质量为⊿m,则微元流体的平均密度为⊿m/⊿V,微元流体的平均密度随体积的变化如图2所示。
当微元体积⊿V从十分小逐渐增大,趋向一个特定的微元体积V时,流体的平均密度逐渐趋向一个极限值,且不再随微元体积的继续增大而发生变化。
当微元体积⊿V比δV小时,这时微元体积内所包含的流体分子数目是那样少,以致流体分子因为其无规章的热运动,进入或离开微元体积的流体分子数目已足以引起该微元体积内流体平均密度的随机波动。
流动传热及传质的控制方程
缺点: ①数学模化的全面和准确性需要不断提高:
Ⅰ、物理问题的数学模型是否正确(回流问题还是边界层问题, 稳态还是非稳态),否则,数值算法的改进没有意义。
Ⅱ、所有物性数据要可靠,否则减少数值误差的努力毫无意义。 ②真实再现某些过程的代价也是极其昂贵的或不可能;(用于气象,
石油) ③有些迫不得已的简化是致命的或大大降低其价值; ④计算结果准确性仍需接受实验或精确解检验。(如对有代表性点
把原来在空间与时间坐标中连续的物理量的场速度场温度场浓度场等用一系列有限个离散点节点上的值的集合来代替通过一定的原则建立起这些离散点上变量值之间关系的代数方程称为离散方程求解所建立起来的代数方程以获得所求解变量的近似值
流动与传热的数值计算
§1 绪论
1.1 引言 1、传热、传质与流体流动的重要性
工程设备(如结晶器,中间包,钢包及锅炉,高炉等) 内部流体流动及热交换过程,自然环境中的污染问题,暴 风雨雪,河流泛滥及着火过程中出现的热、质传递,流动 起着重要作用。 2、对过程估计和认识的必要性
一.质量守恒方程(连续性方程)
1.理论依据:质量守恒定律 2.数学描述: [单位时间内微元体中流体质量的增加]=[ 同一时间间隔内流入该微
元体的净质量] 3.数学表达式:
?? ? ? ?? u?? ? ?? v?? ? ?? w?? 0
?t ?x
?y
?z
?? ? div(? U) ? 0
?t
or
?? ? ? (? U) ? 0
五.控制方程的通用形式 引入背景:比较四个基本控制方程式,虽因变量各不相同,但它
)
?
div(?
gradu)
?
?p ?x
?
Su
化工原理伯努利方程
伯努利方程流体宏观运动机械能守恒原理的数学表达式。
1738年瑞士数学家D.伯努利在《水动力学──关于流体中力和运动的说明》中提出了这一方程。
它可由理想流体运动方程(即欧拉方程)在定态流动条件下沿流线积分得出;也可由热力学第一定律导出。
它是一维流动问题中的一个主要关系式,在分析不可压缩流体的定态流动时十分重要,常用于确定流动过程中速度和压力之间的相互关系。
方程的形式 对于不可压缩的理想流体,密度不随压力而变化,可得:Zg+22u P +ρ=常数式中Z 为距离基准面的高度;P 为静压力;u 为流体速度;ρ为流体密度;g 为重力加速度。
方程中的每一项均为单位质量流体所具有的机械能,其单位为N ·m/kg ,式中左侧三项,依次称为位能项、静压能项和动能项。
方程表明三种能量可以相互转换,但总和不变。
当流体在水平管道中流动时Z 不变,上式可简化为:ρPu +22=常数 此式表述了流速与压力之间的关系:流速大处压力小,流速小处压力大。
对于单位重量流体,取管道的1、2两截面为基准,则方程的形式成为:gu g P Z g u g P Z 2222222111++=++ρρ 式中每一项均为单位重量流体的能量,具有长度的因次,三项依次称为位头、静压头和动压头(速度头)。
对于可压缩理想流体,密度随压力而变化。
若这一变化是可逆等温过程,则方程可写成下式:1211222211ln 22P PP u gZ u gZ ρ++=+若为可逆绝热过程,方程可写为:1211222211ln 22P PP u gZ u gZ ρ++=+式中γ为定压比热容Cp 和定容比热容Cv 之比,即比热容比,也称为绝热指数。
对于粘性流体,流动截面上存在着速度分布,如用平均流速u 表达动能项,应对其乘以动能校正系数d ο。
此外,还需考虑因粘性引起的流动阻力,即造成单位质量流体的机械能损失h f ,若在流体流动过程中,单位质量流体又接受了流体输送机械所做的功W ,在这些条件下,若取处于均匀流段的两截面1和2为基准,则方程可扩充为:α值可由速度分布计算而得, 流体在圆管内作层流流动时α=2;作湍流流动时,α≈1.06。
流体力学三大方程公式及符号含义
流体力学是研究流体运动和力学的学科,涉及流体的运动规律、压力、密度等物理性质。
在流体力学的研究中,三大方程公式是非常重要的理论基础,它们分别是连续方程、动量方程和能量方程。
本文将对这三大方程公式及其符号含义进行详细介绍。
一、连续方程连续方程是描述流体连续性的重要方程,它表达了流体在运动过程中质点的连续性。
连续方程的数学表达式为:\[ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 \]其中,符号和含义说明如下:1.1 ∂ρ/∂t:表示密度随时间的变化率,ρ为流体密度。
1.2 ∇·(ρv):表示流体质量流动率的散度,∇为Nabla算子,ρv为流体的质量流速矢量。
这一方程表明了在运动的流体中,质量是守恒的,即单位体积内的质量永远不会减少,这也是连续方程的基本原理。
二、动量方程动量方程描述了流体运动过程中动量的变化和传递,是流体力学中的核心方程之一。
其数学表达式为:\[ \frac{\partial (\rho \mathbf{v})}{\partial t} + \nabla \cdot (\rho \mathbf{v} \mathbf{v}) = -\nabla p + \nabla \cdot \mathbf{\tau} + \mathbf{f} \]其中,符号和含义说明如下:2.1 ∂(ρv)/∂t:表示动量随时间的变化率。
2.2 ∇·(ρv⃗v):表示动量流动率的散度。
2.3 -∇p⃗:表示流体受到的压力梯度力。
2.4 ∇·τ⃗:表示应力张量的散度,τ为流体的粘性应力张量。
2.5 f⃗:表示单位体积内流体受到的外力。
动量方程描述了流体内部和外部力之间的平衡关系,它是研究流体运动规律和动力学行为的重要方程。
三、能量方程能量方程描述了流体在运动过程中的能量变化规律,包括内能、压力能和动能等能量形式。
化工原理重要单元主要公式汇总
化工原理课程综合温习提纲化工原理重要单元主要公式汇总第1章 流体流动一、机械能衡算方程式 本章内容的核心公式是机械能衡算方程式:g 2ud L g 2u g P Z H g 2u g P Z 22222e 2111⎪⎭⎫ ⎝⎛++++=+++∑ζλρρ (单位:J/N=m ) (1-1)应用公式(1-1)注意以下几点:(1) 稳固流动、不可紧缩性流体、自1-1至2-2的控制体内流体持续。
(2) Z 1、Z 2选择同一水平基准面,通常选择地平面或控制体1-一、2-2中的较低的一个。
(3) P 1、P 2同时以绝对压计或同时以表压计,而且注意单位均统一到N/m 2 。
(4) 自高位槽或高压容器向其他地方输送流体时一般不需要流体输送机械,现在,H e =0 。
(5) 公式中的每一项均是单位流体的能量,每牛顿流体的能量焦耳,形式上的单位是米。
H e 是流体输送机械加给每牛顿流体的能量焦耳数,阻力损失项亦是每牛顿流体的能量损失焦耳数。
(6) 按照所取的1-一、2-2截面的性质,灵活地肯定u 1、u 2的数值。
(7) 阻力损失项中的流速取产生阻力损失的管段上的流速,有时管段不止一段。
(8) 若控制体内的阀门关闭,1-一、2-2截面上的流体能量便再也不有任何关系。
(9) 若在等直径的管段,无流体输送机械,阻力损失能够忽略,(1-1)式变成流体静力学的形式。
应用公式(1-1)可解决以下方面的问题:(1) 在肯定的控制体中,达到必然的流量,肯定流体输送机械加给每牛顿流体的能量焦耳数及功率。
(2) 在肯定的控制体中,达到必然的流量,肯定起始截面1-1的高度或压强。
(3) 在肯定的控制体中,可达到的流量(流速)。
(4) 在肯定的控制体中,达到必然的流量,肯定管径。
公式(1-1)的另两种形式:2ud L 2u P g Z w 2u P g Z 22222e 2111⎪⎭⎫ ⎝⎛++++=+++∑ζλρρ (单位:J/kg )(1-2)ρζλρρρρρ2udL2u P g Z g H 2u P g Z 22222e 2111⎪⎭⎫ ⎝⎛++++=+++∑ (单位:J/m 3=N/m 2) (1-3)因为机械能衡算式中的每一项均是单位流体的能量,故计算流体输送机械的功率时应注意流体的总流量V q (单位:m 3/s)。
流动、传热及传质的控制方程
5、当前数值计算的发展趋势 1)坐标系的选取 正交坐标系:三个方向均是正交的。如飞机头部变化剧烈 部分在飞机表面上,如何划分网格。正交坐标系→非正交坐 标系。(二阶偏导由正交坐标系→非正交坐标系转换要增加 几十项) 2)非稳态现象的模拟仿真 如一钢包浇注过程中,内部钢水流动处于非稳态过程,数值 模拟计算得到,液面下降到200mm处不能下渣。
(U) 0
t
(1)
适用条件:可压缩和可压缩流体,理想流体和实际流体,稳态及非稳态流动。
二.动量守恒方程
1.理论依据:动量守恒定律(牛顿第二定律)
2.数学描述:
[微元体中流体的动量对时间的变化率]=[外界作用在该微元 体上的各种力之和]
3.数学表达式:
DU 1 p f 2U Dt
1)空间坐标一般是双向坐标,而时间坐标则总是单向坐标; 2)空间坐标也可能作为单向坐标;
如在一个坐标方向上有很强的单向流动,则信息的传递或影 响只能是从上游传至下游,即某点上的状态主要受其上游条 件的影响。(如强制对流过程) 3)数学及物理状态上的各自描述:
数学上微分方程中只要有非稳态项→抛物线方程问题→单向 坐标(无论有无空间上的);
div(U ), zz
2
w z
div(U )
xy
yx
u y
v x
,
xz
zx
u z
w x
,
yz
zy
v z
w
y( 4)
其中:μ是动力粘度,λ是第二粘度,一般可取-2/3。
热质交换原理与设备完整版 修订版
JH JD
cf 2
2 可以把对流传热中有关的计算式用于对流传质, 只要 St pr 3 Stm Sc 3 ○
2
2
将对流 传热 计算式 中的 有关 物理参 数及 准则 数用于 对流 传质 中相应 的代 换即 可,如
3 同样可以用类比关系由传热系数 h 计算传质系 t↔c,a↔D, λ ↔D,pr↔Sc,Nu↔Sh,St↔ Stm ○
2 干燥循环的过程: ○ 干燥剂的吸湿和放湿是由干燥剂表面的蒸汽压与环境空气的蒸汽压差造 成的, 当前者较低时, 干燥剂吸湿, 反之放湿, 两者相等时达到平衡, 即既不吸湿也不放湿;
完整的干燥循环由吸附过程,脱附过程和冷却过程构成 3 影响吸收效果的因素:A 除湿剂的结构 B 除湿剂的选择 ○ 第六章
4. 喷淋室的热交换效率系数: 1 (第一热交换效率或全热交换效率) 1 1
ts 2 tw 2 热湿 ts1 tw1
交换越不完善,1 越小( t s1 、 t s 2 空气始终态的湿球温度, t w1 、t w 2 水的始终态的湿球温度) 喷淋室的接触系数 2 (第二热交换效率或通用热交换效率) 2 1 变化完善程度。 5. 喷淋式计算的主要原则: 该喷淋室能达到的1 应等于空气处理过程所需要的1 , 该喷淋室 能达到的 2 应等于空气处理过程需要的 2 ,该喷淋室喷出的水能吸收(或放出)的热量应 等于空气失去(或得到)的热量。
数 hm 3、对流传质过程的准则数:施密特准则数: Sc
v ,运动黏度与物体扩散系数之比;表 Di
示物性对对流传质的影响,速度与浓度边界层的相对宣乌特数 Sh
hml ;斯坦顿数: Di
St
a hm ;刘伊斯准则: Le , 表示温度分布和浓度分布关系的相互关系,体现传热和传 D u
传热与传质最全的计算
qmhCph ( T1 – T2) = qmcCpc ( T2′– T1′)
其中qmh和qmc分别为热、冷流体的质量流量,kg/s或kg/h
二、传热速率方程
1、总传热速率微分方程
dA
通过换热器中任一微元面积dA的间壁 T
两侧流体的传热速率方程,可以仿照对
dφ
流传热速率方程写出:
dφ = K ( T – T′) dA = K △T dA
i1 i
ri
ΔT=T1 –Tn+1
Φ
r2
r1
r3
T1
r4
T2
T4 T3
5、保温层的临界半径
t1----保温层内表面温度;tf----环境温度 r1、r2----分别为保温层内外壁半径; λ---为保温材料的导热系数 α----为对流传热系数;L---为管长
t1 t f
t1 t f
R1 R2
▪ 蒸发:碱液蒸发、PVC干燥
传热在生产中的应用
▪ 2、化工设备和管道的保温(保冷),以减少热量 (冷量)损失。
▪ 保温:如蒸汽管道、热水管道。 ▪ 保冷:-35℃盐水、7℃水管道 ▪ 3、生产中热能的合理利用,废热回收。 ▪ 废热利用:氯化氢合成热用于溴化锂及采暖、转
化反应热用于溴化锂机组
▪ 研究传热的目的
传热的分类
▪ 传热的分类
间歇传热
按连续性
分类:
连续传热 非稳态传热:传热速率常数,
按与时间 的关系
稳态传热:传热速率=常数,
LOGO
一、热量传递的三种基本方式
根据传热的机理不同,热量传递的基本方式分为三种:
导热 对流 热辐射
1、热传导(又称导热)
燃烧学-五-多组分反应流体守恒方程
Cp
''' , S h0 q f ,i m ,i
六、守恒标量的概念
1.简单化学反应模型 化学反应:燃料和氧化剂消失,产生二氧化碳和 水蒸气,燃气温度升高并发出热量。 假设: (1)燃料和氧化剂以化学恰当比进行单步不可逆 反应,生成单一的燃烧产物 1kg燃料+ kg 氧化剂 (1 )kg 产物 (a) (2)各组分的传输特性相同,但可以随空间位置 而变化(每处每参数相等,但可不均匀);
vi vi ,diff v
组分总的质量通量等于对流通量和扩散通量之和,即
组分总的质量流量
q
'' m ,i
q m9;'
扩 散 通 量
m ,i ,diff
对流通量
wi v wi vi ,diff
将(c)式代入式(a),得
( c)
代入分子输运的费克扩散定律,得
H MW H H O MW H O
( wi ) '' ''' qm ,i qm ,i t
由
'' m ,i
组分
i 的质量守恒方程的一般矢量形式为
,
'' m
i 1,2,, N
,得
( a)
q wi vi q v
混合物质量平均速度
v wi vi
(b)
组分速度等于质量平均速度叠加上扩散(布朗运动)速度
( pA) x ( pA) x x qm [(v) xx (v) x ]
上式除以 A,并取极限 x
x 0
,得
或
dp '' dv qm dx dx dp dv vx dx dx
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
例:炉子结构参数、热工操作参数→热工过程参数→炉子的生 产指标 (质量,产量,热耗,寿命)
人类对过程的认识来自于对它们的掌握和控制的需求, 以最小的代价进行高效率的产品生产,长期的发展已创造出 一套定量估计过程的方法。 3、什么是预测或估计? 估计和预测是人们利用已有的知识和手段对未知过程或现象 进行过程的超前了解的方法。(如设计中间包,结构外形已 定,如何设置内部控制元件使之达到理想的环流轨迹,找出 参量间的关系。)
2.2 数值模拟步骤 1、建立反映工程问题或物理过程本质的数学模型
即建立反映问题各量间关系微分方程及相应定解条件。 数学模型:即基本定律(三大基本守恒定律及组元守恒定律)+ 模型(湍流,燃烧,辐射,多相) 基本方程特性具有: 1)封闭性:表明方程组可解; 2)非线性:不能直接求解,只能迭代求解; 3)偶和性:方程的复杂性不能用解析法求只能用数值法求; 4)形式相同:可建立求解的通用程序,基本规律的普遍性是通用 程序的物理基础,而基本方程形式的一致性则是通用程序的数学 基础。
缺点: ①数学模化的全面和准确性需要不断提高:
Ⅰ、物理问题的数学模型是否正确(回流问题还是边界层问题, 稳态还是非稳态),否则,数值算法的改进没有意义。
Ⅱ、所有物性数据要可靠,否则减少数值误差的努力毫无意义。 ②真实再现某些过程的代价也是极其昂贵的或不可能;(用于气象,
石油) ③有些迫不得已的简化是致命的或大大降低其价值; ④计算结果准确性仍需接受实验或精确解检验。(如对有代表性点
⑤具有模拟理想条件的能力: Ⅰ、如钢包吹氩时及镇静时,夹杂物上浮过程; Ⅱ、连铸二冷动态控制; Ⅲ、如实际有一个绝热表面,计算中很易做到; Ⅳ、创造极端条件(高温,高压,高空或无重力影响),而核电
站或化工厂泄露,水坝崩塌和城市大火等过程的实验则几乎无 法进行,计算程序中可改变几何、物理和化学因素及改变环境 条件以加快研究速度,拓宽研究参数变化范围。
流动与传热的数值计算
§1 绪论
1.1 引言 1、传热、传质与流体流动的重要性
工程设备(如结晶器,中间包,钢包及锅炉,高炉等) 内部流体流动及热交换过程,自然环境中的污染问题,暴 风雨雪,河流泛滥及着火过程中出现的热、质传递,流动 起着重要作用。 2、对过程估计和认识的必要性
过程表现为生产环节的一个阶段,其中存在着传热、传 质、流体流动及化学反应等现象,它是支配过程发展的最 基本的现象和子过程。
稳态问题→椭圆型方程问题→双向坐标 物理状态:具有一个单向坐标的流动→边界层型的流动问题 全部都是双向坐标的流动→回流流动问题
4)计算方法上的含义: 如可用一个单向坐标来规定一个给定的状态,那么就有可能大 大节省计算机的存储量与时间。
如: ① 一个不稳态的二维热传导问题:给定 t 时刻Τ分布求t+Δt 时刻
1)空间坐标一般是双向坐标,而时间坐标则总是单向坐标; 2)空间坐标也可能作为单向坐标;
如在一个坐标方向上有很强的单向流动,则信息的传递或影 响只能是从上游传至下游,即某点上的状态主要受其上游条 件的影响。(如强制对流过程) 3)数学及物理状态上的各自描述:
数学上微分方程中只要有非稳态项→抛物线方程问题→单向 坐标(无论有无空间上的);
4、预测的方法 1)实验研究
多数是在相似原理指导下进行的模拟实验,实验结 果是在一定实验条件下,忽略处于次要地位因素下 外推到全比例设备上去的,故降低了模型试验结果 的效能,且高温时测量困难及测量仪表也有误差。
2)数值计算方法 理论计算方法指的是用数学模型对过程或现象进行预测的方 法。数值计算方法利用大型计算机将微分方程→代数方程进 行求解.
2.1 计算流体动力学CFD(Computational Fluid Danamics)
CFD基本思想:把原来在空间与时间坐标中连续的物理量的场 (速度场,温度场,浓度场等),用一系列有限个离散点 (节点)上的值的集合来代替,通过一定的原则建立起这些 离散点上变量值之间关系的代数方程(称为离散方程),求 解所建立起来的代数方程以获得所求解变量的近似值。
优点: ①成本低:资金投入低(设备,测量手段),人力(不用实验
人员) ②资料信息全面完备:空间上每一位置的信息,每一时刻的计
算结果均可得到:Ф=f(x,y,z,t),并且可作为补充实验资料也 是很有价值的。 ③速度快:设计者可在几天时间内研究几十种不同方案。
④具有模拟真实条件的能力:不需缩小或冷态,计算机程序可模 拟任何大小尺寸,任何温度、速度、压力的过程。(如钢包吹 氩)
的Τ场,计算机内存只要供这两个Τ场即可。对于不同的Δt,可 反复使用同样的存储空间;计算上需同时处理的因变量只有一 个二维温度数组,与未来的温度值无关,节省计算时间。 ②一个二维边界层问题,只要给定上游沿垂直流线方向上因变量 值分布,即可顺次求出下游的,只需一维的计算机存储量。
§2 流动与传热数值模拟的基础知识
3)移动边界现象的模拟仿真 如一钢包内部钢水液面不断下降,铸造过程缩孔,结晶器液 面波动。
4)多相多组元的耦合仿真 流体中有固体颗粒,气泡。每相都需要求解方程,且液相与 气相、固相间还发生化学反应。
5)传热、传质、流体流动及燃料燃烧的耦合仿真
1.2 坐标的性质
1、 自变量:Ф=f(x,y,z,t),Ф=f(x,y,t), Ф=f(x,y,z) 2、坐标的选取
三维m×n×g,网格节点数多,如能减为二维,一维可使节 点数目大大减少 。 1)如连铸过程:板坯、方坯可按二维,薄板坯可按一维传热计 算; 2)圆管内轴对称流动:三维 x×y×z 可按二维 r×z 3、单向与双向坐标 定义:在一个坐标上的给定位置处的条件与其两侧的条件变化 有关,则为双向坐标,只与一侧有)
5、当前数值计算的发展趋势 1)坐标系的选取 正交坐标系:三个方向均是正交的。如飞机头部变化剧烈 部分在飞机表面上,如何划分网格。正交坐标系→非正交坐 标系。(二阶偏导由正交坐标系→非正交坐标系转换要增加 几十项) 2)非稳态现象的模拟仿真 如一钢包浇注过程中,内部钢水流动处于非稳态过程,数值 模拟计算得到,液面下降到200mm处不能下渣。