---
title: Shape-Adaptive Reduced Basis Methods
url: https://www.emergentmind.com/topics/shape-adaptive-reduced-basis-method-rbm
type: topic
---

# Shape-Adaptive Reduced Basis Methods

Searching arXiv for recent and foundational papers on shape-adaptive reduced basis methods and related geometric RBM formulations.
{"query":"shape-adaptive reduced basis method geometry parameterized boundary integral reduced basis arXiv", "max_results": 10}
{"query":"2403.10933 Reduced Basis Method for the Elastic Scattering by Multiple Shape-Parametric Open Arcs in Two Dimensions", "max_results": 5}
Searching by exact arXiv id and title to verify the primary source and gather related work.
{"query":"id:2403.10933", "max_results": 5}
Shape-adaptive reduced basis methods (RBMs) are model-order reduction procedures for parametrized PDE and boundary-integral problems in which geometric variation is part of the parameter space and the reduced approximation is constructed to remain effective under shape change. In the literature represented by Pinto–Henríquez’s treatment of elastic scattering by multiple open arcs, by reduced-basis optimization of plasmonic nanoparticle arrays, by rigorous source–mask optimization in 3D computational lithography, and by divergence-conforming reduced models for geometric Navier–Stokes flow, shape adaptivity is realized through parameter-to-geometry maps, reference-domain pullbacks, Piola transforms, or basis functions that depend explicitly on the current shape parameters [2403.10933], [2508.04198], [1011.2902], [1807.11866].

## 1. Geometric parametrization and the meaning of shape adaptivity

A central ingredient is an explicit parametrization of geometry. For multiple elastic cracks in two dimensions, one assumes \(M\) disjoint open arcs \(\Gamma_1,\dots,\Gamma_M\), with each arc \(\Gamma_j\) given as the image of the reference interval \(\hat S=[-1,1]\) under a smooth parameter-to-geometry map
\[
X_j(\cdot;\mu_j):\hat S\to\mathbb R^2,\qquad s\mapsto X_j(s;\mu_j),
\]
where \(\mu_j\in\mathbb R^{P_j}\), and the full parameter is \(\mu=(\mu_1,\dots,\mu_M)\in\mathcal P\subset\mathbb R^P\), \(P=\sum_j P_j\). A common choice expands the shape around a nominal straight reference curve:
\[
X_j(s;\mu_j)=\xi_j(s)+\sum_{p=1}^{P_j}\mu_{j,p}\phi_{j,p}(s),\qquad s\in[-1,1].
\]
The geometry-dependent normal \(n_j(s;\mu_j)\) and metric factor \(|\partial_s X_j(s;\mu_j)|\) then enter the operators directly [2403.10933].

For plasmonic nanoparticle arrays, the configuration of \(M\) analytic nanoparticles is parameterized by
\[
w=(w_1,\dots,w_M)\in W\subset\mathbb R^{MP},
\]
where each \(w_m\) encodes the shape, for instance ellipse semi-axes, orientation, and center. In the 3D source–mask optimization setting, the parameter vector is split into geometry and illumination:
\[
\mu=(\mu_{\mathrm{geom}},\mu_{\mathrm{src}})\in D=D_{\mathrm{geom}}\times D_{\mathrm{src}}\subset\mathbb R^7,
\]
with \(\mu_{\mathrm{geom}}=(w,h,d)\) for mask geometry and \(\mu_{\mathrm{src}}=(r_1,r_2,b_1,b_2)\) for the quadrupole source. In the incompressible-flow formulation, geometric dependence is expressed by a smooth invertible map \(\chi_\mu:\Omega(\mu)\to\hat\Omega\), or equivalently \(F_\mu=\chi_\mu^{-1}:\hat\Omega\to\Omega(\mu)\), between a reference domain and the physical domain [2508.04198], [1011.2902], [1807.11866].

Within these formulations, “shape-adaptive” does not denote a single algorithmic template. In some cases the reduced basis itself depends explicitly on the current geometry; in others the basis is fixed on a reference domain and transported to the physical domain by a geometry-dependent mapping. This suggests that shape adaptivity is best understood as a structural property of the reduced model rather than a specific discretization choice.

## 2. High-fidelity formulations for parametrized geometry

The full-order problems differ by application class, but each exposes geometry in a way that enables reduced modeling.

For elastic scattering by multiple open arcs, the exterior Lamé boundary-value problem is transformed by single-layer and double-layer potentials into a coupled system of boundary integral equations on \(\Gamma=\bigcup_j \Gamma_j\):
\[
\begin{pmatrix}
\tfrac12 I + K(\mu) & -V(\mu)\\[6pt]
W(\mu) & \tfrac12 I - K'(\mu)
\end{pmatrix}
\begin{pmatrix}
\phi\\
\psi
\end{pmatrix}
=
\begin{pmatrix}
f_D(\mu)\\
f_N(\mu)
\end{pmatrix}.
\]
For \(x\) on \(\Gamma_i(\mu)\) and \(y\) on \(\Gamma_j(\mu)\),
\[
[V(\mu)\sigma]_i(s)
=\sum_{j=1}^M\int_{-1}^1
G\bigl(X_i(s;\mu_i),X_j(t;\mu_j)\bigr)\sigma_j(t)|\partial_tX_j(t;\mu_j)|\,dt,
\]
\[
[K(\mu)\sigma]_i(s)
=\sum_{j=1}^M\int_{-1}^1
\partial_{n_y}G\bigl(X_i(s;\mu_i),X_j(t;\mu_j)\bigr)\sigma_j(t)|\partial_tX_j(t;\mu_j)|\,dt,
\]
with analogous expressions for the hypersingular operator \(W(\mu)\) and the adjoint \(K'(\mu)\). Here \(G\) is the elastic fundamental solution [2403.10933].

For plasmonic nanoparticles under the TM approximation, the forward problem is written as a parameterized layer-potential system for densities \(\phi,\varphi\in L^2(\partial D)\):
\[
\begin{pmatrix}
\mathcal S_w^{k_c} & -\,\mathcal S_w^{k_m}\\[1ex]
\dfrac1{\varepsilon_c}\bigl(-\tfrac12I+(\mathcal K_w^{k_c})^*\bigr) &
-\dfrac1{\varepsilon_m}\bigl(\tfrac12I+(\mathcal K_w^{k_m})^*\bigr)
\end{pmatrix}
\begin{pmatrix}
\phi\\
\varphi
\end{pmatrix}
=
\begin{pmatrix}
u^i(x(t;w))\\
\tfrac1{\varepsilon_m}\partial_\nu u^i(x(t;w))
\end{pmatrix},
\]
with parameterized single-layer and Neumann–Poincaré operators \(\mathcal S_w^k\) and \((\mathcal K_w^k)^*\) [2508.04198].

In 3D computational lithography, the truth approximation is a finite-element discretization of time-harmonic Maxwell equations in a unit cell \(\Omega(\mu)\). Variationally, one seeks \(E(\cdot;\mu)\in V(\mu)\) such that
\[
a_\mu(E,v)=f_\mu(v)\qquad \text{for all }v\in V(\mu),
\]
where
\[
a_\mu(E,v)=\int_{\Omega(\mu)}\bigl[\mu_r^{-1}\,\mathrm{curl}\,E\cdot \mathrm{curl}\,v^*
-k_0^2\,\varepsilon_r\,E\cdot v^*\bigr]\,dx.
\]
Discretization with edge elements yields the algebraic system
\[
A(\mu)\,u(\mu)=b(\mu),\qquad u(\mu)\in\mathbb C^N.
\]
In divergence-conforming Navier–Stokes reduction, the high-fidelity problem is first posed on the reference domain and discretized so that the velocity snapshots are pointwise divergence-free; the geometric dependence then enters the pulled-back forms and the Piola transport [1011.2902], [1807.11866].

A recurring theme is that the geometric parameters appear inside kernels, coefficients, domain maps, or test spaces. That feature is precisely what necessitates shape-aware reduction.

## 3. Construction of reduced spaces

The reduced spaces in shape-adaptive RBM are not uniform across applications.

In Pinto–Henríquez’s multiple-arc scattering framework, the offline stage begins by ignoring inter-arc coupling and solving \(M\) independent single-arc boundary-integral equations over a training set \(\{\mu^1,\dots,\mu^{N_{\rm train}}\}\subset\mathcal P\). On arc \(j\), the high-fidelity density is
\[
u_j(\mu^i)=[\phi_j(\cdot;\mu^i),\psi_j(\cdot;\mu^i)].
\]
Collecting all such solutions yields the snapshot set \(S=\{u(\mu^i)\}_{i=1}^{N_{\rm train}}\subset\mathcal X\). Proper Orthogonal Decomposition (POD) is then applied: with correlation matrix \(C_{kl}=\langle u(\mu^k),u(\mu^l)\rangle_{\mathcal X}\) and dominant eigenpairs \((\lambda_n,v_n)\), the POD basis functions are
\[
\phi_n=\frac1{\sqrt{\lambda_n}}\sum_{i=1}^{N_{\rm train}}(v_n)_i\,u(\mu^i),\qquad n=1,\dots,N.
\]
These basis functions are orthonormal in \(\mathcal X\) and capture the largest energy of the snapshot set [2403.10933].

In source–mask optimization, the reduced space is built by a Greedy algorithm over a dense training set \(\Xi_{\rm train}\subset D\). One initializes with a truth solution \(u(\mu_1)\), repeatedly computes reduced approximations \(u_n(\mu)\) and estimators \(\Delta_n(\mu)\), selects
\[
\mu_{n+1}=\arg\max_{\mu\in\Xi_{\rm train}}\Delta_n(\mu),
\]
and enlarges the basis by the newly computed truth solution. The resulting space is
\[
U_N=\mathrm{span}\{u(\mu_1),\dots,u(\mu_N)\}\subset\mathbb C^N.
\]
This is described as a self-adaptive RB construction for variable geometrical parameters [1011.2902].

The plasmonic absorber formulation develops a more explicit form of shape adaptation. Spectral theory for the Laplace NP-operator on each analytic particle \(D_m=D_m(w_m)\) provides an orthogonal eigenbasis
\[
(\mathcal K_{w_m})^*[\psi_{m,i}^{\pm}]=\pm\alpha_{m,i}\psi_{m,i}^{\pm},\qquad
\alpha_{m,i}=\tfrac12 e^{-2i\rho_m},
\]
where \(\rho_m=\cosh^{-1}(a_m/c_m)\) for an ellipse of semi-axes \((a_m,b_m)\) and focal distance \(c_m\). Retaining the first \(N\) eigenmodes \(\{\psi_{m,1},\dots,\psi_{m,N}\}\), one approximates
\[
\phi_m(t)\approx \sum_{i=1}^N c_{m,i}^{\phi}\psi_{m,i}(t;w_m),\qquad
\varphi_m(t)\approx \sum_{i=1}^N c_{m,i}^{\varphi}\psi_{m,i}(t;w_m).
\]
Since the eigenfunctions depend explicitly on the shape parameters \(w_m\), the resulting basis “adapts” as the optimizer moves in \(w\) [2508.04198].

For steady Navier–Stokes flow, shape adaptivity is achieved differently. The reduced basis is formed from divergence-conforming high-fidelity velocity snapshots \(\{\hat{\boldsymbol u}_0^h(\mu_i)\}_{i=1}^N\) on the reference domain, followed by a POD in the \(H^1\)-seminorm with covariance matrix
\[
C^v_{ij}=\bigl(\hat{\boldsymbol u}_0^h(\mu_i),\hat{\boldsymbol u}_0^h(\mu_j)\bigr)_{H^1(\hat\Omega)}.
\]
The reduced basis functions are
\[
\hat{\boldsymbol u}^r_j
=\frac{1}{\sqrt{\lambda_j}}\sum_{i=1}^N v_i^j\,\hat{\boldsymbol u}_0^h(\mu_i),
\qquad j=1,\dots,M,
\]
and are transported to the physical domain by the Piola transform, which preserves divergence-freeness for every \(\mu\) [1807.11866].

## 4. Offline–online decomposition and affine compression

The practical value of RBM depends on separating expensive parameter-independent work from inexpensive parameter-dependent work.

For multiple open arcs, the principal difficulty is nonaffinity: the operators \(V(\mu)\) and \(K(\mu)\) depend on \(\mu\) through the geometry inside the kernel. A modified Empirical Interpolation Method (EIM) is applied to each kernel
\[
\kappa_{ij}(s,t;\mu)=G\bigl(X_i(s;\mu_i),X_j(t;\mu_j)\bigr)
\]
to obtain
\[
\kappa_{ij}(s,t;\mu)\approx \sum_{q=1}^Q \theta_{ij,q}(\mu)\,\kappa_{ij,q}(s,t).
\]
Hence
\[
V(\mu)\approx \sum_{q=1}^Q \Theta_q^V(\mu)\,V_q,\qquad
K(\mu)\approx \sum_{q=1}^Q \Theta_q^K(\mu)\,K_q,
\]
where \(\mu^0\) is a nominal parameter used to fix the metric in the basis operators. With \(\Phi\in\mathcal X\times N\) denoting the matrix of POD basis vectors, the reduced system is
\[
A^N(\mu)x^N(\mu)=f^N(\mu),\qquad
A^N(\mu)=\sum_{q=1}^Q \Theta_q(\mu)\,[\Phi^T A_q\Phi],\qquad
f^N(\mu)=\Phi^T f(\mu).
\]
The offline complexity for assembling and compressing the \(Q\) basis operators is \(O(Q\cdot M\cdot N_h^2)\), while the POD costs \(O(N_{\rm train}^2\cdot M\cdot N_h)\). Online, evaluating the coefficients costs \(O(Q\cdot P)\), assembling \(A^N(\mu)\) costs \(O(Q\cdot N^2)\), and solving the reduced system by a direct method costs \(O(N^3)\), giving an online cost of order \(O(Q\cdot N^2+N^3+Q\cdot P)\) [2403.10933].

In the plasmonic setting, the offline–online split is organized around singularity extraction. The singular part \(\mathbf M_1(w)\) is wavelength-independent and assembled offline from NP-eigenfunction captures, while the smooth remainder \(\mathbf M_2(w,\lambda)\) is assembled online by trapezoidal rule for each \(\lambda\). The forward reduced system reads
\[
\mathbf M(w,\lambda)\mathbf c=\mathbf f(w,\lambda),\qquad
\mathbf M(w,\lambda)=\mathbf M_1(w)+\mathbf M_2(w,\lambda),
\]
and the adjoint problem has the analogous decomposition
\[
\mathbf T(w,\lambda)\mathbf d=\mathbf g(w,\lambda),\qquad
\mathbf T(w,\lambda)=\mathbf T_1(w)+\mathbf T_2(w,\lambda).
\]
By expanding in NP-eigenfunctions, the \(O(\log|t-s|)\) singularity is isolated in closed form [2508.04198].

In the Maxwell and Navier–Stokes formulations, the standard RBM affine decomposition is used. For lithography,
\[
A(\mu)=\sum_{q=1}^{Q_a}\theta_a^q(\mu)A_q,\qquad
b(\mu)=\sum_{q=1}^{Q_b}\theta_b^q(\mu)b_q,
\]
so that the reduced matrices \(A_N^q=U_N^H A_q U_N\) and vectors \(b_N^q=U_N^H b_q\) are precomputed offline, and the online assembly of
\[
A_N(\mu)=\sum_{q=1}^{Q_a}\theta_a^q(\mu)A_N^q,\qquad
b_N(\mu)=\sum_{q=1}^{Q_b}\theta_b^q(\mu)b_N^q
\]
costs \(O(N^2Q_a+NQ_b)\), followed by an \(O(N^3)\) solve. For geometric incompressible flow, the pulled-back forms admit affine expansions such as
\[
a(u,w;\mu)=\sum_{i=1}^{N_a}\theta_i^a(\mu)\,\hat a_i(u,w),\qquad
c(u,v,w;\mu)=\sum_{k=1}^{N_c}\theta_k^c(\mu)\,\hat c_k(u,v,w),
\]
and the online cost is \(\mathcal O(M^2N_a+M^3N_c)\), independent of the high-fidelity dimension [1011.2902], [1807.11866].

## 5. Certification, stability, and methodological distinctions

A defining property of RBM in these works is the presence of computable error surrogates or preserved structural constraints.

For elastic scattering by shape-parametric open arcs, the reduced residual is
\[
r_N(\mu)=f(\mu)-A(\mu)u^N(\mu),\qquad u^N(\mu)=\Phi x^N(\mu),
\]
and the certified error bound is
\[
\|u(\mu)-u^N(\mu)\|_{\mathcal X}
\le
\frac{\|r_N(\mu)\|_{\mathcal X'}}{\beta_{\rm LB}(\mu)}
=:\Delta_N(\mu),
\]
where \(\beta_{\rm LB}(\mu)\) is a lower bound for the inf–sup constant of the full operator and can be certified offline by a Successive Constraint Method. In the 3D Maxwell setting, the residual is
\[
r_N(v;\mu)=f_\mu(v)-a_\mu(u_N(\mu),v)\qquad \text{for all } v\in V,
\]
and if \(a_\mu(v,v)\ge \alpha(\mu)\|v\|_V^2\), then
\[
\|u(\mu)-u_N(\mu)\|_V
\le
\Delta_N(\mu)
:=
\frac{\|r_N(\cdot;\mu)\|_{V'}}{\alpha(\mu)}.
\]
In practice one uses a cheaply computable lower bound \(\alpha_{\rm LB}(\mu)\), and the estimator drives the Greedy basis generation [2403.10933], [1011.2902].

The plasmonic absorber paper defines residual-based estimators
\[
r_F(\mathbf c)=\|\mathbf M\mathbf c-\mathbf f\|_2,\qquad
r_A(\mathbf d)=\|\mathbf T\mathbf d-\mathbf g\|_2,
\]
and states that standard RBM bounds of the form
\[
\|u-u_N\|\le \frac{r_F}{\sigma_{\min}(\mathbf M_1)}
\]
provide certification. In the Navier–Stokes setting, the main structural issue is not an a posteriori bound but stability under geometric variation: the Piola transform guarantees that \(\nabla_x\cdot \Piola(\mu)\hat{\boldsymbol u}_j^r=0\) for all \(\mu\), and pressure is recovered a posteriori at negligible extra cost [2508.04198], [1807.11866].

Several common misconceptions are clarified by these constructions. Shape adaptivity does not require repeated remeshing of a full-order model: in the nanoparticle optimization framework, no remeshing or velocity-field perturbation is required because shape sensitivities enter only through derivatives of the boundary operators. Nor does geometry parametrization necessarily destroy incompressibility: in the divergence-conforming approach, exact solenoidality is retained under geometric transformations by the Piola transform. Similarly, in the multi-arc scattering problem, the offline stage need not resolve the full coupled configuration; the basis is built from decoupled single-arc problems, while the inter-arc interaction is represented later through affine compression [2508.04198], [1807.11866], [2403.10933].

## 6. Applications and reported numerical behavior

The reported numerical results show that shape-adaptive RBM is used both for rapid many-query simulation and for optimization loops with geometric parameters.

For elastic scattering by multiple open arcs, representative tests with \(M=3\) arcs and \(P=10\) shape parameters per arc show that the maximum reduced-basis error
\[
\mathrm{err}(N)=\max_{\mu\in\mathcal P_{\rm test}}\|u(\mu)-u^N(\mu)\|_{L^2(\Gamma)}
\]
decays exponentially in \(N\), reaching \(10^{-6}\) accuracy by \(N\approx 30\). With \(N=30\), \(Q=50\), \(M=3\), and \(P=30\), the online solve takes \(O(10^{-2}\,\mathrm s)\), versus \(O(10\,\mathrm s)\) for a full BIE. As \(M\) increases, the RB approach scales nearly linearly in \(M\), while the full BIE scales like \(O((M\cdot N_h)^3)\); the bound \(\Delta_N(\mu)\) remains tight even as \(P\) increases [2403.10933].

For broadband absorber design with plasmonic nanoparticles, the online stage with \(M\sim 100\) and \(N\sim 10\) is reduced by two orders of magnitude relative to classical BEM. In single- and multiple-particle tests, the RBM achieves relative errors below \(10^{-6}\) in extinction cross-sections and adjoint fields, whereas a 200-point Nyström reference gives \(10^{-2}\)–\(10^{-1}\). Online solutions of both forward and adjoint problems are \(\sim 100\times\) faster than classical BEM at comparable accuracy, and the projected gradient-descent loop converges in \(\mathcal O(10^3)\) iterations for broadband targets over \([150,550]\) nm [2508.04198].

In 3D source–mask optimization, the truth problem has \(N\approx 340\,000\) degrees of freedom and one full solution takes approximately \(10\,000\) s. A reduced basis with \(N=40\) snapshots yields \(H(\mathrm{curl})\)-norm error \(10^{-6}\) and diffraction-order output error below \(10^{-4}\) at \(N=40\), with an online cost of approximately \(1\) s per solve, corresponding to a speed-up of about \(10^4\). An SMO run with approximately \(10^3\) forward solves is reduced from approximately \(10^7\) s to approximately \(10^3\) s, and the optimal parameters enlarge the process window by about \(30\%\) compared with mask-only optimization [1011.2902].

For steady flow around a NACA0015 airfoil, using \(M_v=M=30\) divergence-free velocity modes yields relative velocity error \(\|u_{\rm HF}-u_{\rm RB}\|_{H^1}/\|u_{\rm HF}\|_{H^1}\approx 10^{-4}\) and relative pressure error \(\|p_{\rm HF}-p_{\rm RB}\|_{L^2}/\|p_{\rm HF}\|_{L^2}\approx 10^{-3}\), uniformly over sampled angles of attack. The divergence-conforming Piola-based RB runs in roughly \(10^{-3}\) s per query, versus on the order of \(10^{-2}\) s for a Taylor–Hood supremizer-stabilized RB of comparable accuracy, and remains exactly divergence-free to machine precision for every \(\mu\) [1807.11866].

Taken together, these results indicate that shape-adaptive RBM is not confined to a single physics domain. It appears in boundary-integral scattering, nanoparticle inverse design, rigorous electromagnetic lithography, and incompressible flow, with the common objective of replacing expensive geometry-dependent solves by low-dimensional, rapidly evaluable, and in several cases certified surrogates.

Source: https://www.emergentmind.com/topics/shape-adaptive-reduced-basis-method-rbm