|January 13, 2014, 13:57||
SimpleFoam time step continuity error
Join Date: Jan 2014
Posts: 1Rep Power: 0
Firstly, I'm pretty new to OF so please forgive me if I'm doing something really stupid.
I'm trying to run simpleFoam on a structured O-grid mesh of an aerofoil, generated using this tool: http://hvirvel.dk/airfoilmesher/ with coordinates that I've put in as a .dat file. However, when I run I'm getting an error in which my time step continuity errors are exploding, as is the bounding nuTilda value.
The error message is as follows:
Time = 51
smoothSolver: Solving for Ux, Initial residual = 0.545925, Final residual = 0.0316233, No Iterations 40
smoothSolver: Solving for Uy, Initial residual = 0.529254, Final residual = 0.0427123, No Iterations 38
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#3 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:?
#5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#7 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#11 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
Floating point exception (core dumped)
I have searched and tried the following:
- Setting gradSchemes in fvSchemes to cellLimited 1. This then causes the simulation to not explode, but it does not converge; it diverges very slowly. After time = 5000 it was still nowhere near converged.
- Running with turbulence off. This obviously removes the nuTilda exploding, but the time step continuity errors still explode.
- Using a coarser mesh. This speeds up the error message appearing, but still doesn't fix the problem.
Using an unstructured mesh using gmsh with all the boundary conditions in the 0 folder the same works fine, and converges no problem.
I can't work out what exactly is causing the problem, and any help or suggestions would be greatly appreciated.
My files are all in the zip file here: https://www.dropbox.com/s/hcyfj9yby0...structured.zip
Thanks in advance for any help anyone can give.
|Thread||Thread Starter||Forum||Replies||Last Post|
|mixerVesselAMI2D's mass is not balancing||sharonyue||OpenFOAM Running, Solving & CFD||6||June 10, 2013 09:34|
|pisoFoam with k-epsilon turb blows up - Some questions||Heroic||OpenFOAM Running, Solving & CFD||26||December 17, 2012 04:34|
|Ansys Fluent 13.0 UDF compilation problem in Window XP (32 bit)||Yogini||Fluent UDF and Scheme Programming||7||October 3, 2012 07:24|
|checking the system setup and Qt version||vivek070176||OpenFOAM Installation||22||June 1, 2010 12:34|
|Could anybody help me see this error and give help||liugx212||OpenFOAM Running, Solving & CFD||3||January 4, 2006 19:07|