Please use this identifier to cite or link to this item: http://hdl.handle.net/2122/8273
Authors: La Spina, G.* 
De' Michieli Vitturi, M.* 
Title: High resolution finite volume central schemes for a compressibile two-phase model
Journal: SIAM journal on scientific computing 
Series/Report no.: /34 (2012)
Publisher: Society for Industrial and Applied Mathematics
Issue Date: 2012
DOI: 10.1137/12087089X
Keywords: High-resolution central schemes
MUSCL-Hancock method
theory of thermodynamically compatible system of conservation laws
compressible two-phase
Subject Classification04. Solid Earth::04.08. Volcanology::04.08.04. Thermodynamics 
05. General::05.01. Computational geophysics::05.01.05. Algorithms and implementation 
05. General::05.05. Mathematical geophysics::05.05.99. General or miscellaneous 
Abstract: A modi_cation of the Kurganov, Noelle, Petrova central-upwind scheme [A. Kurganov et al., SIAM J. Sci. Comput., 23 (2001), pp. 707{740] for hyperbolic systems of conservation laws is presented. In this work, the numerical scheme is applied to a single-temperature model for compressible two-phase ow with pressure and velocity relaxations [E. Romenski et al., J. Sci. Comput., 42 (2010), pp. 68{95]. The system of governing equations of this model are expressed in conservative form, which is the necessary condition to use a central scheme. The numerical scheme presented is not based on the complete characteristic decomposition, but only on the information about the local speeds of propagation given by the maximum and minimum eigenvalue of the Jacobian of the uxes. We propose to use the numerical ux formulation of the central-upwind scheme in conjunction with a second-order reconstruction of the primitive variables and the MUSCL-Hancock method, where the boundary extrapolated values are evolved by half time step before the computation of the numerical uxes. To investigate the accuracy and robustness of the proposed scheme, two 1D Riemann-problems of an air/water mixture and a 2D shock-bubble-interaction problem are presented. Furthermore, a detailed comparison with the second order GFORCE scheme and the _rst order Lax-Friedrichs scheme is shown. To integrate the source terms an operator splitting approach is used and, under suitable conditions, it is shown that this integration can be computed analytically.
Appears in Collections:Article published / in press

Files in This Item:
File Description SizeFormat Existing users please Login
PP_SIAM-JSC_LaSpina_etal_2012.pdf1.89 MBAdobe PDF
Show full item record

WEB OF SCIENCETM
Citations

13
checked on Feb 10, 2021

Page view(s) 10

348
checked on Mar 16, 2024

Download(s) 50

63
checked on Mar 16, 2024

Google ScholarTM

Check

Altmetric