Hybrid numerical method applied to transient hygrothermal analysis in an annular cylinder