yade-users team mailing list archive
Mailing list archive
Re: [Question #698948]: Stability of advection modeling by ThermalEngine
Question #698948 on Yade changed:
Zoheir Khademian posted a new comment:
Sorry for my delays. I try to give your comments enough thoughts and do investigation before replying.
I tried both released version  and development version  and still
observed the instabilities meaning the particle temp goes above the max
boundary temp while cell temp is still within the range.
Regarding the diffusivity coefficient kA/L, in the MWE I uploaded last, I have reduced kA/L by reducing fluidk and fluidConductionAreaFactor to arbitrarily low values. I also tried your values (thermal.minimumFluidCondDist=rMean and thermal.fluidConductionAreaFactor=0.1) but still the particle temp is unstable.
Finally I looked at the thermal.cpp and noticed that
applyTempDeltaToSolids  adds pore temp differences to the solid temp
for no apparent reason (particle temp is updated here ). I
deactivated applyTempDeltaToSolids and recompiled the source. I ran the
MWE and the results are stable now.
Let me know if applyTempDeltaToSolids must not be deactivated or I am
missing something here
You received this question notification because your team yade-users is
an answer contact for Yade.