10 research outputs found
Convergence of numerical schemes for a conservation equation with convection and degenerate diffusion
The approximation of problems with linear convection and degenerate nonlinear diffusion, which arise in the framework of the transport of energy in porous media with thermodynamic transitions, is done using a θ-scheme based on the centred gradient discretisation method. The convergence of the numerical scheme is proved, although the test functions which can be chosen are restricted by the weak regularity hypotheses on the convection field, owing to the application of a discrete Gronwall lemma and a general result for the time translate in the gradient discretisation setting. Some numerical examples, using both the Control Volume Finite Element method and the Vertex Approximate Gradient scheme, show the role of θ for stabilising the scheme
Convergence of the the kinetic hydrostatic reconstruction scheme for the Saint Venant system with topography
We prove the convergence of the hydrostatic reconstruction scheme with kinetic numerical flux for the Saint Venant system with Lipschitz continuous topography. We use a recently derived fully discrete sharp entropy inequality with dissipation, that enables us to establish an estimate in the inverse of the square root of the space increment ∆x of the L 2 norm of the gradient of approximate solutions. By Diperna's method we conclude the strong convergence towards bounded weak entropy solutions
Convergence of the kinetic hydrostatic reconstruction scheme for the Saint Venant system with topography
International audienceWe prove the convergence of the hydrostatic reconstruction scheme with kinetic numerical flux for the Saint Venant system with Lipschitz continuous topography. We use a recently derived fully discrete sharp entropy inequality with dissipation, that enables us to establish an estimate in the inverse of the square root of the space increment ∆x of the L 2 norm of the gradient of approximate solutions. By Diperna's method we conclude the strong convergence towards bounded weak entropy solutions
A multi well-balanced scheme for the shallow water MHD system with topography
International audienceThe shallow water magnetohydrodynamic system involves several families of physically relevant steady states. In this paper we design a well-balanced numerical scheme for the one-dimensional shallow water magnetohydrodynamic system with topography, that resolves exactly a large range of steady states. Two variants are proposed with slightly different families of preserved steady states. They are obtained by a generalized hydrostatic reconstruction algorithm involving the magnetic field and with a cutoff parameter to remove singularities. The solver is positive in height and semi-discrete entropy satisfying, which ensures the robustness of the method
A 5-wave relaxation solver for the shallow water MHD system
International audienceThe shallow water magnetohydrodynamic system describes the thin layer evolution of the solar tachocline. It is obtained from the three dimensional incompressible magnetohydrodynamic system similarly as the classical shallow water system is obtained from the incompressible Navier-Stokes equations. The system is hyperbolic and has two additional waves with respect to the shallow water system, the Alfven waves. These are linearly degenerate, and thus do not generate dissipation. In the present work we introduce a 5-wave approximate Riemann solver for the shallow water magnetohydrodynamic system, that has the property to be non dissipative on Alfven waves. It is obtained by solving a relaxation system of Suliciu type, and is similar to HLLC type solvers. The solver is positive and entropy satisfying, ensuring its robustness. It has sharp wave speeds, and does not involve any iterative procedure
Modelling and entropy satisfying relaxation scheme for the nonconservative bitemperature Euler system with transverse magnetic field
The present paper concerns the study of the nonconservative bitem-perature Euler system with transverse magnetic field. We firstly introduce an underlying conservative kinetic model coupled to Maxwell equations. The nonconservative bitemperature Euler system with transverse magnetic field is then established from this kinetic model by hydrodynamic limit. Next we present the derivation of a finite volume method to approximate weak solutions. It is obtained by solving a relaxation system of Suliciu type, and is similar to HLLC type solvers. The solver is shown in particular to preserve positivity of density and internal energies. Moreover we use a local minimum entropy principle to prove discrete entropy inequalities, ensuring the robustness of the scheme
Convergence of numerical schemes for a conservation equation with convection and degenerate diffusion
The approximation of problems with linear convection and degenerate nonlinear diffusion, which arise in the framework of the transport of energy in porous media with thermodynamic transitions, is done using a θ-scheme based on the centred gradient discretisation method. The convergence of the numerical scheme is proved, although the test functions which can be chosen are restricted by the weak regularity hypotheses on the convection field, owing to the application of a discrete Gronwall lemma and a general result for the time translate in the gradient discretisation setting. Some numerical examples, using both the Control Volume Finite Element method and the Vertex Approximate Gradient scheme, show the role of θ for stabilising the scheme
Modelling and entropy satisfying relaxation scheme for the nonconservative bitemperature Euler system with transverse magnetic field
The present paper concerns the study of the nonconservative bitem-perature Euler system with transverse magnetic field. We firstly introduce an underlying conservative kinetic model coupled to Maxwell equations. The nonconservative bitemperature Euler system with transverse magnetic field is then established from this kinetic model by hydrodynamic limit. Next we present the derivation of a finite volume method to approximate weak solutions. It is obtained by solving a relaxation system of Suliciu type, and is similar to HLLC type solvers. The solver is shown in particular to preserve positivity of density and internal energies. Moreover we use a local minimum entropy principle to prove discrete entropy inequalities, ensuring the robustness of the scheme