eddy dissipation concept + infinitely fast chemistry
dear all!
i implemented the eddy dissipation concept combustion model with infinitely fast chemistry according to the paper of Magnussen: "The Eddy Dissipation Concept A Bridge Between Science And Technology", Lisbon, 2005, based on fireFoam, using a keps turbulence model. the code runs and the flow field looks plausible at first sight (in 3D). but when integrating the reaction rate (or the heat release due to burnt fuel) over the whole computational domain, I can see that the reaction rate is overestimated, i.e. more fuel burns as is injected. I tried to increase spatial resolution as well as to change the value for the turbulent dissipation rate at fuel inlet but always observe the same behavior. Interestingly the reaction rate matches quiet well with the injected amount of fuel for the 2D tutorial case applying keps turbulence model. Anybody experienced similar problems?? I would greatly appreciate your comments!! Regards, Aram 
dear all!
i figured out that the problem is not the combustion model. additional to the convective input of fuel (prescribed by mass flow), fuel enters the domain by diffusion due to the sharp gradient of fuel mass fraction fu at inlet. first i tried to add a long fuel inlet pipe but as a result of low inlet velocities the combustion process "travels" back to the inlet causing the same problem as before. Hence, i decided to set the diffusion coefficient alphaEff in the transport equation of fu (fuhsEqn.H) to zero at inlet, what solved the problem: Code:
alphaEffZeroDiff = turbulence>alphaEff(); Code:
const label patchI = patch().index(); Thx in advance, Aram 
Dear mabinty,
I am working with the solver since a few month and I discovered the problem very lately. So I will see to confirm it as soon as possible. In fact, I still haven't found a possibility to compute the diffusive flux. It looks like this is still unsolved: http://lists.cfdonline.com/pipermai...ne/004327.html also: http://www.extendproject.de/userfo...libraryupdate Cheers, Michael 
Sorry for the double post due to a missing reload, questions already raised above (Apart from deactivating completely the diffusion, is there any way to compute the diffusive mass flow through a patch?).

1 Attachment(s)
hey Michael,
my problem was not to compute diffusive fluxes. for a variable phi and a diffusion coefficient alphaEff i d try it like that Code:
volScalarField dDiff = cheers, aram 
All times are GMT 4. The time now is 21:15. 