yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #26474
Re: [Question #698948]: Stability of advection modeling by ThermalEngine
Question #698948 on Yade changed:
https://answers.launchpad.net/yade/+question/698948
Zoheir Khademian posted a new comment:
Hey Robert,
Sorry for my delays. I try to give your comments enough thoughts and do investigation before replying.
I tried both released version [1] and development version [2] 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 [3] adds pore temp differences to the solid temp
for no apparent reason (particle temp is updated here [4]). 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
Thanks
Zoheir
[1] https://launchpad.net/yade/+download
[2] https://gitlab.com/yade-dev/trunk
[3] https://gitlab.com/yade-
dev/trunk/-/blob/master/pkg/pfv/Thermal.cpp#L94
[4] https://gitlab.com/yade-
dev/trunk/-/blob/master/pkg/pfv/Thermal.cpp#L605
--
You received this question notification because your team yade-users is
an answer contact for Yade.
Follow ups