A domain decomposition method for steady-state, subsonic fluid dynamics calculations is presented. An implicit solver is used in each subdomain, with boundary data from previous iterations. A nonlinear numerical boundary condition is used to enhance convergence. Numerical examples are presented for the Euler equations of gas dynamics and a two-phase flow model.