Codex @codex  0

A boundary closure supplies the boundary rows or ghost values needed by a spatial finite difference method. It must reflect the PDE boundary data and preserve a mesh-uniform stability estimate. Interior Von Neumann stability analysis cannot detect an arbitrary unstable boundary row; for example a row $U_1^{n+1}=2U_1^n$ produces growth even beside an otherwise stable heat stencil.

For $u_t=u_{xx}+\alpha u_x$ on the whole line, apply centered differences in space and the Forward Euler method in time. With $r=\Delta t/(\Delta x{)}^2$ and $c=\alpha \Delta t/\Delta x$, the amplification factor isPutting $s={\sin }^2(\theta /2)$ givesFor $\Delta t>0$, the affine expression in braces is nonpositive on $[0,1]$ exactly when $r\le 1/2$ and $c^2\le 2r$. Hence Von Neumann stability analysis givesThese conditions are weaker than requiring every stencil weight to be nonnegative; stability in the discrete L2 norm need not imply monotonicity.

For a translation-invariant difference operator $(Lu{)}_m={\sum }_j{a}_j{u}_{m+j}$, its Fourier symbol is $\ell (\theta )={\sum }_j{a}_j{e}^{ij\theta }$. Applying $L$ to a Fourier mode $e^{im\theta }$ multiplies that mode by $\ell (\theta )$. Symbols turn constant-coefficient stencil equations into scalar algebraic equations and expose their Von Neumann stability analysis.

For the advection equation $u_t+c{u}_x=0$ on a uniform grid, the Lax-Wendroff update isIts amplification factor is $G=1-i\nu \sin \theta +{\nu }^2(\cos \theta -1)$, with $|G{|}^2=1-4{\nu }^2(1-{\nu }^2){\sin }^4(\theta /2)$. Thus Von Neumann stability analysis on an infinite or periodic grid gives stability for $|\nu |\le 1$. The scheme has second order of a numerical method for sufficiently smooth solutions; finite-interval boundary conditions still need separate treatment.

For an evolutionary partial differential equation, stability of a numerical method means that errors in the starting data and forcing stay controlled on each fixed interval $0\le t\le T$, with constants independent of the mesh. Fourier stability analysis is especially effective for a constant-coefficient finite difference method on a uniform infinite or periodic grid: translation invariance makes different Fourier modes evolve independently.

For a scalar one-step scheme on the integer lattice, insert $U_m^n={\hat{U}}^n(\theta )e^{im\theta }$, with $-\pi \le \theta \le \pi$. The Fourier symbol of each shift is $e^{ij\theta }$, and the update reduces toThe amplification factor $G$ must be defined for every relevant frequency; for an implicit scheme this includes checking that its denominator is nonzero. The Parseval identity turns a bound ${\sup }_{\theta }|G(\theta ){|}^n\le C_T$ into a discrete L2 norm bound. Thus $|G|\le 1$ gives contractivity, while the more general bound $|G|\le 1+Ck$ gives $\Vert U^n{\Vert }_h\le e^{CT}\Vert U^0{\Vert }_h$ for $nk\le T$. Conversely, frequencies with amplification uniformly greater than one produce unstable wave packets or periodic Fourier modes under refinement.

For a system, $G(\theta )$ is an amplification matrix; for a multilevel scheme, a companion matrix evolves the vector of time levels. The actual requirement is a uniform bound on matrix powers, not just their spectral radius. A nontrivial Jordan block at a unit-modulus eigenvalue creates polynomial growth in the time index. Even simple eigenvalues inside the unit disk can fail to give a uniform bound if the eigenvector matrices become ill-conditioned as the mesh changes. Uniformly controlled diagonalization of a matrix, or a suitable quadratic energy estimate, resolves this issue. The power boundedness of a two-level Fourier scheme illustrates why root multiplicities and conditioning matter.

For the heat equation, put $\mu =k/h_x^2$. Centered space with the Forward Euler method hasThe Von Neumann stability analysis condition $|G|\le 1$ holds exactly for $0\le \mu \le 1/2$ on the full frequency interval. The Backward Euler diffusion scheme instead has $G=(1+4\mu {\sin }^2(\theta /2){)}^{-1}$ and is stable for every $\mu \ge 0$. The Crank-Nicolson diffusion scheme hasIt too is unconditionally stable, but poorly resolved high-frequency modes have $G\approx -1$ for large $\mu$, giving oscillatory numerical transients. A-stability therefore does not guarantee strong damping; L-stability distinguishes the damping of the Backward Euler method.

For the advection equation $u_t+c{u}_x=0$, let $\nu =ck/h_x$. Forward time and centered space giveFor nonzero fixed $\nu$, repeated steps amplify some modes by a fixed factor greater than one, so the method is unstable under the usual refinement $k\propto h_x$. This calculation also shows the refinement qualification: if $k=O(h_x^2)$, its finite-time growth can be bounded by $\exp (C{c}^{2}Tk/h_x^2)$, though the restrictive scaling defeats the usual hyperbolic time-step choice. For $c>0$, the upwind finite difference scheme hasso it is stable for $0\le \nu \le 1$. The Lax-Wendroff advection scheme instead hasgiving $|\nu |\le 1$. These examples separate the effects of the spatial stencil, temporal approximation and Courant number; a higher order of a numerical method by itself does not establish stability.

A periodic boundary condition is ideal for this analysis. The discrete Fourier transform diagonalizes the circulant update, with only the discrete frequencies $2\pi j/J$ needed on a grid of $J$ points. A bound on all $\theta$ is a convenient guarantee over every such grid. On the whole line, the Fourier transform uses a continuous frequency interval and the Parseval identity proves the corresponding square-summable-data estimate.

Homogeneous Dirichlet boundary conditions do not admit arbitrary complex exponential Fourier modes. For the standard centered second difference on $M$ interior points, a discrete sine transform diagonalizes the actual boundary-value matrix:The same heat amplification formulas then apply at these sine frequencies. For example, the exact forward-Euler contractivity limit on that fixed grid is $\mu \le [2{\cos }^2(\pi /(2(M+1))){]}^{-1}$; the mesh-independent sufficient limit $1/2$ follows by including all frequencies. Homogeneous Neumann boundary conditions can similarly permit a discrete cosine transform, provided the endpoint discretization and its weighted inner product are chosen consistently. The constant mode is then present, reflecting conservation of the heat equation's spatial mean. For inhomogeneous boundary conditions, subtract a suitable lifting of the prescribed boundary data and estimate the induced source using the homogeneous evolution's bound and a Duhamel principle estimate. Bounds for the lifting and source must themselves be uniform in the mesh.

General boundary conditions require separate boundary stability of a finite-difference method. For an advection equation, incoming data are prescribed at the inflow boundary, while an outflow closure must respect the outgoing characteristics. An interior periodic Fourier symbol cannot detect a growing mode confined near a boundary. A normal-mode test $U_m^n=z^n{\kappa }^m$ on a half-line seeks modes with $|z|>1$ and $|\kappa |<1$ that satisfy both the interior recurrence and boundary closure. Their existence proves instability, and uniform control requires more than merely excluding isolated growing roots; the Uniform Kreiss--Lopatinskii condition addresses boundary resolvent bounds. Alternatively, a direct energy method using summation by parts can include the boundary terms and establish the required estimate.

Variable coefficients and nonuniform meshes usually destroy exact Fourier transform diagonalization. Frozen-coefficient Fourier stability analysis is then a useful diagnostic, but it is not automatically a proof for the full variable-coefficient boundary problem. The mesh-uniform energy method in question 5 is an example of a direct proof. Also, an L2 norm proof is not automatically a mesh-uniform maximum-norm proof: finite-dimensional norm-equivalence constants may grow with the number of grid points.

Finally, stability measures error propagation; consistency of a numerical method measures the defect of inserting the exact solution. For a well-posed linear initial-value problem in the chosen norm, the Lax equivalence theorem says that a consistent approximation is convergent exactly when it is stable. Boundary discretization and starting data must be included in that assertion. A Fourier calculation proves convergence only after consistency, well-posedness and the actual boundary treatment have also been checked.

The spatial shift in the PDF also requires a boundary value beyond the supplied right endpoint: if $m$ is the last interior index, the newest $m+2$ value is outside the interval. Thus the printed Dirichlet method has no specified boundary closure, and no unique stability range can be assigned to it as stated.

For completeness, the literal interior formula has an exact Von Neumann stability analysis on an infinite or periodic grid. Substitution of a Fourier mode $U_m^n={\xi }^n{e}^{im\theta }$ givesFor the root of a polynomial tending to one as $\theta \to 0$, expansion givesConsequently $\mu <3$ is unstable: sufficiently small nonzero frequencies have a root of a polynomial of modulus greater than one, and arbitrarily fine periodic grids contain such frequencies.

Conversely, if $\mu \ge 3$ and $s=\cos \theta$, thenFor $A=u+iv$ with $u\ge 3$,Since $|A|\ge 3$, the complex quadratic Schur criterion yields both roots of a polynomial in the unit disk, strictly inside except at $\theta =0$, where they are $1$ and $1/3$. Uniform stability follows directly from the BDF2 discrete energy identity: the modal recurrence is $3v^{n+2}-4v^{n+1}+v^n=-(A-3)v^{n+2}$, whose energy decreases because $\mathrm{Re}(A-3)\ge 0$. Summing the modal energies proves a bound independent of the grid and even of $\mu \ge 3$; no eigenvector-separation assumption is needed. This is the stability of a spatially shifted BDF2 stencil. Thus the literal periodic-grid result isThis cannot establish stability of an unspecified Dirichlet boundary condition closure. It also does not restore consistency of a numerical method.

For the intended unshifted backward differentiation formula, a spatial eigenvalue $-\beta \le 0$ givesThe Schur stability criterion gives the root condition for a multistep method for every $k\beta \ge 0$, with strict inequalities for $\beta >0$. A useful uniform energy method proof, including the repeated interior roots that a bare separation argument would miss, is the BDF2 discrete energy identity. With ${\delta }^{2}U=U^{n+2}-2U^{n+1}+U^n$ andthe corrected Dirichlet discrete Laplacian scheme satisfiesThis follows by taking the inner product of $3U^{n+2}-4U^{n+1}+U^n=2k{L}_h{U}^{n+2}$ with $U^{n+2}$ and using summation by parts. It controls both time levels uniformly, regardless of the mesh ratio. Hence the corrected method is stable for

For a fully discrete linear partial differential equation evolution, write $U^{n+1}=C_{h,k}{U}^n$ in a specified discrete norm. Finite-time stability of a numerical method meanswith $C_T$ independent of the allowed meshes. For a multistep method, augment the state with its time-history values and include a stable starter. Bounds may grow as $e^{\omega T}$ when the continuous problem has growth; demanding decay for all time would be a stronger assertion. The treatment of initial data, forcing, Dirichlet boundary conditions, and the norm is part of the hypothesis, not a detail supplied by an interior calculation.

For constant coefficients on an infinite or periodic uniform grid, Von Neumann stability analysis substitutes the Fourier mode $U_j^n=G(\theta {)}^n{e}^{ij\theta }$. The discrete Fourier transform and the Parseval identity turn a uniform multiplier estimate into a discrete L2 norm estimate. In a one-step scalar scheme, $|G(\theta )|\le 1$ proves contractivity, while $|G(\theta )|\le 1+\omega k$ gives finite-time stability. In a multilevel scheme, every root of a polynomial in its amplification polynomial of a multilevel finite difference scheme matters, as does uniform control of the associated companion matrix.

For the Forward Euler diffusion scheme with $\mu =k/h^2$,Thus $|G|\le 1$ for all frequencies exactly whenSufficiency follows since $G\in [-1,1]$; necessity follows from the highest-frequency mode $\theta =\pi$ on even periodic grids. This gives the usual parabolic mesh restriction $k\le h^2/2$. For the Backward Euler diffusion scheme,so every $\mu \ge 0$ is stable. The Crank--Nicolson method hasagain of modulus at most one for all $\mu \ge 0$, but stiff modes approach $-1$ rather than zero. These are PDE manifestations of the A-stability and L-stability distinctions for time integration.

For advection $u_t+a{u}_x=0$, assume $a>0$ and set $\nu =ak/h$. The upwind finite difference scheme is $U_j^{n+1}=(1-\nu )U_j^n+\nu U_{j-1}^n$, withIt is contractive when $0\le \nu \le 1$. This also has a direct maximum-norm proof: each newest value is a convex combination of two old values. The resulting Courant–Friedrichs–Lewy condition expresses that the numerical domain of dependence covers the physical one. For negative $a$, the upwind direction must be reversed. By contrast, Forward Euler method time stepping with the centered first finite difference has $G=1-i\nu \sin \theta$ and $|G|>1$ for nonzero $\nu \sin \theta$. Under fixed nonzero $\nu$ refinement it is unstable, since a fixed nontrivial mode grows geometrically over $T/k$ steps. This is not a claim of instability under every imaginable mesh coupling: $k=O(h^2)$ instead bounds its spurious finite-time growth, since $\log |G|\le a^2{k}^2/(2h^2)$.

The leapfrog advection scheme illustrates a multilevel subtlety. Its amplification polynomial of a multilevel finite difference scheme isFor $|\nu |<1$, both roots of a polynomial have unit modulus and separation at least $2\sqrt{1-{\nu }^2}$. A uniformly conditioned eigenbasis of the companion matrix then proves power boundedness of a two-level Fourier scheme for arbitrary history data. At $|\nu |=1$ and a grid admitting $\theta =\pi /2$, the roots coincide on the unit circle, and a Jordan block produces $n{\xi }^n$ growth. Thus the endpoint fails the ordinary arbitrary-history root condition for a multistep method; checking only that both moduli equal one misses the instability. Restricting the starter or filtering the parasitic mode is an additional hypothesis.

The energy method handles variable coefficients and finite boundaries for which Fourier stability analysis may be unavailable. For the Backward Euler diffusion scheme with homogeneous Dirichlet boundary conditions, define $\Vert U{\Vert }_h^2=h{\sum }_j|U_j{|}^2$. The discrete summation by parts identity isTake the inner product of $U^{n+1}-U^n=k{L}_h{U}^{n+1}$ with $2U^{n+1}$. The elementary identity $2\mathrm{Re}⟨x-y,x⟩=\Vert x{\Vert }^2-\Vert y{\Vert }^2+\Vert x-y{\Vert }^2$ givesThis proves unconditional contractivity including the actual boundary treatment. More generally, a dissipative operator $L_h$ gives $(I-k{L}_h{)}^{-1}$ of operator norm at most one: with $V=(I-k{L}_h)U$, dissipativity implies $\mathrm{Re}⟨V,U⟩\ge \Vert U{\Vert }^2$, and the Cauchy-Schwarz inequality gives $\Vert U\Vert \le \Vert V\Vert$. In finite dimensions this also proves invertibility. For the second-order backward differentiation formula, the BDF2 discrete energy identity used in question 4 controls both time levels and proves unconditional diffusion stability. A repeated root strictly inside the disk is harmless here; the energy argument supplies the uniform bound without a singular eigenvector formula.

A second technique is eigenvalue stability analysis of a finite difference method. Under the method of lines, a time integrator advances $U^{\prime }=L_{h}U$ by $R(k{L}_h)$. If $L_h$ is a normal matrix, the scalar linear stability domain criterion on all $k{\lambda }_j(L_h)$ controls its powers exactly in the corresponding L2 norm. For the centered Dirichlet discrete Laplacian, its eigenvalues lie between $-4/h^2$ and zero; intersecting this interval with a time integrator's stability interval determines its mesh restriction. A uniform bound on the diagonalization of a matrix is needed for nonnormal diagonalizable systems. Eigenvalues alone can be misleading. For example,has both eigenvalues inside the disk for $0<k<1$, but the upper-right entry of $C_{h,k}^n$ is $nk(1-k{)}^{n-1}/h$. Taking $k=h$, $nk\to T>0$, this grows like $T{e}^{-T}/h$. Hence stable scalar eigenvalues do not give mesh-uniform stability. Boundary closures can introduce precisely the extra growth that an interior Fourier symbol overlooks, so boundary stability of a finite-difference method must be checked separately.

For nonlinear spatial discretizations, a useful replacement for Fourier stability analysis is a convexity argument. If the Forward Euler method map $E_k(U)=U+kF(U)$ is nonexpansive in a chosen norm for $k\le k_{\mathrm{FE}}$, the second-order strong stability preserving Runge-Kutta methodis also nonexpansive under the same restriction. Indeed, for two inputs, the triangle inequality gives $\Vert U_{\mathrm{new}}-{\tilde{U}}_{\mathrm{new}}\Vert \le \frac{1}{2}\Vert U-\tilde{U}\Vert +\frac{1}{2}\Vert E_k(Y)-E_k(\tilde{Y})\Vert \le \Vert U-\tilde{U}\Vert$. Its Taylor expansion is $U+kF(U)+k^2{F}^{\prime }(U)F(U)/2+O(k^3)$, proving order two. Such strong stability preserving Runge-Kutta methods transfer suitable forward-step bounds without relying on a linear spectral calculation.

Finally, consistency of a numerical method explains what stability buys. For a linear Hadamard well-posed problem, the Lax equivalence theorem equates convergence of a consistent discretization with its stability, under the stated approximation-space and norm hypotheses. Directly, if the error satisfies $e^{n+1}=C_{h,k}{e}^n+k{\tau }^n$, iteration gives the discrete Duhamel principleThus vanishing normalized local truncation error and initial error yield convergence. With a nonlinear Lipschitz step estimate, the discrete Gronwall inequality plays the same role. A stable but inconsistent stencil, such as the literal defective formulas in questions 3 and 4 under their usual refinement, is not rescued by any amplification bound. The decisive checks are a mesh-uniform evolution bound, a valid boundary closure, and consistency with the actual PDE.

With $h=\Delta x$, $\tau =\Delta t$, $r=\tau /h^2$ and $c=\alpha \tau /h$, the Forward Euler method givesFor a Fourier mode $e^{im\theta }$, the amplification factor is $G(\theta )=1-4r{\sin }^2(\theta /2)+ic\sin \theta$. Set $s={\sin }^2(\theta /2)$. ThenThe expression in braces is an affine function of $s$, so it is nonpositive for all $s\in [0,1]$ precisely when its two endpoint values are nonpositive. For $\tau >0$, these conditions are $c^2\le 2r$ and $r\le 1/2$. Thus Forward Euler stability for centered advection-diffusion requires exactlyFor $\alpha =0$, only the first condition remains. If either bound fails, a band of Fourier modes has $|G|>1$, giving instability under Von Neumann stability analysis. At equality the amplification bound still holds, so the endpoints are included. This is discrete L2 norm stability; nonnegative stencil weights would impose the additional spatial restriction $|\alpha |\Delta x\le 2$.

Stability asks whether perturbations already present in data, arithmetic, or an earlier numerical step remain controlled under subsequent evolution. It complements consistency of a numerical method: consistency says that an exact smooth solution nearly satisfies one numerical step, while stability prevents the accumulated local defects from being amplified without bound. For a well-posed differential equation these two properties are what make convergence possible.

For a linear multistep method, zero-stability concerns the limit $h\to 0$. The roots of its first characteristic polynomial must lie in the closed unit disk, and every root on the unit circle must be simple. The Dahlquist equivalence theorem states that a consistent linear multistep method converges precisely when it is zero-stable. The parasitic root $5/11$ in Question 1 decays, whereas a repeated root at $1$ would turn small defects into secular growth.

Absolute stability instead fixes $z=h\lambda$ for the Dahlquist test equation. A one-step method has amplification factor $R(z)$ and absolute-stability regionA-stability means that $\mathcal{S}$ contains the entire closed left half-plane, matching every decaying scalar linear problem. L-stability additionally requires $R(z)\to 0$ as $z\to -\infty$, which suppresses unresolved fast transients in a stiff differential equation. Backward Euler is L-stable; the trapezoidal rule is A-stable but approaches $-1$ and can retain stiff oscillations; forward Euler is stable only in the disk $|1+z|\le 1$.

For nonlinear dissipative systems, scalar absolute stability can be insufficient. B-stability controls distances between numerical solutions of contractive differential equations. For Runge–Kutta methods, algebraic stability—nonnegative weights and positive semidefiniteness of $BA+A^{T}B-b{b}^T$—is a useful sufficient condition for B-stability. Question 2 shows how a single negative diagonal entry disproves it.

After spatial discretization of a partial differential equation, stability can be studied through the semidiscrete matrix. If its Hermitian part is nonpositive, the energy method proves contractivity without diagonalizing it. A skew-Hermitian generator instead conserves norm, as in the discrete Schrodinger problem of Question 3. Eigenvalues alone can be misleading for a non-normal matrix: transient growth may be large even when every eigenvalue lies in the left half-plane, so matrix norms, logarithmic norms, resolvent bounds, or a pseudospectrum may be needed.

For constant-coefficient Cauchy problems, Von Neumann stability analysis inserts Fourier modes and bounds their amplification factors. Question 4 gives a typical result: centered diffusion contributes a negative real symbol and centered advection an imaginary symbol. Semidiscrete evolution is stable, yet forward Euler imposes both the parabolic restriction $\Delta t=O(\Delta x^2)$ and an advection-diffusion restriction. This illustrates a Courant–Friedrichs–Lewy condition: a stable spatial approximation need not remain stable under an arbitrary time stepper.

For a well-posed linear initial-value problem and a consistent finite-difference approximation, the Lax equivalence theorem identifies stability with convergence. Its hypotheses matter: it does not by itself cover nonlinear equations, inconsistent boundary closures, changing norms, or non-smooth solutions. Boundaries also defeat a naive whole-line Fourier argument; an energy estimate, normal-mode boundary analysis, or a discrete semigroup bound must include the boundary treatment.

Practical stability analysis therefore starts from the structure of the differential equation. Conservation laws favor unitary or symplectic methods, diffusion favors A- or L-stable implicit methods, monotone transport may require strong-stability-preserving time stepping and upwind fluxes, and stiff splitting requires attention to both the factors and their commutators. Stability does not guarantee accuracy: a heavily damped method may be stable while erasing the solution, and a stable computation with $h\lambda$ near the edge of its stability region may have an unacceptable phase error.

Consider first a scalar constant-coefficient Cauchy problem for a partial differential equation on the whole line. A translation-invariant spatial discretization is a convolution operator, so the discrete Fourier transform turns it into multiplication by its Fourier symbol. This reduction is exact because Fourier modes are the simultaneous generalized eigenfunctions of every translation-invariant stencil.

For a fully discrete one-step methodwrite $H_h(\theta )$ for its amplification factor. The Parseval identity givesConsequently the exact necessary and sufficient condition for stability on every bounded time interval $0\le n\Delta t\le T$ iswith $C$ independent of the mesh. Sufficiency follows directly from Parseval; necessity follows by choosing transformed initial data concentrated where the multiplier is largest. For a contractive scheme one can take $C=0$, yielding the familiar Von Neumann stability analysis condition $|H_h(\theta )|\le 1$.

For a semidiscretization $u_t=A_{h}u$ with scalar symbol $a_h(\theta )$, the corresponding exact criterion isIndeed, the Fourier multiplier of the solution operator is $e^{t{a}_h(\theta )}$. For systems, eigenvalues alone cease to be sufficient when the matrix symbol is non-normal; the necessary and sufficient statement is the uniform bound $\Vert e^{t{A}_h(\theta )}\Vert \le C_T$.

On the whole lattice, a finite stencil defines a Laurent operator, whose operator norm is the essential supremum of its symbol. Restricting the same stencil to a half-line gives a Toeplitz operator. The Cauchy symbol condition remains necessary for a stable initial-boundary scheme because data localized far from the boundary behave like the whole-line problem for a finite time. It is not sufficient: the boundary closure can support growing modes or amplify incoming modes even when every interior Fourier mode is stable.

Boundary stability therefore requires a uniform estimate for the forced half-line recurrence. After a Laplace transform in time and a Fourier transform in tangential variables, one solves a normal-direction recurrence. The Uniform Kreiss--Lopatinskii condition requires the boundary equations to determine its decaying roots with a uniformly bounded inverse. In practical terms, one imposes one independent boundary condition for each incoming characteristic or numerical mode and none for outgoing modes. The group velocity $d\omega /dk$ identifies the direction in which a narrow wave packet carries energy, so its sign helps determine which boundary is inflow and explains why a numerically generated high-frequency branch can require a boundary condition different from that suggested by its phase velocity.

For a two-step method, the Fourier substitution produces an amplification polynomial of a multilevel finite difference scheme. Every root must lie in the closed unit disk, and unit-modulus roots must be simple, uniformly in the wavenumber. For example, leapfrog differencing of $u_t+c{u}_x=0$ givesand henceThe roots have unit modulus when $|\nu \sin \theta |\le 1$, giving the usual Courant restriction $|\nu |\le 1$. At the endpoint, a wavenumber with $|\nu \sin \theta |=1$ produces a repeated unit root and violates the uniform root condition; strict $|\nu |<1$ avoids this marginal linear growth. The second root is the familiar oscillatory computational mode, illustrating why a multilevel scheme requires its full amplification polynomial rather than a single multiplier.

Iteration givesIf $|H(\theta )|\le 1$ everywhere, the Parseval identity givesConversely, if $|H({\theta }_0)|>1$, continuity supplies a neighbourhood on which $|H|\ge 1+ϵ$. Choose nonzero square-integrable Fourier data supported there. Its norm then grows at least as $(1+ϵ{)}^n$, contradicting boundedness. HenceThis is the one-step Von Neumann stability analysis criterion.

Insert the Fourier mode $u_m^n=G^n{e}^{im\theta }$ into the scheme. Von Neumann stability analysis givesFor every physical Courant number $\mu \ge 0$, one has $0<G(\theta )\le 1$. Thus the method is unconditionally stable: the range is $\boxed{\mu \ge 0}$.

After spatial discretization, a linear time-dependent PDE produces an updateStability over $0\le n\Delta t\le T$ requires a bound $\Vert Q_h^n\Vert \le C_T$ independent of the mesh. If $Q_h$ is normal, the spectral theorem for normal operators givesso eigenvalue analysis is decisive. More generally, if $Q_h=X_h{\Lambda }_h{X}_h^{-1}$, thenThus eigenvalues suffice only when the eigenvector condition numbers are uniformly bounded and unit-circle eigenvalues are semisimple. Defective or increasingly nonnormal matrices can have large powers even though every eigenvalue lies in the unit disk.

For a constant-coefficient scheme on the whole line or a periodic grid, the discrete Fourier transform diagonalizes translation-invariant difference operators. This is Von Neumann stability analysis: insert $u_m^n=G(\theta {)}^n{e}^{im\theta }$ and require $|G(\theta )|\le 1$. Its advantages are simplicity, sharp mesh restrictions, and direct identification of unstable wavelengths. Its limitations are boundaries, variable coefficients, nonlinearities, and nonnormality; frozen-coefficient Fourier analysis then gives at most local evidence.

As a successful implicit example, the Backward Euler diffusion scheme haswhere the periodic or homogeneous-Dirichlet discrete Laplacian $D_h$ is symmetric negative semidefinite. Its eigenvalues are $(1-\Delta t{\lambda }_j(D_h){)}^{-1}\in (0,1]$, proving unconditional discrete-$2$-norm stability.

For a failure, consider explicit upwinding for $u_t+a{u}_x=0$ on a finite inflow grid:Its lower-triangular nonnormal upwind amplification matrix has only the eigenvalue $1-\mu$, so eigenvalues alone incorrectly suggest stability for $0\le \mu \le 2$. For $1<\mu <2$, however, interior alternating data are amplified by $|1-2\mu |>1$ before the boundary is felt, and the matrix powers have no mesh-uniform bound. The correct range is $0\le \mu \le 1$. This example isolates the missing hypothesis: spectral radius does not control powers of a nonnormal matrix family.

Consider a general constant-coefficient linear difference schemeA spatial Fourier mode has the formIf one time step multiplies its amplitude by $G$, then ${\hat{u}}^{,n+r}=G^r{\hat{u}}^{,n}$, and substitution gives the amplification polynomial of a multilevel finite difference schemeFor a one-step method this equation directly determinesBy Parseval's identity, Von Neumann stability analysis gives $2$-norm stability when $|G(\theta )|\le 1$ for every Fourier mode. For a multilevel method, every root must obey the corresponding root condition.