Good afternoon everyone, I decided to generate a transfinite mesh with gmsh and use the standard routine that you can find in this tutorial: Defining subdomains for different materials — FEniCSx tutorial to get tagged physical group. The weirdest thing is if you put a break point inside the last py…

Remove the physical marker on the center point of the facet. It creates an extra point in your mesh that shouldn’t be there (im not sure why gmsh creates an extra point though). Also added create connectivity: import pathlib import gmsh from dolfinx import io from meshio import Mesh as MeshIO fro…

Many thanks @dokken for your answer! Everything works fine now. I will try to use Open Cascade and divide my geometry in multiple part, so this point must be here ! Have a good evening!

Hello I exported a mesh from gmsh, but I have problems to imposse a boundary condicition Note, the file .msh it was created from gmsh software import gmsh import numpy as np from mpi4py import MPI from petsc4py import PETSc from dolfinx.fem import (Constant, Function, FunctionSpace, …

Without having the mesh or error message with stack trace, I see no way of giving further guidance.

Traceback (most recent call last): File "/home/gerard/Dropbox/FEniCSx/turbulent_model.py", line 44, in <module> bcu_inflow = dirichletbc(u_inlet, locate_dofs_topological(V, 1, facet_tags)) File "/usr/lib/petsc/lib/python3/dist-packages/dolfinx/fem/bcs.py", line 81, in locate_dofs_topological _entiti…

Thanks ! Yes, I’ve already used the gmsh.model.mesh.embed (and tag vertices) command for the non-transfinite mesh and it works fine. But with transfinite mesh I got the same error. Normally, a node exists as I set the good parity for the number of vertical cells, but the position is slightly differe…

If the type of facet_tags is mesh.MeshTagsMetaClass. You need to add this line before bcu_inflow features = facet_tags.find(int(your tag number for this DC condition)) bcu_inflow = dirichletbc(u_inlet, locate_dofs_topological(V, 1, features)) Have a look to the tutorial I’ve posted in the first co…

I’ve an other question. Is the order in topology.create_connectivity is important ? It’s to create the Agency List that we can find in topology.connectivity ?

[image] LGerard: Traceback (most recent call last): File "/home/gerard/Dropbox/FEniCSx/turbulent_model.py", line 44, in <module> bcu_inflow = dirichletbc(u_inlet, locate_dofs_topological(V, 1, facet_tags)) File "/usr/lib/petsc/lib/python3/dist-packages/dolfinx/fem/bcs.py", line 81, in locate_dof…

[image] flo1: I’ve an other question. Is the order in topology.create_connectivity is important ? It’s to create the Agency List that we can find in topology.connectivity ? You need to call topology.create_connectivity before calling certain functions in DOLFINx. Most functions will throw err…

Tranfinite gmsh mesh and read_meshtags lead to unexpected crash

Errors

dokken February 26, 2024, 8:54pm 4

You can embed the point: How to set 'boundary condition' inside the boundary? - #11 by dokken

Topic		Replies	Views
Need help converting GMSH to FEniCS installation	61	6664	September 28, 2022
Unable to find entity in map, 3D mesh gmsh Legacy DOLFIN mesh	5	70	April 2, 2025
Unable to import extruded geometry mesh mesh dolfinx	5	108	September 30, 2025
Fail to read tagged facets from mesh generated with GMSH mesh io , gmsh	4	112	March 1, 2025
Read meshtag from xdmf created by <mesh generation with gmsh> I/O mesh , dolfinx , gmsh	2	187	July 20, 2024

Tranfinite gmsh mesh and read_meshtags lead to unexpected crash

Related topics