@bestucan n 谢谢教授的解答

@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