CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Running icoFoam. Problems with a 2D Simulation

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

LinkBack Thread Tools Display Modes
Old   October 2, 2014, 23:23
Default Running icoFoam. Problems with a 2D Simulation
New Member
Esteban Vargas
Join Date: Jul 2014
Location: Medellin, Colombia
Posts: 8
Rep Power: 5
vbesteban is on a distinguished road

I'm running this simple problem in openfoam, with a no solution founded.

The problem is the flow of water around a cylinder, I've decided to do it in 2d, because the next simulation will be the flow around a sphere in which I will have to find drag and lift coefficients but in 3D. The inlet velocity is 0.1 m/s, and the pressure is 0.

The error is:
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/"
#2 in "/lib/x86_64-linux-gnu/"
#3 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/"
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/icoFoam"
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/icoFoam"
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/icoFoam"
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/icoFoam"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/"
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/icoFoam"
Excepción de coma flotante (`core' generado)

I have already change the inlet and outlet types of BC, but the problem still going. I did the mesh in gmsh, and exported with no problem to openFoam. I put here the all information. I really appreciate your time, and the way in which you can help me.

The files are here:

Again I really appreciate your time. I'm using OF 2.3 and postprocessing in paraView.

vbesteban is offline   Reply With Quote

Old   October 3, 2014, 02:21
Senior Member
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,794
Rep Power: 30
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym

well... icoFoam is applicable only for non-turbulent cases. While in your case Re = 0.1*D/1e-6. Visually I can estimate diameter of the cylinder as 0.25 of channel width, so it's 2800*0.25 = 700. So final Re = 7e6, it's rather high and solver diverges.
alexeym is offline   Reply With Quote

Old   October 9, 2014, 13:33
Default Problems with time step continuity!
New Member
Esteban Vargas
Join Date: Jul 2014
Location: Medellin, Colombia
Posts: 8
Rep Power: 5
vbesteban is on a distinguished road
Hi Alexey.

Well I did the change and I'm using simpleFoam to solve the case, but I have another problem again.

PHP Code:
Time 1

:  Solving for UxInitial residual 1, Final residual 0.0608736No Iterations 2
:  Solving for UyInitial residual 1, Final residual 0.0635718No Iterations 2
:  Solving for pInitial residual 1, Final residual 0.0224656No Iterations 10
time step continuity errors 
sum local 1.2504e-07, global = 4.29951e-21cumulative 4.29951e-21
:  Solving for epsilonInitial residual 1, Final residual 2.40018e-13No Iterations 2
:  Solving for kInitial residual 1, Final residual 1.59607e-10No Iterations 2
bounding k
min0 max0.616825 average0.500641
0.08 s  ClockTime 1 s


:  Solving for UxInitial residual 0.393844, Final residual 0.0236749No Iterations 2
:  Solving for UyInitial residual 0.400707, Final residual 0.0234863No Iterations 2
:  Solving for pInitial residual 0.994541, Final residual 0.0102574No Iterations 2
time step continuity errors 
sum local 3.13913e-05, global = 6.02984e-21cumulative 1.03293e-20
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/"
#2   in "/lib/x86_64-linux-gnu/"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/"
#4  void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/"
#5   at kEpsilon.C:0
#6  Foam::incompressible::RASModels::kEpsilon::correct() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/"
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/simpleFoam"
#8  __libc_start_main in "/lib/x86_64-linux-gnu/"
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/simpleFoam"
Excepción de coma flotante (`core' generado) 
And I don't get any solution at all. I've changed the epsilon values as some people indicates in the forum and its didn't work.

The problem is a flow around a cylinder...its seems easy to perform!

I really appreciate that somebody can help me!

Thanks alexey and all who support me.

The new files are here.
vbesteban is offline   Reply With Quote


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
[BEGINNER] Problems running tutorials (incompressible) in OF 2.xx Jaro OpenFOAM Running, Solving & CFD 8 March 16, 2016 08:27
CFX thermo-fluid 2D simulation problems maryliz CFX 6 November 1, 2011 00:26
GUI crash and simulation engine still running RPJones FLOW-3D 2 November 9, 2010 09:18
Replacing mesh while running a simulation akultane CFX 1 November 15, 2009 14:46
Problems about multi swirl tubes simulation using Jason FLUENT 4 May 9, 2008 00:33

All times are GMT -4. The time now is 00:07.