[yb8119]

Greetings

I am new to OpenFOAM and I am simulating the flow through canopies.

The version of the OpenFOAM I am using is V3.0.1, and the solver I use is pisoFoam.

The B.C. sets
inlet|outlet|front|back U cyclicAMI p cyclicAMI
lowerWall U fixed 0 p zeroGradient
upperWall U symmetry p symmetry

(lowerWall -- red
outlet -- green
back -- blue
in the attachment)

I put all cyclic boundaries inside one processor by using this in decomposeParDict

method          scotch;
preservePatches (inlet_ outlet_ front_ back_);

and to create inlet flow, I put this in the fvOptions

momentumSource1
{
    type            vectorSemiImplicitSource;
    active          yes;

    vectorSemiImplicitSourceCoeffs
    {
        selectionMode   all;
        volumeMode      absolute;
        injectionRateSuSp
        {
            U           ((1.34 0 0) 0);
        }
    }
}

Then I ran the simulation and misconvergence occurred within several hundred steps.

I am not sure what caused this problem. I tried to use this to create inlet flow

momentumSource
{
    type            meanVelocityForce;
    active          yes;

    meanVelocityForceCoeffs
    {
        selectionMode   all;

        fieldNames      (U);
        Ubar            (1.34 0 0);
    }
}

But it did not help. I also tried to switch the B.C. of the upperWall to slip, not helpful either.

I am wondering what I did wrong on the configuration of the simulation, any suggestion or idea will be appreciated!

Thanks in advance!

Best Regards
Ben