In brief the steps I am using are:
* Generating the appropriate FiniteElement's for fluid, structure,
and mesh motion
* Adding these elements to the hp::FECollection
* Setting the active fe index for the cells over the domain
* Counting the dofs per component with the
DoFTools::count_dofs_per_component method
When this is finished, it recognizes 0 dofs for the mesh motion (shown
below), as the active fe index for all those cells are associated with
the fluid.
If the motion FE is defined on the same part of the domain as the fluid,
then you need to define only two kinds of FESystem objects:
- one for the structure
- one for the fluid/motion components
You'd then have only to possible fe_index values, but you can still
subdivide the linear system into three blocks by splitting by vector
components if you want.
Best
W.
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
