FreeCAD / src /Mod /Fem /femobjects /element_fluid1D.py

Upload folder using huggingface_hub

985c397 verified 10 months ago

13.4 kB

# ***************************************************************************
# *   Copyright (c) 2016 Ofentse Kgoa <kgoaot@eskom.co.za>                  *
# *   Copyright (c) 2018 Bernd Hahnebach <bernd@bimstatik.org>              *
# *   Based on the FemElementGeometry1D by Bernd Hahnebach                  *
# *                                                                         *
# *   This file is part of the FreeCAD CAx development system.              *
# *                                                                         *
# *   This program is free software; you can redistribute it and/or modify  *
# *   it under the terms of the GNU Lesser General Public License (LGPL)    *
# *   as published by the Free Software Foundation; either version 2 of     *
# *   the License, or (at your option) any later version.                   *
# *   for detail see the LICENCE text file.                                 *
# *                                                                         *
# *   This program is distributed in the hope that it will be useful,       *
# *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
# *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
# *   GNU Library General Public License for more details.                  *
# *                                                                         *
# *   You should have received a copy of the GNU Library General Public     *
# *   License along with this program; if not, write to the Free Software   *
# *   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  *
# *   USA                                                                   *
# *                                                                         *
# ***************************************************************************

__title__ = "FreeCAD FEM _element fluid 1D document object"
__author__ = "Ofentse Kgoa"
__url__ = "https://www.freecad.org"

## @package element_fluid1D
#  \ingroup FEM
#  \brief element fluid 1D object

from . import base_femelement


class ElementFluid1D(base_femelement.BaseFemElement):
    """
    The element_fluid1D object
    """

    Type = "Fem::ElementFluid1D"

    known_fluid_types = ["Liquid"]
    # "Gas", "Open Channel" are not implemented in ccx writer
    # known_fluid_types = ["Liquid", "Gas", "Open Channel"]
    known_liquid_types = [
        "PIPE MANNING",
        "PIPE ENLARGEMENT",
        "PIPE CONTRACTION",
        "PIPE INLET",
        "PIPE OUTLET",
        "PIPE ENTRANCE",
        "PIPE DIAPHRAGM",
        "PIPE BEND",
        "PIPE GATE VALVE",
        "LIQUID PUMP",
        "PIPE WHITE-COLEBROOK",
    ]
    known_gas_types = ["NONE"]
    known_channel_types = ["NONE"]

    def __init__(self, obj):
        super().__init__(obj)

        obj.addProperty(
            "App::PropertyEnumeration", "SectionType", "FluidSection", "select fluid section type"
        )
        obj.setPropertyStatus("SectionType", "LockDynamic")

        obj.addProperty(
            "App::PropertyEnumeration",
            "LiquidSectionType",
            "LiquidSection",
            "select liquid section type",
        )
        obj.setPropertyStatus("LiquidSectionType", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "ManningArea",
            "LiquidManning",
            "set area of the manning fluid section",
        )
        obj.setPropertyStatus("ManningArea", "LockDynamic")

        obj.addProperty(
            "App::PropertyLength",
            "ManningRadius",
            "LiquidManning",
            "set hydraulic radius of manning fluid section",
        )
        obj.setPropertyStatus("ManningRadius", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloat",
            "ManningCoefficient",
            "LiquidManning",
            "set coefficient of manning fluid section",
        )
        obj.setPropertyStatus("ManningCoefficient", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "EnlargeArea1",
            "LiquidEnlargement",
            "set initial area of the enlargement fluid section",
        )
        obj.setPropertyStatus("EnlargeArea1", "LockDynamic")
        obj.addProperty(
            "App::PropertyArea",
            "EnlargeArea2",
            "LiquidEnlargement",
            "set enlarged area of enlargement fluid section",
        )
        obj.setPropertyStatus("EnlargeArea2", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "ContractArea1",
            "LiquidContraction",
            "set initial area of the contraction fluid section",
        )
        obj.setPropertyStatus("ContractArea1", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "ContractArea2",
            "LiquidContraction",
            "set contracted area of contraction fluid section",
        )
        obj.setPropertyStatus("ContractArea2", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloat",
            "InletPressure",
            "LiquidInlet",
            "set inlet pressure for fluid section",
        )
        obj.setPropertyStatus("InletPressure", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloat",
            "OutletPressure",
            "LiquidOutlet",
            "set outlet pressure for fluid section",
        )
        obj.setPropertyStatus("OutletPressure", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloat",
            "InletFlowRate",
            "LiquidInlet",
            "set inlet mass flow rate for fluid section",
        )
        obj.setPropertyStatus("InletFlowRate", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloat",
            "OutletFlowRate",
            "LiquidOutlet",
            "set outlet mass flow rate for fluid section",
        )
        obj.setPropertyStatus("OutletFlowRate", "LockDynamic")

        obj.addProperty(
            "App::PropertyBool",
            "InletPressureActive",
            "LiquidInlet",
            "activates or deactivates inlet pressure for fluid section",
        )
        obj.setPropertyStatus("InletPressureActive", "LockDynamic")

        obj.addProperty(
            "App::PropertyBool",
            "OutletPressureActive",
            "LiquidOutlet",
            "activates or deactivates outlet pressure for fluid section",
        )
        obj.setPropertyStatus("OutletPressureActive", "LockDynamic")

        obj.addProperty(
            "App::PropertyBool",
            "InletFlowRateActive",
            "LiquidInlet",
            "activates or deactivates inlet flow rate for fluid section",
        )
        obj.setPropertyStatus("InletFlowRateActive", "LockDynamic")

        obj.addProperty(
            "App::PropertyBool",
            "OutletFlowRateActive",
            "LiquidOutlet",
            "activates or deactivates outlet flow rate for fluid section",
        )
        obj.setPropertyStatus("OutletFlowRateActive", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "EntrancePipeArea",
            "LiquidEntrance",
            "set the pipe area of the entrance fluid section",
        )
        obj.setPropertyStatus("EntrancePipeArea", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "EntranceArea",
            "LiquidEntrance",
            "set the entrance area of the entrance fluid section",
        )
        obj.setPropertyStatus("EntranceArea", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "DiaphragmPipeArea",
            "LiquidDiaphragm",
            "set the pipe area of the diaphragm fluid section",
        )
        obj.setPropertyStatus("DiaphragmPipeArea", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "DiaphragmArea",
            "LiquidDiaphragm",
            "set the diaphragm area of the diaphragm fluid section",
        )
        obj.setPropertyStatus("DiaphragmArea", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "BendPipeArea",
            "LiquidBend",
            "set pipe area of the bend fluid section",
        )
        obj.setPropertyStatus("BendPipeArea", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloat",
            "BendRadiusDiameter",
            "LiquidBend",
            "set ratio of bend radius over pipe diameter of the bend fluid section",
        )
        obj.setPropertyStatus("BendRadiusDiameter", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloat",
            "BendAngle",
            "LiquidBend",
            "set bend angle of the bend fluid section",
        )
        obj.setPropertyStatus("BendAngle", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloat",
            "BendLossCoefficient",
            "LiquidBend",
            "set loss coefficient of the bend fluid section",
        )
        obj.setPropertyStatus("BendLossCoefficient", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "GateValvePipeArea",
            "LiquidGateValve",
            "set pipe area of the gate valve fluid section",
        )
        obj.setPropertyStatus("GateValvePipeArea", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloat",
            "GateValveClosingCoeff",
            "LiquidGateValve",
            "set closing coefficient of the gate valve fluid section",
        )
        obj.setPropertyStatus("GateValveClosingCoeff", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloatList",
            "PumpFlowRate",
            "LiquidPump",
            "set the pump characteristic flow rate of the pump fluid section",
        )
        obj.setPropertyStatus("PumpFlowRate", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloatList",
            "PumpHeadLoss",
            "LiquidPump",
            "set the pump characteristic head loss of the pump fluid section",
        )
        obj.setPropertyStatus("PumpHeadLoss", "LockDynamic")

        obj.addProperty(
            "App::PropertyArea",
            "ColebrookeArea",
            "LiquidColebrooke",
            "set pipe area of the colebrooke fluid section",
        )
        obj.setPropertyStatus("ColebrookeArea", "LockDynamic")

        obj.addProperty(
            "App::PropertyLength",
            "ColebrookeRadius",
            "LiquidColebrooke",
            "set hydraulic radius of the colebrooke fluid section",
        )
        obj.setPropertyStatus("ColebrookeRadius", "LockDynamic")

        obj.addProperty(
            "App::PropertyLength",
            "ColebrookeGrainDiameter",
            "LiquidColebrooke",
            "set grain diameter of the colebrooke fluid section",
        )
        obj.setPropertyStatus("ColebrookeGrainDiameter", "LockDynamic")

        obj.addProperty(
            "App::PropertyFloat",
            "ColebrookeFormFactor",
            "LiquidColebrooke",
            "set coefficient of the colebrooke fluid section",
        )
        obj.setPropertyStatus("ColebrookeFormFactor", "LockDynamic")

        obj.addProperty(
            "App::PropertyEnumeration", "GasSectionType", "GasSection", "select gas section type"
        )
        obj.setPropertyStatus("GasSectionType", "LockDynamic")

        obj.addProperty(
            "App::PropertyEnumeration",
            "ChannelSectionType",
            "ChannelSection",
            "select channel section type",
        )
        obj.setPropertyStatus("ChannelSectionType", "LockDynamic")

        # set property default values
        obj.SectionType = ElementFluid1D.known_fluid_types
        obj.SectionType = "Liquid"
        obj.LiquidSectionType = ElementFluid1D.known_liquid_types
        obj.LiquidSectionType = "PIPE INLET"
        obj.GasSectionType = ElementFluid1D.known_gas_types
        obj.GasSectionType = "NONE"
        obj.ChannelSectionType = ElementFluid1D.known_channel_types
        obj.ChannelSectionType = "NONE"
        obj.ManningArea = 10.0
        obj.ManningRadius = 1.0
        obj.ManningCoefficient = 0.0015  # has units of s/mm^(1/3)
        obj.EnlargeArea1 = 10.0
        obj.EnlargeArea2 = 20.0
        obj.ContractArea1 = 20.0
        obj.ContractArea2 = 10.0
        obj.EntrancePipeArea = 20.0
        obj.EntranceArea = 20.0
        obj.DiaphragmPipeArea = 20.0
        obj.DiaphragmArea = 20.0
        obj.BendPipeArea = 20.0
        obj.BendRadiusDiameter = 1.0
        obj.BendAngle = 0.0
        obj.BendLossCoefficient = 0.0
        obj.GateValvePipeArea = 20.0
        obj.GateValveClosingCoeff = 0.125
        obj.PumpFlowRate = [0, 1.04e-04, 2.08e-4, 3.13e-4, 4.17e-4]
        obj.PumpHeadLoss = [30, 29.17, 26.67, 23.33, 18.33]
        obj.ColebrookeArea = 20.0
        obj.ColebrookeRadius = 1.0
        obj.ColebrookeGrainDiameter = 0.0025
        obj.ColebrookeFormFactor = 1.0
        obj.InletPressure = 1.0
        obj.OutletPressure = 1.0
        obj.InletFlowRate = 1.0
        obj.OutletFlowRate = 1.0
        obj.InletPressureActive = True
        obj.OutletPressureActive = True
        obj.InletFlowRateActive = False
        obj.OutletFlowRateActive = False

	# ***************************************************************************
	# * Copyright (c) 2016 Ofentse Kgoa <kgoaot@eskom.co.za> *
	# * Copyright (c) 2018 Bernd Hahnebach <bernd@bimstatik.org> *
	# * Based on the FemElementGeometry1D by Bernd Hahnebach *
	# * *
	# * This file is part of the FreeCAD CAx development system. *
	# * *
	# * This program is free software; you can redistribute it and/or modify *
	# * it under the terms of the GNU Lesser General Public License (LGPL) *
	# * as published by the Free Software Foundation; either version 2 of *
	# * the License, or (at your option) any later version. *
	# * for detail see the LICENCE text file. *
	# * *
	# * This program is distributed in the hope that it will be useful, *
	# * but WITHOUT ANY WARRANTY; without even the implied warranty of *
	# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
	# * GNU Library General Public License for more details. *
	# * *
	# * You should have received a copy of the GNU Library General Public *
	# * License along with this program; if not, write to the Free Software *
	# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 *
	# * USA *
	# * *
	# ***************************************************************************

	__title__ = "FreeCAD FEM _element fluid 1D document object"
	__author__ = "Ofentse Kgoa"
	__url__ = "https://www.freecad.org"

	## @package element_fluid1D
	# \ingroup FEM
	# \brief element fluid 1D object

	from . import base_femelement


	class ElementFluid1D(base_femelement.BaseFemElement):
	"""
	The element_fluid1D object
	"""

	Type = "Fem::ElementFluid1D"

	known_fluid_types = ["Liquid"]
	# "Gas", "Open Channel" are not implemented in ccx writer
	# known_fluid_types = ["Liquid", "Gas", "Open Channel"]
	known_liquid_types = [
	"PIPE MANNING",
	"PIPE ENLARGEMENT",
	"PIPE CONTRACTION",
	"PIPE INLET",
	"PIPE OUTLET",
	"PIPE ENTRANCE",
	"PIPE DIAPHRAGM",
	"PIPE BEND",
	"PIPE GATE VALVE",
	"LIQUID PUMP",
	"PIPE WHITE-COLEBROOK",
	]
	known_gas_types = ["NONE"]
	known_channel_types = ["NONE"]

	def __init__(self, obj):
	super().__init__(obj)

	obj.addProperty(
	"App::PropertyEnumeration", "SectionType", "FluidSection", "select fluid section type"
	)
	obj.setPropertyStatus("SectionType", "LockDynamic")

	obj.addProperty(
	"App::PropertyEnumeration",
	"LiquidSectionType",
	"LiquidSection",
	"select liquid section type",
	)
	obj.setPropertyStatus("LiquidSectionType", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"ManningArea",
	"LiquidManning",
	"set area of the manning fluid section",
	)
	obj.setPropertyStatus("ManningArea", "LockDynamic")

	obj.addProperty(
	"App::PropertyLength",
	"ManningRadius",
	"LiquidManning",
	"set hydraulic radius of manning fluid section",
	)
	obj.setPropertyStatus("ManningRadius", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloat",
	"ManningCoefficient",
	"LiquidManning",
	"set coefficient of manning fluid section",
	)
	obj.setPropertyStatus("ManningCoefficient", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"EnlargeArea1",
	"LiquidEnlargement",
	"set initial area of the enlargement fluid section",
	)
	obj.setPropertyStatus("EnlargeArea1", "LockDynamic")
	obj.addProperty(
	"App::PropertyArea",
	"EnlargeArea2",
	"LiquidEnlargement",
	"set enlarged area of enlargement fluid section",
	)
	obj.setPropertyStatus("EnlargeArea2", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"ContractArea1",
	"LiquidContraction",
	"set initial area of the contraction fluid section",
	)
	obj.setPropertyStatus("ContractArea1", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"ContractArea2",
	"LiquidContraction",
	"set contracted area of contraction fluid section",
	)
	obj.setPropertyStatus("ContractArea2", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloat",
	"InletPressure",
	"LiquidInlet",
	"set inlet pressure for fluid section",
	)
	obj.setPropertyStatus("InletPressure", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloat",
	"OutletPressure",
	"LiquidOutlet",
	"set outlet pressure for fluid section",
	)
	obj.setPropertyStatus("OutletPressure", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloat",
	"InletFlowRate",
	"LiquidInlet",
	"set inlet mass flow rate for fluid section",
	)
	obj.setPropertyStatus("InletFlowRate", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloat",
	"OutletFlowRate",
	"LiquidOutlet",
	"set outlet mass flow rate for fluid section",
	)
	obj.setPropertyStatus("OutletFlowRate", "LockDynamic")

	obj.addProperty(
	"App::PropertyBool",
	"InletPressureActive",
	"LiquidInlet",
	"activates or deactivates inlet pressure for fluid section",
	)
	obj.setPropertyStatus("InletPressureActive", "LockDynamic")

	obj.addProperty(
	"App::PropertyBool",
	"OutletPressureActive",
	"LiquidOutlet",
	"activates or deactivates outlet pressure for fluid section",
	)
	obj.setPropertyStatus("OutletPressureActive", "LockDynamic")

	obj.addProperty(
	"App::PropertyBool",
	"InletFlowRateActive",
	"LiquidInlet",
	"activates or deactivates inlet flow rate for fluid section",
	)
	obj.setPropertyStatus("InletFlowRateActive", "LockDynamic")

	obj.addProperty(
	"App::PropertyBool",
	"OutletFlowRateActive",
	"LiquidOutlet",
	"activates or deactivates outlet flow rate for fluid section",
	)
	obj.setPropertyStatus("OutletFlowRateActive", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"EntrancePipeArea",
	"LiquidEntrance",
	"set the pipe area of the entrance fluid section",
	)
	obj.setPropertyStatus("EntrancePipeArea", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"EntranceArea",
	"LiquidEntrance",
	"set the entrance area of the entrance fluid section",
	)
	obj.setPropertyStatus("EntranceArea", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"DiaphragmPipeArea",
	"LiquidDiaphragm",
	"set the pipe area of the diaphragm fluid section",
	)
	obj.setPropertyStatus("DiaphragmPipeArea", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"DiaphragmArea",
	"LiquidDiaphragm",
	"set the diaphragm area of the diaphragm fluid section",
	)
	obj.setPropertyStatus("DiaphragmArea", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"BendPipeArea",
	"LiquidBend",
	"set pipe area of the bend fluid section",
	)
	obj.setPropertyStatus("BendPipeArea", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloat",
	"BendRadiusDiameter",
	"LiquidBend",
	"set ratio of bend radius over pipe diameter of the bend fluid section",
	)
	obj.setPropertyStatus("BendRadiusDiameter", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloat",
	"BendAngle",
	"LiquidBend",
	"set bend angle of the bend fluid section",
	)
	obj.setPropertyStatus("BendAngle", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloat",
	"BendLossCoefficient",
	"LiquidBend",
	"set loss coefficient of the bend fluid section",
	)
	obj.setPropertyStatus("BendLossCoefficient", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"GateValvePipeArea",
	"LiquidGateValve",
	"set pipe area of the gate valve fluid section",
	)
	obj.setPropertyStatus("GateValvePipeArea", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloat",
	"GateValveClosingCoeff",
	"LiquidGateValve",
	"set closing coefficient of the gate valve fluid section",
	)
	obj.setPropertyStatus("GateValveClosingCoeff", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloatList",
	"PumpFlowRate",
	"LiquidPump",
	"set the pump characteristic flow rate of the pump fluid section",
	)
	obj.setPropertyStatus("PumpFlowRate", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloatList",
	"PumpHeadLoss",
	"LiquidPump",
	"set the pump characteristic head loss of the pump fluid section",
	)
	obj.setPropertyStatus("PumpHeadLoss", "LockDynamic")

	obj.addProperty(
	"App::PropertyArea",
	"ColebrookeArea",
	"LiquidColebrooke",
	"set pipe area of the colebrooke fluid section",
	)
	obj.setPropertyStatus("ColebrookeArea", "LockDynamic")

	obj.addProperty(
	"App::PropertyLength",
	"ColebrookeRadius",
	"LiquidColebrooke",
	"set hydraulic radius of the colebrooke fluid section",
	)
	obj.setPropertyStatus("ColebrookeRadius", "LockDynamic")

	obj.addProperty(
	"App::PropertyLength",
	"ColebrookeGrainDiameter",
	"LiquidColebrooke",
	"set grain diameter of the colebrooke fluid section",
	)
	obj.setPropertyStatus("ColebrookeGrainDiameter", "LockDynamic")

	obj.addProperty(
	"App::PropertyFloat",
	"ColebrookeFormFactor",
	"LiquidColebrooke",
	"set coefficient of the colebrooke fluid section",
	)
	obj.setPropertyStatus("ColebrookeFormFactor", "LockDynamic")

	obj.addProperty(
	"App::PropertyEnumeration", "GasSectionType", "GasSection", "select gas section type"
	)
	obj.setPropertyStatus("GasSectionType", "LockDynamic")

	obj.addProperty(
	"App::PropertyEnumeration",
	"ChannelSectionType",
	"ChannelSection",
	"select channel section type",
	)
	obj.setPropertyStatus("ChannelSectionType", "LockDynamic")

	# set property default values
	obj.SectionType = ElementFluid1D.known_fluid_types
	obj.SectionType = "Liquid"
	obj.LiquidSectionType = ElementFluid1D.known_liquid_types
	obj.LiquidSectionType = "PIPE INLET"
	obj.GasSectionType = ElementFluid1D.known_gas_types
	obj.GasSectionType = "NONE"
	obj.ChannelSectionType = ElementFluid1D.known_channel_types
	obj.ChannelSectionType = "NONE"
	obj.ManningArea = 10.0
	obj.ManningRadius = 1.0
	obj.ManningCoefficient = 0.0015 # has units of s/mm^(1/3)
	obj.EnlargeArea1 = 10.0
	obj.EnlargeArea2 = 20.0
	obj.ContractArea1 = 20.0
	obj.ContractArea2 = 10.0
	obj.EntrancePipeArea = 20.0
	obj.EntranceArea = 20.0
	obj.DiaphragmPipeArea = 20.0
	obj.DiaphragmArea = 20.0
	obj.BendPipeArea = 20.0
	obj.BendRadiusDiameter = 1.0
	obj.BendAngle = 0.0
	obj.BendLossCoefficient = 0.0
	obj.GateValvePipeArea = 20.0
	obj.GateValveClosingCoeff = 0.125
	obj.PumpFlowRate = [0, 1.04e-04, 2.08e-4, 3.13e-4, 4.17e-4]
	obj.PumpHeadLoss = [30, 29.17, 26.67, 23.33, 18.33]
	obj.ColebrookeArea = 20.0
	obj.ColebrookeRadius = 1.0
	obj.ColebrookeGrainDiameter = 0.0025
	obj.ColebrookeFormFactor = 1.0
	obj.InletPressure = 1.0
	obj.OutletPressure = 1.0
	obj.InletFlowRate = 1.0
	obj.OutletFlowRate = 1.0
	obj.InletPressureActive = True
	obj.OutletPressureActive = True
	obj.InletFlowRateActive = False
	obj.OutletFlowRateActive = False