[flowwwww]

Hello everyone,


I am stuggling with my simulation and i would like to know if one of you can give me a hand.
I have a gemetry that consist in two ducts that inyect air to a chamber and i would like to study the distribution of temperature inside.
i have started with a more simply geometry that just have one inlet duct. In this case, the simulation ran without any problem and i obtained some reasonable results.
inyector simple 1.png un inyector.png


However, when i add the other duct to the geometry and i tried to simulate with se same boundary conditions to the second inlet (i created a new patch), the simulation can not end. The results stops and i get the next.
inyectordoble1.png inyector doble2.png
The only change is in the domain, that is adapted to the size of the two inlet ducts



If someone can help me i really appreciate it.

Best Regards.

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 1e5;

boundaryField
{
    door  //domain
    {
        type            fixedValue;
        value           $internalField;
    }
    walls  //geometry
    {
        type            zeroGradient;
    }
     box  //domain
    {
        type            zeroGradient;
    }
    inlet  //geometry
    {
        type            zeroGradient;
    }
    inlet1  //geometry
    {
        type            zeroGradient;
    }

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 450;

boundaryField
{

    door
    {
        type            inletOutlet;
        value           uniform 450;
        inletValue      uniform 450;
    }
    walls
    {
        type            zeroGradient;
    }
 box
    {
        type            zeroGradient;
    }

    inlet
    {
        type            fixedValue;
        value           uniform 610;
    }
    inlet1
    {
        type            fixedValue;
        value           uniform 635;
    }

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
     door
    {
           type            inletOutlet;
        inletValue      uniform (0 0 0);
        value           $internalField;
    }    
    walls
    {
        type            noSlip;
    }
     box
    {
        type            noSlip;
    }
    inlet
    {
        type            fixedValue;
        value           uniform (0 -6.49 0);
    }
    inlet1
    {
        type            fixedValue;
        value           uniform (0 -6.79 0);
       
    }

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default             steadyState;
}

gradSchemes
{
    default             Gauss linear;
}

divSchemes
{
    default             none;

    div(phi,U)          bounded Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U)))))      Gauss linear;
    div(phi,e)          bounded Gauss upwind;
    div(phi,epsilon)    bounded Gauss upwind;
    div(phi,k)          bounded Gauss upwind;

    div(phid,p)         Gauss upwind;
    div(phi,Ekp)        bounded Gauss upwind;
    div((phi|interpolate(rho)),p)  Gauss upwind;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}


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

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-08;
        relTol          0.1;
        smoother        GaussSeidel;
        nCellsInCoarsestLevel 20;
    }

    Phi
    {
        $p;
    }

    "(U|e|k|epsilon)"
    {
        solver          GAMG;
        tolerance       1e-08;
        relTol          0.1;
        smoother        GaussSeidel;
        nCellsInCoarsestLevel 20;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 2;
    pMinFactor      0.1;
    pMaxFactor      2;
    transonic       no;
    consistent      yes;
    pRefCell 0; 
    pRefValue 1e5;

    residualControl
    {
        p               1e-3;
        U               1e-4;
        e               1e-3;

        // possibly check turbulence fields
        "(k|epsilon|omega)" 1e-3;
    }
}

potentialFlow
{
    nNonOrthogonalCorrectors 10;
}


relaxationFactors
{
    fields
    {
        p               0.8;
    }
    equations
    {
        p               0.8;
        U               0.9;
        e               0.8;
        k               0.9;
        epsilon         0.9;
    }
}

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