Solving 1D non linear differential equation on complex geometry

Geometries can be loaded using meshio, see Transitioning from mesh.xml to mesh.xdmf, from dolfin-convert to meshio for many examples on how to load the mesh. For locally conservative methods, discontinuous Galerkin can be used, see for instance: https://bitbucket.org/fenics-project/dolfin/src/master/python/demo/undocumented/dg-advection-diffusion/demo_dg-advection-diffusion.py