Я работаю над связанной моделью DEM-CFD. Я использую Fipy для CFD, а моя геометрия - цилиндр. Как я могу импортировать созданную сетку в Gmsh в Fipy?
Мне нужно сгенерировать цилиндрическую геометрию и дискретизировать ее путем создания треугольной сетки в 3D в Fipy. Я попытался найти команду с помощью Fipy для создания трехмерной треугольной сетки в цилиндрической геометрии, но я не смог найти ничего подходящего для 3D. Я создал программу Mesh in Gmsh, которая была принята Fipy на основании его помощи. У меня сейчас есть файл сетки, но я не могу импортировать его в Fipy.
Всякий раз, когда я пытался импортировать его в модель, я сталкивался с этой ошибкой Gmsh version must be >= 2.0
. Тем не менее, моя версия Gmsh была 4.
Я также использовал Gmsh V2.5, но ошибка была той же. Я искал ошибку, и я вижу, что некоторые люди тоже боролись с этим. Для решения проблемы рекомендуется поместить файл Gmsh.exe в папку скриптов. Я сделал это, но у меня есть другая ошибка. Новая ошибка
Windows Error: [Error 32] The process cannot access the file because it is being used by another process
когда команда была mesh = fp.Gmsh3D("Cylinder.msh")
.
Я также искал эту ошибку, и кто-то предложил упомянуть файл геометрии после сетки, mesh = fp.Gmsh3D("Cylinder.msh, Cylinder.geo")
.
Теперь я получаю это
GmshException: Gmsh hasn't produced any cells! Check your Gmsh code.
Я попытался сделать круговую сетку, которая является двухмерной треугольной сеткой в круге (это в справке Fipy), но я получил это
GmshException: Gmsh hasn't produced any cells! Check your Gmsh code
Файлы python scritp, цилиндр.msh, цилиндр.geo и Gmsh.exe находятся в одной папке.
Ошибка:
File "C:\Users\jpethati\AppData\Local\Continuum\anaconda2\lib\site-packages\fipy\meshes\gmshMesh.py", line 809, in read
raise GmshException(errStr)
GmshException: Gmsh hasn't produced any cells! Check your Gmsh code.
Гмш выход:
Info : Running 'gmsh c:\users\jpethati\appdata\local\temp\tmpsw2tzg.geo -3 -nopopup -format msh -o c:\users\jpethati\appdata\local\temp\tmp6gbww9.msh' [Gmsh 4.2.2, 1 node, max. 1 thread]
Info : Started on Fri Mar 29 10:25:25 2019
Info : Reading 'c:\users\jpethati\appdata\local\temp\tmpsw2tzg.geo'...
Info : Done reading 'c:\users\jpethati\appdata\local\temp\tmpsw2tzg.geo'
Info : Meshing 1D...
Info : Done meshing 1D (0 s)
Info : Meshing 2D...
Info : Done meshing 2D (0 s)
Info : Meshing 3D...
Info : Done meshing 3D (0 s)
Info : 0 vertices 0 elements
Info : Writing 'c:\users\jpethati\appdata\local\temp\tmp6gbww9.msh'...
Info : Done writing 'c:\users\jpethati\appdata\local\temp\tmp6gbww9.msh'
Info : Stopped on Fri Mar 29 10:25:25 2019
import fipy as fp
Первая попытка
mesh = fp.Gmsh3D("Cylinder.msh")
Вторая попытка
mesh = fp.Gmsh3D("Cylinder.msh, Cylinder.geo")