2026年大学化工技术（化工热力学）试题及答案

2026年大学化工技术（化工热力学）试题及答案一、选择题（每题2分，共20分）1.对于实际气体的焓变计算，若已知某状态下的剩余焓HR，其与理想气体焓变ΔHig的关系为（）。A.ΔH=ΔHig+HRB.ΔH=ΔHigHRC.ΔH=HRΔHigD.ΔH=HR2.某二元溶液中组分1的活度系数γ1=2.5，若该溶液为正规溶液，则组分2的活度系数γ2应满足（）。（已知x1=0.3，x2=0.7）A.γ2=1.0B.γ2=2.5C.γ2=(x1/x2)×2.5D.γ2=exp[(x1²)(gE/RT)]3.下列关于麦克斯韦关系的表述中，错误的是（）。A.(∂T/∂V)S=-(∂P/∂S)VB.(∂T/∂P)S=(∂V/∂S)PC.(∂S/∂V)T=(∂P/∂T)VD.(∂S/∂P)T=-(∂V/∂T)P4.用RK方程计算373K、5MPa下甲烷的压缩因子Z时，若已知RK方程常数a=0.42748R²Tc².5/Pc，b=0.08664RTc/Pc，其中甲烷的Tc=190.6K，Pc=4.596MPa，则a的单位应为（）。A.Pa·m⁶·K⁰.5/mol²B.Pa·m⁶·K⁰.5/molC.Pa·m³·K⁰.5/mol²D.Pa·m³·K⁰.5/mol5.某纯物质在临界点处的特性为（）。A.(∂P/∂V)T=0，(∂²P/∂V²)T>0B.(∂P/∂V)T=0，(∂²P/∂V²)T=0C.(∂P/∂V)T>0，(∂²P/∂V²)T=0D.(∂P/∂V)T<0，(∂²P/∂V²)T<06.对于理想溶液，下列关系式中不成立的是（）。A.GE=0B.HE=0C.SE=0D.μi=μi+RTlnxiD.μi=μi+RTlnxi7.某气体在298K、100kPa下的逸度系数φ=0.95，若该气体在相同温度下压力升至200kPa时，其逸度f的近似值为（）（假设压力变化范围内热容不变，且可忽略动能和势能变化）。A.190kPaB.195kPaC.200kPaD.185kPa8.用Clausius-Clapeyron方程计算纯液体的饱和蒸气压时，若已知汽化热ΔHv=40kJ/mol，温度T1=373K时P1=101.325kPa，则T2=383K时的P2约为（）（R=8.314J/mol·K）。A.143kPaB.120kPaC.155kPaD.170kPa9.某系统经历一个不可逆绝热过程，其熵变ΔS（）。A.大于0B.等于0C.小于0D.无法确定10.对于二元汽液平衡系统，若满足理想气体（气相）和非理想溶液（液相），则相平衡条件为（）。A.yiP=xiγiPi^sB.yiP=xiPi^sC.yiφiP=xiγiPi^sD.yiφiP=xiPi^s二、填空题（每空1分，共15分）1.热力学基本方程dU=TdS-PdV的适用条件是__________、__________。2.范德华方程的表达式为__________，其中常数a反映了__________，b反映了__________。3.超额吉布斯自由能GE与活度系数的关系为__________（以二元溶液为例）。4.剩余性质的定义是__________，其物理意义是__________。5.纯物质的汽液平衡条件是__________，液液平衡条件是__________。6.用普遍化压缩因子法计算实际气体的体积时，若对比温度Tr=1.2，对比压力Pr=0.8，查得Z=0.85，则实际体积V=__________（用理想气体体积Vig表示）。7.某溶液的超额熵SE=-5J/mol·K，超额焓HE=200J/mol，则超额吉布斯自由能GE=__________。三、简答题（每题6分，共30分）1.简述理想气体与理想溶液的区别，并分别写出其化学势表达式。2.说明为何实际气体的逸度系数φ在低压下趋近于1，而高压下可能大于或小于1。3.比较RK方程与PR方程的优缺点，说明PR方程在实际应用中的改进之处。4.对于二元汽液平衡系统，若液相为非理想溶液（正偏差），气相为理想气体，画出其P-x-y图的典型形状（标注坐标轴、泡点线、露点线及关键点）。5.如何利用热力学基本关系式推导(∂H/∂P)T的表达式？写出推导过程并说明其物理意义。四、计算题（共35分）1.（10分）已知某气体的状态方程为P(V-b)=RT（b为常数），试推导该气体的剩余焓HR和剩余熵SR的表达式（以T、P为变量）。2.（12分）用RK方程计算300K、3MPa下丙烷的摩尔体积（结果保留三位有效数字）。已知丙烷的临界参数：Tc=369.8K，Pc=4.25MPa，R=8.314J/mol·K；RK方程形式为P=RT/(V-b)a/[T^0.5V(V+b)]，其中a=0.42748R²Tc².5/Pc，b=0.08664RTc/Pc。3.（13分）某二元溶液在323K下的超额吉布斯自由能模型为GE/RT=0.5x1x2，其中x1为组分1的摩尔分数。已知组分1和组分2的饱和蒸气压分别为P1^s=80kPa，P2^s=40kPa，气相可视为理想气体。（1）写出组分1和组分2的活度系数γ1、γ2的表达式；（2）计算x1=0.4时的泡点压力P及平衡气相组成y1；（3）判断该溶液是否形成共沸物，说明理由。答案一、选择题1.B2.D3.A4.A5.B6.C7.A8.A9.A10.A二、填空题1.封闭系统、可逆过程（或组成不变的均相系统）2.(P+a/Vm²)(Vm-b)=RT；分子间的吸引力；分子本身的体积3.GE=RT(x1lnγ1+x2lnγ2)4.实际状态与同温同压下理想气体状态的热力学性质之差；衡量实际气体与理想气体的偏差5.μv=μl；μl1=μl2（同一组分在两相中的化学势相等）6.0.85Vig（或Z×Vig）7.GE=HETSE=200-300×(-5)=200+1500=1700J/mol（注：需明确温度，假设T=300K）三、简答题1.理想气体假设分子间无相互作用且分子体积可忽略，其化学势μi=μi^θ(T)+RTln(P/Pθ)；理想溶液假设分子间作用力相等，混合时无体积和焓变，其化学势μi=μi(T,P)+RTlnxi（μi为同温同压下纯组分的化学势）。两者区别：理想气体是单相分子无作用的模型，理想溶液是多组分分子间作用均等的模型。1.理想气体假设分子间无相互作用且分子体积可忽略，其化学势μi=μi^θ(T)+RTln(P/Pθ)；理想溶液假设分子间作用力相等，混合时无体积和焓变，其化学势μi=μi(T,P)+RTlnxi（μi为同温同压下纯组分的化学势）。两者区别：理想气体是单相分子无作用的模型，理想溶液是多组分分子间作用均等的模型。2.低压下，实际气体分子间距大，分子间作用力可忽略，趋近于理想气体（φ→1）；高压下，若分子间吸引力主导（如低温高压），则实际压力低于理想气体（φ<1）；若分子间排斥力主导（如高温高压），则实际压力高于理想气体（φ>1）。3.RK方程引入温度修正项a/T^0.5，改进了范德华方程的精度，适用于非极性气体；PR方程将a项修正为a(T)=a(Tr)[1+k(1-Tr^0.5)]²（k为偏心因子函数），并调整了b的系数（b=0.07780RTc/Pc），对极性和非极性物质的相平衡计算更准确，尤其在液相密度预测上优于RK方程。4.P-x-y图中，横坐标为x1（液相组成），纵坐标为P（压力）。泡点线（P-x）位于下方，露点线（P-y）位于上方，因正偏差（γi>1），泡点线低于理想溶液线（拉乌尔定律线），露点线高于理想溶液线，两线在x1=0（P=P2^s）和x1=1（P=P1^s）处交汇。5.由H=U+PV，微分得dH=TdS+VdP（热力学基本方程）。在恒温下对P求偏导：(∂H/∂P)T=T(∂S/∂P)T+V。根据麦克斯韦关系(∂S/∂P)T=-(∂V/∂T)P，代入得(∂H/∂P)T=V-T(∂V/∂T)P。其物理意义是恒温下压力变化对焓的影响，反映了体积随温度变化的特性（如液体的(∂H/∂P)T≈V，因液体膨胀系数小）。四、计算题1.剩余焓HR=H-Hig，Hig=∫CpdT（与P无关），实际H=U+PV。由状态方程P(V-b)=RT，得V=RT/P+b。dU=TdS-PdV（基本方程），结合dS=(∂S/∂T)VdT+(∂S/∂V)TdV，由麦克斯韦关系(∂S/∂V)T=(∂P/∂T)V=R/(V-b)（因P=RT/(V-b)）。则dU=T[(Cp/V)dT+R/(V-b)dV]PdV=CvdT+[TR/(V-b)P]dV。由状态方程TR/(V-b)=P，故dU=CvdT，即U仅与T有关，H=U+PV=U+RT+Pb（因PV=RT+Pb）。Hig=Uig+PVig=Uig+RT（理想气体PV=RT）。HR=H-Hig=(U+RT+Pb)-(Uig+RT)=(U-Uig)+Pb。因实际气体与理想气体在T相同时U差为0（假设Cv相同），故HR=Pb。同理，SR=S-Sig，Sig=∫Cp/TdTRlnP（理想气体），实际S=∫(∂S/∂T)VdT+∫(∂S/∂V)TdV=∫Cv/TdT+∫R/(V-b)dV=CvlnT+Rln(V-b)+C。Sig=CvlnTRlnP+C（理想气体Vig=RT/P，Vig-b≈Vig，故ln(Vig-b)=ln(RT/P)=lnT+lnRlnP）。SR=S-Sig=Rln(V-b)(-RlnP)=Rln[(V-b)P/RT]（因V-b=RT/P，故(V-b)P/RT=1），因此SR=0。2.计算a和b：a=0.42748×(8.314)²×(369.8)^2.5/4.25×10^6先计算Tc².5=369.8^2×369.8^0.5≈(136752.04)×19.23≈2.63×10^6R²=69.12，故a≈0.42748×69.12×2.63×10^6/4.25×10^6≈0.42748×69.12×2.63/4.25≈0.42748×42.8≈18.36(Pa·m⁶·K^0.5/mol²)b=0.08664×8.314×369.8/4.25×10^6≈0.08664×3075.7/4.25×10^6≈266.5/4.25×10^6≈6.27×10^-5(m³/mol)RK方程：P=RT/(V-b)a/[T^0.5V(V+b)]代入T=300K，P=3×10^6Pa，RT=8.314×300=2494.2J/mol整理方程：3×10^6=2494.2/(V-6.27×10^-5)18.36/[300^0.5V(V+6.27×10^-5)]令V的初始猜测为理想气体体积Vig=RT/P=2494.2/3×10^6≈8.314×10^-4m³/mol（0.8314L/mol），代入计算右边第一项≈2494.2/(8.314×10^-4-6.27×10^-5)=2494.2/(7.687×10^-4)≈3.24×10^6Pa第二项≈18.36/(17.32×8.314×10^-4×8.941×10^-4)≈18.36/(17.32×7.43×10^-7)≈18.36/(1.287×10^-5)≈1.426×10^6Pa右边≈3.24×10^6-1.426×10^6=1.814×10^6<3×10^6，说明V偏小，需增大V。假设V=1.0×10^-3m³/mol（1L/mol），则V-b≈9.373×10^-4m³/mol第一项=2494.2/9.373×10^-4≈2.66×10^6Pa第二项=18.36/(17.32×1.0×10^-3×1.0627×10^-3)=18.36/(17.32×1.0627×10^-6)≈18.36/(1.84×10^-5)≈9.98×10^5Pa右边≈2.66×10^6-9.98×10^5≈1.66×10^6<3×10^6，仍偏小。改用迭代法，设V=5×10^-4m³/mol（0.5L/mol），V-b≈4.373×10^-4m³/mol第一项=2494.2/4.373×10^-4≈5.70×10^6Pa第二项=18.36/(17.32×5×10^-4×5.627×10^-4)=18.36/(17.32×2.813×10^-7)=18.36/(4.87×10^-6)≈3.77×10^6Pa右边≈5.70×10^6-3.77×10^6=1.93×10^6<3×10^6，仍小。可能初始假设错误，实际气体在3MPa、300K（Tr=300/369.8≈0.81）下处于液相区，需用液相体积计算。查丙烷液相摩尔体积约为1.5×10^-4m³/mol（0.15L/mol），代入：V=1.5×10^-4m³/mol，V-b≈8.73×10^-5m³/mol第一项=2494.2/8.73×10^-5≈2.86×10^7Pa（远大于3MPa），说明V应更大。正确方法：将RK方程整理为Z^3Z^2+(ABB^2)ZAB=0，其中A=aP/(R²T^2.5)，B=bP/(RT)计算A=18.36×3×10^6/(8.314²×300^2.5)=5.508×10^7/(69.12×5.196×10^6)=5.508×10^7/(3.59×10^8)=0.153B=6.27×10^-5×3×10^6/(8.314×300)=188.1/2494.2≈0.0754方程：Z³Z²+(0.153-0.0754-0.0754²)Z-0.153×0.0754=0计算系数：0.153-0.0754=0.0776，0.0754²≈0.00568，故中间项=0.0776-0.00568=0.0719方程：Z³Z²+0.0719Z-0.0115=0试Z=0.8：0.512-0.64+0.0575-0.0115=-0.082Z=0.9：0.729-0.81+0.0647-0.0115=-0.0278Z=0.95：0.857-0.9025+0.0683-0.0115=0.0113用线性插值：Z=0.94时，0.8306-0.8836+0.0676-0.0115=-0.0069；Z=0.95时0.0113，故Z≈0.945则V=ZRT/P=0.945×2494.2/3×10^6≈0.945×8.314×10^-4≈7.86×10^-4m³/mol（0.786L/mol）3.（1）GE/RT=0.5x1x2，由GE=RT(x1lnγ1+x2lnγ2)，对x1求偏导得lnγ1=(∂(GE/RT)/∂x1)x2=0.5x2²（因x2=1-x1，d(x1x2)/dx1=x2-x1，故导数为0.5(1-2x1)？错误，正确方法：GE/RT=0.5x1(1-x1)=0.5x1-0.5x1²，lnγ1=(∂(GE/RT)/∂x1)T,P,x2=0.5-x1（当x2=1-x1时，d(x1x2)/dx1=1-2x1，故GE/RT=0.5x1x2，∂(GE/RT)/∂x1=0.5x2-0.5x1=0.5(1-2x1)？不，正确的偏导应为保持x2=1-x1，所以d(x1x2)=x2dx1+x1dx2=x2dx1-x1dx1=(x2-x1)dx1，因此∂(GE/RT)/∂x1=0.5(x2-x1)=0.5(1-2x1)。但更简