driftFluxFoam: Separated control of solution of diffusion from bounded advection

This commit is contained in:
Henry
2014-02-24 11:18:30 +00:00
parent 01703a904d
commit e0987483a2
7 changed files with 59 additions and 25 deletions

View File

@ -57,7 +57,7 @@
- fvm::laplacian(mut/rho, alpha1) - fvm::laplacian(mut/rho, alpha1)
); );
alpha1Eqn.solve(); alpha1Eqn.solve(mesh.solver("alpha1Diffusion"));
phiAlpha += alpha1Eqn.flux(); phiAlpha += alpha1Eqn.flux();
} }

View File

@ -47,7 +47,7 @@ runTimeModifiable yes;
adjustTimeStep on; adjustTimeStep on;
maxCo 1; maxCo 5;
maxDeltaT 1; maxDeltaT 1;

View File

@ -19,8 +19,12 @@ solvers
{ {
"alpha1.*" "alpha1.*"
{ {
nAlphaCorr 1; nAlphaCorr 3;
nAlphaSubCycles 4; nAlphaSubCycles 1;
MULESCorr yes;
nLimiterIter 5;
alphaApplyPrevCorr yes;
solver smoothSolver; solver smoothSolver;
smoother symGaussSeidel; smoother symGaussSeidel;
@ -29,6 +33,15 @@ solvers
minIter 1; minIter 1;
} }
alpha1Diffusion
{
solver PCG;
preconditioner DIC;
tolerance 1e-6;
relTol 0;
minIter 1;
}
p_rgh p_rgh
{ {
solver GAMG; solver GAMG;
@ -49,7 +62,8 @@ solvers
"(U|k|epsilon)" "(U|k|epsilon)"
{ {
solver smoothSolver; solver PBiCG;
preconditioner DILU;
smoother symGaussSeidel; smoother symGaussSeidel;
tolerance 1e-7; tolerance 1e-7;
relTol 0.1; relTol 0.1;
@ -65,20 +79,16 @@ solvers
PIMPLE PIMPLE
{ {
nCorrectors 2; momentumPredictor no;
nCorrectors 3;
nNonOrthogonalCorrectors 0; nNonOrthogonalCorrectors 0;
} }
relaxationFactors relaxationFactors
{ {
fields
{
}
equations equations
{ {
"U.*" 1; ".*" 1;
"k.*" 1;
"epsilon.*" 1;
} }
} }

View File

@ -47,8 +47,7 @@ runTimeModifiable yes;
adjustTimeStep yes; adjustTimeStep yes;
maxCo 0.5; maxCo 1;
maxAlphaCo 0.5;
maxDeltaT 0.05; maxDeltaT 0.05;

View File

@ -19,13 +19,27 @@ solvers
{ {
"alpha1.*" "alpha1.*"
{ {
nAlphaCorr 1; nAlphaCorr 2;
nAlphaSubCycles 3; nAlphaSubCycles 1;
MULESCorr yes;
nLimiterIter 5;
alphaApplyPrevCorr yes;
solver smoothSolver; solver smoothSolver;
smoother symGaussSeidel; smoother symGaussSeidel;
tolerance 1e-6; tolerance 1e-6;
relTol 0; relTol 0;
minIter 1;
}
alpha1Diffusion
{
solver PCG;
preconditioner DIC;
tolerance 1e-6;
relTol 0;
minIter 1;
} }
"rho.*" "rho.*"
@ -71,7 +85,7 @@ solvers
PIMPLE PIMPLE
{ {
momentumPredictor yes; momentumPredictor no;
nCorrectors 3; nCorrectors 3;
nNonOrthogonalCorrectors 0; nNonOrthogonalCorrectors 0;

View File

@ -33,11 +33,11 @@ writeInterval 50;
purgeWrite 0; purgeWrite 0;
writeFormat ascii; writeFormat binary;
writePrecision 6; writePrecision 6;
writeCompression compressed; writeCompression uncompressed;
timeFormat general; timeFormat general;

View File

@ -19,8 +19,12 @@ solvers
{ {
"alpha1.*" "alpha1.*"
{ {
nAlphaCorr 1; nAlphaCorr 5;
nAlphaSubCycles 4; nAlphaSubCycles 1;
MULESCorr no; //yes;
nLimiterIter 5;
alphaApplyPrevCorr yes;
solver smoothSolver; solver smoothSolver;
smoother symGaussSeidel; smoother symGaussSeidel;
@ -29,6 +33,15 @@ solvers
minIter 1; minIter 1;
} }
alpha1Diffusion
{
solver PCG;
preconditioner DIC;
tolerance 1e-6;
relTol 0;
minIter 1;
}
p_rgh p_rgh
{ {
solver GAMG; solver GAMG;
@ -65,15 +78,13 @@ solvers
PIMPLE PIMPLE
{ {
nCorrectors 2; momentumPredictor no;
nCorrectors 3;
nNonOrthogonalCorrectors 0; nNonOrthogonalCorrectors 0;
} }
relaxationFactors relaxationFactors
{ {
fields
{
}
equations equations
{ {
".*" 1; ".*" 1;