def pressure(self, temperature, volume, params):
"""
Returns pressure [Pa] as a function of temperature [K] and volume[m^3]
EQ B7
"""
return bm.birch_murnaghan(params['V_0']/volume, params) + \
self.__thermal_pressure(temperature,volume, params) - \
self.__thermal_pressure(300.,volume, params)
def volume(self, pressure,temperature,params):
"""
Returns volume [m^3] as a function of pressure [Pa] and temperature [K]
EQ B7
"""
func = lambda x: bm.birch_murnaghan(params['V_0']/x, params) + \
self.__thermal_pressure(temperature, x, params) - \
self.__thermal_pressure(300., x, params) - pressure
V = opt.brentq(func, 0.5*params['V_0'], 1.5*params['V_0'])
return V
def pressure(self, temperature, volume, params):
return bm.birch_murnaghan(params['V_0']/volume, params) + \
self.__thermal_pressure(temperature,volume, params) - \
self.__thermal_pressure(300.,volume, params)
def pressure(self, temperature, volume, params):
return (
bm.birch_murnaghan(params["V_0"] / volume, params)
+ self.__thermal_pressure(temperature, volume, params)
- self.__thermal_pressure(300.0, volume, params)
)
def pressure(T, V, params):
return bm.birch_murnaghan(params['ref_V']/V, params) + \
mgd_thermal_pressure(T,V, params) - \
mgd_thermal_pressure(300,V, params)
def volume(pressure,T,params):
func = lambda x: bm.birch_murnaghan(params['ref_V']/x, params) + \
mgd_thermal_pressure(T, x, params) - \
mgd_thermal_pressure(300, x, params) - pressure
V = opt.brentq(func, 0.5*params['ref_V'], 1.5*params['ref_V'])
return V
