[mihaipruna]

Please help!

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : rhoSimplecFoam
Date   : Aug 08 2012
Time   : 10:49:13
Host   : "cadnexus-cae"
PID    : 29421
Case   : /home/cadnexus/OpenFOAM/cadnexus-2.1.1/run/sduct_div/sduct_100_inlet_vel_baseline
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 mesh for time = 0


SIMPLE: convergence criteria
    field p     tolerance 1e-05
    field U     tolerance 1e-05
    field omega     tolerance 1e-05

Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<sutherlandTransport<specieThermo<hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1         0.85034;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.85616;
    Prt             1;
    gamma1          0.5532;
    gamma2          0.4403;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    c1              10;
}


Starting time loop

Reading surface description:
    zNormal

--> FOAM Warning : 
    From function Foam::List<Foam::tetIndices> Foam::polyMeshTetDecomposition::faceTetIndices(const polyMesh&, label, label)
    in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 565
    No base point for face 160047, 4(62180 36360 67554 53019), produces a valid tet decomposition.
--> FOAM Warning : 
    From function Foam::List<Foam::tetIndices> Foam::polyMeshTetDecomposition::faceTetIndices(const polyMesh&, label, label)
    in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 565
    No base point for face 160047, 4(62180 36360 67554 53019), produces a valid tet decomposition.
--> FOAM Warning : 
    From function Foam::List<Foam::tetIndices> Foam::polyMeshTetDecomposition::faceTetIndices(const polyMesh&, label, label)
    in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 565
    No base point for face 169137, 4(32651 49488 69117 53693), produces a valid tet decomposition.
--> FOAM Warning : 
    From function Foam::List<Foam::tetIndices> Foam::polyMeshTetDecomposition::faceTetIndices(const polyMesh&, label, label)
    in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 565
    No base point for face 169137, 4(32651 49488 69117 53693), produces a valid tet decomposition.
--> FOAM Warning : 
    From function Foam::List<Foam::tetIndices> Foam::polyMeshTetDecomposition::faceTetIndices(const polyMesh&, label, label)
    in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 565
    No base point for face 174804, 4(24053 54113 70087 54114), produces a valid tet decomposition.
--> FOAM Warning : 
    From function Foam::List<Foam::tetIndices> Foam::polyMeshTetDecomposition::faceTetIndices(const polyMesh&, label, label)
    in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 565
    No base point for face 174804, 4(24053 54113 70087 54114), produces a valid tet decomposition.
--> FOAM Warning : 
    From function Foam::List<Foam::tetIndices> Foam::polyMeshTetDecomposition::faceTetIndices(const polyMesh&, label, label)
    in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 565
    No base point for face 179330, 4(24272 54485 70817 54484), produces a valid tet decomposition.
--> FOAM Warning : 
    From function Foam::List<Foam::tetIndices> Foam::polyMeshTetDecomposition::faceTetIndices(const polyMesh&, label, label)
    in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 565
    No base point for face 179330, 4(24272 54485 70817 54484), produces a valid tet decomposition.
Time = 1

GAMG:  Solving for Ux, Initial residual = 1, Final residual = 0.0181092, No Iterations 2
GAMG:  Solving for Uy, Initial residual = 1, Final residual = 0.0169639, No Iterations 2
GAMG:  Solving for Uz, Initial residual = 1, Final residual = 0.026404, No Iterations 2
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.096998, No Iterations 19
time step continuity errors : sum local = 474.441, global = -56.4332, cumulative = -56.4332
rho max/min : 1.28594 0.5
GAMG:  Solving for h, Initial residual = 1, Final residual = 0.0497605, No Iterations 2
GAMG:  Solving for omega, Initial residual = 0.123066, Final residual = 0.00214602, No Iterations 2
GAMG:  Solving for k, Initial residual = 1, Final residual = 0.0881728, No Iterations 1
ExecutionTime = 177.11 s  ClockTime = 180 s

Time = 2

GAMG:  Solving for Ux, Initial residual = 0.511533, Final residual = 0.0391677, No Iterations 1
GAMG:  Solving for Uy, Initial residual = 0.541348, Final residual = 0.0441529, No Iterations 1
GAMG:  Solving for Uz, Initial residual = 0.588956, Final residual = 0.0252274, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0901925, Final residual = 0.00839962, No Iterations 6
time step continuity errors : sum local = 323.523, global = -51.6731, cumulative = -108.106
rho max/min : 1.28886 0.5
GAMG:  Solving for h, Initial residual = 0.958599, Final residual = 0.0217326, No Iterations 2
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::hPsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#4  Foam::hPsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#5  
 in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/rhoSimplecFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  
 in "/opt/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/rhoSimplecFoam"
Floating point exception (core dumped)
cadnexus@cadnexus-cae:~/OpenFOAM/cadnexus-2.1.1/run/sduct_div/sduct_100_inlet_vel_baseline$

[mihaipruna]

checkmesh results:

Exec   : checkMesh
Date   : Aug 08 2012
Time   : 11:02:04
Host   : "cadnexus-cae"
PID    : 29636
Case   : /home/cadnexus/OpenFOAM/cadnexus-2.1.1/run/sduct_div/sduct_100_inlet_vel_baseline
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

fileName::stripInvalid() called for invalid fileName constant(copy).tar.gz
    For debug level (= 2) > 1 this is considered fatal
Aborted (core dumped)
cadnexus@cadnexus-cae:~/OpenFOAM/cadnexus-2.1.1/run/sduct_div/sduct_100_inlet_vel_baseline$ checkMesh
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : checkMesh
Date   : Aug 08 2012
Time   : 11:03:06
Host   : "cadnexus-cae"
PID    : 29643
Case   : /home/cadnexus/OpenFOAM/cadnexus-2.1.1/run/sduct_div/sduct_100_inlet_vel_baseline
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:           6376565
    faces:            16789124
    internal faces:   15762710
    cells:            5249544
    boundary patches: 12
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     4568654
    prisms:        143215
    wedges:        0
    pyramids:      0
    tet wedges:    1
    tetrahedra:    0
    polyhedra:     537674

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                  
    topYmax             0        0        ok (empty)                        
    bottomYmin          0        0        ok (empty)                        
    inletXmin           0        0        ok (empty)                        
    outletXmax          0        0        ok (empty)                        
    rightZmax           0        0        ok (empty)                        
    leftZmin            36672    41286    ok (non-closed singly connected)  
    SDuctOutlet         29668    30254    ok (non-closed singly connected)  
    vol1face2           605591   690845   ok (non-closed singly connected)  
    vol1face3           0        0        ok (empty)                        
    vol1face4           0        0        ok (empty)                        
    vol1face5           313384   317020   ok (non-closed singly connected)  
    SDuctInlet          41099    41814    ok (non-closed singly connected)  

Checking geometry...
    Overall domain bounding box (-0.100006 -0.0814332 0) (2.05019 0.911797 0.353048)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-2.48235e-15 2.82511e-16 -1.21331e-14) OK.
    Max cell openness = 3.4866e-16 OK.
    Max aspect ratio = 8.79826 OK.
    Minumum face area = 3.80834e-08. Maximum face area = 0.0025653.  Face area magnitudes OK.
    Min volume = 7.99893e-12. Max volume = 0.000127949.  Total volume = 0.257158.  Cell volumes OK.
    Mesh non-orthogonality Max: 57.0097 average: 9.26743
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 3.91413 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

[mihaipruna]

Image of a solved case (using simpleFoam)

[Tobi]

you still have that problem? If yes , I ll have a look at your case.

Tobi

[mihaipruna]

as of now, yes, no matter what I do I am experiencing crashes rather soon in the solution.
Thank you.

[Tobi]

doublepost / sorry

[Tobi]

Your uploaded files are working.
I dont know what you are wanna do exaclty couse your files are in caos
And I think you are not using OF in a terminal arent you? Couse all headers are away.

Well i can solve that block with rhoSimplecFoam

and I am not sure if you declare BC in k and omega like that:

    inletXmin
{
    type            freestream;
    freestreamValue      uniform 12.500000;
    value           uniform 12.500000;
}

Never seen that!
But I have never used the freeStream BC.
I am confused, where is your problem?

Well maybe there is a problem of your BC for p and U

time step continuity errors : sum local = 474.441, global = -56.4332, cumulative = -56.4332

[mihaipruna]

Tobi, that patch is not considered because of the internal flow condition which meshes inside the duct and thus eliminates the patches outside that domain.

[Tobi]

Originally Posted by mihaipruna

Tobi, that patch is not considered because of the internal flow condition which meshes inside the duct and thus eliminates the patches outside that domain.

okay I understand
But your case is working, so whats the problem? I just wanted to help you but I think you dont need my help at all.
Hope you get it work and have a nice sunday.

Tobi

[mihaipruna]

hi Tobi, I included the streamlines I got using the incompressible solver simpleFoam in order to showcase the geometry and desired flow. Unfortunately, to perform a realistic and accurate simulation at the inlet speeds I used - 100m/s, a compressible flow solver is required. and that's where it crashes.

[Tobi]

hi,

okay you used a simpleFoam solver to test if its working. Hmmm ... but your continuity crashes too.
Look at your values. Maybe you can have a look at the first time step and see where your problem zones are.

Can you uplaod your case, if yes, I c an have a look at it.
Tobi