Add a SpinerEosDependsRhoSie constructor that uses an existing eos object to generate a spiner eos#632
Add a SpinerEosDependsRhoSie constructor that uses an existing eos object to generate a spiner eos#632
Conversation
| dependsRhoSie_.P(j, i) = P; | ||
|
|
||
| // Bulk modulus will be computed by calcBMod_() after derivatives are populated | ||
| dependsRhoSie_.bMod(j, i) = robust::EPS(); |
There was a problem hiding this comment.
We can check that the eos provides an isentropic bulk modulus. Probably should use that. If it provides Ks, gamma and Cv should we derive the derivatives from those? Probably.
There was a problem hiding this comment.
modulo thermodynamic consistency...
There was a problem hiding this comment.
Are you suggesting that you prefer finite difference for derivatives and then evaluation of the properties from those? The real question I am thinking about right now, is what is the minimum we expect in the source_eos. I'm inclined to allow (for the dependsRhoSIE case) the source_eos to provide as little as PofRE and TofRE to fill in the tables, taking advantage of ks, gruneisen, and cv it it also has those methods.
There was a problem hiding this comment.
Honestly, I'm not sure. The advantage of finite differences is that they are guaranteed to be robust to thermodynamic consistency issues. The disadvantage is that they are never that accurate.
The worst possible scenario is using finite differences and then trying to use them in thermodynamic derivatives to compute something else. I ran into this in #526 where I didn't have a good mechanism to get P-T derivatives (hence #598). The finite differences there when combined into quantities we want (like a sound speed) could result in pretty large errors.
Some of the formulas from Menikoff and Plohr using these quantities are pretty complex. I'm not sure anybody has looked at the effects of thermodynamic inconsistencies on the accuracy of quantities computed from these quantities. I'd naively think that
If I had to pick, I think I'd agree with you on using
608002d to
522f37e
Compare
…nd an object representing grid parameters
522f37e to
cc4ee8f
Compare
PR Summary
This is currently very DRAFy. Almost completely generated by AI and I haven't fully reviewed the work myself, so if you don't have time to look at it right now, that's fine.
I had asked it to emulate sesame2spiner in terms of parameters for grid generation. Seems somewhat heavyweight, but there are defaults.
PR Checklist
make formatcommand after configuring withcmake.plan_historiesfolder, with a filename the same as the MR number.If preparing for a new release, in addition please check the following: