HYTOS

# coding: utf-8

"""

    This is calculation module

"""

import sys

import os

from AppDocData import *

import math

class Calculation:

    def __init__(self, hmb):

        self.initUnits()

        self._hmb = hmb

        if self._hmb.phase_type == 'Vapor':

            self.calculation_Vapor()

        elif self._hmb.phase_type == 'Liquid':

            self.calculation_Liquid()

        elif self._hmb.phase_type == 'Mixed':

            self.calculation_Mixed()

    def initUnits(self):

        activeDrawing = AppDocData.instance().activeDrawing

        for attr in activeDrawing.attrs:

            if attr[0] == 'Units':

                self.flowrate_mass_unit = attr[1]['Flowrate_Mass']                

                self.flowrate_volume_unit = attr[1]['Flowrate_Volume']

                self.density_unit = attr[1]['Density']

                self.viscosity_unit = attr[1]['Viscosity']                

                self.velocity_unit = attr[1]['Velocity']

                self.pressure_unit = attr[1]['Pressure']

                self.length_unit = attr[1]['Length']

                self.pipe_diameter_unit = attr[1]['Pipe_Diameter']                

                self.roughness_unit = attr[1]['Roughness']                

                self.temperature_unit = attr[1]['Temperature']

    def getBarometricPressure(self):

        barometric_pressure = None

        if self.pressure_unit == 'kg/cm2':

            barometric_pressure = 1.033

        elif self.pressure_unit == 'bar':

            barometric_pressure = 1.01325

        elif self.pressure_unit == 'psi':

            barometric_pressure = 14.7

        elif self.pressure_unit == 'mmHg':

            barometric_pressure = 760

        elif self.pressure_unit == 'kPa':

            barometric_pressure = 101.325

        elif self.pressure_unit == 'MPa':

            barometric_pressure = 0.101325

        return barometric_pressure

    def calculation_Vapor(self):  

        # (1) density 입력

        Ref_baro = self.getBarometricPressure()

        press2 = 100 # To쪽 Equipment의 Pressure

        if self.pressure_unit == 'kg/cm2':

            press2 = press2 + Ref_baro

        elif self.pressure_unit == 'psi':

            press2 = press2 / 14.7 * 1.033 + Ref_baro

        elif self.pressure_unit == 'atm':

            press2 = press2 * 1.033 + Ref_baro

        elif self.pressure_unit == 'bar':

            press2 = press2 / 1.013 * 1.033 + Ref_baro

        elif self.pressure_unit == 'mmHg':

            press2 = press2 / 760 * 1.033 + Ref_baro

        elif self.pressure_unit == 'kPa':

            press2 = press2 / 101.325 * 1.033 + Ref_baro

        elif self.pressure_unit == 'MPa':

            press2 = press2 / 0.101325 * 1.033 + Ref_baro

        temp = self._hmb.temperature

        # temp를 kalvin으로 맞추기

        if self.temperature_unit == '℃':

            temp = temp + 273.15

        elif self.temperature_unit == '℉':

            temp = (temp - 32) / 1.8 + 273.15

        # 'MW 가져오기

        mw = self._hmb.molecular_weight

        # 'Z 가져오기

        z = self._hmb.compress_factor

        # 밀도 계산

        density2 = press2 * mw / 0.08206 / temp / z / 1.033

        # '밀도 입력

        if self.density_unit == 'kg/m3':

            density2 = density2

        elif self.density_unit == 'lb/ft3':

            density2 = density2 * 0.062428

        # (2) static P를 계산

        # 1. density 받음

        if self.density_unit == 'kg/m3':

            density2 = density2

        elif self.density_unit == 'lb/ft3':

            density2 = density2 * 16.0185

        # 2. elevation 받음

        el1 = 200 # From쪽 Equipment의 Elevation

        el2 = 100 # To쪽 Equipment의 Elevation

        if self.length_unit == 'm':

            el1 = el1

            el2 = el2

        elif self.length_unit == "in":

            el1 = el1 * 0.0254

            el2 = el2 * 0.0254

        elif self.length_unit == "ft":

            el1 = el1 * 0.3048

            el2 = el2 * 0.3048

        elif self.length_unit == "yd":

            el1 = el1 * 0.9144

            el2 = el2 * 0.9144

        elif self.length_unit == "mm":

            el1 = el1 * 0.001

            el2 = el2 * 0.001

        elif self.length_unit == 'mile':

            el1 = el1 * 1609.34

            el2 = el2 * 1609.34

        # 3. static head 계산 

        # 'atm으로 계산된 dp

        stat_dp = (el2 - el1) * density2 / 1000 * 9.80665 / 101.325 

        # 4. 압력 유닛에 맞춰 뿌리기

        if self.pressure_unit == 'kg/cm2':

            stat_dp = stat_dp * 1.033

        elif self.pressure_unit == 'psi':

            stat_dp = stat_dp * 14.7

        elif self.pressure_unit == 'atm':

            stat_dp = stat_dp

        elif self.pressure_unit == 'bar':

            stat_dp = stat_dp * 1.013

        elif self.pressure_unit == 'mmHg':

            stat_dp = stat_dp * 1.013

        elif self.pressure_unit == 'kPa':

            stat_dp = stat_dp * 101.325

        elif self.pressure_unit == 'MPa':

            stat_dp = stat_dp * 0.101325

        # '(3) pressure drop 계산

        # 'vapor인 경우 압력강하를 계산해주고 넣어주는 모듈

        # '입력된 vapor의 압력과 밀도를 가지고 끊어서 계산하는 factor를 적용했을때 length가 얼마나 나오는지 판별하는 것.

        ida = self._hmb.inside_pipe_size

        if self.pipe_diameter_unit == 'in':

            ida = ida * 0.0254

        elif self.pipe_diameter_unit == 'mm':

            ida = ida / 1000

        # 'mass 가져오기 기준 kg/h

        mass = self._hmb.flowrate_mass

        if self.flowrate_mass_unit == 'kg/h':

            mass = mass

        elif self.flowrate_mass_unit == 'g/min':

            mass = mass * 60 / 1000

        elif self.flowrate_mass_unit == 'lb/h':

            mass = mass * 0.453592

        elif self.flowrate_mass_unit == 't/h':

            mass = mass * 1000

        # 'g 구하기 (kg/m2/s)

        g = mass / 3600 / (3.1415 * ida ** 2 / 4)

        # '속도 계산 (m/s)

        velocity = 4 * mass / density2 / 3.1415 / ida ** 2 / 3600

        # 'k 가져오기

        k = self._hmb.specific_heat_ratio

        # 'Mach Number 계산

        mach = velocity / ((k * 9.80665 * 847.28 * temp / mw) ** 0.5)

        # '부피유량 계산

        volume = mass / density2

        # '현재 volume은 m3/h임. 부피유량 단위에 맞춰 뿌려줌

        if self.flowrate_volume_unit == 'm3/h':

            volume = round(volume, 3)    

        elif self.flowrate_volume_unit == 'l/min':

            volume = round(volume / 60 * 1000, 3)

        elif self.flowrate_volume_unit == 'ft3/h':

            volume = round(volume * 35.3147, 3)

        elif self.flowrate_volume_unit == 'USgpm':

            volume = round(volume * 4.40287, 3)

        elif self.flowrate_volume_unit == 'BPSD':

            volume = round(volume * 150.955, 3)

        viscosity = self._hmb.viscosity

        # ' viscosity 유닛 변환 (모두 kg/m.s로 바꿀것임)

        if self.viscosity_unit == 'kg/m.sec':

            viscosity = viscosity

        elif self.viscosity_unit == 'cP':

            viscosity = viscosity * 0.001

        elif self.viscosity_unit == 'kg/m.h':

            viscosity = viscosity / 3600

        elif self.viscosity_unit == 'lb/ft.sec':

            viscosity = viscosity * 1.48816

        # 'density case에 따라 re계산

        if self.density_unit == 'kg/m3':

            reynolds = ida * velocity * density2 / viscosity

        elif self.density_unit == 'lb/ft3':

            reynolds = ida * velocity * (density2 * 16.0185) / viscosity

        rough = self._hmb.roughness

        # 'roughness 를 m로 바꿔줌

        if self.roughness_unit == 'm':

            rough = rough

        elif self.roughness_unit == 'ft':

            rough = rough * 0.3048

        elif self.roughness_unit == 'in':

            rough = rough * 0.0254

        elif self.roughness_unit == 'mm':

            rough = rough * 0.001

        # ' reynolds수에 따라 Fanning/Chen friction factor 계산

        if reynolds <= 2100:

            f = 4 * 16 / reynolds

        else:

            af = math.log(rough / ida / 3.7 + (6.7 / reynolds) ** 0.9) / math.log(10)

            f = (-2 * (math.log(rough / 3.7 / ida - 5.02 / reynolds * af) / math.log(10))) ** (-2)

        # '속도와 마하 수 입력

        if self.velocity_unit == 'm/s':

            self._hmb.velocity = round(velocity, 3)

        elif self.velocity_unit == 'ft/s':

            self._hmb.velocity = round(velocity * 3.28084, 3)

        self._hmb.reynolds = round(mach, 5)

        press1est = press2 / 0.95 #HY_Calc_sht.Cells(2, 1)

        density1est = press1est * mw / 0.08206 / temp / z / 1.033

        # ' 라인 길이 도출 (m)

        estlength = (((press1est ** 2 - press2 ** 2) * mw / g ** 2 / 0.08206 / 1.033 ** 2 / temp / z * 101325) - (2 * math.log(density1est / density2))) * ida / f

        # '라인 길이를 유닛에 맞춰서 변환

        if self.length_unit == 'm':

            estlength = estlength

        elif self.length_unit == 'in':

            estlength = estlength * 39.3701

        elif self.length_unit == 'ft':

            estlength = estlength * 3.28084

        elif self.length_unit == 'yd':

            estlength = estlength * 1.09361

        elif self.length_unit == 'mile':

            estlength = estlength * 0.000621371

        elif self.length_unit == 'mm':

            estlength = estlength * 1000

        if estlength > self._hmb.equivalent_length:

            # '주어진 길이를 가지고 압력 강하를 계산하는 모듈이 들어가야함 (safe)

            # ' length 여유 있음. 한번에 isothermal 식으로 계산

            # ' 끊은 결과 length가 여유가 있음. 그래서 단순 계산하는 모듈

            length = self._hmb.equivalent_length

            # 'length 를 m로 변환

            if self.length_unit == 'm':

                length = length

            elif self.length_unit == 'in':

                length = length * 0.0254

            elif self.length_unit == 'ft':

                length = length * 0.3048

            elif self.length_unit == 'yd':

                length = length * 0.9144

            elif self.length_unit == 'mile':

                length = length * 1609.34

            elif self.length_unit == 'mm':

                length = length * 0.001

            # 'Newton's Rule에 따라

            press1 = self.vap_discharge_1(press2, g, mw, temp, f, z, ida, length)

            line_drop_press = 0

            if self.pressure_unit == 'kg/cm2':

                line_drop_press = press1 - press2

            elif self.pressure_unit == 'psi':

                line_drop_press = (press1 - press2) / 1.033 * 14.7

            elif self.pressure_unit == 'atm':

                line_drop_press = (press1 - press2) / 1.033

            elif self.pressure_unit == 'bar':

                line_drop_press = (press1 - press2) / 1.033 * 1.013

            elif self.pressure_unit == 'mmHg':

                line_drop_press = (press1 - press2) / 1.033 * 760

            elif self.pressure_unit == 'kPa':

                line_drop_press = (press1 - press2) / 1.033 * 101.325

            elif self.pressure_unit == 'MPa':

                line_drop_press = (press1 - press2) / 1.033 * 0.101325

            # 'Hole Calculation

            # 'cv와 pump가 다 없을때 (Hole Case)

        else:

            # '끊어서 계산하는 모듈이 들어가야함

            # 'length 여유 없음. 압력 재계산->밀도 재계산->re, f 재계산->압력 재계산 체계로 가야함

            dev = 1

    def vap_discharge_1(self, p2, g, mw, temp, f, z, id, L_eq):

        R = 0.08206

        Tmp_Const = 101325 * mw / (1.033 ** 2 * (g ** 2) * R * temp * z)

        Tmp_Const2 = id / f

        if self.pressure_unit == 'kg/cm2':

            const_iteration = 5

        elif self.pressure_unit == 'psi':

            const_iteration = 75

        elif self.pressure_unit == 'atm':

            const_iteration = 5

        elif self.pressure_unit == 'bar':

            const_iteration = 5

        elif self.pressure_unit == 'mmHg':

            const_iteration = 3800

        elif self.pressure_unit == 'kPa':

            const_iteration = 506

        elif self.pressure_unit == 'MPa':

            const_iteration = 0.5

        x2 = p2 + const_iteration

        delta_x_past = 0    

        loop = True

        while loop:

            X = x2        

            f_x = (Tmp_Const2 * (Tmp_Const * (X ** 2 - p2 ** 2) - 2 * math.log(X / p2))) - L_eq        

            df_x = 2 * Tmp_Const2 * (Tmp_Const * X - 1 / X)        

            x2 = X - (f_x / df_x)    

            delta_x = x2 - X        

            if abs(delta_x_past) == abs(delta_x):

                break

            delta_x_past = delta_x                            

            if delta_x == 0 or x2 < 0:

                loop = False

        '''

        Do

            X = x2        

            f_x = (Tmp_Const2 * (Tmp_Const * (X ^ 2 - p2 ^ 2) - 2 * Log(X / p2))) - L_eq        

            df_x = 2 * Tmp_Const2 * (Tmp_Const * X - 1 / X)        

            x2 = X - (f_x / df_x)    

            delta_x = x2 - X        

            If Abs(delta_x_past) = Abs(delta_x) Then Exit Do        

            delta_x_past = delta_x                            

        Loop Until delta_x = 0 Or x2 < 0

        '''

        return x2

    def calculation_Liquid(self):

        self._hmb.velocity = 2.889

    def calculation_Mixed(self):

        self._hmb.velocity = 2.889

        self.actionCalculation.triggered.connect(self.calculation)

    def calculation(self):

        from AppDocData import AppDocData

        from Calculation import Calculation

        try:

            appDocData = AppDocData.instance()

            if appDocData.imgName is None or appDocData.activeDrawing is None:

                self.showImageSelectionMessageBox()

                return

            hmbs = appDocData.activeDrawing.hmbTable._hmbs

            if hmbs is not None:            

                try:

                    self.progress = QProgressDialog(self.tr("Please wait for a while"), self.tr("Cancel"), 0, 100, self) if not hasattr(self, 'progress') else self.progress

                    self.progress.setWindowModality(Qt.WindowModal)

                    self.progress.setAutoReset(True)

                    self.progress.setAutoClose(True)

                    self.progress.setMinimum(0)

                    self.progress.resize(600,100)

                    self.progress.setWindowTitle(self.tr("Calculate data..."))

                    self.progress.show()

                    maxValue = len(hmbs)

                    self.progress.setMaximum(maxValue)

                    for hmb in hmbs:   

                        self.progress.setValue(self.progress.value() + 1)

                        if hmb.phase_type:

                            Calculation(hmb)                

                        QApplication.processEvents()

                    self.load_HMB()        

                finally:

                    self.progress.setValue(self.progress.maximum())

                    self.progress.hide()

        except Exception as ex:

            message = 'error occured({}) in {}:{}'.format(ex, sys.exc_info()[-1].tb_frame.f_code.co_filename, sys.exc_info()[-1].tb_lineno)

            self.addMessage.emit(MessageType.Error, message)