[marcosgutierrez]

Hi everyone,

I am trying to simulate two opposed injector nozzle spray in a rectangular chamber.

The solver sprayFoam runs but, only one spray is displayed in Paraview.

Here is the code and what I added is commented.

Can you tell me please how to insert the code to simulate 2 injectors in openFOAM 211 CAELinux 2013 with sprayFoam and injectionModel: coneNozzleInjection?

Thank you for any advice.

Marcos

Below the code:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.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 yes;
cellValueSourceCorrection on;

sourceTerms
{
schemes
{
rho explicit 1;
U explicit 1;
Yi explicit 1;
hs explicit 1;
}
}

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

integrationSchemes
{
U Euler;
T analytical;
}
}


constantProperties
{
parcelTypeId 1;

rhoMin 23.23; //was 1e-15
TMin 310; //was 200
pMin 6100000; //was 1000
minParticleMass 1e-15;

rho0 834.197; //was 1000
T0 310; //was 320
Cp0 4187;

youngsModulus 1e9;
poissonsRatio 0.35;

epsilon0 1;
f0 0.5;
Pr 0.7;
Tvap 310; //was273
Tbp 335.3;

constantVolume false;
}


subModels
{
particleForces
{
sphereDrag;
}

injectionModel coneNozzleInjection;

dispersionModel none;

patchInteractionModel standardWallInteraction;

heatTransferModel RanzMarshall;

compositionModel singlePhaseMixture;

phaseChangeModel liquidEvaporationBoil;

surfaceFilmModel none;

atomizationModel none;

breakupModel ReitzDiwakar; // ReitzKHRT;

stochasticCollisionModel none;

radiation off;

coneNozzleInjectionCoeffs
{
SOI 0;
massTotal 7.7e-5; //was 6.0e-6
parcelBasisType mass;
injectionMethod disc;
flowType flowRateAndDischarge;
outerDiameter 4.0e-4; //was 1.9e-4;
innerDiameter 0;
duration 1.6e-3; //was 1.25e-3;
position ( 0 0.4995 0 );
direction ( 0 -1 0 );
parcelsPerSecond 20000000;
flowRateProfile table
(
(0.0000 0.044118)
(0.0001 0.058824)
(0.0002 0.044118)
(0.0003 0.049020)
(0.0004 0.053992)
(0.0005 0.058824)
(0.0006 0.058824)
(0.0007 0.053992)
(0.0008 0.044118)
(0.0009 0.049020)
(0.0010 0.049020)
(0.0011 0.053992)
(0.0012 0.049020)
(0.0013 0.049020)
(0.0014 0.058824)
(0.0015 0.049020)
(0.0016 0.049020)

); //the whole flow rate was modified

Cd constant 0.8; //was 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;
}
}
}

standardWallInteractionCoeffs
{
type rebound;
}

RanzMarshallCoeffs
{
BirdCorrection true;
}

singlePhaseMixtureCoeffs
{
phases
(
liquid
{
C7H16 1; //was C7H16
}
);
}

liquidEvaporationBoilCoeffs
{
enthalpyTransfer enthalpyDifference;

activeLiquids ( C7H16 ); //was C7H16
}

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

/************************************************** ***********
a new submodel is added with another position of the nozzle
************************************************** ***********/

{
particleForces
{
sphereDrag;
}

injectionModel coneNozzleInjection;

dispersionModel none;

patchInteractionModel standardWallInteraction;

heatTransferModel RanzMarshall;

compositionModel singlePhaseMixture;

phaseChangeModel liquidEvaporationBoil;

surfaceFilmModel none;

atomizationModel none;

breakupModel ReitzDiwakar; // ReitzKHRT;

stochasticCollisionModel none;

radiation off;

coneNozzleInjectionCoeffs
{
SOI 0;
massTotal 7.7e-5; //was 6.0e-6
parcelBasisType mass;
injectionMethod disc;
flowType flowRateAndDischarge;
outerDiameter 4.0e-4; //was 1.9e-4;
innerDiameter 0;
duration 1.6e-3; //was 1.25e-3;
position ( 0 0.0005 0 );
direction ( 0 1 0 );
parcelsPerSecond 20000000;
flowRateProfile table
(
(0.0000 0.044118)
(0.0001 0.058824)
(0.0002 0.044118)
(0.0003 0.049020)
(0.0004 0.053992)
(0.0005 0.058824)
(0.0006 0.058824)
(0.0007 0.053992)
(0.0008 0.044118)
(0.0009 0.049020)
(0.0010 0.049020)
(0.0011 0.053992)
(0.0012 0.049020)
(0.0013 0.049020)
(0.0014 0.058824)
(0.0015 0.049020)
(0.0016 0.049020)

); //the whole flow rate was modified

Cd constant 0.8; //was 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;
}
}
}

standardWallInteractionCoeffs
{
type rebound;
}

RanzMarshallCoeffs
{
BirdCorrection true;
}

singlePhaseMixtureCoeffs
{
phases
(
liquid
{
C7H16 1; //was C7H16
}
);
}

liquidEvaporationBoilCoeffs
{
enthalpyTransfer enthalpyDifference;

activeLiquids ( C7H16 ); //was C7H16
}

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
{}


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

[marcosgutierrez]

SOLVED!

Hi every one,

I found the solution:

coneNozzleInjection only works fro 1 Injector or 1 nozzle hole.

You have to change the injectionModel to coneInjection

below the corrected sprayCloudProperties

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.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 yes;
cellValueSourceCorrection on;

sourceTerms
{
schemes
{
rho explicit 1;
U explicit 1;
Yi explicit 1;
hs explicit 1;
}
}

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

integrationSchemes
{
U Euler;
T analytical;
}
}


constantProperties
{
parcelTypeId 1;

rhoMin 23.23; //was 1e-15
TMin 310; //was 200
pMin 6100000; //was 1000
minParticleMass 1e-15;

rho0 834.197; //was 1000
T0 310; //was 320
Cp0 4187;

youngsModulus 1e9;
poissonsRatio 0.35;

epsilon0 1;
f0 0.5;
Pr 0.7;
Tvap 310; //was273
Tbp 335.3;

constantVolume false;
}


subModels
{
particleForces
{
sphereDrag;
}

//injectionModel coneNozzleInjection;

injectionModel coneInjection;

dispersionModel none;

patchInteractionModel standardWallInteraction;

heatTransferModel RanzMarshall;

compositionModel singlePhaseMixture;

phaseChangeModel liquidEvaporationBoil;

surfaceFilmModel none;

atomizationModel none;

breakupModel ReitzDiwakar; // ReitzKHRT;

stochasticCollisionModel none;

radiation off;

//coneNozzleInjectionCoeffs
coneInjectionCoeffs
{
SOI 0;
massTotal 7.7e-5; //was 6.0e-6
parcelBasisType mass;
injectionMethod disc;
flowType flowRateAndDischarge;
outerDiameter 4.0e-4; //was 1.9e-4;
innerDiameter 0;
duration 1.6e-3; //was 1.25e-3;
//position ( 0 0.4995 0 );
//direction ( 0 -1 0 );

positionAxis //added
( //added
(( 0 0.4995 0 ) ( 0 -1 0 )) //added
(( 0 0.0005 0 ) ( 0 1 0 )) //added
); //added

parcelsPerInjector 100000; //added
parcelsPerSecond 100000; //modified was 10000000
flowRateProfile table
(
(0.0000 0.044118)
(0.0001 0.058824)
(0.0002 0.044118)
(0.0003 0.049020)
(0.0004 0.053992)
(0.0005 0.058824)
(0.0006 0.058824)
(0.0007 0.053992)
(0.0008 0.044118)
(0.0009 0.049020)
(0.0010 0.049020)
(0.0011 0.053992)
(0.0012 0.049020)
(0.0013 0.049020)
(0.0014 0.058824)
(0.0015 0.049020)
(0.0016 0.049020)

); //the whole flow rate was modified

Cd constant 0.8; //was 0.9

Umag constant 22.8; //added

thetaInner constant 0.0;
thetaOuter constant 10.0;

sizeDistribution
{
type RosinRammler;

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

standardWallInteractionCoeffs
{
type rebound;
}

RanzMarshallCoeffs
{
BirdCorrection true;
}

singlePhaseMixtureCoeffs
{
phases
(
liquid
{
C7H16 1; //was C7H16
}
);
}

liquidEvaporationBoilCoeffs
{
enthalpyTransfer enthalpyDifference;

activeLiquids ( C7H16 ); //was C7H16
}

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
{}


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