其他三个都能修改，修改定压比热总会报错no user-defined functions of type udf-type-specific-heat have been loaded，望大佬赐教 有偿
#include "udf.h"
#include <math.h>
/*----------- 用户可调参数 -----------*/
#define PHI 0.05 // 纳米颗粒体积分数(5%)
#define N 3.0 // 形状因子（球形颗粒n=3）
#define RHO_P 2200 // SiO2密度 kg/m³
#define CP_P 745 // SiO2比热容 J/kg-K
#define K_P 1.4 // SiO2导热系数 W/m-K
/*--------- 基液CO₂属性函数(8.2MPa) ----------*/
real co2_density(real T) {
/* 密度多项式拟合(273-1000K, 8.2MPa) */
return 682.1 + 1.873 * T - 0.00412 * pow(T, 2) + 3.15e-6 * pow(T, 3); // kg/m³
}
real co2_cp(real T) {
return 890.4 + 2.15 * T - 0.00582 * pow(T, 2) + 4.27e-6 * pow(T, 3); // J/kg-K
}
real co2_conductivity(real T) {
/* 导热系数拟合 */
return 0.0148 + 2.25e-5 * T - 1.1e-8 * pow(T, 2); // W/m-K
}
real co2_viscosity(real T) {
/* 粘度Andrade方程 */
return 1.62e-5 * exp(1350 / (T - 45.3)); // Pa·s
}
/*----------- 纳米流体物性UDF -----------*/
DEFINE_PROPERTY(nano_density, c, t)
{
real T = C_T(c, t);
real rho_bf = co2_density(T);
return PHI * RHO_P + (1 - PHI) * rho_bf; // 混合规则
}
DEFINE_PROPERTY(nano_cp, c, t)
{
// real T = C_T(c, t);
// real rho_bf = co2_density(T);
//real cp_bf = co2_cp(T);
// real numerator = PHI * RHO_P * CP_P + (1 - PHI) * rho_bf * cp_bf;
// real denominator = PHI * RHO_P + (1 - PHI) * rho_bf;
//return numerator / denominator; // 质量加权平均
return 1000.0;
}
DEFINE_PROPERTY(nano_conductivity, c, t)
{
real T = C_T(c, t);
real k_bf = co2_conductivity(T);
real numerator = K_P + (N - 1) * k_bf - (N - 1) * PHI * (k_bf - K_P);
real denominator = K_P + (N - 1) * k_bf + PHI * (k_bf - K_P);
return (numerator / denominator) * k_bf; // 修正Maxwell模型
}
DEFINE_PROPERTY(nano_viscosity, c, t)
{
real T = C_T(c, t);
real mu_bf = co2_viscosity(T);
return mu_bf / pow((1 - PHI), 2.5); // Brinkman模型
}