---
title: Minimum Flow Problem (MFP) Overview
url: https://www.emergentmind.com/topics/minimum-flow-problem-mfp
type: topic
---

# Minimum Flow Problem (MFP) Overview

Searching arXiv for relevant papers on the Minimum Flow Problem, including both continuum/PDE and graph/min-cost flow interpretations.
arxiv_search(query="minimum flow problem Beckmann pedestrian flow minimum flow problem", max_results=5, sort_by="relevance")
arxiv_search(query="minimum cost flow interior point Laplacian generalized minimum cost flow", max_results=10, sort_by="relevance")
The Minimum Flow Problem (MFP) denotes distinct but related optimization problems whose common structure is the determination of an admissible flow under conservation laws, bounds, and an objective functional. In the literature considered here, one formulation is a Beckmann-type continuous flow problem used as the correction phase of a prediction–correction model for pedestrian motion, where a density \(\rho\in[0,1]\) is transported by a desired velocity and then corrected by an optimization problem that enforces admissibility and boundary flux conditions [2302.11315]. Another formulation interprets MFP as the problem of finding the smallest feasible \(s\)–\(t\) flow, or the minimum total flow subject to feasibility, and treats it as a special case of Minimum Cost Flow (MCF), so that modern interior-point methods for MCF apply directly to MFP after a standard reduction [1312.6713]. The term is therefore not attached to a single canonical mathematical object; its meaning depends on whether the ambient setting is a continuum PDE model or a bounded flow problem on a directed graph.

## 1. Terminological scope and problem classes

A common source of ambiguity is that “Minimum Flow Problem” is used in more than one precise sense. In the continuous crowd-motion setting, the problem is posed on a bounded domain \(\Omega\subset\mathbb{R}^N\) with \(N=2\), with boundary split as \(\partial\Omega=\Gamma_N\cup\Gamma_D\), and the unknowns are a density field and a correction flow field. The objective is not merely to minimize a graph-theoretic path cost, but to compute a flux \(\Phi\) and a corrected density \(\rho\) that satisfy a divergence constraint, capacity bounds \(0\le \rho\le 1\), and mixed boundary conditions, while optimizing a convex cost with a boundary gain term [2302.11315].

In graph optimization, by contrast, MFP is often used for the variant where one seeks the smallest \(s\)–\(t\) flow that meets edge lower bounds, or the minimum total amount of flow subject to feasibility. In that interpretation, the problem is a special case of Minimum Cost Flow: one introduces a unit cost on a single “value” edge or otherwise encodes flow value as a linear cost, and minimizing total cost coincides with minimizing the flow value [1312.6713].

These two usages share a variational backbone. Both impose conservation constraints, both admit dual variables with complementary slackness or complementarity structure, and both use convex optimization techniques. They differ, however, in the geometry of the state space, the meaning of the flow variable, the type of admissibility constraint, and the numerical linear algebra required for solution.

## 2. Continuum MFP in prediction–correction pedestrian flow

In the pedestrian-flow formulation, \(\rho=\rho(t,x)\in[0,1]\) is the macroscopic pedestrian density, \(V\) is the spontaneous velocity driving the prediction step, and \(W\) is the “patch” velocity generated by the MFP in the correction step to resolve congestion and enforce admissibility [2302.11315]. The overall model is a prediction–correction scheme.

The prediction phase evolves the density over a short sub-step \([t_k,t_{k+1/2}]\) by the continuity equation
\[
\partial_t \rho + \nabla\cdot(\rho V)=0
\quad \text{in } (t_k,t_{k+1/2})\times\Omega,
\qquad
\rho(t_k,\cdot)=\rho_k.
\]
The desired velocity is obtained from an eikonal potential \(\varphi\), interpreted as expected travel time to \(\Gamma_D\), through
\[
\|\nabla\varphi(x)\|=f(x)\quad \text{in }\Omega,
\qquad
\varphi=0\quad \text{on }\Gamma_D,
\qquad
V(x)=-\frac{\nabla\varphi(x)}{\|\nabla\varphi(x)\|}.
\]
Here \(f(x)>0\) is a given continuous “spontaneous speed,” possibly space-dependent or congestion-aware in variants.

Given the predicted density \(\tilde\rho:=\rho(t_{k+1/2},\cdot)\), the correction phase computes a corrected density and a correction flow by solving the MFP
\[
\inf_{(\Phi,\rho)}\left\{
\int_\Omega F(x,\Phi(x))\,dx
-
\int_{\Gamma_D} g(x)\,\Phi\cdot n\,ds
\;:\;
0\le \rho\le 1,\;
-\nabla\cdot\Phi=\tilde\rho-\rho\ \text{in }\Omega,\;
\Phi\cdot n=\eta\ \text{on }\Gamma_N
\right\}.
\]
In the time-discrete form over a step of size \(\tau>0\),
\[
-\tau\,\nabla\cdot\Phi=\rho_{k+1/2}-\rho
\quad \text{in }\Omega,
\qquad
\Phi\cdot n=\eta\quad \text{on }\Gamma_N.
\]

Two prototypical choices of the flow cost \(F(x,\xi)\) are emphasized. The quadratic choice
\[
F(x,\xi)=\tfrac12\|\xi\|^2
\]
is described as Wasserstein-2 type and connects to Maury–Roudneff–Santambrogio’s \(W_2\)-gradient flow correction. The \(1\)-homogeneous choice
\[
F(x,\xi)=k(x)\,\|\xi\|
\]
with \(k(x)\ge 0\) continuous is granular or Wasserstein-1 type, induces sandpile-like dynamics, and links to \(W_1\)-gradient flows.

The boundary data are integral to the formulation. The term
\[
-\int_{\Gamma_D} g(x)\,\Phi\cdot n\,ds
\]
models door-exit “gain,” with \(g(x)\ge 0\) rewarding flow crossing \(\Gamma_D\). The normal trace on \(\Gamma_D\) is free and optimized, whereas on \(\Gamma_N\) it is prescribed by \(\eta\), which acts as an imposed inflow rate at entrances.

## 3. Duality, complementarity, and PDE structure

For the granular case \(F(x,\xi)=k(x)\|\xi\|\), the correction phase has an explicit dual structure [2302.11315]. Under the technical assumption that there exists \(\tilde g\in W^{1,\infty}(\Omega)\) such that \(\tilde g=g\) on \(\Gamma_D\) and \(\|\nabla\tilde g(x)\|\le k(x)\) a.e., and \(\eta\in W^{-1/s,s}(\Gamma_N)\) for some \(s>1\), the admissible dual potentials are
\[
\mathcal{K}_k
:=
\left\{
p\in W^{1,\infty}(\Omega)
:\;
p=g\ \text{on }\Gamma_D,\;
\|\nabla p(x)\|\le k(x)\ \text{a.e. in }\Omega
\right\}.
\]
The primal minimum equals the dual maximum:
\[
\inf_{(\Phi,\rho)}\Big\{
\int_\Omega k\|\Phi\|\,dx
-
\int_{\Gamma_D} g\,\Phi\cdot n\,ds
\Big\}
=
\max_{p\in\mathcal{K}_k}
\left\{
\int_\Omega \tilde\rho\,p\,dx
+
\tau\int_{\Gamma_N}\eta\,p\,ds
-
\int_\Omega p^{+}\,dx
\right\},
\]
where \(p^{+}=\max(p,0)\).

The hard density constraint enters the dual through the sign-graph relation
\[
\rho\in \mathrm{Sign}^{+}(p)
\quad\Longleftrightarrow\quad
0\le \rho\le 1,\qquad p(1-\rho)=0.
\]
Equivalently, \(p>0\) implies \(\rho=1\), \(p<0\) implies \(\rho=0\), and \(p=0\) allows any \(\rho\in[0,1]\). In the KKT-like optimality conditions, introducing the patch \(W\) and a Lagrange multiplier \(m\ge 0\) for the slope constraint yields
\[
W=m\nabla p,\qquad
p\ge 0,\qquad
\|\nabla p\|\le k,\qquad
m(\|\nabla p\|-k)=0,
\]
together with
\[
0\le \rho\le 1,\qquad
p(\rho-1)=0.
\]
The last relation is the capacity complementarity condition, and the paper states that \(W\) acts only in saturated zones.

The time-discrete correction step can then be written as
\[
\rho_{k+1}-\tau\,\nabla\cdot\big(\rho_{k+1}W\big)=\rho_{k+1/2}
\quad \text{in }\Omega,
\]
with \(p=g\) on \(\Gamma_D\) and \(\Phi\cdot n=\eta\) on \(\Gamma_N\). As \(\tau\to 0\), one formally obtains the continuous-time limiting evolution
\[
\partial_t \rho + \nabla\cdot\big(\rho\,(V-W)\big)=0,
\qquad
W=m\nabla p,
\]
with
\[
p\ge 0,\quad
0\le\rho\le 1,\quad
p(\rho-1)=0,\quad
\|\nabla p\|\le k,\quad
m(\|\nabla p\|-k)=0.
\]

The quadratic case has a different but simpler optimality structure. When
\[
F(x,\xi)=\tfrac12\|\xi\|^2,
\]
the optimality conditions reduce to
\[
\Phi=\nabla p,\qquad
\rho\in\mathrm{Sign}^{+}(p),
\]
and the evolution becomes
\[
\frac{\partial \rho}{\partial t}+\nabla\cdot\big(\rho(V-\nabla p)\big)=0,
\qquad
p\ge 0,\quad 0\le\rho\le 1,\quad p(\rho-1)=0.
\]
The cited work identifies this with the well-known \(W_2\)-projection approach, whereas the \(1\)-homogeneous case aligns with \(W_1\)-gradient flows and granular dynamics. This suggests that the continuum MFP interpolates between a diffusion-like pressure projection and a sandpile-type redistribution mechanism, depending on the choice of \(F\).

## 4. Numerical realization and observed behaviors

The prediction equation is discretized by a finite-volume method on an \(m\times n\) Cartesian grid of cell size \(h\) [2302.11315]. With explicit Euler in time and upwind fluxes, the transport scheme is
\[
\frac{\rho_{i,j}^{k+1/2}-\rho_{i,j}^k}{\tau}
+\frac{1}{h}\left[\rho_{i,j}^k V^x_{i+1/2,j}-\rho_{i-1,j}^k V^x_{i-1/2,j}\right]
+\frac{1}{h}\left[\rho_{i,j}^k V^y_{i,j+1/2}-\rho_{i,j-1}^k V^y_{i,j-1/2}\right]
=0.
\]
The CFL condition is
\[
\max_{i,j}\|V_{i,j}\|\,\tau/h<1/2.
\]

The eikonal problem is computed through a convex variational problem of the form
\[
\min_u \mathcal{F}_h(u)+\mathcal{G}_h(\nabla_h u),
\]
with \(\mathcal{F}_h(u)=-h^2\sum_{i,j}u_{i,j}\) when \(u=0\) on \(\Gamma_D\), and \(\mathcal{G}_h=\mathbb{I}_{B(0,f)}\). Chambolle–Pock iterations are then used, with proximal maps given by projection onto the ball of radius \(f\) for \(\nabla_h u\) and enforcement of \(u=0\) on \(\Gamma_D\).

For the correction phase, the discrete primal problem is
\[
\min_{(\rho,\Phi)}\;
\mathcal{F}_h(\rho,\Phi)+\mathcal{G}_h(\Lambda_h(\rho,\Phi)),
\]
where
\[
\mathcal{F}_h(\rho,\Phi)
=
h^2\sum_{i,j}\tau\,k_{i,j}\|\Phi_{i,j}\|
+
I_{[0,1]}(\rho),
\qquad
\Lambda_h(\rho,\Phi)=\rho-\tau\,\mathrm{div}_h\Phi.
\]
The Chambolle–Pock updates take the form
\[
\rho^{\ell+1}=\mathrm{Proj}_{[0,1]}\big(\rho^\ell-\beta\,\bar p^\ell\big),
\qquad
\Phi^{\ell+1}=\mathrm{shrink}\big(\Phi^\ell-\beta\,\nabla_h\bar p^\ell;\,\beta\tau k\big),
\]
\[
v^{\ell+1}=p^\ell+\alpha\rho^{\ell+1}-\alpha\,\mathrm{div}_h(\Phi^{\ell+1}),
\qquad
p^{\ell+1}=v^{\ell+1}-\alpha\,\mathrm{Proj}_{\{\|\nabla_h p\|\le k\}}\big(v^{\ell+1}/\alpha\big),
\]
\[
\bar p^{\ell+1}=2p^{\ell+1}-p^\ell,
\qquad
\alpha\beta\|\Lambda_h\|^2<1,
\qquad
\theta=1.
\]
For \(\theta=1\) and \(\alpha\beta\|\Lambda_h\|^2<1\), the iterates converge to a saddle point.

The numerical experiments serve primarily to identify qualitative regimes. In a one-room evacuation with constant \(f\), density moves along \(V\) towards the door, correction prevents exceeding \(\rho=1\), and evacuation proceeds smoothly. With a non-constant spontaneous speed
\[
f(x)=\exp\!\big(-3[(x-1/2)^2+(y-1/2)^2]\big),
\]
the crowd avoids the central slow bump, deflecting around it. With
\[
f(x)=|\cos(3x+5y)|+0.2
\]
and a source inflow on the left boundary, successive regions of small \(f\) produce successive congestion zones, and the system reaches equilibrium when inflow equals outflow.

The comparison between homogeneous \(W_1\)/granular correction and quadratic \(W_2\)/Laplacian correction uses the average density at exits over time and the norms \(\|\rho_1-\rho_2\|_{L^\infty}\) and \(\|\rho_1-\rho_2\|_{L^2}\). Both models exhibit qualitatively similar patterns, but the \(W_1\)/granular correction often evacuates faster. In a configuration with a rectangular obstacle near the door, adding an obstacle slows down evacuation and may trap pedestrians, contrasting with some microscopic findings where an obstacle can regulate flow.

The same source also emphasizes limitations. The continuous-time granular model is challenging, uniqueness is nontrivial and an open topic, and in the homogeneous case the optimal flux may be a vector-valued Radon measure rather than \(L^2\). Extensions explicitly listed include heterogeneous agents, anisotropy, time-dependent \(f\) or door policies \(g(t,x)\), feedback coupling \(f=f(p)\), and multiscale interactions.

## 5. Graph-theoretic MFP as a minimum cost flow problem

In the graph setting, the standard directed network is \(G=(V,E)\) with \(|V|=n\), \(|E|=m\), node–edge incidence matrix \(A\in\mathbb{R}^{n\times m}\), node supplies and demands \(b\in\mathbb{R}^n\), edge costs \(c\in\mathbb{R}^m\), lower bounds \(\ell\in\mathbb{R}^m\), and upper bounds \(u\in\mathbb{R}^m\) [1312.6713]. The primal MCF formulation is
\[
\text{minimize } c^T f
\qquad
\text{subject to } Af=b,
\qquad
\ell\le f\le u.
\]
The assumptions stated are \(\ell_e\le u_e\) for all \(e\), nonempty interior, and finite bounds. The paper notes that \(A\) has rank \(n-1\) in a connected graph.

The dual introduces node potentials \(\pi\in\mathbb{R}^n\) and edge duals \(s^+,s^-\in\mathbb{R}^m\):
\[
\text{maximize } \pi^T b + s^{-T}\ell - s^{+T}u
\]
subject to
\[
A^T\pi + s^- - s^+ + c = 0,
\qquad
s^-\ge 0,\quad s^+\ge 0.
\]
The reduced costs are
\[
r_e:=c_e+\pi_{\mathrm{tail}(e)}-\pi_{\mathrm{head}(e)}=s^+_e-s^-_e.
\]
The KKT conditions consist of primal feasibility, dual feasibility, and complementary slackness
\[
s^-_e(f_e-\ell_e)=0,
\qquad
s^+_e(u_e-f_e)=0
\quad \text{for all } e.
\]

Within this framework, MFP is obtained as a special case of MCF. One reduction is to introduce a circulation formulation with an added “value” edge \(e^*\) from \(t\) to \(s\), with cost \(c_{e^*}=1\), capacity \(u_{e^*}=\infty\), and \(b\equiv 0\). Any circulation sending \(F\) units from \(s\) to \(t\) through \(e^*\) incurs cost \(F\), so minimizing cost is the same as minimizing flow value. Another reduction sets all costs to \(0\) except on arcs leaving \(s\) or via an auxiliary variable. In both cases, the LP remains of the form \(\min c^T x\) subject to flow conservation and bound constraints, so the same MCF algorithms apply without modification.

The parameter \(U\) is the “width” or magnitude parameter entering logarithmic factors; in the paper’s tables, \(U\) denotes “the maximum absolute value of capacities and costs.” Costs can be negative, lower bounds are fully supported, and the method is weakly polynomial.

## 6. Interior-point algorithms, complexity bounds, and comparative perspective

The graph-theoretic paper develops a new interior-point method and applies it to maximum flow, minimum cost flow, and lossy generalized minimum cost flow [1312.6713]. For MCF, and thus for MFP via reduction, the stated running time is
\[
\tilde O\!\left(m\sqrt n\,\log^{O(1)}(U/\epsilon)\right)
\]
work and
\[
O\!\left(\sqrt n\,\log^{O(1)}(U/\epsilon)\right)
\]
parallel depth.

The method uses coordinate-wise barriers on the interval \((\ell_e,u_e)\). If one bound is finite, the barrier is logarithmic; if both are finite, the barrier is
\[
\phi_e(x)=-\log\cos(a_e x+b_e),
\]
with parameters chosen so that \(\phi_e\to+\infty\) at both endpoints and \(\phi_e\) is \(1\)-self-concordant. For weights \(w\in\mathbb{R}_{>0}^m\) and parameter \(t>0\), the penalized objective is
\[
f_t(x,w):=t\,c^T x+\sum_e w_e\phi_e(x_e),
\]
and for fixed \(w\), the path
\[
x_w(t):=\arg\min_{A^T x=b} f_t(x,w)
\]
connects an interior weighted center to the LP optimum as \(t\to\infty\).

The projected Newton step is defined using the diagonal Hessian \(D(x)=\mathrm{diag}(\phi''(x))\) and the \(D\)-weighted orthogonal projector onto \(\ker(A^T)\). A mixed centrality measure
\[
\delta_t(x,w)
\]
combines the Hessian norm and coordinatewise multiplicative change. The stated novelty is that the algorithm adapts the weights during the run and thereby achieves \(O(\sqrt n\cdot \mathrm{polylog})\) iterations, which is strictly better than the \(m\)-based bound suggested by the best general barrier for box constraints.

The weight function is itself defined by a convex program in \(w\). The resulting weights satisfy
\[
\|w\|_1=\mathrm{rank}(A)+\beta m\approx \Theta(n),
\qquad
\max_e w_e\le 1+\beta\le 2.
\]
To avoid the expense of computing \(w(x)\) exactly at every step, the algorithm maintains \(w\) close to \(w(x)\) by a “chasing 0” game on \(z=\log(w)-\log(w(x))\), using a potential \(\Phi_\mu(z)=\sum(e^{\mu z_i}+e^{-\mu z_i})\).

When specialized to MCF, the fundamental linear system is
\[
L:=A D(x)^{-1} A^T,
\]
which is a weighted graph Laplacian. Each centering iteration therefore reduces to \(O(1)\) solves in Laplacian-like systems, together with weight updates and normal-force stabilization. Using Spielman–Peng nearly-linear Laplacian solvers, the method becomes parallelizable with the depth and work bounds above.

The comparative claims are explicit. Goldberg–Rao achieved
\[
O\!\left(m\min\left(n^{2/3},m^{1/2}\right)\log\left(n^{2}/m\right)\log U\right)
\]
for exact maximum flow on directed graphs with integer capacities. Even–Tarjan and Karzanov achieved
\[
O(\min\{m^{3/2},mn^{2/3}\})
\]
for dense unit-capacity directed graphs. Mądry obtained
\[
O(m^{10/7})
\]
for maximum flow on uncapacitated directed graphs. Daitch–Spielman obtained
\[
\tilde O(m^{3/2})
\]
for MCF and generalized MCF using dual interior-point methods with Laplacian and M-matrix solvers. The new algorithm improves the iteration count to \(\tilde O(\sqrt n)\) and thus the total work to \(\tilde O(m\sqrt n)\).

The stopping criterion is expressed through a duality gap bound:
\[
c^T x_t(w)-\mathrm{OPT}\le \|w\|_1/t.
\]
For standard MCF with integral inputs, taking \(\epsilon\approx 1/\mathrm{poly}(mU)\) and rounding yields an exact integral solution. For generalized lossy MCF, the result is \(\epsilon\)-approximate. If MFP is defined instead as minimizing some other linear functional of the flow under the same conservation and bound constraints, the method still applies: one replaces \(c\) by the appropriate vector and retains the same algorithmic framework.

The two literatures thus present MFP as a family of optimization problems rather than a single formalism. In one branch, MFP is a continuum correction mechanism coupling eikonal-driven transport to a convex admissibility projection with granular or diffusive character. In the other, it is a bounded linear flow problem reducible to minimum cost flow and solvable by interior-point methods with Laplacian-based linear algebra. The shared conceptual core is flow selection under convex constraints; the differences lie in whether the state variable is a density in a domain or a vector of edge flows in a graph.

Source: https://www.emergentmind.com/topics/minimum-flow-problem-mfp