不建议开启`体积力面插值`算法



  • 目前OpenFOAM基金会,ESI-OpenCFD(OpenFOAMv1706)下的体积力面插值代码提出的本源在于使得较大体积力下的计算压力无震荡且稳定。然而,目前这个算法并不能够提供合理的数据,和实验甚至通常的算法偏离极大。个人觉得主要有以下缺陷:

    1. 速度处理过于激进:在UEqn.H中将时间导数项移动到压力方程,虽然从求解方程的角度是合理的。但是目前这种方法是否是数据偏离的原因暂不明确;
    2. 力的处理任意性:力是不具有守恒特性的,因此基金会仅仅将力插值到面上来处理,再一次,虽然从求解方程的角度是合理的。但是目前这种方法是否是数据偏离的原因暂不明确;

    我不建议开启体积力面插值算法,目前这个算法需要大量的验证和提高。个人主观来说,发布的这个算法计算的是完全错误的结果

    同时,合理的计算体积力需要更好的算法。



  • @李东岳不建议开启`体积力面插值`算法 中说:

    力的处理任意性:力是不具有守恒特性的,因此基金会仅仅将力插值到面上来处理,再一次,虽然从求解方程的角度是合理的。但是目前这种方法是否是数据偏离的原因暂不明确;

    保守力是可以移到压力方程中的,整到面上没啥不对。

    估计处理的时候体积力的旋度没整成0吧。可以测试一下体积力的旋度的量级。



  • 某些体积力的旋度必然不为0.


登录后回复
 

与 CFD 中国 的连接断开,我们正在尝试重连,请耐心等待