← Back to team overview

yade-users team mailing list archive

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