Hi,
I am trying to implement a coupled Stokes-Darcy problem using Step-46
tutorial as an example. So, essentially I would want to use something like
the following
FESystem<dim> fe_stokes(FE_Q<dim>(2), dim, // Stokes pair
FE_Q<dim>(1), 1,
FE_Nothing<dim>(dim), 1, // Darcy pair
FE_Nothing<dim>(), 1,
FE_Nothing<dim>(), 1); // Lagrange multiplier
FESystem<dim> fe_darcy(FE_Nothing<dim>(), dim,
FE_Nothing<dim>(), 1,
FE_RaviartThomas<dim>(0), 1,
FE_DGQ<dim>(0), 1,
FE_DGQ<dim>(0), 1);
hp::FECollection<dim> fe;
fe.push_back(fe_stokes);
fe.push_back(fe_darcy);
//...
dof_handler.distribute_dofs (fe);
But this results in the error message:
--------------------------------------------------------
An error occurred in line <936> of file
</home/eldar/TmpDeal/source/fe/fe.cc> in function
dealii::FiniteElementDomination::Domination
dealii::FiniteElement<<anonymous>, <anonymous>
>::compare_for_face_domination(const dealii::FiniteElement<<anonymous>,
<anonymous> >&) const [with int dim = 2; int spacedim = 2]
The violated condition was:
false
Additional information:
You are trying to use functionality in deal.II that is currently not
implemented. In many cases, this indicates that there simply didn't appear
much of a need for it, or that the author of the original code did not have
the time to implement a particular case. If you hit this exception, it is
therefore worth the time to look into the code to find out whether you may
be able to implement the missing functionality. If you do, please consider
providing a patch to the deal.II development sources (see the deal.II
website on how to contribute).
Stacktrace:
-----------
#0 /home/eldar/Libs/DealII_8.5.1/lib/libdeal_II.g.so.9.0.0-pre:
dealii::FiniteElement<2,
2>::compare_for_face_domination(dealii::FiniteElement<2, 2> const&) const
#1 /home/eldar/Libs/DealII_8.5.1/lib/libdeal_II.g.so.9.0.0-pre:
dealii::FESystem<2,
2>::compare_for_face_domination(dealii::FiniteElement<2, 2> const&) const
>From what I could see the implementation of dealii::FiniteElement<2,
2>::compare_for_face_domination(dealii::FiniteElement<2, 2> const&) const
consists of Assert(false, ...) so for some reason what I am trying to do
was not meant to be done this way.
Could you please suggest what else could be done, or, maybe, there is just
another way of approaching coupled problems like this one using Deal.II?
Thank you!
--
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].
For more options, visit https://groups.google.com/d/optout.
