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
- 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处不能下渣。