[Solved] Don't get it running...

[Elektrix]

Hello everybody,

this is my first post here. At the moment I'm doing my first steps with Code_Saturne 3.3. My problem is that the solver stops right at the beginning of the calculation.

My system is Debian Jessie AMD64 with Code_Saturne 3.3 installed from the Debian repositories.

I want to do some calculation about what happens in a water storage tank.

So here is what I did:

I constructed the tank in Salome 7.4.0 (only the water volume for starting) and meshed it inside Salome with Netgen 1D2D3D and then saved the mesh in a med file. The mesh consists of the volume mesh and a 2D-submesh that contains the surface. For the beginning I have constant temperature in the whole tank to see if the calculation works.

The surface mesh is the boundary region.

I've tried serveral different settings but no matter what I do I always get a

Code: Select all

SIGFPE signal (floating point exception) intercepted!

In an older post I read that this could be due to errors in the mesh but Code_Saturne says:

Code: Select all

  Criterion 1: Orthogonality:
    Number of bad cells detected: 0 -->   0 %

  Criterion 2: Offset:
    Number of bad cells detected: 0 -->   0 %

  Criterion 3: Least-Squares Gradient Quality:
    Number of bad cells detected: 2 -->   0 %

  Criterion 4: Cells Volume Ratio:
    Number of bad cells detected: 0 -->   0 %

  Criterion 5: Guilt by Association:
    Number of bad cells detected: 0 -->   0 %

 Warning:
 --------
    Mesh quality issue has been detected

    The mesh should be re-considered using the listed criteria.

    The calculation will run but the solution quality may be degraded...

But the calculation doesn't run

Can anybody give me further information about this?

I attach the listing for further details (renamed it to listing.txt).

Thanks in advance

Elektrix

[Yvan Fournier]

Hello,

The stacktrace at the end of your "listing" file indicates you have a floating point exception in a "power" operator in the computation of physical properties, so there may be something wrong in your physical property definitions (for example computing the power of a negative value).

You should check your physical property definitions (if you set some properties as "user defined", make sure you also edit the property definitions themselves).

Regards,

Yvan

[Elektrix]

Hello Yvan,

thank you very much for your quick and helpful answer

I had several errors and so the error messages changed while correcting the mistakes:

The Thermal model was set to Temperature Celsius and the initial temperature was 300, corrected it to Kelvin
The formula for the density was wrong (typo), replaced it for test by the constant value, now works also with the corrected formula.
The name of the submesh for the boundary had spaces, replaced them with "_".

Now the calculation is running and I get a result in paravis.

Thanks again very much!

Elektrix