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 ...
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 />
.