Lead sphere implosion results(Pu239 bare core next?)
I want to share the results i got with my hydrocode on the spherical implosion of a 10cm bare lead sphere by a 1Mbar initial external pressure. The 2 snapshots correspond to the pre-reflection and the post-reflection of the converging shock wave. We can see pressures of over 32g/cm^3 near the center.
This is just a small proof of concept where i wanted to test my implementation against this case which is well documented. The obvious next step would be to do the same but with a bare Pu239 sphere. The problem is i don't have the EOS data and can't seem to find it on the internet for obvious reasons. The SESAME library is out of question since i'm not affiliated with a research institution. If i manage to solve the Pu239 implosion and any of you are interested i might consider sharing the code so you guys can experiment with different initial pressures, core radius, initial density etc... A major big step would be to couple the implosion code with a neutronics module to experiment with criticality insertion calculations but that's not an easy task.