[jyotir]

Hi,

I am facing difficulties in simulating an axial compressor using 'rhoPimpleFoam' in OpenFOAM v9. I am a beginner in OpenFOAM.
I am not able to get any meaningful solution. I am pasting the mass flow time history I am getting from the simulation.

I think
1) MRF is not working for my case as I see hardly any change in the solution if I set 'omega' to zero.
2) My solver/schemes settings are probably incorrect

I am including here my case files. I will be grateful for any help in this regard.

I am including my directory structure also.

Thank you,
Jyoti

log (truncated)

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 9
Exec   : rhoPimpleFoam -parallel
Date   : Nov 17 2021
Time   : 13:43:00
Host   : "aselab-server"
PID    : 202923
I/O    : uncollated
Case   : /home/aselab/OpenFOAM/aselab-9/run/axcomp
nProcs : 12
Slaves : 
11
(
"aselab-server.202926"
"aselab-server.202927"
"aselab-server.202928"
"aselab-server.202929"
"aselab-server.202930"
"aselab-server.202931"
"aselab-server.202932"
"aselab-server.202933"
"aselab-server.202934"
"aselab-server.202935"
"aselab-server.202936"
)

Pstream initialised with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
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


Reading thermophysical properties

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

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
    model           kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}

Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity
Creating field dpdt

Creating field kinetic energy K

Creating MRF zone list from MRFProperties
    creating MRF zone: MRF1
No fvModels present
No fvConstraints present
Courant Number mean: 0.000586994 max: 0.725888
surfaceFieldValue patchFlowRateOutlet:
    total faces  = 41184
    total area   = 0.133121


surfaceFieldValue patchFlowRateInlet:
    total faces  = 41184
    total area   = 0.182797



Starting time loop

Courant Number mean: 0.000586994 max: 0.725888
deltaT = 1.13766e-08
Time = 1.13766e-08

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 7.76942e-06, No Iterations 15
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 7.67355e-06, No Iterations 15
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 7.12107e-06, No Iterations 15
smoothSolver:  Solving for h, Initial residual = 1, Final residual = 7.73572e-06, No Iterations 3
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.017194, No Iterations 1
GAMG:  Solving for p, Initial residual = 2.5631e-05, Final residual = 9.98724e-07, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.06168e-08, global = 1.83108e-09, cumulative = 1.83108e-09
GAMG:  Solving for p, Initial residual = 0.00035044, Final residual = 5.05704e-06, No Iterations 1
GAMG:  Solving for p, Initial residual = 5.62303e-06, Final residual = 1.67461e-07, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.53955e-09, global = -2.29958e-10, cumulative = 1.60112e-09
smoothSolver:  Solving for omega, Initial residual = 0.000548122, Final residual = 4.7661e-06, No Iterations 3
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 5.81245e-06, No Iterations 7
ExecutionTime = 12.89 s  ClockTime = 13 s

Courant Number mean: 0.00279069 max: 2.06948
deltaT = 5.4972e-09
Time = 1.68738e-08

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.0415336, Final residual = 4.11692e-07, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0429697, Final residual = 4.24706e-07, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.000414766, Final residual = 6.38689e-07, No Iterations 1
smoothSolver:  Solving for h, Initial residual = 0.0659807, Final residual = 7.81129e-07, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0009249, Final residual = 2.6245e-06, No Iterations 1
GAMG:  Solving for p, Initial residual = 1.82371e-06, Final residual = 9.40284e-09, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 9.82654e-11, global = -7.90478e-12, cumulative = 1.59322e-09
GAMG:  Solving for p, Initial residual = 2.2142e-05, Final residual = 1.64108e-07, No Iterations 1
GAMG:  Solving for p, Initial residual = 1.65298e-07, Final residual = 1.65298e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.71618e-09, global = -5.57808e-10, cumulative = 1.03541e-09
smoothSolver:  Solving for omega, Initial residual = 0.000147196, Final residual = 1.22326e-07, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.0037913, Final residual = 4.34495e-07, No Iterations 1
ExecutionTime = 17.27 s  ClockTime = 17 s

Courant Number mean: 0.0013461 max: 1.06947
deltaT = 5.13999e-09
Time = 2.20138e-08

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.0211019, Final residual = 1.46182e-07, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.02252, Final residual = 1.55663e-07, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.000233115, Final residual = 2.6705e-07, No Iterations 1
smoothSolver:  Solving for h, Initial residual = 0.00219986, Final residual = 2.08177e-07, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.000658396, Final residual = 1.14394e-06, No Iterations 1
GAMG:  Solving for p, Initial residual = 9.70627e-07, Final residual = 9.70627e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.19003e-08, global = -1.46945e-10, cumulative = 8.88463e-10
GAMG:  Solving for p, Initial residual = 1.15099e-05, Final residual = 5.62891e-08, No Iterations 1
GAMG:  Solving for p, Initial residual = 5.65247e-08, Final residual = 5.65247e-08, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.9081e-10, global = -1.58205e-10, cumulative = 7.30258e-10
smoothSolver:  Solving for omega, Initial residual = 4.59091e-05, Final residual = 3.17886e-08, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.00571728, Final residual = 3.07121e-07, No Iterations 1
bounding k, min: -0.000354255 max: 84.7546 average: 0.0585542
ExecutionTime = 21.59 s  ClockTime = 22 s

.
.
.
. TRUNCATED
.
.
.

Courant Number mean: 0.00103208 max: 0.999906
deltaT = 4.49553e-09
Time = 5.58449e-05

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 5.7813e-05, Final residual = 1.31039e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 5.73147e-05, Final residual = 1.25815e-09, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1.37477e-05, Final residual = 3.00424e-10, No Iterations 1
smoothSolver:  Solving for h, Initial residual = 4.72346e-05, Final residual = 1.11587e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 6.23986e-05, Final residual = 1.76731e-11, No Iterations 1
GAMG:  Solving for p, Initial residual = 2.82336e-11, Final residual = 2.82336e-11, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.71983e-12, global = -1.50626e-14, cumulative = 3.02554e-08
GAMG:  Solving for p, Initial residual = 6.71047e-09, Final residual = 6.71047e-09, No Iterations 0
GAMG:  Solving for p, Initial residual = 6.71047e-09, Final residual = 6.71047e-09, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.84106e-10, global = 4.33091e-13, cumulative = 3.02559e-08
smoothSolver:  Solving for omega, Initial residual = 5.73841e-06, Final residual = 5.73841e-06, No Iterations 0
smoothSolver:  Solving for k, Initial residual = 9.98794e-06, Final residual = 9.98794e-06, No Iterations 0
ExecutionTime = 45920.2 s  ClockTime = 45922 s

Courant Number mean: 0.00103208 max: 0.999909
deltaT = 4.49553e-09
Time = 5.58494e-05

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 5.78063e-05, Final residual = 1.31022e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 5.73081e-05, Final residual = 1.258e-09, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1.37471e-05, Final residual = 3.00413e-10, No Iterations 1
smoothSolver:  Solving for h, Initial residual = 4.72299e-05, Final residual = 1.11579e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 6.24109e-05, Final residual = 1.76701e-11, No Iterations 1
GAMG:  Solving for p, Initial residual = 2.82316e-11, Final residual = 2.82316e-11, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.71905e-12, global = -1.48986e-14, cumulative = 3.02558e-08
GAMG:  Solving for p, Initial residual = 6.714e-09, Final residual = 6.714e-09, No Iterations 0
GAMG:  Solving for p, Initial residual = 6.714e-09, Final residual = 6.714e-09, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.84448e-10, global = 4.3293e-13, cumulative = 3.02563e-08
smoothSolver:  Solving for omega, Initial residual = 5.73839e-06, Final residual = 5.73839e-06, No Iterations 0
smoothSolver:  Solving for k, Initial residual = 9.98862e-06, Final residual = 9.98862e-06, No Iterations 0
ExecutionTime = 45924.1 s  ClockTime = 45926 s

Courant Number mean: 0.00103208 max: 0.999911
deltaT = 4.49553e-09
Time = 5.58539e-05

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 5.77997e-05, Final residual = 1.31006e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 5.73015e-05, Final residual = 1.25785e-09, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1.37464e-05, Final residual = 3.00399e-10, No Iterations 1
smoothSolver:  Solving for h, Initial residual = 4.72252e-05, Final residual = 1.11571e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 6.24232e-05, Final residual = 1.76675e-11, No Iterations 1
GAMG:  Solving for p, Initial residual = 2.82297e-11, Final residual = 2.82297e-11, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.7183e-12, global = -1.4714e-14, cumulative = 3.02563e-08
GAMG:  Solving for p, Initial residual = 6.72012e-09, Final residual = 6.72012e-09, No Iterations 0
GAMG:  Solving for p, Initial residual = 6.72012e-09, Final residual = 6.72012e-09, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.85133e-10, global = 4.3279e-13, cumulative = 3.02567e-08
smoothSolver:  Solving for omega, Initial residual = 5.73838e-06, Final residual = 5.73838e-06, No Iterations 0
smoothSolver:  Solving for k, Initial residual = 9.98931e-06, Final residual = 9.98931e-06, No Iterations 0
ExecutionTime = 45928 s  ClockTime = 45930 s

Courant Number mean: 0.00103208 max: 0.999914
deltaT = 4.49553e-09
Time = 5.58584e-05

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 5.7793e-05, Final residual = 1.30989e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 5.72949e-05, Final residual = 1.2577e-09, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1.37458e-05, Final residual = 3.00387e-10, No Iterations 1
smoothSolver:  Solving for h, Initial residual = 4.72205e-05, Final residual = 1.11558e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 6.24355e-05, Final residual = 1.76647e-11, No Iterations 1
GAMG:  Solving for p, Initial residual = 2.82278e-11, Final residual = 2.82278e-11, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.71753e-12, global = -1.45445e-14, cumulative = 3.02567e-08
GAMG:  Solving for p, Initial residual = 6.72174e-09, Final residual = 6.72174e-09, No Iterations 0
GAMG:  Solving for p, Initial residual = 6.72174e-09, Final residual = 6.72174e-09, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.85224e-10, global = 4.32636e-13, cumulative = 3.02571e-08
smoothSolver:  Solving for omega, Initial residual = 5.73836e-06, Final residual = 5.73836e-06, No Iterations 0
smoothSolver:  Solving for k, Initial residual = 9.99e-06, Final residual = 9.99e-06, No Iterations 0
ExecutionTime = 45931.8 s  ClockTime = 45934 s

Courant Number mean: 0.00103209 max: 0.999917
deltaT = 4.49553e-09
Time = 5.58629e-05

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 5.77864e-05, Final residual = 1.30973e-09, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 5.72883e-05, Final residual = 1.25755e-09, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1.37451e-05, Final residual = 3.00375e-10, No Iterations 1
smoothSolver:  Solving for h, Initial residual = 4.72158e-05, Final residual = 1.11547e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 6.24477e-05, Final residual = 1.76617e-11, No Iterations 1
GAMG:  Solving for p, Initial residual = 2.82259e-11, Final residual = 2.82259e-11, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.71676e-12, global = -1.44152e-14, cumulative = 3.02571e-08
GAMG:  Solving for p, Initial residual = 6.72583e-09, Final residual = 6.72583e-09, No Iterations 0
GAMG:  Solving for p, Initial residual = 6.72583e-09, Final residual = 6.72583e-09, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.8564e-10, global = 4.32441e-13, cumulative = 3.02575e-08