发一个新手学习的简单气泡模拟的帖子
用VOF模拟气泡上升过程,可能结果和欧拉模型差很多,这里不做比较了。前期只需划分一个长方形面网格就好,高度长度自己感觉,合适就行,上面那个边设置成压力出口。选用vof模型(记着打开Implicit Body Force),k-ε,非定长记住设置重力加速度为-9.81!设置材料,从材料库里添加water-liquid设置相,主相为water,第二相为air(相就是那个Phases)边界条件其实什么也不用设置然后直接全局初始化接着在菜单栏上选Adapt---Region,shape选择circle然后填入圆的圆心半径,点mark,这样就把气泡区域给标记出来了接着把标记出来的气泡区域填上空气,这一步就是patch,点patch然后把phase选成air,选中Vomule Fraction,再把出来的value填上1,在registers to patch中选中你刚才标记的区域,点击mark,这一步就完成了接着设置动画(这一步就不说了)然后时间步,迭代次数,就能开始计算了!!这次时间步选的大了,和上次计算还不一样,应该注意选小点,本来分裂开的两个气泡应该同时出去的。我也是个新手,希望共同学习进步。这个例子献给正在摸索的童鞋们!动画比较: 第一个动画步长为0.1s,第二个为0.3s,分析结果就不一样~ 没考虑表面张力吧 没有呀,考虑的话实在interaction里设置吧,这个一般设置多少呢? 我刚刚又用欧拉模型算了下,不太一样呀。 考虑张力就不那么容易破裂了 恩,的确是,不过最终还是破裂成两个了,不知道模拟的对不。刚刚又用恒定的表面张力计算了下,但是感觉差别不大! 除了要考虑表面张力之外,时间步长也不能太大,否则必破无疑。 在初始化的时候,如果考虑更精细一些,可以对计算域的压力按不同高度进行初始化。 随着与水面的距离的减小,外界压力越来越小,而气泡体积若过大的话,破裂也是有可能的。 表面张力系数必须设,一般设成0.6. 第二体积分数你用的是什么差分格式。geo-construction,还是施主受主。都会有一定影响
页:
[1]
2