[Raza Javed]

Hello Everyone,


I am using chtMultiRegionSimpleFoam and my Openfoam version is 4.1.


I have created geometry in Salome and and then imported the UNV file to Openfoam.


My geometry consists of a pipe and three rectangular plates, and all these things are in one box.



The rectangular plates are acting as a heat source which I made using fvOption.


Now, the fluid will flow from the pipe, and due to these hot rectangular plates, the temperature of the fluid will change inside the pipe.


At the inlet, I need to put the flow-rate of 5 litres/minutes.



At the inltet, temperature could be room temperature.


Now, My task is to see the temperature change at the outlet.



I have put some boundary conditions for Velocity, Temperature and Pressure, at inlet, outlet and at the walls of pipe.


But my solver is not converging.



I am not getting which boundary conditions to use.


I have tried putting flowRateInletVelocity, but I don't know which boundary conditions to put at other boundaries in this case?


I am putting my geometry and boundary conditions with this post.


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


Thank you

boundary
{
    inlet
    {
        type            patch;
    }
    outlet
    {
        type            patch;
    }

}
U
{
    internalField   uniform (0 1e-3 0);

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

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

T
{
    internalField   uniform 300;

    boundaryField
    {
        inlet
        {
            type            fixedValue;
            value           uniform 450;//$internalField;
            
        }

        outlet
        {
            type            inletOutlet;
            value           $internalField;
            inletValue      $internalField;
        }

        "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            fixedFluxPressure;//zeroGradient;
            value           uniform 0;
        }

        outlet
        {
            type            fixedValue;
            value           uniform 0;
        }

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

p
{
    internalField   uniform 0;

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