Using the Soave-Redlich-Kwong (SRK), or any other equation of state, to determine a specific volume requires...

   

Transcribed Image Text:

Using the Soave-Redlich-Kwong (SRK), or any other equation of state, to determine a specific volume requires a trial and error calculation. You will use Matlab/Octave or Python to do this. The goal of this problem is to determine v (m/kg) for CO at: 200 K and 6.8 atm, 300K and 6.8 atm, 300 K and 21.5 atm, 300 K and 50.0 atm. The SRK equation is as follows: P = @= b = Tr M RT v-b ax v(v + b) 0.427RT2 Pe 0.08664RT Pe (1 + ((0.48508 +1.55171 -0.15613)(1-T05)) T Te (1) (2) (3) (4) (5) Where R is the gas constant,Pe, Te are the critical pressure and temperature, v is the specific volume and w is the acentric factor. 1. Starting with equation 1, derive the following expression for the SRK equation: f(v)= P- RTv + (aa-bP-bRT) v-aab=0 2. Write a MATLAB or Python program to take as inputs a species identifier (e.g CO), the critical temperature, critical pressure, Pitzer acentric factor, and the temperatures and pressures at which v is to be calculated and Use a root finding procedure in Matlab/Octave to determine v for each of the conditions. Using the Soave-Redlich-Kwong (SRK), or any other equation of state, to determine a specific volume requires a trial and error calculation. You will use Matlab/Octave or Python to do this. The goal of this problem is to determine v (m/kg) for CO at: 200 K and 6.8 atm, 300K and 6.8 atm, 300 K and 21.5 atm, 300 K and 50.0 atm. The SRK equation is as follows: P = @= b = Tr M RT v-b ax v(v + b) 0.427RT2 Pe 0.08664RT Pe (1 + ((0.48508 +1.55171 -0.15613)(1-T05)) T Te (1) (2) (3) (4) (5) Where R is the gas constant,Pe, Te are the critical pressure and temperature, v is the specific volume and w is the acentric factor. 1. Starting with equation 1, derive the following expression for the SRK equation: f(v)= P- RTv + (aa-bP-bRT) v-aab=0 2. Write a MATLAB or Python program to take as inputs a species identifier (e.g CO), the critical temperature, critical pressure, Pitzer acentric factor, and the temperatures and pressures at which v is to be calculated and Use a root finding procedure in Matlab/Octave to determine v for each of the conditions.