Commit 98000d2e authored by Salko Jr, Robert's avatar Salko Jr, Robert
Browse files

Add ability to set htc under-relaxation

Description:
Also move the entrainment model setter into an already existing setModel method.

Gitlab Ticket # - XXXX
parent 7032c04b
Loading
Loading
Loading
Loading
+12 −14
Original line number Diff line number Diff line
@@ -179,6 +179,9 @@ class Model:
        # Single-phase turbulent mixing coefficient
        me.beta = 0.037

        # Under-relaxation coefficients
        me.betaHTC = None

    @property
    def is_nodal(self):
        """ Access private is_nodal attr """
@@ -202,15 +205,23 @@ class Model:
        """
        me.title = str(title)

    def setModelOptions(me, beta=None):
    def setModelOptions(me, beta=None, enableEntrainment=None, betaHTC=None):
        """ Set various modeling options

        Any model that can be changed via this method has a default.
        This only needs to be called if the user wishes to change a deafult modeling option.

        Args:
           beta (float): Single-phase turbulent mixing coefficietn
           enableEntrainment (bool): Set to True to enable entrainment/de-entrainment of droplets
           betaHTC (float): Under-relaxation coefficient for fluid/solid heat transfer coefficient
        """

        if beta is not None:
            me.beta = beta
        if enableEntrainment is not None:
            me.enableEntrainment = enableEntrainment
        me.betaHTC = betaHTC

    def setFrictionModel(me, model):
        """ Sets the friction model option
@@ -230,19 +241,6 @@ class Model:
            assert 'C' in model['options']
        me.frictionModel = model

    def setModel(me, enableEntrainment=None):
        """ Use to change modeling options.

        Any model that can be changed via this method has a default.
        This only needs to be called if the user wishes to change a deafult modeling option.

        Args:
           enableEntrainment (bool): Set to True to enable entrainment/de-entrainment of droplets

        """
        if enableEntrainment is not None:
            me.enableEntrainment = enableEntrainment

    def setFluidProperties(me, fluidprops):
        """Set the fluid property table to use.

+2 −0
Original line number Diff line number Diff line
@@ -150,6 +150,8 @@ def writeDeck(model, filename):
        group1Data.append("**         A           B           C\n")
        group1Data.append("{:12.4e}{:12.4e}{:12.4e}\n".format(
            model.frictionModel['options']['A'], model.frictionModel['options']['B'], model.frictionModel['options']['C']))
    if model.betaHTC is not None:
        group1Data.append("{{beta_htc}} {:15.3e}\n".format(model.betaHTC))

    # The Card Group 2 data to be written to the deck
    group2Data = []
+1 −0
Original line number Diff line number Diff line
@@ -25,6 +25,7 @@ ctf model
*Card 1.4
**GTP(1)   VFRAC(3)  GTP(2) VFRAC(4)  GTP(3) VFRAC(5)  GTP(4) VFRAC(6)
     air     0.0001
{beta_htc}       5.000e-02
***********************************************************************************************
*GROUP 2 - Channel Description
***********************************************************************************************
+1 −1
Original line number Diff line number Diff line
@@ -61,7 +61,7 @@ def main():
                         Tcool_linf_abs=1e-3, Tsolid_l2_abs=1e2, Tsolid_linf_abs=1e99, vl_l2_abs=1e2, vl_linf_abs=1e-3,
                         vv_l2_abs=1e2, vv_linf_abs=1e-3, vd_l2_abs=1e2, vd_linf_abs=1e-3,
                         dtmin=1e-7, dtmax=1e-4)
    model.setModel(enableEntrainment=True)
    model.setModelOptions(enableEntrainment=True, betaHTC=0.05)
    model.setConductorModelOptions(nc=0)

    model.generateModel()