Or, you could use as a stable pair of finite elements the FE_Q(k)-FE_DGP(k-1) pair:
D. Boffi and L. Gastaldi. On the quadrilateral Q2-P1 element for the Stokes problem. Int. J. Numer. Meth. Fluids, 39 (2002), 1001-1011. This pair of finite elements is known to produce consistently better results for discontinuous pressures (after all, p is only supposed to be in L2). L. > On 7 Apr 2017, at 23:51, Sumedh Yadav <sumedhyada...@gmail.com> wrote: > > > I would like to add an idea I intend to apply but right now am clueless of > how to apply. > In the earlier post I mentioned issue of 'spurious currents'. In the past > researchers have tried quite a few approaches to handle this. One of them > involves reconstruction of computed pressure field. In essence if you look > carefully at the images I shared in the last post you can see the failure of > taylor-hood elements to handle the discontinuity in pressure which in turn > makes the velocity field oscillatory/spurious. Reconstruction of pressure is > done only in the interface region by interpolating the local normal direction > values of pressure inside and outside this interface region. I am able to > construct the normal field necessary for this local reconstruction. But I > need to traverse to neighboring cells in both directions of this normal > vector in the interface cells. I am clueless to how this can be done in > deal.ii. The normal field image - > > > > -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to dealii+unsubscr...@googlegroups.com. > For more options, visit https://groups.google.com/d/optout. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
