CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Time step continuity error

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 12, 2015, 06:05
Default Time step continuity error
  #1
Member
 
Michele Lopez
Join Date: Oct 2015
Location: Bari, IT
Posts: 46
Rep Power: 3
lpz_michele is on a distinguished road
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


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


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


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

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


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


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

// ************************************************************************* //
lpz_michele is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Compile OF 2.3 on Mac OS X .... the patch gschaider OpenFOAM Installation on Windows, Mac and other Unsupported Platforms 225 August 25, 2015 19:43
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 09:31
CGNS lib and Fortran compiler manaliac Main CFD Forum 2 November 29, 2010 07:25
user subroutine error CFDUSER CFX 2 December 9, 2006 07:31
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


All times are GMT -4. The time now is 10:12.