View Issue Details

IDProjectCategoryView StatusLast Update
0003228OpenFOAMBugpublic2019-04-30 19:28
Reporterrmn_990Assigned Tohenry 
PrioritynormalSeverityminorReproducibilityhave not tried
Status closedResolutionunable to reproduce 
Product Version 
Fixed in Version 
Summary0003228: MPPICFoam/DPMFoam Wrong number of parcels on boundaries
DescriptionI did a simulation with DPMFoam and I considered different wall condition for each boundary( escape, rebound, etc.).
The illogical thing is that I just inject 550,000 particles in 5 seconds but after running, I can see that more than injected particles, I have some particles which go outside of the system or stick on one of the boundaries.
This is the log file and it can show the distribution of particles on each boundary after 314s:


[QUOTE]Courant Number mean: 0.111537 max: 0.968023
Interface Courant Number mean: 0 max: 0
deltaT = 0.00666667
Time = 314.193

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
GAMG: Solving for kinematicCloud:alpha, Initial residual = 0.0014182, Final residual = 0.00109688, No Iterations 1000
Cloud: kinematicCloud
    Current number of parcels = 287167
    Current mass in system = 1.57878e-10
    Linear momentum = (2.97358e-15 -1.83809e-15 1.12724e-14)
   |Linear momentum| = 1.1802e-14
    Linear kinetic energy = 6.06593e-16
    Injector model1:
      - parcels added = 549986
      - mass introduced = 3.02371e-10
    Parcel fate: system (number, mass)
      - escape = 262819, 1.44492e-10
    Parcel fate: patch WALL_IN (number, mass)
      - escape = 0, 0
      - stick = 0, 0
    Parcel fate: patch WALL_OUT (number, mass)
      - escape = 0, 0
      - stick = 0, 0
    Parcel fate: patch LUMEN_WALL (number, mass)
      - escape = 0, 0
      - stick = 944675672, 5.19376e-07
    Parcel fate: patch OUTLET (number, mass)
      - escape = 262819, 1.44483e-10
      - stick = 0, 0
    Parcel fate: patch INLET (number, mass)
      - escape = 0, 0
      - stick = 0, 0
    Min cell volume fraction = 0
    Max cell volume fraction = 2.13323e-05
    Min dense number of parcels = 5.5605

Continous phase-1 volume fraction = 1 Min(alphac) = 0.999979 Max(alphac) = 1
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0.013389, Final residual = 0.0114096, No Iterations 100
Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1
smoothSolver: Solving for alpha.water, Initial residual = 0.0115916, Final residual = 0.0115598, No Iterations 100
Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1
DICPCG: Solving for p_rgh, Initial residual = 2.09388e-06, Final residual = 9.41566e-08, No Iterations 4
time step continuity errors : sum local = 5.14021e-11, global = -8.3902e-12, cumulative = -6.26016e-09
DICPCG: Solving for p_rgh, Initial residual = 3.24005e-06, Final residual = 1.60079e-07, No Iterations 7
time step continuity errors : sum local = 8.73903e-11, global = -1.33005e-11, cumulative = -6.27346e-09
DICPCG: Solving for p_rgh, Initial residual = 1.38674e-06, Final residual = 9.51116e-09, No Iterations 442
time step continuity errors : sum local = 5.19234e-12, global = 3.77542e-15, cumulative = -6.27346e-09
ExecutionTime = 40261.9 s ClockTime = 40514 s

Courant Number mean: 0.111537 max: 0.968022
Interface Courant Number mean: 0 max: 0
deltaT = 0.00666667
Time = 314.2[/QUOTE]

This is KinematicCloudProperties file:
[QUOTE]/*--------------------------------*- C++ -*----------------------------------*\
  ========= |
  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
   \\ / O peration | Website: [url]https://openfoam.org[/url]
    \\ / 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 off;

    maxCo 1.0;

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

    averagingMethod dual;

    integrationSchemes
    {
        U Euler;
    }

    sourceTerms
    {
        schemes
        {
            U semiImplicit 0.00078;
        }
    }
}

constantProperties
{
    rho0 1050;
    alphaMax 0.9;
}

subModels
{
    particleForces
    {
        ErgunWenYuDrag
        {
            alphac alphac;
        }
        gravity;
    }

    injectionModels
    {
        model1
        {
            type patchInjection;
            parcelBasisType fixed;
            SOI 0;
            patch INLET;
            duration 5;
        nParticle 1;
            parcelsPerSecond 110000;
        massTotal 2.84e-11;
            U0 (0 0 0);
            flowRateProfile constant 2.5e-8;
            sizeDistribution
            {
        type fixedValue;
                fixedValueDistribution
                {
                    value 0.000001;
                }
        }
        }
    }

    dispersionModel none;

    patchInteractionModel localInteraction;

    localInteractionCoeffs
    {
        patches
        (
            WALL_IN
            {
                type rebound;
                e 0.97;
                mu 0.09;
            }
            WALL_OUT
            {
                type rebound;
                e 0.97;
                mu 0.09;
            }
            LUMEN_WALL
            {
                type stick;
            }
            OUTLET
            {
                type escape;
            }
            INLET
            {
                type rebound;
                e 0.97;
                mu 0.09;
            }
        );
    }

    heatTransferModel none;

    surfaceFilmModel none;

    packingModel implicit;

    explicitCoeffs
    {
        particleStressModel
        {
            type HarrisCrighton;
            alphaPacked 0.6;
            pSolid 10.0;
            beta 2.0;
            eps 1.0e-7;
        }
        correctionLimitingMethod
        {
            type absolute;
            e 0.9;
        }
    }

    implicitCoeffs
    {
        alphaMin 0.0001;
        rhoMin 1.0;
        applyLimiting true;
        applyGravity false;
        particleStressModel
        {
            type HarrisCrighton;
            alphaPacked 0.6;
            pSolid 5.0;
            beta 2.0;
            eps 1.0e-2;
        }
    }

    dampingModel none;

    isotropyModel stochastic;

    stochasticCoeffs
    {
        timeScaleModel
        {
            type isotropic;
            alphaPacked 0.6;
            e 0.9;
        }
    }

    stochasticCollisionModel none;

    radiation off;
}


cloudFunctions
{}
[/QUOTE]
TagsNo tags attached.

Activities

rmn_990

2019-04-29 23:34

reporter   ~0010438

If the context is not understandable, you can find them in these two files:

kinematicCloudProperties (3,682 bytes)
/*--------------------------------*- 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 off;

    maxCo           1.0;

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

    averagingMethod dual;

    integrationSchemes
    {
        U               Euler;
    }

    sourceTerms
    {
        schemes
        {
            U           semiImplicit 0.00078;
        }
    }
}

constantProperties
{
    rho0            1050;
    alphaMax        0.9;
}

subModels
{
    particleForces
    {
        ErgunWenYuDrag
        {
            alphac alphac;
        }
        gravity;
    }

    injectionModels
    {
        model1
        {
            type            patchInjection;
            parcelBasisType fixed;
            SOI             0;
            patch       INLET;
            duration        5;
	    nParticle       1;
            parcelsPerSecond 110000;
	    massTotal       2.84e-11;
            U0              (0 0 0);
            flowRateProfile constant 2.5e-8;
            sizeDistribution
            {
		type        fixedValue;
                fixedValueDistribution
                {
                    value   0.000001;
                }
	    }
        }
    }

    dispersionModel none;

    patchInteractionModel localInteraction;

    localInteractionCoeffs
    {
        patches
        (
            WALL_IN
            {
                type rebound;
                e    0.97;
                mu   0.09;
            }
            WALL_OUT
            {
                type rebound;
                e    0.97;
                mu   0.09;
            }
            LUMEN_WALL
            {
                type stick;
            }
            OUTLET
            {
                type escape;
            }
            INLET
            {
                type rebound;
                e    0.97;
                mu   0.09;
            }
        );
    }

    heatTransferModel none;

    surfaceFilmModel none;

    packingModel implicit;

    explicitCoeffs
    {
        particleStressModel
        {
            type HarrisCrighton;
            alphaPacked 0.6;
            pSolid 10.0;
            beta 2.0;
            eps 1.0e-7;
        }
        correctionLimitingMethod
        {
            type absolute;
            e 0.9;
        }
    }

    implicitCoeffs
    {
        alphaMin 0.0001;
        rhoMin 1.0;
        applyLimiting true;
        applyGravity false;
        particleStressModel
        {
            type HarrisCrighton;
            alphaPacked 0.6;
            pSolid 5.0;
            beta 2.0;
            eps 1.0e-2;
        }
    }

    dampingModel none;

    isotropyModel stochastic;

    stochasticCoeffs
    {
        timeScaleModel
        {
            type isotropic;
            alphaPacked 0.6;
            e 0.9;
        }
    }

    stochasticCollisionModel none;

    radiation off;
}


cloudFunctions
{}


// ************************************************************************* //
kinematicCloudProperties (3,682 bytes)
log_sample.txt (7,731 bytes)
Courant Number mean: 0.111537 max: 0.968023
Interface Courant Number mean: 0 max: 0
deltaT = 0.00666667
Time = 314.187

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
GAMG:  Solving for kinematicCloud:alpha, Initial residual = 0.00141744, Final residual = 0.00109698, No Iterations 1000
Cloud: kinematicCloud
    Current number of parcels       = 287172
    Current mass in system          = 1.57881e-10
    Linear momentum                 = (2.96341e-15 -1.85579e-15 1.12751e-14)
   |Linear momentum|                = 1.18048e-14
    Linear kinetic energy           = 6.06622e-16
    Injector model1:
      - parcels added               = 549986
      - mass introduced             = 3.02371e-10
    Parcel fate: system (number, mass)
      - escape                      = 262814, 1.44489e-10
    Parcel fate: patch WALL_IN (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch WALL_OUT (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch LUMEN_WALL (number, mass)
      - escape                      = 0, 0
      - stick                       = 944625979, 5.19349e-07
    Parcel fate: patch OUTLET (number, mass)
      - escape                      = 262814, 1.4448e-10
      - stick                       = 0, 0
    Parcel fate: patch INLET (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Min cell volume fraction        = 0
    Max cell volume fraction        = 2.14854e-05
    Min dense number of parcels     = 5.6004

Continous phase-1 volume fraction = 1  Min(alphac) = 0.999979  Max(alphac) = 1
PIMPLE: iteration 1
smoothSolver:  Solving for alpha.water, Initial residual = 0.0134216, Final residual = 0.0114306, No Iterations 100
Phase-1 volume fraction = 1  Min(alpha.water) = 1  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 1  Min(alpha.water) = 1  Max(alpha.water) = 1
smoothSolver:  Solving for alpha.water, Initial residual = 0.0116186, Final residual = 0.0115812, No Iterations 100
Phase-1 volume fraction = 1  Min(alpha.water) = 1  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 1  Min(alpha.water) = 1  Max(alpha.water) = 1
DICPCG:  Solving for p_rgh, Initial residual = 2.08234e-06, Final residual = 9.29223e-08, No Iterations 4
time step continuity errors : sum local = 5.07285e-11, global = 7.93747e-12, cumulative = -6.26446e-09
DICPCG:  Solving for p_rgh, Initial residual = 3.23116e-06, Final residual = 1.61229e-07, No Iterations 6
time step continuity errors : sum local = 8.80189e-11, global = 1.2696e-11, cumulative = -6.25177e-09
DICPCG:  Solving for p_rgh, Initial residual = 1.38829e-06, Final residual = 9.62392e-09, No Iterations 443
time step continuity errors : sum local = 5.25393e-12, global = -2.9407e-15, cumulative = -6.25177e-09
ExecutionTime = 40255.8 s  ClockTime = 40508 s

Courant Number mean: 0.111537 max: 0.968023
Interface Courant Number mean: 0 max: 0
deltaT = 0.00666667
Time = 314.193

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
GAMG:  Solving for kinematicCloud:alpha, Initial residual = 0.0014182, Final residual = 0.00109688, No Iterations 1000
Cloud: kinematicCloud
    Current number of parcels       = 287167
    Current mass in system          = 1.57878e-10
    Linear momentum                 = (2.97358e-15 -1.83809e-15 1.12724e-14)
   |Linear momentum|                = 1.1802e-14
    Linear kinetic energy           = 6.06593e-16
    Injector model1:
      - parcels added               = 549986
      - mass introduced             = 3.02371e-10
    Parcel fate: system (number, mass)
      - escape                      = 262819, 1.44492e-10
    Parcel fate: patch WALL_IN (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch WALL_OUT (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch LUMEN_WALL (number, mass)
      - escape                      = 0, 0
      - stick                       = 944675672, 5.19376e-07
    Parcel fate: patch OUTLET (number, mass)
      - escape                      = 262819, 1.44483e-10
      - stick                       = 0, 0
    Parcel fate: patch INLET (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Min cell volume fraction        = 0
    Max cell volume fraction        = 2.13323e-05
    Min dense number of parcels     = 5.5605

Continous phase-1 volume fraction = 1  Min(alphac) = 0.999979  Max(alphac) = 1
PIMPLE: iteration 1
smoothSolver:  Solving for alpha.water, Initial residual = 0.013389, Final residual = 0.0114096, No Iterations 100
Phase-1 volume fraction = 1  Min(alpha.water) = 1  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 1  Min(alpha.water) = 1  Max(alpha.water) = 1
smoothSolver:  Solving for alpha.water, Initial residual = 0.0115916, Final residual = 0.0115598, No Iterations 100
Phase-1 volume fraction = 1  Min(alpha.water) = 1  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 1  Min(alpha.water) = 1  Max(alpha.water) = 1
DICPCG:  Solving for p_rgh, Initial residual = 2.09388e-06, Final residual = 9.41566e-08, No Iterations 4
time step continuity errors : sum local = 5.14021e-11, global = -8.3902e-12, cumulative = -6.26016e-09
DICPCG:  Solving for p_rgh, Initial residual = 3.24005e-06, Final residual = 1.60079e-07, No Iterations 7
time step continuity errors : sum local = 8.73903e-11, global = -1.33005e-11, cumulative = -6.27346e-09
DICPCG:  Solving for p_rgh, Initial residual = 1.38674e-06, Final residual = 9.51116e-09, No Iterations 442
time step continuity errors : sum local = 5.19234e-12, global = 3.77542e-15, cumulative = -6.27346e-09
ExecutionTime = 40261.9 s  ClockTime = 40514 s

Courant Number mean: 0.111537 max: 0.968022
Interface Courant Number mean: 0 max: 0
deltaT = 0.00666667
Time = 314.2

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
GAMG:  Solving for kinematicCloud:alpha, Initial residual = 0.00141875, Final residual = 0.00109776, No Iterations 1000
Cloud: kinematicCloud
    Current number of parcels       = 287162
    Current mass in system          = 1.57876e-10
    Linear momentum                 = (2.96312e-15 -1.82337e-15 1.128e-14)
   |Linear momentum|                = 1.18044e-14
    Linear kinetic energy           = 6.06554e-16
    Injector model1:
      - parcels added               = 549986
      - mass introduced             = 3.02371e-10
    Parcel fate: system (number, mass)
      - escape                      = 262824, 1.44495e-10
    Parcel fate: patch WALL_IN (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch WALL_OUT (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch LUMEN_WALL (number, mass)
      - escape                      = 0, 0
      - stick                       = 944725374, 5.19403e-07
    Parcel fate: patch OUTLET (number, mass)
      - escape                      = 262824, 1.44485e-10
      - stick                       = 0, 0
    Parcel fate: patch INLET (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Min cell volume fraction        = 0
    Max cell volume fraction        = 2.34655e-05
    Min dense number of parcels     = 6.11655


log_sample.txt (7,731 bytes)

rmn_990

2019-04-30 19:00

reporter   ~0010445

I am using OpenFoam 18.12 and the solver which I used for this simulation is MPPICInterFoam

henry

2019-04-30 19:28

manager   ~0010446

OpenFoam 18.12 is not a valid version it is a fork of OpenFOAM which is not released or supported by the OpenFOAM Foundation. Please upgrade to OpenFOAM-6 or better OpenFOAM-dev and repeat the tests and report with details of how to reproduce the problem.

Issue History

Date Modified Username Field Change
2019-04-29 23:24 rmn_990 New Issue
2019-04-29 23:34 rmn_990 File Added: kinematicCloudProperties
2019-04-29 23:34 rmn_990 File Added: log_sample.txt
2019-04-29 23:34 rmn_990 Note Added: 0010438
2019-04-30 19:00 rmn_990 Note Added: 0010445
2019-04-30 19:28 henry Assigned To => henry
2019-04-30 19:28 henry Status new => closed
2019-04-30 19:28 henry Resolution open => unable to reproduce
2019-04-30 19:28 henry Note Added: 0010446