64
(POD) provides a solution to this problem via the singular value decomposition (SVD)
(Liang et al., 2002).
To obtain the POD basis we first take “snapshots” of the voltage and nonlinear
terms at specific (usually equally-spaced) time points during the simulation. For
convenience we denote the nonlinear term by
N(v(t),w(t)) ≡ (Φ(w(t))e).v(t) -Φ(w(t))Ei.
(3.10)
We save the values of v and N(v, w) at times tɪ, t2, ∙ ∙ ∙ > tn, where tj = j∆t. The
snapshots are stored column-wise in matrices
V = [v(tι) v(t2) ∙ ∙ ∙ v(tτι)]
F = [N(v(t1),w(tι)) N(v(t2),w(t2)) ∙∙∙ N(v(tn),w(tn))].
The matrix V ∈ Rλz×" will be used to build the POD basis, while F ∈ R7v×" will be
used in §3.2.2.
We begin with the SVDs of the snapshot matrices
V = UΣXτ, F = WXYτ,
where UτU = /, XτX = I, WτW7 = I, YτY = I, and Σ and Λ are diagonal
and non-negative. These diagonal elements are ordered in a descending fashion. We