mechy |
August 1, 2013 05:56 |
PisoFoam and PimpleFoam for my problem can not converge
4 Attachment(s)
I want to solve flow around a cylinder and thin plate
the thin plate is placed horizontally after cylinder
I have test many different settings in fvSchemes and fvSolutions
but all of them give uncorrected results
I will be so grateful if anybody can help me
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default Euler;
}
gradSchemes
{
default leastSquares ;
// grad(p) leastSquares ;
// grad(U) leastSquares ;
}
divSchemes
{
default none;
div(phi,U) Gauss limitedLinearV 1;//Gauss linearUpwind grad(U);//
div(phi,k) Gauss limitedLinear 1;
div(phi,omega) Gauss limitedLinear 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;//Gauss linear limited 0.5;//
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
// pcorr ;
p;
}
// ************************************************************************* //
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.3 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 1e-5;
relTol 0;
smoother GaussSeidel;
nPreSweeps 1;
nPostSweeps 2;
cacheAgglomeration true;
directSolveCoarsest true;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 400;
mergeLevels 1;
minIter 1;
}
pFinal
{
solver GAMG;
tolerance 1e-9;
relTol 0;
smoother GaussSeidel;
nPreSweeps 1;
nPostSweeps 2;
cacheAgglomeration true;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 400;
mergeLevels 1;
minIter 1;
}
"(U|k|omega)"
{
// solver PBiCG;
// preconditioner DILU;
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-5;
relTol 0;
nSweeps 1;
minIter 1;
}
"(UFinal|kFinal|omegaFinal)"
{
// solver PBiCG;
// preconditioner DILU;
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-6;
relTol 0.0;
nSweeps 1;
minIter 1;
}
cellMotionU
{
solver GAMG;
tolerance 1e-6;
relTol 1e-3;
smoother GaussSeidel;
cacheAgglomeration true;
directSolveCoarsest true;
nCellsInCoarsestLevel 20;
agglomerator faceAreaPair;
mergeLevels 1;
}
}
PIMPLE
{
nOuterCorrectors 3;
nCorrectors 2;
nNonOrthogonalCorrectors 3;
pRefCell 0;
pRefValue 0;
correctPhi no;
// momentumPredictor yes;
}
PISO
{
nOuterCorrectors 2;
nCorrectors 2;
nNonOrthogonalCorrectors 3;
pRefCell 0;
pRefValue 0;
// ddtPhiCorr no;
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.7;
k 0.7;
omega 0.7;
}
}
potentialFlow
{
nNonOrthogonalCorrectors 0;
}
cache
{
grad(U);
}
// ************************************************************************* //
|