1/**
2 * 计算基于温度和压力的液态乙烯密度。
3 * 
4 * @param {number} temperatureK - 温度（开尔文，有效范围：104K到282K）
5 * @param {number} pressureBar - 压力（bar，有效范围：1到100 bar）
6 * @returns {number} 液态乙烯的密度（kg/m³）
7 */
8function calculateEthyleneDensity(temperatureK, pressureBar) {
9    // 乙烯的常数
10    const A = 700;
11    const Tc = 283.18;  // 临界温度（K）
12    const n = 0.29683;
13    const B = 0.8;
14    const kappa = 0.00125;  // 等温压缩系数（MPa⁻¹）
15    const P_ref = 0.1;  // 参考压力（MPa，1 bar）
16    
17    // 将压力从bar转换为MPa
18    const pressureMPa = pressureBar / 10;
19    
20    // 计算参考压力下的密度
21    const rhoRef = A * Math.pow(1 - temperatureK/Tc, n) - B * temperatureK;
22    
23    // 应用压力修正
24    const rho = rhoRef * (1 + kappa * (pressureMPa - P_ref));
25    
26    return rho;
27}
28
29// 示例用法
30const temp = 200;  // K
31const pressure = 50;  // bar
32const density = calculateEthyleneDensity(temp, pressure);
33console.log(`液态乙烯在${temp}K和${pressure} bar下的密度: ${density.toFixed(2)} kg/m³`);
34

1' Excel VBA 函数用于液态乙烯密度计算
2Function EthyleneDensity(TemperatureK As Double, PressureBar As Double) As Double
3    ' 乙烯的常数
4    Dim A As Double: A = 700
5    Dim Tc As Double: Tc = 283.18 ' 临界温度（K）
6    Dim n As Double: n = 0.29683
7    Dim B As Double: B = 0.8
8    Dim kappa As Double: kappa = 0.00125 ' 等温压缩系数（MPa⁻¹）
9    Dim P_ref As Double: P_ref = 0.1 ' 参考压力（MPa，1 bar）
10    
11    ' 将压力从bar转换为MPa
12    Dim PressureMPa As Double: PressureMPa = PressureBar / 10
13    
14    ' 计算参考压力下的密度
15    Dim rho_ref As Double: rho_ref = A * (1 - TemperatureK / Tc) ^ n - B * TemperatureK
16    
17    ' 应用压力修正
18    EthyleneDensity = rho_ref * (1 + kappa * (PressureMPa - P_ref))
19End Function
20
21' 在Excel单元格中的用法：
22' =EthyleneDensity(200, 50)
23

1function density = ethyleneDensity(temperatureK, pressureBar)
2    % 计算基于温度和压力的液态乙烯密度
3    %
4    % 输入:
5    %   temperatureK - 温度（开尔文，有效范围：104K到282K）
6    %   pressureBar - 压力（bar，有效范围：1到100 bar）
7    %
8    % 输出:
9    %   density - 液态乙烯的密度（kg/m³）
10    
11    % 乙烯的常数
12    A = 700;
13    Tc = 283.18;  % 临界温度（K）
14    n = 0.29683;
15    B = 0.8;
16    kappa = 0.00125;  % 等温压缩系数（MPa⁻¹）
17    P_ref = 0.1;  % 参考压力（MPa，1 bar）
18    
19    % 将压力从bar转换为MPa
20    pressureMPa = pressureBar / 10;
21    
22    % 计算参考压力下的密度
23    rho_ref = A * (1 - temperatureK/Tc)^n - B * temperatureK;
24    
25    % 应用压力修正
26    density = rho_ref * (1 + kappa * (pressureMPa - P_ref));
27end
28
29% 示例用法
30temp = 200;  % K
31pressure = 50;  % bar
32density = ethyleneDensity(temp, pressure);
33fprintf('液态乙烯在%gK和%g bar下的密度: %.2f kg/m³\n', temp, pressure, density);
34

1#include <iostream>
2#include <cmath>
3
4/**
5 * 计算基于温度和压力的液态乙烯密度。
6 * 
7 * @param temperatureK 温度（开尔文，有效范围：104K到282K）
8 * @param pressureBar 压力（bar，有效范围：1到100 bar）
9 * @return 液态乙烯的密度（kg/m³）
10 */
11double calculateEthyleneDensity(double temperatureK, double pressureBar) {
12    // 乙烯的常数
13    const double A = 700.0;
14    const double Tc = 283.18;  // 临界温度（K）
15    const double n = 0.29683;
16    const double B = 0.8;
17    const double kappa = 0.00125;  // 等温压缩系数（MPa⁻¹）
18    const double P_ref = 0.1;  // 参考压力（MPa，1 bar）
19    
20    // 将压力从bar转换为MPa
21    double pressureMPA = pressureBar / 10.0;
22    
23    // 计算参考压力下的密度
24    double rho_ref = A * pow(1.0 - temperatureK/Tc, n) - B * temperatureK;
25    
26    // 应用压力修正
27    double rho = rho_ref * (1.0 + kappa * (pressureMPA - P_ref));
28    
29    return rho;
30}
31
32int main() {
33    double temp = 200.0;  // K
34    double pressure = 50.0;  // bar
35    double density = calculateEthyleneDensity(temp, pressure);
36    
37    std::cout << "液态乙烯在" << temp << "K和" 
38              << pressure << " bar下的密度: " << density << " kg/m³" << std::endl;
39    
40    return 0;
41}
42

1public class EthyleneDensityCalculator {
2    /**
3     * 计算基于温度和压力的液态乙烯密度。
4     * 
5     * @param temperatureK 温度（开尔文，有效范围：104K到282K）
6     * @param pressureBar 压力（bar，有效范围：1到100 bar）
7     * @return 液态乙烯的密度（kg/m³）
8     */
9    public static double calculateEthyleneDensity(double temperatureK, double pressureBar) {
10        // 乙烯的常数
11        final double A = 700.0;
12        final double Tc = 283.18;  // 临界温度（K）
13        final double n = 0.29683;
14        final double B = 0.8;
15        final double kappa = 0.00125;  // 等温压缩系数（MPa⁻¹）
16        final double P_ref = 0.1;  // 参考压力（MPa，1 bar）
17        
18        // 将压力从bar转换为MPa
19        double pressureMPA = pressureBar / 10.0;
20        
21        // 计算参考压力下的密度
22        double rhoRef = A * Math.pow(1.0 - temperatureK/Tc, n) - B * temperatureK;
23        
24        // 应用压力修正
25        double rho = rhoRef * (1.0 + kappa * (pressureMPA - P_ref));
26        
27        return rho;
28    }
29    
30    public static void main(String[] args) {
31        double temp = 200.0;  // K
32        double pressure = 50.0;  // bar
33        double density = calculateEthyleneDensity(temp, pressure);
34        
35        System.out.printf("液态乙烯在%.1fK和%.1f bar下的密度: %.2f kg/m³%n", 
36                          temp, pressure, density);
37    }
38}
39