###############################################################################
# WaterTAP Copyright (c) 2021, The Regents of the University of California,
# through Lawrence Berkeley National Laboratory, Oak Ridge National
# Laboratory, National Renewable Energy Laboratory, and National Energy
# Technology Laboratory (subject to receipt of any required approvals from
# the U.S. Dept. of Energy). All rights reserved.
#
# Please see the files COPYRIGHT.md and LICENSE.md for full copyright and license
# information, respectively. These files are also available online at the URL
# "https://github.com/watertap-org/watertap/"
#
###############################################################################
"""1D reverse osmosis examples"""
from pyomo.environ import ConcreteModel
from idaes.core import FlowsheetBlock
from idaes.core.util.scaling import calculate_scaling_factors
from watertap.examples.flowsheets.full_treatment_train.model_components import (
property_models,
)
from watertap.unit_models.reverse_osmosis_1D import (
ReverseOsmosis1D,
ConcentrationPolarizationType,
MassTransferCoefficient,
PressureChangeType,
)
from watertap.examples.flowsheets.full_treatment_train.util import (
solve_block,
check_dof,
)
[docs]def build_RO(m, base="TDS", level="simple", name_str="RO"):
"""
Builds a 1DRO model at a specified level (simple or detailed).
Requires prop_TDS property package.
"""
if base not in ["TDS"]:
raise ValueError(
"Unexpected property base {base} for build_RO" "".format(base=base)
)
prop = property_models.get_prop(m, base=base)
if level == "simple":
raise ValueError(
"Unexpected RO level {level} for build_RO" "".format(level=level)
)
elif level == "detailed":
# build unit
setattr(
m.fs,
name_str,
ReverseOsmosis1D(
property_package=prop,
has_pressure_change=True,
pressure_change_type=PressureChangeType.calculated,
mass_transfer_coefficient=MassTransferCoefficient.calculated,
concentration_polarization_type=ConcentrationPolarizationType.calculated,
transformation_scheme="BACKWARD",
transformation_method="dae.finite_difference",
finite_elements=10,
),
)
blk = getattr(m.fs, name_str)
# specify unit
blk.area.fix(50)
blk.A_comp.fix(4.2e-12)
blk.B_comp.fix(3.5e-8)
blk.mixed_permeate[0].pressure.fix(101325)
blk.channel_height.fix(1e-3)
blk.spacer_porosity.fix(0.97)
blk.N_Re[0, 0].fix(500)
else:
raise ValueError(
"Unexpected argument {level} for level in build_RO" "".format(level=level)
)
def solve_RO(base="TDS", level="simple"):
m = ConcreteModel()
m.fs = FlowsheetBlock(dynamic=False)
property_models.build_prop(m, base="TDS")
build_RO(m, base=base, level=level)
# specify feed
property_models.specify_feed(m.fs.RO.feed_side.properties[0, 0], base="TDS")
m.fs.RO.feed_side.properties[0, 0].pressure.fix(50e5)
# scaling
calculate_scaling_factors(m)
# initialize
m.fs.RO.initialize(optarg={"nlp_scaling_method": "user-scaling"})
m.fs.RO.display()
check_dof(m)
solve_block(m)
m.fs.RO.report()
return m
if __name__ == "__main__":
# solve_RO(base='TDS', level='simple')
solve_RO(base="TDS", level="detailed")