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
