Hi, I noticed that in OF 1.
I noticed that in OF 1.5 the implementation of chemistryModel::tc() changed:
OF 1.4 says
if(RR_[i][celli] < -SMALL)
and OF 1.5:
for(label celli=0; celli<nCells; celli++)
scalar rhoi = rho_[celli];
scalar Ti = thermo_.T()[celli];
scalar pi = thermo_.p()[celli];
scalar cSum = 0.0;
for(label i=0; i<Ns_; i++)
scalar Yi = Y_[i][celli];
c[i] = rhoi*Yi/specieThermo_[i].W();
cSum += c[i];
const reaction& R = reactions_[i];
R, c, Ti, pi, pf, cf, lRef, pr, cr, rRef
scalar sr = R.rhs()[s].stoichCoeff;
t[celli] += sr*pf*cf;
t[celli] = Nr*cSum/t[celli];
Could anyone explain why this change has been made and how the new implementation works?
Thanks in advance. Markus
This was some time ago so its
This was some time ago so its not entirely clear to me...
The first definition is simply an estimation of the shortest time for a species to be driven to zero.
The second is a bit more coupled to the PaSR combustion model in diesel(Engine)Foam.
If you write the reaction rates in terms of a positive, P, and negative, N, term and linearize the negative term around the concentration c.
omega = P - N = P - c/tau_c
you get the definition of tau_c.
This yields a chemical timescale for each species and
If I remember correctly the calculated tau_c is an averaged value over all species and reactions.
|All times are GMT -4. The time now is 23:16.|