CFD中文网

跟进一下，这是不同边界条件下的压力云图。可以直观的看到入口处压力场出现震荡。
c908d44b-b2db-44fb-91fb-276fddda4220-image.png

我也尝试从数值角度说一下我的看法。计算求解的是NS方程，对于入口处速度和压力均为定值的初始条件，可能会导致结果不满足NS方程。造成所谓的过度约束。而入口出口分别指定速度和压力为固定值，则会避免过度约束的现象发生。

欢迎各位大佬批评指正！

@xpqiu 好的，谢谢您！

@Mina-Lee 好的，谢谢啦。我刚刚试了很久终于试出来了，确实是少给了参数。

学OpenFOAM至今一直水在某大佬的OpenFOAM千人群中，个人认为还是比较活跃的之一(毕竟也水到50多级了)。
但是今天早上一起床，发现我被群内某管理员踢掉了：
qq_pic_merged_1629032541800.jpg

我寻思我也没发什么违禁言论，也没做什么违法群规的事情啊，且前天我还在群里解答疑惑:134:
关于这个问题我私聊了一下管理员，结果发现我问完他直接给我删了？:136:
30bfc521-cae5-4276-8a6b-89c6da167351-Screenshot_20210815_210758_com.tencent.mobileqq.jpg
145f5c52-5a8e-47e9-9341-302bf79191c7-Screenshot_20210815_210832_com.tencent.mobileqq_e.jpg

我不是很懂这波操作，按理说踢人不应该有个规范吗，还是按管理员想踢谁就踢谁吗？

PS：不知道这里是否允许对这类问题进行讨论，侵权麻烦管理员大大直接删帖吧

@赵一铭 整个程序我用的SIMPLE算法，所以连续性方程改成了压力修正方程来求解压力，不能直接求解出密度，怎么办？

取平均之后计算很快。 如果你需要瞬态流场以及瞬态的各项数据，你就用瞬态。 如果你只是想得到一个平均流场，那就用平均。

是的，在新版除了createFields.H没了，fvCFD.H都没了。变量都在class里面初始化，比如VoFSolver。H：

class VoFSolver : public fluidSolver { protected: // Phase properties //- The compressible two-phase mixture autoPtr<VoFMixture> mixturePtr_; VoFMixture& mixture_; //- Name of the alpha convection scheme const word divAlphaName; //- Phase-fraction flux Courant number scalar alphaCoNum; // Kinematic properties //- Velocity field volVectorField U_; //- Volumetric flux field surfaceScalarField phi_;

首先，把计算域特点介绍一下：本算例是为比较圆柱绕流问题与三维螺旋结构绕流的不同。
本研究中的三维螺旋结构是由三个圆螺旋拉伸形成的结构，有点类似于麻花（但是麻花是两个圆柱缠绕的）：
0_1510012974209_p2.PNG
将流体域与该结构做一个布尔运算，得到的流体域的结构如下图所示：
0_1510013093878_p3.PNG
其中左侧为速度入口，使用的是LES模型要求壁面的Y+=1
————————————————分界线——————————————————
圆柱绕流的网格画分比较简单，可以通过O型网格得到较好的结果。
但是对于新结构的网格可能不太适合，不能一通到底的使用O型剖分。
下面是我尝试用workbench 中的Meshing画分的网格的俯视图：利用的是sweep方法。
0_1510013238698_4.PNG
现在想利用ICEM来画分，是想能得到更好的网格？请问大家有没有什么想法

https://stackoverflow.com/questions/30304794/gnuplot-command-with-multiple-with-arguments

大佬这是输出系数矩阵了吗

@wdhj 试过了没有？结果怎么样呢，请给反馈

chemkin原来使用的输入文件中的传输数据chemkinToFOAM识别不了，需要将组分的传递数据用openFOAM的格式自己写一个transportProperties

标准的限制方法是保证插值后网格顶点（vertice）的值不会出现极大值或极小值（比周围网格的最大值大或比周围网格的最小值小）因为梯度插值最容易出现越界值的就是在顶点处，从而用一个标量来scale梯度gradient；另外一种方法是保证网格的面心midponts of cell faces处不出现极大值或极小值即可，就认为满足了插值后不存在越界的问题。我想的是关于第二种方法，面心处插值后不出现最大值或最小值不能代表顶点处不出现最值，对吗

谢谢，我刚把这个导入到程序中，但结果仍然不理想，解出得P'偏大，在修正速度时会导致速度变成负的，而且是进口速度得负几十倍，所以总感觉压力方程和速度方程不匹配，肿么破？？:crying:

@cfd-china link text!![alt text]0_1492654992588_QQ截图20170420102158.png 这是这个文献上说的在固壁上的条件；还有一篇文献直接在壁面设置为0了link text0_1492655090308_QQ截图20170420101901.png

请问拉格朗日的TomiyamaLift如何进行编写

TomiyamaLift { sigma 0.072； }

我这个写法openfoam报错了

是的，我有的时候也会出现这个情况。点pre-Mesh，体网格丢失，不知道还有没有其他人遇到这个问题？

买nas，养pt

ffmpeg -i input.MP4 -c:v libx264 -preset slow -crf 40 -c:a copy output.mp4