Welcome, Guest
Username: Password: Remember me

TOPIC: Boundary conditions file does not match geometry file

Boundary conditions file does not match geometry file 11 years 8 months ago #7961

  • Rho
  • Rho's Avatar
Hello,

When I want to run a simulation, TELEMAC-2D displays an error message (see the attached file)

Can you help me?

thank you

Rho
Attachments:
The administrator has disabled public write access.

Boundary conditions file does not match geometry file 11 years 8 months ago #7962

  • jmhervouet
  • jmhervouet's Avatar
Hello,

There is obviously something wrong in your mesh. The element number printed is impossible. Try to see the mesh with a post-processor.

Regards,

JMH
The administrator has disabled public write access.

Boundary conditions file does not match geometry file 11 years 8 months ago #7963

  • Rho
  • Rho's Avatar
Hello,

I do not understand because normally when something is wrong in a mesh Fudaa displays an error message. But this is not the case here.

I do not know where to look for the error in my element mesh knowing the number printed is actually non-existent ...

Do you have any other suggestions?

Rho
The administrator has disabled public write access.

Boundary conditions file does not match geometry file 11 years 8 months ago #7964

  • jmhervouet
  • jmhervouet's Avatar
Hello,

The other possibility is that the message is correct and that the boundary conditions file is not compatible with the mesh.

The fact that Fudaa works is not a proof, it can display a mesh that is not really a finite element mesh. For example two triangles connected by one point is not a finite element mesh: nothing can flow from one triangle to the other, and the connection point belongs to both triangles and is influenced by both, it would mean that for any physical quantity it has one value for a triangle and another value for the second. It could be that somewhere you have such a situation, a triangle connected to the rest of the mesh by a single point ? If I look at the algorithm that complains, it seems that no element sees a given boundary point.

If you know Fortran you can take subroutine elebd.f (library bief) in your FORTRAN FILE and where the message is written (line with WRITE(LU,810) IEL) you can add a line :

WRITE(LU,*) 'K1=',K1,' N1=',N1

it will give you the guilty boundary point number, and its global number, so that you can see where it is. I'll add this extra information for everybody in the next version, it will be easier to locate a problem.


Regards,

JMH
The administrator has disabled public write access.
Moderators: pham

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.