Multiple species and UDS source terms

Ale

Hi to all,

I wanted to check the right way to write a source term in fluent, so I solved a simple problem of laminar non-premixed combustion, using the Multiple Species method and ticking off only the option Volumetric Reaction.

In the same simulation, I attempted to calculate the profiles of the species by UDF.

The source term I have written for methane is the following:

DEFINE_SOURCE(p_source_ch4, c, t, dS, eqn) {

double source_ch4;

source_ch4=-16/C_R(c,t)*(2.119e11*exp(-2.027e8/8314/C_T(c,t))*pow((C_UDSI(c,t,0)*C_R(c,t)/16),0.2)*pow((C_UDSI(c,t,1)*C_R(c,t)/32),1.3));

return source_ch4; }

C_UDSI(c,t,0) is methane mass fraction; C_UDSI(c,t,1) is oxygen mass fraction. Other source terms for other species are similar. The kinetic expression is the one used by Fluent for the methane-air mixture.

The problem is that contours of UDS are different from the corresponding mass fraction contours calculated by Fluent.

I used for diffusivity the constant-dilute-appx (2.88e-5) and I set to the same value for UDS Diffusivity.

Has anyone idea of the reason of this discrepancy? Where's my mistake?

Thank you very much!

Ale