CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Floating Point Exception - wrong boundaries or general PC problem? – OF 1.6 extend - (http://www.cfd-online.com/Forums/openfoam/113812-floating-point-exception-wrong-boundaries-general-pc-problem-1-6-extend.html)

A.Wendy February 27, 2013 04:50

Floating Point Exception - wrong boundaries or general PC problem? – OF 1.6 extend -
 
Hi Everybody,

Currently I am working on my masterthesis for a simulation of a gyrocopter including autorotation.

For that I started to modify the mixerGGI library to include a threedimensional autorotation.

I think I worked it out but when I try to start my computation I always get a floating point exception. After some time steps, when I restart the calculation I am able to compute again a number of time steps and so on.
This happens also with the not modified mixerGGI library.

I have no idea where the fault is. I think my boundaries are OK. Maybe my computer has just a general problem? I’ve got no idea.
Maybe someone can have a look or try to run my case with OF 1.6 extend and pimpleDyMFoam if the same error occurs?

Best wishes

Andreas

Link to the complete case: http://ubuntuone.com/0aasoEIvPToraOAPetJJie

0/U:
Code:

  /*--------------------------------*- C++ -*----------------------------------*\
  | =========                |                                                |
  | \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
  |  \\    /  O peration    | Version:  1.6-ext                                |
  |  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
  |    \\/    M anipulation  |                                                |
  \*---------------------------------------------------------------------------*/
  FoamFile
  {
      version    2.0;
      format      ascii;
      class      volVectorField;
      location    "0";
      object      U;
  }
  // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
  dimensions      [0 1 -1 0 0 0 0];
 
  internalField      uniform (10 0 0);
                                                                               
  boundaryField                                                               
  {                                                                           
                                                               
      outsideSlider                                                             
      {                                                                       
                  type                      ggi;                                                       
                  value                    uniform (10 0 0);
      }                                                                       
                                                               
      insideSlider                                                                 
      {                                                                       
                  type                      ggi;
                  value                    uniform (10 0 0);
      }
 
      outlet                                                                             
      {                                                                       
          type                              freestream;
                  freestreamValue            uniform (10 0 0);                           
      }                                                                       
                                                                               
      inlet                                                               
      {
          type                              freestream;
                  freestreamValue            uniform (10 0 0);
      }
 
      sides                                                               
      {                                                                       
          type                              freestream;
                  freestreamValue            uniform (10 0 0);
      }
 
      ".*STL"                                                                         
      {                                                                       
          type            movingWallVelocity;
          value                            uniform (0 0 0);
      }
 
 
 
  }

0/p:
Code:

/*--------------------------------*- C++ -*----------------------------------*\
  | =========                |                                                |
  | \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
  |  \\    /  O peration    | Version:  1.6-ext                                |
  |  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
  |    \\/    M anipulation  |                                                |
  \*---------------------------------------------------------------------------*/
  FoamFile
  {
      version    2.0;
      format      ascii;
      class      volScalarField;
      location    "0";
      object      p;
  }
  // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  // density normalized pressure!!! p/rho
  dimensions      [0 2 -2 0 0 0 0];
 
  internalField  uniform 0;
 
  boundaryField                               
  {                                                                           
                                                               
      outsideSlider                                                             
      {                                                                       
                  type                      ggi;                                                       
                  value                    uniform 0;
      }                                                                       
                                                               
      insideSlider                                                                 
      {                                                                       
                  type                      ggi;                                                       
                  value                    uniform 0;
      }
 
      outlet                                                                             
      {                                                                       
          type            freestreamPressure;
      }                                                                       
                                                                               
      inlet                                                               
      {                                                                       
          type            freestreamPressure;
      }
 
      sides                                                               
      {                                                                       
          type            freestreamPressure;
      }
 
      ".*STL"                                                                         
      {                                                                       
          type            zeroGradient;
      }
 
  }
 
  // ************************************************************************* //

Code:

  Boundary-file:
  /*--------------------------------*- C++ -*----------------------------------*\
  | =========                |                                                |
  | \\      /  F ield        | OpenFOAM Extend Project: Open source CFD        |
  |  \\    /  O peration    | Version:  1.6-ext                              |
  |  \\  /    A nd          | Web:      www.extend-project.de                |
  |    \\/    M anipulation  |                                                |
  \*---------------------------------------------------------------------------*/
  FoamFile
  {
      version    2.0;
      format      ascii;
      class      polyBoundaryMesh;
      location    "0/polyMesh";
      object      boundary;
  }
  // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
  24
  (
      outlet
      {
          type            patch;
          nFaces          1600;
          startFace      4362923;
      }
      inlet
      {
          type            patch;
          nFaces          1600;
          startFace      4364523;
      }
      sides
      {
          type            patch;
          nFaces          6400;
          startFace      4366123;
      }
      outsideSlider
      {
          type            ggi;
          nFaces          13632;
          startFace      4372523;
          shadowPatch    insideSlider;
          zone            outsideZone;
          bridgeOverlap  false;
      }
      defaultFaces
      {
          type            empty;
          nFaces          0;
          startFace      4386155;
      }
      insideSlider
      {
          type            ggi;
          nFaces          35144;
          startFace      4386155;
          shadowPatch    outsideSlider;
          zone            insideZone;
          bridgeOverlap  false;
      }
      m1_STL
      {
          type            wall;
          nFaces          5880;
          startFace      4421299;
      }
      m2_STL
      {
          type            wall;
          nFaces          5849;
          startFace      4427179;
      }
      m3_STL
      {
          type            wall;
          nFaces          5849;
          startFace      4433028;
      }
      m4_STL
      {
          type            wall;
          nFaces          5772;
          startFace      4438877;
      }
      m5_STL
      {
          type            wall;
          nFaces          5849;
          startFace      4444649;
      }
      m6_STL
      {
          type            wall;
          nFaces          5849;
          startFace      4450498;
      }
      m7_STL
      {
          type            wall;
          nFaces          5772;
          startFace      4456347;
      }
      m8_STL
      {
          type            wall;
          nFaces          5772;
          startFace      4462119;
      }
      m9_STL
      {
          type            wall;
          nFaces          5801;
          startFace      4467891;
      }
      p1_STL
      {
          type            wall;
          nFaces          5880;
          startFace      4473692;
      }
      p2_STL
      {
          type            wall;
          nFaces          5849;
          startFace      4479572;
      }
      p3_STL
      {
          type            wall;
          nFaces          5849;
          startFace      4485421;
      }
      p4_STL
      {
          type            wall;
          nFaces          5772;
          startFace      4491270;
      }
      p5_STL
      {
          type            wall;
          nFaces          5849;
          startFace      4497042;
      }
      p6_STL
      {
          type            wall;
          nFaces          5849;
          startFace      4502891;
      }
      p7_STL
      {
          type            wall;
          nFaces          5772;
          startFace      4508740;
      }
      p8_STL
      {
          type            wall;
          nFaces          5772;
          startFace      4514512;
      }
      p9_STL
      {
          type            wall;
          nFaces          5801;
          startFace      4520284;
      }
  )
 
  // ************************************************************************* //



All times are GMT -4. The time now is 19:08.