Faça o gráfico do potencial termodinâmico de Planck da água 

Obtenha a entropia, a energia interna e o volume específico da água a diferentes temperaturas e pressões.

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 a relação para gerar valores discretos e realizar interpolação 3D para o potencial de Planck ou a entropia livre 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]];

Faça o gráfico da função interpolada como uma função de temperatura e pressão.

mostre o input completo de Wolfram Languageoculte o input

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