水のプランク熱力学ポテンシャルをプロットする 

さまざまな温度と圧力における水のエントロピー，内部エネルギー，比容積を得る．

In[1]:=

X temperature = UnitConvert[ Quantity[Table[x, {x, 20., 400, 20}], "DegreesCelsius"], "Kelvins"]; pressure = Quantity[Table[10^x, {x, 3, 6, 0.2}], "Pascals"]; opt = {"Temperature" -> temperature, "Pressure" -> pressure}; s = ThermodynamicData["Water", "Entropy", opt]; u = ThermodynamicData["Water", "InternalEnergy", opt]; v = ThermodynamicData["Water", "SpecificVolume", opt];

関係 を使って離散値を生成し，プランクのポテンシャルあるいはギブズの自由エントロピーに対する3D補間を実行する．

In[2]:=

X ifun = Interpolation[Flatten[Table[{{temperature[[j]], pressure[[i]]}, s[[i, j]] - pressure[[i]]*(v[[i, j]]/temperature[[j]]) - u[[i, j]]/temperature[[j]]}, {i, Length[pressure]}, {j, Length[temperature]}], 1]];

補間された関数を，温度と圧力の関数としてプロットする．

完全なWolfram言語入力を表示する入力を隠す

In[3]:=

X fun[t_?NumericQ, p_?NumericQ] := QuantityMagnitude[ Quiet[ifun[Quantity[t, "DegreesCelsius"], Quantity[p, "Pascals"]] ] ]; Plot3D[fun[tt, p*10^3]/1000, {tt, 20, 400}, {p, 1, 10^3}, AxesLabel -> {Style["Deg. Celsius", 14], Style["KPa.", 14], Style["KJ/(kg.K)", 14]}, Mesh -> None, ImageSize -> Medium, Filling -> Axis]

Out[3]=