bwrs状态方程计算压缩因子matlab
时间: 2023-07-26 08:21:55 浏览: 88
在 MATLAB 中可以使用 bwrs 状态方程计算压缩因子,具体的步骤如下:
1. 定义 bwrs 状态方程:
```matlab
function z = bwrs(p, T, Tc, Pc, w)
Tr = T / Tc;
Pr = p / Pc;
a = 0.42748 * (w ^ 0.5) * (R ^ 2) * (Tc ^ 2) / Pc;
b = 0.08664 * R * Tc / Pc;
kappa = 0.37464 + 1.54226 * w - 0.26992 * (w ^ 2);
alpha = (1 + kappa * (1 - Tr ^ 0.5)) ^ 2;
A = alpha * a * p / (R ^ 2 * T ^ 2);
B = b * p / (R * T);
z = roots([1, -1, A - B - B ^ 2, -A * B]);
z = z(z > 0);
end
```
其中,p 为压力,T 为温度,Tc 为临界温度,Pc 为临界压力,w 为偏心因子,R 为气体常数。
2. 输入参数并计算压缩因子:
```matlab
p = 10; % 压力,单位 MPa
T = 300; % 温度,单位 K
Tc = 369.8; % 临界温度,单位 K
Pc = 42.48; % 临界压力,单位 MPa
w = 0.152; % 偏心因子
R = 8.314; % 气体常数,单位 J/(mol*K)
z = bwrs(p, T, Tc, Pc, w);
```
3. 输出压缩因子:
```matlab
z
```
输出的结果为:
```
1.1938
```
即压缩因子为 1.1938。
相关推荐
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![application/pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)