Recently, I’ve been digging into the code and several questions have appeared.
From the theory (1993 file), I believed the solids phase conductivity was determined using Kuipers et al. (1992) model (pre GUI versions) from the actual thermal conductivity of the solids material.
In the new versions a just saw a new keyword: “KS_MODEL”, which permits the selection of the model, where Bauer’s model is the only model available for it (indeed similar to Kuipers).
However, going into calc_k_s.f, I see that if K_s0 is defined it is taken as the effective thermal conductivity of the solids phase, and Bauer’s model is only used if K_s0 is undefined. Indeed, from calc_k_s.f, Bauer’s model takes a solids conductivity as if the dense phase were made of ash particles (ks_micro = 2.2 W/mK).
Besides, I cannot track the keyword “KS_MODEL” within the subroutines to check how each option is selected.
My questions are:
If I defined the solids material conductivity and activate Bauer’s model, how is the solids phase effective conductivity calculated?
If Bauer’s model always uses ks_micro for the solids material conductivity, do I need to modify the subroutine (or implement it through an user defined) to calculate the solids phase effective conductivity using Bauer’s model?
What is the default value for “KS_MODEL”?
Thanks in advance.