Coding in openfoam is quite trivial at least to me, you have plenty of example source code things becomes difficult (for me) when you really have to implement a separate algorithm like additional Newton Raphson loop, but adding terms to equation should not be that hard right? Or atleast I would think so, at last it will take some time. More general pattern recognition than coding.