[obiscolly50]

Hello all,
I modified a simplefoam solver to include temperature, porous media and lagrange particle tracking. I'm trying to solve the flow for a steady-state case in parallel, but i'm experiencing a strange problem. I'm using the flow field from an already converged case as the initial 0 field. After introducing the particles at the inlet patch they remain roughly in the same position for subsequent iterations.
Below is the kinematicCloudProperties file i'm using which was derived from the steady-state particle solver simpleReactingParcelFoam particle properties file.
I'm supsecting two issues:-
1) The particles are stuck just before the first porous zone. I implemented the porosity as a momentum sink in the UEqn not as a physical barrier and didn't think this would be an issue since the particle should take the path of least resistance.
2) I'm setting something wrong in my kinematicCloudProperties file. I've not used steady-state particle tracking before and just used the file fro simpleReactingParcelFoam with very little modification.
Anyone with any idea of what could be wrong. I'm using the openFoam v6.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      kinematicCloudProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solution
{
    active          true;
    coupled         false;
    transient       no;
    
    calcFrequency   10;  //added
    maxTrackTime    5; //added
    cellValueSourceCorrection off;
    maxCo           0.3; 

    interpolationSchemes
    {
        rho             cell;
        U               cellPoint;
        mu              cell;
        DUcDt           cell;
    }

    integrationSchemes
    {
        U               Euler;
	T               analytical;
    }

    sourceTerms
    {
	resetOnStartup  no; //added
        schemes
        {
            U           semiImplicit 1;
        }
    }
}

constantProperties
{
    rho0            1200; //was 8800
    youngsModulus   1.3e5;
    poissonsRatio   0.35;
}

subModels
{
    particleForces
    {
        sphereDrag;
        
        gravity;
        
        pressureGradient
        {
            U U;
        }

    }

    injectionModels
    {
	model1
	{
            type             patchInjection;
	    massFlowRate     0.8e-04;
            parcelBasisType  mass;
            patchName        influent;
	    parcelsPerSecond 1000;
            duration         1;  //set to 1 for steady state
            U0               (-60 0 0);
            flowRateProfile  constant 1;
            //nParticle        1;
            //SOI              10;

            sizeDistribution
            {
                type        fixedValue;
                fixedValueDistribution
                {
                    value   15e-6;
                }
            }
        }
    }

    dispersionModel none;

    patchInteractionModel  localInteraction; //was standardWallInteraction

    localInteractionCoeffs
    {
        patches
        (
            "effluent|soldischarge"
            {
                type escape;
            }
            ".*"
            {
                type rebound;
                e    0.97;
                mu   0.09;
            }
        );
    }

    surfaceFilmModel none;

    stochasticCollisionModel none;
    
    collisionModel none;    

    pairCollisionCoeffs
    {
        maxInteractionDistance  0.00007;

        writeReferredParticleCloud no;

        pairModel pairSpringSliderDashpot;

        pairSpringSliderDashpotCoeffs
        {
            useEquivalentSize   no;
            alpha               0.12;
            b                   1.5;
            mu                  0.52;
            cohesionEnergyDensity 0;
            collisionResolutionSteps 12;
        };
        
        wallModel wallSpringSliderDashpot;

        wallSpringSliderDashpotCoeffs
        {
            useEquivalentSize no;
            collisionResolutionSteps 12;
            youngsModulus   1e10;
            poissonsRatio   0.23;
            alpha           0.12;
            b               1.5;
            mu              0.43;
            cohesionEnergyDensity 0;
        };
    }
}

cloudFunctions
{
    particleTracks1
    {
        type            particleTracks;
        trackInterval   5;
        maxSamples      1000000;
        resetOnWrite    yes;
    }
}

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

Below is the picture of the particles stuck at the porous zone.