Solving a system of coupled PDEs in fenicsx

[dokken]

@Mamadou

right_m_dofs = locate_dofs_geometrical(Q, right)
bc_m_right = dirichletbc(PETSc.ScalarType(0.12), right_m_dofs, Q)



top_m_dofs = locate_dofs_geometrical(Q, top)
bc_m_top = dirichletbc(PETSc.ScalarType(0.12), top_m_dofs, Q)

right_T_dofs = locate_dofs_geometrical(V, right)
bc_T_right = dirichletbc(PETSc.ScalarType(80.), right_T_dofs, V)

top_T_dofs = locate_dofs_geometrical(V, top)
bc_T_top = dirichletbc(PETSc.ScalarType(80.), top_T_dofs, V)

#Creation of dirichlet conditions
#bcm = [bc_m_right, bc_m_top]     #bc for unknown m
#bcT = [bc_T_right, bc_T_top]     #bc for unknown T

bcs = [bc_m_right, bc_m_top, bc_T_right, bc_T_top]  #dirichlet boundary conditions

None of the boundary conditions are created correctly to work for a problem with Mixed elements. See:
https://jsdokken.com/dolfinx-tutorial/chapter3/component_bc.html
Or Solving Stokes equation with NewtonSolver - #4 by dokken