13.07.2015 Views

Assessment and Future Directions of Nonlinear Model Predictive ...

Assessment and Future Directions of Nonlinear Model Predictive ...

Assessment and Future Directions of Nonlinear Model Predictive ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Close-Loop Stochastic Dynamic Optimization 309Fig. 2. Back-<strong>of</strong>f strategy within moving horizon; back-<strong>of</strong>f from active constraintsthe set-point trajectory would be infeasible. By introducing a back-<strong>of</strong>f from theactive constraints in the optimization, the region <strong>of</strong> the set-point trajectory ismoved inside the feasible region <strong>of</strong> the process to ensure, on the one h<strong>and</strong>, feasibleoperation, <strong>and</strong> to operate the process, on the other h<strong>and</strong>, still as closelyto the true optimum as possible. By this means, the black-marked area in Fig.2 illustrates the corrected bounds ỹ max <strong>of</strong> the hard constraints. Here, it shouldhowever be noted that due to the severe bound at the computation <strong>of</strong> the previoushorizon, the initial value at t 0 is rather far away from the constraint limitin the feasible area. Thus, in the first intervall <strong>of</strong> the current moving horizon,the bound is set at the original physical limit to avoid infeasibility. The back-<strong>of</strong>fadjustment starts from the second interval, i.e. from the time point on, where thenext re-optimization begins. The size <strong>of</strong> ỹ max strongly depends on parametricuncertainty, disturbances, <strong>and</strong> the deviation by measurement errors. Thus, theconstraints in (8) within the moving horizon (8 intervals) are now reformulatedas follows with j =2, ..., 8, α =0.5, ˜T max =4K <strong>and</strong> ˜T ad, max =3K:T (j) ≤ 356 K − ˜T max · α (j−2) ; T ad (j) ≤ 500 K − ˜T ad, max · α (j−2) (7)The decision variable is the cooling flow rate. In order to test robustnesscharacteristics <strong>of</strong> the controler, the performances <strong>of</strong> the open-loop nominal solution,the nominal NMPC, <strong>and</strong> the NMPC with the proposed adaptive back-<strong>of</strong>fapproach are compared under different disturbances, namely: catalyst activitymismatch <strong>and</strong> fluctuations <strong>of</strong> the reactor jacket cooling fluid temperature. Additionally,all measurements are corrupted with white noise e.g. component amount8% <strong>and</strong> temperature 2%.3.1 Dynamik Real–Time OptimizationThe size <strong>of</strong> the dynamic operating region around the optimum (see Fig. 2 right)is affected by fast disturbances. These are, however, efficiently buffered by theproposed regulatory NMPC-based approach. On the other h<strong>and</strong>, there are, infact, slowly time-varying non-zero mean disturbances or drifting model parameterswhich change the plant optimum with time. Thus, a online re-optimizationi.e. dynamic real-time optimization (D-RTO) may be indispensable for an optimaloperation. When on-line measurement gives access to the system state,

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!