zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

twoPhaseEulerFoam: Corrected potential conservation issue at inlet-outlet BCs

Browse Source 
Updated MULES to ensure the limiter is 1 at inlet BCs
This commit is contained in:
Henry
2013-08-15 17:03:52 +01:00
parent b9f9598725
commit 13e1addc35
4 changed files with 12 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H
View File

@ -92,6 +92,12 @@
)
);
alphaPhic1.boundaryField() = min
(
phi1.boundaryField()*alpha1.boundaryField(),
alphaPhic1.boundaryField()
);
MULES::explicitSolve
(
geometricOneField(),

5
src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
View File

@ -485,11 +485,12 @@ void Foam::MULES::limiter
const labelList& pFaceCells =
mesh.boundary()[patchi].faceCells();
const scalarField& phiBDPf = phiBDBf[patchi];
const scalarField& phiCorrPf = phiCorrBf[patchi];
forAll(lambdaPf, pFacei)
{
// Limit outlet faces only
if (phiBDPf[pFacei] > 0)
if ((phiBDPf[pFacei] + phiCorrPf[pFacei]) > SMALL*SMALL)
{
label pfCelli = pFaceCells[pFacei];
@ -862,7 +863,7 @@ void Foam::MULES::limiterCorr
forAll(lambdaPf, pFacei)
{
// Limit outlet faces only
if (phiCorrPf[pFacei] > 0)
if (phiCorrPf[pFacei] > SMALL*SMALL)
{
label pfCelli = pFaceCells[pFacei];

5
tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/0/U.particles
View File

@ -28,9 +28,8 @@ boundaryField
outlet
{
type pressureInletOutletVelocity;
phi phi.particles;
value $internalField;
type fixedValue;
value uniform (0 0 0);
}
walls

3
tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/0/alpha.particles
View File

@ -6027,8 +6027,7 @@ boundaryField
{
inlet
{
type fixedValue;
value uniform 0;
type zeroGradient;
}
outlet
{