上海安世亚太汇智科技股份有限公司欢迎您,咨询热线:021—58403100
Resources 资源中心
Hot Resources / 热门资源 More+
2020 - 11 - 25
作者:王晨希  上海安世亚太高频电磁高级技术专家本文共计2216字,阅读时间预计7分钟编者按作者分析了美国空军一号抵御核爆的秘密,包括如何抵御冲击波、核辐射以及电磁脉冲。在抵御电磁脉冲上,美国空军一号用了三层保护机制,可见对飞机这种精密电子飞行器杀伤力更大的当属电磁脉冲。大名鼎鼎的美国空军一号,美国总统的专属座机,被誉为世界上最安全的飞机。这架不同于任何传统客机的总统专机,拥有世界上最顶...
2020 - 10 - 21
作者:陈康 上海交通大学本文共计3111字,阅读时间预计10分钟编者按作者介绍了壁膜现象、欧拉壁膜模型与部分浸润效应的内在机理,并依此建模进行案例测试分析,以此探索不同壁膜接触角下与质量通量下壁膜形态演变规律的差异,希望可以启发更多的思考探索。介绍壁膜广泛存在于日常生活中,并且在工业生产中也扮演重要角色,如汽车车窗的除雾、冰箱设计、食品冷藏技术中冷凝壁膜的杀菌等应用。壁膜是由液滴撞击到固体壁面上形...
2020 - 10 - 13
作者:吴华春 上海安世亚太结构应用工程师本文共计1559字,阅读时间预计5分钟编者按在ANSYS LS-DYNA中,经常遇到Sliding Energy异常,作者在本文中分析了计算模型常出现的2种情况——初始穿透导致的Sliding Energy异常和由于SEGMENT接触导致的Sliding Energy异常,并给出了对应的解决方法。问题提出负的滑移能存在2种情况:■ 第一种情况:GLSTAT中...
2020 - 09 - 26
翻译:赵亚  上海安世亚太流体应用工程师本文共计1320字,阅读时间预计5分钟编者按上周我们谈到多相流的分类及仿真方法,如果大家记忆有些模糊的话,点此链接来回忆上篇内容吧!接下来作者会进行案例分析和行业解决方案解析,展示ANSYS在多相流仿真领域的运用。多相流仿真展望即使当今多相流仿真功能已经十分强大,但是工程师仍在努力推进多相流仿真的边界。他们将融入并模拟更多的物理模型:物化反应、运动...
2020 - 09 - 23
翻译:赵亚  上海安世亚太流体应用工程师本文共计1785字,阅读时间预计6分钟编者按作者分析了多相流的分类、仿真方法及未来发展方向,深入剖析了真实的产品案例和行业解决方案,并提出ANSYS仿真工具可以快速、高效地解决产品设计在多相流仿真中遇到的挑战,为企业的产品设计和生产节约成本。无论是设计高超声速运载工具的除冰系统、进行血液酶测试,还是输送和熔化稀有金属粉末化合物以进行增材制造,或是为...
2020 - 09 - 16
作者:陈诗佳  上海交通大学-巴黎高科卓越工程师学院谢文丽  上海安世亚太流体应用工程师本文共计1755字,阅读时间预计6分钟编者按作者使用Flownex对垃圾处理中心通风系统的局部管线进行模拟,对敞开式处理车间和封闭罐装设备的环境通风系统做出整体评估,并提出改进建议,同时对管线上的风机选型提出明确要求。研究目的随着国内环境三废处理规范要求越来越明细化,各种环保处理工...

服务热线 / Service Hotline

021-58403100

中国 • 上海 • 地址

上海市静安区石门一路288号兴业太古汇一座1201~1204室

热线电话:021-58403100

邮箱地址:sh.marketing@peraglobal.com

客户微信:PeraHZtech

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真
单份价格:
市场价:

份数: - + 件(库存件)

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

介绍

作者:陈康 上海交通大学

本文共计3111字,阅读时间预计10分钟


编者按

作者介绍了壁膜现象、欧拉壁膜模型与部分浸润效应的内在机理,并依此建模进行案例测试分析,以此探索不同壁膜接触角下与质量通量下壁膜形态演变规律的差异,希望可以启发更多的思考探索。

介绍

壁膜广泛存在于日常生活中,并且在工业生产中也扮演重要角色,如汽车车窗的除雾、冰箱设计、食品冷藏技术中冷凝壁膜的杀菌等应用。

壁膜是由液滴撞击到固体壁面上形成的。液滴撞击壁面后的情况有以下四种:

■ 附着(stick):液滴以很小的动能撞击壁面并近似保持球形;

■ 反弹(rebound):液滴改变速度,相对完整地离开壁面;

■ 铺展(spread):液滴以中等动能撞击壁面并铺展为壁膜;

■ 飞溅(splash):液滴的一部分留在壁膜中,另一部分以一些更小尺寸的小液滴离开壁面。


薄膜假设:壁膜厚度远小于壁面的曲率半径,壁膜在厚度上的属性是一致的,且壁膜流动平行于壁面。

基于以上假设下的壁膜模型分为基于场的欧拉壁膜模型和基于粒子的拉格朗日壁膜模型。本文采用欧拉壁膜模型。


基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图1. 壁膜模型示意图

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图2. 液滴与壁面相互作用决策图


欧拉壁膜的质量、动量、能量守恒

质量守恒

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

等式左边:非稳态项和对流项;

等式右边:单位面积下的质量源项,如液滴收集、壁膜分离、壁膜脱落、相变等行为下,需更新壁膜质量源项。


动量守恒

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

等式左边:非稳态项和对流项

等式右边:

→ 第一项:扩散项:其中,压力=气流压力+垂直于壁膜的重力分量+表面张力

→ 第二项:重力源项:平行于壁膜的重力分量

→ 第三&四项:净粘性切应力源项:气流-壁膜之间&壁面与壁膜之间

→ 第五项:动量源项:液滴的收集与分离

→ 第六项:表面力源项:壁膜的表面张力、壁膜与壁面间的接触角


能量守恒

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

等式左边:非稳态项和对流项

等式右边:

→ 第一项:气流-壁膜间与壁膜-壁面间的净热通量

→ 第二项:  :液体从气流中撞击到壁面的能量源项

→ 第三项:相变能量源项: 是蒸发/冷凝的质量变化率,   是相变潜热值


控制方程的应用与求解

■ 质量守恒和动量守恒:欧拉壁膜模型的核心基础

■ 能量守恒:涉及壁膜的传热需要启用

■ 由于欧拉壁膜模型适用的液膜很薄,故采用润滑近似(平行流动),因此这些控制方程是在平行于壁面的局部当地坐标中求解的。


壁膜速度的获取

壁膜速度可以假设为两个组成:外部流动产生的切应力驱动壁膜速度、由于重力产生的重力驱动壁膜速度。

切应力驱动-速度线性变化

由气液交界面处外部气体流动侧,切应力与壁膜侧切应力平衡,得出外部流动产生的切应力驱动壁膜速度:

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图3. 切应力驱动下壁膜速度变化

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真
基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真


重力驱动-速度抛物线变化

由于重力与壁面切应力的平衡,壁膜速度呈现抛物线变化,壁膜速度在气液交界面的梯度为0,得出重力驱动壁膜速度:

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图4. 重力驱动驱动下壁膜速度变化

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真
基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真


不同应用场合下的壁膜速度

■ 在有外部流体流动的重力场下,如湿空气在冷壁面的冷凝结露问题中,壁膜速度同时受气液界面的切应力与壁膜重力驱动。

■ 在无外部流体流动的重力场下,如壁膜生长流动问题中,壁膜速度完全受重力驱动。

■ 在有外部流动的无重力场下,如航天器燃料储液箱问题中,壁膜速度完全受气液界面切应力驱动。


部分浸润效应

壁膜的部分润湿效应描述了诸如溪流、干点和膜破裂等行为,部分浸润现象广泛存在于工业生产与日常生活中,探索其内在作用机理对完善壁膜建模测试过程尤为重要。

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图5. 壁膜流动形态演变示意图


干区和湿润区是由壁面和壁膜的接触线划分的。

■ 壁膜湿润面积分数方程:

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

h:壁膜厚度

hc:临界壁膜厚度(fluent中1.e-7 mm)

■ 部分浸润效应形成机理:液膜与壁面间接触线的表面力限制壁膜扩散流动

■ 表面力的影响因素:壁膜表面张力σ 、壁膜流体/壁面的接触角 θw

■ 表面力方程:

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

σ:表面张力

β:表征实验与模型的差异的经验系数,范围:0-10

w:壁膜湿润面积分数

θw:接触角


测试案例

如图6所示,建立三维模型,尺寸为8x100x200mm。其中流体域工质为空气,进出口速度设为0(此三维流动模型下只为欧拉壁膜构建,下文不再赘述流体流动的设置)。侧面设为欧拉壁膜模型,分为模拟壁膜进口的壁膜入口边界与模拟壁膜形态演变的壁膜流动壁面。

测试内容

建立壁膜流动模型,研究在重力场中不同流量下,壁膜的流动形态。

 

测试目的

探讨部分浸润效应下壁膜流动的内在机理,分析欧拉壁膜模型与部分浸润效应的条件设置思路,通过对比工况,研究不同输入参数(接触角、质量流量)对壁膜流态模拟结果的影响

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图6. 壁膜流动模型


测试流程

网格设置

采用六面体网格,网格尺寸为2x2x2mm,其中壁膜壁面中第一层网格设置为壁膜入口边界(2x80mm),其余部分设为壁膜流动壁面。

物理模型

激活能量方程、层流模型、欧拉壁膜模型。其中:


■ 欧拉壁膜模型的模型与材料设置中,需激活求解壁膜动量选项,并在动量选项中打开重力选项与压力梯度选项,在此基础上才能打开扩展项和表面张力项;在壁膜材料中选取水作为工质,设置其表面张力为0.07194。


■ 欧拉壁膜模型的时间离散与空间离散选项中,为提高计算精度,壁膜的时间、连续与动量离散均选择二阶离散。壁膜最大厚度与壁面第一层网格高度相关,设为2mm。


■ 打开质量与动量方程耦合求解和弯曲平滑,默认情况下,先计算壁膜的质量方程,再计算壁膜的动量方程。而当壁膜特征量——壁膜高度与速度是耦合时(如波浪表面的自由下落壁膜),需采用同时计算质量与动量方程的耦合求解。耦合求解主要应用于弯曲壁膜和壁膜表面张力计算的场合。

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图7. 壁膜流动模型网格

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图8. 欧拉壁膜模型子模型与材料设置

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图9. 欧拉壁膜模型时间离散与空间离散选项

壁膜壁面条件设置

■ 壁膜边界入口设置

在壁膜壁面上方处设置长宽为2x80mm区域为壁膜入口边界条件,给定壁膜质量通量。对比工况设置质量通量分别为5、10、20、50、100、200kg/(m2·s),对比研究不同流量下,壁膜流动形态演变规律。

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图10. 欧拉壁膜模型壁膜边界入口


■ 壁膜流动壁面设置

◆ 将壁膜壁面除壁膜入口边界外区域设为壁膜流动壁面。

◆ 壁膜部分浸润效应的条件设置:激活初始条件下壁膜接触角选项,设置壁膜与壁面间的平均接触角(对比工况分别设置为75°和15°,分别表征疏水与亲水表面)、相对标准误差(人为设定为13.5%)与表征实验与模型差异的经验系数(人为设定为1)。

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

▲ 图11. 欧拉壁膜模型壁膜流动壁面


测试结果

不同质量通量下壁膜的形态演变

从下列图中可以明显看出不同质量通量下壁膜形态演变的差异:


■ 给定壁膜/壁面接触角为75°下,质量通量为5kg/(m2·s)时,壁膜的表面张力占优,壁膜倾向于分离收缩,壁膜由上方入口受重力驱动沿着壁面往下流动,在流动初始便出现壁膜分离,形成溪流现象;


■ 随着质量通量的增加,当质量通量为10kg/(m2·s),壁膜厚度与速度明显增加,流动规律仍然呈现初始壁膜分离并形成溪流现象,但在流动下半段开始出现溪流合并现象;


■ 当质量通量增加为20kg/(m2·s)时,壁膜流动现象变得更为复杂,随着时间的推移,连续壁膜的分离点明显下移,并且溪流流动过程中出现溪流合并再分离的复杂流动现象;


■ 当质量通量增加为50kg/(m2·s)时,壁膜的质量通量占优,壁膜倾向于铺展流动,壁膜流动开始出现转折性变化,随着时间推移,壁膜流动整体呈现连续流动现象。


■ 当质量通量增加为100、200kg/(m2·s)时,壁膜连续流动的规律愈发明显,并且在同一水平线上,在大接触角下的疏水性壁面上的壁膜厚度从左到右先变小再变大,即在左右俩侧边缘处的壁膜最厚。此现象的机理是在壁膜质量通量占优的情况下,壁膜整体呈现连续流动的规律,但由于壁面的疏水性,壁膜仍然存在收缩的趋势,故在边缘处的壁膜最厚并且流动速度最大。


不同接触角下壁膜的形态差异

如图14所示,给定质量通量为:20 kg/(m2·s),对比研究时间为300ms,壁膜/壁面接触角分别为75°和15°,图为壁膜流动的形态差异。


壁膜在疏水性和亲水性表面上的流动分别呈现壁膜分离溪流流动和连续流动的规律。此外,联系上文,对比壁膜在疏水性表面和亲水性表面上的连续流动现象,发现在同一水平线上壁膜在亲水性表面上呈现俩侧厚度小、中间厚度大的规律。

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

壁膜质量流:20 kg/(m2·s)

接触角:75°

时间:300 ms

基于部分浸润效应的欧拉壁膜流动形态演变模拟仿真

壁膜质量流:20 kg/(m2·s)

接触角:75°

时间:300 ms

▲ 图14. 同接触角下壁膜的形态差异


结语

本文简单介绍了欧拉壁膜模型与部分浸润效应的基本原理与建模思路。

此壁膜模型的仿真模拟可应用于提高相变换热器中冷凝器的液体收集效率(有效降低压缩机的液击危害)、填料塔中液体分布器的优化设计(提高液体分布均匀性与填料传质效率)等应用场景。

壁膜模型还包含更多的子模型与相关应用,如可应用于喷淋冷却、灭火场景的拉格朗日壁膜模型建模、可应用于湿空气在冷壁面冷凝与再烘干过程的组分输运与欧拉壁膜模型的耦合建模。


仅借此文抛砖引玉,希望启发更多的思考探索,让模拟仿真工作来源于生活生产,服务于生活生产。


*本文版权归上海安世亚太所有,

如需转载,请与我们联系。

视频
文件下载
推荐 / Products
2020 - 11 - 25
作者:王晨希  上海安世亚太高频电磁高级技术专家本文共计2216字,阅读时间预计7分钟编者按作者分析了美国空军一号抵御核爆的秘密,包括如何抵御冲击波、核辐射以及电磁脉冲。在抵御电磁脉冲上,美国空军一号用了三层保护机制,可见对飞机这种精密电子飞行器杀伤力更大的当属电磁脉冲。大名鼎鼎的美国空军一号,美国总统的专属座机,被誉为世界上最安全的飞机。这架不同于任何传统客机的总统专机,拥有世界上最顶尖的防御系统,包括防核弹冲击系统、反导系统、诱饵系统,内部设施包括总统生活区、飞行指挥中心以及大名鼎鼎的The nuclear football核弹密码箱等。▲ 空军一号空军一号有两架,飞行呼号都是Special Air Mission 28000和29000,当美国任何一个空中管制员听到这个呼号后都明白美国总统就在上空。空军一号曾一度是国家机密,用以在平常以及核战争时期充当美军的最高指挥中心。▲ 空军一号内部那么作为世界上安保最全面的飞机,面对末日核打击,该如何保护总统及机上人员安全?下面为大家揭开空军一号抵御核爆的秘密。▲ 核弹密码箱首先需要说明的是,无论空军一号的防护多么完备,发生核爆时都不可能在原爆点存活。因为现今人类发明的最强隔热物质能够隔绝的温度不过几千摄氏度,而核爆原爆点几公里内的温度可以达到上千万摄氏度,热核武器甚至可以达到上亿摄氏度,所有的物质都会瞬间气化。空军一号对核爆的防护主要是针对的是冲击波、核辐射以及电磁脉冲。 抵御冲击波今年8月份,在黎巴嫩首都贝鲁特发生了巨大爆炸,强大的冲击波造成了多栋房屋受损,玻璃震碎,港口完全被炸毁,距爆炸地点10多公里的总统府建筑也因此受损。与贝鲁特硝酸铵爆炸产生的冲击波相比,核爆产生的冲击波威力要强大得多,覆盖范围也要大非常多。▲ 贝鲁特爆炸冲击波▲ 核冲击波 为了抵抗...
2020 - 10 - 21
作者:陈康 上海交通大学本文共计3111字,阅读时间预计10分钟编者按作者介绍了壁膜现象、欧拉壁膜模型与部分浸润效应的内在机理,并依此建模进行案例测试分析,以此探索不同壁膜接触角下与质量通量下壁膜形态演变规律的差异,希望可以启发更多的思考探索。介绍壁膜广泛存在于日常生活中,并且在工业生产中也扮演重要角色,如汽车车窗的除雾、冰箱设计、食品冷藏技术中冷凝壁膜的杀菌等应用。壁膜是由液滴撞击到固体壁面上形成的。液滴撞击壁面后的情况有以下四种:■ 附着(stick):液滴以很小的动能撞击壁面并近似保持球形;■ 反弹(rebound):液滴改变速度,相对完整地离开壁面;■ 铺展(spread):液滴以中等动能撞击壁面并铺展为壁膜;■ 飞溅(splash):液滴的一部分留在壁膜中,另一部分以一些更小尺寸的小液滴离开壁面。薄膜假设:壁膜厚度远小于壁面的曲率半径,壁膜在厚度上的属性是一致的,且壁膜流动平行于壁面。基于以上假设下的壁膜模型分为基于场的欧拉壁膜模型和基于粒子的拉格朗日壁膜模型。本文采用欧拉壁膜模型。▲ 图1. 壁膜模型示意图▲ 图2. 液滴与壁面相互作用决策图欧拉壁膜的质量、动量、能量守恒质量守恒等式左边:非稳态项和对流项;等式右边:单位面积下的质量源项,如液滴收集、壁膜分离、壁膜脱落、相变等行为下,需更新壁膜质量源项。动量守恒等式左边:非稳态项和对流项等式右边:→ 第一项:扩散项:其中,压力=气流压力+垂直于壁膜的重力分量+表面张力→ 第二项:重力源项:平行于壁膜的重力分量→ 第三&四项:净粘性切应力源项:气流-壁膜之间&壁面与壁膜之间→ 第五项:动量源项:液滴的收集与分离→ 第六项:表面力源项:壁膜的表面张力、壁膜与壁面间的接触角能量守恒等式左边:非稳态项和对流项等式右边:→ 第一项:气流-壁膜间与壁...
2020 - 10 - 13
作者:吴华春 上海安世亚太结构应用工程师本文共计1559字,阅读时间预计5分钟编者按在ANSYS LS-DYNA中,经常遇到Sliding Energy异常,作者在本文中分析了计算模型常出现的2种情况——初始穿透导致的Sliding Energy异常和由于SEGMENT接触导致的Sliding Energy异常,并给出了对应的解决方法。问题提出负的滑移能存在2种情况:■ 第一种情况:GLSTAT中的滑移能在初始状态就出现负值,如图-1所示。▲ 图-1 第一种能量异常图■ 第二种情况:CAE工程师在利用ANSYS LS-DYNA进行计算时,很多时候会碰到图-2所示的状况,严重的能量不平衡。内能、滑移能在计算前期表现还算正常,但是到计算后期,内能、动能、滑移能严重偏离,导致整个模型能量异常,计算结果不可信。▲ 图-2 第二种能量异常图解决方法第一种情况的解决方法(初始穿透)在仿真模型中出现初始穿透,导致初始的滑移能为负值。仿真穿透模型如图-3所示。产生机理:LS-DYNA在进行接触搜索时,发现存在网格穿透问题,需要把穿透部分的节点移到不穿透的位置,这可以在LS-DYNA软件中设置。如图-4所示,当把初始穿透清除后,其滑移能和整个模型能量表现正常(如图-5所示)。▲图-3 存在初始穿透模型▲ 图-4 清除初始穿透模型▲ 图-5 第一种能量正常图初始穿透解决方法有以下几种① 手动调整单元节点  在进行模型网格划分和设置零件厚度时必须确认是否存在干涉,即必须考虑壳单元的接触厚度。如果发生穿透,可以从计算出的message或d3hsp文件中对关键字“initial penetrations”进行搜索,找到相关单元,然后调整单元节点,消除初始穿透。在一些专门的前处理软件中,例如ANSA,HYPERMESH、SpaceClaim在提交计算前对模型进行...
地址:中国·上海·浦东新区平家桥路36号晶耀前滩5号楼901室
电话:+86 0755-2955 6666
传真:+86 0755-2788 8009
邮编:330520
Copyright ©2005 - 2013 上海安世亚太汇智科技股份有限公司
犀牛云提供企业云服务