xfygogo
/
phasicFlow
mirror of https://github.com/PhasicFlow/phasicFlow.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

zResolution is added to cylinder

This commit is contained in:
hamidrezanorouzi 2022-11-24 23:13:15 +03:30
parent 27e3e24cc3
commit 02f7f9af20
1 changed files with 23 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
utilities/geometryPhasicFlow/cylinderWall/cylinderWall.C
View File

@ -1,4 +1,5 @@
#include "cylinderWall.H" #include "cylinderWall.H"
#include "line.H"
bool pFlow::cylinderWall::readCylinderWall(const dictionary& dict) bool pFlow::cylinderWall::readCylinderWall(const dictionary& dict)
@ -9,9 +10,29 @@ bool pFlow::cylinderWall::readCylinderWall(const dictionary& dict)
auto radius2 = dict.getVal<real>("radius2") ; auto radius2 = dict.getVal<real>("radius2") ;
int32 resolution = dict.getValOrSet("resolution", 24 ); int32 resolution = dict.getValOrSet("resolution", 24 );
int32 zResolution = dict.getValOrSet("zResolution", 1);
return createCylinder(p1, p2, radius1, radius2, resolution); triangles_.clear();
triangles_.reserve(2*resolution*zResolution);
line cylAxis(p1, p2);
auto lp1 = p1;
auto dt = static_cast<real>(1.0/zResolution);
real t = 0;
for(int32 i=0; i<zResolution; i++)
{
t += dt;
auto lp2 = cylAxis.point(t);
if(!createCylinder(lp1, lp2, radius1, radius2, resolution))
return false;
lp1 = lp2;
}
return true;
} }
bool pFlow::cylinderWall::createCylinder( bool pFlow::cylinderWall::createCylinder(
@ -26,9 +47,6 @@ bool pFlow::cylinderWall::createCylinder(
real L = zAx.length(); real L = zAx.length();
// number of wall elements will be twice numDiv
triangles_.clear();
triangles_.reserve(2 * numDiv);
realx3Vector r1P(numDiv + 1), r2P(numDiv + 1); realx3Vector r1P(numDiv + 1), r2P(numDiv + 1);
real dTheta = 2 * Pi / numDiv; real dTheta = 2 * Pi / numDiv;