Конвертация сетки из Саломы в OpenFOAM

[nFirth 0]

Здравствуйте. Случилась такая вот неприятность, требуется конвертировать сетку из Саломы в OpenFOAM 1)создаю файл .unv 2) нибираю в терминале

ideasUnvToFoam Mesh_1.unv

Получаю в ответ ошибку сегментирования

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : ideasUnvToFoam Mesh_1.unv
Date   : Aug 13 2012
Time   : 11:25:42
Host   : "DN"
PID    : 4495
Case   : /home/dn/OpenFOAM/dn-2.1.1/PD
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Processing tag:2411
Starting reading points at line 3.
Read 48337 points.

Processing tag:2412
Starting reading cells at line 96680.
First occurrence of element type 11 for cell 1 at line 96681
First occurrence of element type 44 for cell 1402 at line 100884
First occurrence of element type 41 for cell 4157 at line 106394
Read 0 cells and 47681 boundary faces.

Processing tag:2467
Starting reading patches at line 196248.
For group 2 named NewBoundary trying to read 47681 patch face indices.

Sorting boundary faces according to group (patch)
0: NewBoundary is #0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigSegv::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/libc.so.6"
#3  
in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/ideasUnvToFoam"
#4  __libc_start_main in "/lib/libc.so.6"
#5  
in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/ideasUnvToFoam"
Ошибка сегментирования

Очень надеюсь на вашу помощь, с уважением.

Даниил.

[HFL 34]

На исходную сетку в I-DEAS Universal (судя по всему это поверхностная сетка) можно посмотреть ?

[nFirth 0]

PD.zip

Да конечно вот он.

[Vova 419]

PD.zip

Да конечно вот он.

Openfoam требует объёмную сетку, а у вас вроде только поверхносная. В этом может и проблема.

[HFL 34]

Openfoam требует объёмную сетку, а у вас вроде только поверхносная. В этом может и проблема.

ideasUnvToFoam конвертер по идее работает с любыми сетками (2D-3D),но в данном случае это просто больная сетка (смешанная, четырехугольники-треугольники) от которой он ломается

face 0 in patch 0 does not have neighbour cell face: 4(8 7 1400 1401

/*---------------------------------------------------------------------------*\

| =========				 |												 |

| \\	  /  F ield		 | OpenFOAM: The Open Source CFD Toolbox		   |

|  \\	/   O peration	 | Version:  1.6								   |

|   \\  /	A nd		   | Web:	  www.OpenFOAM.org					  |

|	\\/	 M anipulation  |												 |

\*---------------------------------------------------------------------------*/

Build  : 1.6-53b7f692aa41

Exec   : ideasUnvToFoam Mesh_1.unv

Date   : Aug 16 2012

Time   : 15:16:07

Host   : ACER

PID	: 20414

Case   : /home/hfl/OpenFOAM/run/Test

nProcs : 1

SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).



// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Create time



Processing tag:2411

Starting reading points at line 3.

Read 48337 points.



Processing tag:2412

Starting reading cells at line 96680.

First occurrence of element type 11 for cell 1 at line 96681

First occurrence of element type 44 for cell 1402 at line 100884

First occurrence of element type 41 for cell 4157 at line 106394

Read 0 cells and 47681 boundary faces.



Processing tag:2467

Starting reading patches at line 196248.

For facegroup 2 named NewBoundary trying to read 47681 patch face indices.



Sorting boundary faces according to group (patch)

Constructing mesh with non-default patches of size:

	NewBoundary 47681







face 0 in patch 0 does not have neighbour cell face: 4(8 7 1400 1401)#0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"

#1  Foam::error::abort() in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"

#2  Foam::Ostream& Foam::operator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/opt/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/ideasUnvToFoam"

#3  Foam::polyMesh::facePatchFaceCells(Foam::List<Foam::face> const&, Foam::List<Foam::List<int> > const&, Foam::List<Foam::List<Foam::face> > const&, int) const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"

#4  Foam::polyMesh::polyMesh(Foam::IOobject const&, Foam::Xfer<Foam::Field<Foam::Vector<double> > > const&, Foam::List<Foam::cellShape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::List<Foam::word> const&, Foam::word const&, Foam::word const&, Foam::List<Foam::word> const&, bool) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"

#5  main in "/opt/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/ideasUnvToFoam"

#6  __libc_start_main in "/lib/libc.so.6"

#7  _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122





	From function polyMesh::facePatchFaceCells(const faceList& patchFaces,const labelListList& pointCells,const faceListList& cellsFaceShapes,const label patchID)

	in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 125.



FOAM aborting

По опыту из саломовых сеток нормально конвертятся сетки построенные NETGEN 1-2(-3) и Mefisto (треугольные сетки) но OpenFOAM на них считает медленней.

OpenFOAM же действительно считает только на объёмных сетках ...

PS: Я смог затащить вашу сетку в OpenFOAM (предварительно сделав её объёмной через extrusion по оси z в той же самой Salome) только открыв её предварительно в gmsh и сохранив в .msh

Изменено 16 августа 2012 пользователем HFL

[nFirth 0]

Спасибо большое, значит сделаю геометрию объемной, натяну нормальную сетку и попробую запустить.

[vl 327]

Можете попробовать через STL с помощью snappyhex сгенерить сетку

[nFirth 0]

Отлично! Все работает продавил геометрию и сетку разбил грамотно на группы, и все посчиталось.