[decah]

Hi all

I am relatively inexperienced with Lagrangian simulations and am noticing that when running DPMFoam with the nonSphereDrag particle force selected no particles enter the simulation. The simulation just runs without particles.

I am injecting particles into a tank with water whose velocity field is solved in simpleFoam then mapped with mapFields. Particles are successfully injected when ErgunWenYuDrag or sphereDrag are chosen but I would like to simulate non-spherical particles. Here is my kinematicCloudProperties file:

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

solution
{
    active          true;
    coupled         true; 
    transient       yes; 
    cellValueSourceCorrection off;

    interpolationSchemes
    {
        rho.water         cell;
        U.water           cellPoint;
        mu.water          cell;
    }

    integrationSchemes
    {
        U               Euler;
    }

    sourceTerms
    {
        schemes
        {
            U semiImplicit 1;
        }
    }
}

constantProperties
{
    parcelTypeId 1;

    rhoMin          1e-15;
    minParcelMass   1e-15;

    rho0            2650;
    youngsModulus   6e8;  
    poissonsRatio   0.4; 
    constantVolume  false;

    alphaMax        0.9;
}

subModels
{
    particleForces
    {
       nonSphereDrag
	{
            alphac alpha.water;
            phi 0.75;
	};
      gravity;
    }

    injectionModels
    {
        
    
  model1
        {
            type            patchInjection;
            massTotal       0.375;
            SOI             0;
            parcelBasisType mass;
            patchName       inlet;
            duration        1;
            parcelsPerSecond 1000;
            U0              (0.025 0 0);
            flowRateProfile constant 1;
         
  	    sizeDistribution 
            {
             type        RosinRammler;
	    	RosinRammlerDistribution 
		{
		minValue 2e-6; 
		maxValue 100e-6;
	    	d 40e-6; 
	    	n 3; 
	       }        
            } 
	}    

}
dispersionModel none;

    patchInteractionModel localInteraction;

    localInteractionCoeffs
    {
        patches
        (
            fixedWalls
            {
                type rebound;
                e    0.1; 
                mu   0.5; 
            }
            bottom
            {
                type stick;
            }
            atmosphere
            {
                type rebound;
                e    0.1;
                mu   0.1;
            }
            outlet
            {
                type escape;
            }
	    inlet
            {
                type escape;
            }
        );
    }

    StandardWallInteractionCoeffs
    {
        type rebound;
        e    0.85;
        mu   0.1;
    }

    heatTransferModel none;

    surfaceFilmModel none;

    collisionModel none;

    stochasticCollisionModel none;

    radiation off;
}


cloudFunctions
{}


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

If anyone has any ideas or has used nonSphereDrag before any help would be appreciated!


Declan