---
title: Localized Reduced Basis Methods
url: https://www.emergentmind.com/topics/localized-reduced-basis-methods
type: topic
---

# Localized Reduced Basis Methods

Searching arXiv for recent and foundational papers on localized reduced basis methods to ground the article in current literature.
Localized reduced basis methods are reduced basis or localized model order reduction techniques in which approximation is built from local spaces rather than a single global reduced space. In the literature, locality may refer to overlapping or non-overlapping subdomains in physical space, coarse elements or multiscale patches, or neighborhoods in parameter or solution space selected by clustering or anisotropic distance. The shared objective is to replace a global reduced space that becomes inefficient for multiscale coefficients, high contrast, moving discontinuities, or parameterized geometries by a family of smaller spaces that can be combined, selected, or enriched adaptively [1710.02104] [1410.3253] [1204.3846].

## 1. Forms of locality

A standard reduced basis formulation starts from a parametrized variational problem
\[
a(u(\mu),v;\mu)=f(v;\mu)\qquad \forall v\in V,
\]
and replaces the high-dimensional solution by a Galerkin approximation in a low-dimensional space. Localized variants keep this projection structure but alter the architecture of the reduced space. Instead of one global space, they employ local spaces attached either to subdomains, coarse cells, patches, or parameter regions. This suggests that “localized” denotes a family of approximation strategies rather than a single construction principle.

| Notion of locality | Typical construction | Representative papers |
|---|---|---|
| Spatial or patch-based | Local spaces on subdomains, overlapping patches, or coarse cells; global coupling by Galerkin, DG, or Schwarz ideas | [1710.02104], [1410.3253], [1401.7173] |
| Parameter-space or solution-space | Local spaces selected from nearby snapshots or clusters in parameter/snapshot space | [1204.3846], [2212.11934], [2510.10624] |
| Two-scale or optimization-integrated | Local reduced spaces embedded in multiscale or trust-region algorithms | [2211.09607], [2303.03074] |

A common misconception is that localization must mean domain decomposition in physical space. The parameter-local literature instead defines a local approximation space directly in parameter space, sometimes with a fixed online dimension \(N\), and chooses for each query parameter only the nearest snapshots according to an anisotropic metric [1204.3846]. By contrast, multiscale and domain-decomposition formulations localize in the computational domain and assemble a global approximation from locally supported contributions [1410.3253].

## 2. Spatial localization and multiscale constructions

In overlapping-domain formulations, the domain is covered by subdomains \(\{\omega_i\}_{i=1}^{N_D}\) with local spaces
\[
O_i := H^1_0(\omega_i),
\]
extended by zero outside \(\omega_i\). Starting from \(\widetilde V_0=\operatorname{span}\{0\}\), residual-based online enrichment computes the reduced solution \(\widetilde u_n\), forms the residual \(R_n(\cdot)=f(\cdot)-a(\widetilde u_n,\cdot)\), selects the patch with maximal local dual norm \(\|R_n\|_{O_i'}\), solves a local correction problem in that patch, and enriches the global reduced space by one localized function. A globally coupled variant instead solves on \(\widetilde V_n\oplus O_i\) for every patch and is optimal among all one-function enrichments drawn from local fine spaces [1710.02104].

In localized multiscale settings, locality is built into the trial space itself. The reduced basis localized orthogonal decomposition introduces a quasi-interpolation \(I_H\), the fine-scale kernel
\[
W_h:=\{v_h\in V_h\mid I_H(v_h)=0\},
\]
and local correctors \(Q_{h,k}^K(\Phi_z;\mu)\) on coarse patches \(U_k(K)\). The localized multiscale space is
\[
V_{H,k}^{\rm ms}(\mu)=\mathrm{span}\{\Phi_z+Q_{h,k}(\Phi_z;\mu)\mid z\in\mathcal N_H\},
\]
with \(Q_{h,k}(\Phi_z;\mu)=\sum_{K\subset\omega_z}Q_{h,k}^K(\Phi_z;\mu)\). This yields basis functions with local support and sparse global systems while retaining \(H^1\)-accuracy of order \(H\) when \(k\ge C_g|\log H|\) [1410.3253].

The localized reduced basis multiscale method (LRBMS) adopts a coarse partition \(\mathcal T_H\), local reduced spaces \(V_{\mathrm{red}}^T\subset V_h^{k,T}\) on each coarse element \(T\), and a global broken reduced space
\[
V_{H,\mathrm{red}}=\bigoplus_{T\in\mathcal T_H} V_{\mathrm{red}}^T,
\]
coupled by a generalized SWIPDG formulation. This construction combines numerical multiscale discretization and reduced basis compression, and it is explicitly designed to reduce the computational complexity with respect to both the multiscale parameter \(\varepsilon\) and the online parameter \(\mu\) [1401.7173].

## 3. Parameter-space localization and moving geometries

Parameter-local reduced basis methods define, for each query \(\mu\), a local parameter ball
\[
B_\mu=\{\tilde\mu\in\mathcal P\mid d(\mu,\tilde\mu)\le r(\mu)\},
\]
choose exactly \(N\) snapshots in \(B_\mu\), and form a local space \(\mathbb W_\mu\). The distance \(d\) is not fixed a priori: it is built from an empirically approximated Hessian \(H(\mu)\), an anisotropic metric tensor \(M(\mu)=V|\Lambda|V^T\), and a trapezoidal approximation of the geodesic distance. The method gives explicit control over online cost because the online reduced dimension is the user-chosen \(N\), independent of the total number \(K\) of stored snapshots [1204.3846].

For parameterized geometries on unfitted meshes, the main difficulty is that the active set of degrees of freedom changes with \(\mu\). One approach embeds every domain \(\Omega(\mu)\) in a parameter-independent background domain \(\Omega_0\), defines a background space \(V_{h,0}\), and extends each snapshot by zero on inactive basis functions, producing vectors \(\widehat{\mathbf u_h}(\mu)\in\mathbb R^{\mathcal N_{h,0}}\). Locality is then introduced by \(k\)-means clustering in parameter space,
\[
\mathcal P=\bigcup_{k=1}^{N_c}\mathcal P_k,
\]
followed by local POD on each cluster and nearest-centroid selection online. In this setting, localization is explicitly parameter-based rather than spatial [2212.11934].

For crack problems with moving discontinuities, the snapshots are first mapped to a common reference configuration,
\[
\tilde{\mathbf u}(\mu)=\mathbf u_h(\mu)\circ \boldsymbol\tau(\mu),
\]
then partitioned by fuzzy c-means into clusters \(\mathbb S_k\), reduced by local POD into bases \(\mathbb V_k\), and selected online by a classifier \(\hat{\mathcal C}:\mathcal P\to\{1,\dots,N_c\}\). This construction is tailored to Heaviside jumps and crack-tip singularities represented by XIGA enrichments, and it is explicitly presented as localization in solution or parameter space rather than domain decomposition [2510.10624].

## 4. Hyper-reduction, operator compression, and non-intrusive variants

Localized reduced spaces alone do not remove the online cost if the operators remain nonaffine. For unfitted geometries, one route is matrix and vector DEIM on the extended background mesh:
\[
\widehat{\mathbf A}(\mu)\approx \sum_{q=1}^{Q_a}\theta_q^a(\mu)\widehat{\mathbf A_q},\qquad
\widehat{\mathbf f}(\mu)\approx \sum_{q=1}^{Q_f}\theta_q^f(\mu)\widehat{\mathbf f_q},
\]
with POD on operator snapshots, greedy selection of “magic points,” and coefficient reconstruction from sampled entries. To avoid intrusive online operator assembly, the DEIM coefficients are interpolated by radial basis functions, and localization is applied cluster-wise to both DEIM and reduced bases [2212.11934].

A unified unfitted framework on parameterized domains combines a deformation-based map from a reference configuration to \(\Omega(\mu)\), aggregated unfitted finite element spaces, classical and tensor-train reduced bases, and localized dictionaries of reduced subspaces and hyper-reduction approximations obtained via matrix DEIM. In that formulation, the deformation map provides a common reference configuration, while localization produces separate dictionaries for reduced subspaces and for MDEIM approximations [2508.15320].

Localized operator approximation also appears in nonlinear finite volume models. For a spatially resolved three-dimensional lithium-ion battery model, a localized reduced basis scheme is built from non-conforming local approximation spaces and a localized empirical operator interpolation for the model’s nonlinearities, indicating the feasibility of localization beyond linear finite element settings [1606.05070].

A different path is non-intrusive reduction. In the crack-modeling framework, the local reduced coordinates
\[
\mathcal Z_k:\mathcal P^k\to\mathbb R^{N_k},\qquad \mu\mapsto \mathbb V_k^T\tilde{\mathbf u}(\mu),
\]
are approximated by cluster-wise artificial neural networks \(\hat{\mathcal Z}_k\), so that online evaluation requires classification, network inference, reduced reconstruction, and inverse mapping, without intrusive reduced operator assembly [2510.10624]. This suggests that localized reduced basis methods can be paired either with intrusive hyper-reduction or with data-driven coefficient surrogates.

## 5. Error control, online enrichment, and stability

One of the strongest theoretical results in the area is exponential convergence of online enrichment on overlapping domains. For the residual-based algorithm,
\[
\|u-\widetilde u_{n+1}\|_a \le c\,\|u-\widetilde u_n\|_a,\qquad c<1,
\]
and hence
\[
\|\widetilde u_n-u\|_a \le c^n\|u\|_a.
\]
The globally coupled enrichment strategy is moreover optimal among all single-function enrichments drawn from local fine spaces [1710.02104].

For LRBMS, error control is based on conservative flux reconstruction rather than the classical global residual estimator. The resulting localized estimator decomposes into a nonconformity term, a residual term, and a diffusive-flux term,
\[
\eta_{\mathrm{nc}}^T,\qquad \eta_{\mathrm{r}}^T,\qquad \eta_{\mathrm{df}}^T,
\]
all computable locally on each coarse element. This estimator is rigorous with respect to the weak solution and drives adaptive online enrichment through a localized Solve–Estimate–Mark–Refine loop [1501.05202].

For parabolic problems, elliptic reconstruction yields offline/online decomposable a posteriori bounds on both model reduction and full approximation error. The localized RB method is thereby extended to scalar parabolic evolution equations with true error certification, including discretization and model reduction contributions in one bound [1606.09216].

Localized reduced basis methods also have a close relationship with domain decomposition. The Localized Reduced Basis Additive Schwarz method interprets local residual-corrector problems as subdomain solves in additive Schwarz, and the resulting scheme can be viewed as a locally adaptive multi-preconditioning scheme for CG [2103.10884].

Stability, however, is not automatic. In ArbiLoMod for time-harmonic Maxwell’s equations in \(H(\mathrm{curl})\), the reduced inf-sup constant can drop under localized Galerkin projection, producing error spikes and motivating future work on localized test spaces for non-coercive problems [1606.06206]. By contrast, for unfitted Stokes and Navier–Stokes problems on deformed domains, supremizer enrichment is adapted to the reference configuration, and the paper states that the supremizer operator can be defined on the reference configuration without loss of stability [2508.15320].

## 6. Applications, reported performance, and current directions

The application range in the cited literature includes stationary heat conduction with high contrast and channels, parametric multiscale diffusion, Poisson and linear elasticity on trimmed or embedded geometries, time-harmonic Maxwell equations, pore-scale lithium-ion battery models, crack mechanics with XIGA, and PDE-constrained optimization with localized trust-region surrogates [1606.05070] [2303.03074] [2212.11934].

| Setting | Reported outcome | Paper |
|---|---|---|
| 2D Poisson on trimmed domains | Global \(N=182\), \(Q_a=349\), \(Q_f=124\); local max \(N=35\), max \(Q_a\approx 66\), max \(Q_f\approx 23\); global online \(\sim 879\) ms vs local \(\sim 85\)–95 ms | [2212.11934] |
| Quarter-ring linear elasticity | Local RB speedup \(\sim 130\times\) vs FOM; global \(\sim 16.5\times\) | [2212.11934] |
| Crack ROMs on splines/XIGA | 2D edge crack: \(0.03\) s vs FOM \(10\) s, speedup \(\approx 333\times\); 3D edge crack: \(0.0423\) s vs FOM \(35\) s, speedup \(\approx 827\times\) | [2510.10624] |
| Large-scale PDE-constrained optimization | Relaxed TR-LRBM used \(2\) DG-MsFEM evaluations, \(506\) LRBM evaluations, and \(294\) local patch evaluations, versus \(259\) DG-MsFEM evaluations for BFGS with DG-MsFEM | [2303.03074] |

The broader literature positions localized reduced basis methods at the intersection of reduced basis methods, domain decomposition, multiscale discretization, and hyper-reduction. Representative neighboring frameworks named explicitly in the data include LRBMS, GMsFEM, CEM-GMsFEM, ArbiLoMod, Localized DEIM, hp-RB methods, multi-domain empirical interpolation, and the Localized Orthogonal Decomposition [1410.3253] [1710.02104] [1606.06206].

Several open issues recur. One is certification: some geometry-localized and non-intrusive approaches report numerical error and singular value decay but do not implement explicit certified a posteriori estimators [2212.11934] [2510.10624]. Another is stability for non-coercive or saddle-point systems, where localized test-space design and supremizer strategies remain essential [1606.06206] [2508.15320]. A further direction is optimization and multiscale control: adaptive localized reduced basis methods have been embedded in trust-region algorithms, and the thesis literature presents TSRBLOD and TR-LRB methods that avoid FEM evaluations of the involved systems during the optimization loop [2211.09607]. This suggests that the field has evolved from localized approximation spaces toward localized, certified, and optimization-aware reduced models.

Source: https://www.emergentmind.com/topics/localized-reduced-basis-methods