CFD中文网

青

@李东岳上面的blcokMesh确实可以，但是求解器不行，后面发现是系统没装gmp插件，有些内容编译不了，安装之后就可以了

青

@Amadeus 找到问题了，系统没装gmp插件，有些内容编译不了，安装之后就可以了

青

gcc版本也换过，结果都是一样

zf@zf-Lenovo-XiaoXinAir-14-ACN-2021:~/OpenFOAM/OpenFOAM-v2406$ blockMesh 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2406                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _630d60de3b-20240620 OPENFOAM=2406 version=v2406
Arch   : "LSB;label=32;scalar=64"
Exec   : blockMesh
Date   : Dec 25 2024
Time   : 18:21:00
Host   : zf-Lenovo-XiaoXinAir-14-ACN-2021
PID    : 372228
I/O    : uncollated
Case   : /home/zf/OpenFOAM/OpenFOAM-v2406
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



--> FOAM FATAL ERROR: (openfoam-2406)
cannot find file "/home/zf/OpenFOAM/OpenFOAM-v2406/system/controlDict"

    From virtual Foam::autoPtr<Foam::ISstream> Foam::fileOperations::uncollatedFileOperation::readStream(Foam::regIOobject&, const Foam::fileName&, const Foam::word&, bool) const
    in file global/fileOperations/uncollatedFileOperation/uncollatedFileOperation.C at line 629.

FOAM exiting

zf@zf-Lenovo-XiaoXinAir-14-ACN-2021:~/OpenFOAM/OpenFOAM-v2406$ icoFoam
找不到命令 “icoFoam”，但可以通过以下软件包安装它：
sudo apt install openfoam

青

最近在一台Ubuntu 22.04机子上想装一个openfoam-v2406，编译完之后也没报错，且blockMesh可以正常调用，但是求解器一个都不能用，看了一下编译后的platforms文件夹，bin里面一个求解器都没有，这是什么情况，重新编译了两次了，都是这个问题。以前在很多台电脑都装过，但是第一次遇到这种问题，望大神指教！

青

@派大星huangxiaojie

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/

forces1-gas
{
    type          forces;

    libs          ("libforces.so");

    writeControl  timeStep;
    timeInterval  1;

    log           yes;
    patches     (impeller_66 impeller_67);
    rho           rhoInf;     // Indicates incompressible
    log           true;
    rhoInf        700;          // Redundant for incompressible
    phase    gas;
     p p;
    CofR          (0 0.35 0);    // Rotation around centre line of propeller
    pitchAxis     (0 5 0);
}

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

把每个相的扭矩值都分别处理，然后相加即可，跟实验值较为吻合

青

@李东岳好的，谢谢老师

青

@李东岳另外老师还想请教个问题，scalarTransportFoam中库朗数还需要小于1吗，还是说时间步长可以大一些让库朗数大一些

青

@李东岳我的意思是把multiphaseEulerFoam算出来的U.particles和phi.particles一起复制过来作为scalarTransportFoam的初值，然后把原先phi.particles中的边界条件修改一下就可以了

青

@李东岳老师，问题解决了，要把固体的通量场也作为初值给scalarTransportFoam，并且要把phi里面得fixedValue类型改为calculated就可以了，谢谢老师

青

@李东岳或者说是不是多相流中的单相速度场并不能直接拿到scalarTransportFoam中使用，因为单相速度场没有相含率信息，例如一个网格的固含率为0，但是其还是会有一个固相速度，然而scalarTransportFoam仅读取速度场无法知道这一点，会不会导致质量不守恒的问题。

青

@李东岳老师，我继续计算了4个小时也还是没有变化，出口一直没有浓度，而且T的迭代次数已经都是0了，就是说已经稳态了，可是这样如何得到停留时间分布？

青

青

@李东岳老师，因为我想固体颗粒在液体中应该没有扩散，所以我只考虑了对流，扩散系数设了0，是这个问题吗

青

@李东岳进口固定值1，其他wall和出口都是零梯度

青

并且算着算着迭代次数逐渐变为0了

青

首先使用multiphaseEulerFoam模拟了一个气液固三相连续搅拌釜，想获得固体相在釜中的停留时间分布，所以想将固体的速度场和通量场给到scalarTransportFoam中，从入口进入示踪剂以获得固体停留时间分布。这时遇到一个问题，整个釜的空时大概4个小时，但是计算大概1000s之后整个釜内的示踪剂浓度分布就不变化了，出口的浓度也达到稳定并且极低，请问这是什么问题？附上我的设置

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,T)     Gauss linearUpwind grad(T);
}

laplacianSchemes
{
    default         none;
    laplacian(DT,T) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}


solvers
{
    T
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-6;
        relTol          0;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 1;
}
relaxationFactors
{

    equations
    {
        ".*"         1;
    }
}

青

@tidedrinker 你好，我发现of9中可以在force文件中单独定义相名称来提取相分数和对应的密度，输出的力和力矩如下：

forces forces write:
    sum of forces:
        pressure : (-10.813 0 -5.80003)
        viscous  : (-0.000457782 -0.000707157 0.000145593)
        porous   : (0 0 0)
    sum of moments:
        pressure : (-0.00441956 4.53963 -0.00703327)
        viscous  : (-3.21179e-05 7.86609e-05 -2.42533e-05)
        porous   : (0 0 0)

我的搅拌轴是y轴，所以我的扭矩是否就是4.53963N·M，但是这和实验值相差甚远，实验值大概只有0.06N·M，这是为什么？是不是还需要对这个力矩进行下一步处理才可以得到扭矩？需要除以密度嘛？

青

@tidedrinker 但是还有个问题，我在已经有了计算结果之后使用这个force进行后处理 -postProcess时，会报错说没有找到U和p，并不会自动读取U.*和p_rgh，请问您是怎么在多相流中使用的？

Reading field p_rgh

--> FOAM Warning :
    From function int main(int, char**)
    in file ./createFields.H at line 34
    Reading "/mnt/e/Desktop/limited-zeroG/constant/phaseProperties" from line 16 to line 144
    Pressure limits, pMin and pMax, are now read from "/mnt/e/Desktop/limited-zeroG/system/fvSolution/PIMPLE"
forces forces:
    Not including porosity effects


--> FOAM FATAL ERROR:
Could not find U, p

    From function void Foam::functionObjects::forces::initialise()
    in file forces/forces.C at line 205.

FOAM exiting

青

@tidedrinker 感谢，我现在在用of9的multiphaseEulerFoam，也就是说rho和rhoinf随便填一个数都没有影响是吗

青

各位大佬好，最近在做多相搅拌釜，想要提取桨叶的扭矩来计算搅拌功率。想要使用function中的force功能来提取桨叶上的扭矩，但是发现这个功能好像只适用于单相流。因为force文件需要输入rho，但是在多相流中由于相分布并不均匀且在瞬态变化，无法获得桨叶附近的rho值，所以该方法是否无法计算多相流的扭矩？

forces
{
    type          forces;

    libs          ("libforces.so");

    writeControl  timeStep;
    timeInterval  1;

    log           yes;

    patches       ("impeller*");
    rho           rhoInf;     // Indicates incompressible
    log           true;
    rhoInf        1;          // Redundant for incompressible

    CofR          (0 0 0);    // Rotation around centre line of propeller
    pitchAxis     (0 1 0);
}

另外拜读李东岳老师的博士论文发现李老师在论文中提取出了液-液两相流中的搅拌扭矩，说明多相流的扭矩是可以提取的，但是具体如何实现在论文中并未赘述，能否有大神帮忙解答，感谢！

CFD中文网

青山

帖子