Добавление членов в уравнение энергии для твердого тела

Здравствуйте, ну только если переписать с нуля все уравнения с помощью User Defined Scalar.

А какой дополнительный член вы хотите? Это не источник?

Здравствуйте!

Всё-таки это очень похоже на источник!

Вам нужно просто сообразить, какую математику использовать для нахождения производной плотности по времени.

Производную источника по переменной уравнения энергии (по температуре) можно занулить. Всё равно будет работать

Пример добавления теплоты конденсации в уравнение энергии в однофазной постановке:

DEFINE_SOURCE(my_energy_source,c,t,dS,eqn)
{
    real source;    /* vapor source */
    real sigma = 1.0;    /* evaporation efficiency */
    real constant = sigma *
        sqrt( mw_vapor / (2.0 * M_PI* UNIVERSAL_GAS_CONSTANT * C_T(c,t)));    /* condensation relaxation time, 1/s */
    real diameter_drop = 2.0e-5; /* 10-50 mkm */
    real latent_heat = 2258000.0; /* latent_heat, J/kg */

    
    real X_vap = C_YI(c, t, 0);  /* mass fraction of vapor */
    real Y_vap = molefrac(X_vap, mw_vapor, mw_air);  /* mole fraction, alpha vapor */
    real X_drop = C_YI(c, t, 1);  /* mass fraction of droplets */

    real area_density = 6.0/diameter_drop*X_drop*C_R(c,t)/1000.0; 

    real P_operating = RP_Get_Real("operating-pressure");
    real P_sat = H2O_Saturation_Pressure(C_T(c,t));

    real P_vap = P_operating*Y_vap;

    if (C_UDMI(c, t, UDM_RH) > 101.0)
    {
    
    if (X_drop > 0.0)
    {    source = -latent_heat*area_density*Y_vap*constant*P_sat*((100.0-C_UDMI(c, t, UDM_RH))/100.0);
    dS[eqn] = 0.0;    }
    
    
    else
    {source = 1.0e-7*latent_heat;
    dS[eqn] = 0.0;}

    }

    else 
  {
    source = 0.0;
    dS[eqn] = 0.0;
  }
      
      C_UDMI(c, t, UDM_RH_new) = source;
    return source;
}