ChEn-3170: Computational Methods in Chemical Engineering Fall 2018 UMass Lowell; Prof. V. F. de Almeida 08Nov2018
$ \newcommand{\Amtrx}{\boldsymbol{\mathsf{A}}} \newcommand{\Bmtrx}{\boldsymbol{\mathsf{B}}} \newcommand{\Cmtrx}{\boldsymbol{\mathsf{C}}} \newcommand{\Mmtrx}{\boldsymbol{\mathsf{M}}} \newcommand{\Imtrx}{\boldsymbol{\mathsf{I}}} \newcommand{\Pmtrx}{\boldsymbol{\mathsf{P}}} \newcommand{\Qmtrx}{\boldsymbol{\mathsf{Q}}} \newcommand{\Lmtrx}{\boldsymbol{\mathsf{L}}} \newcommand{\Umtrx}{\boldsymbol{\mathsf{U}}} \newcommand{\xvec}{\boldsymbol{\mathsf{x}}} \newcommand{\yvec}{\boldsymbol{\mathsf{y}}} \newcommand{\zvec}{\boldsymbol{\mathsf{z}}} \newcommand{\avec}{\boldsymbol{\mathsf{a}}} \newcommand{\bvec}{\boldsymbol{\mathsf{b}}} \newcommand{\cvec}{\boldsymbol{\mathsf{c}}} \newcommand{\rvec}{\boldsymbol{\mathsf{r}}} \newcommand{\norm}[1]{\bigl\lVert{#1}\bigr\rVert} \DeclareMathOperator{\rank}{rank} \DeclareMathOperator{\abs}{abs} $
your name
¶Context | Points |
---|---|
Precision of the answer | 80% |
Answer Markdown readability | 10% |
Code readability | 10% |
Import the reaction rate constant data from the course repository data/modified-k_x_T_10pts.dat
and make a plot of $k$ versus $T$.
'''1.1 Read and plot experimental data'''
'''1.2 Build the matrix system'''
rank(A) = 3
'''1.3 Compute LS values of the parameters in the modified Arrhenius expression'''
x_vec = [ 3.89438185e+00 2.31978194e+03 -1.47323394e-01] k_0 = 4.91257e+01 [1/s] E_a = 2.31978e+03 [J/mol] alpha = -1.47323e-01 []
'''1.3 Plot of the transformed plane data and fit'''
Is the fit a straight line?
'''1.3 Plot the data and fit in the Cartesian plane'''
'''1.4 Goodness of fit'''
phi(x_LS) = 2.105e-03 sqrt(phi(x_LS)) = 4.589e-02 goodness of fit = 1.451e-02
'''1.4 Contour plot of objective function'''
Import the reaction rate constant data from the course repository data/modified-k_x_T_30pts.dat
and make a plot of $k$ versus $T$.
'''2.1 Read and plot experimental data'''
'''2.2 Build the matrix system'''
rank(A) = 3
'''2.3 Compute LS values of the parameters in the modified Arrhenius expression'''
x_vec = [7.61614166e-01 4.54575077e+02 2.82450412e-01] k_0 = 2.14173e+00 [1/s] E_a = 4.54575e+02 [J/mol] alpha = 2.82450e-01 []
'''2.3 Plot of the transformed plane data and fit'''
Is the fit a straight line?
'''2.3 Plot the data and fit in the Cartesian plane'''
'''2.4 Goodness of fit'''
phi(x_LS) = 9.042e-03 sqrt(phi(x_LS)) = 9.509e-02 goodness of fit = 1.736e-02
'''2.4 Contour plot of objective function'''
Import the reaction rate constant data from the course repository data/modified-k_x_T_90pts.dat
and make a plot of $k$ versus $T$.
'''3.1 Read and plot experimental data'''
'''3.2 Build the matrix system'''
rank(A) = 3
'''3.3 Compute LS values of the parameters in the modified Arrhenius expression'''
x_vec = [ 0.32736356 269.68912137 0.34493353] k_0 = 1.38731e+00 [1/s] E_a = 2.69689e+02 [J/mol] alpha = 3.44934e-01 []
'''3.3 Plot of the transformed plane data and fit'''
Is the fit a straight line?
'''2.3 Plot the data and fit in the Cartesian plane'''
'''2.4 Goodness of fit'''
phi(x_LS) = 3.046e-02 sqrt(phi(x_LS)) = 1.745e-01 goodness of fit = 1.840e-02
'''2.4 Contour plot of objective function'''
Make a table and compare the values of the modified Arrhenius parameters for each Assignment 1 to 3. Compare the values of the norm of the residual and goodness of fit for each data set. Which data set fits the model best? Explain.
2. Would you say the experimental data is consistent? Explain. 3. Based on the values of the parameters computed for the experimental data, what advice would you give to a company designing a chemical reaction system based on this experimental data?Best fit?
Data consistency?
Design advice?