29.04.2014 Views

Diploma thesis as pdf file - Johannes Kepler University, Linz

Diploma thesis as pdf file - Johannes Kepler University, Linz

Diploma thesis as pdf file - Johannes Kepler University, Linz

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

Chapter 2<br />

Interior point method for the<br />

numerical solution of the obstacle<br />

problem<br />

2.1 The central path<br />

The idea of the interior point methods is to replace the constrained minimization<br />

problem by a sequence of unconstrained minimization problems, for solving of<br />

which we can use Newton’s method. An objective functional of the unconstrained<br />

problem we generate by adding barrier functional to the objective functional of the<br />

original constrained problem. Barrier function serves <strong>as</strong> barrier against leaving of<br />

the elements the fe<strong>as</strong>ible region K. Each of the problems in the sequence corresponds<br />

to the objective functional J κ (v) depending on the nonnegative penalty<br />

parameter κ. For our problem we construct J κ (v) in the following way:<br />

∫<br />

J κ (v) = J(v) − κ ln(v − ψ)dω.<br />

Ω<br />

We extend the definition of the ln to the whole real domain axis by setting lnz =<br />

−∞,<br />

z ≤ 0. Then barrier function approaches infinity <strong>as</strong> the elements from the<br />

interior approach the boundary. Thus we obtained the family of the following<br />

unconstrained minimization problem:<br />

11

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

Saved successfully!

Ooh no, something went wrong!