Abstract
Chromatographic processes can be modeled by nonlinear, convection-dominated partial differential equations, together with nonlinear relations: the adsorption isotherms. In this paper we consider the nonlinear equilibrium dispersive (ED) model with adsorption isotherms of Langmuir type. We show that, in this case, the ED model can be written as a system of conservation laws when the dispersion coefficient vanishes. We also show that the function that relates the conserved variables and the physically observed concentrations of the components in the mixture is one to one and it admits a global inverse, which cannot be given explicitly, but can be adequately computed.As a result, fully conservative numerical schemes can be designed for the ED model in chromatography. We propose a Weighted-Essentially-non-Oscillatory second order IMEX scheme and describe the numerical issues involved in its application. Through a series of numerical experiments, we show that our scheme gives accurate numerical solutions which capture the sharp discontinuities present in the chromatographic fronts, with the same stability restrictions as in the purely hyperbolic case.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call