I am working on the following Brinkman–Forchheimer-Darcy problem
The deficiency can be circumvented with the Brinkman–Forchheimer-Darcy law for flows in closely packed media, which leads to the in the attached model.
I am looking for a Matlab code for solving the Brinkman–Forchheimer-Darcy equations by the finite element method .