[oldirtybrindley]

Not sure what's going on. Can't get simpleFoam to run

Time = 1

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.081410493, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.061596729, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.091684187, No Iterations 3
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0059667553, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.40650934, Final residual = 0.0012374181, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.036894522, Final residual = 0.00027957842, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0062564286, Final residual = 4.8770264e-05, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0013606213, Final residual = 1.2695641e-05, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00035198339, Final residual = 3.117196e-06, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00010876533, Final residual = 1.0124916e-06, No Iterations 3
GAMG:  Solving for p, Initial residual = 3.7107364e-05, Final residual = 3.2700493e-07, No Iterations 3
GAMG:  Solving for p, Initial residual = 1.4011258e-05, Final residual = 1.2603484e-07, No Iterations 3
GAMG:  Solving for p, Initial residual = 5.5796704e-06, Final residual = 3.7850784e-08, No Iterations 6
GAMG:  Solving for p, Initial residual = 2.3340101e-06, Final residual = 6.442118e+119, No Iterations 1000
time step continuity errors : sum local = 1.2632859e+123, global = 4.4873194e+101, cumulative = 4.4873194e+101
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:?
#9  Foam::incompressible::RASModels::kOmegaSST::correct() at ??:?
#10  
 at ??:?
#11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12  
 at ??:?
Floating point exception (core dumped)

[oldirtybrindley]

[

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.0-f5222ca19ce6
Exec   : checkMesh
Date   : Mar 07 2014
Time   : 20:04:09
Host   : "AVA-411296"
PID    : 26818
Case   : /home/kathyb/OpenFOAM/kathyb-2.3.0/run/SplitTipTesting/AOA/simplesolvefixed(donotuse)
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           30927
    faces:            72559
    internal faces:   60772
    cells:            21783
    faces per cell:   6.1208741
    boundary patches: 5
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     15749
    prisms:        1429
    wedges:        0
    pyramids:      0
    tet wedges:    128
    tetrahedra:    2
    polyhedra:     4475
    Breakdown of polyhedra by number of faces:
        faces   number of cells
            4   1186
            5   814
            6   740
            7   218
            8   90
            9   678
           10   86
           11   47
           12   401
           13   3
           14   6
           15   150
           18   56

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
                   Patch    Faces   Points                  Surface topology
                   inlet        0        0                        ok (empty)
                  outlet        0        0                        ok (empty)
                   walls       40       67  ok (non-closed singly connected)
               frontBack        0        0                        ok (empty)
                    Foil    11747    13554  ok (non-closed singly connected)

Checking geometry...
    Overall domain bounding box (-0.020220271 -0.010812593 -0.00023150215) (0.020260984 0.010818352 0.15183358)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-5.8212676e-17 4.5742309e-17 1.775036e-17) OK.
    Max cell openness = 2.9738872e-16 OK.
    Max aspect ratio = 21.38308 OK.
    Minimum face area = 6.6009305e-09. Maximum face area = 1.3716404e-05.  Face area magnitudes OK.
    Min volume = 1.1250083e-12. Max volume = 2.2923998e-08.  Total volume = 2.8560307e-05.  Cell volumes OK.
    Mesh non-orthogonality Max: 62.160689 average: 13.151348
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 2.6576963 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

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

solvers
{
    p
    {
        solver           GAMG;
        tolerance        0.02;
        relTol           0;
        smoother         GaussSeidel;
        nPreSweeps       0;
        nPostSweeps      2;
        cacheAgglomeration on;
        agglomerator     faceAreaPair;
        nCellsInCoarsestLevel 10;
        mergeLevels      1;
    }


    U
    {
        solver           PBiCG;
	preconditioner	 DILU;
        tolerance        1e-6;
        relTol		0.1;
    }

    k
    {
        solver           PBiCG;
        preconditioner	 DILU;
        tolerance        1e-6;
        relTol		0.1;
    }

    omega
    {
                solver           PBiCG;
preconditioner	 DILU;
        tolerance        1e-6;
        relTol		0.1;
    }

	cellDisplacement
	{
	solver	GAMG;
	tolerance	1e-5;
	relTol		0;
	smoother	GaussSeridel;
	cacheAgglomeration	true;
	nCellsInCoursestLevel	10;
	agglomerator	faceAreaPair;
	mergeLevels	1;

	}


}

SIMPLE
{
    nNonOrthogonalCorrectors 0;

    pReffCell 0;
    pRefValue 0;
    residualControl
    {
        p 1e-4;
	U 1e-4;
    }
}

potentialFlow
{
    nNonOrthogonalCorrectors 10;
}

relaxationFactors
{
    fields
    {
        p               0.3;
    }
    equations
    {
        U               0.7;
        k               0.7;
        omega           0.7;
    }
}

cache
{
    grad(U);
}

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

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

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
    grad(U)	    Gauss linear;
    grad(k)    	    cellLimited Gauss linear 1;
    grad(omega)    	    cellLimited Gauss linear 1;
}

divSchemes
{
    default         none;
    div(phi,U)      bounded Gauss upwind;;
    div(phi,k)      bounded Gauss upwind;
    div(phi,omega)  bounded Gauss upwind;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited corrected 0.5;
    laplacian(DomegaEff,omega)		Gauss linear limited 0.33;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.33;
}

fluxRequired
{
    default         no;
    p;
}

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

[Lieven]

Hi olddirtybrindly,

Do you get the same crash when you run the case with pisofoam?

Cheers,

L