exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Bugfix: update of conjugate heat solvers and tutorial settings

Browse source
This commit is contained in:
Hrvoje Jasak 2016-04-18 10:16:29 +01:00
parent 07b473621f
commit a946215a8d
11 changed files with 68 additions and 79 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
applications/solvers/coupled/conjugateHeatFoam/Make/options
View file

@ -19,4 +19,5 @@ EXE_LIBS = \
-lradiation \
-lbasicThermophysicalModels \
-lspecie \
-lconjugateHeatTransfer
-lconjugateHeatTransfer \
-llduSolvers

3
applications/solvers/coupled/conjugateHeatSimpleFoam/Make/options
View file

@ -19,4 +19,5 @@ EXE_LIBS = \
-lradiation \
-lbasicThermophysicalModels \
-lspecie \
-lconjugateHeatTransfer
-lconjugateHeatTransfer \
-llduSolvers

2
src/coupledMatrix/coupledFvMatrices/coupledFvMatrix/coupledFvMatrix.C
View file

@ -58,7 +58,7 @@ Foam::word Foam::coupledFvMatrix<Type>::coupledPsiName() const
if (rowI < matrices.size() - 1)
{
cpn += "+";
cpn += "_";
}
}

3
src/coupledMatrix/coupledFvMatrices/coupledFvScalarMatrix/coupledFvScalarMatrix.C
View file

@ -70,7 +70,8 @@ coupledSolverPerformance coupledFvMatrix<scalar>::solve
saveDiag.set(rowI, new scalarField(curMatrix.diag()));
// HR 17/Feb/2013
// Need to be able to compare references to support hacks such as in jumpCyclic
// Need to be able to compare references to support hacks
// such as in jumpCyclic
// psi.set(rowI, new scalarField(curMatrix.psi()));
psi.set(rowI, &curMatrix.psi());
source.set(rowI, new scalarField(curMatrix.source()));

3
tutorials/coupled/conjugateHeatFoam/conjugateCavity/0/solid/materials
View file

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

1
tutorials/coupled/conjugateHeatFoam/conjugateCavity/system/controlDict
View file

@ -59,7 +59,6 @@ functions
K kappaEff;
}
heatFlux
{
// Type of functionObject

13
tutorials/coupled/conjugateHeatFoam/conjugateCavity/system/fvSchemes
View file

@ -42,18 +42,17 @@ divSchemes
laplacianSchemes
{
default none;
laplacian((1|A(U)),p_rgh) Gauss linear uncorrected;
laplacian(kappaEff,T) Gauss linear uncorrected;
laplacian(gammaRad,G) Gauss linear uncorrected;
laplacian(DkEff,k) Gauss linear uncorrected;
laplacian(DepsilonEff,epsilon) Gauss linear uncorrected;
laplacian(nuEff,U) Gauss linear uncorrected;
laplacian((1|A(U)),p_rgh) Gauss linear corrected;
laplacian(kappaEff,T) Gauss linear corrected;
laplacian(gammaRad,G) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(nuEff,U) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(kappaEff) linear;
}

54
tutorials/coupled/conjugateHeatFoam/conjugateCavity/system/fvSolution
View file

@ -23,10 +23,9 @@ solvers
minIter 0;
maxIter 1000;
tolerance 1e-25;
tolerance 1e-7;
relTol 0.01;
};
}
p_rghFinal
{
solver PCG;
@ -34,32 +33,29 @@ solvers
minIter 0;
maxIter 1000;
tolerance 1e-25;
tolerance 1e-7;
relTol 0;
};
}
U
{
solver PBiCG;
solver BiCGStab;
preconditioner DILU;
minIter 1;
maxIter 1000;
tolerance 1e-25;
relTol 0.01;
};
T+T
tolerance 1e-7;
relTol 0;
}
T_T
{
solver BiCG;
solver BiCGStab;
preconditioner Cholesky;
minIter 0;
maxIter 1000;
tolerance 1e-25;
relTol 0.1;
};
tolerance 1e-7;
relTol 0;
}
G
{
solver PCG;
@ -67,31 +63,29 @@ solvers
minIter 1;
maxIter 1000;
tolerance 1e-25;
relTol 0.1;
};
tolerance 1e-7;
relTol 0;
}
epsilon
{
solver PBiCG;
solver BiCGStab;
preconditioner DILU;
minIter 1;
maxIter 100;
tolerance 1e-25;
relTol 0.01;
};
tolerance 1e-7;
relTol 0;
}
k
{
solver PBiCG;
solver BiCGStab;
preconditioner DILU;
minIter 1;
maxIter 100;
tolerance 1e-25;
relTol 0.01;
};
tolerance 1e-7;
relTol 0;
}
}
PISO

1
tutorials/coupled/conjugateHeatSimpleFoam/conjugateCavity/system/controlDict
View file

@ -55,7 +55,6 @@ functions
K kappaEff;
}
heatFlux
{
// Type of functionObject

13
tutorials/coupled/conjugateHeatSimpleFoam/conjugateCavity/system/fvSchemes
View file

@ -42,18 +42,17 @@ divSchemes
laplacianSchemes
{
default none;
laplacian((1|A(U)),p_rgh) Gauss linear uncorrected;
laplacian(kappaEff,T) Gauss linear uncorrected;
laplacian(gammaRad,G) Gauss linear uncorrected;
laplacian(DkEff,k) Gauss linear uncorrected;
laplacian(DepsilonEff,epsilon) Gauss linear uncorrected;
laplacian(nuEff,U) Gauss linear uncorrected;
laplacian((1|A(U)),p_rgh) Gauss linear corrected;
laplacian(kappaEff,T) Gauss linear corrected;
laplacian(gammaRad,G) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(nuEff,U) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(kappaEff) linear;
}

51
tutorials/coupled/conjugateHeatSimpleFoam/conjugateCavity/system/fvSolution
View file

@ -23,32 +23,29 @@ solvers
minIter 0;
maxIter 1000;
tolerance 1e-25;
tolerance 1e-7;
relTol 0.01;
};
}
U
{
solver PBiCG;
solver BiCGStab;
preconditioner DILU;
minIter 1;
maxIter 1000;
tolerance 1e-25;
tolerance 1e-7;
relTol 0.01;
};
T+T
}
T_T
{
solver BiCG;
solver BiCGStab;
preconditioner Cholesky;
minIter 0;
maxIter 1000;
tolerance 1e-25;
relTol 0.1;
};
tolerance 1e-7;
relTol 0.001;
}
G
{
solver PCG;
@ -56,31 +53,29 @@ solvers
minIter 1;
maxIter 1000;
tolerance 1e-25;
tolerance 1e-7;
relTol 0.1;
};
}
epsilon
{
solver PBiCG;
solver BiCGStab;
preconditioner DILU;
minIter 1;
maxIter 100;
tolerance 1e-25;
tolerance 1e-7;
relTol 0.01;
};
}
k
{
solver PBiCG;
solver BiCGStab;
preconditioner DILU;
minIter 1;
maxIter 100;
tolerance 1e-25;
tolerance 1e-7;
relTol 0.01;
};
}
}
SIMPLE
@ -92,11 +87,11 @@ SIMPLE
relaxationFactors
{
p_rgh 0.7;
U 0.3;
T 0.5;
k 0.6;
epsilon 0.6;
U 0.7;
p_rgh 0.3;
T 0.7;
k 0.7;
epsilon 0.7;
}
// ************************************************************************* //