CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > 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 Display Modes
Old   October 23, 2008, 08:48
Default Dear users, I have read man
  #1
Senior Member
 
stephane sanchi
Join Date: Mar 2009
Posts: 300
Rep Power: 9
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
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Command line tool for finding the bounding box of an STL file terrybarnaby OpenFOAM Mesh Utilities 7 November 11, 2013 02:13
SnappyHexMesh plane surface bounding box podallaire OpenFOAM Bugs 22 August 26, 2009 07:34
Bounding problem in running rasinterfoam qtian OpenFOAM Running, Solving & CFD 4 June 30, 2008 22: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 08:55


All times are GMT -4. The time now is 03:14.