Fast and Robust Computation of the Matrix Absolute Value Function. Application to Roe Solver for the Numerical Simulation of Two-Phase Flow Models
- Commissariat a l'Energie Atomique - CEA (France)
- Ecole Centrale Paris, Grande Voie des Vignes, F-92295 Chatenay-Malabry Cedex (France)
The application of the generalized Roe scheme to the numerical simulation of two-phase flow models requires a fast and robust computation of the absolute value of the system matrix. In several models such as the two-fluid model or a general multi-field model, this matrix has a non trivial Eigen-structure and the Eigen decomposition is often ill conditioned. We give two general algorithms avoiding the diagonalization process: an iterative computation, which turns out to be an exact computation, and an interpolation algorithm which is faster and can handle the case of complex eigenvalues. The knowledge of the characteristic polynomial gives us an easy access to the eigenvalues but however, the iterative scheme can be used with only estimates of the eigenvalues, using for example Gershgorin's disk localization. We finally show some numerical results of two-fluid model simulations involving interfacial pressure and virtual mass force models. (authors)
- Research Organization:
- The ASME Foundation, Inc., Three Park Avenue, New York, NY 10016-5990 (United States)
- OSTI ID:
- 20997024
- Resource Relation:
- Conference: 14. International conference on nuclear engineering (ICONE 14), Miami - Florida (United States), 17-20 Jul 2006; Other Information: Country of input: France
- Country of Publication:
- United States
- Language:
- English
Similar Records
A Riemann solver for single-phase and two-phase shallow flow models based on relaxation. Relations with Roe and VFRoe solvers
Singular value decomposition utilizing parallel algorithms on graphical processors