[sammy1]

Hi all,

I've been running some cases in sprayFoam to simulate the breakup and dispersion of water droplets in different air flows. My first case was a liquid jet injected into a transverse crossflow that worked quite nicely, however, my second case is a nozzle spray into a vertical plume of air and is currently experiencing issues during parcel injection. All of my cases are 2D.

I have attached a zipped case file below. Ordinarily I would inject the parcels once the gaseous phase has been allowed to develop, however, due to the constraints of file-size for uploading the case, I've just got the injection starting from time 0. Once injection begins, the solver runs for a few iterations and then stops and hangs on the lines:

Solving 2-D cloud sprayCloud

Cloud: sprayCloud injector: model1
    Added 120 new parcels

The solver sits in this state until I ctrl+c to exit.

I've gone back to run my jet in crossflow solver, and it is still working. I have also copied the ./constant/sprayCloudProperties from the crossflow solver into the plume case and the issue comes up, so I believe this error could be arising due to my boundary conditions in ./0? I've printed my sprayCloudProperties below.

If anyone could take some time to look at the case it would be greatly appreciated. Apologies if I've missed anything with regards to forum rules/post content - this is my first post ever! Please let me know if I can provide any information that might help someone give me a hand.
Cheers,
Sam

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

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

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

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

    integrationSchemes
    {
        U               Euler;
        T               analytical;
    }
}


constantProperties
{
    T0              298;

    // place holders for rho0 and Cp0
    // - reset from liquid properties using T0
    rho0            1000;
    Cp0             4187;

    constantVolume  false;
}


subModels
{
    particleForces
    {
        sphereDrag;
		gravity;
    }

    injectionModels
    {
        model1
        {
            type            coneNozzleInjection;
            SOI             0.25000;
            massTotal       1.4e-4;
            parcelBasisType mass;
            injectionMethod disc;
            flowType        constantVelocity;
			UMag			18;
            outerDiameter   1.5e-4;
            innerDiameter   0;
            duration        0.1;
            position        (0.1 0.275 0);
            direction       (0 -1 0);
            parcelsPerSecond 20000000;
            flowRateProfile table
				(
				(0				0)
				(0.25			1.4e-3)
				(0.35			0)
				(0.36			0)
				);

            Cd              constant 0.9;

            thetaInner      constant 0.0;
            thetaOuter      constant 30.0;

            sizeDistribution
            {
                type        RosinRammler;

                RosinRammlerDistribution
                {
                    minValue        1e-06;
                    maxValue        0.00015;
                    d               0.00015;
                    n               3;
                }
            }
        }
    }

    dispersionModel none;

    patchInteractionModel standardWallInteraction;

    heatTransferModel RanzMarshall;

    compositionModel singlePhaseMixture;

    phaseChangeModel liquidEvaporationBoil;

    surfaceFilmModel none;

    atomizationModel none;

    breakupModel    ReitzDiwakar; // ReitzKHRT;

    stochasticCollisionModel none;

    radiation       off;

    standardWallInteractionCoeffs
    {
        type            escape;
    }

    RanzMarshallCoeffs
    {
        BirdCorrection  true;
    }

    singlePhaseMixtureCoeffs
    {
        phases
        (
            liquid
            {
                H2O               1;
            }
        );
    }

    liquidEvaporationBoilCoeffs
    {
        enthalpyTransfer enthalpyDifference;

        activeLiquids    ( H2O );
    }

    ReitzDiwakarCoeffs
    {
        solveOscillationEq yes;
        Cbag            6;
        Cb              0.785;
        Cstrip          0.5;
        Cs              10;
    }

/*
    ReitzKHRTCoeffs
    {
        solveOscillationEq yes;
        B0              0.61;
        B1              40;
        Ctau            1;
        CRT             0.1;
        msLimit         0.2;
        WeberLimit      6;
    }
*/
    TABCoeffs
    {
        y0              0;
        yDot0           0;
        Cmu             10;
        Comega          8;
        WeCrit          12;
    }
}


cloudFunctions
{}


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