sprayFoam aachenBomb FATAL ERROR when running combustion
1 Attachment(s)
I've already submitted this to the OpenCFD Mantis bug report, but I thought I would post this here also, to see if anyone else has experienced this issue...
I'm running OpenFOAM-2.1.1 (but this also happens with version 2.1.0). I can run the sprayFoam/aachenBomb tutorial just fine. Then I turn on chemistry (in constant/chemistryProperities I change "chemistry off;" to "chemistry on;"), and run the case from scratch, and it also runs just fine. Then I turn on chemistry (as above) and combustion ( in constant/combustionProperities I change "active false;" to "active true;"), and the case experiences a "FOAM FATAL ERROR" in speciesThermoI.H after about 0.000769 seconds of simulated time. I've attached a copy of the log file (log.txt) with most of the intermediate steps deleted so it isn't so large. If anyone else has experienced this problem, please comment. If anyone knows what the problem is and how to fix it, please post a reply. Thanks. Phil :) |
I have no resolution, but I can just say that I see exactly the same in several setups. I could not pinpoint the issue yet. Only clue so far is that is happens as soon as T gets in range 1350-1400 Kin all cases. therm.dat differs only from other version (1.6-ext) on species not in use for this tutorial.
I'll get back if I find something. K |
2 Attachment(s)
Ok... did some digging. Attached are plots generated with this code:
Code:
for(label T=300; T<2400; T++) It seems there is a discontinuity in the enthalpy at ~1390 K, which likely cause the non-convergence. The fault is likely in therm.dat which have diffenrent values for the thermo compared to for instance ftp://ftp.technion.ac.il/pub/support...mics/THERM.DAT Changing therm.dat might help here K |
Hi
I placed the therm.dat file in the etc/thermoData folder, after renaming the old therm.dat, and it did not solve the problem... in which folder should i place it... ? in aachenBomb/chemkin...? BR Alex |
Hi Alex,
Quote:
As for the file that should be replaced, you need to check the file "constant/thermophysicalProperties" that is in your case folder! Best regards, Bruno |
Hi wyldckat
I changed the reference in thermophysicalProperties, pasted the file in chemkin folder, after removing the header... and launched the app... but it seesm that the therm.dat file is ill formed anyway... even after removing the header... : Code:
Ag cr REF ELEMENT T 6/12AG 1. 0. 0. 0.S 200.000 1235.080 1000. 1 the corresponding line in etc/thermoData/them.dat is as follows : (CH2O)3 70590C 3H 6O 3 G 0200.00 4000.00 1500.00 1 A slightly different format.... BR Alex |
Mmm... I guess you're meant only to replace the lines that you actual need...
|
Hello,
I have also experienced the problem mentioned above when turning on combustion and chemistry. I modified the original C7H16 in therm.dat with the one linked here: ftp://ftp.technion.ac.il/pub/support...mics/THERM.DAT and rerun the simulation only to receive the same error as before at around 0.00077 seconds. I am not sure whether I should change the coefficients for all the species within the default reaction scheme. I am using OpenFOAM-2.1.1, please could someone confirm whether this problem is present in the GIT version of OpenFOAM. Many Thanks, ris |
Greetings Ris,
Quote:
Best regards, Bruno |
1 Attachment(s)
Hello Bruno,
The case in question is the tutorial aachenbomb located in /tutorials/lagrangian/sprayFoam/. The only changes made to the case were: Turning chemistry on in constant/chemistryProperties and Turning combustion to true in contant/combustionProperties. I have attached an example case with the changes above. Thank you for your assistance, Regards, ris |
Hi Ris,
OK, so after dwelling a bit on this I finally figured out what Kalle was talking about. The following changes seem to work:
Now we come to the big problem: how on Earth is this going to be validated? The values present in OpenFOAM were possibly validated with some project, problem is with which one? Best regards, Bruno PS: the mantis bug report mentioned in the first post is this one: http://www.openfoam.org/mantisbt/view.php?id=574 |
Hi Bruno,
Thank you very much for this workaround, I will make the changes you have specified and update this thread with the outcome. With respect to validating this workaround, I should be able to get access to experimental data regarding diesel injection and will try to see how the simulation will compare. Its not ideal but it will give us an idea on what is happening. Many Thanks, ris |
Sorry Bruno for not clarifying enough what I meant with my post. You've got my idea anyhow.
I just got inspired so I decided to calculate Cp and H @Tcommon manually using the data in therm.dat retrieved July 14 from OpenCFD. You can paste the following into python to see for yourself: Code:
T=1391 How to use thermodata: http://www.me.berkeley.edu/gri_mech/data/nasa_plnm.html I have no clue right now what is wrong here... fishy. K |
The way janafThermo defines enthalpy for mixtures (weighting together their coefficients, see operators+= and operator+ in janafThermoI.H) does not allow for using different Tcommon for different species in the same mixture. Janaf should crash upon this in the case debug switches are activated. I did not manage to trigger that however. In this case; maybe C7H16 data is OK, but Tcommon for that specie (1391) is used also for mixtures containing also other species (in this case all the other species have Tcommon = 1000K). As those species have inconsistent Cp and H at 1391K for their two sets of coefficient, which causes non-convergence there. Indeed, the Cp-figure in my previous post shows reasonable Cp up to 1000K.
Conclusion: be careful not to mix thermodata's with different Tcommons. I know Cantera for instance, will complain when you use Chemkin thermo data with several Tcommons. OpenFOAM does not check this by default, but it causes obviously confusing problems. The other thermo data for C7H16 in this thread have Tcommon=1000K and the issue goes away. K |
Very nice explanation Kalle!
But now that I think of it, this looks familiar... I think this is still an incomplete feature in OpenFOAM which was discussed sometime ago here: http://www.openfoam.org/mantisbt/view.php?id=327 - with any hope, it'll be properly fixed in the next version of OpenFOAM... |
I did not know about that bug report. Thanks for point that out, Bruno! I have updated the other bug report about our findings here.
For now on, I'd guess we should stick to thermo data with only one Tcommon. K |
hello
I used janaf coefficients in openfoam but an error in thermo model displayed.i draw that in matlab and it showed a large discontinuity in 1000K.i used coefficients that come in openfoamwiki and don't know what is wrong.how you used the coeffs?i need air or nitropen and oxygen as a mixture. Thanks for helping. |
For more information about Ehsan's problem, check his more complete thread here: http://www.cfd-online.com/Forums/ope...incorrect.html
Best regards, Bruno |
thankyou wyldckat
|
Number of coeffs in therm.dat
Dear all,
Sorry, if a little bit off topic. Why etc/thermoData/therm.dat includes 15 coeffs for each specie? Why not 14 (7 for low + 7for high temperature)? And why in etc/thermoData/thermData imported only 14 from 15? Is it correct to use 15 coeffs if the OF accepts 14 only? |
All times are GMT -4. The time now is 00:20. |