CFD中文网

@hangsz 非结构有限体积也有高精度格式

@wwzhao 原来是这样，icem软件的六面体导出网格前都要转成非结构，应该就是为了适应fluent中的非结构计算。

偶然看到一系列的网格生成软件的网址：Mesh generator，很好奇NUMECA的HexPress 和snappyHexMesh都是 unstructured hexahedral（六面体） mesh，%(#ff0000)[六面体不都应该是结构网格吗？]

wellSource是常数吗？那mesh.V是什么参数？

@东岳 在 讨论 Fluent计算效率优于OpenFOAM 中说：

Ur

您好！之前测试的时候，壁面附近网格越密，网格高度变低、长宽比也变大，要求的时间步长deltaT就越小。不知道您说的Ur是什么定义？不知道在壁面附近有没有特殊处理，增加它的时间步长。

@cccrrryyy 是的

@宝丁 您好！有改过湍流模型吗？不同模型结果也是有差别的

@东岳 您好，移动了下就不行了。现在是 https://www.jianguoyun.com/p/DcjVq2YQ9O7FBxiv3MIB

@cccrrryyy fluent的网格边界层更密，为1微米，满足y+小于1。

@wwzhao 谢谢，库郎数中的速度指的是特征速度，就用了远场速度，不过用您说的最大速度来预估更合理。还看到东岳老师分享的链接，介绍三维非结构网格的库郎数计算，这个公式中的最小体积可以用画网格软件得到，但不知道通量该如何取值。

@cccrrryyy 算例远场速度是72m/s，调节网格后时间步长提高到了3×10^-8，使用pointwise监测最窄边的长度为2.8×10^-5，使用这些值计算的库朗数是0.077，用rhoCentralFoam计算过程中显示的最大库朗数确是0.53，有点不理解，不知道这个最大库朗数点的位置在哪。

Mean and max Courant Numbers = 0.001871585 0.5321452

@cccrrryyy 您好！请问你的时间步长具体值一般是多少？还有就是网格边界层高度是多高？我的分别是10e-8s和0.0001m。 算的确实慢，难以接受

@wwzhao @cccrrryyy 您好！我还处于入门阶段，一边用OpenFOAM，一边学OpenFOAMwiki上的3weeks 和编程查询，我的fvSchemes和fvSolution文件如下，目前调试经验是在fvSchemes中设置fluxScheme Tadmor;算得比设置Kurganov;要快。我后面再试下各位的设置，再过来交流。

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

fluxScheme          Tadmor;// Kurganov;

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    grad(U)         cellLimited Gauss linear 1;
    grad(nuTilda)   cellLimited Gauss linear 1;
    grad(e)   cellLimited Gauss linear 1;
}

divSchemes
{
    default         none;
    div(tauMC)      Gauss linear;
    div(phi,nuTilda) Gauss linearUpwind limited;
}



laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}


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

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-7;
        relTol          0.01;
        smoother        GaussSeidel;
        cacheAgglomeration off;
        nCellsInCoarsestLevel 20;
        agglomerator    faceAreaPair;
        mergeLevels     1;
    }

    Phi
    {
        $p;
    }


    pFinal
    {
        $p;
        tolerance       1e-06;
        relTol          0;
    }

    "(k|omega|epsilon|nuTilda|U|R|v2|f|kt|kl|e)"    
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        nSweeps         1;
        tolerance       1e-20;
        relTol          0.01;
    }


  "(rho|rhoU|rhoE)"
    {
        solver          diagonal;
    }


}

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

试了下同样的算例用Fluent 19计算，真的是黑匣子。。。时间步长0.2s都不发散。要是能借鉴一些方法到OpenFOAM里面就好了

@东岳 @wwzhao @cccrrryyy 您好！今天看到cfd-online上的链接：https://www.cfd-online.com/Forums/openfoam-solving/123586-simulations-large-time-steps-high-cfl.html 。里面讨论了 simulations with large time steps (high CFL)，帖子比较老了，说原因是OPENFOAM中没有像fluent那样Pressure-velocity coupling的coupled solver。粗略的看了下文献《A fully coupled OpenFOAM® solver for transient incompressible turbulent flows in ALE formulation》，里面说Implicit + coupling 效率要高。

很开心看到大家这么多的观点，下面也说说我最近的%(#ff0000)[个人]看法

• 不管fluent还是OpenFOAM都只是工具，更关键的是所关注的问题及其内在的科学规律、机理。资深级CFDer应该会有判断garbage result的能力，会对结果有个大体的预期，对意想不到的结果也有分析判断的能力，甚至有的可以根据流场看出用什么软件/模型求解的。fluent设计目的是提供一个a black box solver that is robust，OpenFOAM则可以自编代码，处于食物链的上一级，更高一级的就是全部自编的in-house代码了。希望以后可以看到fluent的UDF提供更多像OpenFOAM的可操作空间。

• fluent的非定常计算采用implicit，隐式的求解和时间无关，所以时间步长相比显式要大，不过隐式也存在迭代收敛的问题，这样的问题在OpenFOAM中也会出现，即not converged within 20 iterations。求解方法上用隐式或显式消耗的时间理论上无法比较，工程上一般隐式快一点。在fluent中用隐式进行定常计算，在计算稳定后，设置库郎数100也是可以的，因为定常的最终结果是一定的。但是非定常就不行了，这会牵涉到时间精度的问题。

• 作为新手，我学习OpenFOAM的过程是先将OpenFOAM用成fluent，然后再准备修改OpenFOAM得到属于自己的fluent。 不知道这样合不合理，主要是直接看代码，难以短期内有进展，很多情况下要求的是做错了也要有个错的结果。。。

• 回到标题，在计算三维LES时，相同网格在fluent里面计算的时间步长是10^-6，在OpenFOAM中用了runTimeModifiable 和 adjustTimeStep变为了10^-9~10^-8，不用就报错。不知道是不是就是大家所说的AMG求解器的原因。

@youmengtian 您好！我完整的格式是这样的

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

fluxScheme          Tadmor;// Kurganov;

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
         default cellMDLimited Gauss linear 1.0;
}

divSchemes
{
    default         none;
    div(tauMC)      Gauss linear;
    div(phi,nuTilda) Gauss linearUpwind limited;
}



laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}


snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}

@adcpk 在 关于rhoCentralFoam计算高超声速问题时温度出现负值的原因 中说：

divSchemes

您好！divSchemes用的是Euler

ddtSchemes
{
    default         Euler;
}

@adcpk 您好！插值格式用的是linear

interpolationSchemes
{
    default         linear;
}

@东岳 这是具体的算例文件 https://www.jianguoyun.com/p/Df7zvZYQ9O7FBxjyusIB 调试了rhoCentralFoam的设置，都计算不了。