CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   rhoSImplecFoam still having issues! (http://www.cfd-online.com/Forums/openfoam-solving/105762-rhosimplecfoam-still-having-issues.html)

mihaipruna August 8, 2012 11:02

rhoSImplecFoam still having issues!
 
3 Attachment(s)
Please help! :(


Code:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  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 August 8, 2012 11:04

checkmesh results:

Code:

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 August 8, 2012 11:06

Image of a solved case (using simpleFoam)
http://cadnexus.com/images/stories/sductslide.jpg

Tobi August 10, 2012 13:22

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

Tobi

mihaipruna August 10, 2012 15:09

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

Tobi August 10, 2012 15:25

doublepost / sorry

Tobi August 10, 2012 15:26

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:

Code:

    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 August 11, 2012 18:04

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 August 12, 2012 04:06

Quote:

Originally Posted by mihaipruna (Post 376622)
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 August 12, 2012 09:41

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 August 12, 2012 10:30

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


All times are GMT -4. The time now is 09:30.