On Fri, Aug 1, 2014 at 11:06 AM, David Knezevic <dkneze...@seas.harvard.edu> wrote: > > On 08/01/2014 12:50 PM, John Peterson wrote: >> >> On Fri, Aug 1, 2014 at 10:28 AM, John Peterson <jwpeter...@gmail.com> >> wrote: >>> >>> What I don't understand is why your original exodus file, which also >>> has num_dim=3, actually displays correctly in paraview! This might >>> take a big longer to track down... >> >> Here's something... your original file references side "5" of the TRI3 >> elements in question: >> elem_ss1 = 1, 17, 33 ; >> side_ss1 = 5, 5, 5 ; >> >> >> while the libmesh-generated version references side "1" >> elem_ss1 = 1, 17, 33 ; >> side_ss1 = 1, 1, 1 ; >> >> So while I don't understand what side "5" even means on a 3-sided >> element, clearly it displays properly in paraview! >> >> Otherwise, all the element numbers in the sets themselves match >> between the two files. >> > > Hi John, > > Thanks for looking into this. For me, the mesh I sent through displays fine > in CUBIT, and the sidesets are correct (see screenshot where I've > highlighted sideset 51 in orange). Then when I use it for a solve in libMesh > (e.g. a Poisson solve with zero Dirichlet BCs), I get BCs on internal edges. > I guess this is because libMesh sets side_ss1 to 1, as you pointed out.
Yeah, that make sense based on what I now know. > So it seems here that the problem is with the mesh that CUBIT is generating, > do you agree? I agree that "side_ss1 = 5, 5, 5" doesn't make any sense, so > I'm not sure what's going on there. Definitely agree. We could possibly put in some kind of workaround if you can figure out the "logic" behind these sideset IDs coming from Cubit. (If we mod it by n_sides or something, do we get the right ID?) -- John ------------------------------------------------------------------------------ Want fast and easy access to all the code in your enterprise? Index and search up to 200,000 lines of code with a free copy of Black Duck Code Sight - the same software that powers the world's largest code search on Ohloh, the Black Duck Open Hub! Try it now. http://p.sf.net/sfu/bds _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
