CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

OnlineLibrary in chemistryModelC Tommaso development

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 7, 2007, 09:19
Default Hi.. this message is particul
  #1
zurdo
Guest
 
Posts: n/a
Hi..
this message is particularly for Tommasso, who I belive implemented the chemistryModel.C

Cold you give me a brief explanation about this onlinelibrary as I am having errors using the combustion model with chemkin files.

There is a complete set of info as: tolerance,logT,cleanAll, etc in the online distionary, which i dont know what they mean?

Some tips are welcome.

Thanks.

Sergio
  Reply With Quote

Old   August 7, 2007, 14:58
Default Dear Sergio, onlineLibrary
  #2
Member
 
Tommaso Lucchini
Join Date: Mar 2009
Posts: 83
Rep Power: 9
lucchini is on a distinguished road
Dear Sergio,

onlineLibrary is used to tabulate the complex chemistry to reduce the computational time. It is a constant approximation version of the ISAT (in-situ adaptive tabulation) algorithm developed by Pope.

If you want to use the tabulated chemistry, the onlineProperties dictionary should be written as follows:

onlineProperties
{
tauStar off;

online on;

tolerance 1.0e-5;

logT off;

maxElements 100000;

cleanAll off;

scaleFactor
{
otherSpecies 1.0;
Temperature 1700;
Pressure 1.0e+10;
}

scaleFactorSolution
{
otherSpecies 1.0;
Temperature 1700;
Pressure 1.0e+10;
}

}

if you want to switch it off,

onlineProperties
{
tauStar off;

online off;

tolerance 1.0e-5;

logT off;

maxElements 100000;

cleanAll off;

scaleFactor
{
otherSpecies 1.0;
Temperature 1700;
Pressure 1.0e+10;
}

scaleFactorSolution
{
otherSpecies 1.0;
Temperature 1700;
Pressure 1.0e+10;
}

}

In the ISAT algorithm chemical source terms are tabulated into a binary tree and are used to approximate the chemical source terms of the nearest composition according to a specified tolerance.

Explanation of the terms:

tauStar: switch it off, it will be used in the future;
logT: if this option is switched on, the logarithm of the temperature will be tabulated instead of the temperature
maxElements: maximum number of elements to be tabulated (up to 5e6 is fine, then you can seriously run into memory troubles)
cleanAll: if switched on, it will clean the library after any time step. If you want to exploit the full capabilities of the method, switch it off
scaleFactor: tabulated data consist of chemical compositions, temperature and pressure. You can scale chemical compositions, temperature and pressure to get more points found
solutionScaleFactor: the same but used to extend the validity range of a tabulated composition.

For further details, please refer to the Pope's paper (you can find them on the web also) and to the code implementation. In particular
- chemistrySolve.C
- chemistryOnlineLibrary.C
should be fine to start looking at them.

If you need more help or you like to setup a possible collaboration, please contact me.

That's all for the moment, have fun with OpenFOAM and combustion.

Tommaso
lucchini is offline   Reply With Quote

Old   August 8, 2007, 06:30
Default Dear Tommaso, Thanks for y
  #3
zurdo
Guest
 
Posts: n/a
Dear Tommaso,

Thanks for your replay.
I am or I was working extensively in fire and combustion. Actually, I did not use much detailed chemistry but mostly Flame Surface Density approaches and Thickened Flame in LES.
I was in the seminar in Zagreb and I congratulate you for you work in internal comb engines.
This is my first time using PsR model for combustion and this shortcut called ISAT. I am printing out some papers and see if they are useful for me.

These are my lines of work:

1) I am particularly interested in autoignition. I believe that it is important for IC as well?.

I am using a 3 step reaction CH4-O2 mechanism. Originally developed for P up to 100 bar.

2) DDD(Detonation to Deflagration). We are not sure about the combustion model to choose. In Open FOAM we could find: the progress variable approach (FSD) using wrinkling factor and burning turbulent velocities and PsR (using ISAT). I dont see any other combustion model.
In this DDD project, I am doing some mesh refinement by scalar quantity. I made some progress but it takes time and I don't have it. So, far I managed to do the dynamic refinement but there is an error when when in the unrefinement. The idea is to follow the flame front at the shock and combustion wave. This is samething could be interesting for you as to refine where the flame front it is.

3) Fires, I need to implement a new combustion model for fire in Open FOAM and new radiation model.

Thanks for your time. and if it is possible to find some common ground I am happy to any collaboration. If you want what I have done so far for the refinement you can have a look. Just let me know.

Thank

Sergio Ferraris.
  Reply With Quote

Old   December 19, 2007, 12:19
Default Hello Tommaso, a few questi
  #4
Senior Member
 
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 182
Rep Power: 9
markusrehm is on a distinguished road
Hello Tommaso,

a few questions regarding the ISAT-tabulation:

I put your suggestions in the chemistryProperties file in the aachenBomb case. It doesn't seem to work. Do I have to switch it on at another place? Is there any output about the tabulation? I use the 1.4.1-dev-version.

Does it run with other solvers like reactingFoam as well?

Is the table printed out to a file?

Regards Markus.
markusrehm is offline   Reply With Quote

Old   December 19, 2007, 19:40
Default hi, where did you find ISAT
  #5
Senior Member
 
Stephan Gerber
Join Date: Mar 2009
Location: Germany
Posts: 118
Rep Power: 9
stephan is on a distinguished road
hi,

where did you find ISAT in the sources - and in which version is it implemented.
i was looking through the sources but couldnt find anything.
thanx
stephan
stephan is offline   Reply With Quote

Old   December 20, 2007, 04:21
Default Stefan, did you get the fil
  #6
Member
 
Tommaso Lucchini
Join Date: Mar 2009
Posts: 83
Rep Power: 9
lucchini is on a distinguished road
Stefan,

did you get the files I sent you by e-mail?
lucchini is offline   Reply With Quote

Old   December 20, 2007, 04:37
Default hi, actually not - would be
  #7
Senior Member
 
Stephan Gerber
Join Date: Mar 2009
Location: Germany
Posts: 118
Rep Power: 9
stephan is on a distinguished road
hi,

actually not - would be nice if you can try it again..
thanx
stephan
p.s. please dont forget the .net instead of .de at the end.
stephan is offline   Reply With Quote

Old   December 20, 2007, 06:11
Default Marcus, please send me an e
  #8
Member
 
Tommaso Lucchini
Join Date: Mar 2009
Posts: 83
Rep Power: 9
lucchini is on a distinguished road
Marcus,

please send me an e-mail with your correct address since nothing of them works.

Thanks a lot,

regards,

Tommaso
lucchini is offline   Reply With Quote

Old   December 20, 2007, 08:38
Default Hello Tommaso, yes I was in
  #9
Senior Member
 
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 182
Rep Power: 9
markusrehm is on a distinguished road
Hello Tommaso,

yes I was in Zagreb. And I got your files. Thank you very much. I thought I saw the files in the dev-Version on the svn-Server. I am already on holiday so i can`t give you more information if it works or not. But next year... ;)

Regards and a blessed Christmas time to all. Markus
markusrehm is offline   Reply With Quote

Old   April 13, 2009, 04:49
Default
  #10
New Member
 
Arun
Join Date: Mar 2009
Posts: 25
Rep Power: 9
arun is on a distinguished road
Hello Tommaso,

I am currently simulating detailed chemical kinetics combustion of diesel spray using dieselEngineFoam. I am using OpenFOAM-1.5.x.

Where can I find the ISAT in the source file?
And how do I "implement/switch on" ISAT in my case setting?
Is ISAT is implemented in dieselEngineFOAM?

Thanks in advance

Arun
arun is offline   Reply With Quote

Old   July 27, 2009, 20:59
Default
  #11
Member
 
edison
Join Date: May 2009
Location: Australia
Posts: 35
Rep Power: 9
Edison_Ge is on a distinguished road
Hello Tommaso and everyone,
I'm doing PDF combustion simulation and now I'm using a DVOD for my chemistry. I find your implementation of the ISAT very interesting and eager to try it. But my concern is like:
1) I'm working on a 4 processor personal workstation. Does this ISAT implementation support parallezation.
2) I found some paper discussing ISAT with reduced chemistry, but I'm working with GRI-3.0 for methane/air. How's the codes' performance?

Hope some familiar with this chemistryonlinelibrary could help me clear these doubts.
Edison_Ge is offline   Reply With Quote

Old   March 17, 2011, 14:17
Default
  #12
Member
 
Ali Shamooni
Join Date: Oct 2010
Posts: 41
Rep Power: 7
Alish1984 is on a distinguished road
Can anyone help me about using chemistryonline class in OF 1.5-dev?
Alish1984 is offline   Reply With Quote

Old   December 26, 2011, 16:52
Smile ISAT ALGORITHM implementation in diselEnginefoam
  #13
New Member
 
Haris
Join Date: Jul 2010
Posts: 8
Rep Power: 8
riseguess is on a distinguished road
Quote:
Originally Posted by arun View Post
Hello Tommaso,

I am currently simulating detailed chemical kinetics combustion of diesel spray using dieselEngineFoam. I am using OpenFOAM-1.5.x.

Where can I find the ISAT in the source file?
And how do I "implement/switch on" ISAT in my case setting?
Is ISAT is implemented in dieselEngineFOAM?

Thanks in advance

Arun
Hi Arun
did you to find answer to your to query. Can you please guide me in this regard as I want to use ISAT algorithm in dieselEngineFoam.

Shahid
riseguess is offline   Reply With Quote

Old   December 26, 2011, 17:09
Default ISAT ALGORITHM IN DIESELeNGIENFOAM
  #14
New Member
 
Haris
Join Date: Jul 2010
Posts: 8
Rep Power: 8
riseguess is on a distinguished road
Quote:
Originally Posted by markusrehm View Post
Hello Tommaso,

a few questions regarding the ISAT-tabulation:

I put your suggestions in the chemistryProperties file in the aachenBomb case. It doesn't seem to work. Do I have to switch it on at another place? Is there any output about the tabulation? I use the 1.4.1-dev-version.

Does it run with other solvers like reactingFoam as well?

Is the table printed out to a file?

Regards Markus.
Hi Markus
did you manage to run the ISAT algorithm with dieselEngineFoam successfully. Can you guide me in this reagrd
Shahid
riseguess is offline   Reply With Quote

Old   December 26, 2011, 17:15
Default
  #15
New Member
 
Haris
Join Date: Jul 2010
Posts: 8
Rep Power: 8
riseguess is on a distinguished road
Quote:
Originally Posted by lucchini View Post
Dear Sergio,

onlineLibrary is used to tabulate the complex chemistry to reduce the computational time. It is a constant approximation version of the ISAT (in-situ adaptive tabulation) algorithm developed by Pope.

If you want to use the tabulated chemistry, the onlineProperties dictionary should be written as follows:

onlineProperties
{
tauStar off;

online on;

tolerance 1.0e-5;

logT off;

maxElements 100000;

cleanAll off;

scaleFactor
{
otherSpecies 1.0;
Temperature 1700;
Pressure 1.0e+10;
}

scaleFactorSolution
{
otherSpecies 1.0;
Temperature 1700;
Pressure 1.0e+10;
}

}

if you want to switch it off,

onlineProperties
{
tauStar off;

online off;

tolerance 1.0e-5;

logT off;

maxElements 100000;

cleanAll off;

scaleFactor
{
otherSpecies 1.0;
Temperature 1700;
Pressure 1.0e+10;
}

scaleFactorSolution
{
otherSpecies 1.0;
Temperature 1700;
Pressure 1.0e+10;
}

}

In the ISAT algorithm chemical source terms are tabulated into a binary tree and are used to approximate the chemical source terms of the nearest composition according to a specified tolerance.

Explanation of the terms:

tauStar: switch it off, it will be used in the future;
logT: if this option is switched on, the logarithm of the temperature will be tabulated instead of the temperature
maxElements: maximum number of elements to be tabulated (up to 5e6 is fine, then you can seriously run into memory troubles)
cleanAll: if switched on, it will clean the library after any time step. If you want to exploit the full capabilities of the method, switch it off
scaleFactor: tabulated data consist of chemical compositions, temperature and pressure. You can scale chemical compositions, temperature and pressure to get more points found
solutionScaleFactor: the same but used to extend the validity range of a tabulated composition.

For further details, please refer to the Pope's paper (you can find them on the web also) and to the code implementation. In particular
- chemistrySolve.C
- chemistryOnlineLibrary.C
should be fine to start looking at them.

If you need more help or you like to setup a possible collaboration, please contact me.

That's all for the moment, have fun with OpenFOAM and combustion.

Tommaso
Hi Tommaso
Can you please guide me how can i use the ISAT algorithm with dieselEngineFoam successfully. I am new to openFoam world and would appreciate your guidance
Shahid
riseguess 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
cfd code development jojo Main CFD Forum 9 April 22, 2007 11:00
Problems with BL development Dan Main CFD Forum 0 October 2, 2006 11:24
CFD software development Jonathan Main CFD Forum 2 August 8, 2006 03:12
Code development terry towelling Main CFD Forum 0 May 18, 2006 04:05
Openfoam development? Hao Meng Main CFD Forum 1 January 9, 2006 13:57


All times are GMT -4. The time now is 08:58.