NivaarExam PrepOfficial exam papers ↗

23-Chem-A6 Process Dynamics and Control · May 2015

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

National Exams / EGBC — May 2015 — 04-Chem-A6 Process Dynamics & Control. Three-hour open-book examination; any non-communicating calculator is permitted. Eight problems are printed and any five constitute a complete paper (each worth 20%); all eight are solved below for completeness. Most parts are quantitative (dynamic modelling, transfer functions, step/pulse/ramp responses, Routh and Nyquist stability, and IMC design); qualitative sketches are drawn as real figures where the paper asks for them.

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.

Problem 7: Nyquist Stability of an Open-Loop-Unstable Process (20%)

Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.

Given. $G_p=\dfrac{1}{s^2-s-2}=\dfrac{1}{(s-2)(s+1)}$ — one RHP pole at $s=+2$ (so $P=1$) and one LHP pole at $s=-1$; open loop $L(s)=k_cG_p$.

Find. (a) the Nyquist verdict at $k_c=1$; (b) the stabilising range of $k_c$.

Approach. Evaluate $L(j\omega)$ at low, mid and high frequency to fix the locus shape, count encirclements of $-1$, apply $Z=N+P$, and cross-check against the characteristic polynomial by the Routh test.

  1. Key points of $L(j\omega)$. $L(j\omega)=\dfrac{k_c}{(j\omega)^2-j\omega-2}=\dfrac{k_c}{-(\omega^2+2)-j\omega}$. At $\omega=0$: $L=\dfrac{k_c}{-2}=-\dfrac{k_c}{2}$ (for $k_c=1$, the point $-0.5$). As $\omega\to\infty$: $L\to0$ (approaching from the upper half-plane, since $\mathrm{Im}\,L=\dfrac{k_c\omega}{(\omega^2+2)^2+\omega^2}>0$). The real part is negative for all $\omega$, so the locus is a small arc confined to the left half-plane, running from $-0.5$ up to the origin.
  2. (a) Encirclement count for $k_c=1$. The locus reaches only to $-0.5$ on the real axis — it stays to the right of the critical point $-1$ and therefore does not encircle it: $N=0$. With $P=1$ open-loop RHP pole, $$Z=N+P=0+1=1\neq0,$$ so there is one closed-loop RHP pole — the loop is unstable.
  3. Direct check ($k_c=1$). Characteristic equation $s^2-s-2+k_c=s^2-s-1=0$ has roots $s=\dfrac{1\pm\sqrt5}{2}=+1.618,\,-0.618$: one root in the RHP — unstable, confirming $Z=1$. ✓
  4. (b) Range of $k_c$ — Routh test. General characteristic equation: $s^2-s+(k_c-2)=0$. The Routh array is $$\begin{array}{c|cc}s^2&1&k_c-2\\ s^1&-1&0\\ s^0&k_c-2&\end{array}$$ The first column contains $-1<0$ for every $k_c$, so there is always a sign change — at least one RHP root regardless of gain.
  5. (b) Conclusion. $$\boxed{\text{No value of }k_c\text{ stabilises the loop.}}$$ Equivalently, the sum of the closed-loop roots equals $+1>0$ (from the $-s$ term, independent of $k_c$), so the roots can never both lie in the LHP. A single proportional gain cannot move both poles left; stabilising this plant requires a controller that adds phase/zeros (e.g. PD or a lead compensator).
ReIm(−1, 0)ω=0 (-0.50)P=1 RHP pole; locus does NOT encircle −1 ⇒ N=0, Z=N+P=1 (unstable)
Problem 7(a): Nyquist plot of $L=1/(s^2-s-2)$ at $k_c=1$. The locus (in the left half-plane) runs from $-0.5$ to the origin and does not encircle $-1$, so $N=0$; with $P=1$, $Z=N+P=1$ — unstable.
ItemResult
Open-loop RHP poles$P=1$ (at $s=+2$)
$L(0)$ at $k_c=1$$-0.5$ (does not reach $-1$)
Encirclements of $-1$ / $Z$$N=0$, $Z=N+P=1$ → unstable
Stabilising range of $k_c$none (Routh first column has $-1$)