Hello Jason,

I am not sure, if FECollection is the correct class to use. Perhabs you want to have look at the FESystem class.

Best Regards,
Markus



Am 01.02.2012 18:14, schrieb Jason Sheldon:
Hi all,

I am working on an FSI solver using deal.ii. Currently I am attempting to set up a partitioned scheme using block(objects) and an hp::FECollection. I have 3 blocks in my system, one for the Fluid, one for the Structure, and one for the Mesh Motion. The equations governing the fluid and the mesh motion need to be defined on the same region of the domain. Currently I do not know of a way to do this this, and only one FESystem gets any dofs defined for it.

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.

(
DEAL:CC-00::Number of DOFS: 188269, (by partition: 188269)
DEAL:CC-00::Number of blocks: 3
DEAL:CC-00::Dofs by Block : 1836+186433+0
DEAL:CC-00::Block 0: 1,1,1,1,0,0,0,0,0:  1836
DEAL:CC-00::Block 1: 0,0,0,0,1,1,1,0,0:  186433
DEAL:CC-00::Block 2: 0,0,0,0,0,0,0,1,1:  0
)

Any assistance with this problem would be appreciated.  Thanks,

--
-Jason Sheldon


_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to