[hjasak]

Have a look at foamCFD web site: I am trying to keep a repository of FOAM PhD Theses and other work I can get my hands on + some additional info.

The problem with most other algorithms are that they haven't been published on their own and various PhD theses are typically the best source of data even if they are a bit out of date. The only reliable and up-to-date stuff is the actual code itself which you get in full source :-)

As for my other work with FOAM, there are references on my private web site. I have tried to keep a lits of all papers published using FOAM, but because of all the other work I need to do this is pretty outdated.

Enjoy,

Hrv

[grtabor]

Quite a lot of info about an LES version of Henry's combustion model is given in a paper in Flow, Turbulence and Combustion :

Large Eddy Simulation of Premixed Turbulent Combustion Using Xgr Flame Surface Wrinkling Model
G. Tabor and H.G. Weller, Flow, Turbulence and Combustion, Volume 72, Number 1 Date: March 2004
Pages: 1 - 27

http://springerlink.metapress.com/(s1e0qlbrkhvn5x451ar0bqrw)/app/home/contribution.asp?referrer=parent&backto=issu e,1,4;journal,8,46;linkingpublicationresults,1:100 237,1
http://www.springeronline.com/sgw/cda/frontpage/0,0,3-0-70-35731861-0,0.html

Other info in a couple of Combustion Symposium papers :

Application of a Flame-Wrinkling LES Combustion Model to a Turbulent Mixing Layer
H.G.Weller, G.Tabor, A.D.Gosman & C.Fureby
In: Proceedings of the 27th Combustion Symposium pp. 899 -- 907 (1998)

Measurements and Large Eddy Simulations of Turbulent Premixed Flame Kernel Growth
I.K.Nwagwe, H.G.Weller, G.Tabor, A.D.Gosman, M.Lawes, C.G.W.Sheppard & R.Wooley
In: Proceedings of the 28th Combustion Symposium (2000)

Also the RANS version is detailed in an IC internal report :


@techreport{Weller:1993,
author = "Weller, H.G.",
title = "{T}he {D}evelopment of a {N}ew {F}lame {A}rea
{C}ombustion {M}odel {U}sing
{C}onditional {A}veraging",
institution = "Imperial College of Science,
Technology and Medicine",
year = 1993,
month = "March",
type = "Thermo-Fluids Section Report",
number = "TF 9307"
}

if this is any use.

I'm in the process of implementing some other combustion models as well - I will pass them around when I'm sure they are working.

Gavin

[pratap]

Hi,

I have question whether in wellers model implemented in XiFoam is it possible to get a flame when you dont specify any ingition but specify the boundary condition of ctilde of 1 and 0 in the outlet and inlet.
This approach is commnly used in other cfd codes. If I want to run a case where I dont want to solve enthalpy euqation (adiabatic conditions) do I still to include a thermo model it is because viscosity is calculated from it.


Gavin you have asked this kind of question earlier did you make that thing to work.

Pratap

[lucchini]

Hi Pratap !!!!!!
yes you can of course, but you have to modify the XiFoam solver as follows:

comment the following line in the bEqn.H

if(ign.ignited())

and then recompile the XiFoam code with the wmake command.
Then you can run XiFoam as you like
Bye
Tommaso

[pratap]

Thanks for your response,

However I have more problems if you consider the bEqn.h

The equation reads like this

fvScalarMatrix bEqn
(
fvm::ddt(rho, b)
+ mvConvection->fvmDiv(phi, b)
+ fvm::div(phiSt, b, "div(phiSt,b)")
- fvm::Sp(fvc::div(phiSt), b)
- fvm::laplacian(turbulence->muEff(), b)
);

The term which responsible for the reaction rate is the follwing
fvm::div(phiSt, b, "div(phiSt,b)")
- fvm::Sp(fvc::div(phiSt), b)

The wellers model close of the reaction rate is
phiSt*|grad(b)|. So does this means that

fvm::div(phiSt, b, "div(phiSt,b)")
- fvm::Sp(fvc::div(phiSt), b)=phiSt*|grad(b)|

If this the case can you tell me why do we need to decompose the phiSt*|grad(b)| into a im0plicit source term fvm::Sp(fvc::div(phiSt) and the other term.

I also have a general question. When I make my new class application do I have to compile on 32 and 64 bit machine separately before I can use them on 64 and 32 bit machine.

Thanks in advance

Pratap

[vijayakumar]

hi i am new to openfoam, i am solving one combustion problem.. block mesh dict is shown below.. i used XIFOAM as solver... while running the case the errors are

Creating turbulence model

Selecting turbulence model kEpsilon
#0 Foam::error::printStack(Foam:stream&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::compressible::turbulenceModels::kEpsilon::kE psilon(Foam::GeometricField<do uble,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#4 Foam::compressible::turbulenceModel::adddictionary ConstructorToTable<foam::compr essible::turbulencemodels::kepsilon>::New(Foam::Ge ometricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#5 Foam::compressible::turbulenceModel::New(Foam::Geo metricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#6 main in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"

vertices
(
(0 0 0)
(0 0.5 0)
(14 0 0)
(14 0.5 0)
(180 0 0)
(180 0.5 0)
(0 22.5 0)
(180 22.5 0)
(0 0 1)
(0 0.5 1)
(14 0 1)
(14 0.5 1)
(180 0 1)
(180 0.5 1)
(0 22.5 1)
(180 22.5 1)
);

blocks
(
hex (0 2 3 1 8 10 11 9) (14 1 1) simpleGrading (1 1 1)
hex (2 4 5 3 10 12 13 11) (166 1 1) simpleGrading (1 1 1)
hex (1 5 7 6 9 13 15 14) (180 22 1) simpleGrading (1 1 1)
);

edges
(
);

patches
(
patch inlet1
(
(1 9 14 6)
)
patch inlet2
(
(0 8 9 1)
)
patch outlet
(
(5 13 15 7)
)
wall topwall
(
(6 14 15 7)
(1 9 11 3)
)
wall bottomwall
(
(0 8 10 2)
(2 10 12 4)
)
empty frontandback
(
(8 9 11 10)
(0 1 3 2)
(9 14 15 13)
(1 6 7 5)
(10 11 13 12)
(2 3 5 4)
)
);

mergePatchPairs
(
);
please can u tell where i am going wrong... i am doing methane combustion.

[vijayakumar]

hi i am new to openfoam, i am solving one combustion problem.. block mesh dict is shown below.. i used XIFOAM as solver... while running the case the errors are

Creating turbulence model

Selecting turbulence model kEpsilon
#0 Foam::error::printStack(Foam:stream&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::compressible::turbulenceModels::kEpsilon::kE psilon(Foam::GeometricField<do uble,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#4 Foam::compressible::turbulenceModel::adddictionary ConstructorToTable<foam::compr essible::turbulencemodels::kepsilon>::New(Foam::Ge ometricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#5 Foam::compressible::turbulenceModel::New(Foam::Geo metricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#6 main in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"

vertices
(
(0 0 0)
(0 0.5 0)
(14 0 0)
(14 0.5 0)
(180 0 0)
(180 0.5 0)
(0 22.5 0)
(180 22.5 0)
(0 0 1)
(0 0.5 1)
(14 0 1)
(14 0.5 1)
(180 0 1)
(180 0.5 1)
(0 22.5 1)
(180 22.5 1)
);

blocks
(
hex (0 2 3 1 8 10 11 9) (14 1 1) simpleGrading (1 1 1)
hex (2 4 5 3 10 12 13 11) (166 1 1) simpleGrading (1 1 1)
hex (1 5 7 6 9 13 15 14) (180 22 1) simpleGrading (1 1 1)
);

edges
(
);

patches
(
patch inlet1
(
(1 9 14 6)
)
patch inlet2
(
(0 8 9 1)
)
patch outlet
(
(5 13 15 7)
)
wall topwall
(
(6 14 15 7)
(1 9 11 3)
)
wall bottomwall
(
(0 8 10 2)
(2 10 12 4)
)
empty frontandback
(
(8 9 11 10)
(0 1 3 2)
(9 14 15 13)
(1 6 7 5)
(10 11 13 12)
(2 3 5 4)
)
);

mergePatchPairs
(
);
please can u tell where i am going wrong... i am doing methane combustion.

[gschaider]

Hi Vijayakumar!

Almost certainly the problem is that you have intial and/or boundary conditions with the value 0 for epsilon and/or k

Something that is also certain is that nobody will answer your questions anymore if you continue to post the same question to multiple (possibly unrelated) threads

Bernhard

PS: on the positive side: you provided a stack-trace, which made the diagnosis of the problem possible. So things are improving. If you could only try to post a question to only 1 thread ....

[vijayakumar]

hi Bernhard sir i will post my problem to this thread itself. i solved this problem in fluent, but in foam from SIX weeks i am trying to solve the above combustion problem, in XIFOAM solver some field variables i can't define, how these variables to be defined (P, T,k, epsilon, tu, ft, fu, b, xi, su) my inlet conditions are... air inlet 0.5 m/s at 300k. fuel inlet 80m/s at 300k, wall tempr is 300k,
i am defining k like this
k = 3/2*(0.5)2 = 0.375 m2 /s2 for air inlet

and epsilon as
&epsilon;=Cµ^0.75*k^1.5/L

i am giving outlet pr as 0pa
i am defining u as 0.5m/s for airinlet and 80 m/s for fuel inlet ... while running the case it is showing error message
Creating turbulence model

Selecting turbulence model kEpsilon
#0 Foam::error::printStack(Foam::-Ostream&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::compressible::turbulenceModels::kEpsilon::kE psilon(Foam::GeometricField<do > const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#4 Foam::compressible::turbulenceModel::adddictionary ConstructorToTable<foam::compr >::New(Foam::GeometricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#5 Foam::compressible::turbulenceModel::New(Foam::Geo metricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#6 main in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"


please tell where i am going wrong, how to correct it.

[joe_]

Hello everybody!

I am working with the engineFoam tutorial trying to simulate the ignition process. My problem is about these ignitionSites in the constant/combustionProperties file.
For each ignitionSite the user has to set values for:

location
diameter
start
duration
strength

I understand the meaning and what dimensions are used in the first four, but not the one with strength. What kind of measure is it? Energy? What dimension?
I really need this information and I have been searching through the forum. I have also read "The Development of a New Flame Area Combustion Model Using Conditional Averaging" report by Henry Weller that this model is built on without success.

So if you know anything about it or where I can find information I would be very happy.

[gschaider]

Hi Vijayakumar!

I have zero experience with XiFoam so I can only guess here.

Your inlet k-epsilon at the inlet seems fine, but check whether k or epsilon are 0 somewhere else (on the internalField or some othere boundary). Usually that is the problem when the k-eps fails. Fluent protects you from such stuff but OpenFOAM does not (it is CFD for responsible adults)

On the other hand: I was on the impression that XiFoam is a compressible solver. So I'm amazed that it didn't fail while creating the thermophysical model with an absolute pressure of 0 on some boundary

Bernhard

[vijayakumar]

hi i am new to openfoam, i am solving one combustion problem.. block mesh dict is shown below.. i used XIFOAM as solver... while running the case the errors are coming


vertices
(
(0 0 0)
(0 0.5 0)
(14 0 0)
(14 0.5 0)
(180 0 0)
(180 0.5 0)
(0 22.5 0)
(180 22.5 0)
(0 0 1)
(0 0.5 1)
(14 0 1)
(14 0.5 1)
(180 0 1)
(180 0.5 1)
(0 22.5 1)
(180 22.5 1)
);

blocks
(
hex (0 2 3 1 8 10 11 9) (14 1 1) simpleGrading (1 1 1)
hex (2 4 5 3 10 12 13 11) (166 1 1) simpleGrading (1 1 1)
hex (1 5 7 6 9 13 15 14) (180 22 1) simpleGrading (1 1 1)
);

edges
(
);

patches
(
patch inlet1
(
(1 9 14 6)
)
patch inlet2
(
(0 8 9 1)
)
patch outlet
(
(5 13 15 7)
)
wall topwall
(
(6 14 15 7)
(1 9 11 3)
)
wall bottomwall
(
(0 8 10 2)
(2 10 12 4)
)
empty frontandback
(
(8 9 11 10)
(0 1 3 2)
(9 14 15 13)
(1 6 7 5)
(10 11 13 12)
(2 3 5 4)
)
);

mergePatchPairs
(
);
please can u tell where i am going wrong... i am doing methane combustion.
i am trying to solve the above combustion problem, in XIFOAM solver some field variables i can't define, how these variables to be defined (P, T,k, epsilon, tu, ft, fu, b, xi, su) my inlet conditions are... air inlet 0.5 m/s at 300k. fuel inlet 80m/s at 300k, wall tempr is 300k,
i am defining k like this
k = 3/2*(0.5)2 = 0.375 m2 /s2 for air inlet

and epsilon as
&epsilon;=Cµ^0.75*k^1.5/L

i am giving outlet pr as 0pa
i am defining u as 0.5m/s for airinlet and 80 m/s for fuel inlet ... while running the case it is showing error message


Foam::error::printStack(Foam:stream&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::compressible::turbulenceModels::kEpsilon::kE psilon(Foam::GeometricField<do uble,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#4 Foam::compressible::turbulenceModel::adddictionary ConstructorToTable<foam::compr essible::turbulencemodels::kepsilon>::New(Foam::Ge ometricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#5 Foam::compressible::turbulenceModel::New(Foam::Geo metricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleTur bulenceModels.so"
#6 main in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"

[viji]

hi
how and were to define boundary condition of separate inlets for air and fuel in a combustion problem.. i am working in XIfoam

[viji]

Exec : XiFoam /home/openfoam14/OpenFOAM/openfoam14-1.4.1/run/OPprashant 111222333
Date : Feb 25 2009
Time : 16:42:04
Host : shani
PID : 5828
Root : /home/openfoam14/OpenFOAM/openfoam14-1.4.1/run/OPprashant
Case : 111222333
Nprocs : 1
Create time

Create mesh for time = 0

Reading combustion properties

Found ignition cells:

214
(
13
3
4
5
6
7
8
9
10
11
12
14
15
16
17
18
19
20
21
22
23
24
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
1278
1279
1280
1281
1282
1446
1447
1448
1449
1450
1451
1452
1453
1454
1455
1456
1457
1458
1459
1460
1461
1627
1628
1629
1630
1631
1632
1633
1634
1635
1636
1637
1638
1639
1640
1808
1809
1810
1811
1812
1813
1814
1815
1816
1817
1818
1819
1991
1992
1993
1994
1995
1996
)


Ignition on

Reading environmentalProperties
Reading thermophysical properties

Selecting thermodynamics package hhuMixtureThermo<homogeneousmixture<consttransport <speciethermo<hconstthermo<per fectgas>>>>>
min(b) = 1

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model kEpsilon
Creating field DpDt

#0 Foam::error::printStack(Foam:stream&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<foam::fvspatchfield,>(Foam::Geometric Field<double,>&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"
#5 Foam::tmp<foam::geometricfield<double,> > Foam::operator/<foam::fvspatchfield,>(Foam::GeometricField<double ,> const&, Foam::tmp<foam::geometricfield<double,> > const&) in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"
#6 main in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/openfoam14/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/XiFoa m"

i am solving combustion problem using XIFoam solver the above error is coming, how to overcome that error....