[petsc-users] DMPlex and Gmsh

[Matthew Knepley]

On Wed, Nov 8, 2023 at 4:50 PM Blaise Bourdin <bourdin at mcmaster.ca> wrote:

> Hi,
>
> I think that you need to use the magical keyword “
> -dm_plex_gmsh_mark_vertices” for that
>

I try to describe the options here:

  https://petsc.org/main/manualpages/DMPlex/DMPlexCreateGmsh/

   Thanks,

      Matt


> Blaise
>
> On Nov 8, 2023, at 1:13 PM, Sharan Roongta <s.roongta at mpie.de> wrote:
>
> Caution: External email.
>
> Dear Petsc team,
>
> I want to load a .msh file generated using Gmsh software into the DMPlex
> object. There are several things I would want to clarify, but I would like
> to start with “Physical tags”.
>
> If I have defined “Physical Points”, “Physical Surface”, and “Physical
> Volume” in my .geo file, I get the physical tags in the “.msh” file.
> When I load this mesh in DMPlex, and view the DM:
>
> call DMView(globalMesh, PETSC_VIEWER_STDOUT_WORLD,err_PETSc)
>   CHKERRQ(err_PETSc)
>
> This is the output I get:
>
> DM Object: n/a 1 MPI process
>   type: plex
> n/a in 3 dimensions:
>   Number of 0-cells per rank: 14
>   Number of 1-cells per rank: 49
>   Number of 2-cells per rank: 60
>   Number of 3-cells per rank: 24
> Labels:
>   celltype: 4 strata with value/size (0 (14), 6 (24), 3 (60), 1 (49))
>   depth: 4 strata with value/size (0 (14), 1 (49), 2 (60), 3 (24))
>   Cell Sets: 1 strata with value/size (8 (24))
>   Face Sets: 6 strata with value/size (2 (4), 3 (4), 4 (4), 5 (4), 6 (4),
> 7 (4))
>
> I was expecting to get the “Node Sets” or “Vertex Sets” also. Is my
> assumption wrong?
>
> If yes, then how can one figure out the boundary nodes and their tags
> where I want to apply certain boundary conditions?
> Currently we apply boundary conditions on faces, therefore “Face Sets” was
> enough. But now we want to apply displacements on certain boundary nodes.
>
> I have also attached the .geo and .msh file (hope you can open it)
> The Petsc version I am using is 3.18.6.
>
>
> Thanks and Regards,
> Sharan Roongta
>
> -----------------------------------------------
> Max-Planck-Institut für Eisenforschung GmbH
> Max-Planck-Straße 1
> D-40237 Düsseldorf
>
> Handelsregister B 2533
> Amtsgericht Düsseldorf
>
> Geschäftsführung
> Prof. Dr. Gerhard Dehm
> Prof. Dr. Jörg Neugebauer
> Prof. Dr. Dierk Raabe
> Dr. Kai de Weldige
>
> Ust.-Id.-Nr.: DE 11 93 58 514
> Steuernummer: 105 5891 1000
> -------------------------------------------------
> Please consider that invitations and e-mails of our institute are
> only valid if they end with …@mpie.de.
> If you are not sure of the validity please contact rco at mpie.de
>
> Bitte beachten Sie, dass Einladungen zu Veranstaltungen und E-Mails
> aus unserem Haus nur mit der Endung …@mpie.de gültig sind.
> In Zweifelsfällen wenden Sie sich bitte an rco at mpie.de
>
>
>
> ------------------------------
> -------------------------------------------------
> Stay up to date and follow us on LinkedIn, Twitter and YouTube.
>
> Max-Planck-Institut für Eisenforschung GmbH
> Max-Planck-Straße 1
> D-40237 Düsseldorf
>
> Handelsregister B 2533
> Amtsgericht Düsseldorf
>
> Geschäftsführung
> Prof. Dr. Gerhard Dehm
> Prof. Dr. Jörg Neugebauer
> Prof. Dr. Dierk Raabe
> Dr. Kai de Weldige
>
> Ust.-Id.-Nr.: DE 11 93 58 514
> Steuernummer: 105 5891 1000
>
>
> Please consider that invitations and e-mails of our institute are
> only valid if they end with …@mpie.de.
> If you are not sure of the validity please contact rco at mpie.de
>
> Bitte beachten Sie, dass Einladungen zu Veranstaltungen und E-Mails
> aus unserem Haus nur mit der Endung …@mpie.de gültig sind.
> In Zweifelsfällen wenden Sie sich bitte an rco at mpie.de
> -------------------------------------------------
> <cube_check.msh><cube_check.geo>
>
>
> —
> Canada Research Chair in Mathematical and Computational Aspects of Solid
> Mechanics (Tier 1)
> Professor, Department of Mathematics & Statistics
> Hamilton Hall room 409A, McMaster University
> 1280 Main Street West, Hamilton, Ontario L8S 4K1, Canada
> https://www.math.mcmaster.ca/bourdin | +1 (905) 525 9140 ext. 27243
>
>

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20231108/3074f974/attachment-0001.html>