Source code for idaes.surrogate.helmet.AncillaryEquations

##############################################################################
# Institute for the Design of Advanced Energy Systems Process Systems
# Engineering Framework (IDAES PSE Framework) Copyright (c) 2018-2019, by the
# software owners: The Regents of the University of California, through
# Lawrence Berkeley National Laboratory,  National Technology & Engineering
# Solutions of Sandia, LLC, Carnegie Mellon University, West Virginia
# University Research Corporation, et al. All rights reserved.
#
# Please see the files COPYRIGHT.txt and LICENSE.txt for full copyright and
# license information, respectively. Both files are also available online
# at the URL "https://github.com/IDAES/idaes-pse".
##############################################################################
"""
Modeling for saturated densities and vapor pressure
"""

import sys
from idaes.surrogate import alamopy
import importlib

from . import DataImport
from . import DataManipulation


molecule = ""
max_time = 1000



global DLfun, DVfun, PVfun
DLfun, DVfun, PVfun = None, None, None

DataToWrite = []


[docs]def DL(): "ALAMO regression of Saturated Liquid Density" global DLfun DataImport.DL(molecule) Values = DataImport.DLValues xval, zval = [], [] for x in Values: manVal = DataManipulation.DL(x) xval.append(manVal[0]) zval.append(manVal[1]) res = alamopy.alamo( xval, zval, zlabels=["%sDL" % molecule], almname="%sDL" % molecule, monomialpower=(-1, 1, 2, 3), expfcns=1, savetrace=True, savescratch=True, ) if res is None: raise Exception("Model does not exist for Saturated Liquid Density") DLfun = importlib.import_module("%sDL" % molecule, "f")
[docs]def getDL(): """ Imports the regressed saturated liquid density function """ global DLfun if DLfun is None: DLfun = importlib.import_module("%sDL" % molecule, "f") return DLfun
[docs]def DV(): """ ALAMO regression of saturated vapor density """ global DVfun DataImport.DV(molecule) Values = DataImport.DVValues xval, zval = [], [] for x in Values: manVal = DataManipulation.DV(x) xval.append(manVal[0]) zval.append(manVal[1]) res = alamopy.alamo( xval, zval, zlabels=["%sDV" % molecule], almname="%sDV" % molecule, monomialpower=(1, 2, 3, 4, 5, 6), savetrace=True, ) # print(res['model'], res['ssr'], res['R2']) if res is None: raise Exception("Model does not exist for Saturated Vapor Density") DVfun = importlib.import_module("%sDV" % molecule, "f")
[docs]def getDV(): """ Imports the regressed saturated vapor density function """ global DVfun if DVfun is None: DVfun = importlib.import_module("%sDV" % molecule, "f") return DVfun
[docs]def PV(): "ALAMO regression of vapor pressure" global PVfun DataImport.PV(molecule) Values = DataImport.Values xval, zval = [], [] for x in Values: manVal = DataManipulation.PV(x) xval.append([manVal[0], manVal[1]]) zval.append(manVal[2]) res = alamopy.alamo( xval, zval, zlabels=["%sPV" % molecule], almname="%sPV" % molecule, monomialpower=(1, 2, 3, 4, 5, 6), ) # print(res['model'], res['ssr'], res['R2']) if res is None: raise Exception("Model does not exist for Vapor Pressure") PVfun = importlib.import_module("%sPV" % molecule, "f")
[docs]def getPV(): """ Imports the regressed vapor pressure function """ global PVfun if PVfun is None: PVfun = importlib.import_module("%sPV" % molecule, "f") return PVfun