[Raza Javed]

Hello Everyone,


I made a very simple geometry in Salome, a rectangular pipe, and I want to model the fluid flow from it.


I try to explain my problem as clear as I can.


My geometry has 2 regions. (one is pipe and other is fluid). You can see the geometry below in the figure. (The yellow region is the fluid and the outer green region is pipe)



I imported this geometry to Openfoam using ideasUnvToFoam. I am using chtMultiRegionSimpleFoam. After changing all the boundary conditions using changeDictionaryDict. When I run the solver, I get the following error:

Solving for fluid region fluid
#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::laminar<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > >::nuEff() const at ??:?
#7  Foam::linearViscousStress<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > >::divDevRhoReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const at ??:?
#8  ? at ??:?
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? at ??:?
Floating point exception (core dumped)

It shows that when it starts solving the fluid region, it gives the error. The problem is that I couldn't even understand what is the error?



I am attaching my changeDictionaryDict from system/region_name for your reference:


PIPE

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      changeDictionaryDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

T
{
    internalField   uniform 300;
    pipeInlet
    {
        type            wall;

    }
    pipeOutlet
    {
        type            wall;

    }
    defaultFaces
    {
        type            wall;

    }
    boundaryField
    {
        "pipeInlet"
        {
            type            fixedValue;
            value           uniform 300;
        }
        "pipeOutlet"
        {
            type            fixedValue;
            value           uniform 300;
        }
        "defaultFaces"
        {
            type            zeroGradient;

        }

        "pipe_to_.*"
        {
            type            compressible::turbulentTemperatureCoupledBaffleMixed;
            Tnbr            T;
            kappaMethod     solidThermo;
            kappaName       none;
            value           uniform 300;
        }
    }
}

And the boundary file from constant/region_name/polymesh is given below:

FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/pipe/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

4
(
    pipeInlet
    {
        type            patch;
        nFaces          36;
        startFace       4896;
    }
    pipeOutlet
    {
        type            patch;
        nFaces          36;
        startFace       4932;
    }
    defaultFaces
    {
        type            patch;
        nFaces          1108;
        startFace       4968;
    }
    pipe_to_fluid
    {
        type            mappedWall;
        inGroups        1(wall);
        nFaces          860;
        startFace       6076;
        sampleMode      nearestPatchFace;
        sampleRegion    fluid;
        samplePatch     fluid_to_pipe;
    }
)

// ************************************************************************* //

And for fluid region, the corresponding files are given below:


FLUID
system/region_name/changeDictionaryDict

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      changeDictionaryDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
boundary
{
    fluidInlet
    {
        type            wall;
    }
    fluidOutlet
    {
        type            wall;
    }
}
U
{
    internalField   uniform (0 0 0.001);

    boundaryField
    {
        fluidInlet
        {
            type            fixedValue;
            value           uniform (0 0 0.001);
        }

        fluidOutlet
        {
            type            inletOutlet;
            inletValue      uniform (0 0 0);
        }

        ".*"
        {
            type            noSlip;
        }
    }
}

T
{
    internalField   uniform 300;

    boundaryField
    {
        fluidInlet
        {
            type            fixedValue;
            value           uniform 300;
        }

        fluidOutlet
        {
            type            inletOutlet;
            inletValue      uniform 300;
        }

        ".*"
        {
            type            zeroGradient;
            value           uniform 300;
        }

        "fluid_to.*"
        {
            type            compressible::turbulentTemperatureCoupledBaffleMixed;
            Tnbr            T;
            kappaMethod     fluidThermo;
            value           uniform 300;
        }
    }
}


epsilon
{
    internalField   uniform 0.01;

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

        maxX
        {
            type            inletOutlet;
            inletValue      uniform 0.01;
        }

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

k
{
    internalField   uniform 0.1;

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

        maxX
        {
            type            zeroGradient;
            value           uniform 0.1;
        }

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


p_rgh
{
    internalField   uniform 0;

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

        fluidOutlet
        {
            type            fixedValue;
            value           uniform 0;
        }

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

p
{
    internalField   uniform 0;

    boundaryField
    {
        ".*"
        {
            type            calculated;
            value           uniform 0;
        }
    }
}

// ************************************************************************* //

and Constant/region_name/polymesh/boundary

FoamFile
{
    version     0.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/fluid/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


3
(
    fluidInlet
    {
        type            wall;
        nFaces          32;
        startFace       5723;
    }

    fluidOutlet
    {
        type            wall;
        nFaces          30;
        startFace       5755;
    }

    fluid_to_pipe
    {
        type            mappedWall;
        inGroups        1 ( wall );
        nFaces          860;
        startFace       5785;
        sampleMode      nearestPatchFace;
        sampleRegion    pipe;
        samplePatch     pipe_to_fluid;
    }

)

I shall be thankful if someone can guide, what exactly this error is? and How can I rectify it?


Thank you