感觉有人手把手带着学OPenfoam是每个openfoam初学者的梦想。。所以最近国内有没有相关机构或个人组织培训？:big_mouth:

好的，谢谢李老师！

请问有谁测试过OpenFOAM的并行效率？在使用酷睿或至强的工作站或集群上并行最多可以达到单核性能的多少倍？我用I7-4790和I7-3920XM测试的结果，用阿姆达尔定律推测最高是2.5倍不到，而且线程数超过物理核心之后就是负优化了，感觉不太可能，无法解释那些上亿网格是怎么处理的。现在正在鼓动老板至少搞一台高性能工作站，同等价位下是主频优先还是线程优先？优先考虑酷睿还是至强？

0_1512907622643_2.jpg
fluent中设置是这样的：
Boundary A, B and C, D
are set as periodic boundary conditions. Boundary E is set as no slip
wall. Boundary F is set as opening pressure boundary.
不知道在OF怎样设置好

FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 100000; boundaryField { inlet { type cyclic; } gasinlet { type cyclic; } outlet1 { type cyclic; } outlet2 { type cyclic; } left { type cyclic; } right { type cyclic; } TWall { type fixedFluxPressure; } atmosphere { type totalPressure; p0 uniform 100000; U U; phi phi; rho rho; psi none; gamma 1; value uniform 100000; } } dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type cyclic; } gasinlet { type cyclic; } outlet1 { type cyclic; } outlet2 { type cyclic; } left { type cyclic; } right { type cyclic; } TWall { type fixedValue; value uniform (0 0 0); } atmosphere { type pressureInletOutletVelocity; value uniform (0 0 0); } }

体积项用setField，但这样不对，求大神帮忙:big_mouth:

@hurricane007 Hi，对您说的学习笔记特别感兴趣，方便的话可否发送一贴学习交流一下，我邮箱是zhhu5634@Colorado.EDU，谢谢！

@nuaawqf 今天测试了dbnsTurbFoam自带的naca0012算例，最大CFL数在100依然能保持稳定，最大自适应步长7.14e-5。更大的CFL目前暂时没有测试。感觉使用隐式时间推进之后求解器的鲁棒性提高了不少。。

果然是边界条件的问题。现在已经能够运行了，等待比较结果中。。。

是不是我的问题没有问到点子上？请各位赐教。

@李东岳 问题已解决，直接在.H中添加了一段简单的二分法求解方程的程序，谢谢。

@sibo 其实大约都是基于虚拟化技术来的。我是给docker装的vncserver，虚拟一个主机，在的形式windows通过vnc访问。这样server和client

@青瓷碎了夏天 面积分可以分两部分，第一步先使用面两边的体变量插值出面上的物理量，比如一个面k，它左右两边为i，j两个控制体，则P_k=interpolate（P_i,P_j）。这个值相当于面上的平均值，与面矢量内积就能够得到面积分了。

这个问题使用fluent算貌似同样会出现这个问题，其实这时候已经出现了部分空化，处理不当产生负压。

@100yearsalone 恩，最近没过来逛逛咱们这个论坛。由于OpenFOAM官方版本一直以pressure-based为主导，所以会产生很多迷惑。

fluent 里的porous 边界如何在openfoam中实现

大家好，之前我发了一个帖子，大意是基于interFoam，在连续性方程中添加源项实现质量的输入和输出，进而造波：
http://www.cfd-china.com/topic/943/添加源项后相方程求解越界-请问这种情况有人遇到过吗

经过多次调整，质量源造波的程序已经基本上完成了。我按照上贴子的公式，把源项加入相方程中并开启了MULES修正，同时在压力方程中也添加了源项，这样求解的时候，质量源区域内的alpha值（默认为1，源区位于水下）大概在1.0*（质量输入）-0.9*（质量输出）之间波动（因为我的源项是类似于正弦波动的）。

然而，我采用相同的方法修改压力方程，使用动网格求解器interDyMFoam的时候，上述的方法出问题了，形成了类似于空化的现象（不懂空化，希望没说错），即相体积分数变得很小，最直观的描述是在水下质量源区域会形成一个气泡，后面相体积分数的求解救完全错误了。

我仔细看了一下求解器的文件，发现interFoam.C和interDyMFoam.C求解时主要的区别在于动网格求解程序中多了一些通量计算以及压力方程中的部分不同。但是由于之前没看过动网格的求解器，因此还是看的不太懂，请问有做过这方面的大牛吗？可否分享一下经验，为什么网格不变形的时候完全正常，但是用了动网格求解器就不行了。:big_mouth:

感激不尽:sunglasses:

是啊，相信CFD在国内也会有不错的发展！

因OpenFOAM现有格式不满足需求，需要实现新的插值格式，有没有熟悉的人说明一下思路，此外我不清楚MeshObjects文件夹下的 upwindCFCCellToFaceStencilObject和upwindCPCCellToFaceStencilObject的作用和区别是什么？谢谢

surfaceScalarField cPhi("cPhi", c & mesh.sf()); fvScalarMatrix c ( fvm::ddt(c) + fvc::div(cPhi) == source. );

@wwzhao 谢谢

@小龙 filter 搜一下就可以找到