Scribd
Upload
3 views4 pages

Experiment Matlab Model Validation

The document outlines an experiment aimed at determining the Residence Time Distribution (RTD) using the Axial dispersion model and validating it with MATLAB against experimental data from a PFR with pulse input tracer. It describes the theory behind the axial dispersion model, the problem statement, and a detailed MATLAB procedure for calculating and validating the RTD. Results include MATLAB code, plots, mean squared error, and correlation coefficients for model performance assessment.

Uploaded by

Apoorv Upadhyay
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
3 views4 pages

Experiment Matlab Model Validation

The document outlines an experiment aimed at determining the Residence Time Distribution (RTD) using the Axial dispersion model and validating it with MATLAB against experimental data from a PFR with pulse input tracer. It describes the theory behind the axial dispersion model, the problem statement, and a detailed MATLAB procedure for calculating and validating the RTD. Results include MATLAB code, plots, mean squared error, and correlation coefficients for model performance assessment.

Uploaded by

Apoorv Upadhyay
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

Experiment

Title: Determination of Residence Time Distribution using Axial dispersion model and its
validation.

Aim:

 Determination of Residence Time Distribution using Axial dispersion model using


MATLAB
 MATLAB Validation of the model’s output with experimental data obtained from a PFR
using pulse input tracer

Theory:

Models are useful for representing flow in real vessels, for scale up, and for diagnosing poor flow.
We have different kinds of models depending on whether flow is close to plug, mixed, or
somewhere in between. primarily with small deviations from plug flow. There are two models for
this: the dispersion model and the tanks-in-series model. They are roughly equivalent. These
models apply to turbulent flow in pipes, laminar flow in very long tubes, flow in packed beds,
shaft kilns, long channels, screw conveyers, etc.

The axial dispersion model accounts for deviations from ideal plug flow in tubular reactors due to
axial mixing. The residence time distribution (RTD), represented by the E-curve (E(θ) or E(t)),
describes the time spent by fluid elements in the reactor. The dimensionless form of the governing
equation is:

𝝏𝑪 𝑫 𝝏𝟐 𝑪 𝝏𝑪
= ∗ −
𝝏𝜽 𝒖𝑳 𝝏𝒛𝟐 𝝏𝒛

Where:

 C = c/c0 Dimensionless concentration


 θ = t/τ: Dimensionless time (τ = 𝒕̅)
 z = x/L: Dimensionless axial position
 D/uL: Dispersion number (indicates degree of mixing)
 Pe = uL/D: Peclet number (inverse of dispersion number)
For a pulse input in an open system and having dispersion number > 0.01, the solution of
above differential equation is given by following set of equations:

Problem Statement:

The tubular reactor has:

 Volume (V) = __ L
 Feed flow rate (q) = __ LPH
 Dispersion Number (D/uL) = __ (should be > 0.01)

Determine the RTD of the reactor through Axial dispersion model. Validate the findings with
obtained experimental data of the reactor using pulse tracer input. Give the parameters
related to validation.

Procedure/MATLAB Approach:

1. Start
a. Initialize the program
2. Define reactor parameters
a. Ask for volume, feed flow rate and dispersion number of reactor.
b. Calculate mean residence time τ.
3. Generate time domain
a. Define time range as per the experimental time domain.
b. Calculate θ.
4. Calculate E-Curve for D/uL > 0.01
a. For each value of θ, compute E(θ) as per the solution.
b. Store the E(θ) values as Et.
5. Calculate Variance
a. σ(θ) as per the formula.
6. Display Results
a. Display τ.
b. Variance σ(θ)
7. Upload experimental data
a. Upload excel file having experimental values of θ and E(θ).
b. Store these values as θ_e and E(θ)_e.
8. Plot E(θ) vs θ curves
a. Plot theoretical values.
b. Plot experimental values on same graph.
c. Label both axis.
d. Name the curve legends.
e. Name the plot title.
9. Validation parameters
a. Calculate mean squared error (MSE) and display the value (Lower MSE indicates
better model performance).
b. Calculate correlation coefficient and display the result.
c. Conclude correlation as per the following criteria.

10. End
a. Terminate the program.
Results:

 Attach developed MATLAB Code.


 Attach generated plot.
 Attach command window image.
 Mean squared error (MSE): __
 Correlation coefficient: __
 Correlation interpretation: __

Conclusion: ___

Precautions: ___

References:

Levenspiel, O. (2006). Chemical Reaction Engineering, 4th Edition. Wiley India. ISBN: 978-
0471530169

You might also like