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

interFoam/kOmegaSST tank filling with printStackError/Mules

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 12, 2012, 08:39
Default interFoam/kOmegaSST tank filling with printStackError/Mules
  #1
Member
 
Simon Arne
Join Date: May 2012
Posts: 42
Rep Power: 6
simpomann is on a distinguished road
Hey there,

I am simulating the filling of a tank with interfoam and try to use k-Omega-SST model for turbulence. The mesh has around 400.000 cells, checkMesh says okay. Unluckily I encounter these errors a lot (sometimes after 0.25s, with different transport properties after 0.4s):
Code:
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.524686, Final residual = 0.000681733, No Iterations 2
time step continuity errors : sum local = 1.32975e-10, global = 2.51813e-15, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0920533, Final residual = 7.70789e-08, No Iterations 5
time step continuity errors : sum local = 2.09622e-14, global = -1.45744e-18, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 0.000158889, Final residual = 7.80272e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.00244614, Final residual = 5.717e-10, No Iterations 2
bounding k, min: -0.000200721 max: 3.00901e+20 average: 1.31669e+15
ExecutionTime = 24866.4 s  ClockTime = 24916 s

Courant Number mean: 2.01358e-08 max: 0.712635
Interface Courant Number mean: 8.29938e-09 max: 0.282305
deltaT = 2.2569e-14
Time = 0.2550889169617

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.408204, Final residual = 0.000561119, No Iterations 2
time step continuity errors : sum local = 2.68439e-10, global = -5.36369e-15, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0698315, Final residual = 5.11837e-08, No Iterations 5
time step continuity errors : sum local = 3.27413e-14, global = -7.29092e-18, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 0.00329008, Final residual = 2.10387e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.0124975, Final residual = 7.00019e-09, No Iterations 2
ExecutionTime = 24868.6 s  ClockTime = 24919 s

Courant Number mean: 3.59009e-08 max: 3.49031
Interface Courant Number mean: 1.43023e-08 max: 0.550178
deltaT = 5.17294e-15
Time = 0.255088916961705

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.587722, Final residual = 0.000830391, No Iterations 2
time step continuity errors : sum local = 7.64239e-11, global = 4.21113e-16, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0960292, Final residual = 3.83386e-08, No Iterations 6
time step continuity errors : sum local = 5.13861e-15, global = -9.16091e-19, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 0.000484156, Final residual = 1.42813e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.00071846, Final residual = 1.26128e-09, No Iterations 2
ExecutionTime = 24870.9 s  ClockTime = 24921 s

Courant Number mean: 1.32675e-08 max: 0.702258
Interface Courant Number mean: 5.5049e-09 max: 0.202586
deltaT = 5.76223e-15
Time = 0.25508891696171

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.143881, Final residual = 0.000365565, No Iterations 2
time step continuity errors : sum local = 4.46931e-11, global = -4.75218e-15, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0402696, Final residual = 3.67927e-08, No Iterations 6
time step continuity errors : sum local = 4.85566e-15, global = -1.63317e-18, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 8.34399e-05, Final residual = 9.95434e-12, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00146217, Final residual = 2.85304e-10, No Iterations 2
ExecutionTime = 24873.2 s  ClockTime = 24923 s

Courant Number mean: 1.83689e-08 max: 0.546694
Interface Courant Number mean: 7.7516e-09 max: 0.262968
deltaT = 6.60544e-15
Time = 0.255088916961717

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0613232, Final residual = 0.000199917, No Iterations 2
time step continuity errors : sum local = 3.64953e-11, global = 1.40235e-15, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.00862543, Final residual = 7.72936e-08, No Iterations 5
time step continuity errors : sum local = 1.45349e-14, global = 4.97913e-17, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 0.000305148, Final residual = 9.20208e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.00171966, Final residual = 6.46387e-10, No Iterations 2
bounding k, min: -2.15624e+09 max: 2.68344e+21 average: 1.08957e+16
ExecutionTime = 24875.4 s  ClockTime = 24925 s

Courant Number mean: 2.361e-08 max: 0.507769
Interface Courant Number mean: 1.01369e-08 max: 0.318716
deltaT = 7.64614e-15
Time = 0.255088916961725

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.24767, Final residual = 0.000361503, No Iterations 2
time step continuity errors : sum local = 1.34031e-10, global = 3.46138e-15, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0585152, Final residual = 3.09369e-08, No Iterations 5
time step continuity errors : sum local = 1.38894e-14, global = 8.65978e-17, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 0.00217707, Final residual = 1.56868e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00609363, Final residual = 2.66106e-09, No Iterations 2
ExecutionTime = 24877.6 s  ClockTime = 24928 s

Courant Number mean: 3.25704e-08 max: 1.92022
Interface Courant Number mean: 1.36638e-08 max: 0.455131
deltaT = 3.18553e-15
Time = 0.255088916961728

[3] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2   in "/lib64/libc.so.6"
[3] #3  void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #4   at MULES.C:0
[3] #5  Foam::MULES::explicitSolve(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #6  
[3]  in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/interFoam"
[3] #7  __libc_start_main in "/lib64/libc.so.6"
[3] #8  
[3]  at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116
[c103-06:14726] *** Process received signal ***
[c103-06:14726] Signal: Floating point exception (8)
[c103-06:14726] Signal code:  (-6)
[c103-06:14726] Failing at address: 0x35c600003986
[c103-06:14726] [ 0] /lib64/libc.so.6(+0x34e10) [0x7fd91ae43e10]
[c103-06:14726] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x7fd91ae43d95]
[c103-06:14726] [ 2] /lib64/libc.so.6(+0x34e10) [0x7fd91ae43e10]
[c103-06:14726] [ 3] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geometricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEERKT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMeshEEESK_RKT0_RKT1_ddi+0x1068) [0x7fd91ca71fa8]
[c103-06:14726] [ 4] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(+0x7c007e) [0x7fd91ca6507e]
[c103-06:14726] [ 5] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERS8_dd+0x80) [0x7fd91ca65310]
[c103-06:14726] [ 6] interFoam() [0x42d35e]
[c103-06:14726] [ 7] /lib64/libc.so.6(__libc_start_main+0xed) [0x7fd91ae3023d]
[c103-06:14726] [ 8] interFoam() [0x425119]
[c103-06:14726] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 14726 on node c103-06 exited on signal 8 (Floating point exception).
Timestep decreases and this happens. I checked velocities and pressure at the last working timestep, nothing out of hand here.
I guess I have to supply my fvSchemes, fvSolutions and my initial conditions to trace this. My guess is a connection to the k-Omega-SST model. Probably my initial settings are bad, but I calculated them several times following the wiki here and can't manage to get it work.
Case is currently starting laminar, if this works I will update you asap.

U
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (0 0 0);

boundaryField

{
    walls
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    //Durchströmt
    
    inlet
    {
       type            flowRateInletVelocity;
       flowRate          constant 0.000666667;//incompressible solver: [m³/s]
       value           uniform (0 0 0); //placeholder
    
    //type            fixedValue;            //Alternative, Sascha
        //value           (1.36 -2.56 0);
    }
        
     atmosphaere_CATIASTL
    {
        type            zeroGradient;
    }
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    
    efr_CATIASTL
     {
        type            fixedValue;
        value           uniform (0 0 0);
    }

    entl_CATIASTL
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    
    pistole_vorn_CATIASTL
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    
    
    pistole_hinten_CATIASTL
     {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    
    innenteile_CATIASTL
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
       
    flv_CATIASTL
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    
   einsatz_efr_CATIASTL
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    
    verschlusskopf_CATIASTL
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    
}

// ************************************************************************* //
p_rgh

Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField

{
    walls
    {
        type            buoyantPressure;
        value           uniform 0;
    }
    //Durchströmt
    
    inlet
    {
        type            buoyantPressure;
        value           uniform 0;
    }
    
    atmosphaere_CATIASTL
   {
        type            totalPressure;
        p0              uniform 0;
        U               U;
        phi             phi;
        rho             rho;
        psi             none;
        gamma           1;
        value           uniform 0;    
    }
    
    //Alles was umströmt wird

    tank_CATIASTL
     {
        type            buoyantPressure;
        value           uniform 0;
    }
    
    efr_CATIASTL
     {
        type            buoyantPressure;
        value           uniform 0;
    }

    entl_CATIASTL
     {
        type            buoyantPressure;
        value           uniform 0;
    }
    
    pistole_vorn_CATIASTL
      {
        type            buoyantPressure;
        value           uniform 0;
    }
    
   // sd_CATIASTL
    //{
    //    type            buoyantPressure;
    //    value           uniform 0;
    //}
    
    pistole_hinten_CATIASTL
    {
        type            buoyantPressure;
        value           uniform 0;
    }
    
    innenteile_CATIASTL
    {
        type            buoyantPressure;
        value           uniform 0;
    }
    
    
    flv_CATIASTL
    {
        type            buoyantPressure;
        value           uniform 0;
    }
    
       einsatz_efr_CATIASTL
    {
        type            buoyantPressure;
        value           uniform 0;
    }
    
    verschlusskopf_CATIASTL
    {
        type            buoyantPressure;
        value           uniform 0;
    }
    
}

// ************************************************************************* //
k
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0.0286;

boundaryField

{
    walls
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    //Durchströmt
    
    inlet
    {
        type            turbulentIntensityKineticEnergyInlet;
        intensity       0.041;       // 5% turbulent intensity
        value           $internalField;
    }
    
        atmosphaere_CATIASTL
    {
        type            inletOutlet;
        inletValue      $internalField;
    }
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    
    efr_CATIASTL
     {
        type            kqRWallFunction;
        value           $internalField;
    }

    entl_CATIASTL
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    
    pistole_vorn_CATIASTL
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    
    sd_CATIASTL
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    
    pistole_hinten_CATIASTL
     {
        type            kqRWallFunction;
        value           $internalField;
    }
    
    innenteile_CATIASTL
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    
    
    flv_CATIASTL
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    
       einsatz_efr_CATIASTL
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    
    verschlusskopf_CATIASTL
    {
        type            kqRWallFunction;
        value           $internalField;
    }
}


// ************************************************************************* //
Omega
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


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

internalField   uniform 1643.71;

boundaryField

{
    walls
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    //Durchströmt
    
    inlet
    {
        type            fixedValue;
        value           uniform 280.78;
    }
       
        atmosphaere_CATIASTL
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }
    
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    
    efr_CATIASTL
     {
        type            omegaWallFunction;
        value           $internalField;
    }

    entl_CATIASTL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    
    pistole_vorn_CATIASTL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    
    sd_CATIASTL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    
    pistole_hinten_CATIASTL
     {
        type            omegaWallFunction;
        value           $internalField;
    }
    
    innenteile_CATIASTL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    
    
    flv_CATIASTL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    
    einsatz_efr_CATIASTL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    
    verschlusskopf_CATIASTL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
}


// ************************************************************************* //
Nut
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0.00010198;

boundaryField

{
    walls
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    //Durchströmt
    
    inlet
    {
        type            calculated;
        value           uniform 0.00010198;
    }
        
        atmosphaere_CATIASTL
    {
        type            calculated;
        value           $internalField;
    }
    
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    
    efr_CATIASTL
     {
        type            nutkWallFunction;
        value           $internalField;
    }

    entl_CATIASTL
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    
    pistole_vorn_CATIASTL
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    
    sd_CATIASTL
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    
    pistole_hinten_CATIASTL
     {
        type            nutkWallFunction;
        value           $internalField;
    }
    
    innenteile_CATIASTL
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    

    
    flv_CATIASTL
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    
    einsatz_efr_CATIASTL
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    
    verschlusskopf_CATIASTL
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    
}
I guess inlet is self explaining. The only outlet is "atmosphaere". All the other stuff are pipes etc, so basically walls.
If somebody spots a mistake in one of the files please tell me! I am super lost with this.

All help is highly appreciated! Big thanks and best regards!
I hope somebody with more experience can help me to fix this, unluckily it's a bit urgent for my thesis.

Greetings!
Simon
simpomann is offline   Reply With Quote

Old   December 12, 2012, 08:40
Default
  #2
Member
 
Simon Arne
Join Date: May 2012
Posts: 42
Rep Power: 6
simpomann is on a distinguished road
fvSchemes
Code:
ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;            //Alternative:default cellMDLimited Gauss linear 1
}

divSchemes
{ 
    div(rho*phi,U)  Gauss limitedLinearV 1;
    div(phi,alpha)  Gauss vanLeer; 
    div(phirb,alpha) Gauss interfaceCompression;
    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;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p;
    pcorr;
    alpha1;
    p_rgh;
}
fvSolutions
Code:
solvers
{
    p_rgh                //Lösung der Druckgleichung
    {
        solver          GAMG;
        tolerance       1e-07;
        relTol          0.01;
        smoother        DIC;
        nPreSweeps      0;
        nPostSweeps     2;
        nFinestSweeps   2;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 10;
        agglomerator    faceAreaPair;
        mergeLevels     1;
    }

    p_rghFinal                //finaler Loop
    {
        solver          PCG;
        preconditioner
        {
            preconditioner  GAMG;
            tolerance       1e-07;
            relTol          0;
            nVcycles        2;
            smoother        DICGaussSeidel;
            nPreSweeps      2;
            nPostSweeps     2;
            nFinestSweeps   2;
            cacheAgglomeration true;
            nCellsInCoarsestLevel 10;
            agglomerator    faceAreaPair;
            mergeLevels     1;
        }
        tolerance       1e-07;
        relTol          0;
        maxIter         10;
    }

    pcorr                //vor dem Loop gelöste Druckvergleichung für Startwerte des Feldes innerhalb des Iterationsschritts
    {
        solver          PCG;
        preconditioner
        {
            preconditioner  GAMG;
            tolerance       1e-05;
            relTol          0;
            smoother        DICGaussSeidel;
            nPreSweeps      0;
            nPostSweeps     2;
            nFinestSweeps   2;
            cacheAgglomeration false;
            nCellsInCoarsestLevel 10;
            agglomerator    faceAreaPair;
            mergeLevels     1;
        }
        tolerance       1e-05;
        relTol          0;
        maxIter         10;
    }


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

    "(U|k|omega)Final"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-07;
        relTol          0;
    }
}

PIMPLE
{
    momentumPredictor no;
    nCorrectors     2;
    nNonOrthogonalCorrectors 0;
    nAlphaCorr      1;
    nAlphaSubCycles 2;
    cAlpha          1;  // Alternative 0
}

And alpha (dont expect problems here)
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    walls
    {
        type            zeroGradient;
    }

    //Durchströmt
    
    inlet
    {
        type            fixedValue;
        value           uniform 1;
    }
        
        atmosphaere_CATIASTL
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            zeroGradient;
    }
    
    efr_CATIASTL
    {
        type            zeroGradient;
    }

    entl_CATIASTL
    {
        type            zeroGradient;
    }
    
    pistole_vorn_CATIASTL
    {
        type            zeroGradient;
    }
    
    sd_CATIASTL
    {
        type            zeroGradient;
    }
    
    pistole_hinten_CATIASTL
    {
        type            zeroGradient;
    }
    
    innenteile_CATIASTL
    {
        type            zeroGradient;
    }
    

    
    flv_CATIASTL
    {
        type            zeroGradient;
    }
    
       einsatz_efr_CATIASTL
    {
        type            zeroGradient;
    }
    
    verschlusskopf_CATIASTL
    {
        type            zeroGradient;
    }

}

// ************************************************************************* //
simpomann is offline   Reply With Quote

Old   December 12, 2012, 13:11
Default
  #3
Member
 
Simon Arne
Join Date: May 2012
Posts: 42
Rep Power: 6
simpomann is on a distinguished road
Hey,

I checked a lamniar case. Same error but at a different timestep
Code:
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62932e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62932e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.065834, Final residual = 0.000307154, No Iterations 2
time step continuity errors : sum local = 2.85125e-12, global = -6.3837e-15, cumulative = -3.90402e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.02062, Final residual = 3.97336e-08, No Iterations 7
time step continuity errors : sum local = 3.9867e-16, global = -5.09275e-19, cumulative = -3.90402e-05
ExecutionTime = 15548.4 s  ClockTime = 15612 s

Courant Number mean: 2.84517e-08 max: 0.867998
Interface Courant Number mean: 5.5157e-09 max: 0.0944208
deltaT = 4.14157e-15
Time = 0.202400597327766

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0662251, Final residual = 0.000298758, No Iterations 2
time step continuity errors : sum local = 2.77135e-12, global = -4.86472e-15, cumulative = -3.90402e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0204396, Final residual = 3.2508e-08, No Iterations 7
time step continuity errors : sum local = 3.2613e-16, global = -4.25822e-19, cumulative = -3.90402e-05
ExecutionTime = 15550.6 s  ClockTime = 15614 s

Courant Number mean: 2.84501e-08 max: 0.867991
Interface Courant Number mean: 5.51484e-09 max: 0.0943805
deltaT = 3.81716e-15
Time = 0.202400597327769

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0657749, Final residual = 0.00031065, No Iterations 2
time step continuity errors : sum local = 2.88026e-12, global = -6.24636e-15, cumulative = -3.90402e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0206025, Final residual = 3.60495e-08, No Iterations 7
time step continuity errors : sum local = 3.61229e-16, global = -3.50359e-19, cumulative = -3.90402e-05
ExecutionTime = 15552.7 s  ClockTime = 15616 s

Courant Number mean: 2.84444e-08 max: 0.867923
Interface Courant Number mean: 5.6614e-09 max: 0.094348
deltaT = 3.51843e-15
Time = 0.202400597327773

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0661369, Final residual = 0.000300055, No Iterations 2
time step continuity errors : sum local = 2.77963e-12, global = -4.82001e-15, cumulative = -3.90402e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0204266, Final residual = 3.17341e-08, No Iterations 7
time step continuity errors : sum local = 3.17909e-16, global = -3.68428e-19, cumulative = -3.90402e-05
ExecutionTime = 15554.9 s  ClockTime = 15619 s

Courant Number mean: 2.8439e-08 max: 0.867925
Interface Courant Number mean: 5.317e-09 max: 0.094304
deltaT = 3.24307e-15
Time = 0.202400597327776

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0657622, Final residual = 0.00031015, No Iterations 2
time step continuity errors : sum local = 2.87193e-12, global = -6.3359e-15, cumulative = -3.90402e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0205844, Final residual = 3.79572e-08, No Iterations 7
time step continuity errors : sum local = 3.79847e-16, global = -3.64393e-19, cumulative = -3.90402e-05
ExecutionTime = 15557.1 s  ClockTime = 15621 s

Courant Number mean: 2.84311e-08 max: 0.867857
Interface Courant Number mean: 5.31414e-09 max: 0.0942674
deltaT = 2.9895e-15
Time = 0.202400597327779

[3] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2   in "/lib64/libc.so.6"
[3] #3  void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #4   at MULES.C:0
[3] #5  Foam::MULES::explicitSolve(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #6  
[3]  in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/interFoam"
[3] #7  __libc_start_main in "/lib64/libc.so.6"
[3] #8  
[3]  at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116
[c103-06:28794] *** Process received signal ***
[c103-06:28794] Signal: Floating point exception (8)
[c103-06:28794] Signal code:  (-6)
[c103-06:28794] Failing at address: 0x35c60000707a
[c103-06:28794] [ 0] /lib64/libc.so.6(+0x34e10) [0x7f3d0d63ae10]
[c103-06:28794] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x7f3d0d63ad95]
[c103-06:28794] [ 2] /lib64/libc.so.6(+0x34e10) [0x7f3d0d63ae10]
[c103-06:28794] [ 3] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geometricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEERKT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMeshEEESK_RKT0_RKT1_ddi+0x1068) [0x7f3d0f268fa8]
[c103-06:28794] [ 4] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(+0x7c007e) [0x7f3d0f25c07e]
[c103-06:28794] [ 5] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERS8_dd+0x80) [0x7f3d0f25c310]
[c103-06:28794] [ 6] interFoam() [0x42d35e]
[c103-06:28794] [ 7] /lib64/libc.so.6(__libc_start_main+0xed) [0x7f3d0d62723d]
[c103-06:28794] [ 8] interFoam() [0x425119]
[c103-06:28794] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 28794 on node c103-06 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
As far as I understand, there is an error during the alpha calculation. But what can cause it? If the turbulence model is allright, checkMesh says okay.... where can I start searching my answers?
Could it be this part from my fvSolution?
PIMPLE { momentumPredictor no; nCorrectors 2; nNonOrthogonalCorrectors 0; nAlphaCorr 1; nAlphaSubCycles 2; cAlpha 1; // Alternative 0 }

Last edited by simpomann; December 12, 2012 at 13:43.
simpomann is offline   Reply With Quote

Old   February 17, 2014, 18:06
Default
  #4
New Member
 
Matej Muller
Join Date: Oct 2011
Location: Slovenia
Posts: 22
Rep Power: 6
matejmuller is on a distinguished road
Hey Simon!

Have you perhaps solved this issue? I'm also trying to run interFoam with the kOmega model and also don't know which solvers and schemes to use, since the tutorial of kOmegaSST exists only for the LTSInterFoam. Have you figured out the right solvers and schemes?

thanks for the help,

Matej
matejmuller 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
tank filling (water/air), I think I messed up the BCs simpomann OpenFOAM Running, Solving & CFD 1 September 8, 2012 07:09
BCs to the outlet patch (filling a simple tank) Ternox OpenFOAM Pre-Processing 3 February 17, 2012 11:21
Filling Tank with Water leff CFX 5 October 1, 2011 09:30
free-surface liquid filling a fuel tank lee CFX 2 June 4, 2008 21:01
how to set boundary condition for a tank filling w namrata FLUENT 6 May 2, 2007 10:23


All times are GMT -4. The time now is 21:35.