关于twoPhaseEulerFoam的问题 ？

[liliang2009223]

本人在利用twoPhaseEulerFoam求救器求解算例bubbleColumn时，将网格文件blockMeshDict中的计算区域改成了五度角的二维轴对称（类似于nozzleFlow2D中的网格），然后把速度入口改成了X方向速度，重力文件改为沿X负方向，其他边界条件没有动，结果运算不能进行，而且alpha出现严重的越界!，想问下到底是什么原因！谢谢！

[李东岳]

Hi Liliang,

请问是哪个版本的OpenFOAM？能否发网格看一下？以及问题所在的log输出？

[liliang2009223]

不好意思 刚看到！ 我用的of版本为2.4. log文件图片和网格图片如下!

[liliang2009223]

@李东岳 版主可否帮忙下什么原因？ 谢谢！

[李东岳]

@liliang2009223

你好，

网格不是很清楚。太密了。你先稀疏一下，网格线符合物理么？
log错误提示exp出错。有可能是你的exp函数里面是一个很大的数。你自己添加代码了么？

[liliang2009223]

@李东岳 没有动代码！exp是什么函数？

[liliang2009223]

@李东岳 checkMesh网格检查没有问题

[李东岳]

@liliang2009223 指数函数，$e^x$. :)

[李东岳]

@liliang2009223 twoPhaseEulerFoam对网格正交要求很高，checkMesh没问题在计算的时候也有可能有问题。哈。可以发网格出来看看正交性和skewness

[liliang2009223]

@李东岳 把网格检查的输出发上来可以吗？因为网格太大没法全局视图

[liliang2009223]

@李东岳 最大的skewness=0.333308

[liliang2009223]

@李东岳

[李东岳]

@liliang2009223 说：

改成了五度角的二维轴对称（类似于nozzleFlow2D中的网格），然后把速度入口改成了X方向速度，重力文件改为沿X负方向，其他边界条件没有动，结果运算不能进行，而且alpha出现严重的越界!，

算例文件多大，压缩一下，上传。我看下。把没必要的文件删除。太大的话发我邮箱，我调试后给你在这里回复。

[liliang2009223]

文件挺大 我发你邮箱吧 谢谢！

[liliang2009223]

@李东岳 李老师文件已经发你info@dyfluid.com邮箱 谢谢！

[liliang2009223]

@李东岳 穿你的文件我已经 foamCleanTutorial ， 边界条件就把原来的速度，温度和压力的值改了下！

[李东岳]

@liliang2009223 你的入口速度是350m/s?

[liliang2009223]

@李东岳 对 入口速度很大

[liliang2009223]

@李东岳 我那个速度是液相的入口速度 内部区域是气相！可能给你的文件没有改过来！不好意思！

[liliang2009223]

@李东岳 李老师不好意思，我明天再重新传你一下原文件吧！谢谢！