Context Navigation

sqpipe_domain.py

Timestamp:

Dec 8, 2011, 11:38:18 AM (12 years ago)

Author:

paul

Message:

Fixed loss of mass in sqpipe and rearranged code a little

File:

-                      r8265
+                      r8274
         """
+        # Equivalent area A is defined by A = (\rho/\rho_0) A_{pipe}
+        # A_{pipe} = width * top
+        # A = width * height
+        # \rho_0 = 1
+        # Therefore \rho = height/top for pressurised states
+        # For free surface \rho = 1
+        # Coarse approximation of mass M in a cell, assume rho is constant
+        # M = \Delta x * \rho * A
+        #   = \Delta_x * A * height/top (pressurised)
+        #   = \Delta_x * A (un-pressurised)
+        area      = self.quantities['area']
+        height    = self.quantities['height']
+        top       = self.quantities['top']
+        #Arrays
+        a   = area.centroid_values
+        h   = height.centroid_values
+        t   = top.centroid_values
+        M = numpy.where(self.state == 0, self.areas * a, self.areas * a * h/t)
+        return M
+        # The equivalent area is the conserved mass quantity
+        # It is equal to \rho/\rho_0 A where A is the cross sectional
+        # area of the fluid where $\rho_0 is some fixed reference density
+        area = self.quantities['area']
+        a = area.centroid_values
+        return a * self.areas
 # Auxillary methods

Note: See TracChangeset for help on using the changeset viewer.