[FLASH-BUGS] multigamma EOS Version: FLASH 2.1.20020605

[Alan Calder]

Hello again, Markus,

It looks like you did indeed find a bug in the code. Thanks for
the report. The formulae were correct in the eos1d.F90 version,
which is the one that cranks the eos on the data row by row.
But, one of our tests does use the point by point version and
we should have caught that.

I am not sure if the correction made it into the released version, 
but it has been committed to the source repositiory.

Thanks again,

Alan





Alan C. Calder
University of Chicago
ASCI Flash Center
5640 South Ellis Ave.
Chicago, IL  60637

email: calder at flash.uchicago.edu
phone:  (773) 834-3904
fax:  (773) 834-3230
web: http://flash.uchicago.edu/~calder/

On Tue, 3 Dec 2002, Markus Gross wrote:

> Hi!
> 
> The multigamma equation of state that came with my version of FLASH 
> (2.1.20020605) seems to have a bug. I only checked pressure, temperature and 
> internal energy so far and found that:
> 
> pressure for input=2 is wrong (I presume)
> ener     for input=3 is wrong (I presume again)
> 
> I attached my "corrected" version to this email.
> 
> I am a bit careful about stating that there may be a mistake because I wonder 
> why it did not crop up earlier during one of your tests - so it may be my 
> fault. Then again you may be more interested in using the Helmoltz eos anyway 
> and therefore may not notice ...
> 
> In my case the riemann solver bailed out - which may be not surprising if 
> there is no thermodynamic correct state to start with...
> 
> Please let me know if I was right or wrong.
> 
> 
> Regards,
> 
> Markus.
> 
> 
> -- 
> _______________________________________________________________
> 
> Markus Gross AMIMechE BEng (Hons.) Mechanical Engineering
> 
> Heriot Watt University Edinburgh
> Department of Mechanical and Chemical Engineering
> 
> Associate member of the Institution of Mechanical Engineers
> Member of the SPIE
> _______________________________________________________________
> 
> further contact:
> 
> Phone   : +44 (0) 131 449 5111 etx. 4737
> 
> UNiX talk: talk markus at lasersim.mce.hw.ac.uk
> 
> _______________________________________________________________
> 
> "Plans are a place to begin," Grove said. "They rarely deliver
> you to where you expect. Make your plans knowing you are going
> to throw them away."
> 
> _______________________________________________________________
> 
> 
> -----BEGIN PGP PUBLIC KEY BLOCK-----
> Version: GnuPG v1.0.6 (GNU/Linux)
> Comment: For info see http://www.gnupg.org
> 
> mQGiBDt30nARBACUL2XadcLU9YS+3ndnWUUpnrhAD3Ab/TqJP1Obtz7H356+vigu
> 2m3exw2P9fky08SBdTbjOCV5OTmVcRW5F0SNMTIl8gXtGdIcQHCcZxi1GcRJtTBj
> npfRCmZojHr8YgsKDdJ9xupNm+UYWcg1BpFX7xe65kXM4oVwVyuCcUT+VwCg6/Hd
> 6oVoRHm1eq1NdMbCfsc9x78D/2uGQSMJBKewQ/uSPpdEg+bkAQo5Y7MooEeYpLk9
> Iroo4rqxp12k0+bWHWVMdfywI80tivorTIgBTKpCtjLaiKen9E9zECo+2OR4+XRQ
> O6cBKXOm+V2lJOpFYOh500tJHPggshD21jksNPKJr9CM1byeZFw38fFJez+DyShZ
> Gmc6A/9L21jqbVjxvAPwaqe09GEmpT1EYOWtcPYSh2cuo1v6XsGFXfi8tnRSizhH
> xsyIfuOUuj0NKRIcXH60awrvABSa45fvQUefPw0hxSwdCnfc0yEqcD4x/03xTZJS
> 1hUK3T4oTlrXNugaCb59ezciEIh2y/d8pW1DR3RMcgXUSz39tLQhTWFya3VzIEdy
> b3NzIDxtLnMuZ3Jvc3NAaHcuYWMudWs+iFcEExECABcFAjt30nAFCwcKAwQDFQMC
> AxYCAQIXgAAKCRBTnmsr4a/0+e3tAJoDk08MfeqB83FXP6g8QlDyTlotSwCgiKwS
> PPkVAbopX8mahpkrphEaunW5Ag0EO3fSghAIAOeLTMjdFkl/xBvjgD2N0n4DCK4N
> EVfaXmkJK5dbTM3K6ZsJjILZ/TfC6R8oxYPsYKwkwLJJEe4E/h1TFg7ogDYbViW8
> DihOyeN8zKSymLTbbikEJX8aNeCdBRuy87X0pBazDq6//+oTIatS2EsvMyKk0tEZ
> Mh0vTmUJqoeorQeAKH66lPsH2A8QI6u0zv17t4S2aLNBXIbKgOJwsglZ1IwjwowR
> KEGsLsHInNkP4bnlBjd51RMiJT979QMPmNgLGEId9Vi0oijtF/v6J91GtH++G6WE
> HM08DsP/BKjRgO9ivWP7lOHSDCJpGYti5K6g38JAc1H2EoHGLD7eRgI/kR8AAwUH
> /j6NVGLxnFGDHGcBXEb2Y4Wwhf/eUeE8Wit0qqRoNypBuzvWl96GXSPj/saQA/Ri
> XcFWR1G3E4V/zdxYnU/xvvGADj9jKHpZjAs73fJJxLc85gzHLcx00tHcQVA07Ihr
> hao8WY+BtKhkZorYEJEvSDJ9xufjL/o1t40eCsHQ3b9Jh/1lqQ2KT1lhDXAaDo1P
> pE38E2uImC0umIW/UtuKJ1k52nXx2qulfcTjAHd35XidwGreH3aHy3fV+x/+wW1T
> CQBHy27dG7JL3ESzPiKW0w2miWjtd4k5jqJTfgKI4wFUhAShxTjp8M051ArzopRJ
> iX+s7RrHuEXHjh45V1zU5BWIRgQYEQIABgUCO3fSggAKCRBTnmsr4a/0+U/rAJ9o
> lOhPGltVnRncYJPzyIk00H0JDQCfWrlnYDh5cRQmhnHilnrVO5Ykr24=
> =i3eC
> -----END PGP PUBLIC KEY BLOCK-----
> 
>