Hello everyone, I came acro
I came across this inconsistency when trying to run a new combination of thermophysicalProperties:
The problem is that multiComponentMixture::cellMixture() does not calculate the mixture value but only takes the value of the first species.
(I posted a detailed description here: http://www.cfd-online.com/cgi-bin/Op...4276#POST24276 )
The reason seem to be different declarations in hConstThermo.H and janafThermo.H. In janafThermo.H the operators += , -= and *= are defined, in hConstThermo they're not. The cell mixture is calculated in multiComponentMixture.C from:
mixture_ += Y_[n][celli]/speciesData_[n].W()*speciesData_[n];
Changing this to:
mixture_ = mixture_ + Y_[n][celli]/speciesData_[n].W()*speciesData_[n];
solves the error. You have to do the same in patchFaceMixture.
Though this solves the problem, it is actually only a workaround for multiComponentMixture. I think the right solution would be to implement += , -= and *= in hConstThermo as well.
All this now seems quite straightforward but it actually took me well over a month to figure out.
Thanks to Tommaso Lucchini for his help that guided me into the right direction!!
Tommaso, can you have a go at
Tommaso, can you have a go at this - should not be a problem...
Dear Dominik, Thanks for th
Thanks for the detailed bug report. I have now added the missing operators to hConstThermo and eConstThermo and these changes will be included in a patched-release of 1.5 as soon as we have sorted out all the other major issues being reported.
P.S. I don't believe the *= op
P.S. I don't believe the *= operator is needed as it simply calls the inherited operator from the equationOfState.
Dear Henry, I am using the
I am using the similar declaration from hCombusionThermo (refer the first post from Dominik) for calculating mixture properties and didn't realize that there might a problem in reading the values. I am using OpenFOAM-1.4
Can you please enlist the missing operators in hConstThermo and eConstThermo?
These changes are in the files
These changes are in the files on our git repository, you could download that and pull out the files/changes you need.
Dear Henry, Thanks for you
Thanks for your reply. I downloaded the files from the git repo. But, there are no documents which mentions the specific changes in the library. I had been working with 1.4.1 and am a bit reluctant to try 1.5 now. I wish somehow you can help me modify the library.
|All times are GMT -4. The time now is 19:37.|