On 3/11/21 10:14 PM, [email protected] wrote:
I am currently trying to solve an advection-diffusion equation at very, very
high Péclet number. I first implemented a GLS-stabilized continuous Galerkin
solution, it works well, but clearly I think a DG approach would be immensely
better here.
The issue is that the velocity field for my advection-diffusion stems from a
Continuous Galerkin solution with another dof handler, formulated using a
FESystem. Everything work well until them. The issue I have is that the
FEInterfaceValues does not seem to want to work with FEValuesExtractors? It
seems like it cannot use the [] operator to extract my velocity values out of
my 4 components (velocity + pressure).
I am sure i'm forgetting something stupid here, but how should I proceed for
vector valued problem using DG?
No, you're not forgetting anything. FEInterfaceValues is still a pretty new
class and nobody has ever written the equivalent of the FEValuesViews classes
that you get from FEValues when you apply an extractor.
I suspect that the effort would not be huge, and that one could just
copy-paste a lot of code from the existing views classes. If you were
interested in doing that, I'd be happy to discuss what I think is going to be
necessary. But the short answer is really that FEInterfaceValues isn't ready
for vector-valued problems at the moment.
Best
Wolfgang
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/95dbdc73-3b8b-6402-98b9-4894446aa55b%40colostate.edu.
