CFD中文网

@samuel-tu
不知道现在还有没有人有这个问题，我之前也遇到过。题主应该用的是com版本的OpenFOAM（v2006, v2112这种），经过ThirdParty编译安装paraFoam, 不支持python shell，无法使用trace录制python宏，但可以在终端使用pvpython。这个是因为默认编译选项中的withPYTHON是关闭的，需要在编译的时候打开withPYTHON选项。
如果不知道咋弄，直接单独装paraview，或者用org版本的OpenFOAM，自动安装paraFoam不用操心。

@aimer
好像，tecplot不只是surfaceField读不了（tecplot2019可以读，但是值都为零），tensorField也读不了，所以壁面剪切应力，速度梯度这些都导不进去。然而tecplot求的速度梯度在壁面处又是错的。感觉还是OpenFOAM postProcess + paraview比较靠谱
不知道有没有解决方案

@李东岳
李老师好！抱歉回复晚了。
算例文件在这，https://www.jianguoyun.com/p/DejQRSIQ-4XWCRiVid0EIAA，稍微有些大，压缩包16.8MB

各位前辈好，一个simpleFoam算低雷诺数（re≈0.1）绕多圆柱的算例，在网格变化稍大的地方出现明显的压力场不光滑的情况，像马赛克一样。网格正交性还不错，也进行了非正交修正（面梯度离散格式corrected，nNonOrthogonalCorrectors=2或4）。这是正常现象吗？

@波流力 网格差不多长这样

@cccrrryyy 感谢感谢，我在正式计算之前，算了单个圆柱的，和Oseen公式推出来的理论解能够比上。
您说的方法我试试，由于OpenFoam我也只是初学，没有特别系统的学过，您说的relaxationFactors 我还不太清楚，我去好好了解一下。

计算一个简单二维多圆柱绕流的动网格算例，网格数量大概40-50万，fvsolution中的nOuterCorrectors和nCorrecters我尝试了不同的值，发现随着该值设置的越来越大，结果才会收敛，最终分别等于20和10。但是这个条件下，计算时间特别久，一个二维40-50万网格的算例要计算将近90个小时，感觉很不合理。想请教一下各位大佬，我是不是哪里设置出错了？

/*--------------------------------*- 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       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    "(p|pcorr)"
    {
        solver          PCG;
        preconditioner        DIC;
        tolerance       1e-04;
        relTol          0.1;
    }
    "(p|pcorr)Final"
    {
        $p;
        relTol          0;
    }
    U
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-04;
        relTol          0.1;
    }
    UFinal
    {
        $U;
        tolerance       1e-04;
        relTol          0.0;
    }
}


PIMPLE
{
    momentumPredictor   no;
    correctPhi          yes;
    nOuterCorrectors    20;
    nCorrectors         10;
    nNonOrthogonalCorrectors 1;
    ddtCorr             true;

    pRefPoint           (20.0 0.0 0.5);
    pRefValue           0.0;
}

relaxationFactors
{
   fields
   {
   }
   equations
   {
       "U.*"             1;
   }
}
// ************************************************************************* //

/*--------------------------------*- 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       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(U)          Gauss linear;
    div(phi,U)      Gauss linearUpwind grad(U);
    div(phi,k)      Gauss linearUpwind grad(U);
    div(phi,K)      Gauss linearUpwind grad(U);
    div(phi,omega)  Gauss linearUpwind grad(U);
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}


// ************************************************************************* //

感谢感谢

楼主搞定了吗，我最近分析一个圆柱加速运动，也想计算一下其附加质量力

@队长别开枪 您好，我想请问一下，pointwise导入openfoam时，如何设置网格各个面的boundary类型呢

请问一下，有哪位大佬用过datToFoam这个命令吗，*.dat文件需要写成什么格式呢？

运用pimpleFoam，求解一个平板的运动问题，想让整个网格一起运动，程序顺利运行了，但是结果明显是错误。猜测是边界条件设置的问题。想请教一下各位大佬，这种情况下，边界条件如何设置？或者是否有其他方面的问题。

补充记录一下今天的发现，只有在合并压力p的时候会出现段错误，合并其他物理量的时候则不会。

我在尝试使用foam-extend4.0的浸入边界法时，并行也总是出错。后来改用foam-extend4.1的浸入边界法pimpleDyMIbFoam，分块并行的时候没问题，reconstructPar的时候出现段错误

利用foam-extend4.1中的pimpleDyMIbFoam求解器求解一个浸入边界法的问题，物面在网格中做拍打运动。单块网格求解没有报错，但是分块并行之后，在运行reconstructPar合并网格的时候，出现如下错误：
本人初学OpenFOAM，不是很熟悉，这个问题纠结了好久了，拜托各位大佬帮帮忙！

找到了解决办法了，使用命令surfaceConvert geometry.stl geometry.ftr可以将.stl转换为.ftr

您好，请问一下这个问题您找到解决办法了吗

用重叠网格求解器，overPimpleDyMFoam计算了一个平板前移，在物面网格和背景网格之间有非常明显的界线。想请教各位大佬，生成网格的时候，内外层网格之间需要注意什么才可以让两层网格间的界限不那么明显。

@bestucan 太感谢了

谢谢老师，不过这是在Windows10系统下运行paraview时出现的情况。