[GrivalszkiP]

Hi Foamers,

Im trying to simulate a free surface river connection using interIsoFoam.

• Geometry: fictitious from the internet. Two inlet, one outlet. Quite complex.
• Mesh: with snappyHexMesh, checkMesh shows that everything is OK.
• Boundary: separated inlet for water and air phases. Values are also fictitious. variableHeightFlowRateInletVelocity have been tried, failure.
• Turbulence model: k-epsilon.

Unfortunately, it crashes after ~0.01 s. Where is my mistake?

Boundaries:

U

    in1_w
    {
        type            	fixedValue;	
		value				uniform (0.5 0 0);

    }

    in1_a
    {
        type            	fixedValue;
		value				uniform (0.5 0 0);

    }

    in2_w
    {
        type            	fixedValue;
		value				uniform (0 -0.6 0);

    }

    in2_a
    {
        type            	fixedValue;
		value				uniform (0 -0.6 0);

    }

    out
    {
        type            zeroGradient;

    }
	
	wall
	{
        type            fixedValue;
		value			uniform (0 0 0);
	}

	RiversComplete
	{
        type            fixedValue;
		value			uniform (0 0 0);
	}

	atm
	{
        type            slip;
	}

p_rgh

    atm
    {
        type            totalPressure;
        p0              uniform 0;
    }

    ".*"
    {
        type            fixedFluxPressure;
        value           uniform 0;
    }

alpha

    in1_w
    {
        type            fixedValue;
		value           uniform 1; 
    }
	
    in1_a
    {
        type            fixedValue;
		value           uniform 0; 
    }

    in2_w
    {
        type            fixedValue;
		value           uniform 1; 
    }
	
    in2_a
    {
        type            fixedValue;
		value           uniform 0; 
    }

    out
    {
        type            zeroGradient;
		value			uniform 0;
    }

    atm
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    wall
    {
        type            zeroGradient;
    }

    RiversComplete
    {
        type            zeroGradient;
    }

k

    in1_w
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in1_a
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in2_w
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in2_a
    {
        type            fixedValue;
        value           $internalField;
    }

    out
    {
        type            zeroGradient;
    }
	
	wall
    {
        type            kqRWallFunction;
        value           $internalField;
    }

    atm
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

	RiversComplete
    {
        type            kqRWallFunction;
        value           $internalField;
    }

epsilon

    in1_w
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in1_a
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in2_w
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in2_a
    {
        type            fixedValue;
        value           $internalField;
    }
	
    out
    {
        type            zeroGradient;
    }

    wall
    {
        type            epsilonWallFunction;
        value           $internalField;
    }

    atm
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    RiversComplete
    {
        type            epsilonWallFunction;
        value           $internalField;
    }

nut

    in1_w
    {
        type            calculated;
		value			uniform 0;
    }
	
    in1_a
    {
        type            calculated;
		value			uniform 0;
    }
	
    in2_w
    {
        type            calculated;
		value			uniform 0;
    }
	
    in2_a
    {
        type            calculated;
		value			uniform 0;
    }

    wall
    {
        type            nutkWallFunction;
        value           uniform 0;
    }
	
	RiversComplete
    {
        type            nutkWallFunction;
        value           uniform 0;
    }


    out
    {
        type            zeroGradient;
    }
	
    atm
    {
        type            calculated;
		value			uniform 0;
    }