Hi, Daniel,
I see. Thank you very much. So, I think our method may be available as long
as it can separate the face.
Best,
Xiang
On Tuesday, 13 October 2020 at 14:38:29 UTC-7 d.arnd...@gmail.com wrote:
> Xiang,
>
> deal.II doesn't have functionality to create separate cells sharing a face
Hi, Bruno,
Thanks a lot. We will have a try based on the example you gave us.
Best,
Xiang
On Tuesday, 13 October 2020 at 10:49:08 UTC-7 bruno.t...@gmail.com wrote:
> Xiang,
>
> Here is an example on how to do it:
>
>
Xiang,
deal.II doesn't have functionality to create separate cells sharing a face
but only to join them. Hence, whatever input you provide to
Triangulation::create_triangulation() must have these faces already. GridIn
should work for that if the faces are disjoint in the input file.
Best,
Daniel
Xiang,
Here is an example on how to do it:
https://github.com/geodynamics/aspect/blob/master/cookbooks/prescribed_velocity/prescribed_velocity.cc#L297-L303
As you can see, you need to know which dof you need to constrain
Best,
Bruno
Le mar. 13 oct. 2020 à 11:02, 孙翔 a écrit :
>
>
> Hi, Bruno,
Hi, Daniel,
Thank you very much. So, if I want to disjoint the faces, how should I do
using dealII or I need to disjoint them manually? If I manually disjoint
them, should I duplicate the joint nodes and assign them with the new
number, then assemble the adjacent element with the duplicated
Hi, Daniel,
Thank you very much. So, if I want to disjoint the faces, how should I do
using dealII or I need to disjoint them manually? If I manually disjoint
them, should I duplicate the joint nodes and assign them with the new
number, then assemble the adjacent element with the duplicated
Xiang,
Adding to Bruno's response: If the interface you want to apply Dirichlet
boundary conditions to stays the same through the whole simulation, then
treating the corresponding faces as disjoint while creating your mesh might
be a better option.
Setting constraints yourself requires you to