dolfin team mailing list archive

[Alexander Jarosch <alexanj@xxxxx>]

it is of course because in:

a = nu*(dot(grad(U), grad(v)) - P*div(v) + div(U)*q)*dx

I multiply the viscosity nu with the whole equation, which is wrong, so 
the pressure is wrong, but the thing which does not work is an approach like

a = (nu*dot(grad(U), grad(v)) - P*div(v) + nu*div(U)*q)*dx

so that the pressure is fine. I need nu as a function in my form, 
because nu changes over the domain.

Alex