[SalumeH]

Hey everyone.

I'm currently doing simulations of a turbulent flux in an open channel which has a static hydraulic wheel with the solver interFoam. The 2D mesh was done using SALOME and the configuration of the wheel's axis and the depth is shown in the figure. The internal field velocity for water is uniform (0.172028 0 0) and the air should have a velocity equal to zero at all instant.

The assigned boundary conditions are for U:

boundaryField
{
    back
    {
        type            empty;
    }
    bottom
    {
        type            noSlip;
    }
    waterwheel
    {
        type            noSlip;
    }
    front
    {
        type            empty;
    }
    waterinlet
    {
        type            fixedValue;
        value           uniform (0.172028 0 0);
    }
    airinlet
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    atmosphere
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform (0 0 0);
    }
}

p_rgh

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

internalField   uniform 0;

boundaryField
{
    waterinlet
    {
        type            fixedFluxPressure;
	value		uniform 0;
    }

    airinlet
    {
        type            fixedFluxPressure;
	value		uniform 0;
    }


    outlet
    {
	type		prghPressure;
	p		uniform 0;
    }

    bottom
    {
        type            fixedFluxPressure;
	value		uniform 0;
    }

    waterwheel
    {
        type            fixedFluxPressure;
	value		uniform 0;
    }


    atmosphere
    {
        type            totalPressure;
        p0              uniform 0;
    }

    front
    {
        type            empty;
    }

    back
    {
        type            empty;
    }
}

alpha.water

boundaryField
{
    back
    {
        type            empty;
    }
    bottom
    {
        type            zeroGradient;
    }
    waterwheel
    {
        type            zeroGradient;
    }
    front
    {
        type            empty;
    }
    waterinlet
    {
        type            fixedValue;
        value           uniform 1;
    }
    airinlet
    {
        type            fixedValue;
        value           uniform 0;
    }
    atmosphere
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    outlet
    {
        type            zeroGradient;
    }
}

k

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

internalField   uniform 1.1098e-4;

boundaryField
{
    waterinlet
    {
        type            fixedValue;
        value           uniform 1.1098e-4;
    }

    airinlet
    {
        type            fixedValue;
        value           uniform 1.1098e-4;
    }

    bottom
    {
        type            kqRWallFunction;
        value           uniform 0;
    }

    waterwheel
    {
        type            kqRWallFunction;
        value           uniform 0;
    }

   outlet
    {
        type            inletOutlet;
        inletValue      uniform 1.1098e-4;
        value           uniform 1.1098e-4;
    }

    atmosphere
    {
        type            inletOutlet;
        inletValue      uniform 1.1098e-4;
        value           uniform 1.1098e-4;
    }

    front
    {
        type            empty;
    }

    back
    {
        type            empty;
    }

}

omega

internalField   uniform 5.6838;

boundaryField
{
    waterinlet
    {
        type            fixedValue;
        value           uniform 5.6838;
    }

    airinlet
    {
        type            fixedValue;
        value           uniform 5.6838;
    }

    bottom
    {
        type            omegaWallFunction;
        value           uniform 0.002;
    }

    waterwheel
    {
        type            omegaWallFunction;
        value           uniform 0.002;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      uniform 5.6838;
        value           uniform 5.6838;
    }

    atmosphere
    {
        type            inletOutlet;
        inletValue      uniform 5.6838;
        value           uniform 5.6838;
    }

    front
    {
        type            empty;
    }

    back
    {
        type            empty;
    }

}

However, I'm getting results that do not agree with the expected outcome. I've tried with different conditions and with meshes with different number of elements. After a lot of time reviewing everything, I haven't been able to figure out what is the reason behind these results.

So far, I think it may be:

1. Boundary and initial conditions.
2. The mesh.
3. Perhaps, I have to do a modification to the InterFoam code.

Could someone please, verify which of the given reasons is what causing my results? Otherwise, which is the real problem?

Thank you very much to all.