[Raza Javed]

Hello Everyone,


I am using chtMultiRegionSimpleFoam and my OpenFoam version is 4.1.


My geometry has a box, and inside this box, there is a fluid region.
The inlet and outlet of the fluid region is at two faces of the box.


I have to specify a flow rate at the inlet. And according to that I set my boundary conditions like below:

boundary
{
    inlet
    {
        type            patch;
    }
    outlet
    {
        type            patch;
    }
    defaultFaces
    {
        type            patch;
    }
}
U
{
    internalField   uniform (1e-3 0 0);

    boundaryField
    {
        inlet
        {
            type                flowRateInletVelocity;//fixedValue;//pressureInletVelocity;
            volumetricFlowRate  0.066;
            extrapolateProfile  yes;
            value               uniform (1 0 0);
        }

        outlet
        {
            type            zeroGradient;
            //value           uniform (0 0 0);//$internalField;
        }
        defaultFaces
        {
            type            zeroGradient;
        }
        "fluid_to_box"
        {
            type            noSlip;
        }
    }
}

T
{
    internalField   uniform 300;

    boundaryField
    {
        inlet
        {
            type            fixedValue;
            value           $internalField;
            
        }

        outlet
        {
            type            zeroGradient;//inletOutlet;
            value           $internalField;
            //inletValue      $internalField;
        }
        defaultFaces
        {
            type            zeroGradient;
        }

        "fluid_to_box"
        {
            type            compressible::turbulentTemperatureCoupledBaffleMixed;
            Tnbr            T;
            kappaMethod     fluidThermo;
            value           uniform 300;
        }
    }
}


epsilon
{
    internalField   uniform 0.01;

    boundaryField
    {
        inlet
        {
            type            fixedValue;
            value           uniform 0.01;
        }

        outlet
        {
            type            inletOutlet;
            inletValue      uniform 0.01;
        }

        ".*"
        {
            type            epsilonWallFunction;
            value           uniform 0.01;
        }
    }
}

k
{
    internalField   uniform 0.1;

    boundaryField
    {
        inlet
        {
            type            inletOutlet;
            inletValue      uniform 0.1;
        }

        outlet
        {
            type            zeroGradient;
            value           uniform 0.1;
        }

        ".*"
        {
            type            kqRWallFunction;
            value           uniform 0.1;
        }
    }
}


p_rgh
{
    internalField   uniform 0;

    boundaryField
    {
        inlet
        {
            type            zeroGradient;
            value           uniform 0;
        }

        outlet
        {
            type            fixedValue;
            value           uniform 0;
        }

        ".*"
        {
            type            fixedFluxPressure;
            value           uniform 0;
        }
    }
}

p
{
    internalField   uniform 0;

    boundaryField
    {
        ".*"
        {
            type            zeroGradient;
            //value           uniform 0;
        }
    }
}

But when I RUN this solver it gives me the error that "maximum number of iterations exceeded". The error is given below:

--> FOAM FATAL ERROR: 
Maximum number of iterations exceeded

    From function Foam::scalar Foam::species::thermo<Thermo,  Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar  (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar,  Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo,  Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar  (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const)  const [with Thermo =  Foam::hConstThermo<Foam::rhoConst<Foam::specie> >; Type =  Foam::sensibleEnthalpy; Foam::scalar = double;  Foam::species::thermo<Thermo, Type> =  Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie>  >, Foam::sensibleEnthalpy>]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::heRhoThermo<Foam::rhoThermo,  Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie>  >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#3  Foam::heRhoThermo<Foam::rhoThermo,  Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie>  >, Foam::sensibleEnthalpy> > > >::correct() at ??:?
#4  ? at ??:?
#5  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6  ? at ??:?
Aborted (core dumped)

I think there is some problem with my boundary conditions of the fluid.

Log file of solver is also attached.

Can someone please guide me where I am doing mistake?


I shall be very thankful if someone can help me out in this.


Thank you

[AbdulazizAlkandari]

Did you manage to fix the problem? If so can you explain please.

[priyankap]

Originally Posted by AbdulazizAlkandari

Did you manage to fix the problem? If so can you explain please.

First, please RUN checkMesh to check if you have any non-orthogonal faces in your mesh.


For me, I was getting this error due to my mesh because I had non-orthogonal faces in my mesh, I refined my mesh to remove the non-orthogonality and then I didn't get this error.


If you can't refine your mesh, then you should add nonOrthoCorrectors in your fvSolution file. OR, you can also modify your fvSchemes to compensate the non-orthogonality.


You can check the link below about the modification of fvSchemes.


https://cfd.direct/openfoam/user-guide/v6-fvschemes/


Hope this helps.


Best,


Priyanka