Dear developers,
I try to use des_usr_var to output convection source in global energy source array.According to the description, Qcv is the heat transfer between the particle and the fluid, but the output value during my calculation is 0.I don’t know why,I use the dmp solver. solar_03_2022-08-12T045253.996079.zip (297.1 KB)
! Calculate the heat transfer coefficient.
CALL CALC_GAMMA_DES(NP, GAMMA)
if (gamma .ne. 0.0) write(*,*) "********GAMMA", gamma
! Calculate heat transfer for a real particle in the parcel, then scale W times
if(CGDEM) then
GAMMAxSA(NP) = GAMMA* l4Pi*DES_CGP_RPR(NP)*DES_CGP_RPR(NP) *DES_CGP_STW(NP)
else
GAMMAxSA(NP) = GAMMA* l4Pi*DES_RADIUS(NP)*DES_RADIUS(NP)
endif
! Calculate the rate of heat transfer to the particle
Qcv = GAMMAxSA(NP)*(lTg - DES_T_s(NP))
the value of GAMMA cming from CALC_GAMMA_DES is 0, which causes GAMMAxSA and Qcv to both be 0.
So the question is, why does CALC_GAMMA_DES always return 0?