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.