I recently came across this GitHub repository, and I found it quite interesting. After going through the code, I have a couple of questions:
It seems that the Governing Equations quite different from the references I. Steinbach 1999 you provided. I would appreciate it if you could provide some insights or clarification on these points.
Where did the various coefficients used in the code come from? Like u=64 and u_3phase_penalty=3.