View Issue Details

IDProjectCategoryView StatusLast Update
0001058OpenFOAMBugpublic2016-05-04 20:40
Reporteruser649Assigned Tohenry  
PrioritynormalSeveritytweakReproducibilityalways
Status closedResolutionunable to reproduce 
PlatformLinuxOSUbuntuOS Version10.04
Summary0001058: LTSReactingParcelFoam with coneNozzleInjection producing abnormal spray pattern
DescriptionHi!

I would like to report a repeating problem I have with a simple water spray in a cylindrical pipe. I specified a conical spray pattern in reactingCloud1Properties and checked many times to make sure I have the correct parameters. The result I obtained with LTSReactingParcelFoam is not what I expected. The spray does not have the cone shape but instead it appears to be injecting on a 2D plane whereas the domain is 3D.

Attached is the screenshot of cross section of cylinder. As you can see, instead of a cone shape, its point of injection seems to be 2D (on Y-Z plane). I took this screenshot with paraview.

Your help on this is much appreciated.
Steps To ReproduceCreate a cylinder (or any shape) as a domain. Run RAS turbulence model to obtain a solution first, and then run LTSReactingParcelFoam for 200 iterations. The residuals were in the 1e-8 to 1e-9.
TagsNo tags attached.

Activities

user649

2013-10-25 13:37

 

reactingCloud1Properties (4,990 bytes)   
The reactingCloud1Properties code is as follow:

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

solution
{
    active          yes;

    transient       no; // yes;
    calcFrequency   1;
    maxTrackTime    5.0;
    maxCo           0.05;

    coupled         true;
    cellValueSourceCorrection off;

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


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

    integrationSchemes
    {
        U               Euler;
        T               analytical;
    }
}


constantProperties
{
    rho0            1000;
    T0              300;
    Cp0             4100;

    youngsModulus   1e9;
    poissonsRatio   0.35;

    epsilon0        1;
    f0              0.5;

    Tvap            273;
    Tbp             373;
    Pr              0.7;
    LDevol          0;
    hRetentionCoeff 1;

    constantVolume  false;
}


subModels
{
    particleForces
    {
        sphereDrag;
        gravity;
    }

    injectionModels
    {
        model1
        {
            type            coneNozzleInjection;
	    SOI			0;
            massFlowRate      0.00277; // 
  //          positionAxis       ((( 0 0 0.1 ) ( 0 0 1 ))); 
  	    position (0.0 0.0 0.0); // ok
  	    direction (0 0 1); // ok
            parcelBasisType mass;
  	    outerDiameter 0.00625; // ok
  	    innerDiameter 0.0062; // ok
            parcelsPerSecond 100;//100; //DEFAULT: 100
//            parcelsPerInjector 40000;//40000; //DEFAULT: 100
            duration        1; // NOTE: set to 1 for steady state
            injectionMethod disc;
            flowType        constantVelocity;//flowRateAndDischarge;//
            UMag            40.0;
            flowRateProfile constant 1;
            Cd              constant 0.9;

            thetaInner      constant 18.0; //from sprayCloudProperties
            thetaOuter      constant 21.0; //from sprayCloudProperties

            sizeDistribution
            {
                type RosinRammler;
                RosinRammlerDistribution
                {
                    minValue        20e-06; //20 micron
                    maxValue        60e-06; //60 micron
                    d		    40e-06; //40 micron
                    n	            4; //spray parameter
                }
            }
        }
    }

    dispersionModel none;//stochasticDispersionRAS;

    patchInteractionModel standardWallInteraction;

    heatTransferModel RanzMarshall;

    compositionModel singleMixtureFraction;

    phaseChangeModel liquidEvaporation;

    devolatilisationModel none;

    surfaceReactionModel none;

    surfaceFilmModel none; //from sprayCloudProperties

    atomizationModel none; //from sprayCloudProperties

    breakupModel    none; // ReitzKHRT; ReitzDiwakar; from sprayCloudProperties

    stochasticCollisionModel none; //from sprayCloudProperties

    radiation       off;

    standardWallInteractionCoeffs
    {
        type            rebound;
    }

    RanzMarshallCoeffs
    {
        BirdCorrection  off; //off
    }

    singleMixtureFractionCoeffs
    {
        phases
        (
            gas
            {
            }
            liquid
            {
                H2O 1;
            }
            solid
            {
            }
        );
        YGasTot0        0;
        YLiquidTot0     1;
        YSolidTot0      0;
    }

    liquidEvaporationCoeffs
    {
        enthalpyTransfer enthalpyDifference; // latentHeat;

        activeLiquids   ( H2O );
    }
}


cloudFunctions
{
 /*   patchPostProcessing
    {
        maxStoredParcels 100;
        patches         ( outlet );
    }
*/
    particleTracks
    {
        trackInterval   1;
        maxSamples      1000000;
        resetOnWrite    yes;
    }

}


// ************************************************************************* //
reactingCloud1Properties (4,990 bytes)   

henry

2015-01-01 14:27

manager   ~0003446

Could you please supply a test case which reproduces the problem in OpenFOAM-2.3.x?

henry

2016-05-04 20:40

manager   ~0006228

Orphaned report

Issue History

Date Modified Username Field Change
2013-10-25 13:37 user649 New Issue
2013-10-25 13:37 user649 File Added: reactingCloud1Properties
2015-01-01 14:27 henry Note Added: 0003446
2016-05-04 20:40 henry Note Added: 0006228
2016-05-04 20:40 henry Status new => closed
2016-05-04 20:40 henry Assigned To => henry
2016-05-04 20:40 henry Resolution open => unable to reproduce