12.07.2015 Views

Metrics of curves in shape optimization and analysis - Andrea Carlo ...

Metrics of curves in shape optimization and analysis - Andrea Carlo ...

Metrics of curves in shape optimization and analysis - Andrea Carlo ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

2.4.5 Centroid energyWe will now propose another simple example where the above phenomenon isaga<strong>in</strong> evident.Example 2.8 (Centroid-based energy) Let us fix a target po<strong>in</strong>t v ∈ lR n .We recall that avg c (c) is the center <strong>of</strong> mass <strong>of</strong> the curve. The energyE(c) := 1 2 |avg c(c) − v| 2 (2.9)penalizes the distance from the center <strong>of</strong> mass to v. Let <strong>in</strong> the follow<strong>in</strong>g c =avg c (c) for simplicity. The directional derivative <strong>of</strong> E <strong>in</strong> direction h is〈〉DE(c; h) = c − v, D(c)(c; h)where <strong>in</strong> turn (by eqn. (2.7) <strong>and</strong> (2.8))∫D(c)(c; h) = h + (c − c)(D s h · D s c) ds = (2.10)c∫= h − D s c (h · D s c) − (c − c)(h · Dsc) 2 dssuppos<strong>in</strong>g that the curve is planar, thench − D s c (h · D s c) = N(h · N)soDE(c; h) =∫〈c − v, N〉(h · N) − 〈c − v, c − c〉κ(h · N) ds .The H 0 gradient descent flow is then∂c∂t = −∇ H 0E(c) = 〈(v − c), N〉N − κN〈 (c − c), (v − c) 〉 .The first term 〈(v − c) · N〉N <strong>in</strong> this gradient descentflow moves the whole curve towards v.PvLet P := {w : 〈(w − c) · (v − c) 〉 ≥ 0} be the half plane“on the v side” . The second term −κN 〈 (c − c) · (v − c) 〉<strong>in</strong> this gradient descent flow tries to decrease the curvelength out <strong>of</strong> P <strong>and</strong> <strong>in</strong>crease the curve length <strong>in</strong> P , <strong>and</strong>this is ill posed.ccWe will come back to this example <strong>in</strong> Proposition 10.20.17

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

Saved successfully!

Ooh no, something went wrong!