CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Mass loading - InjectionModel DPMFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 14, 2018, 10:28
Default Mass loading - InjectionModel DPMFoam
  #1
New Member
 
Join Date: Sep 2018
Posts: 4
Rep Power: 7
scrubber is on a distinguished road
Dear all.

I have a question related to mass flow rate for injection of particles, using DPMsolver. I aim to have a mass flow rate of 1,58*10^(-5) kg/s.

This is my strategy:
Number of particles = 10 000
Volume = (3/4)*pi*(d/2)^3 where d = 70*10^(-6) m
Particle density = 8800 kg/m^3
Duration of simulation 1s.

Mass loading = N_particles*Volume*Density/ Time_simulation = 1,58*10^(-5) kg/s

Can someone verify my code is correct?
See attached code for kinematicCloudProperties.

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

solution
{
    active          true;
    coupled         true;
    transient       yes;
    cellValueSourceCorrection on; // org: off

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

    integrationSchemes
    {
        U               Euler;
    }

    sourceTerms
    {
        schemes
        {
            U semiImplicit 1;
        }
    }
}

constantProperties
{
    parcelTypeId 1;

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

    rho0            8800;
    youngsModulus   117e6;
    poissonsRatio   0.355;

    constantVolume  false;

    alphaMax        0.99;
}

subModels
{
    particleForces
    {
        sphereDrag//WenYuDrag
        {
            alphac alpha.air;
        }
        gravity;
        }

//////
    injectionModels
    {
	model1
	{
            type             patchInjection;
            patchName        inlet;
            duration         1;
	    parcelsPerSecond 10000;
            massTotal        0;
            parcelBasisType  fixed;
            flowRateProfile  constant 1;
            nParticle        1;
            SOI              0.0;
            U0               (9.39 0 0);
            sizeDistribution
            {
                type        fixedValue;
                fixedValueDistribution
                {
                    value   0.00007;
                }
            }
        }
    }
/////



    dispersionModel none;

    patchInteractionModel localInteraction;

    localInteractionCoeffs
    {
        patches
        (
	    inlet
            {
                type escape;
                e    0.97;
                mu   0.09;
            }    

	    upperWall
            {
                type rebound;
                e    0.97;
                mu   0.09;
            }
            lowerWall
            {
                type rebound;
                e    0.97;
                mu   0.09;
            }
            outlet
            {
                type escape;
                e    0.97;
                mu   0.09;
            }
            frontAndBack
            {
                type rebound;
                e    0.97;
                mu   0.09;
            }
        );
    }

    StandardWallInteractionCoeffs
    {
        type rebound;
        e    0.97;
        mu   0.09;
    }

    heatTransferModel none;

    surfaceFilmModel none;

    collisionModel none;//pairCollision;

    pairCollisionCoeffs
    {
        maxInteractionDistance  0.0025;

        writeReferredParticleCloud no;

        pairModel pairSpringSliderDashpot;

        pairSpringSliderDashpotCoeffs
        {
            useEquivalentSize   no;
            alpha               0.02;
            b                   1.5;
            mu                  0.10;
            cohesionEnergyDensity 0;
            collisionResolutionSteps 12;
        };

        wallModel wallSpringSliderDashpot;

        wallSpringSliderDashpotCoeffs
        {
            useEquivalentSize no;
            collisionResolutionSteps 12;
            youngsModulus   1e8;
            poissonsRatio   0.23;
            alpha           0.01;
            b               1.5;
            mu              0.09;
            cohesionEnergyDensity 0;
        };

        U     U.air;
    }

    stochasticCollisionModel none;

    radiation off;
}


cloudFunctions
{}


// ************************************************************************* //
scrubber is offline   Reply With Quote

Reply

Tags
dpmfoam, mass flow rate

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
mass fraction of species Lilly FLUENT 5 March 13, 2022 18:52
DPMFoam diverged when inject particles with high "nParticle" bijan darbari OpenFOAM Running, Solving & CFD 0 July 30, 2016 04:03
Exit Corrected Mass Flow Rate Mesh Sensitivity Study s__s__s CFX 4 July 20, 2016 12:46
Net mass flow inlet vs outlet Nigui28 FLUENT 1 August 12, 2011 11:09
Mass flux and mass flow rate us Phoenics 10 March 31, 2006 03:06


All times are GMT -4. The time now is 17:16.