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; } } // ************************************************************************* //