CFD Online Discussion Forums - Uncertain k values with multiphaseStabilizedTurbulence

Hi!

I have uncertain k values from my multiphase simulations.
I used the new multiphaseStabilizedTurbulence fvOption to get correct k everywhere. Unfortunatelly, I did not get the k values I expected. On my graphs, one can see average k in specified verticals. The green line shows the result I expect, but If I refine my mesh, the k values are not better, but worse. Why is it? Do I have to use different BCs, or is there a mistake in the fvOptions?
I have a RAS with k-omega SST turbulence closure, and I use kqRWallFunction on the walls. For free surface capturing, I use isoAdvector (interIsoFoam).

Thank you in advance!

Code:

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

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

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

|  \\    /   O peration     | Version:  v1812                                 |

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

|    \\/     M anipulation  |                                                 |

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

FoamFile

{

    version     2.0;

    format      ascii;

    class       dictionary;

    location    "constant";

    object      fvOptions;

}

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



    multiphaseStabilizedTurbulence1

    {

        type            multiphaseStabilizedTurbulence;

        active          yes;



        multiphaseStabilizedTurbulenceCoeffs

        {

            // Optional coefficients

            lambda2         0.1;   // A value of 0 sets the nut correction to 0

            Cmu             0.09;  // from k-epsilon model

            C               1.51;  // from k-omega model

            alpha           1.36;  // 1/Prt

        }

    }





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