Supplementary material from "Analysing diffusion and flow-driven instability using semidefinite programming"

However, verifying the transport-driven instability of reaction–diffusion–advection systems requires checking the Jacobian eigenvalues of infinitely many Fourier modes, which is computationally intractable. To overcome this limitation, this paper proposes mathematical optimization algorithms that determine the stability/instability of reaction–diffusion–advection systems by finite steps of algebraic calculations.

Specifically, the stability/instability analysis of Fourier modes is formulated as a sum-of-squares optimization program, which is a class of convex optimization whose solvers are widely available as software packages. The optimization program is further extended for facile computation of the destabilizing spatial modes. This extension allows for predicting and designing the shape of the concentration gradient without simulating the governing equations.