CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Main CFD Forum (https://www.cfd-online.com/Forums/main/)
-   -   Combusting flow with incompressible code? (https://www.cfd-online.com/Forums/main/107212-combusting-flow-incompressible-code.html)

 michujo September 20, 2012 05:06

Combusting flow with incompressible code?

Hi everyone, I have a question concerning the computation of a reacting flow (combustion).

Do you think it is possible to compute a flame using the SIMPLE method?

What I think is that, assuming steady state (so no acoustics going on) and for low-Mach number flows, the gases are essentially incompressible (density differences only due to temperature changes, since the pressure is almost constant accross the flame), therefore I think that it is possible to compute it with an incompressible code (I'm particularly referring to the methodology proposed in Versteeg & Malalasekera or Patankar).

Cheers,
Michujo.

 yonchong September 20, 2012 11:06

Yes. Most of the combustion CFDs are done using SIMPLE and incompressible.

 leflix September 20, 2012 17:14

Quote:
 Originally Posted by michujo (Post 382721) Hi everyone, I have a question concerning the computation of a reacting flow (combustion). Do you think it is possible to compute a flame using the SIMPLE method?
the use of SIMPLE is not a characteristic of incompressible flows. You can use SIMPLE algorithm even for compressible flows.
The question is more can we use a incompressible solver for combustion problem? If your mach number is below 0.3 I guess we can...

 michujo September 20, 2012 17:20

Hi guys, thanks for your feedback.
You're right, I was referring to the use of an incompressible code for combustion, so you pretty much gave a satisfactory answer :). I still have some questions though. For instance: let's say that density depends on temperature, and I compute it as rho=function(T), would mass still be conserved?

My experience with combustion is using a non-commercial LES compressible solver so I was not familiar with the use of compressible codes for combustion. Do you guys know of any good reference that i could check in this regard?

Thanks again.

Cheers,
Michujo.

 yonchong September 20, 2012 17:28

I would recommend using ideal gas option for calculating density. Most of (if not all of) them should have it.

 michujo September 20, 2012 18:01

Thanks yonchong. Actually I ask so because I have in mind developing a 2D incompressible code of my own, so no commercial. Given that pressure is constant (except for the few kilopascals that would drive the flow), the ideal gas law is reduced to a function solely of temperature, is it not? (although I'll have to deal with the change of molecular weight of the gas for a multispecies flow...). Too many doubts!

Thanks.

Cheers,
Michujo.

 yonchong September 20, 2012 18:21

Pressure may not be changing but still need to calculate density.

 yonchong September 20, 2012 18:25

By the way, if you are trying to write your own CFD code have a look at these posts.

http://www.cfd-online.com/Forums/mai...-language.html

http://www.cfd-online.com/Forums/mai...tml#post376764

 julien.decharentenay September 20, 2012 19:21

Michujo,

Mass conservation in SIMPLE (or any pressure based algorithm) is achieved in the pressure correction step. As long as the pressure correction step includes the density variation associated with variation in temperature, the mass is "conserved" - I used inverted comma here as mass conservation depends on how good the pressure solver is and the tolerance used. As a good practice I would recommend that you implement a mass conservation check.

I would recommend doing a bit of reading on the derivation of the low Mach number approximation for pressure based equation (google dilatable incompressible low Mach number equation may give you to the right answer). Otherwise I can try to dig out the references from my PhD (unsteady simulation of reactive flow at low Mach numbers).

Julien

 michujo September 21, 2012 02:41

Dear all, thanks for your contributions. I will try and check out those references.

Cheers,
Michujo.

 All times are GMT -4. The time now is 17:15.