Calcule los valores termodinámicos desde datos interpolados  

Obtenga los datos de volumen de en una presión constante y una temperatura constante desde ThermodynamicData.

muestre la entrada completa de Wolfram Languageoculte la entrada

In[1]:=

X refTemperature = Quantity[40, "DegreesCelcius"]; refPressure = Quantity[10^5, "Pascals"]; With[{p = Quantity[Table[10^x, {x, 4, 6, 0.01}], "Pascals"]}, constantT = Transpose[{p, ThermodynamicData["CarbonDioxide", "SpecificVolume", {"Pressure" -> p, "Temperature" -> refTemperature}]}]]; With[{t = Quantity[Table[x, {x, 30, 50, 0.1}], "DegreesCelcius"]}, constantP = Transpose[{t, ThermodynamicData["CarbonDioxide", "SpecificVolume", {"Temperature" -> t, "Pressure" -> refPressure}]}]];

Interpole los datos y obtenga y .

In[2]:=

X SpVolAtConstantT = Interpolation[constantT]; SpVolAtConstantP = Interpolation[constantP]; dVdP[x_Quantity] := Derivative[1][SpVolAtConstantT][x]; dVdT[x_Quantity] := Derivative[1][SpVolAtConstantP][x];

Trace las derivadas estimadas.

muestre la entrada completa de Wolfram Languageoculte la entrada

In[3]:=

X Column[{ ListLinePlot[ Table[With[{tt = Quantity[10^p, "Pascals"]}, {tt, dVdP[tt]}], {p, 4, 6, 0.01}], Frame -> True, FrameLabel -> Automatic, PlotLabel -> "\!\(\*SubscriptBox[\((\*FractionBox[\(\[PartialD]V\), \(\ \[PartialD]p\)])\), \(T\)]\)", ImageSize -> Medium], ListLinePlot[ Table[{Quantity[t, "DegreesCelcius"], dVdT[Quantity[t, "DegreesCelcius"]]}, {t, 30, 50, 0.1}] , Frame -> True, FrameLabel -> Automatic, PlotLabel -> "(\!\(\*FractionBox[\(\[PartialD]V\), \ \(\[PartialD]T\)]\)\!\(\*SubscriptBox[\()\), \(p\)]\)", ImageSize -> Medium]}]

Out[3]=

Utilice los datos interpolados para calcular la capacidad calorífica isobárica desde la capacidad calorífica isocórica.

In[4]:=

X Cv = ThermodynamicData["CarbonDioxide", "IsochoricHeatCapacity", {"Temperature" -> refTemperature, "Pressure" -> refPressure}];

In[5]:=

X Cv - UnitConvert[refTemperature, "Kelvins"]* UnitSimplify[((dVdT[refTemperature]^2)/ dVdP[refPressure])*(Quantity[1, "Kilograms"]/ Quantity[1, "Moles"])]

Out[5]=

Compare el resultado calculado con los datos disponibles.

In[6]:=

X ThermodynamicData["CarbonDioxide", "IsobaricHeatCapacity", \ {"Temperature" -> refTemperature, "Pressure" -> refPressure}]

Out[6]=