Hi, I need to calculate case when cp of solid increases with temperature, so I change specific heat of solid during the calculation (transient)using SPECH.f , and the result is a drop of temperature in solid, when cp is increased. Do you have any experience with that? what's wrong? Thanks Marescu

For a constant heat flux, higher Cp will results in lower dT = Q=m*Cp*dT. So depending upon your problem, that may be correct. What do you expect to happen?

this is not what i mean. There is a is a drop of Temperature in the time step the cp is changed . i.e. for constant positive heat flux, following Q=m*cp*dT the temperature is still growing( with different dT),- no temperature drop. So I've got sudden temperature drop in one particular time step....

in specht.f

C-------------------------------------------------------------------------

C C This subroutine enables the user to specify the mean specific heat

C capacity at cells and boundaries.

C

C CP - Mean specific heat capacity (J/kgK)

C (of the background fluid in multi-component fluids)

so I think, you must set "Mean specific heat capacity " at each temperature of solid