Skip to content

Reply To: Curved boundary

#9493
Bran
Participant

Still doesn’t work, only this part, can’t be setWallfunctionboundary.
Hump
Following is how I create the geometry:

STLreader<T> Hump( "Hump.stl", converter.getConversionFactorLength(), 1, 1);
CuboidGeometry3D<T> cuboidGeometry(Hump, converter.getPhysDeltaX(), noOfCuboids);
cuboidGeometry.setPeriodicity(true, false, true);
HeuristicLoadBalancer<T> loadBalancer(cuboidGeometry);
SuperGeometry<T,3> superGeometry(cuboidGeometry, loadBalancer);
prepareGeometry(superGeometry, Hump, converter);
void prepareGeometry(SuperGeometry<T,3>& superGeometry,
                     IndicatorF3D<T>& Hump,
                     UnitConverter<T,DESCRIPTOR> const& converter)
{
  OstreamManager clout(std::cout,"prepareGeometry");
  clout << "Prepare Geometry ..." << std::endl;
  superGeometry.rename(0,2,Hump);
  superGeometry.rename(2,1,{0,1,0});
  superGeometry.clean();
  superGeometry.innerClean();
  superGeometry.checkForErrors();
  superGeometry.print();
  olb::Vector<T, 3> PhyMax = superGeometry.getStatistics().getMaxPhysR(2);
  olb::Vector<T, 3> PhyMin = superGeometry.getStatistics().getMinPhysR(2);
  clout << "Dimension of the channel in meters: x = " << PhyMax[0] - PhyMin[0];
  clout << " ; y = " << PhyMax[1] - PhyMin[1];
  clout << " ; z = " << PhyMax[2] - PhyMin[2] << std::endl;
  clout << "Prepare Geometry ... OK" << std::endl;
}
  • This reply was modified 4 weeks, 1 day ago by Bran.