View Issue Details

IDProjectCategoryView StatusLast Update
0003871OpenFOAMBugpublic2022-08-02 12:16
Reporteragustinvo Assigned Tohenry  
PrioritynormalSeverityminorReproducibilityalways
Status closedResolutionno change required 
PlatformGNU/LinuxOSUbuntuOS Version20.04
Product Version10 
Fixed in Version10 
Summary0003871: Can't compute nu when using transport const and rhoConst
DescriptionWhen running a simulation using rhoSimpleFoam

thermoType
{
    ...
    transport const;
    equationOfState rhoConst;
    ...
}

An error appears in the selected compressibleMomentumTransportModel. The turbulence model cannot calculate nu from mu and rho. I tried with several models but the same error appears.

"""
#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::fluidThermo::nu() const at ??:?
#6 Foam::RASModel<Foam::compressibleMomentumTransportModel>::nu() const at ??:?
#7 ? at compressibleMomentumTransportModels.C:?
#8 Foam::linearViscousStress<Foam::RASModel<Foam::compressibleMomentumTransportModel> >::divDevTau(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const at ??:?
"""
Steps To Reproduce1. Run rhoSimpleFoam with constant mu and rho
2. Get an error whathever the turbulence model is
Additional InformationA test case is attached. Just use blockMesh and rhoSimpleFoam
TagsNo tags attached.

Activities

agustinvo

2022-08-02 09:08

reporter  

channelFlow.zip (11,898 bytes)

henry

2022-08-02 09:23

manager   ~0012700

What messages do you get when you run the case you attached? I get

--> FOAM FATAL ERROR:
cannot find file ".../channelFlow/0/R"

agustinvo

2022-08-02 10:08

reporter   ~0012701

Sorry, I attached a case that uses the SSG model. Here there is a LaunderSharmaKE simulation
channelFlow-2.zip (11,906 bytes)

henry

2022-08-02 12:15

manager   ~0012702

You are using the compressibility-based thermo hePsiThermo with an incompressible equation of state so the density evaluated from the compressibility will be 0.

I could not work with your case as the setup is no convoluted with the data moved around the files so I setup the squareBend with

thermoType
{
    type heRhoThermo;
    mixture pureMixture;
    transport const;
    thermo hConst;
    equationOfState rhoConst;
    specie specie;
    energy sensibleEnthalpy;
}

mixture
{
    specie
    {
        molWeight 28.9;
    }
    equationOfState
    {
        rho 1;
    }
    thermodynamics
    {
        Cp 1005;
        Hf 0;
    }
    transport
    {
        mu 1.82e-05;
        Pr 0.71;
    }
}

and it runs fine.

henry

2022-08-02 12:16

manager   ~0012703

User error

Issue History

Date Modified Username Field Change
2022-08-02 09:08 agustinvo New Issue
2022-08-02 09:08 agustinvo File Added: channelFlow.zip
2022-08-02 09:23 henry Note Added: 0012700
2022-08-02 10:08 agustinvo File Added: channelFlow-2.zip
2022-08-02 10:08 agustinvo Note Added: 0012701
2022-08-02 12:15 henry Note Added: 0012702
2022-08-02 12:16 henry Assigned To => henry
2022-08-02 12:16 henry Status new => closed
2022-08-02 12:16 henry Resolution open => no change required
2022-08-02 12:16 henry Fixed in Version => 10
2022-08-02 12:16 henry Note Added: 0012703