Calculating the dissociation chemistry and equation of state properties of materials [Slides]
Many experimental techniques are accurate enough to reveal glaring weaknesses in the material models used by physics codes to design them. This is requiring ever more stringent equation of state (EOS) development. Formerly, EOS models were built using limited gas gun shock data and a few general physics assumptions. With the ever increasing computer power available, we began using density functional theory (DFT) to approximate other EOS regimes such as the principal isentrope, critical point, and sometimes melt. With Sequoia, and only because of Sequoia’s capacity, we are now calculating a limited range, general EOS table and then deriving the above EOS properties based on that table, often on a shorter time frame and more thoroughly than we used to compute just the EOS sections listed above. The first project was to calculate a limited range EOS for CO 2 and specifically track the dissociated species present at each density and temperature. We would not see all of the possible species due to our small sample size and limited time, however, we should see the dominant species and general trends. This work was done to support a new, chemistry aware EOS model that is still under development at LANL. The second project was done for two reasons. The first is LANL and SNL are developing a new amorphous SiO 2 EOS for use in hydrocodes. There is extensive principal Hugoniot and isentrope data, but very little lower density data. There is also no real estimate of the critical point so we found the approximate value of that as well. The second reason for lower density SiO 2 is to test a new method for calculating the Hugoniot of porous materials. Because aerogel (low density SiO 2 ) is used extensively in many shock experiments, there was a large data set to compare with. As can be seen in the figure, the new method appears to work well and is tracking the porous Hugoniot within the experimental error bars (not shown). Finally, SNL is developing a new EOS for LiF because we were seeing discrepancies with the SESAME table. DFT has some issues calculating EOS values for LiF, so we will make some corrections, however, our first attempt appears to have very good results which confirm several experimentally seen properties, as seen on the previous slide.