CFD中文网

@李东岳 李老师，请问如果求解非守恒型的动量方程，那么对于其中一相为0的情况，是不是就直接求解就好了，不用再做别的处理了，如让相为0的相的速度等于按照模型去求出的速度。

topo.png
我需要用ICEM画multiblock网格，几何类似于上图的模型，我知道三角形Y-block和圆形O-block，那请问这两个如何结合起来呢？在这种结构中它的edge并不共用相同的vertices
第二个问题是icem的multiblock的网格如何光顺？我尝试使用了pre-mesh smooth，但是这个功能很慢，容易卡死。我只知道可以用edit edge调整spline来改变edge形状，如果一个一个edge调是不是有点笨？

@卡洛 请问你是用的哪种方法？

有个方向，叫混沌之类的。详细可以了解下蝴蝶效应和非线性的一些内容。

solidworks导入icem时候有些地方确实有点问题，我现在先用sw把大框架画出来，自己在icem里面再编辑一下小框架，最后再切分一下，有点麻烦。你现在已经解决把，有啥好办法吗

这个潘工的联系方式给我呀，500块买不了上当啊

@bestucan 开的是双精度的，之前用pbm都是自己写UDF的，这次用的自带的，第一次遇到这个问题。。。回头有空搞搞，这感觉是软件自身的问题，我装个低版本的c试试，现在装的是 2012，改明儿装个2010，我还是感觉2010比较稳定。

@cccrrryyy 感谢大佬答疑解惑，也谢谢您的分享:140:

多谢老师解答，问题已经解决了

各位老师，最近在研究双欧拉模型模拟液固两相流动，遇到以下两个疑惑，希望老师给予帮助:xiexie:
①想利用UDF对颗粒相剪切粘性系数和固相压力等参数进行修正（如26fbd023-2fb6-48f6-842a-733347fa08d1-图片.png ），公式中包含颗粒拟温度项，在编写UDF的过程中不知道如何获取该参数？？:134:
是否可以利用颗粒脉动速度C表示的颗粒拟温度7bf0c5a9-9020-484b-ba1e-2ba848de0117-图片.png 该方法中如何获取颗粒脉动速度？？？曾尝试用速度获取宏 C_U(c,t); C_V(c,t); C_W(c,t)得到颗粒速度，所得结果明显错误
②湍流模型中的湍流黏度是否要定义:141: 378f2ea3-f977-4a04-904a-05db21608467-图片.png
k-e湍流模型中湍流黏度定义为e94ab208-99c9-49db-9e65-664f5af4574b-图片.png
在输出湍流黏度时，是默认按照该公式进行计算的吗？？？只有有别于改公式时才需要写UDF定义吗:141:

你是不是算的倒虹吸啊

@veen 谢谢，我学习一下

哈哈哈

多相流的有限体积法跟单相流没什么区别啊

@李东岳 谢谢东哥，就是我把resolveFeatrueAngle调大之后画网格没有捕捉这条特征线，画几何的时候因为在两个不同平面上，这条特征线是存在的

@百善孝为先0626 可以的话还是想麻烦您把我的联系方式告知他，这边看能不能合作一下

@李东岳 有没有其他方式参考赠书呀？五行缺书~

我自己在用TDMA算法，用Matlab

@wsy11 谢谢您~我现在已经可以成功构建网格啦^^

@evensun

我猜你的疑问应该是，为什么类似下面这段代码编译器会报错？

1 pressureReference pressureReference(p, simple.dict()); 2 pressureReference pressureReference(pf, simple.dict());

这两条语句有4个pressureReference，我分别编号为1, 2, 3和4号。

编译器运行到第一行代码时，假设当前函数体内没有声明pressureReference，编译器就会首先去全局域（global ）里寻找pressureReference是什么，当找到pressureReference的声明，发现是个类时，编译器把它（1号）解析为一个类（class）。 编译器继续运行，遇到2号pressureReference，编译器主动构造类型为pressureReference（1号）的对象pressureReference（2号）。
注意：这两个名称一样，编译器会怎么做？当前作用域里的pressureReference（2号）会把全局域（global）里的pressureReference（1号）隐藏掉，1号就发挥不了作用了。也就是说，在此后的作用域里pressureReference只是个对象，只发挥2号pressureReference的作用。 编译器继续运行，遇到第二行的3号pressureReference，此时编译器将这个pressureReference解析为一个对象（2号）。 编译器最后发现4号pressureReference，这是什么，3号对象后面一个4号对象，蜜汁操作，编译器理解不了，编译就失败了。

那如果我一定要让3号pressureReference发挥类（class）的作用（像1号pressureReference一样），有没有办法？
有的。在第二行代码的3号前加作用域解析运算符（::），编译器就会去找全局（global）的pressureReference声明了，像上面的第一步一样，然后把3号解析为类（class），再把4号理解为对象。编译就能通过了。
有兴趣可以试试。

1 pressureReference pressureReference(p, simple.dict()); 2 ::pressureReference pressureReference(pf, simple.dict());