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