---
title: 'LinDistFlow: Linearized Radial Feeder Model'
url: https://www.emergentmind.com/topics/lindistflow-model
type: topic
---

# LinDistFlow: Linearized Radial Feeder Model

LinDistFlow is a linearized, lossless approximation of the AC power flow equations tailored for radial distribution systems. It originates from the nonlinear DistFlow model, exploits the typical radial (tree) topology of distribution systems, and uses squared voltage magnitudes directly in the formulation, which makes voltage drops explicit and preserves linear dependence on active and reactive power variables. In the contemporary literature, it serves both as a foundational linear power flow model for feeder analysis and optimization and as a reference point for multiphase generalizations, uncertainty-aware formulations, loss-compensated variants, and parameter-optimized surrogates [2104.02118][2404.05125].

## 1. Canonical formulation

In radial distribution networks, LinDistFlow is commonly written as a simplified branch-flow model in which active and reactive power balances are enforced along the feeder and squared-voltage drops are expressed linearly in line flows. One standard statement is
\[
\begin{align}
(d_i^p - g_i^p) + \sum_{j \in \mathcal{C}_i} f_j^p &= f_{i}^p, \quad p \in \{P,Q\}, \quad \forall i \in \mathcal{N}^+ \\
u_{A_i} - 2(f_i^P R_i + f_i^Q X_i) &= u_i, \quad \forall i \in \mathcal{N}^+ ,
\end{align}
\]
with \(u_i = v_i^2\), line resistances \(R_i\), reactances \(X_i\), and parent/child relations defined by the feeder tree [1803.04912]. Equivalent branch-oriented forms appear throughout the literature, for example
\[
\begin{aligned}
p_{ij} &= p_{j} + \sum_{k:(j,k) \in \mathcal{L}}p_{jk}, \\
q_{ij} &= q_{j}+\sum_{k:(j,k)\in \mathcal{L}}q_{jk}, \\
w_i - w_j &= 2(r_{ij} p_{ij} + x_{ij} q_{ij}),
\end{aligned}
\]
where \(w_i=v_i^2\) and \((i,j)\) is a radial line [2009.02088].

A compact sensitivity form is also standard:
\[
\mathbf{v} = R \mathbf{p} + X \mathbf{q} + \tilde{\mathbf{v}},
\]
where \(\mathbf{v}\) collects squared voltage magnitudes, \(\mathbf{p}\) and \(\mathbf{q}\) collect active and reactive injections, and the entries of \(R\) and \(X\) are sums of line resistances and reactances over the common path from the slack bus to the buses of interest. For single-phase radial networks,
\[
R_{ij} = \sum_{(\zeta, \xi) \in \mathcal{E}_i \cap \mathcal{E}_j} 2r_{\zeta\xi},\quad
X_{ij} = \sum_{(\zeta, \xi) \in \mathcal{E}_i \cap \mathcal{E}_j} 2x_{\zeta\xi},
\]
with \(\mathcal{E}_i\) denoting the unique path from the substation to node \(i\) [1903.00072].

These forms encode the same modeling idea: branch losses and current-squared terms are omitted, while voltage drops are retained through squared-voltage variables. That representation is what makes LinDistFlow especially useful in OPF, sensitivity analysis, and convex relaxations built for radial feeders [2305.00400].

## 2. Assumptions and modeling regime

The key assumption underpinning LinDistFlow is the neglect of quadratic loss terms, specifically terms proportional to the square of line current. In the nonlinear DistFlow model, terms such as \(r_{ij}\ell_{ij}\), \(x_{ij}\ell_{ij}\), and \((r_{ij}^2+x_{ij}^2)\ell_{ij}\) appear in branch power balance and voltage equations; LinDistFlow drops these terms to obtain a linear approximation [2508.04382][2404.05125].

The model therefore assumes radial topology, modest system losses, and explicit use of squared voltage magnitudes. This is why it differs structurally from classic DC power flow. Unlike the DC power flow model, LinDistFlow is tailored to feeders with higher \(R/X\) ratios, unbalanced conditions, and significant voltage deviations, and it explicitly incorporates squared voltage magnitudes and their evolution along branches [2508.04382]. In many formulations it is also described as a lossless approximation of the AC power flow equations for radial networks, with computational tractability as a primary motivation [1803.04912].

These assumptions delimit its scope. LinDistFlow is most naturally associated with radial, single-phase or balanced equivalent feeders, although later work extends it to unbalanced three-phase systems. It remains attractive because, once losses are ignored, voltage, flow, and sometimes uncertainty propagation constraints become linear or convex. A plausible implication is that LinDistFlow is less a single equation than a modeling family organized around three recurring choices: radial topology, squared-voltage variables, and omission of quadratic loss terms.

## 3. Multiphase and generalized variants

Although the classical model is tied to radial single-phase feeders, later work extends its structure in two directions: multiphase unbalanced modeling and topological generalization. For multi-phase radial networks, the linearization becomes
\[
\mathbf{v}_{\Xi} = R_{\Xi} \mathbf{p}_{\Xi} + X_{\Xi} \mathbf{q}_{\Xi} + \tilde{\mathbf{v}}_{\Xi},
\]
with phase-indexed voltage and injection vectors and sensitivity entries
\[
\frac{\partial v_i^{\varphi}}{\partial p_j^{\phi}} =
2\,\mathrm{Re} \left\{ \overline{Z}_{ij}^{\varphi\phi} \cdot \omega^{\varphi-\phi} \right\}, \quad
\frac{\partial v_i^{\varphi}}{\partial q_j^{\phi}} =
-2\,\mathrm{Im} \left\{ \overline{Z}_{ij}^{\varphi\phi} \cdot \omega^{\varphi-\phi} \right\},
\]
where \(\omega = e^{-j 2\pi/3}\) and \(\overline{Z}_{ij}^{\varphi\phi}\) aggregates mutual impedances over common paths [1903.00072].

A more explicit generalization is the Generalized LinDistFlow (GLDF) model:
\[
\boxed{ \bm{v} = \mathrm{diag}(\bar{\bm{E}})\bm{E} + \mathbf{M}\bm{p} + \mathbf{N}\bm{q} + \mathbf{\Lambda} }
\]
with open-circuit voltage \(\bm{E}:=-\mathbf{Y}_{LL}^{-1}\mathbf{Y}_{LS}\bm{V}_S\), \(\mathbf{Z}=\mathbf{Y}_{LL}^{-1}\), and
\[
\mathbf{M} = 2 \mathrm{Re}\left(\mathrm{diag}(\bm{E}) \overline{\mathbf{Z} \cdot \mathrm{diag}(\bm{E})^{-1}} \right), \quad
\mathbf{N} = -2 \mathrm{Im}\left(\mathrm{diag}(\bm{E}) \overline{\mathbf{Z} \cdot \mathrm{diag}(\bm{E})^{-1}} \right).
\]
The offset \(\mathbf{\Lambda}\) enforces exactness at the linearization point \((\bm{S}^*,\bm{V}^*)\) [2104.02118].

This generalized model handles multiphase systems, applies to generic topologies including radial and meshed networks, allows linearization at arbitrary operating points, and reduces to the classic or multiphase LinDistFlow model for a radial network at zero injection [2104.02118]. When both GLDF and fixed-point linearization are linearized at zero injection, GLDF always yields lower errors in squared voltage magnitude than FPL, for any loading, in radial networks; the paper states this as
\[
\left| \hat{v}_i - v_i \right| \leq \left| \tilde{v}_i - v_i \right|,\quad \forall i.
\]
A related formulation appears in uncertainty analysis, where LinDistFlow is treated as equivalent to the Linear Coupled Power Flow model under radial network conditions:
\[
\begin{bmatrix} p\\ q \end{bmatrix}
=
\begin{bmatrix} G & -B\\ -B & -G \end{bmatrix}
\begin{bmatrix} \delta v\\ \delta \theta \end{bmatrix},
\]
with a closed-form inverse on trees involving effective resistances and reactances [2510.17798].

## 4. Role in optimization, control, and market design

LinDistFlow is central to OPF formulations because it converts AC-feeder constraints into linear or convex constraints that remain sensitive to voltage and reactive power. In data-driven distributionally robust OPF for radial systems, the model enables explicit propagation of nodal injection uncertainty into random voltages and flows, so that chance constraints can be recast as second-order conic constraints; the resulting approach solved systems up to the IEEE 8500-bus network in about \(24\) seconds [1803.04912]. In chance-constrained flexibility-request design for local flexibility markets, LinDistFlow likewise supports SOCP reformulations and yields suboptimality gaps of \(3.2\%\) on a 15-bus case and \(8.25\%\) on an 81-bus real network relative to a stochastic market-clearing benchmark [2110.05983].

The model also underpins distributionally robust formulations that account for renewable dependence. In a LinDistFlow-based radial OPF with a copula-based ambiguity set, a conic reformulation preserves tractability while embedding wind-farm dependencies, and reported cost savings reach up to \(12\%\) relative to traditional DRO [2104.08101]. In local electricity market design, LinDistFlow equations are combined with intrusive generalized Polynomial Chaos to propagate uncertainty through the grid model and produce probabilistic locational marginal prices in a two-stage convex formulation [2510.12318].

In control and distributed computation, LinDistFlow’s sensitivity structure is used to decompose large OPF problems hierarchically. For multi-phase radial networks, the block structure of \(R_{\Xi}\) and \(X_{\Xi}\) enables a hierarchical distributed primal-dual algorithm in which regional coordinators and a central coordinator exchange aggregated quantities rather than full-network sensitivities. On a 4,521-node feeder, the algorithm achieves more than 10-fold acceleration in the speed of convergence compared to a centrally coordinated implementation [1903.00072]. A related hierarchical formulation based on autonomous-grid partitions reports more than tenfold improvement in convergence speed on the same scale of feeder [1809.08624].

LinDistFlow is also used in interface flexibility and topology optimization. For TSO-DSO flexibility domain identification on the IEEE 33-bus radial system, the LinDistFlow formulation is an LP solved with Gurobi and required \(2.3\) seconds for 200 boundary points, compared with \(1798.3\) seconds for AC-OPF and \(666.2\) seconds for DistFlow [2009.02088]. In daily topology reconfiguration for improving fairness in photovoltaic curtailment, day-ahead grid constraints are modeled using LinDistFlow, while the real-time stage uses a first-order Taylor linearization of AC power flow [2403.07853]. In optimization-based control of distributed battery storage, a LinDistFlow-based convex multi-period formulation yields reductions of \(11\) per cent in losses and \(26\) per cent in peak substation power compared to other state-of-the-art algorithms [2409.00424].

A further strand of work uses LinDistFlow for economic interpretation. In single-phase radial LDF-OPF, closed-form marginal values of real and reactive demand are derived from dual variables, and upper bounds are obtained for the change in marginal demand prices when apparent power flow limits become binding [2305.00400].

## 5. Accuracy, limitations, and failure modes

The principal limitation of LinDistFlow is the omission of losses. Because the model ignores quadratic loss terms, all flows are treated as lossless; the only effect of line parameters is on voltage drops [2508.04382]. This makes the approximation most accurate when system losses are modest, but it also creates systematic distortions in operating regions where losses materially affect interface exchange, feasible domains, or voltage profiles.

One documented consequence is that LinDistFlow overestimates capability at the point of common coupling. In flexibility aggregation, the lossless model omits the accumulated line-loss terms \(R^\top L\) and \(X^\top L\), so a flexibility region computed via LinDistFlow overestimates the feasible set at the PCC because, in reality, more power needs to be injected to cover losses [2505.01715]. In TSO-DSO flexibility domain identification, the main discrepancies appear in the negative orthant, because LinDistFlow does not consider losses and therefore overestimates the capability of energy that could be generated from the active distribution network [2009.02088].

A second failure mode is temporal error accumulation. Simulations on the KIT Campus Nord network with real demand and solar data show that line losses are generally underestimated by linear models, and line-loss errors tend to accumulate both at the point of common coupling and over extended time horizons [2508.04382]. In day-ahead scheduling, this can manifest as scheduled ESS state of charge drifting below its intended value because of accumulated undercounting of system losses [2508.04382].

A third issue concerns robustness away from the linearization regime. This is one reason later work compares LinDistFlow with generalized linearizations around nonzero operating points. GLDF and FPL are both exact at the linearization point, but GLDF is reported to be more robust and to exhibit more consistent performance as operating points depart from that point [2104.02118].

Recent work also studies LinDistFlow under uncertain network parameters. Using matrix concentration inequalities, the expected operator norm error of the LinDistFlow or LCPF matrix under independent, bounded uncertainties scales as \(O(\Delta \sqrt{n \log n})\), and the probability of large deviation decays exponentially in \(t^2\) [2510.17798]. This suggests that the model’s approximation error can be bounded probabilistically under bounded and independent line uncertainty, rather than treated only as a deterministic modeling residual.

## 6. Compensation, calibration, and current directions

Several recent directions aim not to replace LinDistFlow, but to preserve its structure while correcting its dominant errors. One approach is explicit loss compensation. In flexibility aggregation for integrated transmission-distribution systems, a system loss compensation method first computes the LinDistFlow-based flexibility set and then corrects each candidate interface point by adding estimated active and reactive losses written as quadratic functions of the PCC exchange:
\[
p^\mathrm{loss}(u) = \frac{1}{2} u^\top H_p u + g_p^\top u + c_p,\qquad
q^\mathrm{loss}(u) = \frac{1}{2} u^\top H_q u + g_q^\top u + c_q.
\]
The compensated set
\[
\mathcal{U}_{\text{SLC}}
=
\left\{
u ~\middle|~
u = u_{\mathrm{LDS}} +
\begin{bmatrix}
p^\mathrm{loss}(u_{\mathrm{LDS}})\\
q^\mathrm{loss}(u_{\mathrm{LDS}})
\end{bmatrix},
~~
u_{\mathrm{LDS}} \in \mathcal{U}_{\mathrm{LDS}}
\right\}
\]
is reported to closely approximate the actual feasible region while preserving data privacy [2505.01715].

Another approach is parameter optimization. Rather than fixing LinDistFlow coefficients at physical line parameters, the Optimized LinDistFlow formulation introduces trainable coefficient and bias parameters:
\[
\mathbf{v} = v_0 \mathbf{1}
+ 2\mathbf{A}^{-1} \mathbf{D}_r \mathbf{A}^{-T} (\mathbf{p} + \boldsymbol{\rho})
+ 2\mathbf{A}^{-1} \mathbf{D}_x \mathbf{A}^{-T} (\mathbf{q} + \boldsymbol{\varrho})
+ \boldsymbol{\gamma}.
\]
The parameters are trained offline with sensitivity information using the Truncated Newton Conjugate-Gradient method to minimize discrepancies relative to the nonlinear DistFlow model [2404.05125]. Reported improvements reach up to \(92\%\) in \(L_{1}\)-norm loss and \(88\%\) in \(L_{\infty}\)-norm loss relative to traditional LinDistFlow, and the optimized approximation is also assessed under topology changes and in hosting-capacity optimization [2404.05125].

LinDistFlow has also been adapted to settings where the load model itself is voltage-dependent. For load-altering attacks with ZIP loads, the combination of LinDistFlow with a ZP approximation yields a matrix-based closed-form expression for squared voltages,
\[
U_{(N-1)\times 1} = (I_{(N-1)\times (N-1)} - \Omega'_{(N-1)\times (N-1)})^{-1} \Omega''_{(N-1)\times 1},
\]
which makes analytical attack-impact calculations possible even though the original ZIP model is nonlinear in bus voltage [2311.06429].

Taken together, these developments indicate that LinDistFlow remains a reference model not because it is exact, but because its linear, squared-voltage structure is unusually adaptable. This suggests that the most active line of research is no longer the use of an unmodified lossless approximation alone, but the systematic refinement of LinDistFlow through generalized operating-point selection, loss compensation, uncertainty quantification, and parameter calibration [2104.02118][2404.05125].

Source: https://www.emergentmind.com/topics/lindistflow-model