Extend range of Water 97 fluid properties in FP module

[GiudGiud]

Motivation

Water97FluidProperties::T_from_p_h_ad() Region 5 is not implemented

the code assumes that the (p, h) plot is continuous across h_from_p_T(pressure, vaporTemperature(pressure)), which is not. Region 4 should be a two phase region

Design

@lingzou suggested the following pseudo code for h_from_pT

if (h > h_from_p_T(p, 273.15) && h < h_l_sat_from_p(p))
  return 1;
else if (h < h_g_sat_from(p))
  return 4;
else if (h < h_from_p_T(p, 1073.15)
  return 2;
else if (h < h_from_p_T(p, 2273.15)
  return 5:

and a bisection method for region 5

Impact

More capable fluid properties More capabilities for studying water cooled systems