The finite element discretization in space, using the standard Galerkin procedure yields the following system of non linear, non symmetric and coupled system of equations.
