“…Taking the recently developed Finite Element Method -Flux Corrected Transport Scheme (FEM-FCT) [1,2] as starting point, we focus on a partial differential-algebraic framework that allows a general coupling procedure where boundary conditions are expressed as constraints and appended by means of Lagrange multipliers. The stabilization approach that we provide makes it possible to model and stabilize a flow problem as part of a more general coupled system that may stem from fluid-structure interaction, from domain decomposition or from other multiphysics models, see, e.g., [3,4,5,6,7,8] for such applications.…”