[mturcios777]

In continuing my exploration of the new lagrangian spray libraries, I've tried to run sprayEngineFoam with a profile I've used before to simulate diesel injection. However, there is no injection at any time during the run.

I think it may have to do with the mapping from crank angle to seconds. In the old dieselSpray libraries there was mapping done with userTimeToTime so that the profile can be specified in crank angle degrees (which is pretty standard), which I haven't been able to find yet in the new libraries. Right now I've specified SOI, duration and profile in CAdegrees and get no injection. Is there another combination I should try? I would think specifying everything from the dictionary in crank angles makes the most sense.

I don't know if this qualifies as feature request or bug report, so I'll leave it here for consideration.

[mturcios777]

Totally forgot I hadn't reported; there was indeed a problem with the crank angle mapping. The version currently on git is working properly.

[NablaDyn]

Dear Foamers,

I'm a bit stuck trying to run sprayEngineFoam with fuel injection.

I set up my case as a mixture of the aachenBomb lagrangian tutorial and the one regarding engineFoam. Both based on OpenFOAM versions 4 and 5.

Annoyingly, no fuel is injected and albeit thourough investigation I did not find a solution, yet. I hope some of you can help me.

So, here's a snippet of the output of OpenFOAM during runtime:

Courant Number mean: 0.00010521857883 max: 0.197815497257
deltaT = 1.20141012565e-06
Crank angle = -179.594698012 CA-deg
deltaZ = 1.09762120504e-07
DICPCG:  Solving for yPsi, Initial residual = 9.03191839263e-09, Final residual = 9.03191839263e-09, No Iterations 0
clearance: 0.204997915028
Piston speed = 0.0913610749238 m/s

Solving 2-D cloud sprayCloud
Cloud: sprayCloud
    Current number of parcels       = 0
    Current mass in system          = 0
    Linear momentum                 = (0 0 0)
   |Linear momentum|                = 0
    Linear kinetic energy           = 0
    model1:
        number of parcels added     = 0
        mass introduced             = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 0, 0
    Mass transfer phase change      = 0
    D10, D32, Dmax (mu)             = 0, 0, 0
    Liquid penetration 95% mass (m) = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.00473860342429, Final residual

The engine motion including flow computation itself is doing as expected. But no parcels are injected. Here's the sprayCloudProperties 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      binary;
    class       dictionary;
    location    "constant";
    object      SprayCloudProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solution
{
    active                      true;
    coupled                     true;
    transient                   true;
    cellValueSourceCorrection   true;
    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              323;
    rho0            1000;
    Cp0             4187;
    constantVolume  false;
}

subModels
{
    particleForces
    {
        sphereDrag;
    }

    injectionModels
    {
        model1
        {
            type                coneNozzleInjection;
            SOI                 -180;
            massTotal           1e-05;
            parcelBasisType     mass;
            injectionMethod     disc;
            flowType            flowRateAndDischarge;
            outerDiameter       2.0e-03;
            innerDiameter       0;
            duration            10;
            position            (0 0 0.2);
            direction           (0 0 -1);
            parcelsPerSecond    1000000;
            flowRateProfile     table
            (
                (-181 0)
                (-180 0.3)
                (30 0.3)
                (31 0)
            );

            Cd          constant 0.9;
            thetaInner  constant 0.0;
            thetaOuter  constant 10.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;
    stochasticCollisionModel    none;
    radiation                   false;
    standardWallInteractionCoeffs
    {
        type    rebound;
    }
    RanzMarshallCoeffs
    {
        BirdCorrection  true;
    }
    singlePhaseMixtureCoeffs
    {
        phases
        (
            liquid
            {
                C7H16   1;
            }
        );
    }

    liquidEvaporationBoilCoeffs
    {
        enthalpyTransfer    enthalpyDifference;
        activeLiquids       ( C7H16 );
    }

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

    TABCoeffs
    {
        y0              0;
        yDot0           0;
        Cmu             10;
        Comega          8;
        WeCrit          12;
    }
}

cloudFunctions
{}

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

I already played around with the duration and start values as well as some other stuff but nothing seems to have an effect.

Any help is appreciated!

Kind regards,

Martin