CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   rhoPimpleFoam k-omega sst crashes channel flow (https://www.cfd-online.com/Forums/openfoam-solving/138928-rhopimplefoam-k-omega-sst-crashes-channel-flow.html)

hxaxtma July 14, 2014 10:05

rhoPimpleFoam k-omega sst crashes channel flow
 
1 Attachment(s)
Hi guys,

I am trying to simulate a channelFlow at Re=5000 with rhoPimpleFoam.
Problem is my simulation crashes immediately:
Code:

Create time

Create mesh for time = 0


PIMPLE: no residual control data found. Calculations will employ 10 corrector loops

Reading thermophysical properties

Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture        pureMixture;
    transport      const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

AMI: Creating addressing and weights between 30912 source faces and 30912 target faces
AMI: Patch source sum(weights) min/max/average = 1, 1, 1
AMI: Patch target sum(weights) min/max/average = 1, 1, 1
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
#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  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&) at ??:?
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#6  Foam::compressible::RASModels::kOmegaSST::F2() const at ??:?
#7  Foam::compressible::RASModels::kOmegaSST::F23() const at ??:?
#8  Foam::compressible::RASModels::kOmegaSST::kOmegaSST(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&, Foam::word const&) at ??:?
...............

You can find my inlet conditions in the attachement

Where is the problem?

Thanks in advance!

vonboett July 15, 2014 05:04

you devide by zero, and I guess your initial zero pressure field is the reason.

hxaxtma July 15, 2014 05:18

Quote:

Originally Posted by vonboett (Post 501526)
you devide by zero, and I guess your initial zero pressure field is the reason.

yes that was the reason, thank you!


All times are GMT -4. The time now is 14:19.