CFD Online Discussion Forums - rhoReactingBuoyantFoam_2 species

Dear fellow Foamers,

I am trying to simulate a leak from a pump (box) using rhoReactingBuoyantFoam. Propane is leaking from a hole inside the pump, as shown in the figure below. The problem is, the velocity in the figure below is 10 m/s and i am seeing some propane coming out, But when i gave 1 m/s velocity at the inlet (leak), there is no propane coming out, there is nothing going on, no velocity or anything. But in the log file of the solver, propane species equation is getting solved (i can see the iterations). I thought may be the velocity is too low, and i tried to increase the tolerance in fvSolutions but no change. I tried with 2 m/s velocity then there is some propane coming out. But why its not coming out for 1m/s. I need results for 1 m/s:confused:. Please can someone shed some light on this. Any ideas are welcome.:)
Images order: pump geometry 1,2, leak 1ms, 2ms, 10ms

Here are the case_files

fvSolutions

Code:

/*--------------------------------*- C++ -*----------------------------------*\

| =========                 |                                                 |

| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |

|  \\    /   O peration     | Version:  4.x                                   |

|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |

|    \\/     M anipulation  |                                                 |

\*---------------------------------------------------------------------------*/

FoamFile

{

    version     2.0;

    format      ascii;

    class       dictionary;

    location    "system";

    object      fvSolution;

}

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



solvers

{

    "rho.*"

    {

        solver          diagonal;

    }



    p

    {

        solver           PCG;

        preconditioner   DIC;

        tolerance        1e-6;

        relTol           0.1;

    }

    p_rgh

    {

        solver           PCG;

        preconditioner   DIC;

        tolerance        1e-6;

        relTol           0.1;

    }



    p_rghFinal

    {

        $p;

        tolerance        1e-6;

        relTol           0.0;

    }

    pFinal

    {

        $p;

        tolerance        1e-6;

        relTol           0.0;

    }



    "(U|h|k|omega|epsilon)"

    {

        solver          PBiCG;

        preconditioner  DILU;

        tolerance       1e-10;

        relTol          0.1;

    }



    "(U|h|k|omega|epsilon)Final"

    {

        $U;

        relTol          0;

    }



    "(Yi|YiFinal)"

    {

        solver          PBiCG;

        preconditioner  DILU;

        tolerance       1e-10;

        relTol          0;

        //$hFinal;

    }



}



PIMPLE

{

    momentumPredictor no;

    nOuterCorrectors  1; //2;

    nUCorrectors     2;

    nCorrectors    1; //1;

    nNonOrthogonalCorrectors 0;

}



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

cheers,
nikhil.