Rigorous Multicomponent Reactive Separations Modelling : Complete Consideration of Reaction-Diffusion Phenomena

This paper gives the first step of the development of a rigorous multicomponent reactive separation model. Such a model is highly essential to further the optimization of acid gases removal plants (CO2 capture, gas treating, etc.) in terms of size and energy consumption, since chemical solvents are conventionally used.Firstly, two main modelling approaches are presented: the equilibrium-based and the rate-based approaches. Secondly, an extended rate-based model with rigorous modelling methodology for diffusion-reaction phenomena is proposed. The film theory and the generalized Maxwell-Stefan equations are used in order to characterize multicomponent interactions. The complete chain of chemical reactions is taken into account. The reactions can be kinetically controlled or at chemical equilibrium, and they are considered for both liquid film and liquid bulk. Thirdly, the method of numerical resolution is described. Coupling the generalized Maxwell-Stefan equations with chemical equilibrium equations leads to a highly non-linear Differential-Algebraic Equations system known as DAE index 3. The set of equations is discretized with finite-differences as its integration by Gear method is complex. The resulting algebraic system is resolved by the Newton- Raphson method. Finally, the present model and the associated methods of numerical resolution are validated for the example of esterification of methanol. This archetype non-electrolytic system permits an interesting analysis of reaction impact on mass transfer, especially near the phase interface. The numerical resolution of the model by Newton-Raphson method gives good results in terms of calculation time and convergence. The simulations show that the impact of reactions at chemical equilibrium and that of kinetically controlled reactions with high kinetics on mass transfer is relatively similar. Moreover, the Fick’s law is less adapted for multicomponent mixtures where some abnormalities such as counter-diffusion take place

The Stefan problem with variable thermophysical properties and phase change temperature

In this paper we formulate a Stefan problem appropriate when the thermophysical properties are distinct in each phase and the phase-change temperature is size or velocity dependent. Thermophysical properties invariably take different values in different material phases but this is often ignored for mathematical simplicity. Size and velocity dependent phase change temperatures are often found at very short length scales, such as nanoparticle melting or dendrite formation; velocity dependence occurs in the solidification of supercooled melts. To illustrate the method we show how the governing equations may be applied to a standard one-dimensional problem and also the melting of a spherically symmetric nanoparticle. Errors which have propagated through the literature are highlighted. By writing the system in non-dimensional form we are able to study the large Stefan number formulation and an energy-conserving one-phase reduction. The results from the various simplifications and assumptions are compared with those from a finite difference numerical scheme. Finally, we briefly discuss the failure of Fourier's law at very small length and time-scales and provide an alternative formulation which takes into account the finite time of travel of heat carriers (phonons) and the mean free distance between collisions.Comment: 39 pages, 5 figure

The pear-shaped fate of an ice melting front

A fluid-structure interaction problem with the melting of water around a heated horizontal circular cylinder is analysed with numerical simulations. Dynamic meshing was used for evolving the flow domain in time as the melting front extended radially outward from the cylinder; a node shuffle algorithm was used to retain mesh quality across the significant mesh deformation. We simulated one case above the density inversion point of water and one case below, yielding pear-shaped melting fronts due to thermal plumes either rising or falling from the cylinder, respectively. Results were compared with previous experimental studies and the melting front profiles matched reasonably well and melting rates were in agreement. We confirm that natural convection plays a significant role in the transport of energy as the melt zone increases, and needs to be considered for accurately modelling phase change under these conditions.Comment: Accepted for the 12th International Conference on CFD in Oil & Gas, Metallurgical and Process Industries. SINTEF, Trondheim, Norway. May 30th - June 1st, 201

Drug diffusion from polymeric delivery devices: a problem with two moving boundaries

An existing model for solvent penetration and drug release from a spherically-shaped polymeric drug delivery device is revisited. The model has two moving boundaries, one that describes the interface between the glassy and rubbery states of polymer, and another that defines the interface between the polymer ball and the pool of solvent. The model is extended so that the nonlinear diffusion coefficient of drug explicitly depends on the concentration of solvent, and the resulting equations are solved numerically using a front-fixing transformation together with a finite difference spatial discretisation and the method of lines. We present evidence that our scheme is much more accurate than a previous scheme. Asymptotic results in the small-time limit are presented, which show how the use of a kinetic law as a boundary condition on the innermost moving boundary dictates qualitative behaviour, the scalings being very different to the similar moving boundary problem that arises from modelling the melting of an ice ball. The implication is that the model considered here exhibits what is referred to as ``non-Fickian'' or Case II diffusion which, together with the initially constant rate of drug release, has certain appeal from a pharmaceutical perspective