[Yassin-K]

Hey guys,
I am trying to run the rivuletPanel tutorial using buoyantReactingFoam and replacing the H2O liquid with liquid nitrogen. However, it shows the following error:

/*---------------------------------------------------------------------------*\ 
  =========                 | 
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox 
   \\    /   O peration     | Website:  https://openfoam.org 
    \\  /    A nd           | Version:  9 
     \\/     M anipulation  | 
\*---------------------------------------------------------------------------*/ 
Build  : 9 
Exec   : buoyantReactingFoam 
Date   : Apr 26 2023 
Time   : 10:26:44 
Host   : "isr028" 
PID    : 31425 
I/O    : uncollated 
Case   : /loc_home/khaled/software/testing/rivuletPanel 
nProcs : 1 
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). 
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (file
ModificationSkew 10) 
allowSystemOperations : Allowing user-supplied system call operations 

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 
Create time 

Create mesh for time = 0 


PIMPLE: No convergence criteria found 


PIMPLE: Operating solver in transient mode with 1 outer corrector 
PIMPLE: Operating solver in PISO mode 


Reading thermophysical properties 

Selecting thermodynamics package  
{ 
    type            heRhoThermo; 
    mixture         multiComponentMixture; 
    transport       sutherland; 
    thermo          janaf; 
    energy          sensibleEnthalpy; 
    equationOfState perfectGas; 
    specie          specie; 
} 

Reading field U 

Reading/calculating face flux field phi 

Creating turbulence model. 

Selecting turbulence model type laminar 
Selecting laminar stress model Stokes 
Creating thermophysical transport model 

Selecting thermophysical transport type laminar 
Selecting default laminar thermophysical transport model unityLewisFourier 
Creating reaction model 

Selecting combustion model none 

Reading g 

Reading hRef 
Calculating field g.h 


Reading pRef 
Reading field p_rgh 

Creating field dpdt 

Creating field kinetic energy K 

No MRF models present 

Creating fvModels from "constant/fvModels" 

Selecting finite volume model type surfaceFilm 
    Name: surfaceFilm 
Selecting surfaceFilmModel thermoSingleLayer 
    Selecting region model functions 
    none 

PIMPLE: No convergence criteria found 


PIMPLE: No corrector convergence criteria found 
        Calculations will do 2 corrections 


PIMPLE: Operating solver in transient mode with 2 outer correctors 


Selecting thermodynamics package  
{ 
    type            heRhoThermo; 
    mixture         pureMixture; 
    properties      liquid; 
    energy          sensibleInternalEnergy; 
} 

#0  Foam::error::printStack(Foam::Ostream&) at ??:? 
#1  Foam::sigFpe::sigHandler(int) at ??:? 
#2  ? in "/lib64/libc.so.6" 
#3  ? in "/lib64/libm.so.6" 
#4  powf64 in "/lib64/libm.so.6" 
#5  Foam::N2::rho(double, double) const at ??:? 
#6  Foam::heRhoThermo<Foam::rhoThermo::composite, Foam::pureMixture<Foam::species::thermo
<Foam::thermophysicalPropertiesSelector<Foam::liquidProperties>, Foam::sensibleInternalEn
ergy> > >::calculate() at ??:? 
#7  Foam::heRhoThermo<Foam::rhoThermo::composite, Foam::pureMixture<Foam::species::thermo
<Foam::thermophysicalPropertiesSelector<Foam::liquidProperties>, Foam::sensibleInternalEn
ergy> > >::heRhoThermo(Foam::fvMesh const&, Foam::word const&) at ??:? 
#8  Foam::rhoThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoThermo::compo
site, Foam::pureMixture<Foam::species::thermo<Foam::thermophysicalPropertiesSelector<Foam
::liquidProperties>, Foam::sensibleInternalEnergy> > > >::New(Foam::fvMesh const&, Foam::
word const&) at ??:? 
#9  Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh c
onst&, Foam::word const&) at ??:? 
#10  Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:? 
#11  Foam::regionModels::surfaceFilmModels::kinematicSingleLayer::kinematicSingleLayer(Fo
am::word const&, Foam::fvMesh const&, Foam::dimensioned<Foam::Vector<double> > const&, Fo
am::word const&, bool) at ??:? 
#12  Foam::regionModels::surfaceFilmModels::thermoSingleLayer::thermoSingleLayer(Foam::wo
rd const&, Foam::fvMesh const&, Foam::dimensioned<Foam::Vector<double> > const&, Foam::wo
rd const&, bool) at ??:? 
#13  Foam::regionModels::surfaceFilmModel::addmeshConstructorToTable<Foam::regionModels::
surfaceFilmModels::thermoSingleLayer>::New(Foam::word const&, Foam::fvMesh const&, Foam::
dimensioned<Foam::Vector<double> > const&, Foam::word const&) at ??:? 
#14  Foam::regionModels::surfaceFilmModel::New(Foam::fvMesh const&, Foam::dimensioned<Foa
m::Vector<double> > const&, Foam::word const&) at ??:? 
#15  Foam::fv::surfaceFilm::surfaceFilm(Foam::word const&, Foam::word const&, Foam::dicti
onary const&, Foam::fvMesh const&) at ??:? 
#16  Foam::fvModel::adddictionaryConstructorToTable<Foam::fv::surfaceFilm>::New(Foam::wor
d const&, Foam::word const&, Foam::dictionary const&, Foam::fvMesh const&) at ??:? 
#17  Foam::fvModel::New(Foam::word const&, Foam::dictionary const&, Foam::fvMesh const&) 
at ??:? 
#18  Foam::fvModels::fvModels(Foam::fvMesh const&) at ??:? 
#19  ? at ??:? 
#20  ? at ??:? 
#21  __libc_start_main in "/lib64/libc.so.6" 
#22  ? at /home/abuild/rpmbuild/BUILD/glibc-2.31/csu/../sysdeps/x86_64/start.S:122 
Gleitkomma-Ausnahme (Speicherabzug geschrieben) 
khaled@isr028:/loc_home/khaled/software/testing/rivuletPanel> ^C 
khaled@isr028:/loc_home/khaled/software/testing/rivuletPanel>  cd /loc_home/khaled/softwa
re/testing/rivuletPanel/0

knowing that I set the inlet temperature to be at 70K to make sure that nitrogen is liquid in this case.

I think OpenFOAM has a probelm obtaining the liquid properties in this case.

Has anybody tried running liquid nitrogen as a film before and can share the case? Any ideas about why it is throwing this error?


Attached my case.


Regards,
Khaled

[ShaneheSSSS]

Hello，have you solved this problem？