Incompressible Navier-Stokes with particular scheme & mixed elements FEniCSx

[dokken]

@TOP1

bcu_inflow = dirichletbc(u_inlet, locate_dofs_topological(V, fdim, ft.find(inlet_marker)))

# Walls
u_nonslip = np.array((0,) * mesh.geometry.dim, dtype=PETSc.ScalarType)
bcu_walls = dirichletbc(u_nonslip, locate_dofs_topological(V, fdim, ft.find(wall_marker)), V)

# Obstacle
bcu_obstacle = dirichletbc(u_nonslip, locate_dofs_topological(V, fdim, ft.find(obstacle_marker)), V)

These boundary conditions are not correct for mixed spaces, see for instance
https://jsdokken.com/dolfinx-tutorial/chapter3/component_bc.html#boundary-conditions