#------------------------------------------------------------------------------- # M/ANEOS improved SiO2 #------------------------------------------------------------------------------- try: izetl = sph.vector_of_int(1, -1) sph.initializeANEOS('/proj/nmovshov_hindmost/collisions/ANEOS/ANEOS.INPUT', 'ANEOS.barf', izetl) etamin, etamax = 0.94, 10 rho0 = 2650 pext, pmin, pmax = 0.0, -1e200, 1e200 # these are actually the defaults SiO2 = sph.ANEOS(0, # Material number 1000, # num rho vals 1000, # num T vals etamin*rho0, # minimum density (kg/m^3) etamax*rho0, # maximum density (kg/m^3) 1.0, # minimum temperature (K) 1.0e4, # maximum temperature (K) units, pext, pmin, pmax) os.system('rm -f ANEOS.barf') del izetl, etamin, etamax, rho0 except: print "ANEOS SiO2 not loaded" pass #------------------------------------------------------------------------------- # A polytropic fluid EOS #------------------------------------------------------------------------------- K = 2e5 # polytropic constant n = 1 # polytropic index mu = 2.2e-3 # mean molecular weight