Scribd Logo
Upload

Welcome to Scribd!

  • 2 views

    Example Cycle Printout

    Uploaded by

    Arnob Das

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    Example Cycle Printout

    Uploaded by

    Arnob Das
    2 views4 pages

    Original Title

    example_cycle_printout

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    2 views4 pages

    Example Cycle Printout

    Uploaded by

    Arnob Das

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 4

    File:example_cycle.

    EES 8/20/2023 4:57:24 PM Page 1


    EES Ver. 10.835: #1867: For use by students and faculty, College of Engineering, University of Oklahoma, Stillwater, OK

    1: {MAE 5353 - Simple Example Cycle


    2: 2022-09-13 - Christian Bach}
    3: {! You may want to update the header and description for your current problem}
    4: "Setup"
    5: $UnitSystem SI MASS DEG bar C kJ
    6: "Conditionals"
    7: $IFNOT PARAMETRICTABLE
    8: T_evap=0 [C]
    9: T_cond = 30 [C] "saturation temperature in condenser"
    10: $ENDIF
    11:
    12: "!Thermodynamic performance of a vapor compression refrigeration cycle."
    13:
    14: "Cycle Specs"
    15: R$ = 'R134a'
    16: DELTAT_superheat=5 [C] "superheat at compressor inlet"
    17: DELTAT_subcool=5 [C] "subcooling at condenser outlet"
    18: T[1] = T_evap+DELTAT_superheat,
    19: T[3] = T_cond-DELTAT_subcool
    20:
    21: "State 1: Entering the compressor"
    22: P[1] = pressure(R$,T=T_evap, x=1)
    23: h[1] = enthalpy(R$,T=T[1], P=P[1])
    24: s[1] = entropy(R$,T=T[1], P=P[1])
    25: v[1]=volume(R$,T=T[1], P=P[1])
    26:
    27: "State 2: Entering the condenser (isentropic compression)"
    28: P[2] = P[3] "neglect pressure drop in condenser"
    29: s_id[2] = s[1] "ideal compression process is isentropic"
    30: h_id[2] = enthalpy(R$,P=P[2], s=s_id[2]) "enthalpy at outlet state assuming isentropic process"
    31: eta=(h_id[2]-h[1])/(h[2]-h[1]) "definition of compressor efficiency"
    32: eta=0.7 "compression isentropic efficiency"
    33: s[2]=entropy(R$,P=P[2], h=h[2]) "entropy at actual compressor outlet state"
    34:
    35: "State 3: Entering the expansion device"
    36: P[3] = pressure(R$,T=T_cond,x=0) "saturation pressure in condenser"
    37: h[3] = enthalpy(R$,T=T[3], P=P[3])
    38: s[3] = entropy(R$,T=T[3], P=P[3])
    39:
    40: "State 4: Entering the evaporator from a expansion valve"
    41: P[4] = P[1]
    42: h[4] = h[3] "isenthalpic expansion in throttle"
    43: s[4]=entropy(R$, h=h[4], P=P[4])
    44: T[4]=temperature(R$, h=h[4], P=P[4])
    45:
    46: "State 5: Same as state 1 for closing the cycle when plotting"
    47: P[5] = P[1]
    48: h[5] = h[1]

    Setup

    Conditionals

    T evap = 0 [C]

    T cond = 30 [C] saturation temperature in condenser

    Thermodynamic performance of a vapor compression refrigeration cycle.


    File:example_cycle.EES 8/20/2023 4:57:24 PM Page 2
    EES Ver. 10.835: #1867: For use by students and faculty, College of Engineering, University of Oklahoma, Stillwater, OK

    Cycle Specs

    R$ = 'R134a'

    DT superheat = 5 [C] superheat at compressor inlet

    DT subcool = 5 [C] subcooling at condenser outlet

    T 1 = T evap + DT superheat

    T 3 = T cond – DT subcool

    State 1: Entering the compressor

    P 1 = P R$ , T = T evap , x = 1

    h 1 = h R$ , T = T 1 , P = P 1

    s 1 = s R$ , T = T 1 , P = P 1

    v 1 = v R$ , T = T 1 , P = P 1

    State 2: Entering the condenser (isentropic compression)

    P 2 = P 3 neglect pressure drop in condenser

    s id,2 = s 1 ideal compression process is isentropic

    h id,2 = h R$ , P = P 2 , s = s id,2 enthalpy at outlet state assuming isentropic process

    h id,2 – h 1
    h = definition of compressor efficiency
    h2 – h1

    h = 0.7 compression isentropic efficiency

    s 2 = s R$ , P = P 2 , h = h 2 entropy at actual compressor outlet state

    State 3: Entering the expansion device

    P 3 = P R$ , T = T cond , x = 0 saturation pressure in condenser

    h 3 = h R$ , T = T 3 , P = P 3

    s 3 = s R$ , T = T 3 , P = P 3

    State 4: Entering the evaporator from a expansion valve

    P4 = P1

    h 4 = h 3 isenthalpic expansion in throttle

    s 4 = s R$ , h = h 4 , P = P 4

    T 4 = T R$ , h = h 4 , P = P 4
    File:example_cycle.EES 8/20/2023 4:57:24 PM Page 3
    EES Ver. 10.835: #1867: For use by students and faculty, College of Engineering, University of Oklahoma, Stillwater, OK

    State 5: Same as state 1 for closing the cycle when plotting

    P5 = P1

    h5 = h1

    SOLUTION
    Unit Settings: SI C bar kJ mass deg
    DTsubcool = 5 [C] DTsuperheat = 5 [C] h = 0.7 [-]
    R$ = 'R134a' Tcond = 30 [C] Tevap = 0 [C]

    No unit problems were detected.

    Arrays Table: Main


    Ti Pi hi hid,i si sid,i vi
    [C] [bar] [kJ/kg] [kJ/kg] [kJ/kg-K] [kJ/kg-K] [m 3/kg]

    1 5 2.93 254.9 0.9476 0.07104


    2 7.706 284.3 275.5 0.9755 0.9476
    3 25 7.706 86.4 0.324
    4 -7.437E-10 2.93 86.4 0.3309
    5 2.93 254.9
    File:example_cycle.EES 8/20/2023 4:57:24 PM Page 4
    EES Ver. 10.835: #1867: For use by students and faculty, College of Engineering, University of Oklahoma, Stillwater, OK

    5x102 R134a

    0.5

    0.6

    8
    0.

    0.

    K
    0.9

    g-
    J/k
    102

    1k
    110°C 160°C 230°C
    59°C
    P [bar]

    101
    3 2
    16°C

    4 1 P[i]

    -21°C
    100

    0.2 0.4 0.6 0.8


    2x10-1
    -100 0 100 200 300 400 500
    h [kJ/kg]

    You might also like