CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

chemical time scale

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

Like Tree1Likes
  • 1 Post By jose_rodrig

LinkBack Thread Tools Display Modes
Old   July 16, 2011, 10:33
Default chemical time scale
José Rodrigues
Join Date: Jun 2010
Location: IN+/IST Lisbon
Posts: 53
Rep Power: 9
jose_rodrig is on a distinguished road

Does anybody knows how OpenFOAM calculates the chemical time scale, tc ?

I am deciphering the source code right now, but it would be nice to have somebody else's work to compare with.

Any help?

cdunn6754 likes this.
jose_rodrig is offline   Reply With Quote

Old   August 16, 2011, 03:26
Join Date: Jun 2011
Posts: 68
Rep Power: 8
megacrout is on a distinguished road
Hi José,

If you´re simply looking for the formula, you´ve got a first answer on slide 16 of the following presentation:

If you´re trying to determine how and where exactly this parameter is computed in OpenFOAM format, I´m afraid I can´t help. I actually looked for it - as I am running simulations with reactingFoam and find the question interesting - but did not quite get it. I guess it would help to have more than just my basic knowledge of C++!

I´d be glad to read your findings here in case you got more than that in the meantime.

megacrout is offline   Reply With Quote

Old   November 19, 2013, 13:03
Default Chemical time scale in PaSR model
New Member
Join Date: Jan 2013
Location: Darmstadt
Posts: 25
Rep Power: 7
Manuel CFD is on a distinguished road
Dear all,

I have a doubt concerning chemical time scale computation in ChemistryModel class. Mainly I'm referring to the tc() function reported in the following:

 if (this->chemistry_)
        forAll(rho, celli)
            scalar rhoi = rho[celli];
            scalar Ti = this->thermo().T()[celli];
            scalar pi = this->thermo().p()[celli];
            scalarField c(nSpecie_);
            scalar cSum = 0.0;

            for (label i=0; i<nSpecie_; i++)
                scalar Yi = Y_[i][celli];
                c[i] = rhoi*Yi/specieThermo_[i].W();
                cSum += c[i];

            forAll(reactions_, i)
                const Reaction<ThermoType>& R = reactions_[i];

                omega(R, c, Ti, pi, pf, cf, lRef, pr, cr, rRef);

                forAll(R.rhs(), s)
                    scalar sr = R.rhs()[s].stoichCoeff;
                    tc[celli] += sr*pf*cf;
            tc[celli] = nReaction*cSum/tc[celli];
Could anyone explain me in details how this chemical time is computed? Which is the physical expression taken as reference? Are there any literature works concerning the evalutation of this term?
I've read a lot of works about PaSR model where they talk about the general overview of this combustion model but in each of them they talk about the need of solving detailed reaction system in order to obtain this time. But I've never found how this has to be done.

Thank you.


Last edited by wyldckat; November 19, 2013 at 18:39. Reason: Added [CODE][/CODE]
Manuel CFD is offline   Reply With Quote


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
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM Running, Solving & CFD 8 March 28, 2016 11:08
cfx time scale control Kordou CFX 2 September 22, 2008 11:07
time scale in brownian force Mehdi Asghari FLUENT 0 May 17, 2008 01:02
Transient simulation not converging skabilan OpenFOAM Running, Solving & CFD 12 September 17, 2007 17:48
AMG versus ICCG msrinath80 OpenFOAM Running, Solving & CFD 2 November 7, 2006 16:15

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