[mbookin]

Does somebody know whether thermoParcel contain MPPIC submodels? Is there a way to check it?

It is said in the code it contains Kinematic parces sub-models and MPPIC parcel is a wraper which add new functionality to Kinematic parcel (see below).

ThermoParcel.H
Description
    Thermodynamic parcel class with one/two-way coupling with the continuous
    phase. Includes Kinematic parcel sub-models, plus:
    - heat transfer
/////////////////////////////////////////////
KinematicParcel.H
Description
    Kinematic parcel class with rotational motion (as spherical
    particles only) and one/two-way coupling with the continuous
    phase.
/////////////////////////////////////////////
MPPICParcel.H
Description
    Wrapper around kinematic parcel types to add MPPIC modelling

I see that Ucorrect is not the same in InjectionChannel tutorial if I substitute kinematicParcell with ThermoParcel there. Thats why parcels interact in different way if I use Kinematic or Thermo clouds. I can add case details if see any interest to this post.

[mbookin]

As I see on the forum thermoParcel includes mppic functionality as now developers simply transfered mppic methods into kinematicParcel and do not use any wrapping.

So now the question is why does kinematicParcel and thermoParcel show different collision behaviour having similar input parameters. Please see images of two clouds Ucorrection, more evolved clouds in coloring by its diameter and logs comparison below.

They are evidently different while cloud reading comparison from logs show that input is similar for thermo and kinematic parcel. I also attach cloud properties files for information.

kineticCloudProperties

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2306                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      kinematicCloudProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

    maxCo           0.3;

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

    averagingMethod dual;//basic;

    integrationSchemes
    {
        U               Euler;
    }

    sourceTerms
    {
        schemes
        {
            U semiImplicit 1;
        }
    }
}

constantProperties
{
    rho0            1000;

    alphaMax        0.9;
}

subModels
{
    particleForces
    {
        sphereDrag;
        //WenYuDrag
        //{
        //    alphac alphac;
        //}
        gravity;
//        interface
//        {
//            C  -10;
//            alpha alpha.water;
//        }
    }

    injectionModels
    {
        /*
            % matlab/octave calculation of flow rate
            alpha=0.1
            U=20
            A=0.01^2
            Q=U*A
            d=650e-6
            v=(4/3)*pi*(d/2)^3
            n=1
            rate=alpha*Q/(v*n)
        */

        model1
        {
            type            patchInjection;
            parcelBasisType fixed;
            patch           lowerInlet;
            U0              (18.7939 6.8404 0);
            nParticle       1;
            parcelsPerSecond 1390885;
            sizeDistribution
            {
                type            normal;
                normalDistribution
                {
                    mu              650e-6;
                    sigma           25e-6;
                    minValue        800e-6;
                    maxValue        1600e-6;
                }
            }
            flowRateProfile constant 1;
            massTotal       1;
            SOI             0;
            duration        60;
        }

        model2
        {
            type            patchInjection;
            parcelBasisType fixed;
            patch           upperInlet;
            U0              (18.7939 -6.8404 0);
            nParticle       1;
            parcelsPerSecond 1390885;
            sizeDistribution
            {
                type            normal;
                normalDistribution
                {
                    mu              650e-6;
                    sigma           25e-6;
                    minValue        500e-6;
                    maxValue        800e-6;
                }
            }
            flowRateProfile constant 1;
            massTotal       0;
            SOI             0;
            duration        60;
        }
    }

    dispersionModel none;

    heatTransferModel RanzMarshall;

    RanzMarshallCoeffs
    {
             BirdCorrection  true;
    }

    patchInteractionModel localInteraction;

    localInteractionCoeffs
    {
        patches
        (
            "(frontAndBack|walls)"
            {
                type rebound;
                e    1;
                mu   0;
            }
            "(lowerInlet|upperInlet|outlet)"
            {
                type escape;
            }
        );
    }

    surfaceFilmModel none;

    packingModel none;

    dampingModel relaxation;

    relaxationCoeffs
    {
        timeScaleModel
        {
            type nonEquilibrium;
            alphaPacked 0.58;
            e 0.9;
        }
    }

    isotropyModel stochastic;

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

    stochasticCollisionModel none;
}

cloudFunctions
{}


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

thermoCloudProperties

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2306                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      kinematicCloudProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

    maxCo           0.3;

    interpolationSchemes
    {
        rho             cell;
        U               cellPoint;
//        mu              cell;
        gradAlpha       cellPoint;
        thermo:mu       cell;
        T               cell;
        Cp              cell;
        kappa           cell;
        p               cell;
    }

    averagingMethod dual;//basic;

    integrationSchemes
    {
        U               Euler;
        T               analytical;
    }

    sourceTerms
    {
        schemes
        {
//            U semiImplicit 1;
            rho             explicit 1;
            U               explicit 1;
//            Yi              explicit 1;
            h               explicit 1;
            radiation       explicit 1;
        }
    }
}

constantProperties
{
    rho0            1000;
    T0              293;
    // Place holders for rho0 and Cp0
    // - reset from liquid properties using T0
    Cp0             4187;
    constantVolume  false;
    alphaMax        0.9;
}

subModels
{
    particleForces
    {
        sphereDrag;

        //WenYuDrag
    //    {
    //        alphac alphac;
    //    }
        gravity;
//        interface
//        {
//            C  -10;
//            alpha alpha.water;
//        }
    }

    injectionModels
    {
        /*
            % matlab/octave calculation of flow rate
            alpha=0.1
            U=20
            A=0.01^2
            Q=U*A
            d=650e-6
            v=(4/3)*pi*(d/2)^3
            n=1
            rate=alpha*Q/(v*n)
        */

        model1
        {
            type            patchInjection;
            parcelBasisType fixed;
            patch           lowerInlet;
            U0              (18.7939 6.8404 0);
            nParticle       1;
            parcelsPerSecond 1390885;
            sizeDistribution
            {
                type            normal;
                normalDistribution
                {
                    mu              650e-6;
                    sigma           25e-6;
                    minValue        800e-6;
                    maxValue        1600e-6;
                }
            }
            flowRateProfile constant 1;
            massTotal       1;
            SOI             0;
            duration        60;
        }

        model2
        {
            type            patchInjection;
            parcelBasisType fixed;
            patch           upperInlet;
            U0              (18.7939 -6.8404 0);
            nParticle       1;
            parcelsPerSecond 1390885;
            sizeDistribution
            {
                type            normal;
                normalDistribution
                {
                    mu              650e-6;
                    sigma           25e-6;
                    minValue        500e-6;
                    maxValue        800e-6;
                }
            }
            flowRateProfile constant 1;
            massTotal       0;
            SOI             0;
            duration        60;
        }
    }

    dispersionModel none;

    patchInteractionModel localInteraction;

    localInteractionCoeffs
    {
        patches
        (
            "(frontAndBack|walls)"
            {
                type rebound;
                e    1;
                mu   0;
            }
            "(lowerInlet|upperInlet|outlet)"
            {
                type escape;
            }
        );
    }

    surfaceFilmModel none;

    heatTransferModel RanzMarshall;

    RanzMarshallCoeffs
    {
         BirdCorrection  true;
    }

    packingModel none;

    radiation       off;

    dampingModel relaxation;

    relaxationCoeffs
    {
        timeScaleModel
        {
            type nonEquilibrium;
            alphaPacked 0.58;
            e 0.9;
        }
    }

    isotropyModel stochastic;

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

    stochasticCollisionModel none;
}

cloudFunctions
{}


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