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 />
Ma<strong>in</strong> results: Auxiliary Problem 1 is well posed<br />
For the <strong>Burgers</strong>-with-Dirac-at-zero <strong>mo<strong>de</strong>l</strong> , we apply the mach<strong>in</strong>ery<br />
<strong>de</strong>veloped for conservation laws with discont<strong>in</strong>uous flux (adapted<br />
entropies, Baiti,Jenssen <strong>and</strong> Audusse,Perthame ; revisited <strong>and</strong><br />
generalized recently by BA., Karlsen, Risebro us<strong>in</strong>g the notion of<br />
admissibility germ ). The outcome is:<br />
– <strong>de</strong>f<strong>in</strong>ition(s) of entropy solutions<br />
– uniqueness, cont<strong>in</strong>uous <strong>de</strong>pen<strong>de</strong>nce (L 1 , L 1 loc<br />
<strong>de</strong>pen<strong>de</strong>nce) exactly as <strong>in</strong> the Kruzhkov <strong>theory</strong><br />
In addition, we f<strong>in</strong>d<br />
with doma<strong>in</strong> of<br />
– a priori L ∞ bounds <strong>and</strong> (more <strong>de</strong>licate) variation bounds<br />
– a strik<strong>in</strong>gly simple numerical method (monotone consistent f<strong>in</strong>ite<br />
volume scheme with a trick at the <strong>in</strong>terface )<br />
– convergence of the numerical scheme, existence .<br />
NB: the Riemann solver at the <strong>in</strong>terface was already <strong>de</strong>scribed by<br />
Lagoutière, Segu<strong>in</strong>, Takahashi , so a Godunov scheme could be<br />
constructed; but we seek to avoid us<strong>in</strong>g the Riemann solver because<br />
it is <strong>in</strong>tricate.