Reply To: Problems of Resolved Particles with Periodic Boundary
OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › General Topics › Problems of Resolved Particles with Periodic Boundary › Reply To: Problems of Resolved Particles with Periodic Boundary
Dear yueyq,
solidBoundaries.push_back( SolidBoundary<T, DESCRIPTOR::d>(
std::make_unique<IndicInverse<T, DESCRIPTOR::d>>(
cuboid, cuboid.getMin() – 5 * converter.getPhysDeltaX(),
cuboid.getMax() + 5 * converter.getPhysDeltaX()), 2, wallContactMaterial));
the above part takes the inverse of the whole computation domain, which means each side is regarded as a wall. If you solely want the bottom to act as a wall, then I suggest that you introduce a IndicatorCuboid
instead and place it at the bottom wall. Be sure to place it exactly where the bottom wall is and give it a sufficient thickness, because the particles will overlap it.
That’s also part of the reason why the SolidBoundary
in the example is extended. If you’d have a wall with a infinitesimal thickness, it would lead to problems during the contact treatment, as we want to calculate the overlap volume.
Best regards,
Jan
- This reply was modified 7 months ago by jan.