POLYMATH Report | NLE |

Explicit Equations | 2016-May-31 |

Calculated values of explicit variables

Variable | Value | |

1 | V | 0.065 |

2 | Pc | 72.9 |

3 | R | 0.08206 |

4 | T | 300 |

5 | omega | 0.225 |

6 | Tc | 304.2 |

7 | a | 3.65392 |

8 | b | 0.029668 |

9 | m | 0.826312 |

10 | alpha | 1.01148 |

11 | P | 96.129 |

12 | Tr | 0.986193 |

13 | Pr | 1.31864 |

14 | z | 0.253814 |

Problem source text

# S. 15* - Complex consecutive calculations
# Gas pressure from SRK equation of state # Verified Solution: P = 96.129, z = 0.253814 #Ref.: Prob. 2.9 in Problem Solving in Chemicalâ€¦ P = R*T/(V-b)-alpha*a/(V*(V+b)) # Pressure (atm) a = 0.42747*R^2*Tc^2/Pc #SRK constant b = 0.08664*R*Tc/Pc #SRK constant alpha = (1+m*(1- sqrt(T/Tc)))^2 #SRK constant alpha
m = 0.48508+1.55171*omega-0.1561*omega^2 #SRK constant m Tc = 304.2 #Critical temperature (K) Pc = 72.9 #Critical pressure (atm) R = 0.08206 #Gas constant (L-atm/gmol-K) omega = 0.225 #Acentric factor T = 300 #Temperature (K) V=0.065 # Volume (L/g-mol) #P = 100 #Pressure (atm) Tr=T/Tc Pr=P/Pc z=P*V/(R*T) |

Matlab formatted problem

Create m file called PolyExplicit.m and paste the following text into it.

% S. 15* - Complex consecutive calculations
% Gas pressure from SRK equation of state % Verified Solution: P = 96.129, z = 0.253814 % Ref.: Prob. 2.9 in Problem Solving in Chemicalâ€¦ %Volume (L/g-mol) V = 0.065; %Critical pressure (atm) Pc = 72.9; %Gas constant (L-atm/gmol-K) R = 0.08206; %Temperature (K) T = 300; %Acentric factor omega = 0.225; %Critical temperature (K) Tc = 304.2; %SRK constant a = 0.42747 * R ^ 2 * Tc ^ 2 / Pc; %SRK constant b = 0.08664 * R * Tc / Pc; %SRK constant m m = 0.48508 + 1.55171 * omega - (0.1561 * omega ^ 2); %SRK constant alpha alpha = (1 + m * (1 - sqrt(T / Tc))) ^ 2; %Pressure (atm) P = R * T / (V - b) - (alpha * a / (V * (V + b))); Tr = T / Tc; Pr = P / Pc; z = P * V / (R * T); |

General Settings

number of explicit equations: | 14 |

reporting digits | 8 |