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 -
<https://lh3.googleusercontent.com/-WQV5SbLTHLM/WOgI3JptIjI/AAAAAAAAFRs/9Q2J2STP8YMKzb7T2210RcT6S1Q0z1rfgCLcB/s1600/visit0000.png> -- 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.