FEniCSx or Multiphenicsx for Nonlinear coupled PDE system

Tutorial 03 of multiphenicsx should give you a good starting point on a much simpler problem.

The next release of dolfinx will have codimension 1 support for multiphysics problems. The implementation is completely different from the one in multiphenicsx, and therefore the builtin implementation of dolfinx and the one of multiphenicsx will coexist for the foreseeable future.