Thermodynamische Werte aus interpolierten Daten berechnen 

Ermitteln Sie die Volumendaten für bei konstantem Druck und konstanter Temperatur aus ThermodynamicData.

Die gesamte Wolfram-Language Eingabe zeigenEingabe verbergen

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}]}]];

Interpolieren Sie die Daten und ermitteln Sie und .

In[2]:=

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

Stellen Sie die geschätzten Ableitungen graphisch dar.

Die gesamte Wolfram-Language Eingabe zeigenEingabe verbergen

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

Verwenden Sie die interpolierten Daten, um die isobare Wärmekapazität aus der isochoren Wärmekapazität zu berechnen.

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

Vergleichen Sie das berechnete Ergebnis mit den verfügbaren Daten.

In[6]:=

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

Out[6]=