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