I'm interested in numerically modelling viscoelastic fluids in a 2D channel geometry, driven by a pressure drop. The method consists of solving coupled constitutive equation: (a) Stokes force balance (for incompressible flow) and (b) a polymeric model equation, e.g., Oldroyd-B or similar, which contributes an additional viscoelastic stress into (a). Eqn (a) is subject to no-slip and no-permeation boundary conditions (BCs). The pressure drop is in the x-direction, and the wall-normal direction is y. The details of the viscoelastic component is not necessary for the question - it is mentioned (i) to provide context for my problem, and (ii) as it means that solutions beyond simple Pouseille flow can be found.
Intuitively it seems as though solutions with a net y-velocity component should not allowed as this implies that the sample as a whole is moving vertically up or down in the channel, violating the no-permeation BCs. Vortex solutions where positive vy on one side of the vortex is cancelled out by negative vy on the other are clearly OK.
However I am struggling to show this analytically. Does anyone have any advice for how to show this if true, or can they provide a counter-example if not.
Many thanks!