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

Simulation not converging (core dumped)

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

Like Tree1Likes
  • 1 Post By betakv

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 25, 2016, 06:49
Unhappy Simulation not converging (core dumped)
  #1
New Member
 
Max Vorstadt
Join Date: May 2016
Posts: 28
Rep Power: 10
jet_engine is on a distinguished road
Hi! I have been trying to run simpleFoam on a simulation but apparently, there is some sort of problem (probably regarding the BCs) which makes the simulation not converge and thus, stop running without obtaining results. I'm trying to run a simulation with the k-w SST turbulence model. I am attaching the 0 folder, and in addition a screenshot where all the BCs are detailed. I've tried limiting the Courant number (maxCo 0.5 and setting a variable timestep but it didn't work, and neither did it work when I set an inflow field of 0 instead of one obtained from potentialFoam)

The other screenshots represent the environment so the case is more clear. But basically, a velocity profile enters in the x direction, and it encounters a cylinder and a plate, expecting a wake behind the cylinder. The values of the turbulence parameters have been calculated according to the formulas from the CFD Online wiki (http://www.cfd-online.com/Wiki/Turbu...ary_conditions).

Does anybody see which BC(s) or other problem could be wrong making the simulation not converge?

That's the output of the log.simpleFoam file (the output is Floating point exception (core dumped)):
Code:
Create time

Create mesh for time = 0


SIMPLE: no convergence criteria found. Calculations will run for 500 steps.

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
bounding k, min: 0 max: 0.02535 average: 0
bounding omega, min: 0 max: 162.5 average: 0
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5  Foam::RASModels::kOmegaSST<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::F2() const at ??:?
#6  Foam::RASModels::kOmegaSST<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::F23() const at ??:?
#7  Foam::RASModels::kOmegaSST<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correctNut(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#8  Foam::RASModels::kOmegaSST<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correctNut() at ??:?
#9  Foam::RASModels::kOmegaSST<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::kOmegaSST(Foam::geometricOneField const&, Foam::geometricOneField const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel const&, Foam::word const&, Foam::word const&) at ??:?
#10  Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::adddictionaryConstructorToTable<Foam::RASModels::kOmegaSST<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > >::New(Foam::geometricOneField const&, Foam::geometricOneField const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel const&, Foam::word const&) at ??:?
#11  Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::New(Foam::geometricOneField const&, Foam::geometricOneField const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel const&, Foam::word const&) at ??:?
#12  Foam::TurbulenceModel<Foam::geometricOneField, Foam::geometricOneField, Foam::incompressibleTurbulenceModel, Foam::transportModel>::adddictionaryConstructorToTable<Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > >::NewTurbulenceModel(Foam::geometricOneField const&, Foam::geometricOneField const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel const&, Foam::word const&) at ??:?
#13  ? at ??:?
#14  ? at ??:?
#15  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#16  ? at ??:?
Attached Images
File Type: png Screen Shot 2016-05-25 at 12.47.07.png (43.8 KB, 6 views)
File Type: png Screen Shot 2016-05-25 at 12.30.37.png (14.0 KB, 6 views)
File Type: png Screen Shot 2016-05-25 at 12.30.31.png (41.4 KB, 7 views)
File Type: png Screen Shot 2016-05-25 at 12.46.05.png (81.8 KB, 6 views)
Attached Files
File Type: gz 0.tar.gz (47.3 KB, 5 views)

Last edited by jet_engine; May 25, 2016 at 06:55. Reason: Added the 0 folder
jet_engine is offline   Reply With Quote

Old   May 25, 2016, 08:20
Default
  #2
Member
 
Vojtech Betak
Join Date: Mar 2009
Location: Czech republic
Posts: 34
Rep Power: 18
betakv is on a distinguished road
Dear Max, there is an error in initial and boundary conditon for omega. please try to use value 162.5 instead 0
jet_engine likes this.
betakv is offline   Reply With Quote

Old   May 25, 2016, 08:31
Default
  #3
New Member
 
Max Vorstadt
Join Date: May 2016
Posts: 28
Rep Power: 10
jet_engine is on a distinguished road
Quote:
Originally Posted by betakv View Post
Dear Max, there is an error in initial and boundary conditon for omega. please try to use value 162.5 instead 0
Dear Vojtech,

First of all thank you for your answer! Do you mean the internalField uniform value 0, the values of w at the wall or both? I have tried changing the ones at the walls in addition to the inlet one and the output I obtain is the same.

Edit: Now I've tried setting all them to 162.5 and apparently it runs. Is setting the whole domain's w to the same value physically correct?
Furthermore, is it the same with nut and k?

Last edited by jet_engine; May 25, 2016 at 08:56. Reason: Tested both
jet_engine is offline   Reply With Quote

Reply

Tags
boundary condition, converge, kwsstmodel, simplefoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
icoFoam - Floating point exception (core dumped) File OpenFOAM Running, Solving & CFD 26 March 6, 2021 04:26
If nothing works of converging simulation... miquiita CFX 1 June 7, 2017 22:42
Shall I leave one core free when I run CFD simulation? Pierre1 Main CFD Forum 1 March 12, 2015 04:46
Segmentation fault (core dumped) in Axial_Rotor_MRF sam.ho OpenFOAM Running, Solving & CFD 0 March 17, 2014 01:42
mapField & core dumped sasanghomi OpenFOAM 4 August 23, 2013 10:10


All times are GMT -4. The time now is 12:32.