Reply To: How to add 2d vtk files as geometrical models
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 › How to add 2d vtk files as geometrical models › Reply To: How to add 2d vtk files as geometrical models
July 17, 2024 at 1:11 pm
#8947
jan
Participant
Dear Bobbie,
if your 2D structure can be converted to a VTI, then you can use an approach such as this, which uses the IndicatorBlockData2D:
std::shared_ptr<BlockVTIreader3D<T, T>> vtiReader(
new BlockVTIreader3D<T, T>(vtiFile, "Elevation"));
vtiReader->printInfo();
auto cuboidSample = vtiReader->getCuboid();
T deltaRsample = cuboidSample.getDeltaR() * sourceScale;
Vector<int, 3> extentSample = cuboidSample.getExtent();
Vector<T, 3> originSamplePhys = cuboidSample.getOrigin() * sourceScale;
Vector<T, 3> extentSamplePhys = {deltaRsample * T(extentSample[0] + 0.5),
deltaRsample * T(extentSample[1] + 0.5),
deltaRsample * T(extentSample[2] + 0.5)};
std::shared_ptr<IndicatorBlockData2D<T>> ind(
new IndicatorBlockData2D<T>(vtiReader->getBlockData(), extentSamplePhys,
originSamplePhys, deltaRsample, false));
Best regards,
Jan
