Wolfgang, I follow your suggestion, and may go that route, but I was wondering if it would be possible to define two active indices on a cell. Ideally, I would prefer the fluid and mesh motion to have separate FESystems.
On Wed, Feb 1, 2012 at 3:18 PM, Wolfgang Bangerth <[email protected]>wrote: > > 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/<http://www.math.tamu.edu/%7Ebangerth/> > > -- -Jason Sheldon
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
