Reply To: Curved boundary

[Bran]

Do you mean doing these steps? It will make the whole domain become material 0.

STLreader<T> Hump("Hump.stl", converter.getConversionFactorLength(), 1, 1);
IndicatorLayer3D<T> extendedDomain(Hump, converter.getConversionFactorLength());
CuboidGeometry3D<T> cuboidGeometry(extendedDomain, converter.getPhysDeltaX(), noOfCuboids);
cuboidGeometry.setPeriodicity(true, false, true);
HeuristicLoadBalancer<T> loadBalancer(cuboidGeometry);
SuperGeometry<T,3> superGeometry(cuboidGeometry, loadBalancer);
prepareGeometry(superGeometry, extendedDomain, Hump, converter);

void prepareGeometry(SuperGeometry<T,3>& superGeometry,
                     IndicatorF3D<T>& extendedDomain,
                     STLreader<T>& Hump,
                     UnitConverter<T,DESCRIPTOR> const& converter)
{
  OstreamManager clout(std::cout,"prepareGeometry");
  clout << "Prepare Geometry ..." << std::endl;
  superGeometry.rename(0,2,extendedDomain);
  superGeometry.rename(1,2,Hump);
  superGeometry.clean();
  superGeometry.innerClean();
  superGeometry.checkForErrors();
  superGeometry.print();
  clout << "Prepare Geometry ... OK" << std::endl;
}

This is how my stlfile looks like.

• This reply was modified 3 weeks, 3 days ago by Bran.