@冠竹simplefoam、pimplefoam、rhosimplefoam和rhopimplefoam的层流计算是不是也可以称为准直接数值模拟呢？

在v2012版本中，TEqn.H中有rhoCpPhi的定义为：

const surfaceScalarField rhoCpPhi ( "rhoCpPhi", rhoPhi*(Cp1 - Cp2) + phi*rho2*Cp2 );

这个定义是如何得到的? 如果Cp1 = Cp2，那rhoCpPhi = phi*rho2*Cp2，这个和 phi*rho*Cp有差距。

如果直接把中心点的rhoCp插值到面上和phi相乘应该表示为 fvc::interpolate(rho1*alpha1*Cp1 + rho2*alpha2*Cp2)*phi。

@杨英狄 这个也是计算粒子与壁面的作用，所以就是在particleerosion函数的基础上改的

@李东岳 李老师，fluent和其他机器计算时长的对比有吗

@cccrrryyy 让运维换了个of2106版本，速度有所提升，相当于自己电脑的四分之一速度，离期望值还是低点

@李东岳 烦请东岳老师有空看下，期待您的回复。谢谢！

@bestucan 谢谢大佬，已经解决了

看下tutorials里的Allrun和Allmesh操作步骤

@wenke 你好，你这个问题解决了吗？我也遇到了这种问题，现在也没有解决。

@qing 是这样的，openfoam完全可以这样任意折腾。它是基于文件的，只要能读取、识别到相应的文件就可以进行计算。

我思考了下，

matrix relax 虽然具体实现方式不同，但是本质不变。

你说的是对的。一楼我说的应该有问题（我更正一下）。不仅仅field relax，matrix relax，也会影响计算结果。我当时把matrix relax理解成不会影响计算结果了。

matrix relax是使得计算结果更慢的趋向于真实值，但换来一个对角占优矩阵

field relax是使得计算结果更慢的趋向于真实值

两个都使得计算结果更慢的趋向于真实值。

当时看这个公式的时候：$\phi^{n+1}=\phi^n+\beta(\phi^{n+1}-\phi^n)$，以为$\beta$变小时候，$\phi^{n+1}$还是原来那个值。然而并不是:136: :136: 好尴尬

@李东岳 李老师，您好！现在使用foam-ext 4.0,已经编译完成，计算成功一次，后续计算时一直报错（可以划分网格，计算中读取信息时出现问题）。
bed30e22-ed0d-4079-b2df-86b3002c91bf-image.png
7502d3ba-1a3d-4601-bd2a-3de6d249d09f-image.png

@zhouxu 谢谢，我尝试下。

@李东岳 东岳老师，请问这个速度能不能从10-2下降到10-4或者-6？

这里有介绍：link text

类似的需求：

How do I add some other package to the container?

The easiest way is to modify your Dockerfile and rebuild. Particularly easy if the software you want to add has an Ubuntu deb package.

里面也有 docker file 的示例。看起来这个比较可行

没有经验，一直在 linux 上跑

@李东岳 懂了，谢谢李老师

应该研究下为什么transonic需要提升对角占优

@李东岳 回答一下自己的问题。transonic下的压力方程存在对流项，对流项不利于对角占优矩阵的产生，因此transonic的压力方程要增加方程松弛。

http://dyfluid.com/rhoSimpleFoam.html

@ruibin_tongji
Green-Gauss Cell Based
9f30b62e-ad51-4bd3-ad87-c12ce78f2194-image.png
Least Squares Cell Based
8e14d4c0-98f9-4866-b57a-05a79fe25109-image.png
Green-Gauss Cell Based是最粗暴最直接的，但也许恰好适用于结构网格？我看这个Theory Guide里面是这么说的
38584698-16af-4fa1-8327-194895b34791-image.png

有没有测试过对流项的格式？

@李东岳 几万还是可以算的 几百万的网格我40核的服务器就算不动了

想请教各位，我对OF后处理中力系数forceCoeffs中的几个概念不太能理解，想请教一下各位。

forceCoeffs { type forceCoeffs; libs ( "libforces.so" ); writeControl timeStep; writeInterval 1; patches ( "floatingObject" ); rho rhoInf; log true; rhoInf 1000; liftDir (0 0 1); //升力方向 dragDir (0 0 1); //拖曳力方向 CofR (5 5 2); //用于计算力矩的点 pitchAxis (0 0 1); //pitch方向,围绕转动的轴 magUInf 0; //无限远处来流处的速度 lRef 2; //参考长度 Aref 8; //参考面积 }

首先这里对liftDir升力和dragDir拖曳力的定义是怎样的 ?举个例子，如果模拟的是一个物体在水中匀速运动，升力可否理解为浮力？那这里的拖曳力又是指什么呢？
还是这个例子，pitchAxis解释为围绕转动的轴，这又是指什么？另外如果在静水中运动，那magUInf可以设置为0吗？lRef和Aref是相对谁的参考呢，计算域吗？