[sqek]

Hello

I am modelling large-amplitude thermoacoustic oscillations in FLUENT, using a transient simulation, and the partially premixed combustion model.
I am using the c-equation chemical equilibrium model, with compressibility effects and non-adiabatic treatment, with a c2h4-air mixture.

I would like to calculate the rate of chemical heat release (in W or W/m3) on a cell-by-cell basis (in such a way that I could plot contour maps etc of heat release), in order to see how the phase difference between heat release and pressure varies throughout the combustor.

Is there a way to do this with a field function or UDF please?

Thanks very much