CFD中文网

已经找到了论文Dynamic meshing around fluid-fluid interfaces with applications to droplet tracking in contraction geometries作者的代码，使用OpenFOAM2.1版本测试了一下结果没有问题，代码下载位置，大家可以下来试试。

这个是我运行作者给的算例的结果：

@harvee

@浪迹天大 说是速度和精度要高一些，还可以进行显焓分析等，可以看一下这本论文的第五章。 Implementation and Validation of a Solver for Direct Numerical Simulations of Turbulent Reacting Flows in OpenFOAM

我研究的问题是一个气固两个域耦合问题，在气相中有化学反应。在迭代耦合边界满足温度和热流密度守恒时，化学反应计算出现了问题，所以研究一下这个方法是否可行。昨天刚把OpenFOAM-Extend1.5 + Cantera1.8的编译成功，看能不能改写到extend 3.2版本的。 要是能要到上面论文里的代码就好了，他是2.0版本的，改起来就快了。

@cfd-china 是可以实现的 国外有一个OpenFOAM Extend 1.5 + Cantera 1.8的，不过我还没有调试成功，我用的是Extend 3.2的。

有人在用OpenFOAM + Cantera来计算化学反应吗，我现在在搞这个，有的话交流一下？ 还有请教一个问题，OpenFOAM中化学反应是不是不能是定常的？

@harvee @harvee 不好意思一直没回复 ，最近都没弄了，在搞化学反应一直又问题，你有什么进展吗？cfd-online上那个帖子有回复，我一直也没看。cfd-online二维自适应

@liqiwen 还可以对结构网格进行二次开发， [基于FLUENT软件的导弹意外点火情况下舱室维护仿真平台开发

](http://d.g.wanfangdata.com.cn/Thesis_D429353.aspx)这个论文的3.1.2节有讲用index来对ICEM CFD二次开发的功能。

@wwzhao 再请教一个问题，您是怎么知道这个输出的位置或是这个函数的原理的，我一开始通过输出、调试这个函数并没有搞懂它是怎么回事。很想知道您当时是通过什么途径来弄清楚的？ 感觉好多OpenFOAM的函数分析起来不是很好弄清楚，可能是我对OpenFOAM的原理和C++理解掌握的还不够:sad:

@wwzhao 嗯　看了一下好像是这样，谢谢解答:happy:

@wwzhao Info<<os<<endl;　是我自己写的，是输出到终端，现在想知道os.writeKeyword("neighbourFieldName")<< neighbourFieldName_
<< token::END_STATEMENT << nl;是输出到什么地方去了。

@wwzhao 嗯，我调试了一下，确实向您说的是runTime.write()这个函数触发的。我用Info输出了一下os发现是数字１，按理说是写到一个文件里，但运行后，case目录下并没有输出相关文件。在想是不是要在etc目录下的controlDict中设置一下DebugSwitches开关，不过还是没成功。这个函数的输出可以改成Info输出，显示一些信息，但是很想知道这个函数的输出是怎么回事。

OpenFoam中的边界条件里面都有一个函数write，用于输出一些信息，如oscillatingFixedValue这个边界里面的write函数：
这个函数是怎么触发的呢？

@harvee 嗯 ，这两个我也看过。直接修改边界，目前我也没成功，我看hexRef.C代码里面，没看到有说对于某个特定的边界，可以不切割Z方向的网格，还在研究中。现在还在看代码，如果能成功我会分享的。希望能快点搞出来:laughing:

还有一个14的, ANSYS 14.0 LINUX

@cfd-china
ANSYS 17.0 Linux 你试试.

请教个问题如何上传图片,个人用户有相册集可以把图片上传到服务器上吗?

在论文DYNAMIC MESHING AROUND FLUID-FLUID INTERFACES WITH APPLICATIONS TO DROPLET TRACKING IN CONTRACTION GEOMETRIES第58页最上部有写cAlpha的作用。

The cAlpha keyword specified in the file <case>/system/fvSolution is
a factor that controls the compression of the interface, where 0 corresponds to no
compression. After solving the volume fluid equation, the density and viscosity are
modified using the new values of α.

@cfd-china 装一个Linux版的ANSYS:laughing:

@浪迹天大 The OpenFOAM Technology Primer 这本书动网格一章里面有讲自适应网格的，OpenFOAM求解器大多可以添加自适应，而且步骤比较简单。:happy:

@cfd-china 当时Z方向上的一层网格有被加密码？

我试过采用symmetryPlane和cyclic边界，但是三维方向上网格还是被切了，我现在只想让x和y方向切网格。这个pesudo-2D是指在三维方向上采用类似symmetryPlane和cyclic边界条件处理吗？ 我找了一下 ，好像是这个意思。 谢谢你的回复@administrators

参考论文Dynamic meshing around fluid-fluid interfaces with applications to droplet tracking in contraction geometries附录B 里面的代码，修改OpenFOAM自带的三维自适应网格算法，不过没有成功。我选用了一个最简单的2D导热问题，来调试这个代码，后俩发现在hexRef4.C（这个文件原来是hexRef8.C用于处理三维问题）的函数createInternalFaces里有一个判断if(edgeMidPointI==-1)，这个出现了问题，而论文里并没有写这个，edgeMidPointI这个变量是传递到函数storeMidPointInfo这个函数中，所以现在想改写一下这个函数，但我对这个函数的原理还没有看明白，不知有人实现了OpenFOAM下2D自适应没？ 想交流一下，如果最终能把代码调通，我会把代码也共享一下。