43 void Foam::wedgePolyPatch::initTransforms()
49 patchNormal_ = operator[](0).normal(points);
50 patchNormal_ /=
mag(patchNormal_);
55 sign(patchNormal_.
x())*(
max(
mag(patchNormal_.
x()), 0.5) - 0.5),
56 sign(patchNormal_.
y())*(
max(
mag(patchNormal_.
y()), 0.5) - 0.5),
57 sign(patchNormal_.
z())*(
max(
mag(patchNormal_.
z()), 0.5) - 0.5)
59 centreNormal_ /=
mag(centreNormal_);
63 mag(centreNormal_.
x() + centreNormal_.
y() + centreNormal_.
z())
69 "wedgePolyPatch::wedgePolyPatch(const polyPatch&, "
70 "const fvBoundaryMesh&)"
71 ) <<
"wedge " <<
name()
72 <<
" centre plane does not align with a coordinate plane by "
74 -
mag(centreNormal_.
x()+centreNormal_.
y()+centreNormal_.
z())
78 axis_ = centreNormal_ ^ patchNormal_;
79 scalar magAxis =
mag(axis_);
86 "wedgePolyPatch::initTransforms()"
87 ) <<
"wedge " <<
name()
88 <<
" plane aligns with a coordinate plane." <<
nl
89 <<
" The wedge plane should make a small angle (~2.5deg)"
90 " with the coordinate plane" <<
nl
91 <<
" and the the pair of wedge planes should be symmetric"
92 <<
" about the coordinate plane." <<
nl
93 <<
" Normal of face " << 0 <<
" is " << patchNormal_
94 <<
" , implied coordinate plane direction is " << centreNormal_
99 cellT_ = faceT_ & faceT_;
156 polyPatch(pp, bm, index, newSize, newStart)