4031066a5a
The function now takes care of making the flux relative if the mesh is moving. This is necessary in order to consistently reconstruct the face velocity from the flux and then call correctBoundaryConditions on U.
138 lines
4 KiB
C
138 lines
4 KiB
C
/*---------------------------------------------------------------------------*\
|
|
========= |
|
|
\\ / F ield | foam-extend: Open Source CFD
|
|
\\ / O peration | Version: 4.0
|
|
\\ / A nd | Web: http://www.foam-extend.org
|
|
\\/ M anipulation | For copyright notice see file Copyright
|
|
-------------------------------------------------------------------------------
|
|
License
|
|
This file is part of foam-extend.
|
|
|
|
foam-extend is free software: you can redistribute it and/or modify it
|
|
under the terms of the GNU General Public License as published by the
|
|
Free Software Foundation, either version 3 of the License, or (at your
|
|
option) any later version.
|
|
|
|
foam-extend is distributed in the hope that it will be useful, but
|
|
WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with foam-extend. If not, see <http://www.gnu.org/licenses/>.
|
|
|
|
Application
|
|
icoFoam
|
|
|
|
Description
|
|
Transient solver for incompressible, laminar flow of Newtonian fluids.
|
|
Consistent formulation without time-step and relaxation dependence by
|
|
Jasak and Tukovic.
|
|
|
|
Author
|
|
Hrvoje Jasak, Wikki Ltd. All rights reserved
|
|
|
|
\*---------------------------------------------------------------------------*/
|
|
|
|
#include "fvCFD.H"
|
|
#include "pisoControl.H"
|
|
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
|
|
int main(int argc, char *argv[])
|
|
{
|
|
# include "setRootCase.H"
|
|
# include "createTime.H"
|
|
# include "createMesh.H"
|
|
|
|
pisoControl piso(mesh);
|
|
|
|
# include "createFields.H"
|
|
# include "initContinuityErrs.H"
|
|
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
|
|
Info<< "\nStarting time loop\n" << endl;
|
|
|
|
while (runTime.loop())
|
|
{
|
|
Info<< "Time = " << runTime.timeName() << nl << endl;
|
|
|
|
# include "CourantNo.H"
|
|
|
|
// Time-derivative matrix
|
|
fvVectorMatrix ddtUEqn(fvm::ddt(U));
|
|
|
|
// Convection-diffusion matrix
|
|
fvVectorMatrix HUEqn
|
|
(
|
|
fvm::div(phi, U)
|
|
- fvm::laplacian(nu, U)
|
|
);
|
|
|
|
if (piso.momentumPredictor())
|
|
{
|
|
solve(ddtUEqn + HUEqn == -fvc::grad(p));
|
|
}
|
|
|
|
// Prepare clean 1/a_p without time derivative contribution
|
|
volScalarField rAU = 1.0/HUEqn.A();
|
|
|
|
// --- PISO loop
|
|
|
|
while (piso.correct())
|
|
{
|
|
// Calculate U from convection-diffusion matrix
|
|
U = rAU*HUEqn.H();
|
|
|
|
// Consistently calculate flux
|
|
piso.calcTransientConsistentFlux(phi, U, rAU, ddtUEqn);
|
|
|
|
adjustPhi(phi, U, p);
|
|
|
|
while (piso.correctNonOrthogonal())
|
|
{
|
|
fvScalarMatrix pEqn
|
|
(
|
|
fvm::laplacian
|
|
(
|
|
fvc::interpolate(rAU)/piso.aCoeff(),
|
|
p,
|
|
"laplacian(rAU," + p.name() + ')'
|
|
)
|
|
==
|
|
fvc::div(phi)
|
|
);
|
|
|
|
pEqn.setReference(pRefCell, pRefValue);
|
|
pEqn.solve
|
|
(
|
|
mesh.solutionDict().solver(p.select(piso.finalInnerIter()))
|
|
);
|
|
|
|
if (piso.finalNonOrthogonalIter())
|
|
{
|
|
phi -= pEqn.flux();
|
|
}
|
|
}
|
|
|
|
# include "continuityErrs.H"
|
|
|
|
// Consistently reconstruct velocity after pressure equation
|
|
piso.reconstructTransientVelocity(U, phi, ddtUEqn, rAU, p);
|
|
}
|
|
|
|
runTime.write();
|
|
|
|
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
|
|
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
|
|
<< nl << endl;
|
|
}
|
|
|
|
Info<< "End\n" << endl;
|
|
|
|
return 0;
|
|
}
|
|
|
|
|
|
// ************************************************************************* //
|