mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
Merge branch 'master' of /home/dm4/OpenFOAM/repositories/OpenFOAM-dev
This commit is contained in:
@ -1,36 +0,0 @@
|
|||||||
if (packingLimiter)
|
|
||||||
{
|
|
||||||
// Calculating exceeding volume fractions
|
|
||||||
volScalarField alphaEx(max(alpha1 - alphaMax, scalar(0)));
|
|
||||||
|
|
||||||
// Finding neighbouring cells of the whole domain
|
|
||||||
labelListList neighbour = mesh.cellCells();
|
|
||||||
scalarField cellVolumes(mesh.cellVolumes());
|
|
||||||
|
|
||||||
forAll (alphaEx, celli)
|
|
||||||
{
|
|
||||||
// Finding the labels of the neighbouring cells
|
|
||||||
labelList neighbourCell = neighbour[celli];
|
|
||||||
|
|
||||||
// Initializing neighbouring cells contribution
|
|
||||||
scalar neighboursEx = 0.0;
|
|
||||||
|
|
||||||
forAll (neighbourCell, cellj)
|
|
||||||
{
|
|
||||||
labelList neighboursNeighbour = neighbour[neighbourCell[cellj]];
|
|
||||||
scalar neighboursNeighbourCellVolumes = 0.0;
|
|
||||||
|
|
||||||
forAll (neighboursNeighbour, cellk)
|
|
||||||
{
|
|
||||||
neighboursNeighbourCellVolumes +=
|
|
||||||
cellVolumes[neighboursNeighbour[cellk]];
|
|
||||||
}
|
|
||||||
|
|
||||||
neighboursEx +=
|
|
||||||
alphaEx[neighbourCell[cellj]]*cellVolumes[celli]
|
|
||||||
/neighboursNeighbourCellVolumes;
|
|
||||||
}
|
|
||||||
|
|
||||||
alpha1[celli] += neighboursEx - alphaEx[celli];
|
|
||||||
}
|
|
||||||
}
|
|
||||||
@ -1,36 +0,0 @@
|
|||||||
IOdictionary ppProperties
|
|
||||||
(
|
|
||||||
IOobject
|
|
||||||
(
|
|
||||||
"ppProperties",
|
|
||||||
runTime.constant(),
|
|
||||||
mesh,
|
|
||||||
IOobject::MUST_READ,
|
|
||||||
IOobject::NO_WRITE
|
|
||||||
)
|
|
||||||
);
|
|
||||||
|
|
||||||
scalar preAlphaExp
|
|
||||||
(
|
|
||||||
readScalar(ppProperties.lookup("preAlphaExp"))
|
|
||||||
);
|
|
||||||
|
|
||||||
scalar alphaMax
|
|
||||||
(
|
|
||||||
readScalar(ppProperties.lookup("alphaMax"))
|
|
||||||
);
|
|
||||||
|
|
||||||
scalar expMax
|
|
||||||
(
|
|
||||||
readScalar(ppProperties.lookup("expMax"))
|
|
||||||
);
|
|
||||||
|
|
||||||
dimensionedScalar g0
|
|
||||||
(
|
|
||||||
ppProperties.lookup("g0")
|
|
||||||
);
|
|
||||||
|
|
||||||
Switch packingLimiter
|
|
||||||
(
|
|
||||||
ppProperties.lookup("packingLimiter")
|
|
||||||
);
|
|
||||||
@ -505,7 +505,7 @@ void Foam::cyclicPolyPatch::getCentresAndAnchors
|
|||||||
}
|
}
|
||||||
default:
|
default:
|
||||||
{
|
{
|
||||||
// Assumes that cyclic is planar. This is also the initial
|
// Assumes that cyclic is rotational. This is also the initial
|
||||||
// condition for patches without faces.
|
// condition for patches without faces.
|
||||||
|
|
||||||
// Determine the face with max area on both halves. These
|
// Determine the face with max area on both halves. These
|
||||||
|
|||||||
@ -187,7 +187,7 @@ void Foam::cyclicAMIPolyPatch::calcTransforms
|
|||||||
}
|
}
|
||||||
|
|
||||||
scalar areaError =
|
scalar areaError =
|
||||||
max(errorPos, errorNeg)/(mag(area1) + ROOTVSMALL);
|
min(errorPos, errorNeg)/(mag(area1) + ROOTVSMALL);
|
||||||
|
|
||||||
if (areaError > matchTolerance())
|
if (areaError > matchTolerance())
|
||||||
{
|
{
|
||||||
|
|||||||
Reference in New Issue
Block a user