02.11.2012 Views

Optimisation of Marine Boilers using Model-based Multivariable ...

Optimisation of Marine Boilers using Model-based Multivariable ...

Optimisation of Marine Boilers using Model-based Multivariable ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

36 3. SUMMARY OF CONTRIBUTIONS<br />

The model was <strong>based</strong> on first principles <strong>using</strong> mass and energy balances for describing<br />

the temperature and pressure dynamics. As flow through and efficiency <strong>of</strong> the turbocharger<br />

was required in these balance equations, approximations <strong>of</strong> the turbocharger<br />

turbine and compressor maps partly <strong>based</strong> on physical insight were made. The oxygen<br />

was modelled by describing mole balances for the oxygen in the gas generator and furnace<br />

volumes. The resulting model was <strong>of</strong> sixth order and presented in descriptor form<br />

as:<br />

F(x) dx<br />

= h(x,u,d)<br />

dt<br />

(3.7a)<br />

y = g(x,u,d) (3.7b)<br />

where u = [ ˙mfu1, ˙mfu2] T is the vector <strong>of</strong> inputs being the fuel mass flow to the gas<br />

generator and the furnace. d = [Ta,Tfu,Tm] T is the vector <strong>of</strong> disturbances being<br />

temperatures <strong>of</strong> the inlet air, <strong>of</strong> the fuel and <strong>of</strong> the metal separating the hot flue gas<br />

and the water/steam in the boiler. x = [pgg,Tgg,ω,Tfn,xgg,O2 ,xfn,O2 ]T is the state<br />

vector being pressure in the gas generator, temperature in the gas generator, shaft speed,<br />

temperature in the furnace, and oxygen fraction in the gas generator and the furnace.<br />

Finally, y = [ ˙mfu, ˙ Q,xfn,O2 ]T is the output vector with ˙mfu = ˙mfu1 + ˙mfu2 and<br />

˙Q is the energy transferred to the metal walls. Expanding, (3.7a) reveals the model<br />

structure as:<br />

⎡<br />

⎢<br />

⎣<br />

f11 f12 0 0 0 0<br />

f21 f22 0 0 0 0<br />

0 0 f33 0 0 0<br />

0 0 0 f44 0 0<br />

0 0 0 0 f55 0<br />

0 0 0 0 0 f66<br />

⎤⎡<br />

⎥⎢<br />

⎥⎢<br />

⎥⎢<br />

⎥⎢<br />

⎥⎢<br />

⎥⎢<br />

⎦⎢<br />

⎣<br />

dpgg<br />

dt<br />

dTgg<br />

dt<br />

dω<br />

dt<br />

dTfn<br />

dt<br />

dxgg,O2 dt<br />

dxfn,O2 dt<br />

⎤<br />

⎡<br />

⎥ ⎢<br />

⎥ ⎢<br />

⎥ ⎢<br />

⎥ = ⎢<br />

⎥ ⎢<br />

⎥ ⎢<br />

⎦ ⎣<br />

h1<br />

h2<br />

h3<br />

h4<br />

h5<br />

h6<br />

⎤<br />

⎥<br />

⎦<br />

(3.8)<br />

where the elements fij and hi were found in the model derivation – see [Solberg et al.,<br />

2008c]. F is never singular, hence has a well defined inverse, so (3.7a) can be written<br />

as an ordinary differential equation: ˙x = f(x,u,d) = F −1 (x)h(x,u,d).<br />

Validation against measurements collected from the test setup shows good agreement<br />

between model and measuring data in terms <strong>of</strong> capturing the dynamical behaviour.<br />

However, in terms <strong>of</strong> stationary values these are not represented well by the model<br />

for other outputs than the oxygen level, see Figure 3.5. This was accepted for now<br />

as the main focus was on oxygen control, and constraints on internal states were not<br />

considered.<br />

The performance requirements for the controller <strong>of</strong> the burner were to deliver the requested<br />

fuel flow while keeping the oxygen percentage above 3% and maximise efficiency.<br />

No requirement regarding <strong>of</strong>f-set free tracking <strong>of</strong> the fuel flow reference was<br />

introduced meaning that differences between requested and actual fuel flow should be<br />

handled by e.g. including integral action in the outer controller.

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

Saved successfully!

Ooh no, something went wrong!