Coupled solver update

This commit is contained in:
Hrvoje Jasak 2015-10-28 16:34:23 +00:00
parent 49d85398e8
commit c9af015399
6 changed files with 106 additions and 43 deletions

View file

@ -22,13 +22,13 @@ startTime 0;
stopAt endTime; stopAt endTime;
endTime 300; endTime 100;
deltaT 1; deltaT 1;
writeControl timeStep; writeControl timeStep;
writeInterval 50; writeInterval 20;
purgeWrite 0; purgeWrite 0;
@ -36,7 +36,7 @@ writeFormat ascii;
writePrecision 6; writePrecision 6;
writeCompression uncompressed; writeCompression compressed;
timeFormat general; timeFormat general;
@ -44,4 +44,27 @@ timePrecision 6;
runTimeModifiable yes; runTimeModifiable yes;
functions
(
minMaxU
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name U;
}
minMaxP
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name p;
}
);
// ************************************************************************* // // ************************************************************************* //

View file

@ -28,12 +28,8 @@ gradSchemes
divSchemes divSchemes
{ {
default none; default none;
div(phi,U) Gauss upwind; // div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind; div(phi,U) Gauss linearUpwind Gauss linear;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
div((nuEff*dev(grad(U).T()))) Gauss linear; div((nuEff*dev(grad(U).T()))) Gauss linear;
div((grad(p)|A(U))) Gauss linear; div((grad(p)|A(U))) Gauss linear;
div(U) Gauss linear; div(U) Gauss linear;

View file

@ -42,13 +42,16 @@ solvers
minCoarseEqns 4; minCoarseEqns 4;
nMaxLevels 100; nMaxLevels 100;
scale on; scale on;
coarseSmoother ILU; smoother
fineSmoother ILU; {
smoother ILUCp;
fillInLevel 1;
}
minIter 0; minIter 0;
maxIter 100; maxIter 100;
tolerance 1e-7; tolerance 1e-7;
relTol 0.005; relTol 0.0;
} }
// Segregated // Segregated
@ -86,7 +89,7 @@ fieldBounds
relaxationFactors relaxationFactors
{ {
U 0.95; // U 0.95;
} }
// ************************************************************************* // // ************************************************************************* //

View file

@ -3,7 +3,7 @@
| \\ / F ield | foam-extend: Open Source CFD | | \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.2 | | \\ / O peration | Version: 3.2 |
| \\ / A nd | Web: http://www.foam-extend.org | | \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | | | \\/ M anipulation | For copyright notice see file Copyright |
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
FoamFile FoamFile
{ {

View file

@ -44,4 +44,54 @@ timePrecision 6;
runTimeModifiable yes; runTimeModifiable yes;
functions
(
minMaxU
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name U;
}
minMaxP
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name p;
}
minMaxK
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name k;
}
minMaxEpsilon
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name epsilon;
}
minMaxOmega
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name omega;
}
);
// ************************************************************************* // // ************************************************************************* //

View file

@ -18,38 +18,29 @@ solvers
{ {
Up Up
{ {
solver BiCGStab; solver AMG;
preconditioner Cholesky; cycle V-cycle;
coarseningType AAMG;
norm componentNorm;
normComponent 0;
tolerance 1e-09; nPreSweeps 2;
relTol 0.0; nPostSweeps 2;
groupSize 2;
minIter 1; minCoarseEqns 4;
maxIter 500; nMaxLevels 100;
scale on;
smoother
{
smoother ILUCp;
fillInLevel 1;
} }
// Up minIter 0;
// { maxIter 100;
// solver AMG; tolerance 1e-7;
// cycle V-cycle; relTol 0.001;
// coarseningType AAMG; }
// norm componentNorm;
// normComponent 0;
// nPreSweeps 2;
// nPostSweeps 2;
// groupSize 2;
// minCoarseEqns 4;
// nMaxLevels 100;
// scale on;
// fineSmoother ILU;
// coarseSmoother k ILU;
// minIter 0;
// maxIter 100;
// tolerance 1e-7;
// relTol 0.001;
// }
// Segregated // Segregated
p p