CFD中文网

@东岳 感谢李老师！！是我边界条件设的不合理加上网格密度取值的问题。已经更正。谢谢李老师耐心解答！！不可压的alpha就是groMetricOneField，我想太简单了，以为稳流模型把空气自动过滤掉了。

云图红蓝色是速度，箭头是速度方向（无大小），alpha水没有显示，但是和左侧入口风的界面只差了6个网格，是几乎一样高的。

@东岳

抱歉，我传了两次图都消失了~
看到过博士论文如下图右侧，用kwsst加过，并且成功了。所以我想能不能把稳流模型alpha去掉，但这样会不会造成VOF界面处有问题了呢？OF这么写纯粹为了格式统一还是什么原因呢？

哦好像是定义没有弄明白，miu和nu是不一样的····

@东岳 谢谢李老师，我感觉是因为稳流模型中的alpha导致了空气项的速度老是算不对····而不是其他原因。我一个幼稚的想法，把稳流模型中的alpha去掉···这样是不是会导致空气-水的界面出现问题？···

看到过用一个稳流模型算风-浪的文章，但是文章给出的稳流模型和OF中的感觉有点不一样，就是alpha这块。不知道OF这么写用意是什么···也不清楚我这么一改会不会有别的问题···

OF1906，关于divDevRhoReff(rho,U)的定义，是不是有些问题，下图是OF中代码以及我认为应该对的代码，请问老师是不是这么写的？？

还有就是另一个问题。稳流模型中的alpha和这里的alpha是不是都会导致水-空气都存在的时候空气界面的速度会算不准。一个比较幼稚的想法，如果我删掉alpha，会不会就能算准了···？如果就删掉alpha，会导致空气-水交界面出现问题么····

想在空气项inlet入口加入风速剖面，实现风浪场。现在就做了一个简单的算例想测试一下，但是发现风场有一些奇怪，图片如下。

水深5m，空气高度3m，水槽长度20m。左侧inlet在空气界面加了一个（5,0,0）的速度剖面，但是稳定后发现风向没啥问题，但有一种向上走的趋势。请问各位老师，可能是什么原因导致的呢？OF1906，interFoam；
边界：二维
inlet: U fixedValue（空气界加风场）; p fixedFluxPressure; alpha zeroGradient;
outlet: U zeroGradient; p zeroGradient; alpha variableHeightFlowRate;
bottom:无滑移；
top：U pressureInletOutletVelocity; p totalPressure; alpha inletOutlet;
我感觉边界条件好像没什么问题。暂时没有什么头绪，请各位老师指教~谢谢！！

@东岳 明白了！谢谢老师

@东岳
请问一下老师，为什么将gradSchemes格式改为迎风等就会出现beyond EOF这种错误 ，而只有中心格式不会···就是用的simpleFoam

看了老师~ 是否可以这么理解？UEqn.A是动量方程离散完矩阵装配完后矩阵的对角线值，方程右端的压力梯度项量纲为m·(s-1)，所以UEqn.A实际量纲为s-1。而UEqn.H是这个方程的b，求出HbyA后为了得到面上的phiHbyA，求压力泊松方程。 在simpleFoam的泊松方程中没有找到类似
的这种思想~好像simpleFoam泊松方程求出来的直接是p*，而不是p'

抱歉 倒数二三行是OF上对这两个函数的解释~没有太看懂

请教一下各位老师，请问simple循环中UEqn().A() 和UEqn().H()分别是矩阵的什么呢？A的单位是s-1，H单位是m s-2

对于这两个函数的定义 也不是很了解 谢谢老师~

为了修改并行方便 我把所有代码减少了，即使是做这样的循环（forAll 对某个特定的patch面进行遍历计数），单核没有问题，网格数都是对的上的。但是一并行，就会卡在循环这步，不报错也不继续，请问老师是什么原因导致的呢？

你好老师～我也遇到了类似问题，用forall循环遍历每个网格，并进行累加操作。单核没有问题，多核就卡住，不报错也不继续。请问老师问题解决了吗～

OpenFoam版本为1906

修改位置主要在newMark中~

最近在修改完代码以后出现了并行相关的问题，代码如下：

tempx 和tempw都是vector型，if中的是对某个特定的patch网格面上的量进行操作，并将相关值赋值到tempx和tempw中，returnReduce是想并行将每个process相累加，但是发现编译后单核没有任何问题，但是并行的时候卡在这一步了，不报错也不往下走，如图：

是因为目标patch并没有分到所有process中因而卡住还是什么原因呢？
希望老师能指点一二，祝各位老师新年快乐~

@freevector 谢谢老师~ 我找到了一个方法，用lookup指针去找model_.time().lookupObject<fvMesh>("region0")，您的方法我有个地方不太理解，就是runTime提示没有申明，不太明白什么意思~

请问想要rigidBodyDynamics库中调用到某patch网格面上的信息。现在想到的方法是在库中找相关的继承函数，再类似mesh.boundaryMesh()。但是这个想法没能实现。请问这种思路是否可行。或者OpenFOAM中是否存在类似直接调用constant/polyMesh中网格信息的操作~希望能给出建议~