TELEMAC-MASCARET Forum: Mesh problems (1/1)

Welcome, Guest

TOPIC: Mesh problems

Mesh problems 11 years 7 months ago #8369

konsonaut
OFFLINE
openTELEMAC Guru
Posts: 413
Thank you received: 144

Hello!

I have some problems with running a simulation.
The mesh in SMS format was from another 2-D solver and I should use it for Telemac-2D. The maximum number of points around a node is 16 and I know so the quality of the mesh is not so good but I don't want to generate a new mesh. I generated a boundary conditions file as usual but I get the following error:

READGEO1: TITLE= newSelafin

NUMBER OF ELEMENTS: 353087
NUMBER OF POINTS: 197107
MXPTEL (BIEF) : MAXIMUM NUMBER OF ELEMENTS AROUND A POI
NT: 16
MAXIMUM NUMBER OF POINTS AROUND A POINT
: 16
CORRXY (BIEF):NO MODIFICATION OF COORDINATES

MESH: MESH ALLOCATED

*************************************
* END OF MEMORY ORGANIZATION: *
*************************************

ELEBD: ERROR OF NUMBERING IN THE ELEMENT:2147483647
POSSIBLE REASON:
THE BOUNDARY CONDITION FILE IS NOT
RELEVANT TO THE GEOMETRY FILE

PLANTE: PROGRAM STOPPED AFTER AN ERROR

What is the reason for this error message?

I would be very thankful for any help!

Clemens

The administrator has disabled public write access.

Mesh problems 11 years 7 months ago #8370

c.coulet OFFLINE Moderator Posts: 3722 Thank you received: 1031	Hi Clemens In Telemac system, the maximum number of points around a point is fixed to 10! That's the reason of this message. You could try to modify the sources to change this limit but i'm not sure this will be very easy to do. Maybe you could have a look to the stbtel program but i'm not sure if there is an option to generate a new mesh with an acceptable number of neighbours... Regards
	Christophe The administrator has disabled public write access. The following user(s) said Thank You: konsonaut

Mesh problems 11 years 7 months ago #8388

jmhervouet

Hello Clemens and Christophe,

The variables are MXPTVS (maximum of points around a point in the mesh) and MXELVS (maximum number of elements around a point in the mesh). I am not sure that there is a clear limitation of these numbers as obviously the program does not stop. The only real limitation I remember is for using the frontal product (keyword MATRIX-VECTOR PRODUCT, I think, with presumably 1 : normal and 2 : frontal, not checked). In this case we start by the points with two neighbours, then 3 neighbours, etc., and for every number of neighbours there is a loop, and we stopped at 10 neighbours (but it would be easy to write more loops for larger numbers...).

However it could be that this limitation at 10 is elsewhere in the program, but the program should stop there with a clear message.

So I would say that either the mesh has really something else weird, or we miss a clear stop somewhere before elebd (I do not see any such limitation by 10 in ELEBD, nor in VOISIN which is also a key subroutine at the beginning).

With best regards,

Jean-Michel Hervouet

The administrator has disabled public write access.

Mesh problems 11 years 6 months ago #8393

konsonaut OFFLINE openTELEMAC Guru Posts: 413 Thank you received: 144	Hello, thank you Christophe and Jean-Michel! Indeed the test simulation runs also with maximum number of points around a point in the mesh > 10, in my case 16. Hopefully the simulations will not make any troubles later. The mesh had some other issues, e.g. points which shared more than one outblocked area (holes) and I could resolve these problems with the STBTEL utility. Best regards, Clemens
	The administrator has disabled public write access.