Abstract Winds from black hole accretion flows are ubiquitous. Previous works mainly focus on the launching of wind on the accretion flow scale. It still remains unclear how far the winds can propagate outward and what their large-scale dynamics is. As the first paper of this series, we study the large-scale dynamics of thermal wind beyond accretion scales via analytical and numerical methods. Boundary conditions, which are crucial to our problem, are analyzed and presented based on small-scale simulations combined with observations of winds. Both black hole and galaxy potential are taken into account. For winds originating from hot accretion flows, we find that the wind can reach large scales. The radial profiles of velocity, density, and temperature can be approximated by and , where v r0, ρ 0, and T 0 are the velocity, density, and temperature of winds at the boundary , and γ is the polytropic index. During the outward propagation, the enthalpy and rotational energy compensate for the increase of gravitational potential. For thin disks, we find that because the Bernoulli parameter is smaller, winds cannot propagate as far as the hot winds, but stop at a certain radius where the Bernoulli parameter is equal to the potential energy. Before the winds stop, the profiles of dynamical quantities can also be approximated by the above relations. In this case, the rotational energy alone compensates for the increase in potential energy.