23.01.2014 Views

The Reduced Basis Method Problem Set 1

The Reduced Basis Method Problem Set 1

The Reduced Basis Method Problem Set 1

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

discontinuity interfaces Γ i int ≡ ∂Ω0 ∩ ∂Ω i , i = 1, . . . , 4, where ∂Ω i denotes the boundary of Ω i :<br />

u 0 = u i<br />

−(∇u 0 · ˆn i ) = −k i (∇u i · ˆn i )<br />

}<br />

on Γ i int, i = 1, . . . , 4;<br />

here ˆn i is the outward normal on ∂Ω i . Finally, we introduce a Neumann flux boundary condition<br />

on the fin root<br />

−(∇u 0 · ˆn 0 ) = −1 on Γ root , (2)<br />

which models the heat source; and a Robin boundary condition<br />

−k i (∇u i · ˆn i ) = Bi u i on Γ i ext , i = 0, . . . , 4, (3)<br />

which models the convective heat losses. Here Γ i ext is that part of the boundary of Ωi exposed to<br />

the flowing fluid; note that ∪ 4 i=0Γ i ext = Γ\Γ root.<br />

<strong>The</strong> average temperature at the root, T root (µ), can then be expressed as l O (u(µ)), where<br />

∫<br />

l O (v) = v<br />

Γ root<br />

(recall Γ root is of length unity). Note that l(v) = l O (v) for this problem.<br />

Part 1 - Finite Element Approximation<br />

We saw in class that the reduced basis approximation is based on a “truth” finite element<br />

approximation of the exact (or analytic) problem statement. To begin, we have to show that the<br />

exact problem described above does indeed satisfy the affine parameter dependence and thus fits<br />

into the framework shown in class.<br />

α) Show that u e (µ) ∈ X e ≡ H 1 (Ω) satisfies the weak form<br />

a(u e (µ), v; µ) = l(v), ∀v ∈ X e , (4)<br />

with<br />

a(w, v; µ) =<br />

l(v) =<br />

4∑<br />

∫<br />

k i ∇w · ∇v dA + Bi wv dS,<br />

i=0<br />

∫Ω i Γ\Γ root<br />

∫<br />

v dS.<br />

Γ root<br />

β) Optional: show that u e (µ) is the argument that minimizes<br />

J(w) = 1 2<br />

4∑<br />

∫<br />

k i<br />

i=0<br />

∇w · ∇w dA + Bi<br />

Ω i 2<br />

∫<br />

Γ\Γ root<br />

∫<br />

w 2 dS − w dS (5)<br />

Γ root<br />

over all functions w in X e .<br />

We now consider the linear finite element space<br />

X N = {v ∈ H 1 (Ω)| v| Th ∈ IP 1 (T h ), ∀T h ∈ T h },<br />

2

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

Saved successfully!

Ooh no, something went wrong!