# American Institute of Mathematical Sciences

March  2020, 15(1): 143-169. doi: 10.3934/nhm.2020007

## The Lax-Friedrichs scheme for interaction between the inviscid Burgers equation and multiple particles

 MiraCosta College, 3333 Manchester Avenue, Cardiff-by-the-Sea, CA 92007-1516, USA

Received  July 2019 Revised  October 2019 Published  December 2019

We propose a finite difference method based on the Lax-Friedrichs scheme for a model of interaction between multiple solid particles and an inviscid fluid. The single-particle version has been studied extensively during the past decade. The model studied here consists of the inviscid Burgers equation with multiple nonconservative moving source terms that are singular and account for drag force interaction between the fluid and the particles. Each particle trajectory satisfies a differential equation that ensures conservation of momentum of the entire system. To deal with the singular source terms we discretize a model that associates with each particle an advection PDE whose solution is a shifted Heaviside function. This alternative model is well known but has not previously been used in numerical methods. We propose a definition of entropy solution which directly generalizes the previously defined single-particle notion of entropy solution. We prove convergence (along a subsequence) of the Lax-Friedrichs approximations, and also prove that if the set of times where the particle paths intersect has Lebesgue measure zero, then the limit is an entropy solution. We also propose a higher resolution version of the scheme, based on MUSCL processing, and present the results of numerical experiments.

Citation: John D. Towers. The Lax-Friedrichs scheme for interaction between the inviscid Burgers equation and multiple particles. Networks & Heterogeneous Media, 2020, 15 (1) : 143-169. doi: 10.3934/nhm.2020007
##### References:
 [1] N. Aguillon, F. Lagoutière and N. Seguin, Convergence of finite volume schemes for the coupling between the inviscid Burgers equation and a particle, Math. Comp., 86 (2017), 157-196.  doi: 10.1090/mcom/3082.  Google Scholar [2] N. Aguillon, Riemann problem for a particle-fluid coupling, Math. Models Methods Appl. Sci., 25 (2015), 39-78.  doi: 10.1142/S0218202515500025.  Google Scholar [3] N. Aguillon, Numerical simulations of a fluid-particle coupling, in Finite Volumes for Complex Applications Ⅶ-Elliptic, Parabolic and Hyperbolic Problems. Springer Proceedings in Mathematics & Statistics (eds. J. Fuhrmann, M. Ohlberger M. and C. Rohde), Springer, 78 (2014), 759–767. doi: 10.1007/978-3-319-05591-6_76.  Google Scholar [4] B. Andreianov, K. H. Karlsen and N. H. Risebro, A theory of $L^1$-dissipative solvers for scalar conservation laws with discontinuous flux, Arch. Ration. Mech. Anal., 201 (2011), 27-86.  doi: 10.1007/s00205-010-0389-4.  Google Scholar [5] B. Andreianov, F. Lagoutière, N. Seguin and T. Takahashi, Small solids in an inviscid fluid, Netw. Heterog. Media, 5 (2010), 385-404.  doi: 10.3934/nhm.2010.5.385.  Google Scholar [6] B. Andreianov, F. Lagoutière, N. Seguin and T. Takahashi, Well-posedness for a one-dimensional fluid-particle interaction model, SIAM J. Math. Anal., 46 (2014), 1030-1052.  doi: 10.1137/130907963.  Google Scholar [7] B. Andreianov and N. Seguin, Analysis of a Burgers equation with singular resonant source term and convergence of well-balanced schemes, Discrete Contin. Dyn. Syst., 32 (2012), 1939-1964.  doi: 10.3934/dcds.2012.32.1939.  Google Scholar [8] M. Crandall and A. Majda, Monotone difference approximations for scalar conservation laws, Math. Comp., 34 (1980), 1-21.   Google Scholar [9] H. Holden and N. H. Risebro, Front Tracking for Hyperbolic Conservation Laws, Springer-Verlag, New York, 2002. doi: 10.1007/978-3-642-56139-9.  Google Scholar [10] C. Klingenberg, J. Klotzky and N. Seguin, On well-posedness for a multi-particle fluid model, in Theory, Numerics and Applications of Hyperbolic Problems Ⅱ. HYP 2016. Springer Proceedings in Mathematics & Statistics, (eds. C. Klingenberg and M. Westdickenberg), Springer, 237 (2018), 167–177.  Google Scholar [11] F. Lagoutière, N. Seguin and T. Takahashi, A simple 1D model of inviscid fluid-solid interaction, J. Differ. Equ., 245 (2008), 3503-3544.  doi: 10.1016/j.jde.2008.03.011.  Google Scholar [12] R. J. Leveque, Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, Cambridge, UK, 2002.  doi: 10.1017/CBO9780511791253.  Google Scholar [13] E. Y. Panov, Existence of strong traces for quasi-solutions of multidimensional conservation laws, J. Hyperbolic Differ. Equ., 4 (2007), 729-770.  doi: 10.1142/S0219891607001343.  Google Scholar [14] J. Towers, A fixed grid, shifted stencil scheme for inviscid fluid-particle interaction, Appl. Numer. Math., 110 (2016), 26-40.  doi: 10.1016/j.apnum.2016.08.002.  Google Scholar

show all references

##### References:
 [1] N. Aguillon, F. Lagoutière and N. Seguin, Convergence of finite volume schemes for the coupling between the inviscid Burgers equation and a particle, Math. Comp., 86 (2017), 157-196.  doi: 10.1090/mcom/3082.  Google Scholar [2] N. Aguillon, Riemann problem for a particle-fluid coupling, Math. Models Methods Appl. Sci., 25 (2015), 39-78.  doi: 10.1142/S0218202515500025.  Google Scholar [3] N. Aguillon, Numerical simulations of a fluid-particle coupling, in Finite Volumes for Complex Applications Ⅶ-Elliptic, Parabolic and Hyperbolic Problems. Springer Proceedings in Mathematics & Statistics (eds. J. Fuhrmann, M. Ohlberger M. and C. Rohde), Springer, 78 (2014), 759–767. doi: 10.1007/978-3-319-05591-6_76.  Google Scholar [4] B. Andreianov, K. H. Karlsen and N. H. Risebro, A theory of $L^1$-dissipative solvers for scalar conservation laws with discontinuous flux, Arch. Ration. Mech. Anal., 201 (2011), 27-86.  doi: 10.1007/s00205-010-0389-4.  Google Scholar [5] B. Andreianov, F. Lagoutière, N. Seguin and T. Takahashi, Small solids in an inviscid fluid, Netw. Heterog. Media, 5 (2010), 385-404.  doi: 10.3934/nhm.2010.5.385.  Google Scholar [6] B. Andreianov, F. Lagoutière, N. Seguin and T. Takahashi, Well-posedness for a one-dimensional fluid-particle interaction model, SIAM J. Math. Anal., 46 (2014), 1030-1052.  doi: 10.1137/130907963.  Google Scholar [7] B. Andreianov and N. Seguin, Analysis of a Burgers equation with singular resonant source term and convergence of well-balanced schemes, Discrete Contin. Dyn. Syst., 32 (2012), 1939-1964.  doi: 10.3934/dcds.2012.32.1939.  Google Scholar [8] M. Crandall and A. Majda, Monotone difference approximations for scalar conservation laws, Math. Comp., 34 (1980), 1-21.   Google Scholar [9] H. Holden and N. H. Risebro, Front Tracking for Hyperbolic Conservation Laws, Springer-Verlag, New York, 2002. doi: 10.1007/978-3-642-56139-9.  Google Scholar [10] C. Klingenberg, J. Klotzky and N. Seguin, On well-posedness for a multi-particle fluid model, in Theory, Numerics and Applications of Hyperbolic Problems Ⅱ. HYP 2016. Springer Proceedings in Mathematics & Statistics, (eds. C. Klingenberg and M. Westdickenberg), Springer, 237 (2018), 167–177.  Google Scholar [11] F. Lagoutière, N. Seguin and T. Takahashi, A simple 1D model of inviscid fluid-solid interaction, J. Differ. Equ., 245 (2008), 3503-3544.  doi: 10.1016/j.jde.2008.03.011.  Google Scholar [12] R. J. Leveque, Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, Cambridge, UK, 2002.  doi: 10.1017/CBO9780511791253.  Google Scholar [13] E. Y. Panov, Existence of strong traces for quasi-solutions of multidimensional conservation laws, J. Hyperbolic Differ. Equ., 4 (2007), 729-770.  doi: 10.1142/S0219891607001343.  Google Scholar [14] J. Towers, A fixed grid, shifted stencil scheme for inviscid fluid-particle interaction, Appl. Numer. Math., 110 (2016), 26-40.  doi: 10.1016/j.apnum.2016.08.002.  Google Scholar
Example 8.3. Particle trajectories using basic scheme (upper plot) and MUSCL (lower plot). Both the true (thick line) and approximate (thin line) trajectories are plotted. For the MUSCL scheme (lower plot) the true and approximate trajectories are visually indistinguishable at this level of discretization. $\Delta x = 1.95 \times 10^{-5}$, $\mu = .25$, 102401 time steps
Example 8.4. Basic scheme (left) and MUSCL (right). The horizontal axis represents $x$, and the vertical axis represents $t$. Top level plots: $\Delta x_1 = 3.75 \times 10^{-4}$. Middle level plots: $\Delta x_2 = {1\over 2} \Delta x_1$. Bottom level plots: $\Delta x_3 = {1\over 4} \Delta x_1$. $\mu = .125$ for all plots
Example 8.1. Top: Fluid velocity $u$ at $t = 1$. Exact solution is solid line, with sharp corners. Bottom: Particle position error vs. time. Basic scheme (left plots) and MUSCL scheme (right plots). $\Delta x = .0025$ (dashed line), and $\Delta x = .00125$ (solid line). Both approximations used $\mu = .25$
Example 8.2. Fluid velocity $u$ at $t = 1$. Basic scheme (left plots) and MUSCL scheme (right plots). Exact solution (dashed line) and approximate solution (solid line). Top plots used $\Delta x = .005$, bottom plots used $\Delta x = .000625$. All approximations used $\mu = .25$. A spurious kink is visible. Its magnitude diminishes with grid refinement
Example 8.3. Solution $u$ using basic scheme at $t = .125$ (upper left), and using MUSCL (upper right). True solution (dashed line) and approximate solution (solid line). Both upper plots computed with $\Delta x = .00325$, $\mu = .25$. The lower plots show the error in $u$ in discrete $L^1$ norm as a function of time using the basic scheme (lower left) and MUSCL scheme (lower right). Uses $\Delta x = .00325$ and $\Delta x = .001625$, $\mu = .25$
 [1] Paolo Baiti, Alberto Bressan, Helge Kristian Jenssen. Instability of travelling wave profiles for the Lax-Friedrichs scheme. Discrete & Continuous Dynamical Systems - A, 2005, 13 (4) : 877-899. doi: 10.3934/dcds.2005.13.877 [2] Hermano Frid. Invariant regions under Lax-Friedrichs scheme for multidimensional systems of conservation laws. Discrete & Continuous Dynamical Systems - A, 1995, 1 (4) : 585-593. doi: 10.3934/dcds.1995.1.585 [3] Boris Andreianov, Nicolas Seguin. Analysis of a Burgers equation with singular resonant source term and convergence of well-balanced schemes. Discrete & Continuous Dynamical Systems - A, 2012, 32 (6) : 1939-1964. doi: 10.3934/dcds.2012.32.1939 [4] Navnit Jha. Nonpolynomial spline finite difference scheme for nonlinear singuiar boundary value problems with singular perturbation and its mechanization. Conference Publications, 2013, 2013 (special) : 355-363. doi: 10.3934/proc.2013.2013.355 [5] Tetsuya Ishiwata, Kota Kumazaki. Structure preserving finite difference scheme for the Landau-Lifshitz equation with applied magnetic field. Conference Publications, 2015, 2015 (special) : 644-651. doi: 10.3934/proc.2015.0644 [6] Vanessa Baumgärtner, Simone Göttlich, Stephan Knapp. Feedback stabilization for a coupled PDE-ODE production system. Mathematical Control & Related Fields, 2019, 0 (0) : 0-0. doi: 10.3934/mcrf.2020003 [7] Mohamed Alahyane, Abdelilah Hakim, Amine Laghrib, Said Raghay. Fluid image registration using a finite volume scheme of the incompressible Navier Stokes equation. Inverse Problems & Imaging, 2018, 12 (5) : 1055-1081. doi: 10.3934/ipi.2018044 [8] Ronald E. Mickens. A nonstandard finite difference scheme for the drift-diffusion system. Conference Publications, 2009, 2009 (Special) : 558-563. doi: 10.3934/proc.2009.2009.558 [9] Matúš Tibenský, Angela Handlovičová. Convergence analysis of the discrete duality finite volume scheme for the regularised Heston model. Discrete & Continuous Dynamical Systems - S, 2018, 0 (0) : 0-0. doi: 10.3934/dcdss.2020226 [10] David Bourne, Howard Elman, John E. Osborn. A Non-Self-Adjoint Quadratic Eigenvalue Problem Describing a Fluid-Solid Interaction Part II: Analysis of Convergence. Communications on Pure & Applied Analysis, 2009, 8 (1) : 143-160. doi: 10.3934/cpaa.2009.8.143 [11] Matania Ben–Artzi, Joseph Falcovitz, Jiequan Li. The convergence of the GRP scheme. Discrete & Continuous Dynamical Systems - A, 2009, 23 (1&2) : 1-27. doi: 10.3934/dcds.2009.23.1 [12] Johannes Elschner, George C. Hsiao, Andreas Rathsfeld. An inverse problem for fluid-solid interaction. Inverse Problems & Imaging, 2008, 2 (1) : 83-120. doi: 10.3934/ipi.2008.2.83 [13] Peter Monk, Virginia Selgas. An inverse fluid--solid interaction problem. Inverse Problems & Imaging, 2009, 3 (2) : 173-198. doi: 10.3934/ipi.2009.3.173 [14] Sabine Eisenhofer, Messoud A. Efendiev, Mitsuharu Ôtani, Sabine Schulz, Hans Zischka. On an ODE-PDE coupling model of the mitochondrial swelling process. Discrete & Continuous Dynamical Systems - B, 2015, 20 (4) : 1031-1057. doi: 10.3934/dcdsb.2015.20.1031 [15] Thierry Cazenave, Yvan Martel, Lifeng Zhao. Finite-time blowup for a Schrödinger equation with nonlinear source term. Discrete & Continuous Dynamical Systems - A, 2019, 39 (2) : 1171-1183. doi: 10.3934/dcds.2019050 [16] Youshan Tao, J. Ignacio Tello. Nonlinear stability of a heterogeneous state in a PDE-ODE model for acid-mediated tumor invasion. Mathematical Biosciences & Engineering, 2016, 13 (1) : 193-207. doi: 10.3934/mbe.2016.13.193 [17] Ciro D'Apice, Peter I. Kogut, Rosanna Manzo. On relaxation of state constrained optimal control problem for a PDE-ODE model of supply chains. Networks & Heterogeneous Media, 2014, 9 (3) : 501-518. doi: 10.3934/nhm.2014.9.501 [18] Maria Laura Delle Monache, Paola Goatin. A front tracking method for a strongly coupled PDE-ODE system with moving density constraints in traffic flow. Discrete & Continuous Dynamical Systems - S, 2014, 7 (3) : 435-447. doi: 10.3934/dcdss.2014.7.435 [19] Wei Qu, Siu-Long Lei, Seak-Weng Vong. A note on the stability of a second order finite difference scheme for space fractional diffusion equations. Numerical Algebra, Control & Optimization, 2014, 4 (4) : 317-325. doi: 10.3934/naco.2014.4.317 [20] Wen Li, Song Wang. Pricing American options under proportional transaction costs using a penalty approach and a finite difference scheme. Journal of Industrial & Management Optimization, 2013, 9 (2) : 365-389. doi: 10.3934/jimo.2013.9.365

2018 Impact Factor: 0.871

## Tools

Article outline

Figures and Tables