|
[Sponsors] |
please help me which solver is better for my application |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 22, 2010, 06:29 |
please help me which solver is better for my application
|
#1 |
New Member
usn
Join Date: Jul 2010
Location: Deutschland
Posts: 9
Rep Power: 16 |
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 pressure-pressure 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 |
|
July 22, 2010, 06:54 |
|
#2 |
Member
MSarkar
Join Date: Dec 2009
Posts: 99
Rep Power: 16 |
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 |
|
July 23, 2010, 06:06 |
|
#3 |
New Member
usn
Join Date: Jul 2010
Location: Deutschland
Posts: 9
Rep Power: 16 |
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 |
|
July 26, 2010, 05:26 |
|
#4 |
New Member
usn
Join Date: Jul 2010
Location: Deutschland
Posts: 9
Rep Power: 16 |
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.831e-05 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((rho|A(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*(1|A(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 1e-05; relTol 0; } p { solver GAMG; tolerance 1e-08; 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 1e-05; relTol 0; } h { solver PBiCG; preconditioner DILU; tolerance 1e-06; relTol 0.1; } R { solver PBiCG; preconditioner DILU; tolerance 1e-05; relTol 0; } 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; } omega { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e-07; 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 |
|
July 26, 2010, 05:34 |
|
#5 |
New Member
usn
Join Date: Jul 2010
Location: Deutschland
Posts: 9
Rep Power: 16 |
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.2714e-06, No Iterations 3 DILUPBiCG: Solving for Uy, Initial residual = 0.0729241, Final residual = 4.88299e-06, No Iterations 3 DILUPBiCG: Solving for Uz, Initial residual = 0.0808816, Final residual = 3.54645e-06, 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.21321e-23, Final residual = 2.21321e-23, No Iterations 0 smoothSolver: Solving for k, Initial residual = 0.0534275, Final residual = 2.55456e-05, 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.05299e-06, No Iterations 4 DILUPBiCG: Solving for Uy, Initial residual = 0.0646994, Final residual = 4.51054e-06, No Iterations 3 DILUPBiCG: Solving for Uz, Initial residual = 0.0725841, Final residual = 1.03022e-06, 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:rintStack(Foam::Ostream&) in "/server/appl/Programme/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/server/appl/Programme/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #2 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::calculate() in "/server/appl/Programme/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so" #3 Foam::hPsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::correct() in "/server/appl/Programme/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so" #4 main in "/server/appl/Programme/OpenFOAM/OpenFOAM-1.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/OpenFOAM-1.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 cn-12 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 |
|
July 26, 2010, 10:29 |
Same Problem
|
#6 |
New Member
Alexander
Join Date: May 2010
Posts: 4
Rep Power: 16 |
Same Problem:
Time = 1.65e-05 Courant Number mean: 5.0346e-06 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.81638e-06, No Iterations 4 DILUPBiCG: Solving for Uy, Initial residual = 0.937236, Final residual = 8.94623e-06, No Iterations 4 DILUPBiCG: Solving for Uz, Initial residual = 0.814971, Final residual = 3.0321e-06, No Iterations 4 DILUPBiCG: Solving for e, Initial residual = 0.805664, Final residual = 2.73097e-06, 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/OpenFOAM-1.7.0/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 67. FOAM aborting #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam170/lib/linuxGccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/opt/openfoam170/lib/linuxGccDPOpt/libOpenFOAM.so" #2 Foam::ePsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::eConstThermo<Foam:: perfectGas> > > > >::calculate() in "/opt/openfoam170/lib/linuxGccDPOpt/libbasicThermophysicalModels.so" #3 Foam::ePsiThermo<Foam:ureMixture<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? |
|
July 27, 2010, 04:14 |
|
#7 |
New Member
usn
Join Date: Jul 2010
Location: Deutschland
Posts: 9
Rep Power: 16 |
Hi everyone,
if anyone knows the above problem....please help me |
|
January 8, 2013, 05:16 |
|
#8 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
i have same error message as alexander has had.will it be solved if i make mesh finer?or its irrelevant?
|
|
January 8, 2013, 13:29 |
|
#9 |
Senior Member
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 740
Rep Power: 28 |
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 trial-and-error. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Is Playstation 3 cluster suitable for CFD work | hsieh | OpenFOAM | 9 | August 16, 2015 15:53 |
Working directory via command line | Luiz | CFX | 4 | March 6, 2011 21:02 |
Solver dirtyness | Programmer | Main CFD Forum | 1 | June 20, 2008 09:07 |
why the solver reject it? Anyone with experience? | bearcat | CFX | 6 | April 28, 2008 15:08 |
User Solver in FoamX with a new field variable | vvqf | OpenFOAM Pre-Processing | 2 | October 28, 2005 21:24 |