A compressible mixture model with phase transition
Introduction
In this study, we propose a model for chemically reacting viscous fluid mixtures that may develop a transition between a liquid and a vapor phase. The mixture consists of constituents and is described by partial mass balance equations and a single equation of balance for the barycentric momentum. We exclusively consider isothermal evolutions. To describe the phase transition, we introduce an artificial phase field indicating the present phase by assigning the values 1 and to the liquid and the vapor phase, respectively. Within the transition layer between two adjacent phases, the phase field smoothly changes between 1 and . However, usually the transition layers are very thin leading to steep gradients of the phase field.
This model belongs to the class of diffuse interface models. An alternative model class, that likewise represents phase transitions in fluid mixtures, contains sharp interface models. From the modeling point of view, sharp interface models have a simpler physical basis than diffuse interface models. For this reason, there arises always the non-trivial question if the sharp interface limits of a given diffuse model lead to admissible sharp interface models. The main concern of this paper is a careful discussion of this problem.
While diffuse interface models solve partial differential equations in the transition region, sharp interface models deal with jump conditions across the interface between the phases. Sometimes the jump conditions are mixed with geometric partial differential equations.
For two phases without chemical reactions, our compressible model reduces to an Allen–Cahn/Navier–Stokes type model, which is quite similar to the model derived by Blesgen [1]. Blesgen’s model has been investigated analytically in [2], [3], where existence of strong local-in-time solutions and weak solutions has been shown.
We would like to emphasize that the thermodynamical approaches of Blesgen’s system and our derived model are different. For instance, the phase field variable in Blesgen’s model, which satisfies the Allen–Cahn equation, is a physical quantity, namely the local mass fraction of one phase, whereas in our proposed model simply indicates the present phase. In addition, the function that interpolates between the phases, which will be described later, is linear in Blesgen’s model. This means that, in contrast to our proposed model, in general the equilibria depend on the chosen interpolation function , which implies that Blesgen’s system runs into different equilibria. Our approach with as an artificial phase field variable permits to obtain physical meaningful jump conditions at the interface agreeing with classical laws of thermodynamics.
A modified version of Blesgen’s model can be found in [4]. In contrast to [1] and our introduced model, Witterstein [4] describes a mixture of two compressible fluids, which physically differ, exclusively by different Lamé coefficients which are assumed to depend on the phase field parameter and the mass density. Witterstein’s model [4] also differs in the choice of the free energy, which contains two length scales. Moreover, the minima of the double well potential in the free energy have to be of different heights. This implies that the energy can only be controlled for transition regions with fixed width but not in the sharp interface limit.
Related to our work are diffuse interface models for incompressible and quasi-incompressible fluids. A diffuse interface model of Navier–Stokes–Cahn–Hilliard type for two incompressible, viscous Newtonian fluids, having the same densities, has been introduced by Hohenberg and Halperin in [5]. That model has been modified in several thermodynamically consistent ways such that different densities are allowed, see e.g. [6], [7], [8]. For existence results of strong local-in-time solutions and weak solutions, we refer to [9], [10], [11]. A diffuse interface model for two incompressible constituents which permits the transfer of mass between the phases due to diffusion and phase transitions has been proposed in [12], [13]. The densities of the fluids may be different, which leads to quasi-incompressibility of the mixture.
Our newly introduced diffuse interface model is given by the following system of PDEs for () in , : where is the pressure, the temperature, the atomic mass of constituent , the Boltzmann constant, with . In addition, are the stoichiometric coefficients of possible chemical reactions, the affinities and , and the mobilities.
The work is organized as follows. In the upcoming section we derive the thermodynamically consistent model for multi-component flows with phase transitions and chemical reactions. Section 3 is devoted to the non-dimensionalization, the introduction of two interesting scaling regimes of the system and the setting of asymptotic analysis. Finally, in Sections 4 Sharp interface limit of the dissipative regime, 5 Sharp interface limit of the non-dissipative regime, we determine the sharp interface limits for the two different scaling regimes introduced previously. We would like to emphasize that Section 4.2 contains a conjecture on the incapability of viscous diffuse models to generate viscous sharp models.
Section snippets
Constituents and phases
We consider a fluid mixture consisting of constituents indexed by . The constituents may be subjected to chemical reactions. There are reactions, indexed by , of the general type Thus, there are forward (f) as well as backward (b) reactions. The constants and are positive integers and denotes the stoichiometric coefficient of constituent in the reaction .
The fluid
Asymptotic analysis
To avoid physically impossible scalings, we first nondimensionalize the system (2.22)–(2.25).
Low viscosity case
We start by defining outer solutions in the bulk phases. They are obtained by inserting (3.8) into the scaled equations and comparing the terms order by order. Definition 4.1 A tuple with where , is called an outer solution of the dissipative regime if
Sharp interface limit of the non-dissipative regime
As the non-dissipative regime is rather similar to the dissipative regime treated in the last section we only outline the differences. The equations satisfied by outer solutions are obtained by inserting (3.8) into the scaled equations and comparing the terms order by order. Definition 5.1 A tuple is called an outer solution of the non-dissipative regime if it is an outer solution of the dissipative regime.
Inner solutions are determined from the scaled system by changing coordinates
Acknowledgments
W.D., J.G. and C.K. would like to thank the German Research Foundation (DFG) for financial support of the project “Modeling and sharp interface limits of local and non-local generalized Navier–Stokes–Korteweg Systems”. J.G. was also supported by the EU FP7-REGPOT project “Archimedes Center for Modeling, Analysis and Computation”.
References (23)
A generalization of the Navier–Stokes equations to two-phase flows
J. Phys. D: Appl. Phys. 32
(1999)Strong solutions of the Navier–Stokes equations for a compressible fluid of Allen–Cahn type
Arch. Ration. Mech. Anal.
(2012)- et al.
Analysis of a phase-field model for two-phase compressible fluids
Math. Models Methods Appl. Sci.
(2010) Phase change flows with mass exchange
Adv. Math. Sci. Appl.
(2011)- et al.
Theory of dynamic critical phenomena
Rev. Modern Phys.
(1977) - et al.
Two-phase binary fluids and immiscible fluids described by an order parameter
Math. Models Methods Appl. Sci.
(1996) - et al.
Quasi-incompressible Cahn-Hilliard fluids and topological transitions
R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci.
(1998) - et al.
Thermodynamically consistent frame indifferent diffuse interface models for incompressible two-phase flows with different densities
Math. Models Methods Appl. Sci.
(2012) Existence of weak solutions for a diffuse interface model for viscous, incompressible fluids with general densities
Comm. Math. Phys.
(2009)Strong well-posedness of a diffuse interface model for a viscous, quasi-incompressible two-phase flow
SIAM J. Math. Anal.
(2012)
Existence of weak solutions for a diffuse interface model for two-phase flows of incompressible fluids with different densities
J. Math. Fluid Mech.
Cited by (14)
A new method to discretize a model for isothermal flow with a multi-component equation of state
2023, Journal of Computational and Applied MathematicsCitation Excerpt :The advantage of a phase field model is that the phase transition is handled by the PDE for the phase field quantity. Numerical methods for a sub-model of a phase field model derived by Dreyer, Giesselmann and Kraus [12] are the main focus of this paper. We will introduce both models in the following section.
New insights on the interfacial tension of electrochemical interfaces and the Lippmann equation
2018, European Journal of Applied MathematicsExplicit implicit domain splitting for two phase flows with phase transition
2023, Physics of Fluids