|
[Sponsors] |
November 17, 2021, 11:49 |
Help sought on axial compressor simulation
|
#1 |
New Member
Jyoti Ranjan
Join Date: Oct 2020
Posts: 14
Rep Power: 6 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |