laplacianFoam: Added support for fvOptions

This commit is contained in:
Henry Weller
2016-07-08 11:34:00 +01:00
parent 10ad855d8d
commit c4390f7059

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -30,6 +30,7 @@ Description
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
#include "fvCFD.H" #include "fvCFD.H"
#include "fvOptions.H"
#include "simpleControl.H" #include "simpleControl.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -44,6 +45,7 @@ int main(int argc, char *argv[])
simpleControl simple(mesh); simpleControl simple(mesh);
#include "createFields.H" #include "createFields.H"
#include "createFvOptions.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -55,10 +57,16 @@ int main(int argc, char *argv[])
while (simple.correctNonOrthogonal()) while (simple.correctNonOrthogonal())
{ {
solve fvScalarMatrix TEqn
( (
fvm::ddt(T) - fvm::laplacian(DT, T) fvm::ddt(T) - fvm::laplacian(DT, T)
==
fvOptions(T)
); );
fvOptions.constrain(TEqn);
TEqn.solve();
fvOptions.correct(T);
} }
#include "write.H" #include "write.H"