CFD中文网

@jpzhao 我没有试过cht用这样网格，应该是没有问题的。

@jpzhao 是的，但是交界面网格质量不是很好，建议使用 mergeMeshes，但是tposet好像无法做出光滑的曲面。

@oitocfd 您可以大概说一下在中间遇到了什么困难吗，目前我的大概思路是 fsi 的固体求解就不动，流体这边开发好作为一个类再合并进去

@number44 计算柔性体入水且伴随空化的案例，solid4Foam 加入 overset 功能后不能并行计算，所以并不打算在 solid4Foam 上做算法改进。

@bestucan @number44感谢谢前辈指导， 因为两个版本输出的结果文件格式是一样的，我打算尝试在 fsiFoam 和 overInterDyMFoam 之间做一个接口进行两边数据的更新，类似 Mpcci 这样的工具。

@bestucan 感谢前辈！

改了命名空间，貌似不起作用，是什么原因？

@李东岳 好的，谢谢

本人 c++ 能力有限
fsiFoam 为 foam-extend 4.0 版本上可以成功编译，我需要将 of 2106 上 overInterDyMFoam 作为 fsiFoam 的流体求解器，在把 fsiFoam 升级到 of 2106 当时，发现两个版本头文件的包含、类的名称及定义有诸多不一致，编译时出现多处重定义以及重复声明以及其他的错误，希望有大神们能提供 fsiFoam 升级到 of 2016 上的相关资料、经验或者思路，感激不尽！

@zhouxu 流体材料用的是常数密度，没有用理想流体，所以不需要用绝对压力。

看了一下 N-S 方程，计算前需要将 p 的值除以密度 ，计算完成之后 p 的结果需要乘以密度，简单验证了一下，和可压缩的计算结果基本一致了。

@zhouxu 相对压力为负压

@cccrrryyy 您的意思是说将计算结果中的压力乘以密度对吗，还是说 0 文件当中的 p 需要进行特殊的设置？谢谢。

我在使用 simpleFoam 计算水压时，发现其计算结果比 fluent 小了将近1000倍，使用 rhoSimpleFoam 时，计算结果与 fluent 相近。simpleFoam 已经将 rho 进行了无量纲化，是不是不适用于计算密度不等于 1 的流体材料的压力？或者是大佬们有其他解决方案？谢谢！

@袁宝强 内存多大？

问题解决了，之前是没发现 flowRateOutletVelocity ，把 flowRateInletVelocity 换成 flowRateOutletVelocity 即可

@李东岳 这个网格可简单了，应该不是网格的问题。

    Max skewness = 1.794448 OK.
    Coupled point location match (average 0) OK.
Mesh OK.

大家好，我在使用 simpleFoam 计算外流场时，利用质量流出口作为驱动，计算结果在出口处出现了奇点，以下为最新设置的 U 和 P，供各位参考，奇点的出现原因是什么？是否和边界类型有关？希望能收到大家的回复。

U文件如下

FoamFile
{
    version                                 2.0;
    format                                  ascii;
    class                                   volVectorField;
    location                                0;
    object                                  U;
}
dimensions                              [0 1 -1 0 0 0 0];
internalField                           uniform (0 0 0);
boundaryField
{
    outFlow_car
    {
        type                                    noSlip;
    }
    outFlow_inlet
    {
        type                                    pressureInletUniformVelocity;
        value                                   $internalField;
    }
    outFlow_wall
    {
        type                                    noSlip;
    }
    outFlow_bottom
    {
        type                                    noSlip;
    }
    outFlow_outlet
    {
        type                                    flowRateInletVelocity;
        massFlowRate                            constant -10;
        rho                                     rho;
        rhoInlet                                1;
    }
}

P 文件如下

FoamFile
{
    version                                 2.0;
    format                                  ascii;
    class                                   volScalarField;
    location                                0;
    object                                  p;
}
dimensions                              [0 2 -2 0 0 0 0];
internalField                           uniform 0;
boundaryField
{
    outFlow_car
    {
        type                                    zeroGradient;
    }
    outFlow_inlet
    {
        type                                    fixedValue;
        value                                   $internalField;
    }
    outFlow_wall
    {
        type                                    zeroGradient;
    }
    outFlow_bottom
    {
        type                                    zeroGradient;
    }
    outFlow_outlet
    {
        type                                    zeroGradient;
    }
}

谢谢！

可以试试cartesianMesh，用的是 openMP 并行技术，需要设置一下环境变量（export OMP_NUM_THREADS=n），不过没试过上亿量级的网格划分

@lts您好阿，最近我在使用 simpleFoam 时也遇到了同样的问题，请问您是如何解决的。谢谢！