Example Problem for Experimental Design: Difference between revisions
From charlesreid1
| Line 3: | Line 3: | ||
Problem description: http://files.charlesmartinreid.com/VUQ_Toy_Problem.pdf | Problem description: http://files.charlesmartinreid.com/VUQ_Toy_Problem.pdf | ||
=I/U Map= | =Inputs and Outputs= | ||
There are several (polynomial) response surfaces being fit for the Monte Carlo simulations. This is because there is one response surface for each output or observable. (Technically, these are all part of one large multivariate response surface, but it is easier to think about them as independent response surfaces). | |||
The variables that are included in the response surface analysis are: | |||
* <math>z_i</math> (for <math>i=1,2,3</math>) - the location of measurement of axial concentrations | |||
* <math>\dot{m} = \dot{m}_1 = \dot{m}_2</math> - the mass flowrate of the inlet streams of A and B | |||
* <math>L_{mix}</math> - mixing length (parameter for the mixing model) | |||
* <math>k(T)</math> - reaction rate for the reaction <math>A + B \rightarrow^{k} P</math> | |||
==Variables: I/U Map== | |||
{|class="wikitable" | {|class="wikitable" | ||
| Line 11: | Line 22: | ||
!Logarithmic Scale? | !Logarithmic Scale? | ||
|- | |- | ||
|<math> | |<math>z_1</math> | ||
|<math>0.5 m</math> | |<math>0.5 m</math> | ||
|<math>\pm 0.02 m</math> | |<math>\pm 0.02 m</math> | ||
| Line 17: | Line 28: | ||
|- | |- | ||
|<math> | |<math>z_2</math> | ||
|<math>1.5 m</math> | |<math>1.5 m</math> | ||
|<math>\pm 0.02 m</math> | |<math>\pm 0.02 m</math> | ||
| Line 23: | Line 34: | ||
|- | |- | ||
|<math> | |<math>z_3</math> | ||
|<math>2.5 m</math> | |<math>2.5 m</math> | ||
|<math>\pm 0.02 m</math> | |<math>\pm 0.02 m</math> | ||
| Line 65: | Line 76: | ||
|} | |} | ||
==Response Surfaces== | |||
===Product at Exit Response Surface=== | |||
This response surface maps the response of the mass fraction of product at the exit, <math>y_{P,exit}</math>. | |||
This response surface is a function of several variables: | |||
<math> | |||
y_{P,exit} = y_{P,exit} \left( \dot{m}, k, L_{mix} \right) | |||
</math> | |||
However, the axial location of measurement of product <math>z_i</math> is not considered, because it does not affect the measurement of P at the exit. | |||
===Product at Axial Location Response Surfaces=== | |||
This response surface maps the response of the mass fraction of product at several axial locations, <math>y_{P,z1}, y_{P,z2}, y_{P,z3}</math>. | |||
This response surface is a function of all variables: | |||
<math> | |||
y_{P,zi} = y_{P,exit} \left( \dot{m}, k, L_{mix}, z_i \right) | |||
</math> | |||
(Note, however, that only the corresponding <math>z_i</math> will be a response surface independent variable, since other <math>z_i</math> values have no affect). | |||
===Dealing with Multimodal Variables=== | |||
The variable <math>\dot{m}</math> is an excellent example of modal variables - that is, variables that have certain modes, or distinct sets of values. | |||
In the case of <math>\dot{m}</math>, the two modes are <math>\dot{m} = 1.0 \pm 5%</math> and <math>\dot{m} = 2.0 \pm 5%</math>. | |||
Revision as of 04:53, 29 June 2011
Problem Description
Problem description: http://files.charlesmartinreid.com/VUQ_Toy_Problem.pdf
Inputs and Outputs
There are several (polynomial) response surfaces being fit for the Monte Carlo simulations. This is because there is one response surface for each output or observable. (Technically, these are all part of one large multivariate response surface, but it is easier to think about them as independent response surfaces).
The variables that are included in the response surface analysis are:
- $ z_i $ (for $ i=1,2,3 $) - the location of measurement of axial concentrations
- $ \dot{m} = \dot{m}_1 = \dot{m}_2 $ - the mass flowrate of the inlet streams of A and B
- $ L_{mix} $ - mixing length (parameter for the mixing model)
- $ k(T) $ - reaction rate for the reaction $ A + B \rightarrow^{k} P $
Variables: I/U Map
| Variable Name | Input value (I) | Uncertainty (U) | Logarithmic Scale? |
|---|---|---|---|
| $ z_1 $ | $ 0.5 m $ | $ \pm 0.02 m $ | no |
| $ z_2 $ | $ 1.5 m $ | $ \pm 0.02 m $ | no |
| $ z_3 $ | $ 2.5 m $ | $ \pm 0.02 m $ | no |
| $ \dot{m}_1 $ | $ 1.0 $ | $ \pm 0.05 $ | no |
| $ 2.0 $ | $ \pm 0.10 $ | ||
| $ \dot{m}_2 $ | $ 1.0 $ | $ \pm 0.05 $ | no |
| $ 2.0 $ | $ \pm 0.10 $ | ||
| $ k(T) $ | $ 1 $ | $ 10^{0 \pm 2} $ | yes |
| $ L_{mix} $ | $ 0.3-3.0 $ | yes |
Response Surfaces
Product at Exit Response Surface
This response surface maps the response of the mass fraction of product at the exit, $ y_{P,exit} $.
This response surface is a function of several variables:
$ y_{P,exit} = y_{P,exit} \left( \dot{m}, k, L_{mix} \right) $
However, the axial location of measurement of product $ z_i $ is not considered, because it does not affect the measurement of P at the exit.
Product at Axial Location Response Surfaces
This response surface maps the response of the mass fraction of product at several axial locations, $ y_{P,z1}, y_{P,z2}, y_{P,z3} $.
This response surface is a function of all variables:
$ y_{P,zi} = y_{P,exit} \left( \dot{m}, k, L_{mix}, z_i \right) $
(Note, however, that only the corresponding $ z_i $ will be a response surface independent variable, since other $ z_i $ values have no affect).
Dealing with Multimodal Variables
The variable $ \dot{m} $ is an excellent example of modal variables - that is, variables that have certain modes, or distinct sets of values.
In the case of $ \dot{m} $, the two modes are $ \dot{m} = 1.0 \pm 5% $ and $ \dot{m} = 2.0 \pm 5% $.