In this study, we develop computational models and a methodology for accurate multicomponent flow simulation in underresolved multiscale porous structures [1]. It is generally impractical to fully resolve the flow in porous structures with large length-scale differences due to the tremendously high computational expense. The flow contributions from underresolved scales should be taken into account with proper physics modeling and simulation processes. Using precomputed physical properties such as the absolute permeability, K 0 , the capillary pressure-saturation curve, and the relative permeability, K r , in typical resolved porous structures, the local fluid force is determined and applied to simulations in the underresolved regions, which are represented by porous media. In this way, accurate flow simulations in multiscale porous structures become feasible. To evaluate the accuracy and robustness of this method, a set of benchmark test cases are simulated for both single-component and two-component flows in artificially constructed multiscale porous structures. Using comparisons with analytic solutions and results with much finer resolution resolving the porous structures, the simulated results are examined. Indeed, in all cases, the results successfully show high accuracy with proper input of K 0 , capillary pressure, and K r . Specifically, imbibition patterns, entry pressure, residual component patterns, and absolute/relative permeability are accurately captured with this approach.