please help me which solver is better for my application
Hai everyone,
I am new to OpenFoam...my Geometry (Inlet Duct, constant Inlet Valve position, closed outlet valve and cylinder) is meshed in Icem...mesh size is approximately 5 millionen...now I would like to simulate in OpenFoam but I don't know which solver is suitable for my application...I am searching for compressible steady state Komega sst solver. I tried with rhoPisoFoam but it didn't work....please help me is there any other compressible steady state Komega sst solver available in OpenFoam...and one more thing I would like to give pressurepressure boundary conditions at Inlet and Outlet Inlet at Inlet Duct: total pressure: 0 Pa temperature: 298 k Outlet at Cylinder Bottom: static pressure: 5000 Pa temperature: 299.5 k Thank you in advance 
Hi
You mentioned you want to use steady state compressible flow solver but you used rhoPisoFoam which is a transient solver. However, you can try rhoSimpleFoam that is steady state SIMPLE solver for laminar or turbulent RANS flow of compressible fluid. Hope it will help you. M. Sarkar 
Hi Sarkar,
thank you very much for your reply...ja I know about Piso just I tried to know is there any problem with my mesh...now I am trying with rhoSimpleFoam for sst...soon I give reply whether I succeed or not thank you again 
2 Attachment(s)
Hi Sarkar
I tried with rhoSimple Foam...but not reached to my cfx results....here I am providing my data...please tell me if I am doing anything wrong U: INLET { type zeroGradient; } OUTLET { type zeroGradient; } WALL { type zeroGradient; } p: INLET { //type fixedValue; //value uniform 98038.6; //type zeroGradient; type totalPressure; p0 uniform 1e5; value uniform 1e5; gamma 1.4; psi psi; //rho rho; } OUTLET { type fixedValue; value uniform 95000; } WALL { type zeroGradient; } omega: INLET { type fixedValue; value uniform 2080.55; } OUTLET { type zeroGradient; } WALL { type zeroGradient; } k: INLET { type fixedValue; value uniform 16.977; } OUTLET { type zeroGradient; } WALL { type zeroGradient; } epsilon: INLET { type fixedValue; value uniform 3178.95; } OUTLET { type zeroGradient; } WALL { type zeroGradient; } RAS Properties: Turbulence Properties: RASModel kOmegaSST; simulationType RASModel; turbulence on; printCoeffs on; thermophysicalProperties: thermoType hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>; mixture air 1 28.96 1004.4 0 1.831e05 0.7142; fvSchemes: ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; grad(U) Gauss linear; grad(p) Gauss linear; } divSchemes { div(phi,U) Gauss upwind; div((muEff*dev2(grad(U).T()))) Gauss linear; div(phi,h) Gauss upwind; div(phi,epsilon) Gauss upwind; div(phi,k) Gauss upwind; div(phid,p) Gauss linear; div(U,p) Gauss linear; div(phiU,p) Gauss linear; div(phi,R) Gauss upwind; div(phi,omega) Gauss upwind; div((rho*R)) Gauss linear; div(R) Gauss linear; div(U) Gauss linear; } laplacianSchemes { laplacian(muEff,U) Gauss linear corrected; laplacian(mut,U) Gauss linear corrected; laplacian(alphaEff,h) Gauss linear corrected; laplacian((rhoA(U)),p) Gauss linear corrected; laplacian((rho*rAU),p) Gauss linear corrected; laplacian(DepsilonEff,epsilon) Gauss linear corrected; laplacian(DREff,R) Gauss linear corrected; laplacian(DomegaEff,omega) Gauss linear corrected; laplacian(DkEff,k) Gauss linear corrected; laplacian(1,p) Gauss linear corrected; laplacian((rho*(1A(U))),p) Gauss linear corrected; } interpolationSchemes { default linear; div(U,p) upwind phi; } snGradSchemes { default corrected; } fluxRequired { default no; p ; } fvSolution: rho { solver PCG; preconditioner DIC; tolerance 1e05; relTol 0; } p { solver GAMG; tolerance 1e08; relTol 0.05; smoother GaussSeidel; cacheAgglomeration false; //off > true nCellsInCoarsestLevel 10; // 20 > 10 Fehlermeldung damit behoben agglomerator faceAreaPair; mergeLevels 1; maxIter 500; } U { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0; } h { solver PBiCG; preconditioner DILU; tolerance 1e06; relTol 0.1; } R { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0; } k { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e07; relTol 0.1; } epsilon { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e07; relTol 0.1; } omega { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e07; relTol 0.1; } } SIMPLE { nNonOrthogonalCorrectors 0; pMin pMin [ 1 1 2 0 0 0 0 ] 100; } relaxationFactors { p 0.1; rho 0.05; U 0.3; k 0.05; epsilon 0.05; omega 0.05; h 0.3; R 0.05; } and I attached the results to this message as jpg files please see my data and help me where I am doing wrong....I am new to OpenFoam and also new to CFD thank you in advance 
Hi
and I forgot to tell you that the simulation diverged after some iterations...here is the result... Time = 1299 DILUPBiCG: Solving for Ux, Initial residual = 0.103001, Final residual = 6.2714e06, No Iterations 3 DILUPBiCG: Solving for Uy, Initial residual = 0.0729241, Final residual = 4.88299e06, No Iterations 3 DILUPBiCG: Solving for Uz, Initial residual = 0.0808816, Final residual = 3.54645e06, No Iterations 3 DILUPBiCG: Solving for h, Initial residual = 0.478437, Final residual = 0.0172015, No Iterations 1 GAMG: Solving for p, Initial residual = 0.69215, Final residual = 0.0323442, No Iterations 15 time step continuity errors : sum local = 1241.23, global = 604.855, cumulative = 836.722 rho max/min : 179.449 2205.68 bounding p, min: 3.68828e+09 max: 7.43339e+07 average: 1.19258e+07 smoothSolver: Solving for omega, Initial residual = 2.21321e23, Final residual = 2.21321e23, No Iterations 0 smoothSolver: Solving for k, Initial residual = 0.0534275, Final residual = 2.55456e05, No Iterations 2 bounding k, min: 88298.8 max: 2479.51 average: 5.23806 ExecutionTime = 51147.7 s ClockTime = 51181 s Time = 1300 DILUPBiCG: Solving for Ux, Initial residual = 0.0840612, Final residual = 1.05299e06, No Iterations 4 DILUPBiCG: Solving for Uy, Initial residual = 0.0646994, Final residual = 4.51054e06, No Iterations 3 DILUPBiCG: Solving for Uz, Initial residual = 0.0725841, Final residual = 1.03022e06, No Iterations 4 DILUPBiCG: Solving for h, Initial residual = 0.709417, Final residual = 0.0270025, No Iterations 1 [2] [2] [2] Maximum number of iterations exceeded#0 Foam::error::printStack(Foam::Ostream&) in "/server/appl/Programme/OpenFOAM/OpenFOAM1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/server/appl/Programme/OpenFOAM/OpenFOAM1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #2 Foam::hPsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::calculate() in "/server/appl/Programme/OpenFOAM/OpenFOAM1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so" #3 Foam::hPsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::correct() in "/server/appl/Programme/OpenFOAM/OpenFOAM1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so" #4 main in "/server/appl/Programme/OpenFOAM/OpenFOAM1.6/applications/bin/linux64GccDPOpt/rhoSimpleFoam" #5 __libc_start_main in "/lib64/tls/libc.so.6" #6 _start at ../sysdeps/x86_64/elf/start.S:116 [2] [2] [2] From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const [2] in file /home/dm2/henry/OpenFOAM/OpenFOAM1.6/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 68. [2] FOAM parallel run aborting [2]  MPI_ABORT was invoked on rank 2 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them.   mpirun has exited due to process rank 2 with PID 10612 on node cn12 exiting without calling "finalize". This may have caused other processes in the application to be terminated by signals sent by mpirun (as reported here).  please help me Sarkar or anyone please thank you in advance 
Same Problem
Same Problem:
Time = 1.65e05 Courant Number mean: 5.0346e06 max: 11.2864 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Ux, Initial residual = 0.930813, Final residual = 5.81638e06, No Iterations 4 DILUPBiCG: Solving for Uy, Initial residual = 0.937236, Final residual = 8.94623e06, No Iterations 4 DILUPBiCG: Solving for Uz, Initial residual = 0.814971, Final residual = 3.0321e06, No Iterations 4 DILUPBiCG: Solving for e, Initial residual = 0.805664, Final residual = 2.73097e06, No Iterations 4 > FOAM FATAL ERROR: Maximum number of iterations exceeded From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/opencfd/OpenFOAM/OpenFOAM1.7.0/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 67. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam170/lib/linuxGccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/opt/openfoam170/lib/linuxGccDPOpt/libOpenFOAM.so" #2 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::eConstThermo<Foam:: perfectGas> > > > >::calculate() in "/opt/openfoam170/lib/linuxGccDPOpt/libbasicThermophysicalModels.so" #3 Foam::ePsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::eConstThermo<Foam:: perfectGas> > > > >::correct() in "/opt/openfoam170/lib/linuxGccDPOpt/libbasicThermophysicalModels.so" #4 in "/opt/openfoam170/applications/bin/linuxGccDPOpt/sonicFoam" #5 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6" #6 in "/opt/openfoam170/applications/bin/linuxGccDPOpt/sonicFoam" Why? 
Hi everyone,
if anyone knows the above problem....please help me 
i have same error message as alexander has had.will it be solved if i make mesh finer?or its irrelevant?

The issue is that the extraction of Temperature from enthalpy is failing. Basically we can't find a temperature that makes sense for the enthalpy. I would recommend looking at your temperature field at the few timesteps before the crash. This will give you better information as to what is going on.
This type of error is difficult to diagnose over the Forum as they are case specific, and is subject to trialanderror. 
All times are GMT 4. The time now is 19:46. 