Thanks Roy and Alexander,

I appreciate the insight, I will ask those maintaining PyGmsh about the
order of elements if I can't seem to find the option readily myself, and
then report back here.

Disregarding order for now, if I simply run this as it is, is it assuming
Tri3 elements and a number of elements corresponding to those given in the
mesh?

On Wed, Nov 22, 2017 at 4:36 PM, Alexander Lindsay <alexlindsay...@gmail.com
> wrote:

> I've never used pygmsh, but I know if you're running gmsh from the command
> line, passing the option `-order 2` will create second order elements.
>
> On Wed, Nov 22, 2017 at 2:14 PM, Roy Stogner <royst...@ices.utexas.edu>
> wrote:
>
>>
>> On Wed, 22 Nov 2017, Zack Vitoh wrote:
>>
>>  GmshIO gmsh_io(mesh);
>>>  gmsh_io.read("disk.vtu");
>>>
>>> but--and I am not sure if this is possible--imposing additional nodes on
>>> each element generated in the triangularization generated by gmsh with
>>> nodes such that I have, say, Tri6 elements.
>>>
>>
>> You have Tri3 elements now?  After the read calling
>> mesh.all_second_order();
>> will turn them into Tri6 elements.
>>
>> Is this possible through libmesh, or something I would I have to do
>>> something in gmsh instead?
>>>
>>
>> I'm afraid if you have a disk then this is something you *should* do
>> in gmsh instead if that's possible.  By the time libMesh reads a
>> first-order triangulation, curvature information has been discarded,
>> and we don't try to reconstruct it, so you'll still get straight-sided
>> elements rather than curved-sided elements.  The only way to get the
>> best second-order boundary geometry is to patch it up manually,
>> looping over boundary nodes and "snapping" them to the curve where you
>> know they should be.
>> ---
>> Roy
>>
>>
>> ------------------------------------------------------------
>> ------------------
>> Check out the vibrant tech community on one of the world's most
>> engaging tech sites, Slashdot.org! http://sdm.link/slashdot
>> _______________________________________________
>> Libmesh-users mailing list
>> Libmesh-users@lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/libmesh-users
>>
>
>
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to