multiphase VoF solvers: Ensure the nHatf surfaceField in interfaceProperties is updated following mesh change

This commit is contained in:
Henry Weller
2019-03-19 20:13:40 +00:00
parent f0c3c497e7
commit aaf4a0bfcc
8 changed files with 30 additions and 15 deletions

View File

@ -125,10 +125,10 @@ int main(int argc, char *argv[])
// Make the fluxes relative to the mesh motion
fvc::makeRelative(phi, U);
mixture.correct();
}
mixture.correct();
if (checkMeshCourantNo)
{
#include "meshCourantNo.H"

View File

@ -123,10 +123,10 @@ int main(int argc, char *argv[])
// Make the flux relative to the mesh motion
fvc::makeRelative(phi, U);
mixture.correct();
}
mixture.correct();
if (checkMeshCourantNo)
{
#include "meshCourantNo.H"

View File

@ -115,10 +115,10 @@ int main(int argc, char *argv[])
// Make the flux relative to the mesh motion
fvc::makeRelative(phi, U);
mixture.correct();
}
mixture.correct();
if (checkMeshCourantNo)
{
#include "meshCourantNo.H"

View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -78,6 +78,7 @@ Foam::phaseChangeTwoPhaseMixtures::Kunz::mDotAlphal() const
);
}
Foam::Pair<Foam::tmp<Foam::volScalarField>>
Foam::phaseChangeTwoPhaseMixtures::Kunz::mDotP() const
{
@ -95,7 +96,9 @@ Foam::phaseChangeTwoPhaseMixtures::Kunz::mDotP() const
void Foam::phaseChangeTwoPhaseMixtures::Kunz::correct()
{}
{
phaseChangeTwoPhaseMixture::correct();
}
bool Foam::phaseChangeTwoPhaseMixtures::Kunz::read()

View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -75,6 +75,7 @@ Foam::phaseChangeTwoPhaseMixtures::Merkle::mDotAlphal() const
);
}
Foam::Pair<Foam::tmp<Foam::volScalarField>>
Foam::phaseChangeTwoPhaseMixtures::Merkle::mDotP() const
{
@ -90,7 +91,9 @@ Foam::phaseChangeTwoPhaseMixtures::Merkle::mDotP() const
void Foam::phaseChangeTwoPhaseMixtures::Merkle::correct()
{}
{
phaseChangeTwoPhaseMixture::correct();
}
bool Foam::phaseChangeTwoPhaseMixtures::Merkle::read()

View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -144,7 +144,9 @@ Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer::mDotP() const
void Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer::correct()
{}
{
phaseChangeTwoPhaseMixture::correct();
}
bool Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer::read()

View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -63,6 +63,7 @@ Foam::phaseChangeTwoPhaseMixture::vDotAlphal() const
);
}
Foam::Pair<Foam::tmp<Foam::volScalarField>>
Foam::phaseChangeTwoPhaseMixture::vDotP() const
{
@ -73,6 +74,12 @@ Foam::phaseChangeTwoPhaseMixture::vDotP() const
}
void Foam::phaseChangeTwoPhaseMixture::correct()
{
immiscibleIncompressibleTwoPhaseMixture::correct();
}
bool Foam::phaseChangeTwoPhaseMixture::read()
{
if (immiscibleIncompressibleTwoPhaseMixture::read())

View File

@ -108,10 +108,10 @@ int main(int argc, char *argv[])
// Make the flux relative to the mesh motion
fvc::makeRelative(phi, U);
mixture.correct();
}
mixture.correct();
if (checkMeshCourantNo)
{
#include "meshCourantNo.H"