[jcw]

Hello!
After reading many threads in forum, I am still not sure, which OF solver to use.
I would like to simulate the mixing of 2 gases (air and CH4). Each gas component is entering a separate inlet. At the outlet I would like to evaluate the mixing quality. There is no phase change. Gases enter with different temperatures.
Which OF solver is appropriate? Is there a tutorial? I use OF4.1.
Thanks.

[jcw]

No hints for my question?

[tomf]

Hi,

You could have a look at reactingFoam and it's variants. If you turn off chemistry/combustion you will get just the mixing.

I would suggest looking at either of:

$FOAM_TUTORIALS/combustion/reactingFoam/laminar/counterFlowFlame2D
$FOAM_TUTORIALS/combustion/reactingFoam/laminar/counterFlowFlame2DLTS

depending on whether you want to look at time-accurate or steady-state simulations.

Regards,
Tom

[jcw]

Hello!
Thanks for reply. I actually managed with simpleReactingParcelFoam. Lagrange off, non-reacting.
Regards

[CFDUser29]

Originally Posted by jcw

Hello!
Thanks for reply. I actually managed with simpleReactingParcelFoam. Lagrange off, non-reacting.
Regards

Can you please share you constant and system directories please? I struggle to set a simulation with two gases.

[jcw]

Originally Posted by CFDUser29

Can you please share you constant and system directories please? I struggle to set a simulation with two gases.

Hello! It is a long time ago I used OF. My simulation is done with OF41. Maybe it will help.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \    /   O peration     | Version:  4.0                                   |
|   \  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
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(phid,p)     bounded Gauss upwind;
    div(phi,K)      bounded Gauss linear;
    div(phi,h)      bounded Gauss upwind;
    div(phi,k)      bounded Gauss upwind;
    div(phi,epsilon) bounded Gauss upwind;
    div(phi,omega) bounded Gauss upwind;
    div(phi,Yi_h)   Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear uncorrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         uncorrected;
}

wallDist
{
    method meshWave;
}

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \    /   O peration     | Version:  4.0                                   |
|   \  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver           GAMG;

        tolerance        0;
        relTol           0.05;

        smoother         DICGaussSeidel;



        maxIter          50;
    };

    Phi
    {
        $p;
    }

    "(U|Yi|h|k|omega|epsilon)"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        tolerance       0;
        relTol          0.1;
        maxIter         20;
    }
}

potentialFlow
{
    nNonOrthogonalCorrectors 5;
}

SIMPLE
{
    nCorrectors     2;
    nNonOrthogonalCorrectors 1;
    rhoMin          0.1;
    rhoMax          1.5;
}

relaxationFactors
{
    fields
    {
        p               0.3;
        rho             1;
    }
    equations
    {
        U               0.7;
        h               0.7;
        ".*"            0.7;
    }
}



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