Problem solving a cylinder case 2d. Simple Foam

vbesteban · October 16, 2014, 17:51

Hello.

I'm trying to solve a problem in 2 dimensions. The problem consists in flow around a cylinder. The inlet velocity is 0.1 m/s the pressure is 0, the fluid is water so I've already changed the viscosity in transport properties. But when I run the code this is the result.

And in Parafoam I can't see a full development of the problem. It seems that I don't get any type of result. The fluid is considered as external flow.

I'm using OF2.3 and meshing in gmsh.

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)

I've changed the epsilon values as some people indicates in the forum but this didn't work.

I really appreciate that somebody can help me!

Thanks to all who support me.

The files are here.

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

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
problem with Min/max rho	tH3f0rC3	OpenFOAM	8	July 31, 2019 10:48
a problem with convergence in buoyantSimpleFoam	skuznet	OpenFOAM Running, Solving & CFD	6	November 15, 2017 13:12
conjugate heat transfer in OpenFOAM	skuznet	OpenFOAM Running, Solving & CFD	99	March 16, 2017 06:07
Mass conservation problem with compressibleInterDyMFoam	LESlie	OpenFOAM Running, Solving & CFD	4	July 21, 2015 05:05
pimpleFoam: turbulence->correct(); is not executed when using residualControl	hfs	OpenFOAM Running, Solving & CFD	3	October 29, 2013 09:35

October 16, 2014, 17:51	Problem solving a cylinder case 2d. Simple Foam	#1
vbesteban New Member Esteban Vargas Join Date: Jul 2014 Location: Medellin, Colombia Posts: 8 Rep Power: 12	Hello. I'm trying to solve a problem in 2 dimensions. The problem consists in flow around a cylinder. The inlet velocity is 0.1 m/s the pressure is 0, the fluid is water so I've already changed the viscosity in transport properties. But when I run the code this is the result. And in Parafoam I can't see a full development of the problem. It seems that I don't get any type of result. The fluid is considered as external flow. I'm using OF2.3 and meshing in gmsh. 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) I've changed the epsilon values as some people indicates in the forum but this didn't work. I really appreciate that somebody can help me! Thanks to all who support me. The files are here. https://www.dropbox.com/s/87h6ekxf9p...ce.tar.gz?dl=0[/QUOTE]