songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

GIT: mapPatch.H: unused include files

Browse Source
This commit is contained in:
mattijs 2013-10-14 12:38:24 +01:00
parent 33e8d6f456
commit dc8d0658e6
2 changed files with 0 additions and 172 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

82
src/sampling/include/buildPatch.H
View File

@ -1,82 +0,0 @@
IOdictionary planeToPatchDict
(
IOobject
(
"planeToPatchDict",
runTime.system(),
runTime,
IOobject::MUST_READ_IF_MODIFIED,
IOobject::NO_WRITE
)
);
Switch planeMapping
(
planeToPatchDict.lookup("mappingPlane")
);
Info << "Constructing planePatch" << endl;
cuttingPlane cut(plane(planeToPatchDict), mesh);
primitiveFacePatch planePatch
(
cut.faces(),
cut.points()
);
Info << "Finding target patch" << endl;
word toPatchName
(
planeToPatchDict.lookup("targetPatchName")
);
const fvPatchList& patches = mesh.boundary();
bool targetPatchNamePresent(false);
label targetPatchNumber = -1;
forAll(patches, patchI)
{
if (patches[patchI].name() == toPatchName)
{
targetPatchNamePresent = true;
targetPatchNumber = patchI;
}
}
const labelList& cutCells(cut.cells());
if (!targetPatchNamePresent)
{
targetPatchNumber = 1;
Serr<< "Target patch not present. " << endl;
}
else
{
Serr<< "Target patch name: " << patches[targetPatchNumber].name() << endl;
}
patchToPatchInterpolation interPatch
(
planePatch,
patches[targetPatchNumber].patch()
);
PrimitivePatchInterpolation
<
PrimitivePatch<face, List, const pointField&>
> planePatchInter
(
planePatch
);
PrimitivePatchInterpolation
<
PrimitivePatch<face, List, const pointField&>
> toPatchInter
(
patches[targetPatchNumber].patch()
);

90
src/sampling/include/mapPatch.H
View File

@ -1,90 +0,0 @@
if (planeMapping && targetPatchNamePresent && cut.cut())
{
//Map variables U, phi & k
//Info<< typeInfo("isoLESmodel") << endl;
//Info<< "SGS model" << sgsModel.type() << endl;
scalarField planek(cutCells.size());
scalarField planeNuTilda(cutCells.size());
vectorField planeU(cutCells.size());
forAll(cutCells, cCellsI)
{
if (sgsModel.type() == "SpalartAllmaras")
{
planeNuTilda[cCellsI] = sgsModel.nuTilda()()[cutCells[cCellsI]];
}
else
{
planek[cCellsI] = sgsModel.k()()[cutCells[cCellsI]];
}
planeU[cCellsI] = U[cutCells[cCellsI]];
}
if (sgsModel.type() == "SpalartAllmaras")
{
Info<< "Mapping NuTilda." << endl;
sgsModel.nuTilda()().boundaryField()[targetPatchNumber] ==
interPatch.faceInterpolate(planeNuTilda)();
}
else
{
Info<< "Mapping k." << endl;
sgsModel.k()().boundaryField()[targetPatchNumber] ==
interPatch.faceInterpolate(planek)();
}
U.boundaryField()[targetPatchNumber] ==
interPatch.faceInterpolate(planeU)();
/* sgsModel.k()().boundaryField()[targetPatchNumber] ==
toPatchInter.pointToFaceInterpolate
(
interPatch.pointInterpolate
(
planePatchInter.faceToPointInterpolate
(
planek
)
)
);
U.boundaryField()[targetPatchNumber] ==
toPatchInter.pointToFaceInterpolate
(
interPatch.pointInterpolate
(
planePatchInter.faceToPointInterpolate
(
planeU
)
)
);
*/
scalar Q = sum
(
mesh.Sf().boundaryField()[targetPatchNumber]
& U.boundaryField()[targetPatchNumber]
);
scalar Qbar = sum
(
mesh.Sf().boundaryField()[targetPatchNumber] &
Ubar.value()
);
U.boundaryField()[targetPatchNumber] ==
U.boundaryField()[targetPatchNumber] * (Qbar/Q);
Info<< "Mass flux correction factor: " << (Qbar/Q) << endl;
phi.boundaryField()[targetPatchNumber] ==
(
mesh.Sf().boundaryField()[targetPatchNumber] &
U.boundaryField()[targetPatchNumber]
);
}