0f5f4ed62c
This change lets the solver be looked up from the region database, so
that aspects of the solution algorithm can be accessed by
functionObjects and fvModels and similar. For example, a fluid solver
could be looked up as follows:
const solvers::fluid& s =
mesh().lookupObject<solvers::fluid>(solver::typeName);
161 lines
3.8 KiB
C++
161 lines
3.8 KiB
C++
/*---------------------------------------------------------------------------*\
|
|
========= |
|
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
|
\\ / O peration | Website: https://openfoam.org
|
|
\\ / A nd | Copyright (C) 2023 OpenFOAM Foundation
|
|
\\/ M anipulation |
|
|
-------------------------------------------------------------------------------
|
|
License
|
|
This file is part of OpenFOAM.
|
|
|
|
OpenFOAM 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.
|
|
|
|
OpenFOAM 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 OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
|
|
|
\*---------------------------------------------------------------------------*/
|
|
|
|
#include "compressibleMultiphaseVoF.H"
|
|
#include "CorrectPhi.H"
|
|
#include "geometricZeroField.H"
|
|
#include "fvcDdt.H"
|
|
#include "addToRunTimeSelectionTable.H"
|
|
|
|
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
|
|
|
|
namespace Foam
|
|
{
|
|
namespace solvers
|
|
{
|
|
defineTypeNameAndDebug(compressibleMultiphaseVoF, 0);
|
|
addToRunTimeSelectionTable(solver, compressibleMultiphaseVoF, fvMesh);
|
|
}
|
|
}
|
|
|
|
|
|
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
|
|
|
Foam::solvers::compressibleMultiphaseVoF::compressibleMultiphaseVoF
|
|
(
|
|
fvMesh& mesh
|
|
)
|
|
:
|
|
multiphaseVoFSolver
|
|
(
|
|
mesh,
|
|
autoPtr<multiphaseVoFMixture>
|
|
(
|
|
new compressibleMultiphaseVoFMixture(mesh)
|
|
)
|
|
),
|
|
|
|
mixture
|
|
(
|
|
refCast<compressibleMultiphaseVoFMixture>
|
|
(
|
|
multiphaseVoFSolver::mixture
|
|
)
|
|
),
|
|
|
|
phases(mixture.phases()),
|
|
|
|
p(mixture.p()),
|
|
|
|
pressureReference_
|
|
(
|
|
p,
|
|
p_rgh,
|
|
pimple.dict(),
|
|
false
|
|
),
|
|
|
|
pMin
|
|
(
|
|
"pMin",
|
|
dimPressure,
|
|
mixture
|
|
),
|
|
|
|
K("K", 0.5*magSqr(U)),
|
|
|
|
momentumTransport_
|
|
(
|
|
compressible::momentumTransportModel::New
|
|
(
|
|
rho,
|
|
U,
|
|
rhoPhi,
|
|
mixture
|
|
)
|
|
),
|
|
|
|
momentumTransport(momentumTransport_())
|
|
{
|
|
// Read the controls
|
|
readControls();
|
|
|
|
if (correctPhi)
|
|
{
|
|
rAU = new volScalarField
|
|
(
|
|
IOobject
|
|
(
|
|
"rAU",
|
|
runTime.name(),
|
|
mesh,
|
|
IOobject::READ_IF_PRESENT,
|
|
IOobject::AUTO_WRITE
|
|
),
|
|
mesh,
|
|
dimensionedScalar(dimTime/dimDensity, 1)
|
|
);
|
|
}
|
|
}
|
|
|
|
|
|
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
|
|
|
Foam::solvers::compressibleMultiphaseVoF::~compressibleMultiphaseVoF()
|
|
{}
|
|
|
|
|
|
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
|
|
|
|
void Foam::solvers::compressibleMultiphaseVoF::prePredictor()
|
|
{
|
|
multiphaseVoFSolver::prePredictor();
|
|
|
|
contErr = fvc::ddt(rho)()() + fvc::div(rhoPhi)()();
|
|
|
|
forAll(mixture.phases(), phasei)
|
|
{
|
|
const volScalarField& rho = phases[phasei].thermo().rho();
|
|
contErr.ref() -= fvModels().source(phases[phasei], rho)ρ
|
|
}
|
|
|
|
if (pimple.predictTransport())
|
|
{
|
|
momentumTransport.predict();
|
|
}
|
|
}
|
|
|
|
|
|
void Foam::solvers::compressibleMultiphaseVoF::postCorrector()
|
|
{
|
|
if (pimple.correctTransport())
|
|
{
|
|
momentumTransport.correct();
|
|
}
|
|
}
|
|
|
|
|
|
// ************************************************************************* //
|