Commit 53d50dfd authored by Timo Koch's avatar Timo Koch
Browse files

[simpleh2o][bugfix] Fix liquid enthalpy for incompressible fluid: dh = cp*dT + p/rho

parent d497f147
......@@ -149,7 +149,8 @@ public:
Scalar pressure)
{
static const Scalar tRef = getParam<Scalar>("SimpleH2O.ReferenceTemperature", 293.15);
return liquidHeatCapacity(temperature, pressure)*(temperature - tRef);
return liquidHeatCapacity(temperature, pressure)*(temperature - tRef)
+ pressure/liquidDensity(temperature, pressure);
}
/*!
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment