Test-volField: Corrected equation and test case to run

This commit is contained in:
Henry
2013-09-19 15:09:45 +01:00
parent 663603bf42
commit 2824e86e4f
3 changed files with 11 additions and 22 deletions

View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -22,7 +22,7 @@ License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
Application
volFieldTest
Test-volField
\*---------------------------------------------------------------------------*/
@ -52,8 +52,6 @@ int main(int argc, char *argv[])
mesh
);
//Info<< min(p, p);
Info<< "Reading field U\n" << endl;
volVectorField U
(
@ -85,9 +83,12 @@ int main(int argc, char *argv[])
zeroGradientFvPatchSymmTensorField::typeName
);
//Info<< fvc::div(st) << endl;
solve(fvm::ddt(st) + fvm::div(phi, st) - fvm::laplacian(st));
solve
(
fvm::ddt(st)
+ fvm::div(phi, st)
- fvm::laplacian(dimensionedScalar("D", sqr(dimLength)/dimTime, 1), st)
);
return 0;
}

View File

@ -23,26 +23,22 @@ ddtSchemes
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss linear;
div(phi,st) Gauss linear;
}
laplacianSchemes
{
default none;
laplacian(nu,U) Gauss linear orthogonal;
laplacian((1|A(U)),p) Gauss linear orthogonal;
default Gauss linear orthogonal;
}
interpolationSchemes
{
default linear;
interpolate(HbyA) linear;
}
snGradSchemes

View File

@ -17,15 +17,7 @@ FoamFile
solvers
{
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}
U
st
{
solver PBiCG;
preconditioner DILU;