[Honey]

Hi,

I am willing to create a tapered cylinder in middle of a box with using topoSetDict. I know that is not implemented in OpenFOAM and requires some modification into cylinderToCell.C code located at .../src/meshTools/sets/cellSources/cylinderToCell. Nevertheless, I tried to modify the main function written in cylinderToCell.C from:
void Foam::cylinderToCell::combine(topoSet& set, const bool add) const
{
const vector axis = p2_ - p1_;
const scalar rad2 = sqr(radius_);
const scalar magAxis2 = magSqr(axis);
const pointField& ctrs = mesh_.cellCentres();

forAll(ctrs, cellI)
{
vector d = ctrs[cellI] - p1_;
scalar magD = d & axis;

if ((magD > 0) && (magD < magAxis2))
{
scalar d2 = (d & d) - sqr(magD)/magAxis2;
if (d2 < rad2)
{
addOrDelete(set, cellI, add);
}}}}

to:

void Foam::taperedCylinderToCell::combine(topoSet& set, const bool add) const
{
const vector axis = p2_ - p1_;
const scalar radBase2 = sqr(radiusBase_); //**modified
const scalar radTop2 = sqr(radiusTop_); //**added
const scalar magAxis2 = magSqr(axis);
const pointField& ctrs = mesh_.cellCentres();

forAll(ctrs, cellI)
{
vector d = ctrs[cellI] - p1_;
scalar magD = d & axis;

if ((magD > 0) && (magD < magAxis2))
{
scalar d2 = (d & d) - sqr(magD)/magAxis2;
if ((d2 < radBase2) && (d2 > radTop2)) //**modified
{
addOrDelete(set, cellI, add);
}}}}

Obviously, it generated a hollow cylinder instead of a tapered cylinder I would highly appreciate if anyone here could help me to successfully modify this piece of code!

Thank you in advance