Re: [deal.II] Query about external meshing software

2021-03-09 Thread Zelu Xu
Hello Krishanu: The error is made when you define the line/surface. You don’t have to check the mesh file or each mesh element but the .geo file. For example, here my line loop defines the surface bounded by vertex in the order of [3.7,8,4], but I defined my line as the following. So when I

Re: [deal.II] Query about external meshing software

2021-03-09 Thread krishan...@gmail.com
Hello Jane, Thanks for the response. I understand why there could be an issue with the ordering. But I am using the GUI of Gmsh and the mesh created has about 400 elements. So it is a little difficult to manually check the ordering of the vertices of each element. Do you know if there is any

Re: [deal.II] Query about external meshing software

2021-03-09 Thread Zelu Xu
Hello Krishna: I had this issue before. It just means that when you select mesh entities (vertex) to form a surface, you have to pay attention to the order of which entity goes first, second,... it would be clear if you look at the geometry txt file. Best regards Jane On Tue, Mar 9, 2021 at 3:21

Re: [deal.II] Query about external meshing software

2021-03-09 Thread krishan...@gmail.com
Hello Paras, Thanks for your response. I have been using the GUI of Gmsh to generate the mesh. Do you have any suggestions about what to do in Gmsh to avoid the negative volume issue? Do meshes generated with Abaqus, are able to avoid the above-mentioned issue of negative volume? Thanks in

Re: [deal.II] Query about external meshing software

2021-03-09 Thread Paras Kumar
I use the .inp format of Abaqus, since my code relies on material-IDs which (in my limited experience) are easy to handle with .inp format. I generate meshes using both GMSH (here as you mentioned, one needs to be careful to add curves/surface parts in the correct order to avoid negative volumes)

[deal.II] Query about external meshing software

2021-03-09 Thread krishan...@gmail.com
I am wondering what meshing software produces output files that would work best with deal.ii. I have tried Gmsh, but it seems when this mesh (generated by Gmsh) is imported in deal.ii, it is facing the issue of negative volume fraction of a cell (could be because of the order of nodes written

Re: [deal.II] boundary matrix

2021-03-09 Thread Wolfgang Bangerth
On 3/7/21 10:16 PM, Nikki Holtzer wrote: Thank you for your input! That was my first attempt at tackling this. However, I run into problems with the line for(const auto & face : cell-> face_iterators()). No matter what member function I choose off of the DoFCellAccessor class template