Skip to content

Reply To: Curved boundary

Due to recent bot attacks we have changed the sign-up process. If you want to participate in our forum, first register on this website and then send a message via our contact form.

› Forums › OpenLB › General Topics › Curved boundary › Reply To: Curved boundary

November 7, 2024 at 3:36 am #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 1 year, 2 months ago by Bran.