[javiersanchez]

I have a setup with 4 cylindrical electrodes into a liquid metal, and I must put a positive current in two of them, and a negative current at the others.

The mhdFoam (the solver for mhd systems) requires boundary conditions for U, p and B. The forcing I have is the Lorentz force JxB, but I was thought put this forcing as boundary condition in the electrodes walls, as J=rot(B).

However, I just have found boundary conditions like the Neumann conditions (value of normal gradient at the walls).

Does anyone know if we can put the curl of B value at the walls and the syntaxis in OpenFoam for that?

Thank you in advance.

Have a nice day

Javier