[jet_engine]

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)):

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 ??:?

[betakv]

Dear Max, there is an error in initial and boundary conditon for omega. please try to use value 162.5 instead 0

[jet_engine]

Originally Posted by betakv

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?