[ayhan515]

Hi All,
I'm using dieselEngineFoam with FOAM extend 3.1 to simulate combustion.

Simulation range is -25 +25 crank angle and simulation finished without error.

i use a global reaction (N2 inert)
C7H16 + 11O2 => 7CO2 + 8H2O

SOI =-24 CA
Injected Fuel (C7H16) mass = 12.8189 mg

Mass fraction of O2 and N2 species slightly decreases after the combustion.

@ -25
Total cylinder mass: 0.387525 gram (only air)
mass fraction 0.234 O2 (air) | 0,09068085 gram
0.766 N2 (air) | 0,29684415 gram


@ +25
Total cylinder mass: 0.400344 gram (air+fuel) [mass conservation ok]


mass fraction (from species.-25.dat)
C7H16 2.10186e-14
O2 0.233853 | 0,0936216454 gram
N2 0.765967 | 0,3066502926 gram
CO2 0.000131135
H2O 6.3041e-05


The problem is mass of O2 and N2 species increases after the combustion.
i use below code to write mass fraction of species(species.-25.dat). Maybe something wrong in this code.

<< runTime.value() << token::TAB 

<< thermo.p()[0]  << token::TAB
<< Y[0].name() << token::TAB << Y[0][0]<< token::TAB
<< Y[1].name() << token::TAB << Y[1][1]<< token::TAB
<< Y[2].name() << token::TAB << Y[2][2]<< token::TAB
<< Y[3].name() << token::TAB << Y[3][3]<< token::TAB
<< Y[4].name() << token::TAB << Y[4][4]<< token::TAB
<< thermo.T()[0] << token::TAB

how can i fix this problem?

thanks in advance.

log

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     3.1                                |
|   \\  /    A nd           | Web:         http://www.extend-project.de       |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build    : 3.1-6c0c77dee21a
Exec     : dieselEngineFoamtemmuz2014hrr
Date     : Aug 14 2014
Time     : 15:29:17
Host     : ayhan-virtual-machine
PID      : 24778
CtrlDict : /opt/foam/foam-extend-3.1/etc/controlDict
Case     : /home/ayhan/foam/ayhan-3.1/antorduzsilindir-foam31-multihole-enj-24-swirl-profil-20-derece-grafik-icinhrr
nProcs   : 1
SigFpe   : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create engine time

--> FOAM Warning : 
    From function IOstream::compressionEnum(const word&)
    in file db/IOstreams/IOstreams/IOstream.C at line 74
    bad compression specifier 'off', using 'uncompressed'
Create mesh for time = -25

Selecting engineMesh layered
deckHeight: 0.072
piston position: 0.0639363

Reading thermophysicalProperties
Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics>
Selecting thermodynamics package hsPsiMixtureThermo<reactingMixture<gasThermoPhysics>>
Selecting chemistryReader chemkinReader
Selecting chemistrySolver ode
ODEChemistryModel: Number of species = 5 and reactions = 1
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
    Prt             1;
}

Creating field DpDt

Creating field kinetic energy K


Reading g
Reading combustion properties

Constructing Spray
--> FOAM Warning : 
    From function Cloud<ParticleType>::initCloud(const bool checkClass)
    in file /opt/foam/foam-extend-3.1/src/lagrangian/basic/lnInclude/CloudIO.C at line 124
    Cannot read particle positions file 
    "/home/ayhan/foam/ayhan-3.1/antorduzsilindir-foam31-multihole-enj-24-swirl-profil-20-derece-grafik-icinhrr/-25/lagrangian/defaultCloud"
    assuming the initial cloud contains 0 particles.
Selecting injectorType multiHoleInjector
injectionPressureProfile_.size() = 22, massFlowRateProfile_.size() = 22
end constructor. commonRail den alınan multihole
Selecting atomizationModel off
Selecting dragModel standardDragModel
Selecting evaporationModel standardEvaporationModel
Selecting heatTransferModel RanzMarshall
Selecting wallModel reflect
Selecting breakupModel ReitzKHRT
Selecting collisionModel off
Selecting dispersionModel off
Selecting injectorModel hollowConeInjector
Selecting pdfType RosinRammler
Average Velocity for injector 0: 0.000127159 m/s, injection pressure = 19.9718 bar
Constructing three dimensional spray injection.
Courant Number mean: 0.000346976 max: 0.00280174 velocity magnitude: 5.81392
==startSummary.H yazdiriyor...
Total cylinder mass: 0.387525 gram

Starting time loop

Courant Number mean: 0.000346976 max: 0.00280174 velocity magnitude: 5.81392
Crank angle = -24.975 CA-deg
deltaZ = 7.91565e-06
clearance: 0.00805578
Piston speed = 4.17946 m/s
Piston konumu (mm)=         63.9442
Piston yerdeğiştirme (mm)=     0.00791565
Piston hızı (m/s)=         4.17946
, T = 789.788, p = 1.95045e+06, C7H16 = 0
Solving chemistry
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UEqn.H dosyası tamam 
DILUPBiCG:  Solving for C7H16, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for O2, Initial residual = 0.0030041, Final residual = 3.23546e-16, No Iterations 3
DILUPBiCG:  Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
YEqn.H dosyası tamam 
DILUPBiCG:  Solving for hs, Initial residual = 0.0086605, Final residual = 2.43842e-07, No Iterations 1
EEqn.H dosyası tamam 
DICPCG:  Solving for p, Initial residual = 0.965778, Final residual = 0.0184173, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.204614, Final residual = 0.00605899, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.00679623, Final residual = 0.000255068, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.56393e-08, global = 9.25292e-09, cumulative = 9.25292e-09
pEqn.H dosyası tamam 
DICPCG:  Solving for p, Initial residual = 0.000641585, Final residual = 1.83956e-05, No Iterations 1
DICPCG:  Solving for p, Initial residual = 1.8474e-05, Final residual = 6.7099e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 6.64488e-07, Final residual = 6.64488e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.67414e-11, global = 3.58652e-11, cumulative = 9.28879e-09
pEqn.H dosyası tamam 
DICPCG:  Solving for p, Initial residual = 1.44558e-05, Final residual = 3.33062e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 3.32992e-07, Final residual = 3.32992e-07, No Iterations 0
DICPCG:  Solving for p, Initial residual = 3.32992e-07, Final residual = 3.32992e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.34456e-11, global = 2.47535e-12, cumulative = 9.29126e-09
pEqn.H dosyası tamam 
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00150506, Final residual = 3.05437e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.000300666, Final residual = 8.07714e-09, No Iterations 1
Mean pressure (pa):        1.95275e+06
Mean temperature (K):        815.163
Mean u'(m/s):            2.3296
Hacim (cm-küp)            46.6023
Toplam silindir kütlesi (gram):    0.387525

Number of parcels in system.... | 0
Injected liquid mass........... | 0 mg
Liquid Mass in system.......... | 0 mg
SMD, Dmax...................... | 0 mu, 0 mu
Added gas mass................. | 2.71051e-13 mg
Evaporation Continuity Error... | 2.71051e-13 mg

HRRate = 0 J/s-W

Heat release1 [J]: 0

sumHRRate 0 in -24.975
Heat release2 [J]: 0

ExecutionTime = 0.11 s  ClockTime = 0 s

******************

ourant Number mean: 0.000315112 max: 0.00232653 velocity magnitude: 4.9143
Crank angle = 25 CA-deg
deltaZ = -7.91565e-06
clearance: 0.0080637
Piston speed = -4.17946 m/s
Piston konumu (mm)=         63.9363
Piston yerdeğiştirme (mm)=     -0.00791565
Piston hızı (m/s)=         -4.17946
, T = 733.705, p = 2.34941e+06, C7H16 = 2.10186e-14
Solving chemistry
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
UEqn.H dosyası tamam 
DILUPBiCG:  Solving for C7H16, Initial residual = 0.000109686, Final residual = 4.94516e-16, No Iterations 2
DILUPBiCG:  Solving for O2, Initial residual = 9.89698e-05, Final residual = 1.33232e-15, No Iterations 2
DILUPBiCG:  Solving for CO2, Initial residual = 0.000122413, Final residual = 2.71565e-14, No Iterations 2
DILUPBiCG:  Solving for H2O, Initial residual = 0.000122413, Final residual = 2.71565e-14, No Iterations 2
YEqn.H dosyası tamam 
DILUPBiCG:  Solving for hs, Initial residual = 0.00027299, Final residual = 7.1057e-10, No Iterations 1
EEqn.H dosyası tamam 
DICPCG:  Solving for p, Initial residual = 0.49718, Final residual = 0.018074, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.0209321, Final residual = 0.0014516, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.00145968, Final residual = 9.29612e-05, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.94866e-07, global = -2.41915e-08, cumulative = -1.4224e-05
pEqn.H dosyası tamam 
DICPCG:  Solving for p, Initial residual = 9.84205e-05, Final residual = 6.10055e-06, No Iterations 1
DICPCG:  Solving for p, Initial residual = 6.10029e-06, Final residual = 6.54184e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 6.52666e-07, Final residual = 6.52666e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.36808e-09, global = -6.15307e-10, cumulative = -1.42246e-05
pEqn.H dosyası tamam 
DICPCG:  Solving for p, Initial residual = 6.76295e-07, Final residual = 6.76295e-07, No Iterations 0
DICPCG:  Solving for p, Initial residual = 6.76295e-07, Final residual = 6.76295e-07, No Iterations 0
DICPCG:  Solving for p, Initial residual = 6.76295e-07, Final residual = 6.76295e-07, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.41761e-09, global = -6.15307e-10, cumulative = -1.42252e-05
pEqn.H dosyası tamam 
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00487561, Final residual = 5.09457e-09, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.0016339, Final residual = 3.83589e-09, No Iterations 1
Mean pressure (pa):        2.34247e+06
Mean temperature (K):        1095.94
Mean u'(m/s):            1.05203
Hacim (cm-küp)            46.6481
Toplam silindir kütlesi (gram):    0.400344

Number of parcels in system.... | 0
Injected liquid mass........... | 12.8189 mg
Liquid Mass in system.......... | 0 mg
SMD, Dmax...................... | 0 mu, 0 mu
Added gas mass................. | 12.8189 mg
Evaporation Continuity Error... | 1.22663e-05 mg


ExecutionTime = 594.75 s  ClockTime = 595 s

End