---
title: Semi-analytical Solver (SAS) Overview
url: https://www.emergentmind.com/topics/semi-analytical-solver-sas
type: topic
---

# Semi-analytical Solver (SAS) Overview

Searching arXiv for the cited SAS-related papers and closely related formulations.
Search query: arXiv 1412.0904 semi-analytic solutions power system differential algebraic equations; 1802.04134 differential transformation method power system; 1804.07330 time-power series SAS; 1810.11418 semi-analytical orbit propagation STELA; 2210.03298 natural gas network semi-analytical solution; 1409.1113 Grad-Shafranov semi-analytical solver; 2506.15185 semi-analytical forward solver elastic parameters; 2509.16766 Noh problem semi-analytic solutions.
A semi-analytical solver (SAS) is a solver architecture in which part of a governing model is handled in closed or semi-closed form while the remaining computation is carried out by localized numerical evaluation, reduced auxiliary solves, or short-window propagation. In the power-system formulation of SAS, the semi-analytic solution is a closed-form function of symbolic variables including time, the initial state and the parameters on system operating conditions, and hence able to directly give trajectories on system state variables that are accurate for at least a certain time window [1412.0904]. In orbital mechanics, semi-analytical methods remove fast oscillations analytically and integrate only the slow evolution of mean elements [1810.11418]. Other arXiv uses of the term include multivariate-polynomial cellwise solutions for natural-gas PDEs [2210.03298], eigenfunction-based semi-analytical equilibrium construction for the Grad–Shafranov equation [1409.1113], a direct method of lines that preserves analytical radial dependence in elasticity [2506.15185], and self-similar shock solutions closed by Rankine–Hugoniot algebraic systems for the Noh problem with a black-box equation of state [2509.16766].

## 1. Definitional scope and common structure

The phrase “semi-analytical solver” does not denote a single algorithm. In the cited literature it designates a class of methods that retain an analytic representation for a strategically chosen part of the problem while delegating the rest to recursion, averaging, collocation, linear algebra, or localized numerical integration. In power systems this analytic part is typically a short-time series in time, parameterized by initial conditions and operating-point parameters [1412.0904]. In orbit propagation it is the analytical removal and later recovery of short-periodic terms from osculating elements [1810.11418]. In elasticity it is the analytical solution along one coordinate after discretization in another [2506.15185].

Taken together, these formulations suggest a common structural pattern. First, the governing equations are rewritten so that one component of the dynamics becomes symbolically tractable: a decomposition series, a mean-element system, a polynomial ansatz, a separable eigenexpansion, or a self-similar reduction. Second, the semi-analytical representation is truncated or specialized to a finite basis. Third, the resulting object is evaluated or fitted over finite windows, cells, or regions. Fourth, auxiliary mechanisms control validity: divergence indicators in power-system transient simulation, resonance-retention criteria and perturbation truncation in orbit propagation, collocation and continuity conditions in gas networks, boundedness and eigenvalue selection in elasticity, or Newton solves for Rankine–Hugoniot closure in the Noh problem [1412.0904; 1810.11418; 2210.03298; 2506.15185; 2509.16766].

A plausible implication is that “semi-analytical” in current computational practice is best understood operationally rather than taxonomically: it identifies a solver in which analytic structure is preserved only to the extent that it yields a computational advantage without requiring a fully closed-form global solution.

## 2. Power-system SAS formulations

In transient stability simulation, the baseline model is a nonlinear differential-algebraic system
$$
\dot x(t)=f(x(t),y(t),p), \qquad 0=g(x(t),y(t),p),
$$
with $x(t)\in\mathbb{R}^m$ the dynamic states, $y(t)\in\mathbb{R}^k$ the algebraic variables, and $p$ the fixed parameters such as inertia $H$, damping $D$, machine EMFs, network admittances, and mechanical power inputs [1412.0904]. In practice one often eliminates $y$ via network equations to obtain coupled ODEs. For the classical 2nd-order model of a $K$-machine system,
$$
M\ddot\delta + D\dot\delta = P_m - P_e(\delta),
$$
or, in first-order form,
$$
\dot\delta=\omega, \qquad M\dot\omega + D\omega = P_m - P_e(\delta).
$$

The Adomian-decomposition formulation derives the SAS by expanding both the state trajectory and the nonlinear term in infinite series and then obtaining each term recursively in closed form [1412.0904]. For an ODE of the form $\ddot\delta+a\dot\delta=F(\delta)$, the method applies the Laplace transform, assumes
$$
\delta(t)=\sum_{n=0}^\infty \delta_n(t), \qquad F(\delta(t))=\sum_{n=0}^\infty A_n(\delta_0,\ldots,\delta_n),
$$
with Adomian polynomials
$$
A_n = \frac{1}{n!}\frac{d^n}{d\lambda^n}\left[F\!\left(\sum_{k=0}^\infty \delta_k\lambda^k\right)\right]_{\lambda=0},
$$
and derives recursive relations in the Laplace domain. Truncation after $N$ terms yields
$$
\delta_{\mathrm{SAS}}(t)=\sum_{n=0}^{N-1}\delta_n(t),
$$
which is accurate over a finite time window $T$ that depends on $N$ and on system parameters [1412.0904]. A modified ADM may fold linear damping terms and certain constant parts into the zeroth term, so that each $\delta_n(t)$ becomes a polynomial in $t$ plus possibly sinusoids.

The differential transformation method (DTM) produces the same general short-time objective through a different coefficient calculus [1802.04134]. For an analytic function $x(t)$, the transform and inverse transform are
$$
X(k)=\frac{1}{k!}\left.\frac{d^k x(t)}{dt^k}\right|_{t=0}, \qquad
x(t)=\sum_{k=0}^{\infty}X(k)t^k.
$$
Applied to an index-0 DAE,
$$
\dot\varphi(t)=f(\varphi(t),\psi(t)), \qquad \psi(t)=g(\varphi(t)),
$$
the DTM yields recursions
$$
\Phi^{(k)} = \frac{1}{k}F\!\bigl(\Phi^{(0)},\dots,\Phi^{(k-1)},\Psi^{(0)},\dots,\Psi^{(k-1)}\bigr), \qquad
\Psi^{(k)} = G\!\bigl(\Phi^{(0)},\dots,\Phi^{(k)}\bigr),
$$
followed by the truncated expansions
$$
\varphi(t)=\sum_{k=0}^K \Phi^{(k)} t^k, \qquad \psi(t)=\sum_{k=0}^K \Psi^{(k)} t^k.
$$
The distinguishing claim is computational: the recursive formula avoids generating and storing its complete symbolic expression, which makes both derivation and evaluation of the SAS more efficient especially for multi-machine power systems [1802.04134].

A third power-system variant writes the SAS directly as an $N$th-order power series about $t_0$,
$$
x_{\mathrm{sas}}(t)=\sum_{k=0}^N a_k (t-t_0)^k,
$$
with coefficients obtained by matching the Taylor expansion of $f(x(t))$ and the recursion
$$
a_{k+1}=\frac{b_k}{k+1}, \qquad b_k=b_k(a_0,a_1,\dots,a_k)
$$
[1804.07330]. This framework extends to general DAEs through the dynamic-bus method, in which every PV or load bus is represented as a dynamic bus with voltage series in time, component SASs are evaluated in parallel, the network algebraic equations are solved once at the end of the window, and the voltage-series coefficients are re-fitted over the interval [1804.07330].

## 3. Windowing, offline–online execution, and error control

The power-system literature gives the SAS a distinct offline–online decomposition. In the two-stage SAS-based approach, the offline stage selects symbolic variables such as time, initial states, parameters, and optionally selected admittance matrix entries to cover multiple contingency topologies; applies Multi-stage Modified ADM to derive symbolically $\delta_n(t)$ up to order $N$ for each generator; stores the closed-form expressions; and estimates the maximum length of window of accuracy $R_A$ [1412.0904]. The online stage carries out numerical integration only through the fault-on period to obtain the post-fault initial state, evaluates the SAS for all generators over the first interval, extracts the terminal state, and repeats sequentially until the horizon is covered.

Local validity is central. Because a truncated SAS departs from the true solution beyond some $R_A$, the power-system ADM paper defines a divergence indicator
$$
I_{\mathrm{LOA}}(t)=\frac{d}{dt}[A_{N+1}(t)],
$$
and chooses a threshold $I_{\mathrm{LOA},\max}$ so that $I_{\mathrm{LOA}}(R_A)=I_{\mathrm{LOA},\max}$ [1412.0904]. In a multi-machine case one computes $R_{A,i}$ per generator and then picks $T\le \min_i R_{A,i}$. The time-power-series formulation uses an error-rate upper bound
$$
r(h)=\sup_{t\in[t_0,t_0+h]}\left\|f(x_{\mathrm{sas}}(t))-\dot x_{\mathrm{sas}}(t)\right\|,
$$
with cumulative error bounded by $r(h)h$, and selects the largest $h$ satisfying $r(h)h\le \varepsilon$ [1804.07330]. These are two different mechanisms for the same operational requirement: finite-window reliability.

The DTM formulation makes the offline–online split even more explicit. Offline, one derives once the recursive rules symbolically in terms of the initial state and parameters. At runtime, whenever a new contingency or operating point changes the initial state or parameters, the same recursion rapidly computes the numeric coefficients for the current window [1802.04134]. The paper further studies an optimal order $K$ by balancing increased window length $t_w(K)$ against increased per-window evaluation time, with total time
$$
T_{\mathrm{total}}(K)=\frac{T_{\mathrm{sim}}}{t_w(K)}\,t_{\mathrm{eval}}(K).
$$
The reported example on the IEEE 39-bus system states that, for tolerance $10^{-5}$ pu, $K=12$ and $t_w=0.20$ s were optimal [1802.04134].

A recurrent misconception is to equate SAS with a single long-horizon closed-form trajectory. The cited formulations do not support that interpretation. In these power-system papers, the semi-analytic object is explicitly local, and long simulations are obtained by repeated reinitialization over short windows [1412.0904; 1802.04134; 1804.07330].

## 4. SAS in orbit propagation, gas networks, plasma equilibrium, elasticity, and shock verification

In Earth-orbit propagation, semi-analytical methods seek to retain only the slow evolution of orbital elements by removing fast orbital-periodic oscillations analytically [1810.11418]. The state is written in a convenient element set, for example the non-singular equinoctial elements
$$
E=\{a,h,k,P,Q,\lambda\},
$$
with
$$
h=e\sin(\omega+\Omega),\quad k=e\cos(\omega+\Omega),\quad
P=\sin(i/2)\cos\Omega,\quad Q=\sin(i/2)\sin\Omega,\quad
\lambda=M+\omega+\Omega.
$$
Osculating equations are written in Gauss form, perturbations are expanded in Fourier and Legendre series, short-periodic corrections are separated from mean elements,
$$
\hat E_i = E_i + \sum_{j=1}^\infty \epsilon^j \eta_{i,j}(E,t),
$$
the mean-element rates are obtained by averaging over the fast angle or angles, and the averaged equations are then integrated with a numerical solver [1810.11418]. The computational benefit is that, by eliminating the $\sim 1$-revolution oscillations analytically, the numerical solver may take much larger steps, on the order of days.

In natural-gas transients, SAS is a cellwise polynomial approximation for the isothermal pipeline PDE system
$$
\frac{\partial p(x,t)}{\partial t}+\frac{v^2}{S}\frac{\partial q(x,t)}{\partial x}=0,\qquad
\frac{\partial q(x,t)}{\partial t}+S\frac{\partial p(x,t)}{\partial x}
+\frac{\lambda v^2}{2dS}\frac{q(x,t)|q(x,t)|}{p(x,t)}=0
$$
[2210.03298]. The $(x,t)$ domain is partitioned into uniform cells of size $\Delta L\times \Delta T$, coordinates are normalized to $[0,1]\times[0,1]$, and in each cell the SAS ansatz is
$$
p^{(e,I)}(x,t,s)=\sum_{r=0}^R \sum_{m=0}^M \sum_{n=0}^m p^{(e,I)}_{n,m-n,r} x^n t^{m-n} s^r,
$$
$$
q^{(e,I)}(x,t,s)=\sum_{r=0}^R \sum_{m=0}^M \sum_{n=0}^m q^{(e,I)}_{n,m-n,r} x^n t^{m-n} s^r,
$$
where the coefficients are determined by PDE residual equations, initial-boundary projection at collocation points, and seamless coupling between adjacent cells. Each time slice reduces to a square linear system $Ac=b$, and a reduced scheme, SAS-2, replaces the friction term by a cell-wise linear model,
$$
s\,C_3\frac{q^2}{p}\approx s\,C_3 C_4 q,
$$
to simplify assembly and reduce run time with minor loss of accuracy [2210.03298].

For the Grad–Shafranov equation, the solver becomes semi-analytical after specific constitutive assumptions render the PDE linear [1409.1113]. With pressure proportional to the poloidal flux and $F^2(\psi)$ quadratic in $\psi$, the normalized equation becomes
$$
\frac{\partial^2\bar\psi}{\partial x^2}-\frac1x\frac{\partial\bar\psi}{\partial x}
+\frac{\partial^2\bar\psi}{\partial y^2}
= -\varphi^2\left(\frac{\tilde\beta}{2}x^2+a\bar\psi+b\right),
$$
which is decomposed into a particular solution and a homogeneous part with eigenvalue parameter $s$ [1409.1113]. The full flux is then written as a linear combination of independent basis solutions whose coefficients are fitted under physical constraints such as aspect ratio, elongation, triangularity, X-point geometry, and current-profile conditions.

In heterogeneous and anisotropic elasticity, the semi-analytical forward solver is the direct method of lines [2506.15185]. A star-shaped domain is mapped to log-polar coordinates,
$$
x=e^\rho \tilde r(\phi)\cos\phi,\qquad y=e^\rho \tilde r(\phi)\sin\phi,
$$
only the angular variable is discretized with a periodic $P_1$ finite-element space, and the static Navier equations reduce to a second-order ODE system in $\rho$,
$$
B_2 U''(\rho)+B_1 U'(\rho)+B_0 U(\rho)=0.
$$
Seeking modes of the form
$$
U(\rho)=e^{\rho\gamma}\xi
$$
yields a quadratic eigenvalue problem
$$
[\gamma^2 B_2+\gamma B_1+B_0]\xi=0,
$$
and the bounded solution is reconstructed from the eigenmodes with nonnegative real part [2506.15185]. The paper emphasizes that the eigen-exponents converge to the true stress-singularity exponents and that no special singular elements are required.

For the Noh problem with a black-box equation of state, the semi-analytic object is a self-similar two-region shock solution [2509.16766]. The pre-shock solution is
$$
\rho_R(x,t)=\rho_0\left(1-\frac{u_0 t}{x}\right)^m,\qquad P_R=P_0,\qquad u_R=u_0,
$$
the post-shock state is constant with $u_L=0$, and Rankine–Hugoniot conditions reduce the closure to a nonlinear $3\times3$ system in $(\rho_L,P_L,D)$ or $(\rho_L,e_L,D)$ depending on the EoS representation. A Newton iteration then produces the shock speed and post-shock state, from which the full piecewise semi-analytic solution follows [2509.16766].

## 5. Computational characteristics and reported performance

The principal computational claim of SAS methods is not that they eliminate numerics, but that they shift expensive computation away from the most frequent operations. In the ADM-based power-system solver, once offline derivation is done, online costs are only closed-form evaluations and no iterative solves; the method also exhibits a high degree of parallelism because each generator is independent [1412.0904]. In the DTM formulation, the avoidance of complete symbolic expansion reduces derivation and storage burdens, and coefficients at the same order but different machines are independent, which enables parallel evaluation [1802.04134]. In the time-power-series formulation, the computational advantage comes from evaluating component-level SASs over large windows while solving the network algebraic equations far fewer times [1804.07330].

Several papers report representative benchmarks:

| Domain | Representative setup | Reported outcome |
|---|---|---|
| Power systems (ADM) | IEEE 10-machine, 39-bus; order $N=3$; $T=0.04$ s | online $\approx 4.5$ ms per 4 s simulation in Case A; speedup $\sim 40\times$ real time [1412.0904] |
| Power systems (DTM) | IEEE 10-machine 39-bus; $K=12$; $t_w=0.20$ s | serial SAS time $0.118$ s; ideal 40-core time $\approx 0.00295$ s [1802.04134] |
| Power systems (time-power series) | Polish 2383-bus, 327 machines + 366 motors | ideal component-level parallelism gives total $\approx 6.9$ s; $\approx 7.9\times$ speed-up vs. Euler [1804.07330] |
| Orbit propagation | STELA vs THALASSA over decades | STELA $\sim 2\times$ faster for LEO at $\sim 10$ km accuracy; non-averaged methods comparable or faster in GTO/HEO [1810.11418] |
| Natural gas | 6-node network, $\Delta T=0.2$ s | SAS-2 with $M=3$ gives $\log_{10}(\mathrm{ERR})=-3.13$ in $3.84$ s [2210.03298] |

The orbit-propagation comparison is especially instructive because it isolates the cost of the semi-analytical approximation itself. The paper decomposes total integration error into dynamical, model truncation, short-periodic, and numerical components, and argues that the first three are attributable to distinct approximations required by the method of averaging [1810.11418]. As a result, STELA can be very fast in LEO and MEO regimes, but its attainable accuracy is fundamentally limited by averaging-inherent errors. For the test cases considered, regularized non-averaged methods were found to be up to two times slower than semi-analytical for the LEO orbit, to have comparable speed for the GTO, and to be ten times as fast for the HEO for the same accuracy [1810.11418].

The natural-gas results show a different mechanism. There the gain comes from using a much larger grid cell than the widely used finite difference method. In the single-pipeline example, SAS-2 with $M=4$ reaches error $10^{-3}$ in $0.116$ s, whereas the reference FDM with $\Delta L=1$ m and $\Delta T=10^{-3}$ s takes $0.185$ s to simulate $200$ s [2210.03298]. This is not merely a faster linear solve; it is a different discretization regime enabled by the cellwise polynomial ansatz.

## 6. Limitations, misconceptions, and methodological significance

The literature is explicit that SAS methods have domain-specific limits. In power systems, the SAS is only locally accurate; one must choose the window carefully or adaptively monitor the divergence indicator, and high-order SAS increases expression complexity [1412.0904]. Offline symbolic derivation can also become expensive as the number of symbolic parameters grows. The time-power-series formulation similarly notes that offline derivation can be memory- and CPU-intensive for very large-order series and that network algebraic equations still have to be solved, even if much fewer times [1804.07330]. In the DTM setting, larger order increases both the admissible time window and the evaluation burden, so optimal performance requires order selection rather than monotonically increasing truncation order [1802.04134].

In orbit propagation, the main limitation is structural rather than implementation-specific. Dynamical error, model truncation error, and short-periodic reconstruction error cannot be reduced by simply decreasing the numerical step size; one must derive higher-order averaging, include more harmonics, or switch to non-averaged methods [1810.11418]. Near tesseral or third-body resonances, double averaging is required, and highly eccentric or high-altitude orbits exhibit large long-periodic terms and large lunar perturbations that increase sensitivity to truncation choices.

In the elasticity solver, the method requires a star-shaped domain, piecewise-constant material coefficients in the angular sectors, and solution of a generalized eigenproblem whose cost can dominate when the angular discretization is large [2506.15185]. In the Grad–Shafranov solver, linearity depends on specific profile assumptions for pressure and current functions [1409.1113]. In the Noh-problem solver, the classic two-region structure imposes a geometry-specific restriction: in cylindrical and spherical symmetry one must have $P_0=0$ and $e(\rho\to0,P)\to\text{constant}$ to retain the classic Noh solution structure [2509.16766].

A second misconception is that “semi-analytical” implies “more accurate than numerical integration” in a universal sense. The cited orbit-propagation paper directly contradicts that generalization: semi-analytical solvers may be faster in some regimes, but non-averaged regularized formulations can match or exceed both speed and accuracy in others [1810.11418]. Conversely, a purely numerical baseline is not always the computational optimum when closed-form local propagation, averaging, or coordinatewise analytic resolution is available [1412.0904; 2210.03298; 2506.15185].

The broader significance of SAS methods lies in this trade-off architecture. They preserve enough analytical structure to reduce repeated solves, exploit parallelism, enlarge time or space cells, or encode singular behavior in basis functions, while accepting locality, truncation, or model specialization as the price of efficiency. This suggests that SAS is less a single solver family than a recurrent design principle: isolate a tractable substructure, make it explicit, and use that explicit representation to reorganize the dominant computational cost.

Source: https://www.emergentmind.com/topics/semi-analytical-solver-sas