probes和sample提取wallShearStress问题

李东岳

如果找不到湍流模型，如要挂在求解器上，求解器上有turbulence的库：driftFluxFoam -postProcess -func ... 有几点：

首先多相流领域，湍流模型研究相对不是热点。因为多相流本身就模化的很厉害，很难做到单向NS那种高解析的结果。因此多相流领域湍流模型的研究，再怎么精准也提高不到单向流NS方程那种。当然了肯定有多相流湍流研究的文章，但是看他们的数据的话，偏差30%都有。不过我认为也是可以理解的。
在这个背景下，用driftFluxFoam来算流场，然后依据流场结果+湍流模型来计算壁面剪切力，这里面模化的因素忒大了。就单单说泥沙颗粒一层在壁面网格的处理，有搞出好多种方法。

coolhhh

测试了一下，driftFluxFoam 好像无法直接用wallShearStress这个函数，还是会提示找不到模型

driftFluxFoam -postProcess -func wallShearStress -latestTime

李东岳

@coolhhh 没开虚拟机。log怎么提示。用的什么版本的openfoam

coolhhh

@李东岳
李老师，我测试了下of8和of2206两个版本：
测试算例为: OpenFOAM-8/tutorials/multiphase/driftFluxFoam/RAS/mixerVessel2D
运行完算例后都执行命令：driftFluxFoam -postProcess -func wallShearStress -latestTime

of8版本结果
of2206版本结果
有没有可能是这个算例用的湍流模型 buoyantKEpsilon 有关

RASModel        buoyantKEpsilon;

李东岳

如果直接用postProcess是会出现没有这个专属的湍流模型，一会我试一下

李东岳

我试了一下OpenFOAM-10，这个可以用

dyfluid@dyfluid-virtual-machine:~/tutorials/tutorials-10/multiphase/driftFluxFoam/RAS/mixerVessel2DMRF$ ./Allrun 
Running blockMesh on /home/dyfluid/tutorials/tutorials-10/multiphase/driftFluxFoam/RAS/mixerVessel2DMRF
Running driftFluxFoam on /home/dyfluid/tutorials/tutorials-10/multiphase/driftFluxFoam/RAS/mixerVessel2DMRF
^C
dyfluid@dyfluid-virtual-machine:~/tutorials/tutorials-10/multiphase/driftFluxFoam/RAS/mixerVessel2DMRF$ driftFluxFoam -postProcess -func wallShearStress
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 10-03dd22794e29
Exec   : driftFluxFoam -postProcess -func wallShearStress
Date   : Aug 14 2024
Time   : 15:27:16
Host   : "dyfluid-virtual-machine"
PID    : 96237
I/O    : uncollated
Case   : /home/dyfluid/tutorials/tutorials-10/multiphase/driftFluxFoam/RAS/mixerVessel2DMRF
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: No convergence criteria found


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode


Time = 0s
Reading field p_rgh

Reading field U

Reading/calculating face flux field phi


Reading g

Reading hRef
Calculating field g.h

Reading incompressibleTwoPhaseInteractingMixture

Selecting mixture viscosity model BinghamPlastic
Selecting viscosity model constant
Creating MRF zone list from MRFProperties
    creating MRF zone: MRF1
Selecting relative velocity model simple
Selecting turbulence model type RAS
Selecting RAS turbulence model buoyantKEpsilon
buoyantKEpsilonCoeffs
{
    Cg              0.85;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

No fvModels present
No fvConstraints present
wallShearStress wallShearStress:
    processing all wall patches

wallShearStress wallShearStress write:
    writing object wallShearStress
    min/max(rotor) = (-0.493384 -0.493384 -2.2956e-20), (0.493384 0.493384 3.1335e-20)
    min/max(stator) = (-0 -0 -0), (-0 -0 -0)


Time = 0.25s
Reading field p_rgh

Reading field U

Reading/calculating face flux field phi


Reading g

Reading hRef
Calculating field g.h

Reading incompressibleTwoPhaseInteractingMixture

Selecting mixture viscosity model BinghamPlastic
Selecting viscosity model constant
Creating MRF zone list from MRFProperties
    creating MRF zone: MRF1
Selecting relative velocity model simple
Selecting turbulence model type RAS
Selecting RAS turbulence model buoyantKEpsilon
buoyantKEpsilonCoeffs
{
    Cg              0.85;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

No fvModels present
No fvConstraints present
wallShearStress wallShearStress:
    processing all wall patches

wallShearStress wallShearStress write:
    writing object wallShearStress
    min/max(rotor) = (-0.52359 -0.524005 -4.956e-19), (0.524154 0.523798 4.51473e-19)
    min/max(stator) = (-0.0495746 -0.0495683 -4.97663e-20), (0.0496049 0.0495712 1.89403e-19)


Time = 0.5s
Reading field p_rgh

Reading field U

Reading/calculating face flux field phi


Reading g

Reading hRef
Calculating field g.h

Reading incompressibleTwoPhaseInteractingMixture

Selecting mixture viscosity model BinghamPlastic
Selecting viscosity model constant
Creating MRF zone list from MRFProperties
    creating MRF zone: MRF1
Selecting relative velocity model simple
Selecting turbulence model type RAS
Selecting RAS turbulence model buoyantKEpsilon
buoyantKEpsilonCoeffs
{
    Cg              0.85;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

No fvModels present
No fvConstraints present
wallShearStress wallShearStress:
    processing all wall patches

wallShearStress wallShearStress write:
    writing object wallShearStress
    min/max(rotor) = (-0.573367 -0.573406 -9.16331e-19), (0.573481 0.573189 5.36707e-19)
    min/max(stator) = (-0.0448395 -0.0451926 -5.10002e-20), (0.0447064 0.0447249 1.00214e-19)


End

dyfluid@dyfluid-virtual-machine:~/tutorials/tutorials-10/multiphase/driftFluxFoam/RAS/mixerVessel2DMRF$

coolhhh

再试了下，of9、of10都可以计算，但of8，of2206，of2212，of2306都不行。估计com版本这个问题都没修复

李东岳

@coolhhh

吴建民

@coolhhh 非常感谢老师，我今天装了个of8，完美解决了

吴建民

@吴建民在 probes和sample提取wallShearStress问题中说：

@coolhhh 非常感谢老师，我今天装了个of8，完美解决了

说错了，是of9

CFD中文网

probes和sample提取wallShearStress问题