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

Best regards,
Likun

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

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!

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

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!

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

Yes, I am still interested in FGM, hope you can send me a copy. Thank you
:)

yangdongchao500@126.com

Already sent! :)

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?

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

CombustionAddon's example question
 

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?

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

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:


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

Ok first problem solved!
Anyway if someone would share his experience on the FGM model implemented by Hannes I would be interested in!

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

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

I solved the problem!
 

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.

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

Z-CO2 FGM is not working well
 

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.

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