GEKKO动态热仿真中CoolProp两相区BSpline拟合结果异常咨询

• 实现逻辑差异：GEKKO内置bspline采用全局最小二乘多项式拟合框架，并非Scipybisplrep使用的局部分块过点B样条，对带陡变、平台特征的物性曲面会产生显著龙格过冲，即观测到的温度突跃问题。
• 物性特性不匹配：两相区的T(P,h)映射存在固有非光滑特性——定压下从饱和液到饱和气的焓区间内，温度恒等于对应压力的饱和温度，形成长平台结构，全局样条无足够节点自由度精准贴合这类转折，会在平台边缘出现幅值远超真实范围的过冲。
• 拟合模式差异：GEKKO内置bspline为最小二乘近似拟合，并非严格通过所有采样点的插值，当采样点分布不均匀、存在非光滑结构时，拟合残差会被多项式放大为显著的曲线过冲。
• 冗余建模：两相区温度仅为压力的单值函数，无需构建h→p→T的二维映射，跨相变边界的全局二维插值本身就放大了拟合误差。

1. 按热力学域拆分插值映射，规避跨相变全局拟合 #

不要将过冷液区、两相区、过热区的物性放在同一个插值域中，分别构建独立的拟合块：

• 单相区（过冷液/过热气）：该区域温度对压力、比焓均为光滑连续函数，可直接使用GEKKO原生bspline实现，拟合精度可满足求解要求
• 两相区：跳过二维焓-压力映射，直接构建一维p→饱和温度Tsat(p)的插值关系，从根源上消除焓平台带来的拟合误差

2. 优化插值采样点分布 #

在同一插值域内，加密物性变化梯度大的位置的采样密度，减少全局拟合的过冲幅度：

• 相变边界（饱和液焓hf(p)、饱和气焓hg(p)）两侧额外补充3-5个密采样点，覆盖转折区域
• 保持采样点维度对齐：调用m.bspline前验证数据维度匹配，即z_data.shape == (len(x_data), len(y_data))，避免因矩阵转置导致的拟合偏差
• 降低插值阶数：若仍需在跨区场景下使用全局bspline，将kx、ky从3调整为2，可有效削弱高次多项式的过冲效应

3. 动态仿真场景的光滑域切换 #

由于GEKKO不支持硬判断的分支逻辑，可通过光滑sigmoid函数实现三个热力学域的计算衔接，所有切换函数均自带解析导数，可直接接入GEKKO求解器：

# 光滑切换逻辑示例
def smooth_step(x, x0, k=0.01):
    return 1/(1 + np.exp(-k*(x - x0)))

# 分域计算温度
tau = (1 - smooth_step(h, hf(p)))*tau_subcool \
      + smooth_step(h, hf(p))*(1 - smooth_step(h, hg(p)))*Tsat(p) \
      + smooth_step(h, hg(p))*tau_superheat

• 输入变量的对数变换无法解决该问题：偏差来源是插值方法的实现差异，并非压力、比焓的尺度范围问题
• 单纯调整插值阶数无法完全消除跨相变区的过冲：全局拟合的本质决定了其对带平台曲面的拟合能力上限，必须通过分域建模才能实现两相区的精准映射

内容的提问来源于stack exchange，提问作者Hagobert