TAGGED: laminar-flame, soot
February 15, 2023 at 2:14 pmChen WeiSubscriber
I'm trying to simulate laminar diffusion flame in ANSYS Fluent within radiation and soot models. Within the Moss-Brookes model, I always get an incorrect soot volume fraction profile, even I set the case almost the same, in comparison with a reference paper "On the relative contributions of soot to radiative heat transfer at different oxygen indices in ethylene – O2/CO2 laminar diffusion flames". Some differences are:
- Due to the limited computational force. The mesh is not the fine (near 30,000 for 2D axis-symmetric)
- The thermal and transport properties are obtained from GRI-Mech 3.0, because the properties used in the paper are not mentioned clearly.
- I'm using the ANSYS Fluent Student version, while the reference paper uses 19.1.
The soot volume fraction profile in my case always has a peak in the axis of the flame, while the profile from others' experiment or simulation data should have a peak in the flame's fringe.
I even get a case from someone, who claims that the case could produce a correct soot volume fraction "shape". However, I still get the wrong shape with that case, even if I don't do anything but read the case and run it (300 steps for a stable field and then ignite by patching the field with 1,300 K).
I think I must miss something during the setup. But I really cannot find it out.
February 16, 2023 at 2:26 pmRobAnsys Employee
The version shouldn't matter as Student is exactly the same code, the only difference is there is a cell count limit of 512k.
It's more likely that boundary conditions, model settings or mesh aren't the same, or aren't suitable. If you look at the results how do the other fields compare with the paper? How well converged is the case?
February 20, 2023 at 8:07 amChen WeiSubscriber
Thank you for your response and for taking the time to look into my issue further!
Here is some additional information for you.
In fact, the case I receive is just from the paper’s author. She gave me two files, case.dat and data.dat, but she didn’t give chem.inp, trans.dat and therm.dat. She said she create three files based on GRI-Mech 3.0 by removing the NOx part. So I create the files related to the chemical mechanism by myself based on her instruction.
Following the author's guidance, I attempted to replicate her work using the case.dat files she provided. My replication process included the following steps:
- Importing the case.dat file into ANSYS Fluent.
- Addressing an error that occurred due to the absence of the thermal and transport database by importing the chemical mechanism and database myself.
- Running the simulation without making any modifications to the imported case.
- Performing 1000 steps to achieve a cold stable field.
- Igniting the flame by patching the entire mesh with a temperature of 1300 K.
- Allowing the simulation to run until the flame reaches full coverage.
I was able to compare my results with the reference data provided by the author in the data.dat file. Overall, most of the fields in my simulation are similar to the reference, including some minor species that are crucial for soot formation, such as C2H2, C2H4, OH, and O. However, there are some discrepancies between my soot profiles and the reference data. Specifically, the shape of the soot volume fraction profile is not matching with the reference data, which has a peak in the flame's fringe rather than at the axis of the flame, as mentioned in the last reply.
We have the same Rate_of_Nucleation and Rate_of_Mass_Nucleation profiles, which is reasonable because we share the same temperature and pressure fields. However, there are some differences in the rest of the results. The source term, normalized nuclei concentration, and soot mass fraction all differ between my simulation and the author's simulation. Specifically, my results appear more diffusive compared to hers, which have much sharper profiles.
I suspected that the soot model parameters might be different in the two cases, which could explain the differences in the results. To test this hypothesis, I followed these steps:
- Imported the case.dat and data.dat files.
- Solved only the soot and nuclei equations.
- Set the relaxation factor for both variables to be 1e-5.
- Ran only one step.
If the two cases had different parameter sets, then the source terms would have changed after this step. However, I found that the source terms remained the same.
I am completely perplexed and unable to identify any other factors that could be contributing to the differences between my results and hers. I appreciate your prompt attention to this matter and eagerly anticipate your reply.
March 20, 2023 at 6:50 amChen WeiSubscriber
Do you have any idea about the results?
March 20, 2023 at 9:11 amRobAnsys Employee
Apologies, I missed this. Thanks for the reminder.
If you set the relaxation factors that low not much will change, but the creation term is only part of the calculation: there's removal too. I'm not sure why it's different, but suspect there is a slight difference somewhere in the chemistry.
- You must be logged in to reply to this topic.
Boost Ansys Fluent Simulations with AWS
Computational Fluid Dynamics (CFD) helps engineers design products in which the flow of fluid components is a significant challenge. These different use cases often require large complex models to solve on a traditional workstation. Click here to join this event to learn how to leverage Ansys Fluids on the cloud, thanks to Ansys Gateway powered by AWS.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
- Suppress Fluent to open with GUI while performing in journal file
- Floating point exception in Fluent
- What are the differences between CFX and Fluent?
- Heat transfer coefficient
- Getting graph and tabular data from result in workbench mechanical
- The solver failed with a non-zero exit code of : 2
- Difference between K-epsilon and K-omega Turbulence Model
- Time Step Size and Courant Number
- Mesh Interfaces in ANSYS FLUENT
- error in cfd post
© 2023 Copyright ANSYS, Inc. All rights reserved.