Grafique el potencial termodinámico de Agua de Planck  

Obtenga la entropía, energía interna y el volumen específico del agua a diversas temperaturas y presiones.

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

Use la relación para generar valores discretos y realizar una interpolación en 3D para el potencial de Planck o entropía libre de Gibbs.

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

Grafique la función interpolada como una función de la temperatura y la presión.

muestre la entrada completa de Wolfram Languageoculte la entrada

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