Force Partitioning Method (FPM)

• FPM is a fluid-mechanical framework that partitions pressure-induced loads into additive contributions from distinct physical mechanisms, spatial regions, and vortex structures.
• It employs an auxiliary Laplace potential combined with the Q-criterion to clearly separate rotation-dominated from strain-dominated regions in complex flow fields.
• Applications on low-aspect-ratio wings reveal that cross-span vortices can contribute more to lift than spanwise vortices, challenging traditional aerodynamic assumptions.

Searching arXiv for recent and foundational papers on the Force Partitioning Method and closely related formulations. arxiv_search(query="Force Partitioning Method aerodynamic load partitioning Q criterion influence potential Menon Mittal", max_results=10, sort_by="relevance") arxiv_search(query="(Menon et al., 2022) force partitioning lift low aspect ratio wings", max_results=5, sort_by="relevance") arxiv_search(query="(Zhu et al., 2023) force and moment partitioning method swept wings", max_results=5, sort_by="relevance") Force Partitioning Method (FPM) is a fluid-mechanical framework for quantifying how distinct flow features contribute to pressure-induced aerodynamic loads on an immersed body. In the formulation used for incompressible, vortex-dominated flows, an auxiliary harmonic potential encodes the body geometry and the target load direction, and a projection of the Navier–Stokes equations yields additive decompositions of force by physical mechanism, spatial region, and identified flow structure. In the low-aspect-ratio wing study that made the method central to lift attribution, FPM was used to separate spanwise-oriented and cross-span-oriented vorticity, and also rotation-dominated and strain-dominated regions, thereby providing a quantitative account of how different vortex structures generate positive or negative lift (Menon et al., 2022). Closely related formulations extend the same logic to moments, to PIV-based experiments, and to modal decompositions of the $Q$-field in unsteady loading and aeroacoustic problems (Zhu et al., 2023, Prakhar et al., 9 Jan 2025).

1. Definition and analytical scope

FPM is designed to partition the pressure force into additive components attributable to distinct physical mechanisms, spatial regions, and identified flow structures. In the low-aspect-ratio wing formulation, the method partitions the pressure-induced aerodynamic load into contributions from unsteady acceleration, convection, viscous diffusion, and far-field acceleration; into contributions from regions over the wing, in the near wake, and in the far wake; and into contributions from spanwise-oriented versus cross-span-oriented vortices, as well as from rotation-dominated versus strain-dominated regions (Menon et al., 2022).

The method is rigorous in the sense used in the source literature: the partition is obtained from the governing equations rather than from heuristic vortex labeling. Because the relevant force contribution is expressed as a linear volume integral weighted by an influence potential, the contribution from any masked feature or spatial subdomain is additive, and the sum over mutually exclusive partitions recovers the total vortex-induced force. This additivity is the basis for statements such as $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$, where “outside” denotes the residual from outside the analysis subdomain (Menon et al., 2022).

In comparative terms, the method differs from classical surface-pressure versus viscous decompositions, which do not isolate distinct volumetric features such as leading-edge vortex cores and their surrounding strain fields, and from impulse-theory or Lamb-vector-based approaches, which provide integral or global measures but can be less straightforward to partition by local features in complex three-dimensional flows. The source literature characterizes the present FPM, rooted in Quartappelle–Napolitano, Wu, and Howe and extended by Chang, Zhang, and Menon–Mittal, as enabling exact volumetric attribution of pressure-induced loads to identified structures via $-\rho\int 2Q\,\phi_i\,dV$, separation of rotation and strain contributions, and consistent additivity across subdomains and features in fully three-dimensional vortex-dominated flows (Menon et al., 2022).

2. Mathematical foundation

The method is formulated for incompressible flow governed by the Navier–Stokes equations and continuity,

$$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$

FPM introduces an auxiliary scalar potential field $\phi_i$ that depends on the body geometry and on the force direction $i$. For lift on a wing, $i=2$. The potential solves a Laplace problem in the fluid domain $V_f$,

$\nabla^2 \phi_i = 0 \quad \text{in } V_f,$

with mixed Neumann boundary conditions

$$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$

where $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$0 is the body surface and $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$1 is the outer boundary (Menon et al., 2022).

Projecting the momentum equation onto $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$2 and integrating over $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$3 yields a representation of the surface pressure force

$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$4

In the wing formulation, the resulting terms correspond respectively to unsteady acceleration on the body, convective or vortex-induced effects, viscous diffusion, and far-field flow acceleration. For the stationary wings considered there, the unsteady body term is zero and the far-field term is negligible in the large domain; viscous diffusion contributes modestly, about $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$5 of the total lift, while the convective term dominates the pressure-induced lift (Menon et al., 2022).

The connection to vortical structure follows from the pressure Poisson equation,

$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$6

where $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$7 and $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$8 are the strain-rate and rotation-rate parts of $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$9 and

$-\rho\int 2Q\,\phi_i\,dV$0

Substitution yields the central FPM relation for the vortex-induced component,

$-\rho\int 2Q\,\phi_i\,dV$1

Rotation-dominated regions with $-\rho\int 2Q\,\phi_i\,dV$2 and strain-dominated regions with $-\rho\int 2Q\,\phi_i\,dV$3 therefore induce pressure forces of opposite sign. This relation underlies the direct separation of vortex-core and strain-field contributions in later applications (Menon et al., 2022).

A closely related pressure-loading partition used in modal force partitioning writes

$-\rho\int 2Q\,\phi_i\,dV$4

with

$-\rho\int 2Q\,\phi_i\,dV$5

where $-\rho\int 2Q\,\phi_i\,dV$6 is the acceleration-reaction term, $-\rho\int 2Q\,\phi_i\,dV$7 is the outer-boundary term, and $-\rho\int 2Q\,\phi_i\,dV$8 is the pressure loading induced by viscous diffusion of momentum. For the moderate-to-high Reynolds number flows considered there, $-\rho\int 2Q\,\phi_i\,dV$9 for bodies with negligible rigid-body acceleration, $$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$0 for steady inflow and distant $$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$1, and $$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$2 is small compared with $$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$3 (Prakhar et al., 9 Jan 2025).

3. Partitioning workflow and feature classification

The implementation used for low-aspect-ratio wings proceeds by first solving the auxiliary Laplace problem for $$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$4 on the CFD mesh. Because $$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$5 depends only on geometry and domain, it can be precomputed once per configuration. Velocity gradients are then used to form $$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$6, $$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$7, and the scalar field

$$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$8

To distinguish orientations of vorticity, the study defines the spanwise component of the unit vorticity vector,

$$\rho\,\frac{\partial \boldsymbol{u}}{\partial t} + \rho\,(\boldsymbol{u}\cdot\nabla)\boldsymbol{u} = -\nabla p + \mu\,\nabla^2 \boldsymbol{u}, \qquad \nabla\cdot\boldsymbol{u} = 0.$$9

and partitions the flow into spanwise-oriented vorticity, defined by $\phi_i$0, and cross-span-oriented vorticity, defined as the complementary set. Each orientation class is then split into rotation-dominated and strain-dominated regions through the sign of $\phi_i$1 (Menon et al., 2022).

Spatial attribution is performed by defining cuboidal integration subdomains. In the wing application, each cuboid spans the analysis subdomain in $\phi_i$2 and $\phi_i$3 and has length $\phi_i$4 in $\phi_i$5; the near-body cuboid is centered around mid-chord, and successive cuboids tile the wake downstream in one-chord increments. Partitioned vortex-induced lift is then computed as

$\phi_i$6

for each region $\phi_i$7 and for each feature mask. The force is converted to a lift coefficient through

$\phi_i$8

The study emphasizes that additivity is guaranteed by linearity of the integral and by the use of mutually exclusive masks (Menon et al., 2022).

Moment-oriented extensions use the same architecture with a different influence potential. For spanwise pitching moment about a specified pitch axis, the potential $\phi_i$9 satisfies

$i$0

and the vorticity-induced moment is

$i$1

In this formulation the moment density is $i$2, which permits direct mapping of positive and negative moment contributions from LEV, TV, and related structures (Zhu et al., 2023).

Experimental implementations follow the same logic. In a two-dimensional pitching-wing study, the physics-based force and moment partitioning method computes the instantaneous aerodynamic moment directly from measured velocity fields without reconstructing pressure, using the same $i$3–$i$4 formulation, and supplements it with an added-mass torque obtained from body kinematics and $i$5 (Zhu et al., 2023).

4. Low-aspect-ratio wings: lift decomposition and physical mechanism

In the principal wing study, the geometry is a rectangular flat-plate wing of thickness $i$6 of chord $i$7, with aspect ratios $i$8 and $i$9, angles of attack $i=2$0 and $i=2$1, and incompressible flow at chord-based Reynolds number $i=2$2. The flow solver is a sharp-interface immersed boundary method with second-order spatial accuracy, Adams–Bashforth second-order time stepping, and multigrid solution of the pressure Poisson equation. The computational domain is $i=2$3 with a $i=2$4 grid, and the FPM analysis subdomain is $i=2$5, centered about the wing and extending $i=2$6 downstream; that subdomain captures about $i=2$7 of the vortex-induced lift (Menon et al., 2022).

Within these conditions, $i=2$8 accounts for about $i=2$9 of the total lift; viscous diffusion contributes the remainder, about $V_f$0, while viscous shear contributes a very small negative amount. The principal comparative result is that in three of four cases cross-span-oriented vorticity contributes more to lift than spanwise-oriented vorticity, contrary to the common expectation that leading-edge, spanwise vortices dominate lift (Menon et al., 2022).

Configuration	Spanwise contribution	Cross-span contribution
$V_f$1, $V_f$2	$V_f$3 of total, $V_f$4 of vortex-induced lift	$V_f$5 of total, $V_f$6 of vortex-induced lift
$V_f$7, $V_f$8	$V_f$9 of total, $\nabla^2 \phi_i = 0 \quad \text{in } V_f,$0 of vortex-induced lift	$\nabla^2 \phi_i = 0 \quad \text{in } V_f,$1 of total, $\nabla^2 \phi_i = 0 \quad \text{in } V_f,$2 of vortex-induced lift
$\nabla^2 \phi_i = 0 \quad \text{in } V_f,$3, $\nabla^2 \phi_i = 0 \quad \text{in } V_f,$4	$\nabla^2 \phi_i = 0 \quad \text{in } V_f,$5 of total, $\nabla^2 \phi_i = 0 \quad \text{in } V_f,$6 of vortex-induced lift	$\nabla^2 \phi_i = 0 \quad \text{in } V_f,$7 of total, $\nabla^2 \phi_i = 0 \quad \text{in } V_f,$8 of vortex-induced lift
$\nabla^2 \phi_i = 0 \quad \text{in } V_f,$9, $$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$0	$$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$1 of total, $$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$2 of vortex-induced lift	$$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$3 of total, $$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$4 of vortex-induced lift

The spatial partition clarifies why the balance changes between over-wing and wake regions. Over the wing, lift is largest and is dominated by spanwise-oriented vortex cores, namely the LEV and leading-edge shear layer. Cross-span vorticity in this region, primarily tip vortices, yields net negative lift due to the strain-dominated regions, consistent with classical downwash-induced lift reduction. In the wake, the total vortex-induced lift shows a non-monotonic trend with downstream distance and becomes negative near $$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$5 downstream of mid-chord. The magnitude of this negative wake lift is about $$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$6 for $$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$7, $$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$8, $$\boldsymbol{\hat{n}}\cdot\nabla \phi_i = \begin{cases} n_i & \text{on } B,\ 0 & \text{on } \Sigma, \end{cases}$$9 for $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$00, $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$01, $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$02 for $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$03, $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$04, and $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$05 for $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$06, $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$07, each expressed as a fraction of the peak over-wing lift (Menon et al., 2022).

FPM identifies the negative wake lift as a spanwise-vorticity effect. The positive contribution from spanwise cores drops rapidly due to tilting, while the negative contribution from spanwise strain regions increases because of growth of vortex-induced strain in the wake. By contrast, cross-span vorticity in the wake produces net positive lift: cross-span cores, identified with streamwise vortices and horseshoe legs, strengthen via tilting and stretching, while their strain regions grow less. The balance between rotation and strain in the two orientation classes, governed by vortex tilting and stretching, explains both why cross-span vorticity contributes more to lift in most cases and why spanwise vorticity in the near wake can produce net negative lift (Menon et al., 2022).

5. Extensions to moments, oscillations, and modal structure

The force-to-moment extension is usually termed the Force and Moment Partitioning Method (FMPM). In pitching swept wings, FMPM was used to correlate three-dimensional vortex dynamics with the resultant unsteady aerodynamic moment by solving an influence potential for pitching moment about the mid-chord, mid-span axis and integrating $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$08 over a multi-layer stereoscopic PIV volume. That study reports that LEV generally contributes positive moment, TV increasingly contributes negative moment with higher sweep, and TEV contributes negligibly because it lies near the trailing edge where $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$09. The same partitioning was used to interpret non-monotonic sweep effects on energy transfer and stability boundaries, including an optimal sweep angle of $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$10 for power extraction and destabilization, and added-mass moment peaks of about $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$11, $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$12, and $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$13 of the total peak moment for $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$14, $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$15, and $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$16, respectively (Zhu et al., 2023).

A related two-dimensional experimental study of a sinusoidally pitching NACA 0012 wing in quiescent water used the same $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$17–$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$18 formulation to estimate vortex-induced aerodynamic moment from mid-span PIV. In that case, the method partitions torque into a vorticity-induced contribution,

$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$19

an added-mass term obtained from prescribed kinematics,

$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$20

and leading- and trailing-edge subdomain contributions. The study reports that the FMPM-estimated moment follows the trend of direct measurements but underestimates the magnitude by about $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$21, with likely causes identified as near-wall PIV errors and three-dimensional effects missed by 2C2D measurements (Zhu et al., 2023).

FPM has also been applied to flow-induced vibration of cylinders. In that setting, the method partitions net hydrodynamic force into kinematic, vorticity-induced, viscous, potential, and outer-boundary components, and the vorticity-induced term is further split between a near-body shear-layer region and the wake. The principal physical conclusion is that wake vortex shedding is necessary to initiate oscillations, but it is the vorticity associated with the boundary layer over the cylinder that is responsible for the sustenance of flow-induced vibrations. Quantitatively, the cited cases show positive energy transfer from the shear layer and negative energy transfer from the wake at stationary states (Menon et al., 2020).

A more recent extension is modal force partitioning. When FPM is applied after decomposing the velocity field, the nonlinear dependence of $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$22 on velocity gradients produces inter-modal interaction terms that are difficult to interpret. The proposed remedy is to decompose the $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$23-field directly and evaluate modal forces as

$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$24

Applied to a two-dimensional cylinder, a NACA 0015 airfoil, and a three-dimensional revolving wing, this direct $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$25-field decomposition yielded more interpretable force and aeroacoustic partitions. For the revolving wing, the first six $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$26-modes capture about $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$27 of lift variance, and the three leading modes contribute about $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$28 of total acoustic intensity (Prakhar et al., 9 Jan 2025).

6. Metrics, validation, limitations, and terminological ambiguity

A data-driven variant has introduced two FPM-derived metrics for unsteady lift attribution in an impulsively started rotor blade at span-based Reynolds number $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$29: a $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$30-strength metric,

$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$31

and vortex proximity metrics,

$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$32

These satisfy

$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$33

In that rotor case, vortex-induced lift contributes about $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$34 of total pressure lift, viscous diffusion about $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$35, the acceleration reaction is non-zero only at the first time step owing to the numerical impulsive start, and Coriolis and centrifugal terms do not contribute to the pressure force when the computational domain is sufficiently large. The reported interpretation is that early lift extrema arise from the interplay between how much vortex content is present and how close positive-$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$36 cores and negative-$$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$37 coronae are to the lifting surface (Prakhar et al., 30 Jul 2025).

Across the cited literature, several robustness checks recur. In the low-aspect-ratio wing study, the sum of partitions matches the total vortex-induced lift, and vortex-induced lift plus viscous diffusion agrees with the total pressure-induced lift, with viscous shear very small and negative. In the swept-wing moment study, FMPM moment distributions qualitatively and quantitatively explain measured moment trends across sweep angles, although omitted non-vortex terms can produce magnitude mismatch. In the quiescent pitching-wing study, three phase-averaging strategies yielded similar moments, but near-wall PIV limitations remained a major source of bias (Menon et al., 2022, Zhu et al., 2023, Zhu et al., 2023).

The main limitations stated in the sources are also consistent. FPM depends on accurate velocity-gradient fields, because errors in $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$38 directly affect the load estimate. The auxiliary potential requires a Laplace solve on a sufficiently large domain with appropriate Neumann conditions. For stationary wings the body-acceleration term is zero, but for flapping or pitching wings it must be retained. Orientation thresholds such as $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$39 worked well for the reported low-aspect-ratio wing flows, but other flows may require tighter cones or modal decompositions. Finite measurement volumes in experiments can truncate far-field contributions, and two-dimensional PIV cannot capture cross-span velocity and gradient components (Menon et al., 2022, Zhu et al., 2023, Prakhar et al., 30 Jul 2025).

A recurrent source of confusion is terminological rather than conceptual. In condensed-matter phonon calculations, one paper discusses a “partition of interatomic force constants” and states explicitly that the authors do not coin the term “Force Partitioning Method (FPM).” That work concerns redistribution of supercell force constants in real space and is distinct in scope, equations, and application domain from the fluid-dynamical FPM based on $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$40 and $$C_L^\omega \approx C_L^\text{span} + C_L^\text{cross} + C_L^\text{outside}$$41 (Lee et al., 2020).