[lpz_michele]

Hi everyone!

This is my first time writing on this CFD forum. I'm trying to simulate flowfield around a VAWT, Darrieus type. But something is going wrong with time step continuity error. On the first time step it reaches a global value of e+36.

I'm using AMI to simulate the rotation of the turbine. I've exported the mesh from Pointwise than merged 2 different part with mergeMeshes utility.

Note that:

• pala*, refers to the 3 airfoil of the rotor
• cubo*, refers to the external volume
• All the other b.c. are cyclicAMI

These are all the files i'm using in system

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     pimpleDyMFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         0.7853976;

deltaT          0.0000109083;

writeControl    adjustableRunTime;

writeInterval   0.0000109083;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  no;

maxCo           1;

// ************************************************************************* //

functions
(
forces
{
type forces;
functionObjectLibs ("libforces.so");
outputControl timeStep;
outputInterval 1;
patches (pala1 pala1sopra pala1sotto pala2 pala2sopra pala2sotto pala3 pala3sopra pala3sotto);
// pname p;
// Uname U;
rhoName rhoInf;
log true;
rhoInf 1.205;
CofR (0 0 0);
}
);

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default                     leastSquares;
    grad(p)                     leastSquares;
    grad(U)                    leastSquares;

}

divSchemes
{
    default                     none;

    div(phi,U)                  Gauss linearUpwind grad(U);
    div(phi,k)                  Gauss upwind;
    div(phi,omega)             Gauss upwind;
    div((nuEff*dev(T(grad(U)))))     Gauss linear;
}

laplacianSchemes
{
    default                     Gauss linear corrected;
}

interpolationSchemes
{
    default                     linear;
    interpolate(HbyA)             linear;

}

snGradSchemes
{
    default                     corrected;
}

fluxRequired
{
    default         no;
    p               ;
    pcorr           ;
}


// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    pcorr
    {
        solver GAMG;
    smoother GaussSeidel;
    nPreSweeps 0;
    nPostSweeps 2;
    cacheAgglomeration off;
    agglomerator faceAreaPair;
    nCellsInCoarsestLevel 10;
    mergeLevels 1;

    tolerance 0.02;
    relTol 0;
    }

    p
    {
    $pcorr;
    tolerance 1e-07;
    relTol 0.01;
    }  
    
    pFinal
    {
        $p;
    tolerance     1e-07;
        relTol          0;
    }


    U
    {
    solver        smoothSolver;
    smoother    GaussSeidel;
    tolerance    1e-8;
    relTol        0.1;

    }
    
    UFinal
    {
    solver        smoothSolver;
    smoother    GaussSeidel;
    tolerance    1e-8;
    relTol        0.1;

    }
    
    "(k|omega)"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        tolerance       1e-7;
        relTol          0.1;
    }

    "(k|omega)Final"
    {
        $U;
    tolerance     1e-7;
        relTol          0;
    }

    cellMotionUx
    {
    solver        PCG;
    preconditioner     DIC;
    tolerance    1e-08;
    relTol        0;
    }

}

PIMPLE
{
    correctPhi        no;
    nOuterCorrectors    100;
    nCorrectors         2;
    nNonOrthogonalCorrectors 0;
    
    pRefCell         0;
    pRefValue        0;

residualControl
  {
    "(p|U|k|omega)"
     {
      tolerance 1e-2;
      relTol 0;
      }
  }
}

relaxationFactors
{
    fields
    {
    }
    equations
    {
    "p.*"        0.3;
    "U.*"        0.7;
    "k.*"        0.8;
    "omega.*"        0.8;
    }
}


// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      decomposeParDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

numberOfSubdomains     4;

preservePatches     (fisso fissosopra fissosotto volumelato volumesopra volumesotto);

method              scotch;

distributed         no;

roots             ();


// ************************************************************************* //

And p, U, k, omega in 0

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [ 0 2 -2 0 0 0 0 ];

internalField   uniform 0.375;

boundaryField
{
   
   
    inlet
    {
        type                turbulentIntensityKineticEnergyInlet;
    intensity         0.05;
    value             $internalField;
    }
    
    
    outlet
    {
         type                inletOutlet;
        inletValue          $internalField;
        value               $internalField;      
    }
  
  
    pala1
    {
        type             kqRWallFunction;
        value               $internalField;
    }
    
    
    pala1sopra
    {
        type              kqRWallFunction;
        value               $internalField;      
    }
    
    
    pala1sotto
    {
    type             kqRWallFunction;
        value                  $internalField;
    }
    
    
    pala2
    {
        type               kqRWallFunction;
        value               $internalField;   
    }
    
    
    pala2sopra
    {
        
       type            kqRWallFunction;
       value                   $internalField;
    }
    
    
    pala2sotto
    {
       type               kqRWallFunction;
        value               $internalField;
    }
    
    
    pala3
    {
        type            kqRWallFunction;
        value               $internalField;
    }
    
    
    pala3sopra
    {
        type              kqRWallFunction;
        value               $internalField;    
    }
    
    
    pala3sotto
    {
     type             kqRWallFunction;
     value                   $internalField;
    }
    
      volumelato
    {
        type         cyclicAMI;
    value         $internalField;
     }
     
     
    volumesopra
    {
        type         cyclicAMI;
    value         $internalField;
    }
    
    
    volumesotto
    {
        type         cyclicAMI;
    value         $internalField;
    }
    
    cubodx
    {
        type         zeroGradient;
        
    }
    
    
    cubosopra
    {
        type         zeroGradient;
        
    }
    
    
    cubosotto
    {
        type         zeroGradient;
       
    }
    
    
    cubosx
    {
        type         zeroGradient;
        
    }
    
    
    fisso
    {
        type         cyclicAMI;
    value         $internalField;
    }
    
    
    fissosopra
    {
        type         cyclicAMI;
        value         $internalField;
    }
    
    
    fissosotto
    {
        type         cyclicAMI;
        value         $internalField;
    }
   
}


// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



dimensions      [0 0 -1 0 0 0 0];

internalField   uniform 1118.034;

boundaryField
{
     cubodx
    {
        type         zeroGradient;
        
    }
    
    
    cubosopra
    {
        type         zeroGradient;
        
    }
    
    
    cubosotto
    {
        type         zeroGradient;
       
    }
    
    
    cubosx
    {
        type         zeroGradient;
        
    }
    
    inlet
    {
        type            turbulentMixingLengthFrequencyInlet;
    mixingLength    0.001;
        value           $internalField;  
    }
    
    
    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;      
    }
    
    
    fisso
    {
        type         cyclicAMI;
    value         $internalField;
    }
    
    
    fissosopra
    {
        type         cyclicAMI;
        value         $internalField;
    }
    
    
    fissosotto
    {
        type         cyclicAMI;
        value         $internalField;
    }
    
     volumelato
    {
        type         cyclicAMI;
    value         $internalField;
     }
     
     
    volumesopra
    {
        type         cyclicAMI;
    value         $internalField;
    }
    
    
    volumesotto
    {
        type         cyclicAMI;
    value         $internalField;
    }
    
       pala1
    {
        type         omegaWallFunction;
    value         $internalField;
    }
    
    
    pala1sopra
    {
        type         omegaWallFunction;
    value         $internalField;
    }
    
    
    pala1sotto
    {
    type         omegaWallFunction;
    value         $internalField;
    }
    
    
    pala2
    {
        type         omegaWallFunction;
    value         $internalField;
    }
    
    
    pala2sopra
    {
        
       type        omegaWallFunction;
       value         $internalField;
    }
    
    
    pala2sotto
    {
       type          omegaWallFunction;
       value         $internalField;
    }
    
    
    pala3
    {
        type        omegaWallFunction;
    value         $internalField;
    }
    
    
    pala3sopra
    {
        type         omegaWallFunction;
    value         $internalField;
    }
    
    
    pala3sotto
    {
     type         omegaWallFunction;
     value         $internalField;
    }
    
}


// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     pimpleDyMFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         0.7853976;

deltaT          0.0000109083;

writeControl    adjustableRunTime;

writeInterval   0.0000109083;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  no;

maxCo           1;

// ************************************************************************* //

functions
(
forces
{
type forces;
functionObjectLibs ("libforces.so");
outputControl timeStep;
outputInterval 1;
patches (pala1 pala1sopra pala1sotto pala2 pala2sopra pala2sotto pala3 pala3sopra pala3sotto);
// pname p;
// Uname U;
rhoName rhoInf;
log true;
rhoInf 1.205;
CofR (0 0 0);
}
);

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    cubodx
    {
        type freestream;
    freestreamValue $internalField;
    }
    
    
    cubosopra
    {
        type   freestream;
    freestreamValue $internalField;      
    }
    
    
    
    cubosotto
    {
        type freestream;
    freestreamValue $internalField;       
    }
    
    
    
    cubosx
    {
        type freestream;
    freestreamValue $internalField;        
    }
    
    
    
    fisso
    {
        type cyclicAMI;
    value $internalField;
    }
    
    
    
    fissosopra
    {
        type cyclicAMI;
        value $internalField;
    }
    
    
    
    fissosotto
    {
        type cyclicAMI;
        value $internalField;
    }
    
    
    
    inlet
    {
        type fixedValue;
    value uniform (10.0 0 0);       
    }
    
    
    
    outlet
    {
        type zeroGradient;
        
    }
  
    pala1
    {
        type fixedValue;
    value $internalField;
    }
    
    
    pala1sopra
    {
        type fixedValue;
    value $internalField;
        
    }
    
    
    pala1sotto
    {
    type fixedValue;
    value $internalField; 
    }
    
    
    
    pala2
    {
        type fixedValue;
    value $internalField;
      
    }
    
    
    
    pala2sopra
    {
        
       type    fixedValue;
       value     $internalField;
    }
    
    
    
    pala2sotto
    {
       type  fixedValue;
       value $internalField; 
    }
    
    
    
    pala3
    {
        type    fixedValue ;
    value $internalField;
    }
    
    
    
    pala3sopra
    {
        type fixedValue;
    value $internalField;
      
    }
    
    
    
    pala3sotto
    {
     type fixedValue;
     value $internalField;
    }
    
    
    volumelato
    {
        type cyclicAMI;
    value $internalField;
     }
     
     
    volumesopra
    {
        type cyclicAMI;
    value $internalField;
    }
    
    
    volumesotto
    {
        type cyclicAMI;
    value $internalField;
    }
    
   
}

// ************************************************************************* //