|
[Sponsors] |
March 26, 2024, 18:51 |
Laminar Flow Combustion
|
#1 |
New Member
Asif Mahmud Rayhan
Join Date: Mar 2024
Posts: 2
Rep Power: 0 |
I am trying to re-do a paper work where the following set ups are used for getting the methane air combustion.
Axisymmetric. 1. Viscous Model: Laminar 2. Energy equation: on 3. Radiation: DO 4. Species Transport: Volumetric, Stiff Chemistry, Full multicomponent diffusion, GRI Mech 3.0 imported for transport and thermal data 5. velocity inlet: Fuel - 3m/s, 1 mole fraction CH4 at species air - 2m/s, 0.23 Oxygen & 0.77 Nitrogen mole fraction at species. Both 298 k temperature 6. Pressure Outlet 7. 1000 iteration and residual chosen 1*10^-6 OUTPUT: The reaction happened just at the inlet part and some temperature changes in the inlet and outlet part of the 2D flow region. Again velocity is found on those small part near to the inlet and outlet but the rest of the chamber's fluid's velocity remain zero. Why the flow and reactions are not happened in the middle part of the fluid region? I will appreciate all of your contribution. If you want I can share the paper with you. |
|
March 27, 2024, 17:38 |
|
#2 |
Senior Member
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 676
Rep Power: 21 |
Which solver? And which paper? Steady-state?
Usually, you start such a simulation first without any chemical reactions, only flow and transport of species. This is fast and you can then check that at least this part of the simulation works as expected. Then, you ignite a specific part of the domain by setting the temperature to a high value, like 1500K. You have to carefully choose this region to make sure the flame does not blow off or flash back. (Gas) Radiation and multicomponent diffusion is usually not important in methane-air flames, so you can switch it off in this phase because they are also time-consuming. Once it has converged a bit, you can activate these models. Also, if this is Fluent, then 1000 iterations and R=1e-6 is not enough. You probably need 100,000 iterations and R=1e-10. It might be more than that when you include radiation. Other solvers might be worse. |
|
April 1, 2024, 14:14 |
|
#3 |
New Member
Asif Mahmud Rayhan
Join Date: Mar 2024
Posts: 2
Rep Power: 0 |
I have tried according to the way you were asked. First I tried with air. Then I use CH4-air, and then chemkin GRI 3.0. Even when I am applying the volumetric reactions and Full multicomponent diffustion. It all works and showed convergence.
But when I go to see the contour and velocity vector, I found that the solution is not correct because the flow does not go from inlet to outlet fully. In the middle of the 2D domain the velocity magnitude is shown zero. Some velocity is found at the region close to inlet and outlet. The rest of the part has no flow. |
|
April 4, 2024, 02:56 |
|
#4 |
Senior Member
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 676
Rep Power: 21 |
How is the statement 'it all works' compatible with your results?
When you simulate without combustion, only air flow, what does the result look like then? If you place a monitor in the middle of your domain and observe pressure, velocity and species mass fractions, do you observe that all variables have reached their steady states? |
|
Tags |
combustion model, laminar flow, methane-air, species transport model |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Intuition for why flow follows convex surfaces | lopp | Main CFD Forum | 47 | February 1, 2022 14:14 |
Drag Force Ratio for Flat Plate | Rob Wilk | Main CFD Forum | 40 | May 10, 2020 05:47 |
Will the results of steady state solver and transient solver be same? | carye | OpenFOAM Running, Solving & CFD | 9 | December 28, 2019 06:21 |
Optimize for laminar flow, assume it valid for turbulent flow? | Chander | CFX | 15 | November 6, 2011 06:06 |
Optimize for laminar flow, assume it valid for turbulent flow? | Chander | Main CFD Forum | 2 | October 24, 2011 08:43 |