“…The constraint operator (see (3)) may become: the trace operator under contact conditions [1][2][3], the jump operator for cracks and anticracks [4][5][6], the gradient operator in plasticity [7], the divergence operator under incompressibility conditions [8][9][10], a state-constraint in mathematical programs with equilibrium constraints [11,12], and the like. The constraint problems are related to parameter identification problems (see the theory in References [13][14][15] and application to biological systems in Reference [16]), to inverse problems by the mean of observation data used in mathematical physics [17,18] and in acoustics [19][20][21], to overdetermined and free-boundary problems [22,23]. As an application, in the current paper we focus on the incompressible Brinkman flow problem under a divergence-free condition (see the related modeling of porous medium in References [24,25], well-posedness analysis in Reference [26], and fluid-porous coupling with numerics in References [27,28]).…”