Jean-Paul, The boundary ID output from the four different calculations are as follows (before any boundary IDs are set to 6):
1) 3D, no imported mesh: equal numbers for boundary IDs 0, 1, 2, 3, 4, and 5. This calculation converges to a solution (as expected). 2) 2D, no imported mesh: equal numbers for boundary IDs 0, 1, 2, and 3. This calculation converges to a solution (as expected). 3) 3D, imported mesh: some boundary IDs are 1 and 2, while most boundary IDs are 3. This calculation converges to a solution. 4) 2D, imported mesh: all boundary IDs are 0. This calculation does not converge to a solution. After cycling through all of the face cells and using set_boundary_id(), each calculation (correctly) has 6 as a boundary ID for some number of faces. With the information above and your comments, I take it that boundary ID = 0 means the boundary is traction-free and no dirichlet BC is specified? Assuming this is the case, then it makes sense there would be no convergence. It also means that the 3D calculation with the imported mesh is not the same as the original. I have not yet determined why I cannot read the .vtu output files in Paraview, which is why I did not know this before (I could only tell it converged to a solution). I do find it interesting that the two meshes were not imported with the same boundary IDs. I guess in the future I should assume that all boundary IDs are zero and reset them for each calculation? Also, is there a quick reference guide for what each boundary ID means? I have not yet found such a reference on the deal.ii website. Thank you for your help, James On Monday, April 13, 2020 at 3:26:14 AM UTC-4, Jean-Paul Pelteret wrote: > > It might be that one of the constrained boundaries is not traction free... > > > Correction: This should have read "It might be that one of the constrained > boundaries is *now* traction free..." > > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/f9d48ea6-76ab-4cd9-a5d7-aaedc0d81537%40googlegroups.com.
