17.08.2013 Views

A particle-in-Burgers model: theory and numerics - Laboratoire de ...

A particle-in-Burgers model: theory and numerics - Laboratoire de ...

A particle-in-Burgers model: theory and numerics - Laboratoire de ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

Mo<strong>de</strong>l <strong>and</strong> motivation Auxiliary steps Results h = 0: coupl<strong>in</strong>g h = 0: <strong>de</strong>f<strong>in</strong>ition, uniqueness h = 0: <strong>numerics</strong>, existence The coupled problem<br />

Case of u frozen: evolv<strong>in</strong>g h = h(·)<br />

Proposition (<strong>mo<strong>de</strong>l</strong>l<strong>in</strong>g/“traces” <strong>in</strong>terpretation of the ODE on h(·) )<br />

For every drag force, the ODE <strong>in</strong> the coupled problem writes<br />

mh ′′ <br />

(t) =<br />

−<br />

<br />

(u−) 2 /2−h ′ (t)u−<br />

<br />

(u+) 2 /2−h ′ (t)u+<br />

Notice that the right-h<strong>and</strong> si<strong>de</strong> above is expressed as the difference of the<br />

normal components of the 2D-field (u, u 2 /2) on the curve {x = h(t)} from<br />

the left <strong>and</strong> from the right . Comb<strong>in</strong><strong>in</strong>g this observation with the Green-Gauss<br />

formula, we get the follow<strong>in</strong>g weak formulation of the ODE:<br />

Lemma (second <strong>in</strong>terpretation of the ODE on h(·) )<br />

Let u be a weak solution of the PDE on {x = h(t)}; let h ∈ W 2,∞ (0, T). Then<br />

h(·) verifies the ODE if <strong>and</strong> only if for all ξ ∈ D([0, T)), for all ψ ∈ D(R) such<br />

that ψ ≡ 1 on the set {x ∈ R : ∃t ∈ [0, T] such that h(t) = x}, there holds<br />

T<br />

−m h<br />

0<br />

′ (t)ξ ′ (t)dt = mh ′ T <br />

(0)ξ(0)+ uψξt +<br />

0 R<br />

u2<br />

2 ξψx<br />

<br />

+ u0ψξ(0).<br />

R<br />

<br />

.

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

Saved successfully!

Ooh no, something went wrong!