About MHD equations

You can do as


The penalty formulation (3) will work if the boundary is made of flat pieces. Otherwise you have to lump it, see Boundary condition - #6 by fb77