CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   Sandia Flame D OpenFoam tutorial case - bad result? (https://www.cfd-online.com/Forums/openfoam-programming-development/220375-sandia-flame-d-openfoam-tutorial-case-bad-result.html)

Brandani September 5, 2019 06:31

Sandia Flame D OpenFoam tutorial case - bad result?
 
3 Attachment(s)
Hi there!

I am trying to study combustion modelling within OpenFoam and am currently looking into the Sandia Flame D case - /opt/openfoam6/tutorials/combustion/reactingFoam/RAS/SandiaD_LTS

I copy the case into my own directory and execute provided Allrun script.
The results of the temperature field at the last time step are shown in the first picture (myD).

However (to me) - these results seem to be incorrect. If you google Sandia D flame tutorial in openFoam you get these screenshots that make more sense in terms of results.

For example from the following thesis(http://repozitorij.fsb.hr/9279/1/Del..._diplomski.pdf), page 79 (picture D2) - the results are very different yet according to the information there the conditions are identical to those found in the OpenFoam tutorial case.

Also there is a paper (https://link.springer.com/article/10...494-014-9561-5) looking into Sandia D flames and their OpenFoam simulation looks as follows (D3).

So my question is do I need to work with the case and modify to achieve these more sensible results or is it just a simple question of post-processing -as I assume the in-built openFoam tutorial on this case is correct.

Thank you!

Brandani September 6, 2019 10:37

Update on this -- I installed the most recent dev version of OpenFoam, which contains the updated Sandia Flame D tutorial and produces much better results.

cryabroad September 11, 2019 03:35

Hi Brandani,

Would it be possible for you to share the updated results you mentioned? We are implementing the flamelet/progress variable approach in OpenFOAM, and obviously Flame D is one of the most used benchmark cases. I wonder what the results look like using the existing solvers in OpenFOAM.

Thanks!

Brandani September 11, 2019 05:10

Quote:

Originally Posted by cryabroad (Post 744419)
Hi Brandani,

Would it be possible for you to share the updated results you mentioned? We are implementing the flamelet/progress variable approach in OpenFOAM, and obviously Flame D is one of the most used benchmark cases. I wonder what the results look like using the existing solvers in OpenFOAM.

Thanks!

So just to clarify:

1) What I was saying is that SandiaD flame in OpenFOAM version 6 seemed to produce inaccurate results (I believe there was no ignition with the EDC combustion model - 1st picture in my first post)

2) Then, what I found is that the OpenFOAM development version (https://openfoam.org/version/dev/) + OpenFoam v7 (https://openfoam.org/release/7/) must have had their in-built tutorial on Sandia Flame D updated because there was ignition and results seem to match those pictures (2nd and 3rd) I attached in my first post

3) So install v7 or the dev version of OpenFOAM and you will have a working base case. I did not do any modifications to the case in these versions.

Hope this helps!

Brandani September 11, 2019 05:14

1 Attachment(s)
Quote:

Originally Posted by Brandani (Post 744430)
So just to clarify:

1) What I was saying is that SandiaD flame in OpenFOAM version 6 seemed to produce inaccurate results (I believe there was no ignition with the EDC combustion model - 1st picture in my first post)

2) Then, what I found is that the OpenFOAM development version (https://openfoam.org/version/dev/) + OpenFoam v7 (https://openfoam.org/release/7/) must have had their in-built tutorial on Sandia Flame D updated because there was ignition and results seem to match those pictures (2nd and 3rd) I attached in my first post

3) So install v7 or the dev version of OpenFOAM and you will have a working base case. I did not do any modifications to the case in these versions.

Hope this helps!

This is the T field at last time step produced by the in-built case in the dev version.

FluentStarter August 27, 2021 01:48

Does anyone have an explanation for the results in v6? One gets actually the same "bad" result in v2106.

Amirreza_pro January 7, 2022 11:41

1 Attachment(s)
Hi dear Brandani
I would appreciate if you could answer my question.
Actually I'm using SandiaD_LTS test case, but the time unite which OpenFOAM is using does not make sense me at all. Is the unite of time in second? Because in the default setting, 'end time' is set to 1500, but this time for combustion (injection) is not a reasonable value ! In combustion the time is ususallly in millisecond scale.
Do you have probably any Idea about time units in Openfoam, specially by using this start time and end time and also Delta T.
In attachment you can also find the picture of the ControlDict.orig file.

Thanks alot in advance

kerim January 18, 2022 01:10

Sandia Flame D OpenFoam tutorial case - bad result?
 
All quantities in OpenFOAM should be in SI system of unit.

dlahaye January 18, 2022 13:29

End time 1500 here refers to the maximal number of time steps (intended as number of iterations).

The value to the time step (or time increment) in controlled by the PIMPLE dictionary inside system/fvSolution


All times are GMT -4. The time now is 04:26.