@bestucan

如果是计算时要实时计算边界条件，这个可能符合你的要求：
http://openfoamwiki.net/index.php/Contrib/groovyBC#Usage_of_the_boundary_condition

感谢回复，需要的确实是在计算中实时计算的边界条件，我先了解一下:xinxin:

又仔细看了一下，两个版本应该是差不多的。
相对速度是因为CdRe不是曳力系数 而是曳力系数×雷诺数 所以可以用粘度、直径还有密度换算成相对速度。

@李东岳 老师，内部也在动的，截图可能看不出来，只是这个边界更明显。
75376c62-5086-4608-879f-b7823f7ebaa8-1639391378.png
7a6ca8ef-4d95-4e82-b6e1-77c1203adfab-1639391406(1).png

@bestucan 是的，原因应该是相似的，感谢您的解答。

@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方程笔记》里面。