[FLASH-BUGS] multigamma EOS Version: FLASH 2.1.20020605

[Alan Calder]

Hi Markus,

Thanks for sending a corrected version of the multigamma eos
routine. I will take a look at it and let you know what I
learn.

I should mention, though, that we are about to release FLASH 2.2,
which has several bug fixes from 2.1. I would recommend moving to
that when it is released, which will be later this week.

Thanks for sending your report.

Regards,

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-----
> 
>