|
[Sponsors] |
May 25, 2016, 06:49 |
Simulation not converging (core dumped)
|
#1 |
New Member
Max Vorstadt
Join Date: May 2016
Posts: 28
Rep Power: 10 |
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 ??:? Last edited by jet_engine; May 25, 2016 at 06:55. Reason: Added the 0 folder |
|
May 25, 2016, 08:20 |
|
#2 |
Member
Vojtech Betak
Join Date: Mar 2009
Location: Czech republic
Posts: 34
Rep Power: 18 |
Dear Max, there is an error in initial and boundary conditon for omega. please try to use value 162.5 instead 0
|
|
May 25, 2016, 08:31 |
|
#3 | |
New Member
Max Vorstadt
Join Date: May 2016
Posts: 28
Rep Power: 10 |
Quote:
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 |
||
Tags |
boundary condition, converge, kwsstmodel, simplefoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |