Hi, I think it is rather an issue of implementation to get the terms contributing to the off-diagonal blocks in the mass matrix right. Essentially its the velocity basis functions in an inner product with the bases coming from the pressure space.
Best, Martin On Tue, Feb 1, 2011 at 2:35 PM, Markus Bürg <[email protected]> wrote: > Hello Martin, > > as far as I can oversee this issue, I think there is some mistake in the > derivation of the variational formulation. You cannot add a scalar value to > a tensor. Did you derive the variational formulation in another way than > tutorial step-22 does? > > Best Regards, > Markus > > > > Am 01.02.11 14:18, schrieb Martin Stoll: > > >> Hi, >> >> thank you for pointing this out what I obviously meant was >> "...phi_u[j[ will be of Tensor element and phi_p[i] a scalar..." >> >> This comes from the discretization of a control problem where a Q1 control >> "meets" a Q2 test function. >> >> Sorry about the confusion. >> >> Best, >> Martin >> > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > -- *Martin Stoll* *Postdoctoral Research Fellow* Computational Methods in Systems and Control Theory Max Planck Institute for Dynamics of Complex Technical Systems Sandtorstr. 1 D-39106 Magdeburg Germany Email: [email protected] URL : http://www.mpi-magdeburg.mpg.de/people/stollm Tel :+49 391 6110 384
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
