A numerical framework for computing steady states of structured population models and their stability

Figure 1.  Results-of-the Results of the numerical simulations. a) $a,\,b\,,c$ values satisfying the necessary condition (12), form a 3D surface (blue surface). Steady states of the Sinko-Streifer model only exist on the red line b) Comparison of exact stationary solution (for the point marked with red star in Figure 1a) with approximate stationary solution for $n=100$. c) Absolute error between exact stationary solution and approximate stationary solution decays linearly as the dimension of approximate subspaces $\mathcal{X}_{n}$ increase.

Figure 2.  Existence and stability regions for the steady states of the PBE a) Existence region for the steady states of the PBE forms a wedge like shape. b) Stability region for $b=0.1$, $a\in[0,\,15]$ and $c\in[0,\,5]$. c) Stability region for $b=0.5$, $a\in[0,\,15]$ and $c\in[0,\,5]$. d) Stability region for $b=1.0$, $a\in[0,\,15]$ and $c\in[0,\,5]$. Color bar represents the real part of rightmost eigenvalue of the Jacobian matrix evaluated at each steady state. Yellow regions represents the region for which a positive steady state does not exists.

Figure 3.  An example steady-state solution of the PBE for $b=0.5,\,a=c=1$. b) Steady states for increasing renewal rate and $b=c=1$

Figure 4.  Time evolution of the flocculation model with arbitrary initial conditions. a) Four different initial conditions are chosen close to the steady state. b) Solution of the PBE for those initial conditions at $t=10$. c) Evolution of the total number $M_{0}(t)$ of the flocs for $t\in[0,\,10]$. d) Evolution of the total mass $M_{1}(t)$ of the flocs for $t\in[0,\,10]$.

Figure 5.  Change in zeroth and first moments with increasing dimension of the approximate space $\mathcal{X}_{n}$. a) Change in the total number and the total mass of the flocs with respect to increasing dimension $n$. Dashed red lines and dotted green lines corresponds to the total number and the total mass of the flocs of the steady state for $n=1000$, respectively. b) Steady state solution for $n=100$ and $n=500$.

Figure 6.  Eigenvalues of the Jacobian $J_{\mathcal{F}}(\alpha)$ multiplied by $\Delta x$ for the steady state illustrated in Figure 3a. a) Eigenvalues of the Jacobian plotted in the complex plane for $n=20$. b) Eigenvalues of the Jacobian plotted in the complex plane for $n=50$. c) Eigenvalues of the Jacobian plotted in the complex plane for $n=200$. d) Change in the rightmost eigenvalue for increasing $n$. Dashed black line corresponds to the rightmost eigenvalue of the Jacobian for $n=1000$.