Diverse subsonic initial-boundary-value problems (flows in a closed volume initiated by blowing or suction through permeable walls, flows with continuously distributed sources, viscous flows with substantial heat fluxes, etc.) are considered, to show that they cannot be solved by using the classical theory of incompressible fluid motion which involves the equation div u = 0. Application of the most general theory of compressible fluid flow may not be best in such cases, because then we encounter difficulties in accurately resolving the complex acoustic phenomena as well as in assigning the proper boundary conditions. With this in mind a new non-local mathematical model, where div u ≠ 0 in the general case, is proposed for the simulation of unsteady subsonic flows in a bounded domain with continuously distributed sources of mass, momentum and entropy, also taking into account the effects of viscosity and heat conductivity when necessary. The exclusion of sound waves is one of the most important features of this model which represents a fundamental extension of the conventional model of incompressible fluid flow. The model has been built up by modifying both the general system of equations for the motion of compressible fluid (viscous or inviscid as required) and the appropriate set of boundary conditions. Some particular cases of this model are discussed. A series of exact time-dependent solutions, one- and two-dimensional, is presented to illustrate the model.