private void CalculateEnthalpyFromTemperature(double pValue, double tValue, double moistureContentValue) { //double tBoilingPointOfSolvent = MoistureProperties.GetSaturationTemperature(pValue); double tBoilingPointOfSolution = GetBoilingPoint(pValue); double hValue; if (tValue <= tBoilingPointOfSolution) { hValue = CalculateLiquidEnthalpy(pValue, tValue, moistureContentValue); Calculate(specificEnthalpy, hValue); Calculate(vaporFraction, 0.0); } else if (tValue > tBoilingPointOfSolution) { double hBubble = GetBubblePointEnthalpy(pValue, moistureContentValue); double evapHeat = GetEvaporationHeat(tBoilingPointOfSolution); SteamTable steamTable = SteamTable.GetInstance(); SubstanceStatesAndProperties props = steamTable.GetPropertiesFromPressureAndTemperature(pValue, tBoilingPointOfSolution + 1.0e-6); double hVapor1 = props.enthalpy; props = steamTable.GetPropertiesFromPressureAndTemperature(pValue, tValue); double hVapor2 = props.enthalpy; double cs = GetCpOfAbsoluteDryMaterial(tValue); hValue = hBubble + moistureContentValue * (evapHeat + (hVapor2 - hVapor1)) + (1.0 - moistureContentValue) * cs * (tValue - tBoilingPointOfSolution); Calculate(specificEnthalpy, hValue); Calculate(vaporFraction, moistureContentValue); } }
private double CalculateSpecificEnthalpyVaporOutlet(double tBoilingPointOutlet) { //SteamTable steamTable = SteamTable.GetInstance(); //SubstanceStatesAndProperties props = steamTable.GetPropertiesFromPressureAndTemperature(vaporOutlet.Pressure.Value, tBoilingPointOutlet); //return props.enthalpy; double hValue = 0; if (inlet is DryingMaterialStream) { DryingMaterialStream dmsVaporOutlet = vaporOutlet as DryingMaterialStream; double pValue = vaporOutlet.Pressure.Value; double hBubble = dmsVaporOutlet.GetBubblePointEnthalpy(pValue, 1.0); double tBoilingPoint = dmsVaporOutlet.GetBoilingPoint(pValue, 0); double evapHeat = vaporOutlet.GetEvaporationHeat(tBoilingPoint); SteamTable steamTable = SteamTable.GetInstance(); SubstanceStatesAndProperties props = steamTable.GetPropertiesFromPressureAndVaporFraction(pValue, 1.0); double hVapor1 = props.enthalpy; props = steamTable.GetPropertiesFromPressureAndTemperature(pValue, tBoilingPointOutlet + 1.0e-6); double hVapor2 = props.enthalpy; hValue = hBubble + evapHeat + (hVapor2 - hVapor1); } return(hValue); }