NivaarExam PrepOfficial exam papers ↗

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

Question 7 of 8: Nyquist Stability of an Open-Loop-Unstable First-Order Process

Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)

Notes on this paper

National Exams / EGBC — December 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 First-Order 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{20}{s-3}$ — a single right-half-plane pole at $s=+3$, so the open loop has $P=1$ unstable pole; $L(s)=k_cG_p=\dfrac{20k_c}{s-3}$.

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

Approach. Plot a few points of $L(j\omega)$ to fix the locus, count encirclements of $-1$, apply $Z=P+N$ (with $N$ the clockwise encirclements), and cross-check against the characteristic equation.

  1. Key points ($k_c=1$). $L(j\omega)=\dfrac{20}{j\omega-3}=\dfrac{20(-3-j\omega)}{9+\omega^2}$. At $\omega=0$: $L=-\tfrac{20}{3}=-6.67$ (negative real axis). At $\omega=3$: $L=\dfrac{20}{-3+3j}=-3.33-3.33j$. As $\omega\to\infty$: $L\to0$ from the third quadrant. The full locus ($\omega:-\infty\to\infty$) is a circle centred at $(-3.33,0)$ of radius $3.33$, passing through $-6.67$ and the origin.
  2. (a) Encirclement and stability. The critical point $-1$ lies inside that circle (between $-6.67$ and $0$), so the locus encircles $-1$ once counter-clockwise: $N=-1$. With $P=1$ open-loop RHP pole, $$Z=P+N=1+(-1)=0,$$ so there are no closed-loop RHP poles — the loop is stable at $k_c=1$. (The single RHP pole demands exactly one CCW encirclement, which the locus supplies.)
  3. Direct check. Characteristic equation $1+L=0\Rightarrow s-3+20k_c=0\Rightarrow s=3-20k_c$. For $k_c=1$, $s=-17<0$ — stable, confirming $Z=0$. ✓
  4. (b) Range of $k_c$. The real-axis crossing at $\omega=0$ is $-20k_c/3$. The locus encircles $-1$ (giving $Z=0$) only when this crossing lies to the left of $-1$: $\dfrac{20k_c}{3}>1$. Equivalently the closed-loop pole $s=3-20k_c$ is in the LHP when $$\boxed{k_c>\frac{3}{20}=0.15}.$$ For $0can stabilise this plant — unlike a process with two mismatched RHP contributions — because there is only one pole to drive into the LHP.
ReIm-6-4-2ω=0 (−6.67)ω=3 (−3.33,−3.33)(−1, 0)P=1 RHP pole; locus encircles −1 once CCW ⇒ N=−1, Z=P+N=0 (stable at k₊=1)
Problem 7(a): Nyquist plot of $L=20/(s-3)$ at $k_c=1$ — a circle from $-6.67$ ($\omega=0$) through $(-3.33,-3.33)$ to the origin. It encircles $-1$ once counter-clockwise ($N=-1$); with $P=1$, $Z=P+N=0$, so the closed loop is stable.
ItemResult
Open-loop RHP poles$P=1$ (at $s=+3$)
$L(0)$ at $k_c=1$$-6.67$ (encircles $-1$)
Encirclements / $Z$ at $k_c=1$$N=-1$, $Z=P+N=0$ → stable
Stabilising range$k_c>3/20=0.15$