s2 - process product optimization using design experiments and response surface methodolgy
DESCRIPTION
An intensive practical course mainly for PhD-students on the use of designs of experiments (DOE) and response surface methodology (RSM) for optimization problems. The course covers relevant background, nomenclature and general theory of DOE and RSM modelling for factorial and optimisation designs in addition to practical exercises in Matlab. Due to time limitations, the course concentrates on linear and quadratic models on the k≤3 design dimension. This course is an ideal starting point for every experimental engineering wanting to work effectively, extract maximal information and predict the future behaviour of their system. Mikko Mäkelä (DSc, Tech) is a postdoctoral fellow at the Swedish University of Agricultural Sciences in Umeå, Sweden and is currently visiting the Department of Chemical Engineering at the University of Alicante. He is working in close cooperation with Paul Geladi, Professor of Chemometrics, and using DOE and RSM for process optimization mainly for the valorization of industrial wastes in laboratory and pilot scales.”TRANSCRIPT
![Page 1: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/1.jpg)
Process/product optimization using design of experiments and response surface methodology
M. Mäkelä
Sveriges landbruksuniversitetSwedish University of Agricultural Sciences
Department of Forest Biomaterials and TechnologyDivision of Biomass Technology and ChemistryUmeå, Sweden
![Page 2: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/2.jpg)
Contents
Practical course, arranged in 4 individual sessions:
Session 1 – Introduction, factorial design, first order models
Session 2 – Matlab exercise: factorial design
Session 3 – Central composite designs, second order models, ANOVA,
blocking, qualitative factors
Session 4 – Matlab exercise: practical optimization example on given data
![Page 3: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/3.jpg)
Session 1Introduction
Why experimental design
Factorial design
Design matrix
Model equation = coefficients
Residual
Response contour
![Page 4: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/4.jpg)
Session 2
Factorial design
Research problem
Design matrix
Model equation = coefficients
Degrees of freedom
Predicted response
Residual
ANOVA
R2
Response contour
![Page 5: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/5.jpg)
Research problem
A chemist is interested on the effect of temperature (A), catalyst concentration (B)
and time (C) on the molecular weight of polymer produced
She performed a 23 factorial design
Molecular weight (in order): 2400, 2410, 2315, 2510, 2615, 2625, 2400, 2750
Parameter Low High
Temperature, A (°C) 100 120
Catalyst conc., B (%) 4 8
Time, C (min) 20 30
Myers, Montgomery & Anderson-Cook, Response Surface Methodology, 3rd ed., 2009, 131.
![Page 6: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/6.jpg)
Research problem
Empirical model:
, ,
⋯
In matrix notation:
→
yy⋮y
1 ⋯1 ⋯1 ⋮ ⋮ ⋱ ⋮1 ⋯
bb⋮b
ee⋮e
Measure ChooseUnknown!
![Page 7: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/7.jpg)
Design matrixN:o A B C Resp.
1 -1 -1 -1 2400
2 1 -1 -1 2410
3 -1 1 -1 2315
4 1 1 -1 2510
5 -1 -1 1 2615
6 1 -1 1 2625
7 -1 1 1 2400
8 1 1 1 2750
Build a design matrix with interactions in Matlab
![Page 8: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/8.jpg)
Coefficients
Model coefficients
Least squares calculation
Coefficient significance?
8 model terms
![Page 9: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/9.jpg)
Degrees of freedom
In statistics, degrees of freedom (df) are lost by imposing linear constraints on
a sample
E.g. sample variance:
∑
where ∑ 0. Hence 1 residuals can be used to completely
determine the other.
In RSM, dfs are lost due to the constraints imposed by the coefficients
Residual df with observations and 1 model terms
![Page 10: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/10.jpg)
Coefficients
Remove unimportant model terms
and calculate new coefficients
![Page 11: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/11.jpg)
Predicted response
Calculate predicted response,
Based on X and b
Graphical tools are useful
![Page 12: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/12.jpg)
Residuals
Calculate model residuals
A common way is to scale the residuals
E.g. standardized residual
→ Need an error approximation
![Page 13: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/13.jpg)
Error estimation
Estimated error of predicted response
∑ , where p = k + 1
Standardized residuals
>│3│ a potential outlier
![Page 14: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/14.jpg)
ANOVA
Analysis of variance (ANOVA) allows to statistically test the model
Fisher’s F test
H0: ⋯ 0
H1: 0 for at least one j
Parameter df Sum of squares (SS)
Meansquare (MS)
F-value p-value
Total corrected n-1 SStot MStot
Regression k SSmod MSmod MSmod/MSres
<0.05>0.05
Residual n-k-1 SSres MSres
Lack of fit
Pure error
![Page 15: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/15.jpg)
R2 statistic
Data variability explained by the
model
Compares model and total sum of
squares
R2 = 95.2%
![Page 16: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/16.jpg)
Response contour
For a k=3 design, 2D contours require
one constant factor
![Page 17: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/17.jpg)
Response contour
Use of the model
Optimization within specific coordinates
Prediction of an optimum
Verification
A = 115 and B = 7
1.5.51.25
2645 90
![Page 18: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/18.jpg)
Research problem
The chemist has a possibility to perform three
additional center-points in her factorial design
Molecular weight values: 2800, 2750, 2810
Can this change our view on the behaviour of
the system?
A
BC
![Page 19: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/19.jpg)
Session 2
Factorial design
Research problem
Design matrix
Model equation = coefficients
Degrees of freedom
Predicted response
Residual
ANOVA
R2
Response contour
![Page 20: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/20.jpg)
Nomenclature
Design matrix
Coefficient
Degree of freedom
Prediction
Residual
Outlier
ANOVA
Sum of squares
Mean square
Contour
![Page 21: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/21.jpg)
Contents
Practical course, arranged in 4 individual sessions:
Session 1 – Introduction, factorial design, first order models
Session 2 – Matlab exercise: factorial design
Session 3 – Central composite designs, second order models, ANOVA,
blocking, qualitative factors
Session 4 – Matlab exercise: practical optimization example on given data
![Page 22: S2 - Process product optimization using design experiments and response surface methodolgy](https://reader031.vdocument.in/reader031/viewer/2022020110/559c19db1a28ab22598b46f1/html5/thumbnails/22.jpg)
Thank you for listening!