Unable to Apply no_penetration LCBC

[YasamanNemati]

Hello SfePy Community,

I am working on modeling the effects of veins on a resistivity tool in a borehole. For simplification, I'm starting with a basic cube made from a single material. In my model, I apply a potential to the two ends of the cylinder and aim to enforce a no-penetration boundary condition. However, I'm encountering difficulties in correctly applying this linear combination boundary condition (LCBC).
You can view my full script and mesh here:
code-and-mesh.zip

this is the error I receive:
ValueError: assertion failed!

Could this be an issue with how I am specifying the DOFs or regions for the no penetration LCBC? Or is there something else I'm missing in my LCBC and region definitions? Also, can you please let me know if I need to add anything to my equation in order to consider the defined lcbc? I'm so confused.

Any help troubleshooting would be greatly appreciated.

Thanks!

[rc]

I have not yet had time to look at your code, but note that specifying LCBCs can be tricky, so there is an alternative - a penalty term, see [1]. Maybe that would work for you?

[1] http://sfepy.org/doc-devel/examples/navier_stokes-stokes_slip_bc.html

[rc]

Also compare [2], [3].

[2] http://sfepy.org/doc-devel/examples/multi_physics-biot_npbc.html
[3] http://sfepy.org/doc-devel/examples/multi_physics-biot_npbc_lagrange.html

[YasamanNemati]

Thank you so much. I tried to follow the steps in these examples, but I'm still unable to apply the no_penetration condition for my model. Do you have any other suggestions? I'd really appreciate it.