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

Covergence Problem with simpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 4, 2013, 12:11
Default Covergence Problem with simpleFoam
  #1
Senior Member
 
sfigato's Avatar
 
Marco Longhitano
Join Date: Jan 2013
Location: Aachen
Posts: 103
Rep Power: 4
sfigato is on a distinguished road
Send a message via Skype™ to sfigato
Hi Foamers,

we are simulating a valve with simpleFoam! We have carried out the first simulation by using a first order (upwind) numerical schems for divergence. So, we have switched it with a second order one (Gauss Linear). Here appears many problems with convergence!!

Anyone knows how set properly fvSchemes in order to have a seocnd order accurancy??

Here is our fvSchemes and fvSolutiuon:



ddtSchemes
{
default steadyState;
}

gradSchemes
{
default cellLimited Gauss linear 1;
grad(p) cellLimited Gauss linear 1;
grad(U) cellLimited Gauss linear 1;
grad(omega) cellLimited Gauss linear 1;
grad(k) cellLimited Gauss linear 1;
}

divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV cellMDLimited Gauss linear 1;
div(phi,k) Gauss linearUpwind cellLimited Gauss linear 1;
div(phi,omega) Gauss linearUpwind cellLimited Gauss linear 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear limited 0.5;
laplacian((1|A(U)),p) Gauss linear limited 0.5;
laplacian(DkEff,k) Gauss linear limited 0.5;
laplacian(DepsilonEff,epsilon) Gauss linear limited 0.5;
laplacian(DomegaEff,omega) Gauss linear limited 0.5;
laplacian(DREff,R) Gauss linear limited 0.5;
laplacian(DnuTildaEff,nuTilda) Gauss linear limited 0.5;
}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}



__________________________________________________ ______-

solvers
{
p
{
solver GAMG;
tolerance 1e-08;
relTol 0.05;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration true;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;

}
U
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-08;
relTol 0.1;
nSweeps 1;
}

k
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-08;
relTol 0.01;
nSweeps 1;
}

omega
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-08;
relTol 0.01;
nSweeps 1;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 2;
residualControl
{
p 1e-5;
U 1e-5;
k 1e-5;
omega 1e-5;
}
}


relaxationFactors
{
default 0;
p 0.3;
U 0.7;
k 0.7;
omega 0.7;
nut 0.7;
}
}


Thx

Regards
Marco
sfigato is offline   Reply With Quote

Old   February 6, 2013, 08:31
Default
  #2
Senior Member
 
sfigato's Avatar
 
Marco Longhitano
Join Date: Jan 2013
Location: Aachen
Posts: 103
Rep Power: 4
sfigato is on a distinguished road
Send a message via Skype™ to sfigato
Hi Foamers,

no one has an idea about which problem could be???
Pleaseee answer!!

Regards
Marco
sfigato 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
2D foil with simpleFoam: stability problem Nolwenn OpenFOAM Running, Solving & CFD 0 July 5, 2012 05:16
Trying to run a benchmark case with simpleFoam spsb OpenFOAM 3 February 24, 2012 10:07
SimpleFOAM + SST-Model + problem with convergence A.Devesa OpenFOAM Running, Solving & CFD 0 November 9, 2010 05:43
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13


All times are GMT -4. The time now is 20:19.