[A.Wendy]

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:

  /*--------------------------------*- 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:

/*--------------------------------*- 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;
      }
   
  }
   
  // ************************************************************************* //

  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;
      }
  )
   
  // ************************************************************************* //