[MoeCFD]

Hello All,

I have to admit that my level in OpenFoam is basic. I am going through the learning process. So, without further ado, I am trying to run simpleFoam on a fibrous material (think about it as a porous media). I created the mesh using SnappyHexMesh and and an stl file that I have. It is 3D problem with one inlet and one outlet, below are the boundary conditions that I am using:

For pressure:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{

    fixedWalls_patch20083
    {
        type            zeroGradient;
    }

    left
    {
        type            fixedValue;
		value			uniform 1e-3;
    }

    right
    {
        type            fixedValue;
		value			uniform 0;
    }

    bottom
    {
        type            zeroGradient;
    }

    top
    {
        type            zeroGradient;
    }

    frontAndBack
    {
        type            zeroGradient;
    }
}

For velocity

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (0 0 0);

boundaryField
{

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

    left
    {
        type            zeroGradient;
    }
    right
    {
        type            zeroGradient;
    }
    bottom
    {
        type            fixedValue;
		value			uniform (0 0 0);
    }
    top
    {
        type            fixedValue;
		value			uniform (0 0 0);
    }
    frontAndBack
    {
        type            fixedValue;
		value			uniform (0 0 0);
    }
}

Here is my fvSolution:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers 
{

    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-07;
        relTol          0.1;
    }

    U
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-07;
        relTol          0;
    }

}



/*
    U
    {
        solver           smoothSolver;
        smoother         GaussSeidel;
        tolerance        1e-12;
        relTol           0;
        nSweeps          1;
    }

    p
    {
        solver           GAMG;
        tolerance        1e-9;
        relTol           0.1;
        smoother         GaussSeidel;
        nPreSweeps       0;
        nPostSweeps      2;
        cacheAgglomeration on;
        agglomerator     faceAreaPair;
        nCellsInCoarsestLevel 10;
        mergeLevels      1;
    }
*/



SIMPLE
{
    nNonOrthogonalCorrectors 0;

    residualControl
    {
        p               1e-6;
        U               1e-6;
    }

	pRefCell 0;
	pRefValue 0;
}

relaxationFactors
{
    p               0.15;
    U               0.9;
}



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

When I run the command SimpleFoam to start the simulation based on the inputs shown above, I get the following (see below)
which doesn't show me what went wrong with the simulation. Is it something with the way I defined the boundary conditions or simulation settings? I appreciate your help in advance. Also, I am running the simulation in serial for now. However, I created the Mesh by running SnappyHexMesh in parallel then used reconstruct to put it together.

Build  : 4.1
Exec   : simpleFoam
Date   : Jul 07 2017
Time   : 21:57:43
Host   : 
PID    : 592
Case   : 
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: convergence criteria
    field p      tolerance 1e-06
    field U      tolerance 1e-06

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
No MRF models present

No finite volume options present


Starting time loop

Time = 1

DILUPBiCG:  Solving for Ux, Initial residual = 0.219962, Final residual = 2.60365e-08, No Iterations 10
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 4.89906e-08, No Iterations 11
DILUPBiCG:  Solving for Uz, Initial residual = 0.143777, Final residual = 9.97743e-08, No Iterations 9
moe@DESKTOP:/example$