Skip to main content
Engineering LibreTexts

1.8: Industrial Process Models

  • Page ID
    17748
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    Industrial processes comprise procedures involving exchange of chemical, electrical or mechanical energy to aid in the manufacturing of industrial products. Industrial process models are mathematical models used to describe those processes.

    An industrial process model, in its simplified form, can be represented by a first-order ODE accompanied by a dead-time, i.e., there is a finite time delay between the application of input and the appearance of the process output.

    Let \(\tau\) represent the time constant associated with an industrial process, \(\tau _

    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[3]/span[2]/span, line 1, column 1
    
    \) represent the dead time, and \(K\) represent the process dc gain; then, the simplified industrial process dynamics are represented by the following delay-differential equation: \[\tau \frac{{\rm d}y(t)}{{\rm d}t} +y(t)=Ku(t-t_{d} ).\] Application of the Laplace transform produces the following first-order-plus-dead-time (FOPDT) model of an industrial process: \[G(s)=\frac{Ke^{-\tau _
    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[3]/span[5]/span, line 1, column 1
    
    s} }{\tau s+1} .\]
    The process parameters \(\{ K,\; \tau ,\; \tau _
    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[3]/span[6]/span, line 1, column 1
    
    \}\)
    , respectively, denote the process gain, the process time constant, and the process dead time. These parameters can be identified from the process response to a unit-step input.

    We note that the process model involves a transendental function, \(e^{-\tau _

    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[4]/span[1]/span, line 1, column 1
    
    s}\). For analysis and controller design purposes, this term is replaced by a suitable rational approximation. Typical approximations obtained from Taylor series expansion of the transendental term are given as: \[e^{-\tau _
    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[4]/span[2]/span[1], line 1, column 1
    
    s} \approx 1-\tau _
    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[4]/span[2]/span[2], line 1, column 1
    
    s, e^{-\tau _
    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[4]/span[2]/span[3], line 1, column 1
    
    s} =\frac{1}{e^{\tau _
    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[4]/span[2]/span[4], line 1, column 1
    
    s} } \approx \frac{1}{1+\tau _
    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[4]/span[2]/span[5], line 1, column 1
    
    s} , e^{-{\tau }_ds}=\frac{e^{-{\tau }_ds/2}}{e^
    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[4]/span[2]/span[6], line 1, column 1
    
    \approx \frac{1-{\tau }_ds/2}{1+{\tau }_ds/2}.\]
    The last expression above is called the first-order Pade’ approximation and is popular in idustrial process control applications.

    Example 1.13: An industrial process model

    The process parameters of a stirred-tank bioreactor are given as: \(\{ K,\; \tau ,\; \tau _

    ParseError: invalid DekiScript (click for details)
    Callstack:
        at (Courses/University_of_Arkansas_Little_Rock/Introduction_to_Control_Systems_(Iqbal)/01:_Mathematical_Models_of_Physical_Systems/1.08:_Industrial_Process_Models), /content/body/p[6]/span[1]/span, line 1, column 1
    
    \} =\left\{20,0.5,1\right\}\). The transfer function model of the process is formed as: \(G(s)=\frac{20e^{-s} }{0.5s+1}\).

    By using a first-order Pade’ approximation, a rational transfer function model of the industrial process with delay is obtained as: \(G(s)=\frac{20\left(1-0.5s\right)}{\left(0.5s+1\right)^{2} }\).

    Further, a model of this kind is called non-minimum phase due to the extra phase added by the right-half-plane zero in the numerator (see Chapter 2 for details).


    1.8: Industrial Process Models is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.

    • Was this article helpful?