CFD Online Discussion Forums

Page 1 of 2

Show 40 post(s) from this thread on one page

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - swirlinjector (https://www.cfd-online.com/Forums/openfoam-solving/65032-swirlinjector.html)

az1362f

June 1, 2009 22:03

swirlinjector

Hello to all,

I used swirlinjector for simulating dieselfoam solver in aachenbomb case instead of unitinjector but I face this error:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : dieselFoam
Date : Jun 01 2009
Time : 17:57:59
Host : spray
PID : 25432
Case : /home/openfoam1.5/OpenFOAM/openfoam1.5-1.5/run/tutorials/dieselFoam/aachenBomb
nProcs : 1

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

Create mesh for time = 0

Reading thermophysicalProperties
Selecting thermodynamics package hMixtureThermo<reactingMixture>
Selecting chemistryReader chemkinReader
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
alphah 1;
alphak 1;
alphaEps 0.76923;
muLimiter on;
Lsgs 0.0002;
}

Creating field DpDt

Constructing chemical mechanism
Selecting ODE solver SIBS
chemistryModel::chemistryModel: Number of species = 6 and reactions = 1

Reading environmentalProperties
Reading combustion properties

Constructing Spray
Selecting injectorType swirlInjector
Selecting atomizationModel off
Selecting dragModel standardDragModel
Selecting evaporationModel off
Selecting heatTransferModel RanzMarshall
Selecting wallModel reflect
Selecting breakupModel ReitzKHRT
Selecting collisionModel off
Selecting dispersionModel off
*** glibc detected *** dieselFoam: free(): invalid next size (fast): 0x0853b3f8 ***
======= Backtrace: =========
/lib/libc.so.6[0xb69bc4b6]
/lib/libc.so.6(cfree+0x89)[0xb69be179]
/home/openfoam1.5/OpenFOAM/ThirdParty/gcc-4.3.1/platforms/linux/lib/libstdc++.so.6(_ZdlPv+0x21)[0xb6b85051]
/home/openfoam1.5/OpenFOAM/ThirdParty/gcc-4.3.1/platforms/linux/lib/libstdc++.so.6(_ZNSs4_Rep10_M_destroyERKSaIcE+0x1d )[0xb6b6361d]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZN4Foam10dictionary3addEPNS_5entry Eb+0x52c)[0xb6cff45c]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZN4Foam10dictionary3addERKNS_4word ERKS0_b+0x5c)[0xb6cffeac]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZN4Foam7IOerror4exitEi+0x153)[0xb6cb7a03]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZNK4Foam10dictionary11lookupEntryE RKNS_4wordEb+0xec)[0xb6cfd98c]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZNK4Foam10dictionary6lookupERKNS_4 wordEb+0x2c)[0xb6cfd9cc]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libliquidMixture.so(_ZN4Foam13liquidMixtureC1ERKNS _10dictionaryE+0x821)[0xb70b1061]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libliquidMixture.so(_ZN4Foam13liquidMixture3NewERK NS_10dictionaryE+0x3a)[0xb70b13fa]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libdieselSpray.so(_ZN4Foam5sprayC1ERKNS_21volPoint InterpolationERKNS_14GeometricFieldINS_6VectorIdEE NS_12fvPatchFieldENS_7volMeshEEERKNS4_IdS7_S8_EESE _SE_RKNS_17combustionMixtureERKNS_7PtrListINS_19su therlandTransportINS_12specieThermoINS_11janafTher moINS_10perfectGasEEEEEEEEERKNS_10dictionaryESV_+0 x3f5)[0xb714a925]
dieselFoam[0x805edd4]
/lib/libc.so.6(__libc_start_main+0xe0)[0xb696bfe0]
dieselFoam(_ZNK4Foam11regIOobject11writeObjectENS_ 8IOstream12streamFormatENS1_13versionNumberENS1_15 compressionTypeE+0xcd)[0x805d021]
======= Memory map: ========
08048000-080be000 r-xp 00000000 08:0a 6275793 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam
080be000-080bf000 r--p 00076000 08:0a 6275793 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam
080bf000-080c0000 rw-p 00077000 08:0a 6275793 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam
080c0000-08562000 rw-p 080c0000 00:00 0 [heap]
b6100000-b6121000 rw-p b6100000 00:00 0
b6121000-b6200000 ---p b6121000 00:00 0
b627a000-b63d3000 rw-p b627a000 00:00 0
b63d3000-b6408000 r--s 00000000 08:09 1231203 /var/run/nscd/passwd
b6408000-b640a000 rw-p b6408000 00:00 0
b640a000-b640c000 r-xp 00000000 08:09 1855055 /lib/libutil-2.6.1.so
b640c000-b640e000 rw-p 00001000 08:09 1855055 /lib/libutil-2.6.1.so
b640e000-b640f000 rw-p b640e000 00:00 0
b640f000-b6422000 r-xp 00000000 08:09 1855032 /lib/libnsl-2.6.1.so
b6422000-b6424000 rw-p 00012000 08:09 1855032 /lib/libnsl-2.6.1.so
b6424000-b6426000 rw-p b6424000 00:00 0
b6426000-b645b000 r-xp 00000000 08:0a 6358081 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-pal.so.0.0.0
b645b000-b645c000 r--p 00034000 08:0a 6358081 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-pal.so.0.0.0
b645c000-b645d000 rw-p 00035000 08:0a 6358081 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-pal.so.0.0.0
b645d000-b647f000 rw-p b645d000 00:00 0
b647f000-b64cf000 r-xp 00000000 08:0a 6358083 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-rte.so.0.0.0
b64cf000-b64d0000 r--p 00050000 08:0a 6358083 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-rte.so.0.0.0
b64d0000-b64d2000 rw-p 00051000 08:0a 6358083 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-rte.so.0.0.0
b64d2000-b6571000 r-xp 00000000 08:0a 6214169 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESfilters.so
b6571000-b6573000 r--p 0009e000 08:0a 6214169 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESfilters.so
b6573000-b6574000 rw-p 000a0000 08:0a 6214169 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESfilters.so
b6574000-b6575000 rw-p b6574000 00:00 0
b6575000-b65e4000 r-xp 00000000 08:0a 6358070 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libmpi.so.0.0.0
b65e4000-b65e5000 r--p 0006e000 08:0a 6358070 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libmpi.so.0.0.0
b65e5000-b65ec000 rw-p 0006f000 08:0a 6358070 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libmpi.so.0.0.0
b65ec000-b65f2000 rw-p b65ec000 00:00 0
b65f2000-b6606000 r-xp 00000000 08:0a 6214141 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libz.so
b6606000-b6607000 r--p 00013000 08:0a 6214141 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libz.so
b6607000-b6608000 rw-p 00014000 08:0a 6214141 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libz.so
b6608000-b662e000 r-xp 00000000 08:0a 6214167 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESdeltas.so
b662e000-b662f000 r--p 00026000 08:0a 6214167 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESdeltas.so
b662f000-b6630000 rw-p 00027000 08:0a 6214167 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESdeltas.so
b6630000-b6734000 r-xp 00000000 08:0a 6214185 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libcompressibleLESModels.so
b6734000-b6737000 r--p 00103000 08:0a 6214185 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.Aborted
openfoam1.5@spray:~/OpenFOAM/openfoam1.5-1.5/run/tutorials/dieselFoam/aachenBomb>

what is the problem?
is it related to injector model that i used?I used hollowconeinjector model as injector model.

please help me.

az1362f

June 2, 2009 12:00

swirl injector

niklas

June 3, 2009 02:46

did it work before you changed injector type?

what is the temperature of the gas and liquid ?

az1362f

June 3, 2009 11:26

swirl injector

Hello Dear Dr;

before every things I am thankful for your replies dear.
yes it worked even when I changed injector type to commonrail injector and it still worked,but when I changed injector type to swirlinjector and just copy and paste swirlInjectorProps instead of commonRailInjectorProps in injectorProperties file,it generated this error.
I am confused.
is it related to injector model?should I change injector model to pressure swirl?

thanks
With Best Regard

az1362f

June 3, 2009 11:29

pressure swirl injector model

It seems that I have to define three parameter for pressure swirl injector including:1-ConeAngle 2-ConeInterval 3-maxKv?do I have to define any more parameters?
what are these constant in the pressure swirl model:1-ConeInterval 2-maxKv?
Is there any reference or article that explain this model?

thanks
best regard

niklas

June 4, 2009 02:18

Quote:

Originally Posted by az1362f (Post 218147)

again...
check the source file if you are uncertain about things.
src/lagrangian/dieselSpray/spraySubModels/injectorModel/pressureSwirl/pressureSwirlInjector.C

coneAngle is the umbrellaAngle for the spray and coneInterval is the sprayAngle. I dont know if those terms are clearer or not.
but if you consult the source code the injected parcel will have the angle
calculated like this

Code:

    scalar c = rndGen_.scalar01();

    angle_ = coneAngle_[n]  + 2.0 * coneInterval_[n] * (0.5 - c) ;

maxKv is the injectors maximum allowed discharge coeffients,
which is used to calculate the velocity

Code:

scalar v = kV * sqrt(2.0*deltaPressure/rhoFuel);

az1362f

June 4, 2009 07:52

Quote:

Originally Posted by niklas (Post 218206)

again...
did it work before you changed injector type?

what is the temperature of the gas and liquid ?
[/CODE]

Hello Dear Dr;

before every things I am thankful for your replies dear.
yes it worked even when I changed injector type to commonrail injector and it still worked,but when I changed injector type to swirlinjector and just copy and paste swirlInjectorProps instead of commonRailInjectorProps in injectorProperties file,it generated this error.tempretures is the same as original commonrailprops.
I am confused.I don't know what is my mistake.
is it related to injector model?should I change injector model to pressure swirl?

thanks
With Best Regard

az1362f

June 6, 2009 08:37

swirl injector bug

1 Attachment(s)

Hi

I used swirlinjector for simulating dieselfoam solver in aachenbomb case instead of unitinjector but I face this error:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : dieselFoam
Date : Jun 01 2009
Time : 17:57:59
Host : spray
PID : 25432
Case : /home/openfoam1.5/OpenFOAM/openfoam1.5-1.5/run/tutorials/dieselFoam/aachenBomb
nProcs : 1

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

Create mesh for time = 0

Reading thermophysicalProperties
Selecting thermodynamics package hMixtureThermo<reactingMixture>
Selecting chemistryReader chemkinReader
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
alphah 1;
alphak 1;
alphaEps 0.76923;
muLimiter on;
Lsgs 0.0002;
}

Creating field DpDt

Constructing chemical mechanism
Selecting ODE solver SIBS
chemistryModel::chemistryModel: Number of species = 6 and reactions = 1

Reading environmentalProperties
Reading combustion properties

Constructing Spray
Selecting injectorType swirlInjector
Selecting atomizationModel off
Selecting dragModel standardDragModel
Selecting evaporationModel off
Selecting heatTransferModel RanzMarshall
Selecting wallModel reflect
Selecting breakupModel ReitzKHRT
Selecting collisionModel off
Selecting dispersionModel off
*** glibc detected *** dieselFoam: free(): invalid next size (fast): 0x0853b3f8 ***
======= Backtrace: =========
/lib/libc.so.6[0xb69bc4b6]
/lib/libc.so.6(cfree+0x89)[0xb69be179]
/home/openfoam1.5/OpenFOAM/ThirdParty/gcc-4.3.1/platforms/linux/lib/libstdc++.so.6(_ZdlPv+0x21)[0xb6b85051]
/home/openfoam1.5/OpenFOAM/ThirdParty/gcc-4.3.1/platforms/linux/lib/libstdc++.so.6(_ZNSs4_Rep10_M_destroyERKSaIcE+0x1d )[0xb6b6361d]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZN4Foam10dictionary3addEPNS_5entry Eb+0x52c)[0xb6cff45c]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZN4Foam10dictionary3addERKNS_4word ERKS0_b+0x5c)[0xb6cffeac]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZN4Foam7IOerror4exitEi+0x153)[0xb6cb7a03]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZNK4Foam10dictionary11lookupEntryE RKNS_4wordEb+0xec)[0xb6cfd98c]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZNK4Foam10dictionary6lookupERKNS_4 wordEb+0x2c)[0xb6cfd9cc]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libliquidMixture.so(_ZN4Foam13liquidMixtureC1ERKNS _10dictionaryE+0x821)[0xb70b1061]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libliquidMixture.so(_ZN4Foam13liquidMixture3NewERK NS_10dictionaryE+0x3a)[0xb70b13fa]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libdieselSpray.so(_ZN4Foam5sprayC1ERKNS_21volPoint InterpolationERKNS_14GeometricFieldINS_6VectorIdEE NS_12fvPatchFieldENS_7volMeshEEERKNS4_IdS7_S8_EESE _SE_RKNS_17combustionMixtureERKNS_7PtrListINS_19su therlandTransportINS_12specieThermoINS_11janafTher moINS_10perfectGasEEEEEEEEERKNS_10dictionaryESV_+0 x3f5)[0xb714a925]
dieselFoam[0x805edd4]
/lib/libc.so.6(__libc_start_main+0xe0)[0xb696bfe0]
dieselFoam(_ZNK4Foam11regIOobject11writeObjectENS_ 8IOstream12streamFormatENS1_13versionNumberENS1_15 compressionTypeE+0xcd)[0x805d021]
======= Memory map: ========
08048000-080be000 r-xp 00000000 08:0a 6275793 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam
080be000-080bf000 r--p 00076000 08:0a 6275793 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam
080bf000-080c0000 rw-p 00077000 08:0a 6275793 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam
080c0000-08562000 rw-p 080c0000 00:00 0 [heap]
b6100000-b6121000 rw-p b6100000 00:00 0
b6121000-b6200000 ---p b6121000 00:00 0
b627a000-b63d3000 rw-p b627a000 00:00 0
b63d3000-b6408000 r--s 00000000 08:09 1231203 /var/run/nscd/passwd
b6408000-b640a000 rw-p b6408000 00:00 0
b640a000-b640c000 r-xp 00000000 08:09 1855055 /lib/libutil-2.6.1.so
b640c000-b640e000 rw-p 00001000 08:09 1855055 /lib/libutil-2.6.1.so
b640e000-b640f000 rw-p b640e000 00:00 0
b640f000-b6422000 r-xp 00000000 08:09 1855032 /lib/libnsl-2.6.1.so
b6422000-b6424000 rw-p 00012000 08:09 1855032 /lib/libnsl-2.6.1.so
b6424000-b6426000 rw-p b6424000 00:00 0
b6426000-b645b000 r-xp 00000000 08:0a 6358081 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-pal.so.0.0.0
b645b000-b645c000 r--p 00034000 08:0a 6358081 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-pal.so.0.0.0
b645c000-b645d000 rw-p 00035000 08:0a 6358081 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-pal.so.0.0.0
b645d000-b647f000 rw-p b645d000 00:00 0
b647f000-b64cf000 r-xp 00000000 08:0a 6358083 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-rte.so.0.0.0
b64cf000-b64d0000 r--p 00050000 08:0a 6358083 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-rte.so.0.0.0
b64d0000-b64d2000 rw-p 00051000 08:0a 6358083 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-rte.so.0.0.0
b64d2000-b6571000 r-xp 00000000 08:0a 6214169 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESfilters.so
b6571000-b6573000 r--p 0009e000 08:0a 6214169 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESfilters.so
b6573000-b6574000 rw-p 000a0000 08:0a 6214169 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESfilters.so
b6574000-b6575000 rw-p b6574000 00:00 0
b6575000-b65e4000 r-xp 00000000 08:0a 6358070 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libmpi.so.0.0.0
b65e4000-b65e5000 r--p 0006e000 08:0a 6358070 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libmpi.so.0.0.0
b65e5000-b65ec000 rw-p 0006f000 08:0a 6358070 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libmpi.so.0.0.0
b65ec000-b65f2000 rw-p b65ec000 00:00 0
b65f2000-b6606000 r-xp 00000000 08:0a 6214141 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libz.so
b6606000-b6607000 r--p 00013000 08:0a 6214141 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libz.so
b6607000-b6608000 rw-p 00014000 08:0a 6214141 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libz.so
b6608000-b662e000 r-xp 00000000 08:0a 6214167 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESdeltas.so
b662e000-b662f000 r--p 00026000 08:0a 6214167 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESdeltas.so
b662f000-b6630000 rw-p 00027000 08:0a 6214167 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESdeltas.so
b6630000-b6734000 r-xp 00000000 08:0a 6214185 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libcompressibleLESModels.so
b6734000-b6737000 r--p 00103000 08:0a 6214185 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.Aborted
openfoam1.5@spray:~/OpenFOAM/openfoam1.5-1.5/run/tutorials/dieselFoam/aachenBomb>

error is :

*** glibc detected *** dieselFoam: free(): invalid next size (fast): 0x0853b3f8 ***

or in some case:

*** glibc detected *** dieselFoam: malloc(): memory corruption (fast): 0x0b5a3a70 ***

I changed several injector model but none of them worked.
it seems that dieselfoam has a bug.
could any body help me.
it is my injectorproperties file that I used:

please help me it is urgent.
thanks

az1362f

June 6, 2009 08:41

swirl injector bug

1 Attachment(s)

Hi

I used swirlinjector for simulating dieselfoam solver in aachenbomb case instead of unitinjector but I face this error:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : dieselFoam
Date : Jun 01 2009
Time : 17:57:59
Host : spray
PID : 25432
Case : /home/openfoam1.5/OpenFOAM/openfoam1.5-1.5/run/tutorials/dieselFoam/aachenBomb
nProcs : 1

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

Create mesh for time = 0

Reading thermophysicalProperties
Selecting thermodynamics package hMixtureThermo<reactingMixture>
Selecting chemistryReader chemkinReader
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
alphah 1;
alphak 1;
alphaEps 0.76923;
muLimiter on;
Lsgs 0.0002;
}

Creating field DpDt

Constructing chemical mechanism
Selecting ODE solver SIBS
chemistryModel::chemistryModel: Number of species = 6 and reactions = 1

Reading environmentalProperties
Reading combustion properties

Constructing Spray
Selecting injectorType swirlInjector
Selecting atomizationModel off
Selecting dragModel standardDragModel
Selecting evaporationModel off
Selecting heatTransferModel RanzMarshall
Selecting wallModel reflect
Selecting breakupModel ReitzKHRT
Selecting collisionModel off
Selecting dispersionModel off
*** glibc detected *** dieselFoam: free(): invalid next size (fast): 0x0853b3f8 ***
======= Backtrace: =========
/lib/libc.so.6[0xb69bc4b6]
/lib/libc.so.6(cfree+0x89)[0xb69be179]
/home/openfoam1.5/OpenFOAM/ThirdParty/gcc-4.3.1/platforms/linux/lib/libstdc++.so.6(_ZdlPv+0x21)[0xb6b85051]
/home/openfoam1.5/OpenFOAM/ThirdParty/gcc-4.3.1/platforms/linux/lib/libstdc++.so.6(_ZNSs4_Rep10_M_destroyERKSaIcE+0x1d )[0xb6b6361d]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZN4Foam10dictionary3addEPNS_5entry Eb+0x52c)[0xb6cff45c]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZN4Foam10dictionary3addERKNS_4word ERKS0_b+0x5c)[0xb6cffeac]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZN4Foam7IOerror4exitEi+0x153)[0xb6cb7a03]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZNK4Foam10dictionary11lookupEntryE RKNS_4wordEb+0xec)[0xb6cfd98c]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so(_ZNK4Foam10dictionary6lookupERKNS_4 wordEb+0x2c)[0xb6cfd9cc]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libliquidMixture.so(_ZN4Foam13liquidMixtureC1ERKNS _10dictionaryE+0x821)[0xb70b1061]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libliquidMixture.so(_ZN4Foam13liquidMixture3NewERK NS_10dictionaryE+0x3a)[0xb70b13fa]
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libdieselSpray.so(_ZN4Foam5sprayC1ERKNS_21volPoint InterpolationERKNS_14GeometricFieldINS_6VectorIdEE NS_12fvPatchFieldENS_7volMeshEEERKNS4_IdS7_S8_EESE _SE_RKNS_17combustionMixtureERKNS_7PtrListINS_19su therlandTransportINS_12specieThermoINS_11janafTher moINS_10perfectGasEEEEEEEEERKNS_10dictionaryESV_+0 x3f5)[0xb714a925]
dieselFoam[0x805edd4]
/lib/libc.so.6(__libc_start_main+0xe0)[0xb696bfe0]
dieselFoam(_ZNK4Foam11regIOobject11writeObjectENS_ 8IOstream12streamFormatENS1_13versionNumberENS1_15 compressionTypeE+0xcd)[0x805d021]
======= Memory map: ========
08048000-080be000 r-xp 00000000 08:0a 6275793 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam
080be000-080bf000 r--p 00076000 08:0a 6275793 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam
080bf000-080c0000 rw-p 00077000 08:0a 6275793 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/dieselFoam
080c0000-08562000 rw-p 080c0000 00:00 0 [heap]
b6100000-b6121000 rw-p b6100000 00:00 0
b6121000-b6200000 ---p b6121000 00:00 0
b627a000-b63d3000 rw-p b627a000 00:00 0
b63d3000-b6408000 r--s 00000000 08:09 1231203 /var/run/nscd/passwd
b6408000-b640a000 rw-p b6408000 00:00 0
b640a000-b640c000 r-xp 00000000 08:09 1855055 /lib/libutil-2.6.1.so
b640c000-b640e000 rw-p 00001000 08:09 1855055 /lib/libutil-2.6.1.so
b640e000-b640f000 rw-p b640e000 00:00 0
b640f000-b6422000 r-xp 00000000 08:09 1855032 /lib/libnsl-2.6.1.so
b6422000-b6424000 rw-p 00012000 08:09 1855032 /lib/libnsl-2.6.1.so
b6424000-b6426000 rw-p b6424000 00:00 0
b6426000-b645b000 r-xp 00000000 08:0a 6358081 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-pal.so.0.0.0
b645b000-b645c000 r--p 00034000 08:0a 6358081 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-pal.so.0.0.0
b645c000-b645d000 rw-p 00035000 08:0a 6358081 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-pal.so.0.0.0
b645d000-b647f000 rw-p b645d000 00:00 0
b647f000-b64cf000 r-xp 00000000 08:0a 6358083 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-rte.so.0.0.0
b64cf000-b64d0000 r--p 00050000 08:0a 6358083 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-rte.so.0.0.0
b64d0000-b64d2000 rw-p 00051000 08:0a 6358083 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libopen-rte.so.0.0.0
b64d2000-b6571000 r-xp 00000000 08:0a 6214169 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESfilters.so
b6571000-b6573000 r--p 0009e000 08:0a 6214169 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESfilters.so
b6573000-b6574000 rw-p 000a0000 08:0a 6214169 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESfilters.so
b6574000-b6575000 rw-p b6574000 00:00 0
b6575000-b65e4000 r-xp 00000000 08:0a 6358070 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libmpi.so.0.0.0
b65e4000-b65e5000 r--p 0006e000 08:0a 6358070 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libmpi.so.0.0.0
b65e5000-b65ec000 rw-p 0006f000 08:0a 6358070 /home/openfoam1.5/OpenFOAM/ThirdParty/openmpi-1.2.6/platforms/linuxGccDPOpt/lib/libmpi.so.0.0.0
b65ec000-b65f2000 rw-p b65ec000 00:00 0
b65f2000-b6606000 r-xp 00000000 08:0a 6214141 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libz.so
b6606000-b6607000 r--p 00013000 08:0a 6214141 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libz.so
b6607000-b6608000 rw-p 00014000 08:0a 6214141 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libz.so
b6608000-b662e000 r-xp 00000000 08:0a 6214167 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESdeltas.so
b662e000-b662f000 r--p 00026000 08:0a 6214167 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESdeltas.so
b662f000-b6630000 rw-p 00027000 08:0a 6214167 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libLESdeltas.so
b6630000-b6734000 r-xp 00000000 08:0a 6214185 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libcompressibleLESModels.so
b6734000-b6737000 r--p 00103000 08:0a 6214185 /home/openfoam1.5/OpenFOAM/OpenFOAM-1.Aborted
openfoam1.5@spray:~/OpenFOAM/openfoam1.5-1.5/run/tutorials/dieselFoam/aachenBomb>

error is :

*** glibc detected *** dieselFoam: free(): invalid next size (fast): 0x0853b3f8 ***

or in some case:

*** glibc detected *** dieselFoam: malloc(): memory corruption (fast): 0x0b5a3a70 ***

I changed several injector model but none of them worked.
it seems that dieselfoam has a bug.
could any body help me?
I attach my injectorproperties file that I used.

please help me it is urgent.
thanks

az1362f

June 8, 2009 05:03

need help

any body know how can I fix it?

niklas

June 8, 2009 07:31

1 Attachment(s)

yup, there was a bug in the swirlInjector.

If you have different number of entries in the pressureProfile and massFlowRateProfile it wouldnt work.
I've attached a fix to solve this problem

save the unzipped file in src/lagrangian/dieselSpray/injector/swirlInjector/

then go to src/lagrangian/dieselSpray
and run 'wmake libso'

then you should be able to run dieselFoam

az1362f

June 10, 2009 11:18

need header file

1 Attachment(s)

Hello Dear Dr ;

I am thank you for your kindly attention.

I did every things that you said but I faced error when I try to compile it with wmake libso command,I though that I have to change swirlInjector.h to fix this errors,I tried to do it but I can not completely fixed it.

I want you to please send me swirlInjector.h file,then I can fix it.

regard

niklas

June 11, 2009 02:21

what on earth made you change the .h file??

I thought it was overwhelmingly obvious that you should save it as the .C file.

az1362f

June 11, 2009 07:18

change .h file to include class

I did not change swirlInjector.c file to swirlInjector.h file.
after I copied swirlInjector.c to swirlInjector folder I tried to compiled it with wmake libso but it generated some errors like these :

Making dependency list for source file injector/swirlInjector/swirlInjector.C
SOURCE=injector/swirlInjector/swirlInjector.C ; g++ -m32 -Dlinux -DDP -Wall -Wn
o-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoReposi
tory -ftemplate-depth-40 -I/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/finiteVol
ume/lnInclude -I/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/lagrangian/basic/lnI
nclude -I/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/turbulenceModels/RAS/compre
ssible/lnInclude -I/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/turbulenceModels/
LES/LESdeltas/lnInclude -I/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/turbulence
Models/LES/compressible/lnInclude -I/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/
thermophysicalModels/basic/lnInclude -I/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/s
rc/thermophysicalModels/liquids/lnInclude -I/home/openfoam1.5/OpenFOAM/OpenFOAM-
1.5/src/thermophysicalModels/liquidMixture/lnInclude -I/home/openfoam1.5/OpenFOA
M/OpenFOAM-1.5/src/thermophysicalModels/thermophysicalFunctions/lnInclude -I/hom
e/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/thermophysicalModels/specie/lnInclude -I
/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/thermophysicalModels/combustion/lnIn
clude -I/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/thermophysicalModels/pdfs/ln
Include -IlnInclude -I. -I/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/OpenFOAM/l
nInclude -I/home/openfoam1.5/OpenFOAM/OpenFOAM-1.5/src/OSspecific/Unix/lnInclude
-fPIC -pthread -c $SOURCE -o Make/linuxGccDPOpt/swirlInjector.o
injector/swirlInjector/swirlInjector.C:190: error: prototype for ‘const Foam::ve
ctor Foam::swirlInjector: position(Foam::label) const’ does not match any in cla
ss ‘Foam::swirlInjector’
injector/swirlInjector/swirlInjector.H:161: error: candidates are: virtual Foam:
:vector Foam::swirlInjector: position(Foam::scalar, bool, Foam::scalar, const Fo
am::vector&, const Foam::vector&, const Foam::vector&, Foam::Random&) const
injector/swirlInjector/swirlInjector.H:149: error: virtual const
Foam::vector Foam::swirlInjector: position() const
injector/swirlInjector/swirlInjector.C:205: error: prototype for ‘Foam::vector F oam::swirlInjector: position(Foam::label, Foam::scalar, bool, Foam::scalar, cons t Foam::vector&, const Foam::vector&, const Foam::vector&, Foam::Random&) const’ does not match any in class ‘Foam::swirlInjector’
injector/swirlInjector/swirlInjector.H:161: error: candidates are: virtual Foam: :vector Foam::swirlInjector: position(Foam::scalar, bool, Foam::scalar, const Fo am::vector&, const Foam::vector&, const Foam::vector&, Foam::Random&) const
injector/swirlInjector/swirlInjector.H:149: error: virtual const Foam::vector Foam::swirlInjector: position() const
injector/swirlInjector/swirlInjector.C:240: error: no ‘Foam::label Foam::swirlIn jector::nHoles() const’ member function declared in class ‘Foam::swirlInjector’
injector/swirlInjector/swirlInjector.C:254: error: prototype for ‘const Foam::ve ctor& Foam::swirlInjector::direction(Foam::label, Foam::scalar) const’ does not match any in class ‘Foam::swirlInjector’
injector/swirlInjector/swirlInjector.H:167: error: candidate is: virtual const F oam::vector& Foam::swirlInjector::direction() const
injector/swirlInjector/swirlInjector.C:383: error: no ‘Foam::vector Foam::swirlI njector::tan1(Foam::label) const’ member function declared in class ‘Foam::swirl Injector’
injector/swirlInjector/swirlInjector.C:388: error: no ‘Foam::vector Foam::swirlI njector::tan2(Foam::label) const’ member function declared in class ‘Foam::swirl Injector’
make: *** [Make/linuxGccDPOpt/swirlInjector.o] Error 1

I thought that these errors are generated because I did not redefine the swirlInjector.h file to include new members that you add in swirlInjector.c file like:

Foam::label Foam::swirlInjector::nHoles() const
{
return 1;
}

and it did not declare in swirlInjector.h then the compiler generated this error:

injector/swirlInjector/swirlInjector.C:240: error: no ‘Foam::label Foam::swirlIn jector::nHoles() const’ member function declared in class ‘Foam::swirlInjector’
injector/swirlInjector/swirlInjector.C

then because of this I tried to change the swirlInjector.h file to declare new function and data member, but it still generated errors(last attached file). I want you to send me a suggestion or swirlInjector.h file.

thanks
with regard

niklas

June 11, 2009 07:24

ah, ok.

You are using 1.5.
please update to the latest git 1.5.x version.
the correction will already be there then.

az1362f

June 11, 2009 07:44

where can I downloa it

thanks
Where can I download it?

az1362f

June 14, 2009 04:58

liquid tip penetration

1 Attachment(s)

hi

I implemented these line in dieselfoam solver to save liquid peneteration lenght in every time step,but my result showed an approximately a constant liquid penetration.(for example 10.xxxx)

std :: fstream penet
penet.open("penetration.txt",std::ios::app);
penet << dieselSpray.liquidPenetration(1.0)<<"\t"<<runTime. timeName()<<"\n"<<endl;
penet.close();

what is my mistake?
I used default setting of dieselfoam tutorials.
I attached a penetration result.
please help me.

az1362f

June 14, 2009 15:43

liquid & gas properties

hi

I solved my last problem.
I want to use a diesel fuel with density of 840 kg/m3 and nitrogen (N2) as a chamber gas.
how can I set this two parameters?
I checked C7H16 files (.C & .H) but I don't know some parameters,for example rho_(61.38396836, 0.26211, 540.2, 0.28141),I checked html document but it does not describe these parameters:

NSRDSfunc5 (scalar a, scalar b, scalar c, scalar d)

any body can describe them and say how can I get the description of another input parameters like
pv_(87.829, -6996.4, -9.8802, 7.2099e-06, 2),
hl_(540.20, 499121.791545248, 0.38795, 0, 0, 0),
cp_(540.20, 6.11976102401216, 3137.69909384855, 182.274175063868, -254.530511150515),etc.

can I set chamber gas by changing N2 and O2 files in 0 directory in dieselfoam tutorial directory?
how can I describe a chamber gas with a defined properties like density viscosity,etc?

please help me
regard

az1362f

June 15, 2009 17:13

gas and liquid property

hi

I solve my last problem but still I don't know how can I describe chamber gas?

for example for a chamber full of N2 as gas instead of air.

farbfilm

June 16, 2009 03:21

Hi,

did you check the 'thermophsicalProperties'?? There you can specify an inert gas!
I think, in the '0'-directory you can then define the composition of your gas (with the 'N2' or 'O2'-files etc.)!

niklas

June 16, 2009 06:53

Quote:

Originally Posted by az1362f (Post 219381)

hi

I solve my last problem but still I don't know how can I describe chamber gas?

for example for a chamber full of N2 as gas instead of air.

In the 0-directory you have the N2 and O2 files that contain the mass-fraction information for these species.
If you check them you will see that in the N2-file you have a line that reads like this -
internalField uniform 0.766;

and in the O2-file you have a line that reads like this -
internalField uniform 0.233;

These two species form 'air'.
I will leave it up to you to figure out how to modify it so you only have N2.

az1362f

June 16, 2009 16:20

initial condition

Hi

I know that these (N2 & O2 files in 0 directory) are initial condition for chamber gas.I am not sure about it when I change them to for example N2=1.0 and O2=0.0,is it like that chamber gas is only N2 in all of time steps?
But in thermophysicalProperties file I introduce chemkin file that include species like O2 N2 Co2 H2o,doesn't dieselfoam use these species as chamber gas for next time steps?
which part of code call for gas properties or where can I see gas properties in Src?

I change N2 to 1.0 and O2 to 0.0 but my result did n't effect any more(SMD & liquidpenetration).also I changed N2 to 1.0 and O2 to 1.0 but dieselfoam did not generate an error,How is it?
what is the role of inertSpecie in thermophysicalProperties file ?

Please guid me to select N2 just as a chamber gas.

thanks
regard

niklas

June 17, 2009 03:04

Quote:

Originally Posted by az1362f (Post 219517)

I know that these (N2 & O2 files in 0 directory) are initial condition for chamber gas.I am not sure about it when I change them to for example N2=1.0 and O2=0.0,is it like that chamber gas is only N2 in all of time steps?

No. It is only the initial conditions.
You have to specify boundary conditions and initial conditions for every specie
that you define in chem.inp, but to make life easier there is a file called Ydefault.
Lets say you have NO2, NO, CO, OH defined in the chem.inp.
Normally you would set the initial massfraction to zero for these, but if you have alot of species its quite exhausting just setting everything to zero, thats
where the Ydefault file comes in.
So if you havent got a file for NO2, NO, CO or OH the solver will look for a file called Ydefault and use the intial conditions and boundary conditions from that file.
Take a look in the file
src/thermophysicalModels/combustion/mixtureThermos/mixtures/combustionMixture/combustionMixture.C
if you want to see how it is done.

Quote:

Originally Posted by az1362f (Post 219517)

But in thermophysicalProperties file I introduce chemkin file that include species like O2 N2 Co2 H2o,doesn't dieselfoam use these species as chamber gas for next time steps?
which part of code call for gas properties or where can I see gas properties in Src?

the properties are defined in therm.dat.
all the species will have the laminar viscosity-coefficients hardcoded for the sutherland expression to 1.67212e-6 and 170.672
this is done in
src/thermophysicalModels/combustion/chemistryReaders/chemkinReader/chemkinLexer.L
You can look here for more info on format
http://openfoamwiki.net/index.php/Contrib_dieselFoam

Quote:

Originally Posted by az1362f (Post 219517)

I change N2 to 1.0 and O2 to 0.0 but my result did n't effect any more(SMD & liquidpenetration).also I changed N2 to 1.0 and O2 to 1.0 but dieselfoam did not generate an error,How is it?

I fail to see why changing the gas to just N2 would affect these parameters.
Usually you use N2 in experiments to avoid combustion, but in the CFD world you can just set the chemistry to off in chemistryProperties to turn off all the reactions.
When the massfractions are read everything will be normalized to 1. Just run a few timesteps and check the values and you will see that they will have changed to 0.5.

Quote:

Originally Posted by az1362f (Post 219517)

what is the role of inertSpecie in thermophysicalProperties file ?

since the massfractions have to always add up to 1 you dont have to solve the transport equations for all the species,
the inertSpecie will instead be calculated as 1- sum(Yi)
check the file
applications/solvers/combustion/dieselEngineFoam/YEqn.H

az1362f

June 18, 2009 01:55

diesel Fuel

Hi

thank you very much for your complete answers dear Dr Nordin.

I found that fuel properties that I want to use are close to diesel fuel.does openfoam has properties of these fuel?
I check in SRC and found C12H26,is this fuel used as diesel fuel in common?
I now that C12H26 is one of the component of diesel fuel. is it correct to use it as diesel fuel?

I used C12H26 in thermophysical properties file and chemkin file as second fuel But when I ran diesel fuel I face an error, I check therm.dat file in chemkin directory of case but I didn't find C12H26 coefficient in this file,I am sure that this error is related to this file.
then how can i find these coefficient?
Do you have these coefficient.I will be thankful if you send them for me.

Thanks
regard

niklas

June 22, 2009 06:07

http://garfield.chem.elte.hu/Burcat/burcat.html

az1362f

June 22, 2009 21:25

nasa polynomial

Hi Dear Dr Nordin;

I am thankful for your kindly help.I am so glad that have a friend like you.

I also found a pdf like this:

http://www.docstoc.com/docs/4201307/...ondensed-Phase

Thanks
Regard

az1362f

June 22, 2009 21:29

Kiva Or openFoam

Hi Dear Dr Nordin;

I know that you work with both KIva and openFoam code,I have a question:

Which of them are really suitable for simulating spray simulation?

niklas

June 23, 2009 02:00

Quote:

Originally Posted by az1362f (Post 220123)

Hi Dear Dr Nordin;

I know that you work with both KIva and openFoam code,I have a question:

Which of them are really suitable for simulating spray simulation?

both,
but if you are planning to alter/add something on your own
I know which one I would pick.

az1362f

June 23, 2009 04:31

kiva or openfoam

Hi Dear Dr Nordin;

first by following my last question I am planning to alter/add something on my own,which one I would pick?

second I copy new therm.dat file that you putted for me,but I face this error :

| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : dieselFoam
Date : Jun 23 2009
Time : 22:35:57
Host : spray
PID : 8865
Case : /home/openfoam1.5/OpenFOAM/openfoam1.5-1.5/run/tutorials/dieselFoam/aachenBomb
nProcs : 1

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

Create mesh for time = 0

Reading thermophysicalProperties
Selecting thermodynamics package hMixtureThermo<reactingMixture>
Selecting chemistryReader chemkinReader

C12H26 not found in table. Valid entries are
820
(
CH3CHCL
CL2SISICL
CH3CHCH
SIF2N
A2C2H>2
CH3OCL
CH2CCLOH
.
.
.
.
.
CCL2CCLO
H2SINH3
KO2<S>
A1C2H3*
A1C2H>2
CLCH2OH
H2SISIH2
CCL3CCLO
H3SICH2
H3SICH3
)

in which table I would define C12H26?

thanks
regard

niklas

June 23, 2009 04:55

Its clearly telling you that it cannot find C12H26 in therm.dat
thats because in Burcat's file the name is N-DODECANE

Code:

C12H26            T 5/99C  12H  26    0    0G   200.000  6000.000  1000.000    1

3.70187925E+01  5.54721488E-02-1.92079548E-05 3.08175574E-09-1.84800617E-13    2

-5.26984458E+04-1.61453501E+02 2.13264480E+01-3.86394002E-02 3.99476113E-04    3

-5.06681097E-07 2.00697878E-10-4.22475053E+04-4.85848300E+01-3.49836226E+04    4

...and I would of course use openfoam.

JMP	June 23, 2009 10:35

saving liquid spray penetration

Hi,

I'm a begginer into OpenFoam and C++. Maybe this is an obvious question.
In order to write liquid penetration, did you include those lines into dieselFoam.c or did you create a *.h such as spraySummary.h?
Is required any additional coding in order to use 'ofstream' into dieselfoam.c ,e.g.
#include <iostream>
#include <fstream>
using namespace std;?
Thanks.

Quote:

Originally Posted by az1362f (Post 219236)

az1362f

June 23, 2009 11:23

liquid peneteration

Hello Dear Pastor;

of course you can add these line in both place in spraysummary.h or dieslFoam.c,because dieselFoam call this *.H file in main program but I directly wrote these lines to dieselfoam.c after line :Incloude" spraysummary.h" or you can copy paste the lines inside spraysummary.h to dieselfoam.c and then write this code to dieselfoam.c.
second when you usestd::ios::app, it is enough, you don't need to write using namespace std;

regard

az1362f

June 23, 2009 11:27

C12h26

Hi Dear Dr Nordin

No,I changed Therm.dat(N-DODECANE ->C12H26) before I ran dieselfoam,but it generated this error.
what is the problem?

regard

az1362f

June 23, 2009 11:28

C12h26

Hi Dear Dr Nordin

No,I changed Therm.dat(N-DODECANE ->C12H26) before I ran dieselfoam,but it generated this error.
what is the problem?

regard

JMP	June 23, 2009 11:41

Hi John,

Thanks for your answer. I also add those lines after 'include "spraySummary.H" ' but does not work for me. I get the following compilation error

dieselFoam.C: In function âint main(int, char**)â:
dieselFoam.C:117: error: expected initializer before âstdâ
dieselFoam.C:119: error: âpenetâ no se declarÃ³ en este Ã¡mbito
make: *** [Make/linux64GccDPOpt/dieselFoam.o] Error 1

I'll keep trying...

Quote:

Originally Posted by az1362f (Post 220226)

niklas

June 24, 2009 02:14

Quote:

Originally Posted by az1362f (Post 220229)

Hi Dear Dr Nordin

No,I changed Therm.dat(N-DODECANE ->C12H26) before I ran dieselfoam,but it generated this error.
what is the problem?

regard

You are probably picking up the wrong therm.dat file.
Check thermophysicalProperties and see which file you are using.

CHEMKINThermoFile "~OpenFOAM/thermoData/therm.dat";

az1362f

June 25, 2009 04:29

validating exprimental case

2 Attachment(s)

Hi Dear Dr Nordin;

Thanks.
yes,I make a mistake.
Now I am trying to validate a non-evaporating spray of C7H16 experiment with these characteristic:
max.injection pressure= 60 Mpa
injection duration=2.36 ms
mass of fuel= 9.84 mg
Diameter of the nozzle=0.02 cm
ambient temperature = 293 k
ambient pressure =17 bar

But my result has a large difference with experiment.
I attached a zip file including my diselfoam case setup, my tip penetration results and experiment results.

it seems that spray can not penetrate in chamber,I changed ambient pressure of 17 bar(17e+5) to (1.7e+5) and strangely my result of tip penetration were close to experiment.

I don't know what is my problem.
it seems that ambient pressure should be 1.7e+5.

please help me,I am confused.

Thanks
Regard

niklas

June 25, 2009 09:44

you cant run a 2D simulation like that.

az1362f

June 25, 2009 09:53

correct mesh

hi

what is my mistake?
what is its correction?
how can i describe an axi-symmetric mesh for diesel foam?

thanks
regard

az1362f

June 26, 2009 04:40

3D mesh

1 Attachment(s)

hi Dear Niklas

what is my mistake?
what is its correction?
how can i describe an axi-symmetric mesh for diesel foam?

I ran a 3D simulation of dieselfoam with last characteristic of n-heptane experiment and add its result to a new plot and attached it.

there are two strange things:
1- I saw some fluctuations in liquid penetration curve in 3D cases.
2- likewise liquid penetration is far from experimental result.
what is the reason of these?

I set the initial drop diameter to nozzle diameter according to injector model.is it correct?

I am confused.
Please help me.

thanks
regard

All times are GMT -4. The time now is 20:29.

Page 1 of 2

Show 40 post(s) from this thread on one page