A new numerical method is described for the fully-iterated, conjugate solution of two discrete sub-models, involving (1) a transport network model for heat, moisture, and air flows in a high-permeability, fluid-filled cavity; and (2) a fractured porous medium. The transport network sub-model is an integrated parameter, Computational Fluid Dynamics (CFD) solver, describing the thermal-hydrologic transport processes in the flow channel system of the cavity with laminar or turbulent flow and convective heat and mass transport, using MULTIFLUX.The porous medium submodel, using TOUGH2, is a solver for the heat and mass transport in the