CFD Online Discussion Forums - Running icoFoam. Problems with a 2D Simulation

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - Running icoFoam. Problems with a 2D Simulation (https://www.cfd-online.com/Forums/openfoam-solving/142514-running-icofoam-problems-2d-simulation.html)

Running icoFoam. Problems with a 2D Simulation

Dears.

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::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/icoFoam"
#5
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/icoFoam"
#6
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/icoFoam"
#7
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/icoFoam"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
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:

https://www.dropbox.com/s/zjtkw3vuwe...hp.tar.gz?dl=0
https://www.dropbox.com/s/5b7xtereey...drohp.zip?dl=0

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

Thanks.

Hi,

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.

Problems with time step continuity!

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

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0608736, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0635718, No Iterations 2
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0224656, No Iterations 10
time step continuity errors : sum local = 1.2504e-07, global = 4.29951e-21, cumulative = 4.29951e-21
smoothSolver:  Solving for epsilon, Initial residual = 1, Final residual = 2.40018e-13, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 1.59607e-10, No Iterations 2
bounding k, min: 0 max: 0.616825 average: 0.500641
ExecutionTime = 0.08 s  ClockTime = 1 s

Time = 2

smoothSolver:  Solving for Ux, Initial residual = 0.393844, Final residual = 0.0236749, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.400707, Final residual = 0.0234863, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.994541, Final residual = 0.0102574, No Iterations 2
time step continuity errors : sum local = 3.13913e-05, global = 6.02984e-21, cumulative = 1.03293e-20
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#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/libincompressibleRASModels.so"
#5   at kEpsilon.C:0
#6  Foam::incompressible::RASModels::kEpsilon::correct() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#7  
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/simpleFoam"
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  
 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.

https://www.dropbox.com/s/87h6ekxf9p...ce.tar.gz?dl=0