Abstract
A computational method was investigated to deal with the fluid-structure thermal interaction as well as the mechanical interactions. For this purpose, our multiphase model was improved in order to estimate non-isotropic thermal conductivity in the computational cell including different phases using subcell numerical procedures. The proposed computational method was applied to the natural convection flows driven by temperature difference in the fluid-saturated porous media, consisting of heat conductive solid matrix. As a result of computations, the validity of non-isotopic model was confirmed by comparing the theoretical results. In addition, it was shown that the natural convection flows, which arise in the vertical and horizontal rectangular boxes filled with many spherical objects. are reasonably predicted in the present method.
