skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A positive and entropy-satisfying finite volume scheme for the Baer–Nunziato model

Abstract

We present a relaxation scheme for approximating the entropy dissipating weak solutions of the Baer–Nunziato two-phase flow model. This relaxation scheme is straightforwardly obtained as an extension of the relaxation scheme designed in for the isentropic Baer–Nunziato model and consequently inherits its main properties. To our knowledge, this is the only existing scheme for which the approximated phase fractions, phase densities and phase internal energies are proven to remain positive without any restrictive condition other than a classical fully computable CFL condition. For ideal gas and stiffened gas equations of state, real values of the phasic speeds of sound are also proven to be maintained by the numerical scheme. It is also the only scheme for which a discrete entropy inequality is proven, under a CFL condition derived from the natural sub-characteristic condition associated with the relaxation approximation. This last property, which ensures the non-linear stability of the numerical method, is satisfied for any admissible equation of state. We provide a numerical study for the convergence of the approximate solutions towards some exact Riemann solutions. The numerical simulations show that the relaxation scheme compares well with two of the most popular existing schemes available for the Baer–Nunziato model, namely Schwendeman–Wahle–Kapila'smore » Godunov-type scheme and Tokareva–Toro's HLLC scheme . The relaxation scheme also shows a higher precision and a lower computational cost (for comparable accuracy) than a standard numerical scheme used in the nuclear industry, namely Rusanov's scheme. Finally, we assess the good behavior of the scheme when approximating vanishing phase solutions.« less

Authors:
 [1];  [2];  [3]
  1. CMAP, École Polytechnique CNRS, UMR 7641, Route de Saclay, F-91128 Palaiseau cedex (France)
  2. EDF-R&D, Département MFEE, 6 Quai Watier, F-78401 Chatou Cedex (France)
  3. Université de Lyon, CNRS UMR 5208, Université Lyon 1, Institut Camille Jordan, 43 bd 11 novembre 1918, F-69622 Villeurbanne cedex (France)
Publication Date:
OSTI Identifier:
22622247
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 330; Other Information: Copyright (c) 2016 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; APPROXIMATIONS; COMPARATIVE EVALUATIONS; COMPUTERIZED SIMULATION; CONVERGENCE; DENSITY; ENTROPY; EQUATIONS OF STATE; FLOW MODELS; ISENTROPIC PROCESSES; MATHEMATICAL SOLUTIONS; MULTIPHASE FLOW; NUCLEAR INDUSTRY; NUMERICAL ANALYSIS; PARTIAL DIFFERENTIAL EQUATIONS; RELAXATION; SOUND WAVES; STABILITY; TWO-PHASE FLOW

Citation Formats

Coquel, Frédéric, E-mail: frederic.coquel@cmap.polytechnique.fr, Hérard, Jean-Marc, E-mail: jean-marc.herard@edf.fr, and Saleh, Khaled, E-mail: saleh@math.univ-lyon1.fr. A positive and entropy-satisfying finite volume scheme for the Baer–Nunziato model. United States: N. p., 2017. Web. doi:10.1016/J.JCP.2016.11.017.
Coquel, Frédéric, E-mail: frederic.coquel@cmap.polytechnique.fr, Hérard, Jean-Marc, E-mail: jean-marc.herard@edf.fr, & Saleh, Khaled, E-mail: saleh@math.univ-lyon1.fr. A positive and entropy-satisfying finite volume scheme for the Baer–Nunziato model. United States. doi:10.1016/J.JCP.2016.11.017.
Coquel, Frédéric, E-mail: frederic.coquel@cmap.polytechnique.fr, Hérard, Jean-Marc, E-mail: jean-marc.herard@edf.fr, and Saleh, Khaled, E-mail: saleh@math.univ-lyon1.fr. Wed . "A positive and entropy-satisfying finite volume scheme for the Baer–Nunziato model". United States. doi:10.1016/J.JCP.2016.11.017.
@article{osti_22622247,
title = {A positive and entropy-satisfying finite volume scheme for the Baer–Nunziato model},
author = {Coquel, Frédéric, E-mail: frederic.coquel@cmap.polytechnique.fr and Hérard, Jean-Marc, E-mail: jean-marc.herard@edf.fr and Saleh, Khaled, E-mail: saleh@math.univ-lyon1.fr},
abstractNote = {We present a relaxation scheme for approximating the entropy dissipating weak solutions of the Baer–Nunziato two-phase flow model. This relaxation scheme is straightforwardly obtained as an extension of the relaxation scheme designed in for the isentropic Baer–Nunziato model and consequently inherits its main properties. To our knowledge, this is the only existing scheme for which the approximated phase fractions, phase densities and phase internal energies are proven to remain positive without any restrictive condition other than a classical fully computable CFL condition. For ideal gas and stiffened gas equations of state, real values of the phasic speeds of sound are also proven to be maintained by the numerical scheme. It is also the only scheme for which a discrete entropy inequality is proven, under a CFL condition derived from the natural sub-characteristic condition associated with the relaxation approximation. This last property, which ensures the non-linear stability of the numerical method, is satisfied for any admissible equation of state. We provide a numerical study for the convergence of the approximate solutions towards some exact Riemann solutions. The numerical simulations show that the relaxation scheme compares well with two of the most popular existing schemes available for the Baer–Nunziato model, namely Schwendeman–Wahle–Kapila's Godunov-type scheme and Tokareva–Toro's HLLC scheme . The relaxation scheme also shows a higher precision and a lower computational cost (for comparable accuracy) than a standard numerical scheme used in the nuclear industry, namely Rusanov's scheme. Finally, we assess the good behavior of the scheme when approximating vanishing phase solutions.},
doi = {10.1016/J.JCP.2016.11.017},
journal = {Journal of Computational Physics},
number = ,
volume = 330,
place = {United States},
year = {Wed Feb 01 00:00:00 EST 2017},
month = {Wed Feb 01 00:00:00 EST 2017}
}
  • Numerical methods for the Baer-Nunziato (BN) two-phase flow model have attracted much attention in recent years. In this paper, we present a new gas kinetic scheme for the BN two-phase flow model containing non-conservative terms in the framework of finite volume method. In the view of microscopic aspect, a generalized Bhatnagar-Gross-Krook (BGK) model which matches with the BN model is constructed. Based on the integral solution of the generalized BGK model, we construct the distribution functions at the cell interface. Then numerical fluxes can be obtained by taking moments of the distribution functions, and non-conservative terms are explicitly introduced intomore » the construction of numerical fluxes. In this method, not only the complex iterative process of exact solutions is avoided, but also the non-conservative terms included in the equation can be handled well.« less
  • The structure is investigated of the Baer ideal of a finitely generated algebra of arbitrary finite signature over an arbitrary field or over a Noetherian commutative-associative ring satisfying a system of Capelli identities of order n + 1. It is proved that the length of the Baer chain of ideals in such an algebra is at most n. It is proved that the quotient of this algebra modulo the largest nilpotent ideal is representable.
  • Cited by 2
  • Recently a method for solving the Baer-Kouri-Levin-Tobocman integral equations was presented and applied to several collinear reactive systems. The method is based on using a distortion potential in the unperturbed Hamiltonian so that the translational part of the Green's function involves distorted elastic wave functions. In this work we report on a solution of these equations without a distortion potential so that the Green's function involves free (asymptotic) translational wave functions. The advantages of using these functions is that a certain amount of computer work is saved because some of the integrations can be carried out analytically. A disadvantage ismore » that the rate of convergence with respect to both vibrational and translational basis functions is slower than when the smooth distortion potential is employed. Two diagnostics of the accuracy of the results are found to be the size of the determinant of coefficients of the simultaneous algebraic equations and symmetry of the R-italic matrix. In addition, it is found that if a distortion potential is used, one should choose it so that the resulting perturbation is made small and of least extent possible. This accelerates convergence of the solution with respect to basis size.« less
  • This paper examines a class of explicit finite-difference advection schemes derived along the method of lines. An important application field is large-scale atmospheric transport. The paper therefore focuses on the demand of positivity. For the spatial discretization, attention is confined to conservative schemes using five points per direction. The fourth-order central scheme and the family of {kappa}-schemes, comprising the second-order central, the second-order upwind, and the third-order upwind biased, are studied. Positivity is enforced through flux limiting. It is concluded that the limited third-order upwind discretization is the best candidate from the four examined. For the time integration attention ismore » confined to a number of explicit Runge-Kutta methods of orders two to four. With regard to the demand of positivity, these integration methods turn out to behave almost equally and no best method could be identified. 16 refs., 4 figs., 4 tabs.« less