CFD中文网

@李东岳 谢谢李老师，网格检查过了，没有问题，把cyclicAMI的边界设置为symmetry是可以计算的。刚才按照您的建议试过了，也还是报相同的错误，错误信息如下：

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2112                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _6e1fca0e-20220610 OPENFOAM=2112 patch=220610 version=2112
Arch   : "LSB;label=32;scalar=64"
Exec   : dbnsFoam
Date   : Jun 06 2024
Time   : 18:13:40
Host   : DESKTOP-H49IU2O
PID    : 137843
I/O    : uncollated
Case   : /home/duan/JIA_Xiaotong/ICSFoam-master/tutorials/test_5
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Steady-state analysis detected

pseudoTime: tolerance = 1e-10, relTol = 1e-10

Local timestepping selected

Reading thermophysical properties

Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

AMI: Creating addressing and weights between 18620 source faces and 18620 target faces
AMI: Patch source sum(weights) min:0 max:1.84903 average:0.999732
AMI: Patch target sum(weights) min:0.030709 max:1.15763 average:0.994612
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib/x86_64-linux-gnu/libc.so.6
#3  ? at ??:?
#4  ? at ??:?
#5  Foam::psiThermo::addfvMeshConstructorToTable<Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#6  ? at ??:?
#7  Foam::psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#8  ? in ~/OpenFOAM/duan-v2112/platforms/linux64GccDPInt32Opt/bin/dbnsFoam
#9  ? in /lib/x86_64-linux-gnu/libc.so.6
#10  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#11  ? in ~/OpenFOAM/duan-v2112/platforms/linux64GccDPInt32Opt/bin/dbnsFoam
Floating point exception

老师，有没有可能是因为叶栅模型的旋转半径太大（120m)导致的呢？

我在在设置旋转周期边界时遇到了以下问题：
我使用的是一个openfoam外开发的求解器，求解C3X叶栅的模拟，在使用相同的边界条件时，当设置叶片两个周期性边界为对称边界时，求解一切正常，但将其改为cyclicAMI时则计算无法进行，个人角度感觉是边界条件设置的问题，尝试过在·boundary文件中设置separationVector等，均会出下以下错误：

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Steady-state analysis detected

pseudoTime: tolerance = 1e-10, relTol = 1e-10

Local timestepping selected

Reading thermophysical properties

Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

AMI: Creating addressing and weights between 18620 source faces and 18620 target faces
AMI: Patch source sum(weights) min:0 max:1.84903 average:0.999732
AMI: Patch target sum(weights) min:0.030709 max:1.15763 average:0.994612
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib/x86_64-linux-gnu/libc.so.6
#3  ? at ??:?
#4  ? at ??:?
#5  Foam::psiThermo::addfvMeshConstructorToTable<Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#6  ? at ??:?
#7  Foam::psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#8  ? in ~/OpenFOAM/duan-v2112/platforms/linux64GccDPInt32Opt/bin/dbnsFoam
#9  ? in /lib/x86_64-linux-gnu/libc.so.6
#10  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#11  ? in ~/OpenFOAM/duan-v2112/platforms/linux64GccDPInt32Opt/bin/dbnsFoam
Floating point exception

网格的质量没有问题，boundary处设置的条件如下：

(
    blade
    {
        type            wall;
        inGroups        1(wall);
        nFaces          25088;
        startFace       6512669;
    }
    inlet
    {
        type            patch;
        nFaces          8526;
        startFace       6537757;
    }
    outlet
    {
        type            patch;
        nFaces          8526;
        startFace       6546283;
    }
    Lower_periodicity
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          18620;
        startFace       6554809;
	matchTolerance  0.0001;
        transform       rotational;
        neighbourPatch  Upper_periodicity;
        rotationAxis    (1 0 0);
        rotationCentre  (0 0 0);
    }
    Upper_periodicity
    {
        type            cyclicAMI;
        inGroups        1(cyclicAMI);
        nFaces          18620;
        startFace       6573429;
	matchTolerance  0.0001;
        transform       rotational;
        neighbourPatch  Lower_periodicity;
        rotationAxis    (1 0 0);
        rotationCentre  (0 0 0);
    }
    back
    {
        type            wall;
        inGroups        1(wall);
        nFaces          22363;
        startFace       6592049;
    }
    front
    {
        type            wall;
        inGroups        1(wall);
        nFaces          22363;
        startFace       6614412;
    }
)

同时也尝试过将cyclicAMI改为cyclic，但是两个边界的网格差异过大，有以下问题：

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0



--> FOAM FATAL ERROR: (openfoam-2112 patch=220610)
face 0 area does not match neighbour by 68.9656% -- possible face ordering problem.
patch:Lower_periodicity my area:6.66444e-09 neighbour area:1.36797e-08 matching tolerance:0.1
Mesh face:6554809 fc:(-0.0309194 -0.0589349 119.957)
Neighbour fc:(-0.0196249 0.058833 119.957)
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
Rerun with cyclic debug flag set for more information.

    From void Foam::cyclicPolyPatch::calcTransforms(const primitivePatch&, const pointField&, const vectorField&, const pointField&, const vectorField&)
    in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 215.

FOAM exiting

望大佬帮忙解答下困惑。

@李东岳 李老师，yang2009这篇sci论文里面我没看到有壁面剪应力的相关设置。由于这部分内容是我师兄做的，我先依据硕士论文解释您提出的几个问题，具体的我还需要再查证之后再给您回复。
硕士论文中引用的入流条件也是出自这篇sci文献，有出入的是壁面函数的粗糙高度取值，硕士论文中的粗糙高度是依据文献（Blocken et al. 2007, CFD simulation of the atmospheric boundary layer: wall function problems）计算所得：

壁面剪应力的取值是依据硕士论文中的一些假定提出的，而非yang2009中相关参数。

@李东岳 李老师您好，我们这边之前在fluent里面有验证过该入流条件，具体结果可参见刘安彬的硕士论文《平衡大气边界层的CFD数值模拟》。但是在模拟时上部流场的湍动能并无明显衰减，只是近壁面有衰减，并在施加壁面剪应力之后可有效改善。您可以考虑下顶面边界条件、壁面剪应力这两方面的设置。
下图模拟结果截取自上述硕士论文，未加任何修正

下图是施加壁面剪应力后的模拟结果

你是否比较Upwind格式下的残差曲线和linearUpwind 格式的差别呢？这两方法在OpenFOAM中的收敛性是不一样的。或者设置Porbes探针观察下数据采样点最后的波动情况。我觉的是收敛性影响的。有复杂地形或者建筑群的边界层计算时，OpenFOAM高阶阶方案的收敛性并不好，会影响计算精度@李东岳

请教各位老师两个问题（1）在采用有限体积法对方程进行离散之后，动量控制方程中左侧的惯性力项如何计算？
（2）在高湍流度的流动中，方程右侧的粘性力项很小，流体的受力状态如何理解？（左侧惯性力与右侧压力梯度的平衡吗？广义源项的加入如何影响流动的发展呢？） 谢谢~~~

@东岳 好的，我试一下。谢谢老师~

@东岳 嗯嗯，是在同一个文件夹0/下的

@东岳 李老师，我边界条件代码修改后如下：

    INLET
    {
        type            fixedProfile;
	profile 	tableFile;

	profileCoeffs
	{
            nHeaderLine         0;          // Number of header lines
            refColumn           0;          // Reference column index
            componentColumns    (1 2 3);    // Component column indices
            separator           ",";        // Optional (defaults to ",")
            mergeSeparators     no;         // Merge multiple separators
            file                "inVelocity.txt";//"inletProfiles/inVelocity.csv";
            outOfBounds         clamp;      // Optional out-of-bounds handling
            interpolationScheme linear;     // Optional interpolation scheme
        }
        direction        (0 1 0);
        origin           0;
     }

同时inVelocity.txt数据按照下面格式给出，但还是有同样的报错信息

2.31474,0,0
3.07832,0,0
3.44144,0,0
3.64768,0,0
3.78464,0,0
3.88646,0,0
3.96812,0,0
4.03721,0,0
4.09802,0,0
4.15319,0,0
4.20449,0,0
4.25314,0,0
4.30008,0,0
4.34601,0,0
4.39149,0,0
4.43702,0,0
4.483,0,0
4.5298,0,0

模拟二维大气边界层流动，入口添加ESDU的对数律剖面时，采用如下方式定义，

INLET
   {
       type            fixedProfile;
   profile 	tableFile;

   profileCoeffs
   {
           nHeaderLine         0;          // Number of header lines
           refColumn           1;          // Reference column index
           componentColumns    (1);    // Component column indices
           separator           " ";        // Optional (defaults to ",")
           mergeSeparators     no;         // Merge multiple separators
           file                "inVelocity.txt";//"inletProfiles/inVelocity.csv";
           outOfBounds         clamp;      // Optional out-of-bounds handling
           interpolationScheme linear;     // Optional interpolation scheme
       }
       direction        (0 1 0);
       origin           0;
    }

运行simpleFoam时，提示如下错误

inVelocity.txt文件为

0.899251
3.78464
3.88646
3.96812
4.03721
4.09802 
4.15319
4.20449
4.25314
4.30008 
4.34601 
4.39149	 
4.43702
4.483
4.5298
4.57773
4.62706
4.67807 
4.73099
4.78608 
4.84357
4.90368
4.96665
5.03268
5.10199	 
5.17476	 
5.25119	 
5.33143	 
5.41563	 
5.50394	 
5.59645	 
5.69327	 
5.79447	 
5.90012	 
6.01023

不知道问题在哪里？还请各位指教