Python 程式碼實作熱力學循環效率計算

根據題目要求，我們需要求解η(x)=0.3对应的T1/T2的值。首先，我们需要推导出η(x)的表达式。

根据题目中给出的PV图，我们可以得到以下信息：

• 定温过程：'PV^γ=constant'
• 定压过程：'PV=constant'
• 定容过程：'P/T=constant'

设定温过程的初始状态为1，终态为2，定压过程的初始状态为2，终态为3，定容过程的初始状态为3，终态为1。

根据上述信息，我们可以得到以下关系：

• 定温过程：'T_1V_1^(γ-1)=T_2V_2^(γ-1)'
• 定压过程：'P_2V_2=P_3V_3'
• 定容过程：'P_3/T_3=P_1/T_1'

根据热效率的定义，我们有：

'η=W/Q_h=1-Q_c/Q_h=1-Q_c/Q_h=Q_h-Q_c/Q_h=1-Q_c/Q_h=1-T_3(S_3-S_2)/T_1(S_1-S_4)'

其中，'S_1'、'S_2'、'S_3'、'S_4'分别表示过程1、2、3、4的熵。

根据熵的定义，我们有：

'S=CVln(T)+Rln(V)'

将以上关系带入热效率的定义中，我们可以得到：

'η=1-T_3(CV_3ln(T_3)+Rln(V_3)-CV_2ln(T_2)-Rln(V_2))/T_1(CV_1ln(T_1)+Rln(V_1)-CV_4ln(T_4)-Rln(V_4))'

其中，'V_4'表示过程4的体积，由定容过程的关系'P_3/T_3=P_1/T_1'可以得到'V_4=V_3(T_3/T_1)'

由于题目中给出的是温度的比值，我们可以将上式以温度的比值形式表示：

'η=1-T_3(CV_3ln(T_3)+Rln(V_3)-CV_2ln(T_2)-Rln(V_2))/T_1(CV_1ln(T_1)+Rln(V_1)-CV_4ln(T_4)-Rln(V_4))=T_3-T_2/T_1-T_4'

将以上关系整理后即可得到η(x)的表达式。

接下来，我们可以使用scipy中的数值求解方法来求解η(x)=0.3的解。以下是使用scipy的代码示例：

import numpy as np
from scipy.optimize import fsolve

# 定义η(x)的表达式
def eta(x):
    T1 = 1
    T2 = x
    T3 = x
    T4 = 1
    gamma = 5/3
    CV = 1
    R = 1
    V2 = 1
    V3 = 1
    V4 = V3 * (T3/T1)
    numerator = T3 * (CV * V3 * np.log(T3) + R * np.log(V3) - CV * V2 * np.log(T2) - R * np.log(V2))
    denominator = T1 * (CV * V1 * np.log(T1) + R * np.log(V1) - CV * V4 * np.log(T4) - R * np.log(V4))
    return (T3 - T2) / (T1 - T4) - 0.3

# 使用fsolve求解η(x)=0.3的解
x_solution = fsolve(eta, 1)
T1_T2 = x_solution[0]
print("T1/T2 =", T1_T2)

运行以上代码后，可以得到T1/T2的值。