CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

please help me which solver is better for my application

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 22, 2010, 05:29
Default please help me which solver is better for my application
  #1
New Member
 
usn
Join Date: Jul 2010
Location: Deutschland
Posts: 9
Rep Power: 7
Ger_US is on a distinguished road
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
Ger_US is offline   Reply With Quote

Old   July 22, 2010, 05:54
Default
  #2
Member
 
MSarkar
Join Date: Dec 2009
Posts: 99
Rep Power: 7
msarkar is on a distinguished road
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
msarkar is offline   Reply With Quote

Old   July 23, 2010, 05:06
Default
  #3
New Member
 
usn
Join Date: Jul 2010
Location: Deutschland
Posts: 9
Rep Power: 7
Ger_US is on a distinguished road
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
Ger_US is offline   Reply With Quote

Old   July 26, 2010, 04:26
Default
  #4
New Member
 
usn
Join Date: Jul 2010
Location: Deutschland
Posts: 9
Rep Power: 7
Ger_US is on a distinguished road
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


Attached Images
File Type: jpg OpenFoam.jpg (39.9 KB, 23 views)
File Type: jpg cfx.jpg (36.1 KB, 24 views)
Ger_US is offline   Reply With Quote

Old   July 26, 2010, 04:34
Default
  #5
New Member
 
usn
Join Date: Jul 2010
Location: Deutschland
Posts: 9
Rep Power: 7
Ger_US is on a distinguished road
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
Ger_US is offline   Reply With Quote

Old   July 26, 2010, 09:29
Default Same Problem
  #6
New Member
 
Alexander
Join Date: May 2010
Posts: 4
Rep Power: 7
wise is on a distinguished road
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?
wise is offline   Reply With Quote

Old   July 27, 2010, 03:14
Default
  #7
New Member
 
usn
Join Date: Jul 2010
Location: Deutschland
Posts: 9
Rep Power: 7
Ger_US is on a distinguished road
Hi everyone,

if anyone knows the above problem....please help me
Ger_US is offline   Reply With Quote

Old   January 8, 2013, 05:16
Default
  #8
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,205
Rep Power: 17
immortality is on a distinguished road
i have same error message as alexander has had.will it be solved if i make mesh finer?or its irrelevant?
immortality is offline   Reply With Quote

Old   January 8, 2013, 13:29
Default
  #9
Senior Member
 
mturcios777's Avatar
 
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 727
Rep Power: 18
mturcios777 will become famous soon enough
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.
mturcios777 is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Is Playstation 3 cluster suitable for CFD work hsieh OpenFOAM 9 August 16, 2015 14:53
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
Solver dirtyness Programmer Main CFD Forum 1 June 20, 2008 08:07
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08
User Solver in FoamX with a new field variable vvqf OpenFOAM Pre-Processing 2 October 28, 2005 20:24


All times are GMT -4. The time now is 16:47.