Skip to content

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

#8921
jan
Participant

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.