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 Classification: | 04. 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 | Size | Format | Existing users please Login |
---|---|---|---|---|
PP_SIAM-JSC_LaSpina_etal_2012.pdf | 1.89 MB | Adobe PDF |
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