上海安世亚太汇智科技股份有限公司欢迎您,咨询热线:021—58403100
Resources 资源中心
Hot Resources / 热门资源 More+
2021 - 01 - 20
作者:李桂花 上海安世亚太结构应用工程师本文共计1199字,阅读时间预计4分钟编者按ICT技术已经是现代电子企业常用的测试技术,本文中作者使用ANSYS Sherlock软件进行ICT测试,评估由于过应力导致的风险组件,并根据结果提出改进措施。随着微电子技术和电子组装技术的发展,PCBA(Printed Circuit Board Assembly,印刷电路板组件)布局和组装呈现出复杂化...
2021 - 01 - 13
作者:高征宇 上海安世亚太流体应用工程师本文共计1435字,阅读时间预计5分钟编者按作者由浅入深,通过《数码宝贝》引申出数字孪生的发展历程与生活中的实际运用,分析了数字孪生运用中的重点及难点,最后通过案例展示了数字孪生的使用方式。或许大家有发现数字孪生的概念正慢慢地变得火热起来。不仅变成了众多行业、企业,各种大中小型研讨会追捧的话题,也成为了大家茶余饭后的谈资。数字孪生可以说1970年Apollo...
2021 - 01 - 07
作者:李超峰 上海安世亚太流体应用工程师文章发布:上海安世亚太官方订阅号(搜索:PeraShanghai)联系我们:021-58403100本文共计1371字,阅读时间预计4分钟编者按作者通过一个优化模具通道的案例进行分析和讲解,拆分了该案例的优化步骤,让我们清楚地了解到ANSYS POLYFLOW在挤出工艺中的重要作用。背景挤出成型是聚合物、玻璃及食品加工连续性生产工艺的重要方法。由于聚合物在拉...
2020 - 12 - 16
作者:刘琼  上海安世亚太光学应用工程师文章首发:上海安世亚太官方订阅号(搜索:PeraShanghai)联系我们:021-58403100本文共计1609字,阅读时间预计5分钟编者按《赛博朋克 2077》作为最近大热的一款游戏,因其高还原度场景、酷炫的光影效果和未来感而令无数玩家心驰神往。那么这款游戏到底为什么比其他3D游戏画面更逼真呢?作者在文中对其进行了分析,让你更懂《赛博朋克 2...
2020 - 12 - 09
编译:吴华春 上海安世亚太结构应用工程师本文共计1816字,阅读时间预计6分钟编者按本篇内容原文来自LSTC公司,作者阐述了如何使用LS-PREPOST在LS-DYNA关键字模型中查找并移除交叉边缘和初始穿透。什么是穿透?交叉边缘/ Crossed edges单元的边缘穿过另一个shell单元中间表面,或表面穿过solid单元,在LS-DYNA中不是经典意义上的穿透。但是,在LS-DYN...
2020 - 11 - 25
作者:王晨希  上海安世亚太高频电磁高级技术专家本文共计2216字,阅读时间预计7分钟编者按作者分析了美国空军一号抵御核爆的秘密,包括如何抵御冲击波、核辐射以及电磁脉冲。在抵御电磁脉冲上,美国空军一号用了三层保护机制,可见对飞机这种精密电子飞行器杀伤力更大的当属电磁脉冲。大名鼎鼎的美国空军一号,美国总统的专属座机,被誉为世界上最安全的飞机。这架不同于任何传统客机的总统专机,拥有世界上最顶...

服务热线 / 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
2021 - 01 - 20
作者:李桂花 上海安世亚太结构应用工程师本文共计1199字,阅读时间预计4分钟编者按ICT技术已经是现代电子企业常用的测试技术,本文中作者使用ANSYS Sherlock软件进行ICT测试,评估由于过应力导致的风险组件,并根据结果提出改进措施。随着微电子技术和电子组装技术的发展,PCBA(Printed Circuit Board Assembly,印刷电路板组件)布局和组装呈现出复杂化的趋势。PCB中的器件尺寸微型化、引脚间距紧凑化都使得 PCBA在极大程度上增大了产品失效的可能性。ICT技术ICT( In-Circuit Test),即在线测试技术,作为电子产品印制板组装加工中的重要测试手段,已为大多数现代化电子企业所采用。它使用一系列测试探针和测试夹具,在一个电路板的一面或两面来测试制造过程中的电气连接。每个测试探针施加一个力在一个特定的电路板位置,称为测试点,由设计确定。在测试过程中,所有这些测试点力的联合作用使电路板产生弯曲,如果应力值足够高,焊点可能会失效,PCB板上的元件也将产生超过允许范围的应力。图1  ICT测试仪图2 电路板常见的失效模式:陶瓷电容器的弯曲开裂、垫坑为了提前避免失误和发现问题,在生产过程中需要对ICT步骤进行应力应变测量并监控,以确定产品应力应变处于允许范围内。如果测量值超过了电路板允许的应力应变水平的最大值,将对电路板进行重新布局设计或者调整夹具设计,或者按要求改变流程,使得应力应变数值回到允许范围之内。ANSYS SherlockANSYS Sherlock软件是目前唯一一款基于可靠性物理的电子设计软件,能够在早期设计阶段为电子硬件的元器件、板件和系统级提供快速准确的寿命预测,其提供的ICT分析模块基于有限元算法可以快速对测试点和夹具施加的机械应力进行建模,并根据预测的电路板部件的应力对其进行打分,列出测...
2021 - 01 - 13
作者:高征宇 上海安世亚太流体应用工程师本文共计1435字,阅读时间预计5分钟编者按作者由浅入深,通过《数码宝贝》引申出数字孪生的发展历程与生活中的实际运用,分析了数字孪生运用中的重点及难点,最后通过案例展示了数字孪生的使用方式。或许大家有发现数字孪生的概念正慢慢地变得火热起来。不仅变成了众多行业、企业,各种大中小型研讨会追捧的话题,也成为了大家茶余饭后的谈资。数字孪生可以说1970年Apollo 13事故中拯救了三位宇航员的模拟器是第一个数字孪生的应用。而概念上的起源是2002年Michael Grieves在密西根大学做演讲时在PLM概念模型中首次提出的现实与虚拟空间的概念,之后在其著作中被称之为镜像空间模型(Mirror Space Model),演化至今成为了数字孪生(Digital Twin)。数字孪生是一种超越现实的概念,如今可以被视为一个或多个重要且彼此依赖的装备系统间的数字映射系统。数字孪生是某个特定实体产品/系统的数字表示,能够更深入地洞察该产品的状态、性能和行为;能够实现系统设计、优化、预见性维护、优化工业产品管理。这种方式充分利用了物理模型、传感器更新、运行历史等数据,集成为一个多学科、多物理量、多尺度、多概率的仿真过程,在虚拟空间中完成映射,从而反映相对应的实体装备的全生命周期过程。举个简单点的例子:不知道大家还是否记得《数码宝贝(Digital Monster)》 ,数码暴龙机与数码暴龙的结合就是比较典型的数字孪生早教片了。▲暴龙实体▲暴龙数字体这个听起来可能是养成 + 战斗类游戏,而反应在工业上严谨度会更高。类似的是:你有关心的重要设备(Monster)、各种不同的工况、不同的环境以及他们之间的影响(Buff),我们通过监控虚拟设备的运行情况,结合实际工况/数据的实时采集(HP/MP),能及时预测、发现并解决问题,使得系统良好运转(Fi...
2021 - 01 - 07
作者:李超峰 上海安世亚太流体应用工程师文章发布:上海安世亚太官方订阅号(搜索:PeraShanghai)联系我们:021-58403100本文共计1371字,阅读时间预计4分钟编者按作者通过一个优化模具通道的案例进行分析和讲解,拆分了该案例的优化步骤,让我们清楚地了解到ANSYS POLYFLOW在挤出工艺中的重要作用。背景挤出成型是聚合物、玻璃及食品加工连续性生产工艺的重要方法。由于聚合物在拉伸和挤压过程中表现出复杂的粘弹性行为,因此熔体在模具中和自由端的流动行为非常复杂。目前主要采用实验方法和数学解析的方法来进行研究,前者由于需要反复模具试制、多次试模才能得到合格的产品,且一类经验在其他类型产品中复制性差,造成生成时间周期长,难以满足激烈的市场竞争需求;对于数学解析的方法,随着流变学的发展和计算机技术提高,使这种方法越来越成熟,ANSYS POLYFLOW经过多年的发展,目前已包括丰富的求解功能:可以求解的网格类型多包括三角形、四边形、四面体、六面体等,网格变形工具包括网格叠加技术、网格自适应、网格重构,适应求解网格大变形问题;对于前处理软件有很好的接口可以读入ANSYS meshing、ICEM、GAMBIT网格工具生成的网格;丰富的求解模型包括牛顿流体、广义牛顿流体、粘弹性流体等多种流变模型。因此,POLYFLOW特别适合应用于聚合物成型加工、玻璃工业及食品行业中的材料流动模拟。在塑料产品加工行业应用包括挤出成型、吹塑成型、拉丝、涂覆工艺的流动、传热问题。POLYFLOW在挤出成型中的优化应用在挤出成型中,挤出产品的形状是挤出口模设计的一个重要目标,而挤出胀大现象给具有形状和尺寸要求的模具设计带来了一定难度。另外,熔体流出口模截面速度的差异,同样会造成产品内应力增大,造成产品变形严重。对此,我们可以通过优化模具内部通道,改善出口截面流速的差异,...
地址:中国·上海·浦东新区平家桥路36号晶耀前滩5号楼901室
电话:+86 0755-2955 6666
传真:+86 0755-2788 8009
邮编:330520
Copyright ©2005 - 2013 上海安世亚太汇智科技股份有限公司
犀牛云提供企业云服务