Solving Complex Combustion Challenges
In this tech tip, we will discuss the omnipresence of combustion. Combustion is everywhere: from burners around the world to coal combustion, gas turbines in aircraft engines, and furnaces. Today, we will focus on two specific examples:
- Larger dissimulation of natural gas combustion and CO pollutant emission prediction, based on an example from the DLR in Germany, conducted by my colleague Nassim Ansari.
- A lifted methane-air jet flame in a vitiated cold flow, using data and setup from the University of Berkeley, with simulations conducted by my colleague Pravin Nakhut.
Example 1: Natural Gas Combustion and CO Emission Prediction
This example involves the Pre-Synstra burner from the DLR in Germany. The process involves mixing air and methane fuel, swirling them, and observing the combustion in the chamber to predict CO emissions.
- Geometry Overview: The combustion chamber's geometry includes an air inlet with a swirler for mixing air and fuel.
- Mesh Details: High-quality polyhedral mesh and boundary layers are used to ensure accurate flow resolution.
Before diving into the simulation, a brief theoretical overview:
- We track two scalars of interest: the fuel-air mixing ratio and the flame location.
- Using partially premixed combustion, we employ the C equation to track the flame front and generate a flamelet-generated manifold for flame speed and temperature data.
Simulation setup involves:
- Importing a Chemkin mechanism from Reaction Design.
- Defining boundaries with methane as fuel and air as the oxidizer.
- Generating flamelets to obtain post-flame temperature and flame speed information.
Results include:
- Isosurface of velocity colored by flame temperature, showing a cup shape due to swirling flow.
- Isosurface of the mixture fraction, illustrating the air-to-fuel ratio and recirculation region.
- Animation of the temperature field on the mid-plane of the combustion chamber, capturing flame structure.
- CO mass fraction isosurface, highlighting carbon monoxide as a pollutant.
Example 2: Lifted Methane-Air Jet Flame in a Vitiated Cold Flow
This example features a simpler geometry with a steady simulation. Fuel is injected into a central tube, surrounded by a cold air flow, simulating combustion.
- Key Feature: Simulating a lifted flame, where the flame appears to sit in the middle of the air.
Validation is crucial, and the results show:
- Temperature in the centerline of ANSYS results versus experiments, showing excellent agreement.
- Good prediction of flame lift distance and CO emission, dependent on model and mesh resolution.
Overall, using the right model and mesh resolution provides accurate CO predictions, validated by experimental data.
Hi, for this tech tip we will be talking about combustion and as you can see combustion is everywhere. Think about all the burners that exist around the world. Think about coal combustion, think about gas turbines in aircraft engines for example, think about furnaces etc.
But today I want to focus on two examples. The first one, larger dispersion of natural gas combustion and CO pollutant and emission prediction in an example from the DLR in Germany and my colleague Nassim Ansari did this example.
Another example that uses data and setup from the University of Berkeley will be about a lifted methane air jet flame in a vitiated coal flow and thank you very much for my colleague Pravin Nakhat for doing this simulation.
Now on the first example, larger dispersion of natural gas combustion and CO pollutant and emission prediction. Without further ado, let's go to the first example from the DLR in Germany.
It's called the Pre-Synstra burner and what you can see here is we're going to take air and fuel, the fuel being methane. We're actually going to mix them, make them swirl as well and we'll observe the combustion in the combustion chamber and look at CO emission prediction.
Here is an overview of the geometry of the combustion chamber and what we're seeing now is actually the inlet of air and as you can see we have a swirler and the air will actually swirl and be mixed directly with the fuel.
Here's a little bit more detail of the swirler and the mesh on the swirler and the injection code. What you can also see is the high-quality polyhedral mesh in this plane that I cut in the middle of the combustion chamber.
And if we zoom in we also see the high-quality boundary layer that we put near to the wall to make sure we resolve the boundary layer of the flow correctly. And as usual, high-quality surface mesh as well. Now just a little bit of theory before we start looking at the simulation.
What we'll do here is obviously not track the hundreds of species involved in the combustion, that would be too expensive, but we'll track two scales of interest. We will track the fuel-air mixing ratio, it will give us information like flame speed, post-flame temperature, etc.
And of course, as well, the flame location that we will track and actually the flame surface as you can see on the schematic at the bottom of the picture. Now let's see how that's done. I actually go in the species panel, I will open it. Let me put it back in the center for you to see everything.
We're going to choose partially premixed combustion. We're going to use the C equation, that's what's going to track our flame front. We'll use the flame lead generated manifold to generate all this information we talked about like flame speed, flame temperature, etc.
We set up our operating pressure and what we do is actually we simply import a Kempkin mechanism from reaction design. You see the mechanism is imported and I also import the thermodynamic properties. Then I go to boundaries. What are my boundaries?
Well on fuel I have full methane and on the oxidizer side I have air. So you see that I entered the composition of air. I give the fuel temperature and oxidizer temperature. And then I will generate those flame leads.
Those flame leads that will give me the information I'm looking for on post-flame temperature, flame speed, etc. Now let's generate those and of course what I can do is check the flame that we generated for example temperature as a function of the mixture fraction.
And what you can see here is the temperature to which the mixture fraction makes. As you can see, temperature is dependent on phase couples of course.
The temperature emits the temperature to suffice concentration and the temperature A-.» So you take the temperature and see how this could generate densities. Let's generate that. I badly unload on the test and there are four different functions, the first one is the evening temperature.
Here is the temperature peak when the mixture fraction is actually a stoichiometric mix, which is completely expected. Next thing I have to do is actually take those flamelets, but also include the influence of the turbulent flow. And that's why I generate a PDF table.
And the PDF table will be what contains the key information we're going to use in the simulation for chemistry. Now, I will not go into the detail of running a large eddy simulation computation because you can find that in the turbulence tech tip.
What I will do, however, is show you some of the results. Here you see an isosurface of velocity colored by the flame temperature. And you can see the cup shape of this isosurface. Why? Because the flow is swirling. And actually that creates...
a recirculation region in the middle of the combustion chamber. And that's what gives you this shape. Here is another view. Isosurface of the mixture fraction. Or that's basically an isosurface of the air to fuel ratio. And here you see the same characteristics.
You see the shape that is a little bit like a cup. And that's because of the recirculation region present in the middle of the combustion chamber. And other way to look at it, you could... take for example an animation of the temperature field on the mid-plane of the combustion chamber.
As you can see here. Typical of LES simulation. Those nice simulations with a lot of structure and flame structure that you can capture. And as I told you, we were going to look at prediction of pollutant emission. So here I give you an isosurface of the CO mass fraction.
CO being for carbon monoxide. Which obviously is a pollutant. Now let's go to the second example. Lifted methane air jet flame in a vitiated cold flow. Now the geometry of this example is more simple. And it will be a steady simulation. We are injecting fuel in a tube at the center.
And we have a cold flow of air all around the fuel. And we simulate the combustion in this system. And what is interesting in this simulation. Is that we can actually... simulate what we call a lifted flame. Really when the flame looks like sitting in the middle of the air.
And that's actually something that is not easy to do. Now I just said it is not easy to simulate. So I know what you want to see. You want to see validation. And here is a great validation. For actually looking at the temperature in the centerline of ENSYS results versus experiments.
And as you can see they agree very well. And there is a really good prediction of the flame lift distance. CO emission prediction. Now there are a lot of models available. And results can be dependent upon the grid resolution. But what's important to note here.
Is when you use the right model and the right mesh resolution. You can have a pretty good estimate of the CO prediction. You see the experimental data with the dots. And you see all the results. And you can see your best result in blue. And don't be fooled. Those results are actually extremely good.
And that's actually a very good validation of the simulation.

