|
[Sponsors] |
July 16, 2011, 10:33 |
chemical time scale
|
#1 |
Member
José Rodrigues
Join Date: Jun 2010
Location: IN+/IST Lisbon
Posts: 53
Rep Power: 15 |
Hi,
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? Jose |
|
August 16, 2011, 03:26 |
|
#2 |
Member
Tibo
Join Date: Jun 2011
Posts: 68
Rep Power: 14 |
Hi José,
If you´re simply looking for the formula, you´ve got a first answer on slide 16 of the following presentation: http://www.openfoamworkshop.org/6th_...ist_slides.pdf 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. Tibo |
|
November 19, 2013, 12:03 |
Chemical time scale in PaSR model
|
#3 |
Member
Manuel
Join Date: Jan 2013
Location: Bremen
Posts: 33
Rep Power: 13 |
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: Code:
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]; } } 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. Manuel Last edited by wyldckat; November 19, 2013 at 17:39. Reason: Added [CODE][/CODE] |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Transient simulation not converging | skabilan | OpenFOAM Running, Solving & CFD | 14 | December 16, 2019 23:12 |
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 |
AMG versus ICCG | msrinath80 | OpenFOAM Running, Solving & CFD | 2 | November 7, 2006 15:15 |