Compute valores termodinâmicos a partir de dados interpolados 

Obtenha os dados de volume para a uma pressão e temperatura constantes com ThermodynamicData.

mostre o input completo de Wolfram Languageoculte o input

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 sobre os dados e obtenha e .

In[2]:=

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

Faça o gráfico das derivadas estimadas.

mostre o input completo de Wolfram Languageoculte o input

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]=

Use os dados interpolados para calcular a capacidade térmica isobárica a partir da capacidade térmica 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 o resultado calculado com dados disponíveis.

In[6]:=

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

Out[6]=