Capturing a level of modeling of the flow inside a multistage turbomachine, such as unsteadiness for example, can be done at different levels of detail, either by capturing all deterministic features of the flow with a pure unsteady method or by settling for an approximated solution at a lower computational cost. The harmonic methods stand in this second category. Among them, the “nonlinear harmonic method” (NLHM) from He and Ning [1998, “Efficient Approach for Analysis of Unsteady Viscous Flows in Turbomachines,” AIAA J., 36, pp. 2005–2012] revealed the most efficient. This method consists of solving the fully nonlinear 3D steady problem and a linearized perturbation system in the frequency domain. As it has been shown by the authors that the circumferential variations exhibit a harmonic behavior, it is proposed here to adapt the NLHM to the throughflow model, where the main nonlinear system would be the common throughflow equations and the auxiliary system would give access to the circumferential stresses. As the numerical local explicit impermeability conditions are unsupported by Fourier series, the adaptation of this technique to the throughflow model relies on a reformulation of the blade effect by a smooth force field as in the “immersed boundary method” from Peskin [2002, “The Immersed Boundary Method,” Acta Numerica, 11, pp. 1–39]. A simple example of an inviscid flow around a cylinder will illustrate the preceding developments, bringing back the mean effect of the circumferential nonuniformities into the meridional flow.