zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

solvers::compressibleVoF,compressibleVoF::correctCoNum(): Removed unnecessary nInternalFaces() test

Browse Source 
Resolves bug-report https://bugs.openfoam.org/view.php?id=3929
This commit is contained in:
Henry Weller
2022-11-10 14:30:21 +00:00
parent 745e4ff80d
commit e1c0207695
2 changed files with 20 additions and 32 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
applications/solvers/modules/fluid/compressibleVoF/compressibleVoF.C
View File

@ -45,22 +45,16 @@ void Foam::solvers::compressibleVoF::correctCoNum()
{
fluidSolver::correctCoNum(phi);
alphaCoNum = 0;
scalar meanAlphaCoNum = 0;
const scalarField sumPhi
(
mixture.nearInterface()().primitiveField()
*fvc::surfaceSum(mag(phi))().primitiveField()
);
if (mesh.nInternalFaces())
{
const scalarField sumPhi
(
mixture.nearInterface()().primitiveField()
*fvc::surfaceSum(mag(phi))().primitiveField()
);
alphaCoNum = 0.5*gMax(sumPhi/mesh.V().field())*runTime.deltaTValue();
alphaCoNum = 0.5*gMax(sumPhi/mesh.V().field())*runTime.deltaTValue();
meanAlphaCoNum =
const scalar meanAlphaCoNum =
0.5*(gSum(sumPhi)/gSum(mesh.V().field()))*runTime.deltaTValue();
}
Info<< "Interface Courant Number mean: " << meanAlphaCoNum
<< " max: " << alphaCoNum << endl;

32
applications/solvers/modules/fluid/multiphaseEuler/multiphaseEuler/multiphaseEuler.C
View File

@ -62,31 +62,25 @@ void Foam::solvers::multiphaseEuler::read()
void Foam::solvers::multiphaseEuler::correctCoNum()
{
CoNum = 0;
scalar meanCoNum = 0.0;
scalarField sumPhi
(
fvc::surfaceSum(mag(phi))().primitiveField()
);
if (mesh.nInternalFaces())
forAll(phases, phasei)
{
scalarField sumPhi
sumPhi = max
(
fvc::surfaceSum(mag(phi))().primitiveField()
sumPhi,
fvc::surfaceSum(mag(phases[phasei].phi()))().primitiveField()
);
forAll(phases, phasei)
{
sumPhi = max
(
sumPhi,
fvc::surfaceSum(mag(phases[phasei].phi()))().primitiveField()
);
}
CoNum = 0.5*gMax(sumPhi/mesh.V().field())*runTime.deltaTValue();
meanCoNum =
0.5*(gSum(sumPhi)/gSum(mesh.V().field()))*runTime.deltaTValue();
}
CoNum = 0.5*gMax(sumPhi/mesh.V().field())*runTime.deltaTValue();
const scalar meanCoNum =
0.5*(gSum(sumPhi)/gSum(mesh.V().field()))*runTime.deltaTValue();
Info<< "Courant Number mean: " << meanCoNum
<< " max: " << CoNum << endl;
}