@shouan 这是我的pointMotionUx文件
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: v2006 |
| \ / A nd | Website: www.openfoam.com |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class pointScalarField;
object pointMotionUx;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type uniformFixedValue;
uniformValue table
(
( 0.000000 0.001250 )
( 0.020000 0.001300 )
( 0.040000 0.001500 )
( 0.060000 0.001550 )
( 0.080000 0.001700 )
( 0.100000 0.001900 )
( 0.120000 0.002000 )
( 0.140000 0.002200 )
。。。
。。。
);
}

bottom { type zeroGradient; } reef { type zeroGradient; } outlet { type fixedValue; value uniform 0; } atmosphere { type zeroGradient; } Back { type empty; } front { type empty; }

}

@zhe 好久没来了才看到回复，感谢！！

@袁宝强 我记得可压缩求解器里面有个可压缩系数，psi，不知道是否有用

@队长别开枪 非常感谢您的帮助！

@李东岳 谢谢李哥

不清楚 过几天准备写个弹簧试试看能不能准

@wwq 在 请教关于k和epsilon壁面函数的问题 中说：

https://www.openfoam.com/documentation/guides/latest/doc/guide-bcs-wall-turbulence-epsilonWallFunction.html

非常感谢:146:

@李东岳 谢谢李老师！A fifth-order shock capturing scheme.pdf 李老师，我是想把这个文献里的第二部分模块数值算法植入进去，我昨天上传的是第二部分前面一块，后面那一块主要是进行重构，去计算q_L和q_R，然后通过前面的公式去进行计算，前面那块我看和您发的rhocentralfoam解析中很像，但后面的公式不知道应该怎么植入，没有思路，想请教一下老师:chouchou:

@李东岳 因为换个地方继续算，重新分块的时候就根据最后时刻的网格来划分，所以用scotch 就是每个块基本同样网格数（上下几十浮动），而用simple划就会是那种XYZ方向的拉升

lz解决了吗？我也遇到同样的问题了，高压气流冲击真空环境，使用rhoCentralFoam，算了一段时间之后就出现negative initial temperature t0，改了好久，可以推迟出现的时间，但总会出现，不知道怎么弄了

@李东岳 非常感谢东岳老师的回复，这个问题出现在均匀网格，有的时候甚至同样的位置，但是气域的速度不同，他最后采样的范围也不一样。。。

@turbulence 02525037-654f-4a52-bc51-dcef3cc03688-image.png

@李东岳
感谢李老师的回复，前段时间做了其他的事情，就搁置了这个。哈哈
这个问题解决了，是因为边界条件的设置。
另外请教李老师另外的问题：

OpenFOAM中有求解瞬态不可压多孔介质的求解器吗？ 或者关于二相流中造单纯的流，关于相界面处的边界条件设置问题，李老师有更好的建议吗？
谢谢。:146:

@shiyu 在 BoundaryFoam求解近壁面流动（倾斜壁面） 中说：

已知local坐标系下的值，想转换成global坐标系下，方便给定openfoam的边界值，openfoam里有可以直接利用的tools吗？

没有现成的。不过你写代码的话，也不是特别费事。

http://www.dyfluid.cn/theory.pdf

合适的内容我会更新在《无痛苦NS方程笔记》里面。

最近把理论这部分弄了一下，

捕获.PNG

libs ("libmyfvOptions.so"); 这句不能像一些官方提供的案例那样放在fvoptions文件里，得放在controldict文件下
Screenshot from 2021-12-03 17-31-56.png

@zhouxu 我就是设置的滑移边界，存在槽道流特征

测力时出现很强烈的震荡 想请教下是什么原因导致的 黑色的是想要验证的图（没有震荡），红色的是我计算的。用的Openfoam自带的测力工具。
17c6ea16-98db-4fc2-8b0f-25c32d22a333-image.png
ea932f4f-4999-469d-bdb4-d254f739078e-image.png a83847ec-c340-453c-a6c1-37a0913808f8-image.png

@李东岳 好的，感谢李老师回复，我们服务器系统是centos6.5，好像不太支持新版本:136:

@lyt 难道是后面一定要加两个（）（），才行？这是什么逻辑

@caijinjin 是计算停在这里还是终端会停止，方便的话可以交流一下吗，我微信15098932699

按说，那个 pdf 的大部分勘误记录都会跳转到这里
https://cfd-china.com/topic/2338/无痛苦n-s方程笔记-勘误/38

“错”和“误”不一样，（所以某些理论上无法消除、工具条件限制带来的与事实不符的差距叫“误差”不叫“错差”）（我对你措辞没意见，只是看你挺有互助热情，也“互助”你一下，说给你听听，参考参考。）

可以直接说出笔误位置，细微的笔误难以发现（容易发现也不会误了）。比如：

"接下来定义温度 T 在 (x1, y1, z1) 点的物质导数为移动的无穷小微团通过 (x1, y1, z1) 点的时候"
这句第二个坐标（或者第二个（x1, y1, z1））标号应该是(x2, y2, z2)

我看你之前也发过类似的帖子，可以看完了一遍书，把发现的笔误统一发一下，也方便修改。

欢迎继续:146:

@chaoscfd 我用了 但还是mesh.findcell最快

我猜是这样的，
process 里必须有动网格，这样计算的过程只涉及网格如何动，

如果 process 里没有动网格，那么当动网格运动到该区域，还要造动网格。

至于 process 里没有动网格也能算，那是动网格没有运动到该 process 覆盖的区域。

可以测试，在计算区域内设置动网格，但是不让它动，比如平静的水面。应该可以一直算。

paraview 里可以看每个 process 的区域，算不动的那步估计就是跨区域运动了

顶一下

这个是安装 docker 的问题，不是 OF 的问题 :wocao:
去 docker 相关社区问更容易得到答案:chigua:

@tens 再次感谢您！重新尝试了一下，可以输出了！

应该影响不大。这种颗粒应该不多，不会所有的颗粒都恰好那样。所以对alpha的影响也很小。除非颗粒直径非常大。不过直径非常大的话，那就会引起收敛性的问题。

b75ecf32b61ea077acecbb62002a6c8.png
这是我后来自己的解释

要做简单例子的定量对比，来拿验证方法可靠性的，涉及到水破碎分离的，列如dambreak就很经典 但是那个实验没有涉及到破碎

@fu 在使用laminar（Quasi-DNS）计算的时候，壁面处的网格Y+<1 ,且整个区域网格也要相对较小；当使用湍流模型的时候，壁面处理跟laminar应该不同，使用壁面函数之类的；如果题主使用一套网格应用于两种流动模型，是否会导致计算结果的偏差呢？还请大佬指点

@李东岳 好的，谢谢老师！具体的公式我也在查，文献中没有直接给出

@李东岳 谢谢老师回复~
老师可以帮我看一下这个错误怎么修改吗，按照您之前的一个帖子，定义了一个随时间变化的热流边界~

HEATED { type externalWallHeatFluxTemperature; mode flux; q coded; // Heat flux code #{ scalar q=*this; scalar t = this->db().time().value(); if(t >= 0 && t<=1.5) { q = 180000; } else if(t<=3) { q = 200000; } else if(t<=4.5) { q = 250000; } else if(t<=6) { q = 300000; } else { q = 0; } q.writeEntry("", os); #}; value uniform 300; kappaMethod lookup; kappa "K"; } /home/dyfluid/OpenFOAM/dyfluid-v2006/run/bianhuaQ/0/SOLID/T.boundaryField.HEATED:32:23: error: cannot convert ‘const Foam::PatchFunction1Types::qPatchFunction1ScalarField’ to ‘Foam::scalar’ {aka ‘double’} in initialization /home/dyfluid/OpenFOAM/dyfluid-v2006/run/bianhuaQ/0/SOLID/T.boundaryField.HEATED:33:31: error: ‘const class Foam::PatchFunction1Types::qPatchFunction1ScalarField’ has no member named ‘db’ /home/dyfluid/OpenFOAM/dyfluid-v2006/run/bianhuaQ/0/SOLID/T.boundaryField.HEATED:55:9: error: request for member ‘writeEntry’ in ‘q’, which is of non-class type ‘Foam::scalar’ {aka ‘double’} /home/dyfluid/OpenFOAM/dyfluid-v2006/run/bianhuaQ/0/SOLID/T.boundaryField.HEATED:55:24: error: ‘os’ was not declared in this scope; did you mean ‘cos’?

基本上，不发散能对的上就行 :mihu:

刚盖是不动的盖子? 速度为0不就可以了么

@氕氘-川 批量保存的脚本老报错： ERROR: In C:\bbd\df0abce0\build\paraview\src\paraview\ParaViewCore\VTKExtensions\Default\vtkCSVWriter.cxx, line 82
vtkCSVWriter (000001F6DAFC08C0): Unable to open file: /-0.1150.csv
ERROR: In C:\bbd\df0abce0\build\paraview\src\paraview\ParaViewCore\VTKExtensions\Default\vtkCSVWriter.cxx, line 82
vtkCSVWriter (000001F6DAFCC8C0): Unable to open file: /-0.0650.csv
...

这是我的脚本：
import numpy as np
from paraview import numpy_support as ns
import paraview as pv
import paraview.simple as ps

afoam = FindSource('a.foam')

for y in np.arange(-0.115,0.115,0.05):

# Make a cut of the wing to have a local 2D profile Slice1 = Slice( Input=afoam ) Slice1.SliceType = "Plane" Slice1.Triangulatetheslice = 0 Slice1.SliceOffsetValues = [0.0] Slice1.SliceType.Origin = [0.0, y, 0.0] Slice1.SliceType.Normal = [0.0, 1.0, 0.0] descriptiveStatistics1 = DescriptiveStatistics(Slice1) descriptiveStatistics1.VariablesofInterest = ['T'] descriptiveStatistics1.AttributeMode = 0 Render() writer = CreateWriter(str(y)+".csv",OutputPort(descriptiveStatistics1, 1)) writer.FieldAssociation = "Cells" writer.UpdatePipeline() del writer #SaveData(str(y)+".csv", proxy=OutputPort(descriptiveStatistics1, 1), FieldAssociation='Cells') descriptiveStatistics1.UpdatePipeline()

@wsxfyy 你可以多试试，有的可以改，有的会报错

@李东岳 多谢东岳老师，重新划分网格后可以运行了

0b9f25fa49feec7e6e3148b3f3df445.png
这个意思是说我Foam：：twoPhaseMixtureThermo：：New没定义？
2099f75c499411d7558623e7b73b800.png
这个应该算定义吗？能通过报错中的（.text.startup+0X35ed）找到为啥错吗。。
有遇到过相同问题的朋友吗

试试zeroGradient

@李东岳 好滴好滴，非常感谢，我这就试试:xinxin:

各位前辈中午好，coalChemistryFoam自带case是求解两种物质：煤和石灰石，我想添加新的需要反应的颗粒，请问需要怎么修改求解器呢？谢谢大家。

@李东岳 在 Simple算法中通量phi的计算 中说：

也就是说在边界处phi是不正确的。不知道是我没有找到phi边界修正，

边界处的phi不需要计算，是用户给的。加入速度1m/s，乘以网格边界面积，就是phi啊，并且不变。

如果是pressureInletVelocity 边界呢？ 代码里可以看出pressureInletVelocity 需要用边界上的phi来更新出新的边界速度。如果是这种边界边界上的phi又应该如何确认呢？

@李东岳 李老师，您好，我按照您的建议进行了相关的测试。结果发现当更换为GAMG时计算比PCG慢非常多，如果我在计算一开始取残差与相对残差值较大时，完全不收敛的情况基本上到500步左右就改善了。但是我又发现一个特别奇怪的现象：
首先通过pimpleFoam (nOuterCorrectors=1) 计算到一定步时计算结果收敛很好，残差基本下降至1e-5左右；但是当我在运行时更改nOuterCorrectors=2时，初始残差立即上升至0.4左右，计算异常缓慢，并且继续运行算例残差也下降不了。请问这个可能是什么原因导致的呢？

替代文字

你说上面这个图。这个算力对应的不是这个图。所以不一样。不过如果你增加近气速度的话，应该也是可以的。你加大10倍试试。

您好，请问这个问题您是否有好的解决办法了吗？
强制设为0可能不太对，反映出来的可能不是压实效果。。

@李东岳 东岳老师，这个主要是因为我的是粘弹性流体，在动量方程里边加入了弹性相，不只有雷诺数，还有魏森贝格数，可能是这两个无量纲数导致的:135:

4年过去了，当时发完这个贴就出去玩了。:jingya:

@ruibin_tongji 你看你的log是不是有些迭代次数太多了