Hostname: page-component-8448b6f56d-m8qmq Total loading time: 0 Render date: 2024-04-18T19:29:39.993Z Has data issue: false hasContentIssue false

Fluid displacement between two parallel plates: a non-empirical model displaying change of type from hyperbolic to elliptic equations

Published online by Cambridge University Press:  29 October 2004

M. SHARIATI
Affiliation:
Department of Chemical Engineering, University of Southern California, Los Angeles, CA 90089-1211, USA Present address: Applied Biosystems, Foster City, CA, USA.
L. TALON
Affiliation:
Laboratoire Fluides Automatique et Systèmes Thermiques, Universités P. et M. Curie and Paris Sud, CNRS (UMR 7608) Bâtiment 502, Campus Universitaire, 91405 Orsay Cedex, France
J. MARTIN
Affiliation:
Laboratoire Fluides Automatique et Systèmes Thermiques, Universités P. et M. Curie and Paris Sud, CNRS (UMR 7608) Bâtiment 502, Campus Universitaire, 91405 Orsay Cedex, France
N. RAKOTOMALALA
Affiliation:
Laboratoire Fluides Automatique et Systèmes Thermiques, Universités P. et M. Curie and Paris Sud, CNRS (UMR 7608) Bâtiment 502, Campus Universitaire, 91405 Orsay Cedex, France
D. SALIN
Affiliation:
Laboratoire Fluides Automatique et Systèmes Thermiques, Universités P. et M. Curie and Paris Sud, CNRS (UMR 7608) Bâtiment 502, Campus Universitaire, 91405 Orsay Cedex, France
Y. C. YORTSOS
Affiliation:
Department of Chemical Engineering, University of Southern California, Los Angeles, CA 90089-1211, USA

Abstract

We consider miscible displacement between parallel plates in the absence of diffusion, with a concentration-dependent viscosity. By selecting a piecewise viscosity function, this can also be considered as ‘three-fluid’ flow in the same geometry. Assuming symmetry across the gap and based on the lubrication (‘equilibrium’) approximation, a description in terms of two quasi-linear hyperbolic equations is obtained. We find that the system is hyperbolic and can be solved analytically, when the mobility profile is monotonic, or when the mobility of the middle phase is smaller than its neighbours. When the mobility of the middle phase is larger, a change of type is displayed, an elliptic region developing in the composition space. Numerical solutions of Riemann problems of the hyperbolic system spanning the elliptic region, with small diffusion added, show good agreement with the analytical outside, but an unstable behaviour inside the elliptic region. In these problems, the elliptic region arises precisely at the displacement front. Crossing the elliptic region requires the solution of essentially an eigenvalue problem of the full higher-dimensional model, obtained here using lattice BGK simulations. The hyperbolic-to-elliptic change-of-type reflects the failing of the lubrication approximation, underlying the quasi-linear hyperbolic formalism, to describe the problem uniformly. The obtained solution is analogous to non-classical shocks recently suggested in problems with change of type.

Type
Papers
Copyright
© 2004 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)