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 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

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

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

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::printStack(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::pureMixture<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::pureMixture<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

Same Problem
 

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::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 trial-and-error.