Solve the nonlinear coupled-field problem with FEniCSX

2

Try making a mixed function space with u and y:

v_el = ufl.VectorElement("CG", mesh.ufl_cell(), 1)
mixed_el = ufl.MixedElement([v_el, v_el])
VP = dolfinx.FunctionSpace(mesh, mixed_el)
(v, w) = ufl.TestFunctions(VP)
q = dolfinx.Function(VP)
u, y = ufl.split(q)
1 Like