[kumar.kyutech]

I am trying to simulate capillary flow in a 3D open rectangular microchannel using interFoam. I modified the capillary rise tutorial example to make it for an open channel.

After I run the simulation, the results don't look correct. The fluid is not moving and the free surface is not visible.

Nomenclature:
wallss - open surface
wall - microchannel walls

open.JPG

Can anyone help identify the problem?

Here is my source code:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type inletOutlet;
value uniform 1;
inletValue uniform 1;
}

atmosphere
{
type zeroGradient;
}

walls
{
type constantAlphaContactAngle;
theta0 30;
limit gradient;
value uniform 0;
}

wallss
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

defaultFaces
{
type empty;
}
}


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

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{
inlet
{
type fixedValue;
value uniform 0;
}

atmosphere
{
type fixedValue;
value uniform 0;
}

walls
{
type fixedFluxPressure;
}

wallss
{
type totalPressure;
p0 uniform 0;
}

defaultFaces
{
type empty;
}
}


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

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}

atmosphere
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}

walls
{
type noSlip;
}

wallss
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}

defaultFaces
{
type empty;
}
}


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