gmsh issueshttps://gitlab.onelab.info/gmsh/gmsh/-/issues2024-03-19T02:12:58Zhttps://gitlab.onelab.info/gmsh/gmsh/-/issues/2846merge a STL surface with volumeric domain2024-03-19T02:12:58ZDuo Limerge a STL surface with volumeric domainHello GMSH,
I need some help with my model. I have a fault surface with curvature along both directions, which is in STL format. I first build a large box domain, and then merge this STL into GMSH, using gmsh.merge(stlFile). This STL sh...Hello GMSH,
I need some help with my model. I have a fault surface with curvature along both directions, which is in STL format. I first build a large box domain, and then merge this STL into GMSH, using gmsh.merge(stlFile). This STL should be an internal cut into the box. I am not sure how to do it. I used the following lines:
# Step 2: Load the STL surface file
gmsh.merge(stlFile)
gmsh.model.mesh.classifySurfaces(math.pi, curveAngle=math.pi)
# create a geometry for the discrete curves and surfaces
gmsh.model.mesh.createGeometry()
but I got error saying: Error : Wrong topology of boundary mesh for parametrization
for the createGeometry().
Any idea or suggestion would be appreciated:-)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2844Retrieve the tags after fuse command2024-03-16T12:53:04ZGuido VallifuocoRetrieve the tags after fuse commandGood morning
through the process described in the following script, I do not get any outDimTagsMap from the fuse command. is there a method by which I can retrieve the tags of the child surfaces?
```python
# we have to obtain a wather...Good morning
through the process described in the following script, I do not get any outDimTagsMap from the fuse command. is there a method by which I can retrieve the tags of the child surfaces?
```python
# we have to obtain a wathertight volume
cut_map1, cut_map2 = gmsh.model.occ.cut(wings_volume_dimtags, fuselage_volume_dimtags, -1, True,False)
gmsh.model.occ.synchronize()
children_map3, children_map4 = gmsh.model.occ.fuse(wings_volume_dimtags, fuselage_volume_dimtags, -1, True, True)
gmsh.model.occ.synchronize()
```
if I used surfaces instead of volumes could anything change?https://gitlab.onelab.info/gmsh/gmsh/-/issues/2841Periodic Mesh: Impossible to mesh periodic surface 22024-03-15T08:03:09ZRaffael CasagrandePeriodic Mesh: Impossible to mesh periodic surface 2I try to create a mesh with the following .geo file:
```
SetFactory("OpenCASCADE");
Torus(1) = {0, -0, 0, 0.5, 0.2, 2*Pi};
// Create slice of cake:
Point(10) = {0,0,-2.5,1};
Point(11) = {5,0,-2.5,1};
Point(12) = {5/Sqrt(2),5/Sqrt(2),-2....I try to create a mesh with the following .geo file:
```
SetFactory("OpenCASCADE");
Torus(1) = {0, -0, 0, 0.5, 0.2, 2*Pi};
// Create slice of cake:
Point(10) = {0,0,-2.5,1};
Point(11) = {5,0,-2.5,1};
Point(12) = {5/Sqrt(2),5/Sqrt(2),-2.5,1};
Line(3) = {11, 10};
Line(4) = {10, 12};
Circle(5) = {11, 10, 12};
Curve Loop(2) = {5, -4, -3};
Plane Surface(2) = {2};
Extrude {0, 0, 5} {
Surface{2};
}
BooleanIntersection{ Volume{2}; Delete;}{ Volume{1}; Delete; }
a = 1./Sqrt(2);
Periodic Surface {3} = {1} Affine {a,-a,0,0, a,a,0,0, 0,0,1,0, 0,0,0,1};
```
I tried the `Automatic`, `Delauny` and `Frontal-Delauny Algorithm`, but I always get the error "Impossible to mesh periodic surface 2" plus a few more errors.
Any ideas?https://gitlab.onelab.info/gmsh/gmsh/-/issues/2839independent meshes from the same file2024-03-14T16:13:28ZAndrea Vigliottiindependent meshes from the same fileI have a volume with _matrix_ material with some _inclusions_. Since they have very different typical length and properties I would like to have two separate meshes with very different mesh size and later deal with the interpolation betw...I have a volume with _matrix_ material with some _inclusions_. Since they have very different typical length and properties I would like to have two separate meshes with very different mesh size and later deal with the interpolation between the nodes that lie on the same surface on common boundaries.
One way to do that would be to have two separate script that produce the same model, one will mesh and produce an input file for the _matrix_ only, the other for the _inclusions_ only, and be careful about assigning labels to the surfaces in order to be able to match the nodes ion the two models that lie on common surfaces
I was wondering if there is a way to do this in one go, from the same model produce different meshes for different regions each with its own nodes on the same surface, belonging to the same set
this will also make it very easy to deal with numbering of the nodes and the elements,
thanks,
Andreahttps://gitlab.onelab.info/gmsh/gmsh/-/issues/2836How to mesh several parts and export them in one mesh file2024-03-14T15:28:21ZJohannes ACKVAHow to mesh several parts and export them in one mesh filehello,
i wonder how to mesh several parts and export them in ONE mesh file.
I would like to create them in several models, thus
gmsh.model.add("model_1")
..
gmsh.model.add("model_2")
..
but when i export with
gmsh.write('trial_mesh.med')...hello,
i wonder how to mesh several parts and export them in ONE mesh file.
I would like to create them in several models, thus
gmsh.model.add("model_1")
..
gmsh.model.add("model_2")
..
but when i export with
gmsh.write('trial_mesh.med')
only the mesh of the actual model is exported.
What must i do in addition?
thank you
Johannes_ACKVAhttps://gitlab.onelab.info/gmsh/gmsh/-/issues/2834Extracting connectivities of the reference element2024-03-14T13:51:32ZIgor KuszczakExtracting connectivities of the reference elementHello,
I am trying to recreate the `basix.cell.sub_entity_connectivity` functionality from `fenics-basix` in Gmsh. The function returns the connectivities of a reference element type. Is this possible to using Gmsh Python API? I imagine ...Hello,
I am trying to recreate the `basix.cell.sub_entity_connectivity` functionality from `fenics-basix` in Gmsh. The function returns the connectivities of a reference element type. Is this possible to using Gmsh Python API? I imagine that this could be hardcoded quite easily too but I don't understand the convention used for face numbering.https://gitlab.onelab.info/gmsh/gmsh/-/issues/2833Polar Extrusion2024-03-15T12:00:30ZAxel ThomasPolar ExtrusionHi,
I have a question about geometry extrusion. I have a cylindrical surface on which I would like to add particular patterns (representing sensors) and then like to form a 3D hexahedral mesh from this geometry. Is it possible to extrud...Hi,
I have a question about geometry extrusion. I have a cylindrical surface on which I would like to add particular patterns (representing sensors) and then like to form a 3D hexahedral mesh from this geometry. Is it possible to extrude 2D geometry in the (r, theta) plane?
Thanks in advance,
Axel
![pipe_with_embeded_sensor\_.png](/uploads/9a2e28fe4a653bd7d87589d398cd6b67/pipe_with_embeded_sensor_.png)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2831Intel-based Mac - MPICH - linking fails for libgmsh.dylib, missing _ompi_mpi_...2024-03-11T20:21:34ZTidbit SoftwareIntel-based Mac - MPICH - linking fails for libgmsh.dylib, missing _ompi_mpi_* symbolsI've successfully compiled/configured Gmsh for a number of platforms. I updated today from version 4.10.5 to 4.12.2. I am compiling on an Intel-based Mac running the latest version of macOS, using MPICH as provided by PETSc. Everything c...I've successfully compiled/configured Gmsh for a number of platforms. I updated today from version 4.10.5 to 4.12.2. I am compiling on an Intel-based Mac running the latest version of macOS, using MPICH as provided by PETSc. Everything compiles, but at linktime I get,
```
[ 94%] Linking CXX shared library libgmsh.dylib
ld: warning: -commons use_dylibs is no longer supported, using error treatment instead
Undefined symbols for architecture x86_64:
"_ompi_mpi_char", referenced from:
GmshRemote() in GmshRemote.cpp.o
GmshRemote() in GmshRemote.cpp.o
GmshRemote() in GmshRemote.cpp.o
GmshRemote() in GmshRemote.cpp.o
GmshRemote() in GmshRemote.cpp.o
GmshRemote() in GmshRemote.cpp.o
GmshRemote() in GmshRemote.cpp.o
GmshRemote() in GmshRemote.cpp.o
...
```
and so on for many `_opmi` symbols. At configure time, MPICH libs and binaries were found, OpenMPI was not (which is expected, and good). I assume there's some section in the configuration that's broken, causing a conditional that makes GmshRemote.cpp use the OpenMPI implementation in every case or by mistake (given that MPICH was detected). I haven't started digging yet, but thought it was worth a bug report anyhow.https://gitlab.onelab.info/gmsh/gmsh/-/issues/2830Facets are overlapping(python)2024-03-06T11:26:35ZJOHN PSYCHFacets are overlapping(python)I tried to 3D mesh the goemtry but the following error occurs:
"Invalid boundary mesh (overlapping facets) on surface 8 surface 8"
I tried to use "coherence" as explained in other issues but it didn't work.
[Kelvin_cell_one_body_v2.txt]...I tried to 3D mesh the goemtry but the following error occurs:
"Invalid boundary mesh (overlapping facets) on surface 8 surface 8"
I tried to use "coherence" as explained in other issues but it didn't work.
[Kelvin_cell_one_body_v2.txt](/uploads/9e2082cd9b422ca6ce50d52473f9fd59/Kelvin_cell_one_body_v2.txt)
```
import gmsh
import sys
gmsh.initialize()
gmsh.open("c:\GMSH\Kelvin_cell_one_body v[Kelvin_cell_one_body_v2.txt]
gmsh.model.geo.synchronize()
# Generate mesh:
gmsh.model.mesh.generate(3)
gmsh.model.coherence
# Write mesh data:
gmsh.write("GFG.msh")
gmsh.model.geo.synchronize()
gmsh.model.coherence()
gmsh.fltk.run()
gmsh.finalize()
```https://gitlab.onelab.info/gmsh/gmsh/-/issues/2829FUSE with no overlap (CPP API)2024-03-04T08:28:33ZLars NielsenFUSE with no overlap (CPP API)# GMSH FUSE
Dear GMSH board
I have a task which I would like to solve using gmsh and the cpp API. I have two volumes and I would like to union them and use the output from the stl-file. For this I use the FUSE command. When there is an...# GMSH FUSE
Dear GMSH board
I have a task which I would like to solve using gmsh and the cpp API. I have two volumes and I would like to union them and use the output from the stl-file. For this I use the FUSE command. When there is an overlap between the two volumes the cpp routine works fine. When the two volumes share a surface, my stl-file is empty. Can you help me making the union when the two volumes share a surface?
I have made two boxes and a FUSE in the GUI and that works fine. What am I doing wrong when using the API?
I have uploaded an example.
[Test.cpp](/uploads/2d433a1868b94dae44e74a87c2824781/Test.cpp)
Kind regards N_Bohrhttps://gitlab.onelab.info/gmsh/gmsh/-/issues/28283D Boundary Layer intersects itself2024-03-12T14:47:48ZOle Albers3D Boundary Layer intersects itselfHi, I'm trying to create a 3D Boundary Layer for an aircraft in a farfield.
But when the Boundary Layer is to thick, I get an Error meshing them im 3D (PLC Error: A segment and a facet intersect at point). In the Fltk it looks like this...Hi, I'm trying to create a 3D Boundary Layer for an aircraft in a farfield.
But when the Boundary Layer is to thick, I get an Error meshing them im 3D (PLC Error: A segment and a facet intersect at point). In the Fltk it looks like this happens where the wing and fuselage intersect and the Boundary Layer extrusion intersects itself.
I create the Boundary Layer like in the naca_boundary_layer_3d.py (#2812 (different aircraft this time))
In my case, when I choose 1e-6 as the thickness of the first layer I'll get the Error, when I chose 1e-7 I won't.
I attached both meshes in 2D. Is there anything I can do to prevent this from happening when meshing them 3D?
[BoundaryLayer_1e-6.msh](/uploads/1a2cc2785ce399c3cc8edff32cdf0a97/BoundaryLayer_1e-6.msh)
[BoundaryLayer_1e-7.msh](/uploads/c7fa477f8e723bb1c2901708f1cb3db9/BoundaryLayer_1e-7.msh)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2826When using Compound Curve, mesh size failure is controlled via Transfinite Cu...2024-03-02T12:25:19Zmike sunWhen using Compound Curve, mesh size failure is controlled via Transfinite Curve.I made a simple example, by Compound Surface {1, 2} to see the two surfaces as a composite surface, and also used Compound Curve, in order to better control the size of the mesh, I would like to control the mesh size according to the len...I made a simple example, by Compound Surface {1, 2} to see the two surfaces as a composite surface, and also used Compound Curve, in order to better control the size of the mesh, I would like to control the mesh size according to the length of the curve, I only found Transfinite Curve, according to the length of the setCompound Curve to set the number of control points to control the mesh size, but when the curve is setCompound, the Transfinite command fails, the Transfinite command for the curve that is not setCompound succeeds. I would like to know how should I assign mesh points for this Compound Curve? The mesh generation effect image and geo file are as follows:![1709279213287](/uploads/ea343535354e35a3f1293dee64728c8e/1709279213287.jpg)[compound1.geo](/uploads/b5e15e47fd9b039d1383613dac8ca5db/compound1.geo)
The contents of the geo file are as follows
// Gmsh project created on Fri Mar 01 15:24:45 2024
SetFactory("OpenCASCADE");
//+
Point(1) = {0, 0, 0, 1.0};
//+
Point(2) = {1, 0, 0, 1.0};
//+
Point(3) = {1, 1, 0, 1.0};
//+
Point(4) = {0, 1, 0, 1.0};
//+
Line(1) = {1, 2};
//+
Line(2) = {2, 3};
//+
Line(3) = {3, 4};
//+
Line(4) = {4, 1};
//+
Curve Loop(1) = {4, 1, 2, 3};
//+
Plane Surface(1) = {1};
//+
Point(5) = {1, 1.5, 0, 1.0};
//+
Point(6) = {0, 1.5, 0, 1.0};
//+
Line(5) = {3, 5};
//+
Line(6) = {5, 6};
//+
Line(7) = {6, 4};
//+
Curve Loop(2) = {6, 7, -3, 5};
//+
Plane Surface(2) = {2};
//+
Compound Curve {4, 7};
//+
Compound Curve {2, 5};
//+
Compound Surface {1, 2};
//+
Transfinite Curve {7, 5} = 5 Using Progression 1;
//+
Transfinite Curve {6, 3, 4, 2, 1} = 10 Using Progression 1;https://gitlab.onelab.info/gmsh/gmsh/-/issues/2825mesh transition question with python API2024-03-03T11:52:27ZJohannes ACKVAmesh transition question with python APIhello
see the attachment for my question.
I want to have the 4 bottom volumic parts meshed by extrusion with hexaeders. But the extrude method:
gmsh.model.geo.extrude( [(2,11),(2,12),(2,13),(2,14),], length,0,0, numElements=[1,1,1], ...hello
see the attachment for my question.
I want to have the 4 bottom volumic parts meshed by extrusion with hexaeders. But the extrude method:
gmsh.model.geo.extrude( [(2,11),(2,12),(2,13),(2,14),], length,0,0, numElements=[1,1,1], heights=[ 0.2, 0.5, 1.], recombine=True)
does not work here (meshing fails with "Exception: Could not find extruded node (1.6, -1.732050807568877, -1) in surface 87" ) and would not create coincident nodes along the lateral surfaces.
it seems that it fails because gmsh "wants" to avoid non-coincident nodes on the lateral surfaces. I would have to calculate by hand the values in heights=[] ( the above [ 0.2, 0.5, 1.] is wrong because only estimates). Calculating these values seems impossible because the mesh biasing of the hexahedral part has been on other lines set with something like
setTransfiniteCurve( l1, nb_of_nodes, meshType="Progression", coef=2)
and this methods has a "built-in" rounding which is not transparent to me.
My supposition is that for the extrusion instead of
numElements=[1,1,1], heights=[ 0.2, 0.5, 1.]
there should be a method like the above "nb_of_nodes, meshType="Progression", coef=2"
Does this exist? i didn't find it in the doc (gmsh.py). Or how can this task be resolved?
Many thanks
Johannes_ACKV[mesh_transition_problem.zip](/uploads/45beaa779f5ae28e2f1adbc246c6d092/mesh_transition_problem.zip)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2820Creating inclusion from a STEP file2024-02-28T19:36:53ZAntoine BEAUQUELCreating inclusion from a STEP fileHello all,
I'm contacting you with the following problem: I'm working on the meshing of a geometry defined using a .stp file (example attached).
When meshing the part, I'd like to recover the volumes of the recesses so that I can mesh ...Hello all,
I'm contacting you with the following problem: I'm working on the meshing of a geometry defined using a .stp file (example attached).
When meshing the part, I'd like to recover the volumes of the recesses so that I can mesh them and obtain a conform mesh with the rest of the part.
I assume that this operation is possible, but I haven't found the solution to this problem in the documentation or in other issues.
Thank you for your help and all the work.
[TOLE MINCE.stp](/uploads/36c89242f4cfec9154301150a8d7afa7/TOLE_MINCE.stp)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2819Huge memory consumption when gmsh.initialize() (Julia API) on windows2024-02-28T23:55:23ZTommy AHuge memory consumption when gmsh.initialize() (Julia API) on windowsHi, all. I just tried Julia API and noticed `gmsh.initialize()` may have some problem. When I run this, I can see a huge memory consumption (going up and down repetitively from 0 to 4GB) on the task manager. The initialization is done af...Hi, all. I just tried Julia API and noticed `gmsh.initialize()` may have some problem. When I run this, I can see a huge memory consumption (going up and down repetitively from 0 to 4GB) on the task manager. The initialization is done after about 30 seconds (with no error). I can use gmsh functions with no problem after that but the problem makes it practically unusable. I confirmed this on two of my windows machines (10 & 11). I also tried it on Ubuntu for windows (11) but I didn't see the same problem.
How to reproduce:
```sh
>> winget install julia -s msstore
```
then launch the Julia pkg manager
```sh
>> add Gmsh
```
then run the following code.
```jl
import Gmsh: gmsh
gmsh.initialize()
```
Julia version: 1.10.1
Gmsh version: 0.3.0
<details>
<summary>systeminfo</summary><pre><code class="language-bash">
Host Name:
OS Name: Microsoft Windows 11 Home
OS Version: 10.0.22621 N/A Build 22621
OS Manufacturer: Microsoft Corporation
OS Configuration: Standalone Workstation
OS Build Type: Multiprocessor Free
Registered Owner:
Registered Organization: N/A
Product ID:
System Manufacturer: Dell Inc.
System Model: Inspiron 14 5425
System Type: x64-based PC
Processor(s): 1 Processor(s) Installed.
[01]: AMD64 Family 25 Model 80 Stepping 0 AuthenticAMD ~2301 Mhz
BIOS Version: Dell Inc. 1.12.1, 2023/10/12
Windows Directory: C:\WINDOWS
System Directory: C:\WINDOWS\system32
Boot Device: \Device\HarddiskVolume1
System Locale: ja;Japanese
Input Locale: ja;Japanese
Time Zone: (UTC+09:00) Osaka, Sapporo, Tokyo
Total Physical Memory: 7,528 MB
Available Physical Memory: 2,419 MB
Virtual Memory: Max Size: 18,126 MB
Virtual Memory: Available: 3,255 MB
Virtual Memory: In Use: 14,871 MB
Page File Location(s): C:\pagefile.sys
Domain: WORKGROUP
Logon Server:
Hotfix(s): 4 Hotfix(s) Installed.
[01]: KB5034467
[02]: KB5012170
[03]: KB5034765
[04]: KB5034225
Network Card(s): 1 NIC(s) Installed.
[01]: MediaTek Wi-Fi 6 MT7921 Wireless LAN Card
Connection Name: Wi-Fi
DHCP Enabled: Yes
DHCP Server: 192.168.10.1
Hyper-V Requirements: A hypervisor has been detected. Features required for Hyper-V will not be displayed.
</code></pre>
</details>https://gitlab.onelab.info/gmsh/gmsh/-/issues/2815Coherence / BooleanFragments different outputs across Gmsh versions2024-02-25T16:24:55ZJoão PinheiroCoherence / BooleanFragments different outputs across Gmsh versionsHello,
I've been having issues with the newest version of Gmsh 4.12.2. When i input my .geo files i get weird geometries, which does not happen on the previous version 4.11.1. I can't figure out what I am doing wrong.
Attached are 2 simp...Hello,
I've been having issues with the newest version of Gmsh 4.12.2. When i input my .geo files i get weird geometries, which does not happen on the previous version 4.11.1. I can't figure out what I am doing wrong.
Attached are 2 simple files that hopefully show the problem.
On the first file error_coherence, Volume in BoundingBox does not find any volumes.
On error_fragments, "Error : Could not fix wire in surface" when applying the second BooleanFragments.
[error_coherence.geo](/uploads/2c76b52558266f94e6bef9f8f43a6f10/error_coherence.geo)
[error_fragments.geo](/uploads/8d127e9eb400fb0388b6800e1895b132/error_fragments.geo)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2814Transfinite Curve Progression for two Lines2024-02-25T20:23:52ZMads SkjaerbaekTransfinite Curve Progression for two LinesHi
I’m trying to simulate a nozzle flow into a box for an axisymmetric case. I want to create a structured mesh with higher density near the nozzle (Progression). To get the correct boundary conditions, I must define two lines near the...Hi
I’m trying to simulate a nozzle flow into a box for an axisymmetric case. I want to create a structured mesh with higher density near the nozzle (Progression). To get the correct boundary conditions, I must define two lines near the nozzle. Is there a way to make GMSH behave like the two lines are one when using the transfinite progression command? I have tried the Compound line but can't get it to work. I have attached a picture of my domain. And this is my code:
```
// Gmsh project created on Tue Feb 20 10:09:58 2024
SetFactory("OpenCASCADE");
//+
Point(1) = {0, 400, 0, 1.0};
//+
Point(2) = {1400, 400, 0, 1.0};
//+
Point(3) = {1400, 0, 0, 1.0};
//+
Point(4) = {1400, 6.65, 0, 1.0};
//+
Point(5) = {1425.3, 12.65, 0, 1.0};
//+
Point(6) = {1425.3+22+45+34.95, 12.65, 0, 1.0};
//+
Point(7) = {1425.3+22+45+34.95, 0, 0, 1.0};
//+
Point(8) = {1400+2.53*1, 6.70136001, 0, 1.0};
//+
Point(9) = {1400+2.53*2, 6.99751999, 0, 1.0};
//+
Point(10) = {1400+2.53*3, 7.62848, 0, 1.0};
//+
Point(11) = {1400+2.53*4, 8.554639999, 0, 1.0};
//+
Point(12) = {1400+2.53*5, 9.65, 0, 1.0};
//+
Point(13) = {1400+2.53*6, 10.74536, 0, 1.0};
//+
Point(14) = {1400+2.53*7, 11.67152, 0, 1.0};
//+
Point(15) = {1400+2.53*8, 12.30248, 0, 1.0};
//+
Point(16) = {1400+2.53*9, 12.59864, 0, 1.0};
//+
Point(17) = {0, 0, 0, 1.0};
//+
Point(18) = {0, 40, 0, 1.0};
//+
Point(19) = {1400+25.3, 0, 0, 1.0};
//+
Point(20) = {0, 6.65, 0, 1.0};
//+
Point(21) = {1400, 40, 0, 1.0};
//+
Line(1) = {5, 6};
//+
Line(2) = {6, 7};
//+
Line(3) = {3, 19};
//+
Line(4) = {4, 3};
//+
Line(5) = {1, 2};
//+
Line(6) = {2, 21};
//+
Spline(7) = {4, 8, 9, 10, 11, 12, 13, 14, 15, 16, 5};
//+
Line(8) = {1, 18};
//+
Line(9) = {19, 5};
//+
Line(10) = {7, 19};
//+
Line(11) = {17, 3};
//+
Line(12) = {18, 20};
//+
Line(13) = {21, 4};
//+
Line(14) = {20, 17};
//+
Curve Loop(1) = {5, 6, 13, 4, -11, -14, -12, -8};
//+
Plane Surface(1) = {1};
//+
Curve Loop(2) = {3, 9, -7, 4};
//+
Plane Surface(2) = {2};
//+
Curve Loop(3) = {1, 2, 10, 9};
//+
Plane Surface(3) = {3};
//+
Transfinite Surface {1} = {17, 3, 2, 1};
//+
Transfinite Surface {3} = {19, 7, 6, 5};
//+
Transfinite Curve {8, 12} = 20 Using Progression 1;
//+
Transfinite Curve {6, 13} = 20 Using Progression 0.6;
//+
Recombine Surface {1};
//+
Recombine Surface {3};
//+
Transfinite Curve {4, 14} = 10 Using Progression 1;
//+
Transfinite Curve {11, 5} = 700 Using Progression 0.999;
```
It is line 8 and 12, and 6 and 13, that I want GMSH to consider as one line when using Transfinite Curve
![Screenshot_2024-02-21_153058](/uploads/123deb9b3ada62d165c393b592409c79/Screenshot_2024-02-21_153058.png)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2813about hexahedra mesh generate2024-02-25T21:02:04Zmingyu huangabout hexahedra mesh generateI am currently engaged in a project that necessitates the generation of either all-tetrahedral or all-hexahedral meshes for various geometrical configurations related to GIS. I've observed that Gmsh effortlessly generates tetrahedral ele...I am currently engaged in a project that necessitates the generation of either all-tetrahedral or all-hexahedral meshes for various geometrical configurations related to GIS. I've observed that Gmsh effortlessly generates tetrahedral elements. However, when it comes to generating hexahedral meshes, it seems my only option is to use the subdivisionAlgorithm to convert tetrahedral meshes into hexahedral ones. I also explored using the transfinite algorithm and recombine3DAll, but the recombine option doesn't seem to work, and the transfinite algorithm might not meet my requirements.
Is there an alternative method within Gmsh for generating hexahedral meshes directly? I am particularly interested in employing the frontal-Delaunay for quads or quasi-structured quad as the 2D algorithm, aiming to produce hexahedral meshes resembling a rectangular prism. Thank you for your assistance
[t1.geo](/uploads/b2b71be2a7afe7724628f92d422a4ee8/t1.geo)![p1](/uploads/e72af9b63d7d49b04fc99f78c9f1dd6c/p1.png)https://gitlab.onelab.info/gmsh/gmsh/-/issues/2811Define material properties in mesh (CellData)2024-02-25T21:03:50ZOctavio Castillo ReyesDefine material properties in mesh (CellData)Dear all,
I'm working with gmsh to build a mesh of the earth subsurface. I would to assign some properties to each element of a given physical volume. It is possible? If yes, How can I do it (from a .geo file)?
Thank you for your time ...Dear all,
I'm working with gmsh to build a mesh of the earth subsurface. I would to assign some properties to each element of a given physical volume. It is possible? If yes, How can I do it (from a .geo file)?
Thank you for your time and support
Besthttps://gitlab.onelab.info/gmsh/gmsh/-/issues/2810Gmsh Orthogonality check2024-02-25T21:05:03ZIlaria FicheraGmsh Orthogonality checkHello,
I am using Gmsh to create a tetrahedral volumetric mesh in a room. I need this mesh to be orthogonal. I was searching if Gmsh has already a function that could check if the created mesh is orthogonal. Is there something like that...Hello,
I am using Gmsh to create a tetrahedral volumetric mesh in a room. I need this mesh to be orthogonal. I was searching if Gmsh has already a function that could check if the created mesh is orthogonal. Is there something like that? Or does Gmsh already create an orthogonal mesh if certain settings are considered?
Thank you
Regards,
Ilaria