CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Help sought on axial compressor simulation

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 17, 2021, 11:49
Default Help sought on axial compressor simulation
  #1
New Member
 
Jyoti Ranjan
Join Date: Oct 2020
Posts: 14
Rep Power: 6
jyotir is on a distinguished road
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)
Code:
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  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
Attached Images
File Type: png massflow.png (19.0 KB, 15 views)
File Type: png directory_structure.png (11.5 KB, 17 views)
Attached Files
File Type: gz 0.tar.gz (1.3 KB, 14 views)
File Type: gz constant.tar.gz (1.4 KB, 11 views)
File Type: gz system.tar.gz (1.6 KB, 11 views)
jyotir is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Strange axial velocity distribution in swirling flow aravind vashista Main CFD Forum 0 August 25, 2020 13:08
ActuatorDiskExplicitForce in OF2.1. Help be_inspired OpenFOAM Programming & Development 10 September 14, 2018 12:12
Need help in muti-stage axial compressor simulation using mixing plane pchoopanya FLUENT 0 July 9, 2018 03:51
plottting tangnetial and axial velocity vectors on a plane Shamoon Jamshed Tecplot 2 December 4, 2017 14:55
How to calcualte Axial flux of angular momentum and Axial flux of axial momentum ? nanavati OpenFOAM Running, Solving & CFD 0 November 21, 2014 10:47


All times are GMT -4. The time now is 22:27.