Dirichlet and Boundary conditions in Darcy tutorial.

Hi @fverdugo ... I think that the Dirichlet and Neumann Boundary conditions are flipped in the Darcy's tutorial + RT FEs.

In particular, in the Darcy problem, one imposes the trace of the pressure in the Dirichlet boundary, and the normal component of the flux in the Neumann boundary. Indeed, Dirichlet boundary conditions are natural boundary conditions (i.e., imposed weakly) for Darcy, and Neumann boundary condition are essential boundary conditions for Darcy (i.e., imposed strongly).

See, e.g., pag. 153 of

https://team-pancho.github.io/documents/anIntro2FEM_2015.pdf

Is there anything I am missing here? Do you agree? I can help with the changes required to fix this, if needed.

Locate the Darcy tutorial and its Raviart–Thomas finite-element examples, then compare the boundary-condition descriptions and implementation with the issue's cited FEM reference. Confirm which conditions are labeled or applied incorrectly; done means the tutorial's explanation and affected examples consistently match the intended Darcy formulation.

Written by the indexing model from the issue text.