23-Chem-A6 Process Dynamics and Control · May 2016
Question 7 of 8: Nyquist Stability of an Open-Loop-Unstable Process
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Reference texts: D. E. Seborg, T. F. Edgar, D. A. Mellichamp & F. J. Doyle III, Process Dynamics and Control (4th ed., Wiley) — Laplace-domain modelling, linearisation, transfer functions, feedback stability, frequency response and IMC design; G. Stephanopoulos, Chemical Process Control: An Introduction to Theory and Practice (Prentice Hall) — Routh test, Nyquist criterion and dead-time systems; D. R. Coughanowr & S. E. LeBlanc, Process Systems Analysis and Control (3rd ed., McGraw-Hill) — first-order thermal/level dynamics, step/ramp/pulse response and block-diagram algebra. Standard control conventions (deviation variables; unity valve/sensor gains unless stated) are used throughout. This sitting is an open-book paper of eight problems of which any five constitute a complete exam; all eight are solved here.
Problem 7: Nyquist Stability of an Open-Loop-Unstable Process (20%)
Given. $G_p(s)=\dfrac{100}{s-10}$ — an open-loop-unstable first-order process (one pole at $s=+10$, so $P=1$ RHP pole). Proportional controller $G_c=K_c$; open loop $L(s)=K_c G_p=\dfrac{100K_c}{s-10}$.
Find. (a) Nyquist verdict at $K_c=1,\ 0.01$; (b) the stability limit on $K_c$.
Problem 7: Nyquist plot of $L=100/(s-10)$ at $K_c=1$ — a circle from $L(0)=-10$ to the origin. The critical point $-1$ lies inside, giving one counter-clockwise encirclement; with $P=1$ open-loop RHP pole, $Z=N+P=-1+1=0$, so the loop is stable. At $K_c=0.01$ the circle only reaches $-0.1$, so $-1$ is not encircled ($Z=1$, unstable).
Approach. Apply the Nyquist criterion $Z=N+P$, where $Z$ is the number of closed-loop RHP poles, $P$ the open-loop RHP poles, and $N$ the clockwise encirclements of $-1$ by $L(j\omega)$; because $P=1$, stability ($Z=0$) requires exactly one counter-clockwise encirclement of $-1$. Cross-check against the closed-loop pole directly.
Shape of the Nyquist plot. $L(j\omega)=\dfrac{100K_c}{j\omega-10}=\dfrac{100K_c(-10-j\omega)}{100+\omega^2}$. The real part is always negative and the locus is a circle through the origin ($\omega\to\infty$) and through $L(0)=-10K_c$ ($\omega=0$) — a circle on the real axis from $-10K_c$ to $0$.
(a) $K_c=1$. The circle runs from $-10$ to $0$; the critical point $-1$ lies inside it, so $L$ makes one counter-clockwise encirclement: $N=-1$. Then $Z=N+P=-1+1=0$ — stable. Direct check: $s-10+100K_c=0\Rightarrow s=10-100(1)=\boxed{-90}$ (LHP). ✓
(a) $K_c=0.01$. The circle now runs only from $-10(0.01)=-0.1$ to $0$, so $-1$ lies outside it: no encirclement, $N=0$. Then $Z=N+P=0+1=1$ — unstable (one closed-loop RHP pole). Direct check: $s=10-100(0.01)=\boxed{+9}$ (RHP). ✓
(b) Limiting gain. Encirclement of $-1$ (hence $Z=0$) requires the circle to extend past $-1$, i.e. $10K_c>1$. The boundary is $L(0)=-10K_c=-1$: $$\boxed{K_{c,\lim}=0.1}$$ The system is stable for $K_c>0.1$ and unstable for $K_c<0.1$ — note the reverse sense: an open-loop-unstable process needs enough gain to be stabilised. (Direct: $s=10-100K_c<0\Leftrightarrow K_c>0.1$.)