The goal of the research was to demonstrate the impact of thin porous interfacial transition zones (ITZs) between aggregates and cement matrix on fluid flow in unsaturated concrete caused by hydraulic/capillary pressure. To demonstrate this impact, a novel coupled approach to simulate the two‐phase (water and moist air) flow of hydraulically and capillary‐driven fluid in unsaturated concrete was developed. By merging the discrete element method (DEM) with computational fluid dynamics (CFD) under isothermal settings, the process was numerically studied at the meso‐scale in two‐dimensional conditions. A flow network was used to describe fluid behaviour in a continuous domain between particles. Small concrete specimens of a simplified particle mesostructure were subjected to fully coupled hydro‐mechanical simulation tests. A simple uniaxial compression test was used to calibrate the pure DEM represented by bonded spheres, while a permeability and sorptivity test for an assembly of spheres was used to calibrate the pure CFD. For simplified specimens of the pure cement matrix, cement matrix with aggregate, and cement matrix with aggregate and ITZ of a given thickness, DEM/CFD simulations were performed sequentially. The numerical results of permeability and sorptivity were directly compared to the data found in the literature. A satisfactory agreement was achieved. Porous ITZs in concrete were found to reduce sorption by slowing the capillary‐driven fluid flow, and to speed the full saturation of pores when sufficiently high hydraulic water pressures were dominant.