[iatnart]

Dear all,
I'm new on OpenFoam. I'm training with the models in this link:
https://github.com/iatnart/ac-room
I'm trying to solve the 'buoyantFoam' and I get the message

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 10
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 10-c4cf895ad8fa
Exec : buoyantFoam
Date : Jun 23 2023
Time : 15:04:29
Host : "DESKTOP-H5IDDY0"
PID : 588
I/O : uncollated
Case : /mnt/e/googledrive/05-OpenFOAM/02.Project/01-SGF_LA/00.S6-bouyant
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: Convergence criteria found
p_rgh: tolerance 0.0001
U: tolerance 0.0001
h: tolerance 0.0001
"(k|epsilon|omega)": tolerance 0.001


PIMPLE: Operating solver in steady-state mode with 1 outer corrector
PIMPLE: Operating solver in SIMPLE mode


Reading thermophysical properties

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

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
model kOmegaSST;
turbulence on;
printCoeffs on;
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.555556;
gamma2 0.44;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}

Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity

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

No fvModels present
No fvConstraints present


I don't know to fix this and why I get the problem.

Thank for all,

[Yann]

Hello,

No MRF models present

No fvModels present
No fvConstraints present

Those are usual startup outputs for pretty much any solver. It's only here to inform you there is no MRF model nor fvModels/fvConstraints used in your case

It's perfectly fine, unless you were trying to use one of these functionalities.
These are optional features, your solver does not need it to run you case. If the solver crashes, it's related to something else.

Don't you get an error when it stops?

Regards,
Yann

[iatnart]

Originally Posted by Yann

Hello,

No MRF models present

No fvModels present
No fvConstraints present

Those are usual startup outputs for pretty much any solver. It's only here to inform you there is no MRF model nor fvModels/fvConstraints used in your case

It's perfectly fine, unless you were trying to use one of these functionalities.
These are optional features, your solver does not need it to run you case. If the solver crashes, it's related to something else.

Don't you get an error when it stops?

Regards,
Yann

thank Yann to reply,
the message above was all content when I tried to on Bash ubuntu on window
When I ran on Linux so that's the message:

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 10-c4cf895ad8fa
Exec   : buoyantFoam
Date   : Jun 23 2023
Time   : 19:02:20
Host   : "iatnart"
PID    : 7375
I/O    : uncollated
Case   : /home/iatnart/Documents/fillingRoom
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: Convergence criteria found
        p_rgh: tolerance 0.0001
        U: tolerance 0.0001
        h: tolerance 0.0001
        "(k|epsilon|omega)": tolerance 0.001


PIMPLE: Operating solver in steady-state mode with 1 outer corrector
PIMPLE: Operating solver in SIMPLE mode


Reading thermophysical properties

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

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
    model           kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}

Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity

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

No fvModels present
No fvConstraints present
#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  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::compressibleMomentumTransportModel> >, Foam::compressibleMomentumTransportModel>::F2() const at ??:?
#8  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::compressibleMomentumTransportModel> >, Foam::compressibleMomentumTransportModel>::F23() const at ??:?
#9  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::compressibleMomentumTransportModel> >, Foam::compressibleMomentumTransportModel>::correctNut() at ??:?
#10  ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/buoyantFoam"
#11  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13  ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/buoyantFoam"

[Yann]

OK so basically the solver crashes the moment it tries starting time loop.

The classic Foam::error:: printStack happens when some math get wrong (for instance it ends up trying to divide something by 0) and it doesn't clearly tell you what is wrong.
(more details here: Foam::error::PrintStack)

We can see this block in your log:

#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  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::compressibleMomentumTransportModel> >, Foam::compressibleMomentumTransportModel>::F2() const at ??:?
#8  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::compressibleMomentumTransportModel> >, Foam::compressibleMomentumTransportModel>::F23() const at ??:?
#9  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::compressibleMomentumTransportModel> >, Foam::compressibleMomentumTransportModel>::correctNut() at ??:?

This seems to indicate the problem is related to the turbulence model. So you can double check your setup to see if something related to turbulence seems wrong.

After a quick look at your boundary conditions, I see you defined inletOutlet BC for nut on all your inlets and outlet:

boundaryField
{
    frontAndBack
    {
        type            nutUWallFunction;
        value           uniform 0;
    }
    Machine_heat
    {
        type            nutUWallFunction;
        value           uniform 0;
    }
    Machine_wall
    {
        type            nutUWallFunction;
        value           uniform 0;
    }
    Walls
    {
        type            nutUWallFunction;
        value           uniform 0;
    }
    inlet1
    {
	type		inletOutlet;
	inletValue	$internalField;
	value		$internalField;
    }
    inlet2
    {
	type		inletOutlet;
	inletValue	$internalField;
	value		$internalField;
    }
    outlet1
    {
	type		inletOutlet;
	inletValue	$internalField;
	value		$internalField;
    }
    outlet2
    {
	type		inletOutlet;
	inletValue	$internalField;
	value		$internalField;
    }
}

It should be defined as calculated, because nut is calculated by the turbulence model. (have a look at the tutorials in $FOAM_SCR/heatTransfer/buoyantFoam to have some boundary conditions setup examples)

Your error might be related to the fact you impose nut=0 on your inlets and outlets. Try to replace it by calculated BC and rerun your case.
I cannot guarantee this is the cause of your error and I didn't check the rest of your case setup. If it doesn't solve the issue, repeat the process, and compare your case with some tutorials to spot differences and potential errors.

Good luck!
Yann