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

KEpsilon error

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 14, 2008, 18:37
Default Hi I ran a case with kEpsi
  #1
Member
 
mohd mojab
Join Date: Mar 2009
Posts: 31
Rep Power: 8
mou_mi is on a distinguished road
Hi

I ran a case with kEpsilon model and set the wall for "k" and "epsilon" as a fixedValue for B.C. and OF gave me this error for "k" and "epsilon",


Create mesh for time = 70

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model kEpsilon

Starting time loop

Time = 70.001

Courant Number mean: 0.00103968 max: 0.0494671
DILUPBiCG: Solving for Ux, Initial residual = 5.85686e-05, Final residual = 4.6585e-10, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.000515567, Final residual = 3.09483e-09, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.00109401, Final residual = 7.4057e-09, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.942377, Final residual = 9.56964e-07, No Iterations 273
time step continuity errors : sum local = 2.11503e-12, global = 2.74312e-14, cumulative = 2.74312e-14
DICPCG: Solving for p, Initial residual = 0.000860842, Final residual = 9.75e-07, No Iterations 164
time step continuity errors : sum local = 4.12992e-12, global = -7.69118e-14, cumulative = -4.94806e-14


--> FOAM FATAL ERROR : fixedValue is the wrong epsilon patchField type for wall-functions on patch wall
should be zeroGradient

From function wall-function evaluation
in file /home/OpenFOAM/OpenFOAM-1.4.1/src/finiteVolume/lnInclude/checkPatchFieldTypes.H at line 13.

FOAM exiting


I changed them to "zeroGradient" for wall boundary, after some calculation in the first time step it gave this error.

Create mesh for time = 70

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model kEpsilon

Starting time loop

Time = 70.001

Courant Number mean: 0.00103968 max: 0.0494671
DILUPBiCG: Solving for Ux, Initial residual = 5.85686e-05, Final residual = 4.6585e-10, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.000515567, Final residual = 3.09483e-09, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.00109401, Final residual = 7.4057e-09, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.942377, Final residual = 9.56964e-07, No Iterations 273
time step continuity errors : sum local = 2.11503e-12, global = 2.74312e-14, cumulative = 2.74312e-14
DICPCG: Solving for p, Initial residual = 0.000860842, Final residual = 9.75e-07, No Iterations 164
time step continuity errors : sum local = 4.12992e-12, global = -7.69118e-14, cumulative = -4.94806e-14
#0 Foam::error::printStack(Foam:stream&) in "/usr/local/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/usr/local/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt in "/lib64/tls/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/usr/local/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<foam::fvpatchfield,>(Foam::GeometricF ield<double,>&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&) in "/usr/local/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libincompressibleTurbule nceModels.so"
#5 Foam::tmp<foam::geometricfield<double,> > Foam::operator/<foam::fvpatchfield,>(Foam::tmp<foam::geometricfie ld<double,> > const&, Foam::GeometricField<double,> const&) in "/usr/local/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libincompressibleTurbule nceModels.so"
#6 Foam::turbulenceModels::kEpsilon::correct() in "/usr/local/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libincompressibleTurbule nceModels.so"
#7 main in "/usr/local/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/turbFoam"
#8 __libc_start_main in "/lib64/tls/libc.so.6"
#9 Foam::regIOobject::readIfModified() in "/usr/local/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/turbFoam"
Floating exception




would you help me what happened here?

Thank you
mou
mou_mi is offline   Reply With Quote

Old   July 15, 2008, 03:18
Default Hi, Floating Point Exceptio
  #2
Member
 
Andrew King
Join Date: Mar 2009
Location: Perth, Western Australia, Australia
Posts: 81
Rep Power: 8
andersking is on a distinguished road
Hi,

Floating Point Exception means that some of you numbers are unbounded - Inf (infinity) or NaN (not a number). Sometimes this can be overcome by executing
unset FOAM_SIGFPE
and running your simulation again (this basically stops checking for out of bounds errors).
Most of the time this won't fix the problem, however, it is more likely to be a problem with boundary conditions or the mesh.

Cheers
Andrew
__________________
Dr Andrew King
Fluid Dynamics Research Group
Curtin University
andersking is offline   Reply With Quote

Old   July 15, 2008, 04:00
Default Hi Mou From your error-mess
  #3
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,619
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Mou

From your error-messages, especially #3 and #6, tells me that you are dividing by 0 in your turbulence model. Do you on any of your boundaries, inlet, outlet, etc, set either k or epsilon to 0 in stead of 1e-11?
My experience is that these properties cannot be zero, as you are dividing with them. I also believe that it would be necessary to put k!=0 in the interior, but you probably already have that.

Best regards,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   July 16, 2008, 18:25
Default hi Niels as you mentioned,
  #4
Member
 
mohd mojab
Join Date: Mar 2009
Posts: 31
Rep Power: 8
mou_mi is on a distinguished road
hi Niels

as you mentioned, I set 0 for the k and it was the problem.

Thank you
mou
mou_mi 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
Understanding k from kepsilon markh83 OpenFOAM Post-Processing 3 December 5, 2008 04:42
Add new RASModel kEpsilon modification ivanwhlau OpenFOAM Running, Solving & CFD 3 August 21, 2008 04:36
About kEpsilon turbulence model osimonsimon OpenFOAM Running, Solving & CFD 10 April 24, 2008 02:52
Kepsilon BC and intialization podallaire OpenFOAM Running, Solving & CFD 0 October 17, 2007 20:29
KEpsilon and KOmega larry OpenFOAM Running, Solving & CFD 3 June 29, 2006 01:38


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