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

Bounding epsilon and K with rasInterFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 23, 2008, 09:48
Default Dear users, I have read man
  #1
Senior Member
 
stephane sanchi
Join Date: Mar 2009
Posts: 314
Rep Power: 18
openfoam_user is on a distinguished road
Dear users,

I have read many threads about 'bounding problems for k and epsilon', but without finding how to solve it.

I am trying to compute a free surface calculation with a boat.

The mesh was generated using ICEMCFD Hexa.

The cae crashes after 2-4 iterations with the following message :

MULES: Solving for gamma
Liquid phase volume fraction = 0.499998 Min(gamma) = -2.96063e-12 Max(gamma) = 1.01884
MULES: Solving for gamma
Liquid phase volume fraction = 0.499998 Min(gamma) = -2.96285e-12 Max(gamma) = 1.01883
MULES: Solving for gamma
Liquid phase volume fraction = 0.499998 Min(gamma) = -2.96507e-12 Max(gamma) = 1.01883
MULES: Solving for gamma
Liquid phase volume fraction = 0.499998 Min(gamma) = -2.96729e-12 Max(gamma) = 1.01882
DILUPBiCG: Solving for Ux, Initial residual = 0.611633, Final residual = 8.16394e-07, No Iterations 625
DILUPBiCG: Solving for Uy, Initial residual = 0.673241, Final residual = 8.43185e-07, No Iterations 631
DILUPBiCG: Solving for Uz, Initial residual = 0.927937, Final residual = 8.75002e-07, No Iterations 564
DICPCG: Solving for pd, Initial residual = 1, Final residual = 0.038253, No Iterations 4
DICPCG: Solving for pd, Initial residual = 6.54561e-09, Final residual = 6.54561e-09, No Iterations 0
DICPCG: Solving for pd, Initial residual = 6.54561e-09, Final residual = 6.54561e-09, No Iterations 0
DICPCG: Solving for pd, Initial residual = 1.36844e-07, Final residual = 1.2005e-08, No Iterations 1
DICPCG: Solving for pd, Initial residual = 3.54892e-08, Final residual = 3.54892e-08, No Iterations 0
DICPCG: Solving for pd, Initial residual = 3.54892e-08, Final residual = 3.54892e-08, No Iterations 0
DICPCG: Solving for pd, Initial residual = 1.63417e-07, Final residual = 1.93669e-08, No Iterations 1
DICPCG: Solving for pd, Initial residual = 4.80168e-08, Final residual = 4.80168e-08, No Iterations 0
DICPCG: Solving for pd, Initial residual = 4.80168e-08, Final residual = 4.80168e-08, No Iterations 0
time step continuity errors : sum local = 0.000631323, global = 8.45631e-07, cumulative = 8.45629e-07
DILUPBiCG: Solving for epsilon, Initial residual = 0.870335, Final residual = 32077.1, No Iterations 2001
bounding epsilon, min: -6.62374e+16 max: 1.14934e+17 average: -7.2268e+11
DILUPBiCG: Solving for k, Initial residual = 0.995697, Final residual = 8.46144e-05, No Iterations 2001
bounding k, min: -2.04339e+12 max: 8.30451e+14 average: 2.41778e+11
ExecutionTime = 282.88 s ClockTime = 284 s

Courant Number mean: 0.252986 max: 452593
deltaT = 3.58683e-11
Time = 0.00503247

MULES: Solving for gamma
[1] [3] ##00 Foam::error::printStack(Foam:stream&)Foam::error::printStack(Foam:stream&) in "/shared/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigFpeHandler(int) in "/shared/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigFpeHandler(int)[2] #0 Foam::error::printStack(Foam:stream&)[0] in "/shared/OpenF#O0 AFoam::error::printStack(Foam:stream&)M/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #2 in [5] "#/0share dFoam::error::printStack(Foam:stream&)/OpenFOAM/Ope[4] n#F0O AFoam::error::printStack(Foam:stream&)M-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #2 ?? in "/shared/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[2] #1 Foam::sigFpe::sigFpeHandler(int) in "/shared/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigFpeHandler(int)?? in in "/shared/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[5] #1 Foam::sigFpe::sigFpeHandler(int)"/shared/OpenFOAM in /O"p/elniFbO6A4M/-l1i.b5c/.lsiob./6l"i
n[1] u#x36 4Gvoid Foam::MULES::limiter<foam::onefield,>(Foam::Field< double>&, Foam::oneField const&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int)ccDPOpt/libOpenFOAM.so"
[4] #1 Foam::sigFpe::sigFpeHandler(int) in "/shared/OpenFOAM/OpenFOAM-1.5 in /"l/isbh/alriendu/xO6p4eGncFcODAPMO/pOtp/elniFbOOApMe-n1F.O5A/Ml.isbo/"l
i[2] n#u2x 64GccDPOpt/libOpenFOAM.so"
[0] #2 in "/lib64/libc.so.6"

The error message is much longer. This is only the beginning.

I have tried to play with many parameters (deltaT, maxDeltaT, k and setFields, ...) but without success !!!!!!!!!!!!

Stephane.
openfoam_user 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
[surface handling] Command line tool for finding the bounding box of an STL file terrybarnaby OpenFOAM Meshing & Mesh Conversion 8 November 4, 2016 14:09
SnappyHexMesh plane surface bounding box podallaire OpenFOAM Bugs 22 August 26, 2009 08:34
Bounding problem in running rasinterfoam qtian OpenFOAM Running, Solving & CFD 4 June 30, 2008 23:54
gridding when bounding surfaces are rough? Herb Schilling Main CFD Forum 1 December 1, 2000 06:58
Bounding frames increments Ricky Wong FLUENT 3 April 6, 2000 09:55


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