Reply To: Curved boundary
OpenLB – Open Source Lattice Boltzmann Code › Forums › on 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.
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.