Hi, Makis
The Calculix FSI2 example here is not available for download:Turek-Hron FSI 2 Benchmark
I have modified my discrete scheme, but it is still far from the standard result. Do you have any better suggestions? This is my fvsolution file and fvscheme file.
ddtSchemes
{
default backward;
}
gradSchemes
{
default cellLimited Gauss linear 1;
}
divSchemes
{
default none;
div(phi,U) Gauss linearUpwind grad(U);
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear limited 1;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default limited 1;
}
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 1e-6;
relTol 1e-5;
smoother DICGaussSeidel;
}
pFinal
{
$p;
tolerance 1e-07;
relTol 0;
}
pcorr
{
solver GAMG;
tolerance 1e-5;
relTol 1e-4;
smoother GaussSeidel;
}
pcorrFinal
{
$pcorr;
relTol 0;
}
phi
{
$p;
}
"(U|cellDisplacement)"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-8;
relTol 1e-6;
}
"(U|cellDisplacement)Final"
{
$U;
relTol 0;
}
}
PIMPLE
{
nOuterCorrectors 50;
nCorrectors 2;
nNonOrthogonalCorrectors 1;
//tolerance 1.0e-12;
correctPhi yes;
//relTol 1e-4;
//pisoTol 1e-6;
consistent true;
residualControl
{
U
{
tolerance 1e-5;
relTol 0;
}
p
{
tolerance 1e-5;
relTol 0;
}
}
}
PISO
{
nNonOrthogonalCorrectors 1;
}
potentialFlow
{
nNonOrthogonalCorrectors 1;
}