Mathematical Model
Mathematical Model
Mathematical Model
1. State the modeling objectives and the end use of the model. They determine the
required levels of model detail and model accuracy.
2. Draw a schematic diagram of the process and label all process variables.
3. List all of the assumptions that are involved in developing the model. Try for
parsimony; the model should be no more complicated than necessary to meet the
modeling objectives.
4. Determine whether spatial variations of process variables are important. If so, a
partial differential equation model will be required.
5. Write appropriate conservation equations (mass, component, energy, and so
forth).
6. Introduce equilibrium relations and other algebraic equations (from
thermodynamics, transport phenomena, chemical kinetics, equipment geometry,
etc.).
7. Perform a degrees of freedom analysis (Section 2.3) to ensure that the model
equations can be solved.
8. Simplify the model. It is often possible to arrange the equations so that the
dependent variables (outputs) appear on the left side and the independent
variables (inputs) appear on the right side. This model form is convenient for
computer simulation and subsequent analysis.
9. Classify inputs as disturbance variables or as manipulated variables.
Modeling Approaches
Physical/chemical (fundamental, global)
Model structure by theoretical analysis
Material/energy balances
Heat, mass, and momentum transfer
Thermodynamics, chemical kinetics
Physical property relationships
Model complexity must be determined (assumptions)
Can be computationally expensive (not real-time)
May be expensive/time-consuming to obtain
Good for extrapolation, scale-up
Does not require experimental data to obtain (data required for
validation and fitting)
Conservation Laws
Conservation of Mass
rate of mass
rate of mass rate of mass
(2-6)
accumulation in out
Conservation of Component i
rate of component i i
rate of component
accumulation in
Conservation of Energy
The general law of energy conservation is also called the First Law of
Thermodynamics. It can be expressed as:
Number of
births (West
Germany)
or d V
w1 w2 w (2-2)
dt
d Vx
w1 x1 w2 x2 wx (2-3)
dt
The corresponding steady-state model was derived in Ch. 1 (cf. Eqs. 1-1 and 1-2).
0 w1 w2 w (2-4)
0 w1 x1 w2 x2 wx (2-5)
Equation 2-13 can be simplified by expanding the accumulation term using the chain
rule for differentiation of a product:
d Vx dx dV
V x (2-14)
dt dt dt
Substitution of (2-14) into (2-13) gives:
dx dV
V x w1 x1 w2 x2 wx (2-15)
dt dt
Substitution of the mass balance in (2-12) for dV / dt in (2-15) gives:
dx
V x w1 w2 w w1x1 w2 x2 wx (2-16)
dt
After canceling common terms and rearranging (2-12) and (2-16), a more convenient
model form is obtained:
dV 1
w1 w2 w (2-17)
dt
dx w1 w
x1 x 2 x2 x (2-18)
dt V V
Stirred-Tank Heating Process
Assumptions:
1. Perfect mixing; thus, the exit temperature T is also the temperature of the tank
contents.
2. The liquid holdup V is constant because the inlet and outlet flow rates are equal.
3. The density and heat capacity C of the liquid are assumed to be constant. Thus,
their temperature dependence is neglected.
4. Heat losses are negligible.
1. Changes in potential energy and kinetic energy can be neglected because they are
small in comparison with changes in internal energy.
2. The net rate of work can be neglected because it is small compared to the rates of
heat transfer and convection.
Model Development - I
For a pure liquid at low or moderate pressures, the internal energy is approximately equal
to the enthalpy, U int H , and H depends only on temperature. Consequently, in the
subsequent development, we assume that Uint = H and U int H where the caret (^) means
per unit mass. As shown in Appendix B, a differential change in temperature, dT,
produces a corresponding change in the internal energy per unit mass, dU int ,
Model Development - II
An expression for the rate of internal energy accumulation can be derived from Eqs. (2-
29) and (2-30):
dU int dT
VC (2-31)
dt dt
Note that this term appears in the general energy balance of Eq. 2-10.
Suppose that the liquid in the tank is at a temperature T and has an enthalpy, H .
Integrating Eq. 2-29 from a reference temperature Tref to T gives,
H H ref C T Tref (2-32)
where H ref is the value of H at Tref. Without loss of generality, we assume that H ref 0
(see Appendix B). Thus, (2-32) can be written as:
H C T Tref (2-33)
H i C Ti Tref (2-34)
Substituting (2-33) and (2-34) into the convection term of (2-10) gives:
wH w C Ti Tref w
C T Tref
(2-35)
w s =0.83 106 g hr
H v =600 cal g
C=l cal g o C
w=104 kg hr
=103 kg m3
V=20 m3
V 2 10 4 kg
V 2 104 kg
2hr
w 104 kg hr
dy
2 =-y+6 10-5 u
dt
Dynamic model
y T T
u ws ws
Step 1: t=0 double ws
T(0)=T y(0)=0
u = +0.83 106 g hr
dy
2 =-y+50
dt
y=50 l-e-0.5t
final T=yss +T=50+90=140o C
o
Step 2: maintain T = 140 C / 24 hr
y=50e-0.5t t
(self-regulating, but slow)
y0
Process Dynamics
Process control is inherently concerned with unsteady state behavior (i.e., "transient
response", "process dynamics")
Stirred tank heater: assume a "lag" between heating element temperature Te, and process
fluid temp, T.
heat transfer limitation = heA(Te T)
Energy balances
dT
Tank: wCTi +h e A(Te -T)-wCT=mC
dt
dTe
Chest: Q-h e A(Te -T)=m e Ce
dt
dT dTe
At s.s. 0, 0
dt dt
Specify Q calc. T, Te
2 first order equations 1 second order equation in T
Relate T to Q (Te is an intermediate variable)
mm e Ce d 2 y me Ce m eC e m dy 1
2
y u
h e A e dt heAe wC w dt wC
Fig. 2.2
1
q= h Rv = line resistance
Rv
dh 1
A qi h (2-57)
dt Rv
linear ODE
P p gh
If q = C v P - Pa Pa : ambient pressure
dh
A qi Cgh
v
*
q Ci hv (2-61)
dt
nonlinear ODE
Table 2.2. Degrees of Freedom Analysis
1. List all quantities in the model that are known constants (or parameters that can be
specified) on the basis of equipment dimensions, known physical properties, etc.
2. Determine the number of equations NE and the number of process variables, NV.
Note that time t is not considered to be a process variable because it is neither a
process input nor a process output.
3. Calculate the number of degrees of freedom, NF = NV NE.
4. Identify the NE output variables that will be obtained by solving the process
model.
5. Identify the NF input variables that must be specified as either disturbance
variables or manipulated variables, in order to utilize the NF degrees of freedom.
3 parameters: V , ,C
4 variables: T , Ti , w, Q
1 equation: Eq. 2-36
Thus the degrees of freedom are NF = 4 1 = 3. The process variables are classified as:
1 output variable: T
3 input variables: Ti, w, Q
For temperature control purposes, it is reasonable to classify the three inputs as:
Bioreactions
Fed-Batch Bioreactor