[.bastian]

Hi.

I tried to set up an MRF case to be solved with simpleFoam based on the mixerVessel2D tutorial.
When i try to run simpleFoam the process quits after creating the mesh for time = 0 without an reasonable error.



Can anyone give me a hint what i'm doing wrong?


Here is the log.simpleFoam file:

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 6-47517f2ebb1b
Exec   : simpleFoam
Date   : Nov 25 2019
Time   : 16:34:20
Host   : "892836"
PID    : 254
I/O    : uncollated
Case   : /mnt/c/Users/72/Documents/lnxWork/TSOLA
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


SIMPLE: No convergence criteria found

Reading field p

Reading field U

Reading/calculating face flux field phi

AMI: Creating addressing and weights between 12536 source faces and 6271 target faces
AMI: Patch source sum(weights) min/max/average = 0.0919981, 1.0503, 0.982951
AMI: Patch target sum(weights) min/max/average = 0, 1.00164, 0.983866
AMI: Creating addressing and weights between 7424 source faces and 3337 target faces
AMI: Patch source sum(weights) min/max/average = 0.999738, 1.0005, 1.00006
AMI: Patch target sum(weights) min/max/average = 0.999067, 1.00017, 0.999992
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 0.075 average: 0
bounding epsilon, min: 0 max: 2 average: 2
RAS
{
    RASModel        kEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

Creating MRF zone list from MRFProperties
    creating MRF zone: MRF1
No finite volume options present

[Krao]

With the information provided by you, it is very difficult to add any input without seeing the case setup.

[.bastian]

Originally Posted by Krao

With the information provided by you, it is very difficult to add any input without seeing the case setup.

Mainly I used the tutorial "mixerVessel2D", which runs without any issues, for my setup. fvSolution, fvSchemes and controlDict are taken from it.
Below additional information to my case-setup.
Do you have any particular file in mind that might give the solution?

What surprises me the most is that the termination is done without error message.

system/controlDict:

/*--------------------------------*- 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      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     simpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         500;

deltaT          1;

writeControl    timeStep;

writeInterval   50;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;


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

system/fvSchemes:

/*--------------------------------*- 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;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      bounded Gauss limitedLinearV 1;
    div(phi,k)      bounded Gauss limitedLinear 1;
    div(phi,epsilon) bounded Gauss limitedLinear 25;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}


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

system/fvSolution:

/*--------------------------------*- 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-08;
        relTol          0.05;
        smoother        GaussSeidel;
        nCellsInCoarsestLevel 20;
    }

    U
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        nSweeps         2;
        tolerance       1e-07;
        relTol          0.1;
    }

    k
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        nSweeps         2;
        tolerance       1e-07;
        relTol          0.1;
    }

    epsilon
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        nSweeps         2;
        tolerance       1e-07;
        relTol          0.1;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 2;
  //  pRefCell        0;
  //  pRefValue       0;
}

relaxationFactors
{
    fields
    {
        p               0.3;
    }
    equations
    {
        U               0.5;
        k               0.5;
        epsilon         0.5;
    }
}


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

constant/transportProperties:

/*--------------------------------*- 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    "constant";
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

transportModel  Newtonian;

nu              [0 2 -1 0 0 0 0] 1e-05;

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

constant/MRFProperties:

/*--------------------------------*- 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    "constant";
    object      MRFProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

MRF1
{
    cellZone    rotor;
    active      yes;

    nonRotatingPatches ();

    origin    (0 0 0);
    axis      (0 0 1);
    omega     -753.98; //in rad/sec
}

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

0/U:

/*--------------------------------*- 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       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    rotorInlet
    {
        type            flowRateInletVelocity;
        massFlowRate    6.93;
        value           $internalField;
    }

    rotorOutlet
    {
        type            cyclicAMI;
        value           $internalField;
    }

    rotorWalls
    {
        type            noSlip;
    }

    statorInlet
    {
        type            cyclicAMI;
        value           $internalField;
    }
    statorOutlet
    {
        type            cyclicAMI;
        value           $internalField;
    }

    statorWalls
    {
        type            noSlip;
    }

    voluteInlet
    {
        type            cyclicAMI;
        value           $internalField;
    }
    voluteOutlet
    {
        type            pressureInletOutletVelocity;
        value           $internalField;
    }

    voluteWalls
    {
        type            noSlip;
    }

0/p:

/*--------------------------------*- 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       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    rotorInlet
    {
        type            zeroGradient;
    }

    rotorOutlet
    {
        type            cyclicAMI;
        value           $internalField;
    }

    rotorWalls
    {
        type            zeroGradient;
    }

    statorInlet
    {
        type            cyclicAMI;
        value           $internalField;
    }
    statorOutlet
    {
        type            cyclicAMI;
        value           $internalField;
    }

    statorWalls
    {
        type            zeroGradient;
    }

    voluteInlet
    {
        type            cyclicAMI;
        value           $internalField;
    }
    voluteOutlet
    {
        type            fixedValue;
        value           5e5;
    }

    voluteWalls
    {
        type            zeroGradient;
    }

0/nut:

/*--------------------------------*- 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       volScalarField;
    location    "0";
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
  rotorInlet
  {
    type            calculated;
    value           uniform 0;
  }

  rotorOutlet
  {
      type            cyclicAMI;
      value           $internalField;
  }

  rotorWalls
  {
    type            nutkWallFunction;
    value           uniform 0;
  }

  statorInlet
  {
      type            cyclicAMI;
      value           $internalField;
  }
  statorOutlet
  {
      type            cyclicAMI;
      value           $internalField;
  }

  statorWalls
  {
    type            nutkWallFunction;
    value           uniform 0;
  }

  voluteInlet
  {
      type            cyclicAMI;
      value           $internalField;
  }
  voluteOutlet
  {
      type            zeroGradient;
  }

  voluteWalls
  {
    type            nutkWallFunction;
    value           uniform 0;
  }

}
// ************************************************************************* //

0/k:

/*--------------------------------*- 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       volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
  rotorInlet
  {
    type            turbulentIntensityKineticEnergyInlet;
    intensity       0.05;
    value           uniform 0.075;
  }

  rotorOutlet
  {
      type            cyclicAMI;
      value           $internalField;
  }

  rotorWalls
  {
    type            kqRWallFunction;
    value           uniform 0;
  }

  statorInlet
  {
      type            cyclicAMI;
      value           $internalField;
  }
  statorOutlet
  {
      type            cyclicAMI;
      value           $internalField;
  }

  statorWalls
  {
    type            kqRWallFunction;
    value           uniform 0;
  }

  voluteInlet
  {
      type            cyclicAMI;
      value           $internalField;
  }
  voluteOutlet
  {
      type            zeroGradient;
  }

  voluteWalls
  {
    type            kqRWallFunction;
    value           uniform 0;
  }
}
// ************************************************************************* //

0/epsilon:

/*--------------------------------*- 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       volScalarField;
    location    "0";
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -3 0 0 0 0];

internalField   uniform 2;

boundaryField
{

    rotorInlet
    {
      type            fixedValue;
      value           uniform 0;
    }

    rotorOutlet
    {
        type            cyclicAMI;
        value           $internalField;
    }

    rotorWalls
    {
      type            epsilonWallFunction;
      value           uniform 0;
    }

    statorInlet
    {
        type            cyclicAMI;
        value           $internalField;
    }
    statorOutlet
    {
        type            cyclicAMI;
        value           $internalField;
    }

    statorWalls
    {
      type            epsilonWallFunction;
      value           uniform 0;
    }

    voluteInlet
    {
        type            cyclicAMI;
        value           $internalField;
    }
    voluteOutlet
    {
        type            zeroGradient;
    }

    voluteWalls
    {
      type            epsilonWallFunction;
      value           uniform 0;
    }
}
// ************************************************************************* //

[Krao]

Hi Bastian,

your setup looks ok, now I found out the problem,

AMI: Creating addressing and weights between 12536 source faces and 6271 target faces
AMI: Patch source sum(weights) min/max/average = 0.0919981, 1.0503, 0.982951
AMI: Patch target sum(weights) min/max/average = 0, 1.00164, 0.983866
AMI: Creating addressing and weights between 7424 source faces and 3337 target faces
AMI: Patch source sum(weights) min/max/average = 0.999738, 1.0005, 1.00006
AMI: Patch target sum(weights) min/max/average = 0.999067, 1.00017, 0.999992
Selecting incompressible transport model Newtonian

As it can be seen from the bold text, the problem lies with AMI patch weight, the patch weight is zero somewhere. That means there is no matching between the stationary cell zone and rotor cell zone. I think the mesh should be investigated thoroughly.

But I am not sure why OpenFOAM is not throwing out any error. Go through the following threads, to know more about 0 weight error. One half propeller simulation with MRFSimpleFoam crashing and Periodic cyclic AMI + MRF crashing on pressure

Best wishes,

K. Rao

[.bastian]

Hi Krao,



Thank you for your answer.



I did a complete rework of the case including a better mesh.
Although the weighting has improved, the result remains the same as you can see in the log file.

Any idea what to check next?

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 6-47517f2ebb1b
Exec   : simpleFoam
Date   : Nov 28 2019
Time   : 17:12:12
Host   : "892836"
PID    : 316
I/O    : uncollated
Case   : /mnt/c/Users/72/Documents/lnxWork/TSOLA
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


SIMPLE: No convergence criteria found

Reading field p

Reading field U

Reading/calculating face flux field phi

AMI: Creating addressing and weights between 9430 source faces and 4928 target faces
AMI: Patch source sum(weights) min/max/average = 0.999215, 1.02469, 0.999989
AMI: Patch target sum(weights) min/max/average = 0.990802, 1.00066, 1.00006
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 0.075 average: 0
bounding epsilon, min: 0 max: 2 average: 2
RAS
{
    RASModel        kEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

Creating MRF zone list from MRFProperties
    creating MRF zone: MRF1
No finite volume options present

[Krao]

Hi Bastian,

difficult to add anything, the setup looks decent. Leave log file, before stopping does the terminal shows anything or is it also blank? If I get time on weekend I can go through your setup. If you don't want to share geometry, you can remove the geometry and important info from mesh and you can send the case to me. I can see if I could do something. Or if you think it is not a big secret then you can post it here and anyone can help you.

K. Rao

[.bastian]

Hello Krao, sorry I'm late back.

It seems that the problem was due to a too big imbalance between the number of source faces to target faces.
After I adjusted the mesh again, it worked.

Nevertheless, I have switched to foam-extend. The GGI interfaces have more options and look much more stable and evolved than AMI.

[Krao]

Hi Bastian,

After receiving your feedback about foam-extend, I have though to make a comparison of results obtained using openFoam and foam-extend. I tried MRFsimpleFoam with ggi instead of AMI. I am facing problems during parallel running. My case runs very well in serial, when I decompose and run the case in parallel, my simulation stop after reading 'Reading/calculating face flux field phi'. Did you face similar problem by any chance? Which type of decomposition method are you using? It would be nice to know.

Thank you,
K. Rao

[.bastian]

Originally Posted by Krao

[...] Which type of decomposition method are you using? [...]

Hi Krao,

generally, in the decomposeParDict, you have to specify all interfaces (ggi/mixingPlane) as globalFaceZones.

Additionally, at least for mixing planes, you have to do the partitioning patchConstrained. This means that the interfaces are assigned to a single processor instead of getting "cutted".
Like i mentioned before, i'm not sure if its a only necessary in case of mixing planes or also affects ggi interfaces. Anyway it works for both.

See below my decomposePar dict. Just adapt the names of the interfaces to your case. The numbers after the names in the patchConstrainedCoeffs sets the number of the processor to which they will be assigned to.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.1                                |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      decomposeParDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

nCores    12;

numberOfSubdomains $nCores;

method    patchConstrained;

globalFaceZones
(
    rotorOutletZone
    rotorPeriodicPSZone
    rotorPeriodicSSZone
    statorInletZone
);

patchConstrainedCoeffs
{
    method                         metis;
    numberOfSubdomains    $nCores;
    patchConstraints
    (
        (rotorOutletMP 1)
        (statorInletMP 1)
        (rotorPeriodic_PS 2)
        (rotorPeriodic_SS 2)
    );
}


metisCoeffs
{
    processorWeights
    (
        1
        1
        1
        1
        1
        1
        1
        1
        1
        1
        1
        1
    );
}
// ************************************************************************* //

.bastian