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

rhoSImplecFoam still having issues!

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 8, 2012, 12:02
Default rhoSImplecFoam still having issues!
  #1
Senior Member
 
Mihai Pruna
Join Date: Apr 2010
Location: Boston
Posts: 195
Rep Power: 15
mihaipruna is on a distinguished road
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$
Attached Files
File Type: gz 0 (copy).tar.gz (1.3 KB, 19 views)
File Type: gz constant (copy).tar.gz (1.0 KB, 10 views)
File Type: gz system.tar.gz (3.8 KB, 13 views)
mihaipruna is offline   Reply With Quote

Old   August 8, 2012, 12:04
Default
  #2
Senior Member
 
Mihai Pruna
Join Date: Apr 2010
Location: Boston
Posts: 195
Rep Power: 15
mihaipruna is on a distinguished road
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 is offline   Reply With Quote

Old   August 8, 2012, 12:06
Default
  #3
Senior Member
 
Mihai Pruna
Join Date: Apr 2010
Location: Boston
Posts: 195
Rep Power: 15
mihaipruna is on a distinguished road
Image of a solved case (using simpleFoam)
mihaipruna is offline   Reply With Quote

Old   August 10, 2012, 14:22
Default
  #4
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
you still have that problem? If yes , I ll have a look at your case.

Tobi
Tobi is offline   Reply With Quote

Old   August 10, 2012, 16:09
Default
  #5
Senior Member
 
Mihai Pruna
Join Date: Apr 2010
Location: Boston
Posts: 195
Rep Power: 15
mihaipruna is on a distinguished road
as of now, yes, no matter what I do I am experiencing crashes rather soon in the solution.
Thank you.
mihaipruna is offline   Reply With Quote

Old   August 10, 2012, 16:25
Default
  #6
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
doublepost / sorry
Tobi is offline   Reply With Quote

Old   August 10, 2012, 16:26
Default
  #7
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to 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:

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
Tobi is offline   Reply With Quote

Old   August 11, 2012, 19:04
Default
  #8
Senior Member
 
Mihai Pruna
Join Date: Apr 2010
Location: Boston
Posts: 195
Rep Power: 15
mihaipruna is on a distinguished road
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.
mihaipruna is offline   Reply With Quote

Old   August 12, 2012, 05:06
Default
  #9
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by mihaipruna View Post
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
Tobi is offline   Reply With Quote

Old   August 12, 2012, 10:41
Default
  #10
Senior Member
 
Mihai Pruna
Join Date: Apr 2010
Location: Boston
Posts: 195
Rep Power: 15
mihaipruna is on a distinguished road
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.
mihaipruna is offline   Reply With Quote

Old   August 12, 2012, 11:30
Default
  #11
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to 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
Tobi is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Porous media setup issues in Fluent Bernard Van FLUENT 29 January 26, 2017 05:09
rhoSimplecFoam U-turn j-blindi OpenFOAM Running, Solving & CFD 6 October 29, 2011 18:21
License and Network issues scottneh STAR-CCM+ 2 September 12, 2011 19:27
FLUENT Speed Issues on Cluster cfd23 FLUENT 2 April 4, 2010 00:43
Grid size, convergence issues franzdrs Main CFD Forum 3 June 18, 2009 08:57


All times are GMT -4. The time now is 16:45.