[Weiqiang Liu]

Hi all,

I am doing a methane combustion simulation with fluent. The mechanism is imported in chemkin format. I just followed all settings with literature.

My problem is: The general trends are very similar with literature results. However, the species profiles are not very reasonable. I suspect it's because of my case is not converged yet.

I check the flux report of mass and energy and they are all balanced. however, the residuals are not all below 1e-6 and they are oscillating with small magnitude.

Can anybody who have combustion simulation experience give me some suggestions on how to get full converged results.

Best regards

Weiqiang

[vinerm]

1e-6 is already too small. Even if residuals are below 1e-4, it's more than good enough until the values of species fractions are themselves of that order.

[Weiqiang Liu]

Originally Posted by vinerm

1e-6 is already too small. Even if residuals are below 1e-4, it's more than good enough until the values of species fractions are themselves of that order.

Hi Vinerm,

In my case, the mole fraction of combustion products starts from 2% from inlet and increases to 16% at outlet.

However, in literature the combustion products increases from 0 to 16%, which is more reasonable I think.

also mole fraction of fuel at inlet is smaller than the value I set. they should be the same I think.

Best regards

Weiqiang

[vinerm]

When the post processing is done, even at the inlet boundary, values are plotted based on the values at the centers of the adjacent cells. Therefore, what you are looking at is actually the solution at the center and those values extrapolated to the inlet.

[Weiqiang Liu]

Originally Posted by vinerm

When the post processing is done, even at the inlet boundary, values are plotted based on the values at the centers of the adjacent cells. Therefore, what you are looking at is actually the solution at the center and those values extrapolated to the inlet.

Hi Vinerm,

Do you mean the profile or even contour are node values even if I selected node values in post processing? Because I selected node value when I export data files.

If all values are nodes values instead of boundary value, then how can the author draw species profile starting from 0?

Best regards

Weiqiang

[vinerm]

Node values are not even boundary values, those are node values that could be lying on the boundary. If you specify 0 mass (or mole) fraction for a specie at the inlet, it would be strange to have a rather high mass fraction within such as small distance as half-cell thickness. Reactions take time and flow can travel quite some distance before products are produced. But it all depends on the reaction model being used, flow speed, and flow structure; a highly diffusion controlled system can have high mole fraction of a product close to the inlet.

[Weiqiang Liu]

Originally Posted by vinerm

Node values are not even boundary values, those are node values that could be lying on the boundary. If you specify 0 mass (or mole) fraction for a specie at the inlet, it would be strange to have a rather high mass fraction within such as small distance as half-cell thickness. Reactions take time and flow can travel quite some distance before products are produced. But it all depends on the reaction model being used, flow speed, and flow structure; a highly diffusion controlled system can have high mole fraction of a product close to the inlet.

Hi Vinerm,

In my simulation, chemical reactions are very fast. All fuels are consumed in a very short distance. Can it be a highly diffusion controlled system?

Best regards

Weiqiang

[vinerm]

To look at whether a system is diffusion controlled or not, you need to compare diffusion and convection time-scales. Diffusion time scale is given by diffusion coefficients, such as, thermal diffusivity or mass diffusivity and convection time scale is given by ratio of domain length scale and flow velocity. If diffusion time-scale is much larger than convection, which it usually is, then it is not diffusion controlled. However, if diffusion is much faster, i.e., its time-scale is smaller, then system becomes diffusion controlled. This is possible in laminar or highly recirculating flows.

[Weiqiang Liu]

Originally Posted by vinerm

To look at whether a system is diffusion controlled or not, you need to compare diffusion and convection time-scales. Diffusion time scale is given by diffusion coefficients, such as, thermal diffusivity or mass diffusivity and convection time scale is given by ratio of domain length scale and flow velocity. If diffusion time-scale is much larger than convection, which it usually is, then it is not diffusion controlled. However, if diffusion is much faster, i.e., its time-scale is smaller, then system becomes diffusion controlled. This is possible in laminar or highly recirculating flows.

Hi, Vinerm,

Thanks very much for your patience. It enlightened me a lot.

Best regards

Weiqiang