[petsc-users] How to construct DMPlex of cells with different topological dimension?

袁煕 yuanxi at advancesoft.jp
Sun Oct 31 20:01:58 CDT 2021


| Yuan: can you add a fake Hex over the hanging shell here? Give them zero
stiffness and constrain the hanging nodes that the fake elements create.

I don't think it is a good idea although it would work theoretically.  It
would be a disaster when you try to do calculations over several tens of
millions of elements this way. Maybe I should give up automatic
interpolation but call DMPlexDistribute directly after setting up Dofs upon
each node. Is it possible?

Yuan

2021年11月1日(月) 5:01 Mark Adams <mfadams at lbl.gov>:

>
>
>>>    Is there anything wrong in my above consideration?  Any suggestions?
>>>
>>>         ------------
>>>       /|             /|
>>>      / |           /  |  cell 0: Hex
>>>     /  |          /   |
>>>    ------------/   |
>>>    |   |         |   |
>>>    |   |         |   |   cell 1: Quad
>>>    |   --------|---|------------
>>>    |  /          |  /             /
>>>    | /           | /             /
>>>    |/            |/             /
>>>    -------------------------------------------
>>>                                   cell 2: line
>>>
>>> Much thanks for your help.
>>>
>>
>> If you are solving something where everything is embedded in a volumetric
>> mesh, then there is no problem. However, if you really have
>> the mesh above, where lower dimensional pieces are sticking out of the
>> mesh, then Plex can represent the mesh, but automatic interpolation
>> (creation of edges and faces) will not work.
>>
>
> Yuan: can you add a fake Hex over the hanging shell here? Give them zero
> stiffness and constrain the hanging nodes that the fake elements create.
> Mark
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20211101/fa70c663/attachment.html>


More information about the petsc-users mailing list