下载文档原格式
流体流动数值模拟流体流动现象普遍存在于⾃然界及多种⼯程领域中。
所有这些流动过程都遵循质量守恒、动量守恒、能量守恒和组分守恒等基本物理定律；⽽且流动若处于湍流状态，则该流动系统还要遵守附加的湍流输运⽅程。
本讲座将依据流体运动的特性阐述计算流体动⼒学的相关基础知识及任务；在流体运动所遵循的守恒定律及其数学描述的基础上，介绍数值求解这些基本⽅程的思想及其求解过程。
第⼀节计算流体动⼒学概述计算流体动⼒学（CFD）技术⽤于流体机械部流动分析及其性能预测，具有成本低，效率⾼，⽅便、快捷⽤时少等优点。
近年来随着计算流体⼒学和计算流体动⼒学及计算机技术的发展, CFD技术已成为解决各种流体运动和传热，以及场问题的强有⼒、有效的⼯具，⼴泛应⽤于⽔利、⽔电，航运，海洋，冶⾦，化⼯，建筑，环境，航空航天及流体机械与流体⼯程等科学领域。
利⽤数值计算模拟的⽅法对流体机械的部流动进⾏全三维整机流场模拟，进⽽进⾏性能预测的⽅法越来越⼴泛地被从事流体机械及产品性能取决于各种场特性的设计、科研等科技⼈员所使⽤；过去只有通过实验才能获得的某些结果或结论，现在完全可借助CFD模拟的⼿段来准确地获取。
这不仅既可以节省实验资源,还可以显⽰从实验中不能得到的许多场特性的细节信息。
⼀、什么是计算流体动⼒学计算流体动⼒学(Computational Fluid Dynamics，简称CFD)是通过计算机数值计算和图像显⽰，对包含流体流动和有热传导等相关物理现象的系统所做的分析。
CFD的基本思想可以归结为：把原来在时间域及空间域上连续的物理场（如速度场和压⼒场，以及热⼒场等），⽤⼀系列有限个离散点上变量值的集合来代替；并通过⼀定的原则和规律建⽴起关于这些离散点上的场变量之间关系，从⽽组成这些场变量之间关系的代数⽅程组；然后求解这种代数⽅程组，来获得这些场变量的近似值[1-3]；这就是流动的数值计算。
或者直观地说，通过数值计算中的各种离散⽅法，把描述连续流体运动的控制偏微分⽅程离散成代数⽅程组，由此建⽴该流动的数值模型；再根据问题的具体情况，设定边界条件和初始条件封闭⽅程组；然后通过计算机数值计算求解这种代数⽅程组，从⽽获得描述该流场场变量的某些运动参数的数值解。
cfd数值方法CFD(Computational Fluid Dynamics)数值方法，即计算流体力学数值方法，是通过利用数值计算方法对流体运动进行数值模拟，从而求解流体力学方程的一种方法。
CFD数值方法已经成为了流体力学分析中的重要分支，并且在航空、汽车、船舶、电子、建筑等领域得到了广泛的应用。
CFD数值方法的基本原理是将流体动力学方程组离散化，采用数值方法求解得到流场、温度场、压力场等物理量。
在CFD数值方法中，我们需要对流体运动的连续性、动量和能量守恒等方程进行求解。
这些方程是流体力学方程的基础，在CFD数值方法中有多种不同的求解方法。
其中，最常用的方法是有限体积法(Finite Volume Method，简称FVM)，这种方法将求解区域划分成若干个小体积，对每一小体积应用质量守恒、动量守恒和能量守恒方程进行求解。
在FVM方法中，需要对流体运动中的速度、压力、浓度等物理量进行离散化处理，并通过代数方程求解得到数值解。
除了FVM方法外，还有有限元法(Finite Element Method，简称FEM)、差分法(Finite Difference Method，简称FDM)等数值方法。
这些方法中，FEM方法的应用场景较广，可以对非结构化网格进行求解，其优点体现在对高级复杂结构的求解和可视化方面，但应用在液体/气体流体求解时，计算速度相对慢。
而FDM方法因为其求解速度快、实现简单等特点，在实际工程计算中应用较多。
总的来说，CFD数值方法在流体力学方面的应用发挥了重要作用，为工业生产与科学研究提供了有力支持。
但是，由于计算机性能限制，CFD在求解实际问题时也面临着许多挑战，尤其在复杂流动物理行为的求解中，还需要进一步发展数值技术，提高计算精度和效率。
基于CFD-PBE耦合的连续结晶过程数值模拟
潘立杰;刘伯潭;赵文立;盖晓龙
【期刊名称】《化学工程》
【年(卷),期】2024(52)1
【摘要】通过CFD(计算流体力学)与PBE(粒数衡算方程)耦合对连续结晶过程进行数值模拟。
采用QMOM(积分矩方法)求解PBE。
在考虑成核、生长、聚并和破碎的情况下,通过网格独立性验证确定计算网格数量。
在此基础上研究挡板设置、停留时间、搅拌转速以及搅拌桨类型等因素对CSD(晶体粒径分布)演变的影响,结果表明:挡板的设置增强了结晶器内湍流程度;高转速和低停留时间会获得平均粒径较小的颗粒。
对比3种搅拌桨对CSD演变的影响,轴流式的螺旋桨功耗小,适用于大颗粒晶体结晶过程,而径流式的圆盘涡轮桨适用于微粒结晶过程。
【总页数】7页(P41-47)
【作者】潘立杰;刘伯潭;赵文立;盖晓龙
【作者单位】天津科技大学化学工程与材料科学学院;天津大学化学工程联合国家重点实验室
【正文语种】中文
【中图分类】TQ026.5
【相关文献】
1.异型坯结晶器内钢水流动和凝固过程的耦合数值模拟
2.连铸结晶器内铸坯温度场和应力场耦合过程数值模拟
3.基于多相流耦合过程数值模拟的茶鲜叶离心式连续
脱水设备参数模拟与优化4.基于连续-非连续耦合方法的降雨滑坡数值模拟研究5.漏斗型结晶器内坯壳凝固过程热力耦合数值模拟
因版权原因，仅展示原文概要，查看原文内容请购买。
万方数据万方数据万方数据万方数据万方数据第９期张锴等：流化床内颗粒流体两相流的ＣＦＤ模拟时难以获得颗粒的真实堆积率，因此研究者们需要假设最大颗粒堆积率，如洪若瑜等［４９’５６巧７１采用ｏ．５５，Ｃｈｅｎ等№１４３取ｏ．６０，Ｌｅｔｔｉｅｒｉ等［４５］选Ｏ．６２。
３．１液固体系在Ｏ．５ｍ（高）×０．１ｍ（宽）的二维流化床考察了液（ＩＤｌ＝１０００ｋｇ・ｍ一，产ｌ一１．ｏ×１０－３Ｐａ・ｓ）固（佛＝３０００ｋｇ・ｍ～，或一２．５×１０－３ｍ）体系内网格尺度、时间步长和收敛判据对床层固含率分布特性的影响。
结果表明：（１）从整体来看网格数目和时间步长对床层固含率分布的影响不大，但是从局部放大图可以发现，当网格数目（１０×５０和１５×７５）较少时，平衡时垂直方向上的固含率出现振荡，且１０×５０网格的振荡幅度大于１５×７５的网格，而网格数目（２０×１００和３０×１５０）较多时，床层固含率趋于均匀分布特征；（２）通过对０．０１、Ｏ．００５、０．００１、Ｏ．Ｏ００５ｓ和Ｏ．ｏ００１ｓ时间步长的模拟表明，ｏ．００１ｓ时间步长给出了更适宜的模拟结果；（３）收敛判据取１０一、１０－６和１０＿。
，所得模拟结果几乎完全一致，详细结果见文献［５８］。
３．２气固体系首先采用摄像法考察了图２所示中心孔口为Ｏ．０１０ｍ的２．Ｏｍ（高）×Ｏ．３ｍ（宽）拟二维流化床内射流形成及发展过程、射流穿透深度和射流频率。
实验以常温和常压下的空气为流化介质，ＧｅｌｄａｒｔＢ类物料的玻璃珠（佛＝２５５０ｋｇ・ｍ一，矾一２５０～３００肛ｍ，“ｍｆ一０．０７ｍ・ｓ－１）为固体。
通过对射流气速为７．０７ｍ・ｓ。
１的１２００张图像进图２实验装置流程示意图Ｆｉｇ．２｛ｋｈｅｍａｔｉｃｄｉａｇｒａｍｏｆｅｘｐｅｒｉｍｅｎｔａｌａｐｐａｒａｔｕｓ行逐帧分析，发现当时间为ｏ．０２５ｓ时射流已经形成并开始逐渐长大，到ｏ．１５０ｓ时，该射流在分布器上方脱落形成气泡，并有新的射流产生。
