[ashish.svm]

FSI (fluid structure interaction) problem. Need to generate a mesh around the given fluid mesh. Fluid region mesh will be 2D, solid region 3D

An example of fluid region mesh is shown below
top view (XY)
fluidMesh_XY_view.png
3D view
fluidMesh_3D_view.png

Fluid region mesh is generated from the output of another code and this mesh is in OpenFOAM format. To use this mesh in another meshing software, this can be converted to STL format, but that generates triangle and this mesh information is lost.

Until now I have explored snappyHexMesh option from OpenFOAM but in that case, I will have hard time generating a conforming mesh for fluid and solid regions.
Then I explored SALOME mesh generation software. I converted my fluid mesh to STL format and imported that in SALOME. In this I can generate a mesh but the generated mesh is very fine (as shown below) with a large number of tetrahedral elements.
meshSalomeTetrahedral.jpg

Please suggest me methods to create mesh for this type of problem. Main criteria for the mesh are following:
1) If I can generate mainly hexahedral elements because tetrahedral results in larger number of cells and that will increase the runtime.
2) Finer mesh near the given fluid region and coarse mesh as we move away from the fluid surface.

[peterhess]

Hello!

Well I am not sure if understand the problem you are facing right way...

Anyway, I belive that openFoam dont recognize 2D Mesh, if that the problem you have here.

I suggest trying a 3D Mesh for Fluid. At least one cell in the third direction.

Regards

Peter

[anufagbemi]

I'm having similar problems. i tried creating a mesh around my fluid domain in salome. But I'm having troubles merging both of them together and as a result can't obtain the solid mesh.
Both meshes are unstructured 3D meshes. Please if you were successful in yours let me know.

[peterhess]

Hello!

Post your case please!

Regards

Peter

[ashish.svm]

Originally Posted by anufagbemi

I'm having similar problems. i tried creating a mesh around my fluid domain in salome. But I'm having troubles merging both of them together and as a result can't obtain the solid mesh.
Both meshes are unstructured 3D meshes. Please if you were successful in yours let me know.

Yes.
I used SALOME and then With the "Partition" option, I merged the two domains and then meshed the partitioned geometry.

[anufagbemi]

It's too big to upload. You can send me your email or something so i could send you the trisurface files. (I also separated the solid from the general case folder).
Fluid domain and solid domains are just inversions of the other done with Avizo. I tried to use the snappyHexMesh to mesh them separately then run them. However anytime I try to run it, It tells me this;

--> FOAM FATAL ERROR: 
Master point addressing is not correct

    From function ExtendedGGIInterpolation::masterToSlavePointInterpolate(const Field<Type> pf)
    in file /home/anu/foam/foam-extend-4.0/src/foam/lnInclude/GGIInterpolation.C at line 492.

FOAM aborting

Aborted (core dumped)

I'm guessing the reason for this is that the meshes are not continuous at some points at the interface.
Now I'm attempting to use salome.

N.B the Allrun file has not yet been edited, i.e. you'll need to run the commands manually. Thanks for your help.

Samuel

[anufagbemi]

@ashish.svm It seems the partition is only available in "geometry" mode. My fluid domain is a mesh and is very difficult to read it as a geometry because it has ridiculous amount of faces.

[anufagbemi]

Originally Posted by ashish.svm

Yes.
I used SALOME and then With the "Partition" option, I merged the two domains and then meshed the partitioned geometry.

Do you mind if you share your case file. I want to see how you obtained your interface for both phases.

[ashish.svm]

Originally Posted by anufagbemi

Do you mind if you share your case file. I want to see how you obtained your interface for both phases.

Sorry for late reply. I have attached my case to the following link
https://www.dropbox.com/s/rd6phx4cw5...lex60.hdf?dl=0

[ashish.svm]

Originally Posted by anufagbemi

@ashish.svm It seems the partition is only available in "geometry" mode. My fluid domain is a mesh and is very difficult to read it as a geometry because it has ridiculous amount of faces.

I had the same problem. My fluid domain was also a mesh. But mine was very simple and it consisted of only plane surfaces so I wrote a python code (for SALOME) to create geometry. Code to create geometry is following:

# -*- coding: utf-8 -*-

###
### This file is generated automatically by SALOME v8.3.0 with dump python functionality
###

import sys
import salome

salome.salome_init()
theStudy = salome.myStudy

import salome_notebook
notebook = salome_notebook.NoteBook(theStudy)
sys.path.insert( 0, r'E:/SALOME-8.3.0-WIN64/WORK/DDMGeom')

###
### GEOM component
###

import GEOM
from salome.geom import geomBuilder
import math
import SALOMEDS


geompy = geomBuilder.New(theStudy)

Vertex_c0001 = geompy.MakeVertex(0, 15, 0)
Vertex_c0002 = geompy.MakeVertex(0.5, 15, 0)
Vertex_c0003 = geompy.MakeVertex(-0.5, 15, 0)
Vertex_c0004 = geompy.MakeVertex(1, 15, 0)
Vertex_c0005 = geompy.MakeVertex(-1, 15, 0)
Vertex_c0006 = geompy.MakeVertex(1.192, 15.2432, 0)
Vertex_c0007 = geompy.MakeVertex(1.308, 14.7568, 0)
Vertex_c0008 = geompy.MakeVertex(-1.5, 15.0002, 0)
Vertex_c0009 = geompy.MakeVertex(1.0759, 15.7295, 0)
Vertex_c0010 = geompy.MakeVertex(1.4656, 14.2844, 0)
Vertex_c0011 = geompy.MakeVertex(-1.7481, 15.2504, 0)
Vertex_c0012 = geompy.MakeVertex(-1.7519, 14.7504, 0)
Vertex_c0013 = geompy.MakeVertex(1.4241, 13.8487, 0)
Vertex_c0014 = geompy.MakeVertex(0.95991, 16.2159, 0)
Vertex_c0015 = geompy.MakeVertex(-1.7442, 15.7503, 0)
Vertex_c0016 = geompy.MakeVertex(1.011, 16.684, 0)
Vertex_c0017 = geompy.MakeVertex(-1.7403, 16.2503, 0)
Vertex_c0018 = geompy.MakeVertex(1.2293, 17.1338, 0)
Vertex_c0019 = geompy.MakeVertex(-1.7364, 16.7503, 0)
Vertex_c0020 = geompy.MakeVertex(1.4476, 17.5836, 0)
Vertex_c0021 = geompy.MakeVertex(-1.5361, 17.1526, 0)
Vertex_c0022 = geompy.MakeVertex(1.6659, 18.0335, 0)
Vertex_c0023 = geompy.MakeVertex(1.9597, 18.4269, 0)
Vertex_c0024 = geompy.MakeVertex(2.3447, 18.745, 0)
Vertex_c0025 = geompy.MakeVertex(-1.7558, 14.2504, 0)
Vertex_c0026 = geompy.MakeVertex(2.7661, 19.0113, 0)
Vertex_c0027 = geompy.MakeVertex(-1.7351, 13.7514, 0)
Vertex_c0028 = geompy.MakeVertex(-1.6654, 13.2569, 0)
Vertex_c0029 = geompy.MakeVertex(-1.4192, 13.1625, 0)
Vertex_c0030 = geompy.MakeVertex(-1.8175, 12.8603, 0)
Vertex_c0031 = geompy.MakeVertex(3.2239, 19.2085, 0)
Vertex_c0032 = geompy.MakeVertex(3.4478, 19.5385, 0)
Vertex_c0033 = geompy.MakeVertex(3.4734, 19.0391, 0)
Vertex_c0034 = geompy.MakeVertex(3.4222, 20.0378, 0)
Vertex_c0035 = geompy.MakeVertex(3.3966, 20.5372, 0)
Vertex_c0036 = geompy.MakeVertex(3.371, 21.0365, 0)
Vertex_c0037 = geompy.MakeVertex(3.3455, 21.5359, 0)
Vertex_c0038 = geompy.MakeVertex(3.3179, 22.0351, 0)
Vertex_c0039 = geompy.MakeVertex(3.2846, 22.534, 0)
Vertex_c0040 = geompy.MakeVertex(3.2418, 23.0321, 0)
Vertex_c0041 = geompy.MakeVertex(3.2537, 23.5283, 0)
Vertex_c0042 = geompy.MakeVertex(3.3264, 24.023, 0)
Vertex_c0043 = geompy.MakeVertex(3.399, 24.5177, 0)
Vertex_c0044 = geompy.MakeVertex(3.6281, 24.9241, 0)
Vertex_c0045 = geompy.MakeVertex(4.0138, 25.2424, 0)
Vertex_c0046 = geompy.MakeVertex(4.3994, 25.5607, 0)
Vertex_c0047 = geompy.MakeVertex(4.785, 25.8789, 0)
Vertex_c0048 = geompy.MakeVertex(4.8823, 26.2691, 0)
Vertex_c0049 = geompy.MakeVertex(5.0734, 25.8071, 0)


Normal_c0001 = geompy.MakeVectorDXDYDZ(6.1232e-17, 1, 0)
Normal_c0002 = geompy.MakeVectorDXDYDZ(6.1232e-17, 1, 0)
Normal_c0003 = geompy.MakeVectorDXDYDZ(6.1232e-17, 1, 0)
Normal_c0004 = geompy.MakeVectorDXDYDZ(6.1232e-17, 1, 0)
Normal_c0005 = geompy.MakeVectorDXDYDZ(6.1232e-17, 1, 0)
Normal_c0006 = geompy.MakeVectorDXDYDZ(-0.9727, -0.23207, 0)
Normal_c0007 = geompy.MakeVectorDXDYDZ(-0.9727, -0.23207, 0)
Normal_c0008 = geompy.MakeVectorDXDYDZ(0.00073192, 1, 0)
Normal_c0009 = geompy.MakeVectorDXDYDZ(-0.9727, -0.23207, 0)
Normal_c0010 = geompy.MakeVectorDXDYDZ(-0.9172, -0.39843, 0)
Normal_c0011 = geompy.MakeVectorDXDYDZ(-0.99997, 0.0077962, 0)
Normal_c0012 = geompy.MakeVectorDXDYDZ(-0.99997, 0.0077962, 0)
Normal_c0013 = geompy.MakeVectorDXDYDZ(-0.82545, 0.56447, 0)
Normal_c0014 = geompy.MakeVectorDXDYDZ(-0.9727, -0.23207, 0)
Normal_c0015 = geompy.MakeVectorDXDYDZ(-0.99997, 0.0077962, 0)
Normal_c0016 = geompy.MakeVectorDXDYDZ(-0.89966, 0.43658, 0)
Normal_c0017 = geompy.MakeVectorDXDYDZ(-0.99997, 0.0077962, 0)
Normal_c0018 = geompy.MakeVectorDXDYDZ(-0.89966, 0.43658, 0)
Normal_c0019 = geompy.MakeVectorDXDYDZ(-0.99997, 0.0077962, 0)
Normal_c0020 = geompy.MakeVectorDXDYDZ(-0.89966, 0.43658, 0)
Normal_c0021 = geompy.MakeVectorDXDYDZ(-0.60916, 0.79305, 0)
Normal_c0022 = geompy.MakeVectorDXDYDZ(-0.89966, 0.43658, 0)
Normal_c0023 = geompy.MakeVectorDXDYDZ(-0.67429, 0.73847, 0)
Normal_c0024 = geompy.MakeVectorDXDYDZ(-0.5978, 0.80165, 0)
Normal_c0025 = geompy.MakeVectorDXDYDZ(-0.99997, 0.0077962, 0)
Normal_c0026 = geompy.MakeVectorDXDYDZ(-0.46741, 0.88404, 0)
Normal_c0027 = geompy.MakeVectorDXDYDZ(-0.99588, -0.090682, 0)
Normal_c0028 = geompy.MakeVectorDXDYDZ(-0.98212, -0.18823, 0)
Normal_c0029 = geompy.MakeVectorDXDYDZ(-0.60439, 0.79669, 0)
Normal_c0030 = geompy.MakeVectorDXDYDZ(-0.60439, 0.79669, 0)
Normal_c0031 = geompy.MakeVectorDXDYDZ(-0.32137, 0.94695, 0)
Normal_c0032 = geompy.MakeVectorDXDYDZ(-0.99869, -0.051181, 0)
Normal_c0033 = geompy.MakeVectorDXDYDZ(-0.99869, -0.051181, 0)
Normal_c0034 = geompy.MakeVectorDXDYDZ(-0.99869, -0.051181, 0)
Normal_c0035 = geompy.MakeVectorDXDYDZ(-0.99869, -0.051181, 0)
Normal_c0036 = geompy.MakeVectorDXDYDZ(-0.99869, -0.051181, 0)
Normal_c0037 = geompy.MakeVectorDXDYDZ(-0.99869, -0.051181, 0)
Normal_c0038 = geompy.MakeVectorDXDYDZ(-0.99825, -0.059095, 0)
Normal_c0039 = geompy.MakeVectorDXDYDZ(-0.99726, -0.073937, 0)
Normal_c0040 = geompy.MakeVectorDXDYDZ(-0.99524, -0.09749, 0)
Normal_c0041 = geompy.MakeVectorDXDYDZ(-0.98939, 0.14528, 0)
Normal_c0042 = geompy.MakeVectorDXDYDZ(-0.98939, 0.14528, 0)
Normal_c0043 = geompy.MakeVectorDXDYDZ(-0.98939, 0.14528, 0)
Normal_c0044 = geompy.MakeVectorDXDYDZ(-0.63654, 0.77125, 0)
Normal_c0045 = geompy.MakeVectorDXDYDZ(-0.63654, 0.77125, 0)
Normal_c0046 = geompy.MakeVectorDXDYDZ(-0.63654, 0.77125, 0)
Normal_c0047 = geompy.MakeVectorDXDYDZ(-0.63654, 0.77125, 0)
Normal_c0048 = geompy.MakeVectorDXDYDZ(-0.92407, -0.38223, 0)
Normal_c0049 = geompy.MakeVectorDXDYDZ(-0.92407, -0.38223, 0)


Plane_c0001 = geompy.MakePlane(Vertex_c0001, Normal_c0001, 0.5)
Plane_c0002 = geompy.MakePlane(Vertex_c0002, Normal_c0002, 0.5)
Plane_c0003 = geompy.MakePlane(Vertex_c0003, Normal_c0003, 0.5)
Plane_c0004 = geompy.MakePlane(Vertex_c0004, Normal_c0004, 0.5)
Plane_c0005 = geompy.MakePlane(Vertex_c0005, Normal_c0005, 0.5)
Plane_c0006 = geompy.MakePlane(Vertex_c0006, Normal_c0006, 0.5)
Plane_c0007 = geompy.MakePlane(Vertex_c0007, Normal_c0007, 0.5)
Plane_c0008 = geompy.MakePlane(Vertex_c0008, Normal_c0008, 0.5)
Plane_c0009 = geompy.MakePlane(Vertex_c0009, Normal_c0009, 0.5)
Plane_c0010 = geompy.MakePlane(Vertex_c0010, Normal_c0010, 0.5)
Plane_c0011 = geompy.MakePlane(Vertex_c0011, Normal_c0011, 0.5)
Plane_c0012 = geompy.MakePlane(Vertex_c0012, Normal_c0012, 0.5)
Plane_c0013 = geompy.MakePlane(Vertex_c0013, Normal_c0013, 0.5)
Plane_c0014 = geompy.MakePlane(Vertex_c0014, Normal_c0014, 0.5)
Plane_c0015 = geompy.MakePlane(Vertex_c0015, Normal_c0015, 0.5)
Plane_c0016 = geompy.MakePlane(Vertex_c0016, Normal_c0016, 0.5)
Plane_c0017 = geompy.MakePlane(Vertex_c0017, Normal_c0017, 0.5)
Plane_c0018 = geompy.MakePlane(Vertex_c0018, Normal_c0018, 0.5)
Plane_c0019 = geompy.MakePlane(Vertex_c0019, Normal_c0019, 0.5)
Plane_c0020 = geompy.MakePlane(Vertex_c0020, Normal_c0020, 0.5)
Plane_c0021 = geompy.MakePlane(Vertex_c0021, Normal_c0021, 0.5)
Plane_c0022 = geompy.MakePlane(Vertex_c0022, Normal_c0022, 0.5)
Plane_c0023 = geompy.MakePlane(Vertex_c0023, Normal_c0023, 0.5)
Plane_c0024 = geompy.MakePlane(Vertex_c0024, Normal_c0024, 0.5)
Plane_c0025 = geompy.MakePlane(Vertex_c0025, Normal_c0025, 0.5)
Plane_c0026 = geompy.MakePlane(Vertex_c0026, Normal_c0026, 0.5)
Plane_c0027 = geompy.MakePlane(Vertex_c0027, Normal_c0027, 0.5)
Plane_c0028 = geompy.MakePlane(Vertex_c0028, Normal_c0028, 0.5)
Plane_c0029 = geompy.MakePlane(Vertex_c0029, Normal_c0029, 0.5)
Plane_c0030 = geompy.MakePlane(Vertex_c0030, Normal_c0030, 0.5)
Plane_c0031 = geompy.MakePlane(Vertex_c0031, Normal_c0031, 0.5)
Plane_c0032 = geompy.MakePlane(Vertex_c0032, Normal_c0032, 0.5)
Plane_c0033 = geompy.MakePlane(Vertex_c0033, Normal_c0033, 0.5)
Plane_c0034 = geompy.MakePlane(Vertex_c0034, Normal_c0034, 0.5)
Plane_c0035 = geompy.MakePlane(Vertex_c0035, Normal_c0035, 0.5)
Plane_c0036 = geompy.MakePlane(Vertex_c0036, Normal_c0036, 0.5)
Plane_c0037 = geompy.MakePlane(Vertex_c0037, Normal_c0037, 0.5)
Plane_c0038 = geompy.MakePlane(Vertex_c0038, Normal_c0038, 0.5)
Plane_c0039 = geompy.MakePlane(Vertex_c0039, Normal_c0039, 0.5)
Plane_c0040 = geompy.MakePlane(Vertex_c0040, Normal_c0040, 0.5)
Plane_c0041 = geompy.MakePlane(Vertex_c0041, Normal_c0041, 0.5)
Plane_c0042 = geompy.MakePlane(Vertex_c0042, Normal_c0042, 0.5)
Plane_c0043 = geompy.MakePlane(Vertex_c0043, Normal_c0043, 0.5)
Plane_c0044 = geompy.MakePlane(Vertex_c0044, Normal_c0044, 0.5)
Plane_c0045 = geompy.MakePlane(Vertex_c0045, Normal_c0045, 0.5)
Plane_c0046 = geompy.MakePlane(Vertex_c0046, Normal_c0046, 0.5)
Plane_c0047 = geompy.MakePlane(Vertex_c0047, Normal_c0047, 0.5)
Plane_c0048 = geompy.MakePlane(Vertex_c0048, Normal_c0048, 0.5)
Plane_c0049 = geompy.MakePlane(Vertex_c0049, Normal_c0049, 0.5)


geompy.addToStudy(Vertex_c0001,'Vertex_c0001')
geompy.addToStudy(Vertex_c0002,'Vertex_c0002')
geompy.addToStudy(Vertex_c0003,'Vertex_c0003')
geompy.addToStudy(Vertex_c0004,'Vertex_c0004')
geompy.addToStudy(Vertex_c0005,'Vertex_c0005')
geompy.addToStudy(Vertex_c0006,'Vertex_c0006')
geompy.addToStudy(Vertex_c0007,'Vertex_c0007')
geompy.addToStudy(Vertex_c0008,'Vertex_c0008')
geompy.addToStudy(Vertex_c0009,'Vertex_c0009')
geompy.addToStudy(Vertex_c0010,'Vertex_c0010')
geompy.addToStudy(Vertex_c0011,'Vertex_c0011')
geompy.addToStudy(Vertex_c0012,'Vertex_c0012')
geompy.addToStudy(Vertex_c0013,'Vertex_c0013')
geompy.addToStudy(Vertex_c0014,'Vertex_c0014')
geompy.addToStudy(Vertex_c0015,'Vertex_c0015')
geompy.addToStudy(Vertex_c0016,'Vertex_c0016')
geompy.addToStudy(Vertex_c0017,'Vertex_c0017')
geompy.addToStudy(Vertex_c0018,'Vertex_c0018')
geompy.addToStudy(Vertex_c0019,'Vertex_c0019')
geompy.addToStudy(Vertex_c0020,'Vertex_c0020')
geompy.addToStudy(Vertex_c0021,'Vertex_c0021')
geompy.addToStudy(Vertex_c0022,'Vertex_c0022')
geompy.addToStudy(Vertex_c0023,'Vertex_c0023')
geompy.addToStudy(Vertex_c0024,'Vertex_c0024')
geompy.addToStudy(Vertex_c0025,'Vertex_c0025')
geompy.addToStudy(Vertex_c0026,'Vertex_c0026')
geompy.addToStudy(Vertex_c0027,'Vertex_c0027')
geompy.addToStudy(Vertex_c0028,'Vertex_c0028')
geompy.addToStudy(Vertex_c0029,'Vertex_c0029')
geompy.addToStudy(Vertex_c0030,'Vertex_c0030')
geompy.addToStudy(Vertex_c0031,'Vertex_c0031')
geompy.addToStudy(Vertex_c0032,'Vertex_c0032')
geompy.addToStudy(Vertex_c0033,'Vertex_c0033')
geompy.addToStudy(Vertex_c0034,'Vertex_c0034')
geompy.addToStudy(Vertex_c0035,'Vertex_c0035')
geompy.addToStudy(Vertex_c0036,'Vertex_c0036')
geompy.addToStudy(Vertex_c0037,'Vertex_c0037')
geompy.addToStudy(Vertex_c0038,'Vertex_c0038')
geompy.addToStudy(Vertex_c0039,'Vertex_c0039')
geompy.addToStudy(Vertex_c0040,'Vertex_c0040')
geompy.addToStudy(Vertex_c0041,'Vertex_c0041')
geompy.addToStudy(Vertex_c0042,'Vertex_c0042')
geompy.addToStudy(Vertex_c0043,'Vertex_c0043')
geompy.addToStudy(Vertex_c0044,'Vertex_c0044')
geompy.addToStudy(Vertex_c0045,'Vertex_c0045')
geompy.addToStudy(Vertex_c0046,'Vertex_c0046')
geompy.addToStudy(Vertex_c0047,'Vertex_c0047')
geompy.addToStudy(Vertex_c0048,'Vertex_c0048')
geompy.addToStudy(Vertex_c0049,'Vertex_c0049')


geompy.addToStudy(Normal_c0001,'Normal_c0001')
geompy.addToStudy(Normal_c0002,'Normal_c0002')
geompy.addToStudy(Normal_c0003,'Normal_c0003')
geompy.addToStudy(Normal_c0004,'Normal_c0004')
geompy.addToStudy(Normal_c0005,'Normal_c0005')
geompy.addToStudy(Normal_c0006,'Normal_c0006')
geompy.addToStudy(Normal_c0007,'Normal_c0007')
geompy.addToStudy(Normal_c0008,'Normal_c0008')
geompy.addToStudy(Normal_c0009,'Normal_c0009')
geompy.addToStudy(Normal_c0010,'Normal_c0010')
geompy.addToStudy(Normal_c0011,'Normal_c0011')
geompy.addToStudy(Normal_c0012,'Normal_c0012')
geompy.addToStudy(Normal_c0013,'Normal_c0013')
geompy.addToStudy(Normal_c0014,'Normal_c0014')
geompy.addToStudy(Normal_c0015,'Normal_c0015')
geompy.addToStudy(Normal_c0016,'Normal_c0016')
geompy.addToStudy(Normal_c0017,'Normal_c0017')
geompy.addToStudy(Normal_c0018,'Normal_c0018')
geompy.addToStudy(Normal_c0019,'Normal_c0019')
geompy.addToStudy(Normal_c0020,'Normal_c0020')
geompy.addToStudy(Normal_c0021,'Normal_c0021')
geompy.addToStudy(Normal_c0022,'Normal_c0022')
geompy.addToStudy(Normal_c0023,'Normal_c0023')
geompy.addToStudy(Normal_c0024,'Normal_c0024')
geompy.addToStudy(Normal_c0025,'Normal_c0025')
geompy.addToStudy(Normal_c0026,'Normal_c0026')
geompy.addToStudy(Normal_c0027,'Normal_c0027')
geompy.addToStudy(Normal_c0028,'Normal_c0028')
geompy.addToStudy(Normal_c0029,'Normal_c0029')
geompy.addToStudy(Normal_c0030,'Normal_c0030')
geompy.addToStudy(Normal_c0031,'Normal_c0031')
geompy.addToStudy(Normal_c0032,'Normal_c0032')
geompy.addToStudy(Normal_c0033,'Normal_c0033')
geompy.addToStudy(Normal_c0034,'Normal_c0034')
geompy.addToStudy(Normal_c0035,'Normal_c0035')
geompy.addToStudy(Normal_c0036,'Normal_c0036')
geompy.addToStudy(Normal_c0037,'Normal_c0037')
geompy.addToStudy(Normal_c0038,'Normal_c0038')
geompy.addToStudy(Normal_c0039,'Normal_c0039')
geompy.addToStudy(Normal_c0040,'Normal_c0040')
geompy.addToStudy(Normal_c0041,'Normal_c0041')
geompy.addToStudy(Normal_c0042,'Normal_c0042')
geompy.addToStudy(Normal_c0043,'Normal_c0043')
geompy.addToStudy(Normal_c0044,'Normal_c0044')
geompy.addToStudy(Normal_c0045,'Normal_c0045')
geompy.addToStudy(Normal_c0046,'Normal_c0046')
geompy.addToStudy(Normal_c0047,'Normal_c0047')
geompy.addToStudy(Normal_c0048,'Normal_c0048')
geompy.addToStudy(Normal_c0049,'Normal_c0049')


geompy.addToStudy(Plane_c0001,'Plane_c0001')
geompy.addToStudy(Plane_c0002,'Plane_c0002')
geompy.addToStudy(Plane_c0003,'Plane_c0003')
geompy.addToStudy(Plane_c0004,'Plane_c0004')
geompy.addToStudy(Plane_c0005,'Plane_c0005')
geompy.addToStudy(Plane_c0006,'Plane_c0006')
geompy.addToStudy(Plane_c0007,'Plane_c0007')
geompy.addToStudy(Plane_c0008,'Plane_c0008')
geompy.addToStudy(Plane_c0009,'Plane_c0009')
geompy.addToStudy(Plane_c0010,'Plane_c0010')
geompy.addToStudy(Plane_c0011,'Plane_c0011')
geompy.addToStudy(Plane_c0012,'Plane_c0012')
geompy.addToStudy(Plane_c0013,'Plane_c0013')
geompy.addToStudy(Plane_c0014,'Plane_c0014')
geompy.addToStudy(Plane_c0015,'Plane_c0015')
geompy.addToStudy(Plane_c0016,'Plane_c0016')
geompy.addToStudy(Plane_c0017,'Plane_c0017')
geompy.addToStudy(Plane_c0018,'Plane_c0018')
geompy.addToStudy(Plane_c0019,'Plane_c0019')
geompy.addToStudy(Plane_c0020,'Plane_c0020')
geompy.addToStudy(Plane_c0021,'Plane_c0021')
geompy.addToStudy(Plane_c0022,'Plane_c0022')
geompy.addToStudy(Plane_c0023,'Plane_c0023')
geompy.addToStudy(Plane_c0024,'Plane_c0024')
geompy.addToStudy(Plane_c0025,'Plane_c0025')
geompy.addToStudy(Plane_c0026,'Plane_c0026')
geompy.addToStudy(Plane_c0027,'Plane_c0027')
geompy.addToStudy(Plane_c0028,'Plane_c0028')
geompy.addToStudy(Plane_c0029,'Plane_c0029')
geompy.addToStudy(Plane_c0030,'Plane_c0030')
geompy.addToStudy(Plane_c0031,'Plane_c0031')
geompy.addToStudy(Plane_c0032,'Plane_c0032')
geompy.addToStudy(Plane_c0033,'Plane_c0033')
geompy.addToStudy(Plane_c0034,'Plane_c0034')
geompy.addToStudy(Plane_c0035,'Plane_c0035')
geompy.addToStudy(Plane_c0036,'Plane_c0036')
geompy.addToStudy(Plane_c0037,'Plane_c0037')
geompy.addToStudy(Plane_c0038,'Plane_c0038')
geompy.addToStudy(Plane_c0039,'Plane_c0039')
geompy.addToStudy(Plane_c0040,'Plane_c0040')
geompy.addToStudy(Plane_c0041,'Plane_c0041')
geompy.addToStudy(Plane_c0042,'Plane_c0042')
geompy.addToStudy(Plane_c0043,'Plane_c0043')
geompy.addToStudy(Plane_c0044,'Plane_c0044')
geompy.addToStudy(Plane_c0045,'Plane_c0045')
geompy.addToStudy(Plane_c0046,'Plane_c0046')
geompy.addToStudy(Plane_c0047,'Plane_c0047')
geompy.addToStudy(Plane_c0048,'Plane_c0048')
geompy.addToStudy(Plane_c0049,'Plane_c0049')


Compound_c0001 = geompy.MakeCompound([Plane_c0001, \
Plane_c0002, \
Plane_c0003, \
Plane_c0004, \
Plane_c0005, \
Plane_c0006, \
Plane_c0007, \
Plane_c0008, \
Plane_c0009, \
Plane_c0010, \
Plane_c0011, \
Plane_c0012, \
Plane_c0013, \
Plane_c0014, \
Plane_c0015, \
Plane_c0016, \
Plane_c0017, \
Plane_c0018, \
Plane_c0019, \
Plane_c0020, \
Plane_c0021, \
Plane_c0022, \
Plane_c0023, \
Plane_c0024, \
Plane_c0025, \
Plane_c0026, \
Plane_c0027, \
Plane_c0028, \
Plane_c0029, \
Plane_c0030, \
Plane_c0031, \
Plane_c0032, \
Plane_c0033, \
Plane_c0034, \
Plane_c0035, \
Plane_c0036, \
Plane_c0037, \
Plane_c0038, \
Plane_c0039, \
Plane_c0040, \
Plane_c0041, \
Plane_c0042, \
Plane_c0043, \
Plane_c0044, \
Plane_c0045, \
Plane_c0046, \
Plane_c0047, \
Plane_c0048, \
Plane_c0049, \
])

geompy.addToStudy( Compound_c0001,'Compound_c0001')


Glue_c0001 = geompy.MakeGlueEdges(Compound_c0001, 0.001)
geompy.addToStudy( Glue_c0001,'Glue_c0001')
geompy.ExportIGES(Glue_c0001,'./DFNGeometry.iges')

[mhgolestan]

Originally Posted by anufagbemi

I'm having similar problems. i tried creating a mesh around my fluid domain in salome. But I'm having troubles merging both of them together and as a result can't obtain the solid mesh.
Both meshes are unstructured 3D meshes. Please if you were successful in yours let me know.

Hi @anufagbemi
I want to create a mesh similar to yours (porous structure), but faced the same error. Could you solve this issue?