zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

interFoam/ras/damBreak: updated for predictor-corrector MULES and changes to phi in pEqn

Browse Source
This commit is contained in:
Henry
2014-02-03 18:04:18 +00:00
parent faf867c84c
commit a67dfdb23a
8 changed files with 68 additions and 64 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
tutorials/multiphase/interFoam/ras/damBreak/system/controlDict
View File

@ -47,8 +47,8 @@ runTimeModifiable yes;
adjustTimeStep on;
maxCo 0.2;
maxAlphaCo 0.2;
maxCo 1;
maxAlphaCo 1;
maxDeltaT 1;

7
tutorials/multiphase/interFoam/ras/damBreak/system/fvSchemes
View File

@ -27,14 +27,11 @@ gradSchemes
divSchemes
{
div(rhoPhi,U) Gauss linear;
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss linear;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
div((muEff*dev(T(grad(U))))) Gauss linear;
}
@ -58,7 +55,7 @@ fluxRequired
default no;
p_rgh;
pcorr;
alpha;
alpha.water;
}

46
tutorials/multiphase/interFoam/ras/damBreak/system/fvSolution
View File

@ -17,18 +17,27 @@ FoamFile
solvers
{
alpha.water
"alpha.water.*"
{
nAlphaCorr 1;
nAlphaSubCycles 4;
cAlpha 2;
nAlphaCorr 2;
nAlphaSubCycles 1;
alphaOuterCorrectors yes;
cAlpha 1;
MULESCorr yes;
nLimiterIter 3;
solver PBiCG;
preconditioner DILU;
tolerance 1e-8;
relTol 0;
}
pcorr
{
solver PCG;
preconditioner DIC;
tolerance 1e-10;
tolerance 1e-5;
relTol 0;
}
@ -42,35 +51,34 @@ solvers
p_rghFinal
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
$p_rgh;
relTol 0;
}
"(U|k|epsilon)"
"(U|k|epsilon).*"
{
solver PBiCG;
preconditioner DILU;
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-06;
relTol 0;
}
"(U|k|epsilon)Final"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0;
minIter 1;
}
}
PIMPLE
{
momentumPredictor no;
nOuterCorrectors 1;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
}
relaxationFactors
{
equations
{
".*" 1;
}
}
// ************************************************************************* //

8
tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/0/p_rgh
View File

@ -58,18 +58,18 @@ boundaryField
patchType cyclic;
jump uniform 0;
value uniform 0;
D 700;
D 1000;
I 500;
length 1.05;
length 0.15;
}
porous_half1
{
type porousBafflePressure;
patchType cyclic;
value uniform 0;
D 700;
D 1000;
I 500;
length 1.05;
length 0.15;
}
}

2
tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/controlDict
View File

@ -29,7 +29,7 @@ deltaT 0.001;
writeControl adjustableRunTime;
writeInterval 0.1;
writeInterval 0.05;
purgeWrite 0;

4
tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/createBafflesDict
View File

@ -49,9 +49,9 @@ baffles
{
type porousBafflePressure;
patchType cyclic;
D 700;
D 1000;
I 500;
length 1.05;
length 0.15;
jump uniform 0;
value uniform 0;
}

7
tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/fvSchemes
View File

@ -27,14 +27,11 @@ gradSchemes
divSchemes
{
div(rhoPhi,U) Gauss linear;
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss linear;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
div((muEff*dev(T(grad(U))))) Gauss linear;
}
@ -58,7 +55,7 @@ fluxRequired
default no;
p_rgh;
pcorr;
alpha;
alpha.water;
}

48
tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/fvSolution
View File

@ -17,18 +17,27 @@ FoamFile
solvers
{
alpha.water
"alpha.water.*"
{
nAlphaCorr 1;
nAlphaSubCycles 4;
cAlpha 2;
nAlphaCorr 2;
nAlphaSubCycles 1;
alphaOuterCorrectors yes;
cAlpha 1;
MULESCorr yes;
nLimiterIter 3;
solver PBiCG;
preconditioner DILU;
tolerance 1e-8;
relTol 0;
}
pcorr
{
solver PCG;
preconditioner DIC;
tolerance 1e-10;
tolerance 1e-5;
relTol 0;
}
@ -42,41 +51,34 @@ solvers
p_rghFinal
{
solver PCG;
preconditioner DIC;
tolerance 1e-08;
$p_rgh;
relTol 0;
}
"(U|k|epsilon)"
"(U|k|epsilon).*"
{
solver PBiCG;
preconditioner DILU;
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-06;
relTol 0;
}
"(U|k|epsilon)Final"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-07;
relTol 0;
minIter 1;
}
}
PIMPLE
{
momentumPredictor yes;
momentumPredictor no;
nOuterCorrectors 1;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
}
relaxationFactors
{
"U.*" 1;
"k.*" 1;
"epsilon.*" 1;
equations
{
".*" 1;
}
}
// ************************************************************************* //