CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

FGM combustion model

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree14Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   June 5, 2014, 11:26
Default
  #41
Member
 
Likun
Join Date: Feb 2013
Posts: 52
Rep Power: 5
Likun is on a distinguished road
Send a message via Skype™ to Likun
Quote:
Originally Posted by dongchao yang View Post
Hi, all
Is there any one using this code?
The links of LTTRostockExtensions above seem invalid now...
I'm intrested in FMG model, and want to get the code mentioned above.
Hope some one can help me.

yangdongchao@gmail.com

Best Regards
Dongchao
Are you still interested in this code? I have a copy, if you want I can share it with you.

Best regards,
Likun
Likun is offline   Reply With Quote

Old   June 11, 2014, 04:45
Default
  #42
Member
 
Likun
Join Date: Feb 2013
Posts: 52
Rep Power: 5
Likun is on a distinguished road
Send a message via Skype™ to Likun
Quote:
Originally Posted by Likun View Post
Hi Foamers,

Has anyone succeed in compiling the FGM models in OF 2.3? I tried to do this, but there are some errors:

In file included from lnInclude/SIJPDFthermo.H:314,
from thermophysicalModels/compositionPdf/abstractProgressVariableThermo/abstractProgressVariableThermo.C:31:
lnInclude/SIJPDFthermo.C: In member function 'virtual void Foam::SIJPDFthermo<IndexDriver>::correct()':
lnInclude/SIJPDFthermo.C:1170: error: 'psi_' was not declared in this scope
lnInclude/SIJPDFthermo.C:1206: error: 'mu_' was not declared in this scope
lnInclude/SIJPDFthermo.C:1247: error: 'mu_' was not declared in this scope
lnInclude/SIJPDFthermo.C: In member function 'bool Foam::SIJPDFthermo<IndexDriver>::readPreIntegrated TableFromCache()':
lnInclude/SIJPDFthermo.C:1446: error: 'psi_' was not declared in this scope
lnInclude/SIJPDFthermo.C: In member function 'void Foam::SIJPDFthermo<IndexDriver>::writeCache() const':
lnInclude/SIJPDFthermo.C:1492: error: 'psi_' was not declared in this scope

can anybody help? Where are these errors come from? Is it because of the different OF version, because this CombustionAddons was previously using in OF1.7.
I got a step further, the above errors are because the differences of the basicThermo in the different OF versions. In the OF1.7 the compressibility 'psi' and viscosity 'mu' are defined in the basicThermo, but in OF2.3 they are defined in another file named psiThermo. I could solve the above error by including the header file 'psiThermo.H' in to 'abstractProgressVariableThermo.H' where the 'psi' and 'mu' are required. But the difference on the thermophysical model in the different OpenFOAM versions cause a lot of other errors, I shall try work them out.

Best,
Likun
Likun is offline   Reply With Quote

Old   August 4, 2014, 22:55
Default
  #43
hzw
New Member
 
Join Date: May 2014
Posts: 2
Rep Power: 0
hzw is on a distinguished road
Quote:
Originally Posted by hannes View Post
Hello,

I have already implemented this. The whole code is here:

http://janus.fms.uni-rostock.de/view...ompositionPdf/

Best regards,
Hannes
Hi, Hannes

I would be very appreciated if you can send me the source code of the FGM combustion model implemented in OpenFOAM, as i cannot get into the website you provided above for unknown server error.

Best regards!
hzw is offline   Reply With Quote

Old   August 5, 2014, 02:00
Default
  #44
Member
 
Hannes Kröger
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 96
Rep Power: 10
hannes is on a distinguished road
Dear hzw,

the server "janus" you tried is offline for a while now.
Try to clone the code from the git repo this sourceforge project instead:
http://sf.net/p/ofca

i.e. by
$ cd ~/OpenFOAM
$ git clone http://git.code.sf.net/p/ofca/code CombustionAddons

Regards, Hannes
hannes is offline   Reply With Quote

Old   August 5, 2014, 07:45
Default
  #45
hzw
New Member
 
Join Date: May 2014
Posts: 2
Rep Power: 0
hzw is on a distinguished road
Dear Hannes,

I've got what i need, thank you very much for your hard work. And please do not get bothered if I had any problem to consult.

Best Regards!
hzw is offline   Reply With Quote

Old   August 5, 2014, 11:11
Default
  #46
Member
 
Likun
Join Date: Feb 2013
Posts: 52
Rep Power: 5
Likun is on a distinguished road
Send a message via Skype™ to Likun
Dear Hannes,

Do you have an implementation of FGM model in the new version of OpenFoam (OpeFoam-2.3.0) or in a more recent version than OF-1.6?

Best regrds,
Likun
Likun is offline   Reply With Quote

Old   September 2, 2014, 05:27
Smile
  #47
New Member
 
dongchao yang
Join Date: Oct 2010
Posts: 6
Rep Power: 7
dongchao yang is on a distinguished road
Quote:
Originally Posted by Likun View Post
Are you still interested in this code? I have a copy, if you want I can share it with you.

Best regards,
Likun
Yes, I am still interested in FGM, hope you can send me a copy. Thank you


yangdongchao500@126.com
dongchao yang is offline   Reply With Quote

Old   September 2, 2014, 06:24
Default
  #48
Member
 
Likun
Join Date: Feb 2013
Posts: 52
Rep Power: 5
Likun is on a distinguished road
Send a message via Skype™ to Likun
Quote:
Originally Posted by dongchao yang View Post
Yes, I am still interested in FGM, hope you can send me a copy. Thank you


yangdongchao500@126.com

Already sent!
Likun is offline   Reply With Quote

Old   February 23, 2015, 11:38
Default
  #49
Member
 
Join Date: Feb 2015
Posts: 30
Rep Power: 3
Stefano Puggelli is on a distinguished road
Hi Likun,
I would be very interested in the FGM model for OF 2.3. It would be possible to have a copy of it?
Stefano Puggelli is offline   Reply With Quote

Old   February 24, 2015, 04:07
Default
  #50
Member
 
Likun
Join Date: Feb 2013
Posts: 52
Rep Power: 5
Likun is on a distinguished road
Send a message via Skype™ to Likun
Quote:
Originally Posted by hannes View Post
Dear hzw,

the server "janus" you tried is offline for a while now.
Try to clone the code from the git repo this sourceforge project instead:
http://sf.net/p/ofca

i.e. by
$ cd ~/OpenFOAM
$ git clone http://git.code.sf.net/p/ofca/code CombustionAddons

Regards, Hannes
Hi Stefano,

You can find the code from the links provided by Hannes. But this code is for very old OF version (1.6 or 1.7), and does not work in the new version.

Best,
Likun
Likun is offline   Reply With Quote

Old   March 3, 2015, 10:53
Default CombustionAddon's example question
  #51
New Member
 
xiao yinli
Join Date: Feb 2015
Posts: 1
Rep Power: 0
gravityxiao is on a distinguished road
Dear foamers,
I downloaded the code ofca of Hannes Kroege and built all the solvers with OF1.7.0. All the examples are ok except the SydneyHM1. There are some errors when I run the SydneyHM1 with edcSimpleFoam. I've two questions:
1. In controlDict, there is a line:
libs("libVAPExtensions.so");
what's it for?
When i run it with this line there is a warning: no such file or directory. I cannot find this file in anywhere.
2. The error info is : rhoMax is undefined, when i added the following line:
rhoMax rhoMax [1 -3 0 0 0 0 0]
rhoMin rhoMin [1 -3 0 0 0 0 0]
the error info become :
LHS and RHS + have different dimensions.
Is something wrong with the source code?
gravityxiao is offline   Reply With Quote

Old   March 6, 2015, 09:44
Default
  #52
Member
 
Join Date: Feb 2015
Posts: 30
Rep Power: 3
Stefano Puggelli is on a distinguished road
Hi,
I have compiled all the library and also the solver that are included into the code. However I am having some problems with the compilation of generateILDM. I have installed both cantera 1.7 and octave 3.2.4 but when I try to compile I recive the following error:



/usr/local/include/octave-3.2.4/octave/MArray.h:58: undefined reference to `Array<double>::~Array()'
/usr/local/include/octave-3.2.4/octave/MArray.h:58: undefined reference to `Array<double>::~Array()'
/usr/local/include/octave-3.2.4/octave/MArray.h:58: undefined reference to `Array<double>::~Array()'
Make/linux64GccDPOpt/generationLoop.o: In function `Foam::generationLoop::appendTableDescriptions(Foa m::ListDescription&, Foam::ListDescription&)':
/home/user/OpenFOAM/OpenFOAM_extensions/CombustionAddons/applications/generateILDM/generationLoop.C:290: undefined reference to `Foam::ListDescription::appendEntry(Foam::word const&)'
Make/linux64GccDPOpt/generationLoop.o: In function `~MArray':
/usr/local/include/octave-3.2.4/octave/MArray.h:58: undefined reference to `Array<double>::~Array()'
/usr/local/include/octave-3.2.4/octave/MArray.h:58: undefined reference to `Array<double>::~Array()'
/usr/local/include/octave-3.2.4/octave/MArray.h:58: undefined reference to `Array<double>::~Array()'
/usr/local/include/octave-3.2.4/octave/MArray.h:58: undefined reference to `Array<double>::~Array()'
/usr/local/include/octave-3.2.4/octave/MArray.h:58: undefined reference to `Array<double>::~Array()'
Make/linux64GccDPOpt/main.o.data.rel.ro._ZTVN7Cantera11IdealGasMixE[vtable for Cantera::IdealGasMix]+0x20): undefined reference to `Cantera::State::setMoleFractions(double const*)'
Make/linux64GccDPOpt/main.o.data.rel.ro._ZTVN7Cantera11IdealGasMixE[vtable for Cantera::IdealGasMix]+0x28): undefined reference to `Cantera::State::setMoleFractions_NoNorm(double const*)'
Make/linux64GccDPOpt/main.o.data.rel.ro._ZTVN7Cantera11IdealGasMixE[vtable for Cantera::IdealGasMix]+0x30): undefined reference to `Cantera::State::setMassFractions(double const*)'
Make/linux64GccDPOpt/main.o.data.rel.ro._ZTVN7Cantera11IdealGasMixE[vtable for Cantera::IdealGasMix]+0x38): undefined reference to `Cantera::State::setMassFractions_NoNorm(double const*)'
Make/linux64GccDPOpt/main.o.data.rel.ro._ZTVN7Cantera11IdealGasMixE[vtable for Cantera::IdealGasMix]+0x40): undefined reference to `Cantera::State::setConcentrations(double const*)'
Make/linux64GccDPOpt/main.o.data.rel.ro._ZTVN7Cantera11IdealGasMixE[vtable for Cantera::IdealGasMix]+0x58): undefined reference to `Cantera::Phase::freezeSpecies()'
Make/linux64GccDPOpt/main.o.data.rel.ro._ZTVN7Cantera11IdealGasMixE[vtable for Cantera::IdealGasMix]+0x68): undefined reference to `Cantera::ThermoPhase::duplMyselfAsThermoPhase()'
Make/linux64GccDPOpt/main.o.data.rel.ro._ZTVN7Cantera11IdealGasMixE[vtable for Cantera::IdealGasMix]+0xd0): undefined reference to `Cantera::ThermoPhase::activityConvention() const'
Make/linux64GccDPOpt/main.o.data.rel.ro._ZTVN7Cantera11IdealGasMixE[vtable for Cantera::IdealGasMix]+0xf0): undefined reference to `Cantera::ThermoPhase::getUnitsStandardConc(double *, int, int)'


I would like to know if someone has some ideas to solve these compilation problems.

Best regards
Stefano Puggelli is offline   Reply With Quote

Old   March 12, 2015, 04:43
Default
  #53
Member
 
Join Date: Feb 2015
Posts: 30
Rep Power: 3
Stefano Puggelli is on a distinguished road
Ok finally I am able to compile the whole solver also with generate ILDM ( I don't see the generateFGM that is mentioned in some previous posts... maybe it has been included in the generateILDM but I would like to test it). However when I try to run the generateILDM for the moment I recive the same error posted by Jose:

Quote:
Originally Posted by jose_rodrig View Post
Hi

Ok, I was able to compile finally both generateFGM and generateILDM.

However, when I run either I get the message:

Code:
jrodrig@pup:~/OpenFOAM-1.5.x/src/LTTRostockExtensions/applications/generateILDM/GRI1.2$ generateILDM 
terminate called after throwing an instance of 'std::logic_error'
  what():  basic_string::_S_construct NULL not valid
Aborted
Is this a normal behavior of the program??

How do I use generateILDM??

Regards
Jose

Someone can help me in solving such problem? And furthermore can someone maybe share with us an example of the generateILDMdict that is not included in the code release?

Thanks in advance
Stefano
Stefano Puggelli is offline   Reply With Quote

Old   March 12, 2015, 07:12
Default
  #54
Member
 
Join Date: Feb 2015
Posts: 30
Rep Power: 3
Stefano Puggelli is on a distinguished road
Ok first problem solved!
Anyway if someone would share his experience on the FGM model implemented by Hannes I would be interested in!
Stefano Puggelli is offline   Reply With Quote

Old   March 17, 2015, 10:07
Default
  #55
Member
 
Join Date: Feb 2015
Posts: 30
Rep Power: 3
Stefano Puggelli is on a distinguished road
Hi foamers,
I am trying to use the generateILDM solver released by Hannes Kroger with the rest of the code but I have some problems (maybe someone can help me).
I am running the tutorial test case for the flamelet generation for the Gri30 mechanism provided with the code but I noticed that when I try to create a flamelet with a p=100000Pa T=300K and phi=0.8 (for example but I think that the problem is the temperature level) the solver says me:

Specified enthalpy h=-205470 is below the absolute minimum enthalpy hzero=-8187.86

I understand from the ChemicalSystem.C file what is the enthalpy hzero ( I think it should be the enthalpy of products at the minimum temperature fixed at 290K), but I don't understand how I can solve this problem and why the actual enthalpy with the specified conditions is negative.

Let me know if someone has some ideas about these points... maybe we can discuss together.

Best,
Stefano
Stefano Puggelli is offline   Reply With Quote

Old   April 2, 2015, 07:45
Default
  #56
New Member
 
Jaehyeon Kim
Join Date: Sep 2014
Posts: 4
Rep Power: 3
Ryunjin is on a distinguished road
Hi Stefano,

I have same problem that you have.

(hzero is too high to start downward z loop.)

Did you solve the problem?

Best regards
Ryunjin is offline   Reply With Quote

Old   April 6, 2015, 01:40
Default I solved the problem!
  #57
New Member
 
Jaehyeon Kim
Join Date: Sep 2014
Posts: 4
Rep Power: 3
Ryunjin is on a distinguished road
Hi Stefano,

I solved the problem.

when conservedVaribleLoops variable is 'z' in 'generateIldmDict' file,

Change 'generateIldmDict' file from

defaultConservedVariables 3
(
p 100000
phi 1.0
T 300
);

to

defaultConservedVariables 3
(
p 100000
phi 1.0
Tunburnt 300
);

when conservedVaribleLoops variable is 'z',

mixture (unburnt gas) temperature have to be fixed.

Otherwise, mixture (unburnt gas) enthalpy is fixed at each flamelet.

but when mixture fraction is very small(e.g. z=0), that gas can't be same

as z=stoichiometric gas enthalpy. so error occurs.

I am so sorry. I'm not good at English. I hope you can understand what i'm

saying.
arvindpj likes this.
Ryunjin is offline   Reply With Quote

Old   April 6, 2015, 05:16
Default
  #58
Member
 
Join Date: Feb 2015
Posts: 30
Rep Power: 3
Stefano Puggelli is on a distinguished road
Hi Ryunjin,
I find for the moment the same solution that you reported even if in this way we generate a group of premixed flamelet (based on the mixture fraction value) with different values of enthalpy, while fixing the T temperature the flamelet are generated at the same enthalpy level. I think that it could be a problem when you are interested in generating flamelets for a iso-entalphy test case.
I am now trying to validate the FGM approach using the bunsen-Chen test case (one of the examples) with a direct comparison both with experimental and commercial codes results.
If the results are good I would like to try to implement such approach in actual version of OpenFOAM.
Stefano
Stefano Puggelli is offline   Reply With Quote

Old   May 14, 2015, 07:40
Default Z-CO2 FGM is not working well
  #59
New Member
 
Jaehyeon Kim
Join Date: Sep 2014
Posts: 4
Rep Power: 3
Ryunjin is on a distinguished road
Hi Stefano,

Now i'm solving BunsenChenF2 case. but results are not good solution.

My results are not similar to within result data. (Hannes' result data maybe.)

Moreover, temperature values are not instantanous changing each time.

(Temperature change at main stream part of near inlet only.)

The most bad thing is Z_mean value does not change anything each time. (just freezing.)

I mean results looks almost converged to solution(after calculating a few times , then temperature and mixture fractions are almost freezing.)

I have validated Pitz&Daily 3d case before calculating Bunsen case.

But that case (Pitz&Daily case) works well. So i think something related with mixture fraction have problem.

Please reply anything. thank you.
Ryunjin is offline   Reply With Quote

Old   May 19, 2015, 12:02
Default
  #60
Member
 
Join Date: Feb 2015
Posts: 30
Rep Power: 3
Stefano Puggelli is on a distinguished road
Hi Ryunjin,
I am still working on the bunsen test case since, instead of using the turbulence generator of the University of Rostoc, I prefer to use the mappedInlet boundry condition so some time has been required for the non reactive test case. In the next day I will try the reactive with the PDF solver.
I don't know which could be your problem since I don't know exactly your numerical setup but if you can post some other details we can discuss togheter.
I am also working on the PItz&Dayly case but with the ATF solver updated with the Charlette formulation for the efficency function.
Stefano
Stefano Puggelli is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Liquid rocket engine non-premixed combustion model Erik FLUENT 2 November 28, 2013 15:09
second and final question about flamelet combustion model AdidaKK CFX 0 October 17, 2009 07:25
Combustion Model for CAI/HCCI - AVL FIRE Patryk Main CFD Forum 0 July 24, 2009 16:43
How to model combustion with only methane stream? ligang zheng FLUENT 2 May 1, 2007 09:16
combustion model Hennie van der Westhuizen CD-adapco 7 February 27, 2002 03:10


All times are GMT -4. The time now is 07:26.