Constrained Optimization

I am trying to impose a total energy dissipation constraint on a fluid flow topology optimization problem.

The equation is: ∅=∫(τ:S+α(u.u))dV

Where α is a function of the density variable.

I added the expression (intop1(spf.Qvd+alpha*(u^2+v^2))) in the constraints list in the optimization node and since the expression is in the integral form, I guess it is a global expression.

I evaluate the total energy dissipation after each optimization and the value is always way above the upper bound. How can I impose this constraint on my problem?