FreeFem++ is an open source partial differential equation solver. freefem scripts can solve multiphysics non linear systems in 2D and 3D.